补充3(Matlab程序设计).docVIP

Matlab程序设计 Matlab有两种工作方式： 一种是由键盘逐行输入命令（每行可以有多个），Matlab可以立即执行并显示结果。 另一种是执行一个命令语句集文件，它自动按照文件中排好的命令和语句，顺序执行并显示结果。这种命令语句集文件必须使用扩展名“.m”，因此称为M-文件。可以由一系列命令语句组成，也可从中调用其它M-文件，甚至可以递归调用自身。M-文件中的语句对整个Matlab工作空间里的变量起作用，特别适宜于数据分析、解方程、设计等需要使用大量语句的用户，从键盘逐条输入语句效率太低。 2.4.1 M-文件 M-文件分为两类： 一类是程序文件，它自动执行一长串语句和命令；另一类是函数文件，它把用户自建的函数加入到现存的函数库中。 程序文件和函数文件都是用普通的ASCII码写成，用户可以用熟悉的编辑器或文字处理器创建或修改。 文本编辑器通常用颜色来区分程序内容的类别： 绿色：注释部分，程序并不执行。 黑色：程序主体部分。 红色：属性值的设定。 蓝色：控制流程，如for、if…then等语句。 1 建立新的M-文件 （1）通过菜单栏建立；（2）通过工具栏建立；（3）通过命令窗输入命令edit建立。 2 打开已有的M-文件 （1）菜单中的open命令；（2）工具栏中的打开按钮；（3）在命令窗输入：edit 命令； 2.4.2 数据的输入输出 1 imput函数 a=input(输入一个矩阵) 2 disp函数 a=[1 2 3;4 5 6;7 8 9]; disp(a) 例：求一元二次方程的根，在屏幕上提示输入a、b、c的值，输出计算结果在屏幕上。 a=input(a=);b=input(b=);c=input(c=); d=b*b-4*a*c; x1=(-b+sqrt(d))/2*a;x2=(-b-sqrt(d))/2*a; disp([x1,x2]) 2.4.3 选择结构 1 单分支if语句 if 条件 语句组 end 例 a=input(a=) if a0 disp(a) end 2 双分支语句 if 条件 语句组1 else 语句组2 end 例： x=input(x=); if x==10 y=cos(x+1)+sqrt(x*x+1); else y=x*sqrt(x+sqrt(x)); end disp(y); 3 多分支if语句 if 条件 语句组1 elseif 语句组2 … else 语句组3 end 2.4.4 循环结构 1 循环语句 for i=2:2:100 s=s+1; end disp(s); 也可以写为a=2:2:100;for i=a 2 while语句 while expression statements end 例：s=0;n=1;while n=100? ? s=s+1/n;? ? n=n+1;ends 要注意避免死循环，如果执行中出现死循环，可以用命令Ctrl+C来终止出程序的执行；也可以直接在程序中加入命令break跳出循环。 2 循环嵌套 for n=1:5 for m=5:-1:1 A(n,m)=n^2+m^2; end disp(A) end 2.4.5 函数文件 function [输出形参表]=函数名（输入形参表） 例 计算n的阶乘，定义阶乘函数为M-文件 function [f]=jch(n) f=1; for i=2:n f=f*i; end end 在命令窗输入n=10;jch(n)或jch(10) 注意M-文件的存储名和地址问题。 调用M-文件的函数 function [f]=sjch(n) f=0; for i=1:n f=f+jch(n); end 可采用递归调用 function [f]=dgich(n) if n1 f=dgich(n-1)+jch(n); else f=0; end