MatLab3-程序设计1解析

1、MatLab &数学建模第三讲MatLab程序设计一、M文件M文件可以分为脚本文件(Script和函数文件(Function两种。1脚本文件:(1编写文档：点击MATLAB指令窗口上面最左端的图标，即新建文 件，就可打开MATLAB文件编辑器，像word 一样。用户即可在空白窗口中编写程 序。例如输入如下一段程序area. m的例子% 新建 M-file, area.m%计算一个球的体积r = input(Type radius 输入半径:;Area=pi*rA2;volume=(4/3*pi*rA3;fprintf(半径 The radius is %12.5fn,rfprintf(面积 T

2、he area of a circle is %12.5fn,Areafprintf(体积 The volume of a sphere is %12.5fn,volume1ATLAEkla |iil (in* fe|CfhiVvMl WfWfew3TUE DfrWlW| E 1MATL些却Ft3j t 3 34 F* 曹阳/IM t pictor*Ir- 1 S3 f 4/04/04 -W ILCtvrer*Lr-r i点击编辑调试器工具条图标，在弹出的Windows标准风格的 保存为”对话框中, 选择保存文件夹，键入新编文件名（如area点击【保存】键，就完成了文件保存。（2点击文件编辑器

3、上面工具条中的保存,命名（例如将上面的程序命名为picture , 然后保存。这样的文件就是M-文件:area.m（3运行:i在命令窗口中输入文件名（如上面的area然后执行。ii或直接在文件编辑器上面的工具条中找到debug里面的run （即运行，点击即可（或直接按F5。iii如果发现错误,在写好的程序里直接改正，然后再保存，再运行，注意一定要先保存，后运行iiu lA聿卷*S Tta radiui is 7 OTOSJ面樹 Tlbt u tf a 疋逛打.l& 12. 5t37 住IE Um velwt f shtrt ii 33. SLQ32函数文件：我们经常用到的像sin、cos、ex

4、p这样的一些函数都是 MATLAB 软件自身所带的函数，因此直接应用即可，但有时我们为了解决一些问题需要自己编 写函数。自己编写函数有两个基本要求i必须在MATLAB文件编辑器中编写,也是M-文件。ii函数名和文件名必须相同。例编写函数3222 6.30.05 3.14x x x f x x -+-=+-,计算 3(2(1(2f f f +(1打开MATLAB文件编辑器，即点图标，输入function y= fun1(x % 表示y是x的函数,x是自变量,funl是我们自己命名的函数名y=(xA3 - 2*xA2 + x - 6.3/(xA2 + 0.05*x - 3.14;然后保存。注:在自

5、己编写的函数前都要写上function(关键字，表示这是自己定义的函数。 funl表示函数名，那么最后保存文件也要命名为fun 1。在指令窗口运行以下指令：fun 1(1*fu n1(2+fu n1(3*fu n1(3ans =-12.6023例:计算阶层函数fact.mfunction output = fact(nif n = 1output = 1;return;endoutput = n *fact( n-1;二MATLAB控制流for循环结构For循环允许一组命令以固定的和预定的次数重复。For循环的一般形式是for x = arraycomma ndsend在for和end语句之间

6、的commands按数组中的每一列执行一次。【例】for n=1:10x(n=si n(n *pi/10;end? xx =Columns 1 through 70.3090 0.5878 0.8090 0.9511 1.0000 0.9511 0.8090Columns 8 through 100.5878 0.3090 0.0000For循环可按需要嵌套。n=0:1:10;for i=1:11for j=1:11y(i=si n(n (i;n(j=n(j*10;endendy =Columns 1 through 70 -0.5064 0.9300 -0.8027 -0.1425 -0.9

7、765 -0.5118 Columns 8 through 110.8586 -0.9957 0.9917 0.9287n =1.0e+012 *Columns 1 through 70 0.1000 0.2000 0.3000 0.4000 0.5000 0.6000 Columns 8 through 110.7000 0.8000 0.9000 1.0000可以利用break命令跳出for循环【例】一个简单的for循环示例。for i=1:10; %i 依次取 1,2,10,.x(i=i; %对每个i值,重复执行由该指令构成的循环体，end;x%要求显示运行后数组x的值。x =1 2 3

8、 4 5 6 7 8 9 10while循环结构与For循环以固定次数求一组命令的值相反，While循环以不定的次数求一组语句的值。While循环的一般形式是：while expressi oncomma ndsend只要在表达式里的所有元素为真，就执行while和end语句之间的comma nds。【例】num=0;EPS=1;while (1+EPS1EPS=EPS/2;num=nu m+1;endnumnum =53? EPS=2*EPSEPS =2.2204e-016【例】x = zeros(1,6; % x是一个的零矩阵i = 1;while i = 6,x(i = 1/i;i =

9、i+1;end?x1.0000 0.5000 0.3333 0.2500 0.2000 0.1667可以利用 break命令跳出 while 循环 while循环可按需要嵌套。【例】Fibonacci数组的元素满足Fibonacci规则：ak2 dk ak丨，（k 12 ；且al a2 1。现要求该数组中第一个大于 10000 的元素。a(1=1;a(2=1;i=2; while a(i5 cost=(1-20/100*cost; end cost cost = 200假如有两个选择，If-Else-E nd 结 构是：6if expressi on comma nds evaluated i

10、f True else comma nds evaluated if False eSd 这里，如果表达式为真，则执行第一组命令；如果表达式是假，则执行第二组命 令。当有三个或更多的选择时，If-Else-E nd结构采用形式if expressi on1 comma nds evaluated if expressi on1 is True elseif expressi on2 comma nds evaluated if expressi on2 is True elseif expressi on3 comma nds evaluated if expressi on3 is True

11、 elseif expressi on4 comma nds evaluated if expressi on4 is Truelseif , . else comma nds evaluated ifno other expression is True enc最后的这种形式，只和所碰到的、与第一个真值表达 式相关的命令被执行；接下来的关系表达式不检验，跳过其余的If-Else-E nd结构。而且，最后的else命令可有可无。【例】if rand(10.5 disp(i love you elsedisp(i donot love you end【例】一个简单的分支结构。cost=10;nu

12、mber=12; ifnu mber8 7sums=number*0.95*cost; end, sums sums = 114.0000(例】用 for 循环指令来寻 求 Fibo nacc 数组中第一个大于 10000 的元素。n=100;a=o nes(1, n; for i=3: n a(i=a(i- 1+a(i-2; if a(i=10000 a(i, break; end end i ans = 10946 i = 21 %跳出所在的一级循 环。MATLAB的输入与输出语句输入语句o输入数值x=i nput(please in put anu mber: please in put a nu mber:22 x = 22 输入字符串 x=in put(please in put a string:1,1 please inpLLl a strin莒:this is a string x - this is a siring 输出语句 o 输出显示 命令 自由格式(disp 8disp(23+454-29*4 361 disp(11 22 33; 44 55 66; 77 88 99 11 22 33 44 55 66 77 88 99 disp(this is a siring this is a string格式化输出(fprintf。 fpri