第4章 MATLAB程序设计13685.pptVIP

循环结构 for语句 for语句的格式为： for 循环变量=表达式1:表达式2:表达式3循环体语句end 其中表达式1的值为循环变量的初值，表达式2的值为步长，表达式3的值为循环变量的终值。步长为1时，表达式2可以省略。 例4.7 一个三位整数各位数字的立方和等于该数本身则称该数为水仙花数。输出全部水仙花数。 程序如下： for m=100:999 m1=fix(m/100);%求m的百位数字 m2=rem(fix(m/10),10); %求m的十位数字 m3=rem(m,10);%求m的个位数字 if m==m1*m1*m1+m2*m2*m2+m3*m3*m3 disp(m) end end 例4.8 已知y，当n=100时，求y的值。 程序如下： y=0;n=100; for i=1:n y=y+1/i/i; end y 例4.9 求定积分。 以梯形法为例，程序如下： a=0;b=3*pi; n=1000; h=(b-a)/n; x=a; s=0; f0=exp(-0.5*x)*sin(x+pi/6); for i=1:nx=x+h;f1=exp(-0.5*x)*sin(x+pi/6);s=s+(f0+f1)*h/2;f0=f1; end s 在实际MATLAB编程中，采用循环语句会降低其执行速度，所以前面的程序通常由下面的程序来代替： n=100; i=1:n; f=1./i.^2; y=sum(f) for语句更一般的格式为：for 循环变量=矩阵表达式循环体语句end执行过程是依次将矩阵的各列元素赋给循环变量，然后执行循环体语句，直至各列元素处理完毕。 例4.10 写出下列程序的执行结果。s=0; a=[12,13,14;15,16,17;18,19,20;21,22,23];for k=as=s+k;enddisp(s); while语句while语句的一般格式为：while (条件)循环体语句end其执行过程为：若条件成立，则执行循环体语句，执行后再判断条件是否成立，如果不成立则跳出循环。 例4.11 从键盘输入若干个数，当输入0时结束输入，求这些数的平均值和它们之和。 程序如下： sum=0; n=0; val=input(Enter a number (end in 0):); while (val~=0)sum=sum+val;n=n+1;val=input(Enter a number (end in 0):); end if (n 0)summean=sum/n end 例4.12 根据矩阵指数的幂级数展开式求矩阵指数。 程序如下： X=input(Enter X:); E=zeros(size(X)); F=eye(size(X)); n=1; while norm(F,1)0E=E+F;F=F*X/n;n=n+1; end E expm(X) break语句和continue语句 与循环结构相关的语句还有break语句和continue语句。它们一般与if语句配合使用。 break语句用于终止循环的执行。当在循环体内执行到该语句时，程序将跳出循环，继续执行循环语句的下一语句。 continue语句控制跳过循环体中的某些语句。当在循环体内执行到该语句时，程序将跳过循环体中所有剩下的语句，继续下一次循环。 例4.13 求[100，200]之间第一个能被21整除的整数。 程序如下： for n=100:200 if rem(n,21)~=0continue end break end n 循环的嵌套 如果一个循环结构的循环体又包括一个循环结构，就称为循环的嵌套，或称为多重循环结构。 例4.14 若一个数等于它的各个真因子之和，则称该数为完数，如6=1+2+3，所以6是完数。求[1,500]之间的全部完数。 for m=1:500 s=0; for k=1:m/2 if rem(m,k)==0s=s+k; end end if m==sdisp(m); end end 例4.15 用筛选法求某自然数范围内的全部素数。 程序如下： m=input(m=); p=1:m; p(1)=0; for i=2:sqrt(m)for j=2*i:i:mp(j)=0;end end n=find(p~=0); p(n) 函数文件 function 输出形参表 = 函数名(输入行参表) 注释说明部分（可选） 函数体语句(必须) 其中以 function开始的一行为引导行，表示该 M文件是函数文件。 注释语句段的