第4讲-Matlab程序设计.pptVIP

2、程序控制结构 for n=100:200 if rem(n,21)~=0 continue end break end n %输出结果105 例 求[100，200]之间第一个能被21整除的整数。 第4讲-Matlab程序设计全文共58页，当前为第32页。 2、程序控制结构 如果一个循环结构的循环体又包括一个循环结构，就称为循环的嵌套，或称为多重循环结构。 （4） 循环的嵌套 例 若一个数等于它的各个真因子之和，则称该数为完数，如6=1+2+3，所以6是完数。求[1,500]之间的全部完数。 for m=1:500 s=0; for k=1:m/2 if rem(m,k)==0 s=s+k; end end if m==s disp(m); end end 结果为6，28，496 第4讲-Matlab程序设计全文共58页，当前为第33页。 3、函数文件 function 输出形参表=函数名(输入形参表) 注释说明部分 函数体语句 （1）函数文件的基本结构 以function开头的一行为引导行，表示该M文件是一个函数文件。 函数名的命名规则与变量名相同。 输入形参为函数的输入参数，输出形参为函数的输出参数。 当输出形参多于一个时，则应该用方括号括起来。 函数文件名通常由函数名加上扩展名.m组成，函数名要与函数文件名一致； 第4讲-Matlab程序设计全文共58页，当前为第34页。 3、函数文件 function 输出形参表=函数名(输入形参表) 注释说明部分 函数体语句 （1）函数文件的基本结构 紧随函数文件引导行之后以%开头的第一注释行，一般包括大写的函数文件名和函数功能简要说明，供lookfor关键词查询和help在线帮助时使用； 第一注释行及之后连续的注释行，通常包括函数输入、输出参数的含义及调用格式说明等信息，构成全部在线帮助文本； 与在线帮助文本相隔一空行的注释行，包括函数文件编写和修改的信息，如作者、修改日期、版本等内容，用于软件档案管理。 第4讲-Matlab程序设计全文共58页，当前为第35页。 3、函数文件 function 输出形参表=函数名(输入形参表) 注释说明部分 函数体语句 （1）函数文件的基本结构 如果在函数中插入return语句，则执行到该语句就结束函数的运行，程序流程转到调用该函数的位置； 在函数文件中，也可以不使用return语句，这时被调用的函数执行完成后自动返回。 第4讲-Matlab程序设计全文共58页，当前为第36页。 3、函数文件 例 编写函数文件，求半径为r的圆的面积和周长。 function [s,p]=fcircle(r) %FCIRCLE calculate the area and perimeter of a circle of radii r %r 圆半径 %s 圆面积 %p 圆周长 %2004年7月30日编 s=pi*r*r; p=2*pi*r; 执行下面语句，观察结果 [s,p]=fcircle(10) help fcircle lookfor fcircle 第4讲-Matlab程序设计全文共58页，当前为第37页。 3、函数文件 （2）函数调用 [输出实参表]=函数名(输入实参表) 要注意的是，函数调用时各实参出现的顺序、个数，应与函数定义时形参的顺序、个数一致，否则会出错。 函数调用时，先将实参传递给相应的形参，从而实现参数传递，然后再执行函数的功能。 第4讲-Matlab程序设计全文共58页，当前为第38页。 3、函数文件 例利用函数文件，实现直角坐标(x,y)与极坐标(ρ,θ)之间的转换。 函数文件tran.m： function [rho,theta]=tran(x,y) rho=sqrt(x*x+y*y); theta=atan(y/x); 调用tran.m的命令文件main1.m： x=input(Please input x=:); y=input(Please input y=:); [rho,the]=tran(x,y); rho the 第4讲-Matlab程序设计全文共58页，当前为第39页。 3、函数文件 例 利用函数的递归调用，求n！。 n!本身就是以递归的形式定义的： 显然，求n!需要求(n-1)!，这时可采用递归调用。 function f=factor(n) if