基本结构、循环、线代.pptVIP

MATLAB软件包简介 数学软件的应用 MATLAB使用参考书 1。张志涌等. 精通MATLAB6.5，第1版. 北京航空航天大学出版社， 2003．3 2。薛定宇等. 高等应用数学问题的 MATLAB求解，第1版. 清华大学出版社， 2004．8 3。姜启源等. 数学实验，第1版. 高等教育 出版社，1999．7 MATLAB ----MATrix LABoratory的简写 MathWorks公司 特长：数值计算，也具有符号计算功能（采用MAPLE软件作为核心软件） MATLAB基本命令与结构 基本命令 循环与控制 符号运算 矩阵操作及运算 MATLAB命令窗口 2+2/3 syms a b; int(a^2+b,a) vpa(pi) vpa(pi,100) A=[1 2 3;2 4 5; 3 6 7] 基本命令和结构\basic1.m 求π的近似值 .M文件的编写和运行 1。打开M—文件编辑器编写 不看结果的语句加 ； 注释行前加 %，后面的字符显示为绿色 2。将文件存在自己的工作目录下 3。将命令窗口的工作目录选为自己的目录 4。在命令提示符后输入M文件名即可 MATLAB 语言的变量名规则 由一个字母引导，后面可以为其他字符 区分大小写 Abc ABc 有效 MYvar12, MY_Var12 和 MyVar12_ 错误的变量名 12MyVar, _MyVar12 MATLAB 的保留常量 eps, i, j, pi, NaN, Inf, i=sqrt(-1) lastwarn, lasterr 变量命名规则 MATLAB是区分变量大、小写的 数值型数据结构 双精度数值变量 IEEE标准，64位 (8字节)，11指数位，52数值位和一个符号位 double() 函数的转换 其他数据类型 uint8，常用于图像表示和处理，8位 int8(), int16(), int32(),uint16(), uint32() 直接赋值语句 例1 表示矩阵 基本语句结构 1、2、3分别变成第一二三行第四列，1、2、3、4分别变成一二三四列第四行 冒号表达式 这里，s1、s3 分别是起点、终点值；s2 为步长 例2 用不同的步距生成 (0,pi) 间向量 v5与v1的区别：最后会输出3.1415，而v1输出最后一个数是3.0000 数的基本运算 floor(a) 向下取整 ceil(a) 向上取整 round(a) 四舍五入 factor(n) 质因数分解 gcd(m,n) 最大公约数 lcm(m,n) 最小公倍数 MATLAB 语言流程控制 循环结构 for 结构 while 结构 例3 用循环求解 例4 用循环结构 求最小的 m 输出的s是10011，m是142，即所求 例5 求 %%该程序是用来看运行时间的%% 转移结构 例6 用循环结构 求最大的 m 开关结构 circle.m MATLAB 语言的函数的基本结构 函数调用语句 例6 用循环结构 求最大的 m function [m,s]=findsum(k) s=0; m=0; while(s=k) m=m+1; s=s+m; end; [s,m] findsum(10000) findsum.m 例7 求和 s = 0; for n = 0:1:100; y = factorial(n); s = s+1/y; end; vpa(s,40) function k=factorial(n) if n1 k=n*factorial(n-1); elseif any([0,1]==n) k=1; %%% 0!=1!=1 end appro_e 基本命令和结构\factorial.m \appro_e.m %%阶乘表示成一个函数%% 符号计算(symbolic calculation)工具箱 symbolic1.m 矩阵相关运算 matrix1.m poly_A.m det(A) % 矩阵的行列式 inv(A) % 矩阵的逆 rank(A) % 矩阵的秩 eig(A) % 矩阵的特征值 [v,d] = eig(sym(A)) % 矩阵的特征值v