第4章 MATLAB程序设计精要.pptVIP

*/47 例4.13 求[100，200]之间第一个能被21整除的整数。 程序如下： for n=100:200 if rem(n,21)~=0 continue end break end n */47 4．循环的嵌套 如果一个循环结构的循环体又包括一个循环结构，就称为循环的嵌套，或称为多重循环结构。 例4.14 若一个数等于它的各个真因子之和，则称该数为完数，如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 note：for m=1:500 x=1:fix(m/2); if (rem(m,x)==0)*x==m disp(m) end end */47 例4.15 用筛选法求某自然数范围内的全部素数。 程序如下： m=input(m=); p=1:m; p(1)=0; for i=2:sqrt(m) for j=2*i:i:m p(j)=0; end end n=find(p~=0); %nonzeros(p) 或者 p(p~=0) p(n) 用筛选法求素数的基本思想是： 要找出2~m之间的全部素数， 首先在2~m中划去2的倍数(不包括2)， 然后划去3的倍数(不包括3)， 由于4已被划去， 再找5的倍数 (不包括5)，…， 直到再划去不超过的数的倍数， 剩下的数都是素数。 Note：primes(m) m以内的素数表 isprime(x) 判断x是否为素数 */47 4.3 函数文件 4.3.1 函数文件的基本结构 函数文件由function语句引导，其基本结构为： function 输出形参表=函数名(输入形参表) 注释说明部分 函数体语句 其中以function开头的一行为引导行，表示该M文件是一个函数文件。函数名的命名规则与变量名相同。输入形参为函数的输入参数，输出形参为函数的输出参数。当输出形参多于一个时，则应该用方括号括起来。 Note:函数文件名？ 注释部分的约定？return语句？局部变量？ */47 例4.16 编写函数文件求半径为r的圆的面积和周长。 函数文件如下： function [s,p]=fcircle(r) %FCIRCLE calculate the area and perimeter of a circle of radii r %r 圆半径 %s 圆面积 %p 圆周长 %2013年4月5日编 s=pi*r*r; p=2*pi*r; help fcircle lookfor Note:函数体中一般不用input()和disp(). */47 4.3.2 函数调用 函数调用的一般格式是： [输出实参表]=函数名(输入实参表) 要注意的是，函数调用时各实参出现的顺序、个数，应与函数定义时形参的顺序、个数一致，否则会出错。函数调用时，先将实参传递给相应的形参，从而实现参数传递，然后再执行函数的功能。 Note：传值拷贝，不是传地址。 系统给出的函数往往是矩阵函数，而用户定义的往往是标量（数量）函数。 */47 例4.17 利用函数文件，实现直角坐标(x,y)与极坐标(ρ,θ)之间的转换。 函数文件tran.m： function [rho,theta]=tran(x,y) %cart2pol ,pol2cart 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 在MATLAB中，函数可以嵌套调用，即一个函数可以调用别的函数，甚至调用它自身。一个函数调用它自身称为函数的递归调用。 */47 例4.18 利用函数的递归调用，求n！。 n!本身就是以递归的形式定义的： 显然，求n!需要求(n-1)!，这时可采用递归调用。递归调用函数文件factor.m如下： function f=factor(n)%最好不要取此函数名 i