海南大学《MATLAB基础与应用》课件-第3章程序设计.pptVIP

例3-8 一个三位整数各位数字的立方和等于该数本身则称该数为水仙花数。输出全部水仙花数。 程序如下： 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 例3-9：使用for结构计算1+2+3+…+100 程序如下： clear sum=0; for i=1:100 sum=sum+i; end sum sum = 5050 例3-9 已知 , i为自然数. 当n=100时，求y的值。 程序如下： y=0; n=100; for k=1:n y=y+1/(2*k-1); end y 3.2.3 循环结构(续) 在实际MATLAB编程中，采用循环语句通常会降低其执行速度，所以前面的程序可由下面的程序来代替： n=100; k=1:2:2*n-1; y=sum(1./k); y 3.2.3 循环结构(续) for语句更一般的格式为： for 循环变量=矩阵表达式 循环体语句 end 执行过程是依次将矩阵的各列元素赋给循环变量，然后执行循环体语句，直至各列元素处理完毕。 例3-10 写出下列程序的执行结果。 s=0;k=0; a=[12,13,14;15,16,17;18,19,20;21,22,23]; for k=a s=s+k; end disp(s); 3.2.3 循环结构(续) 2．while语句 while语句的一般格式为： while (条件) 循环体语句 end 其执行过程为：若条件成立，则执行循环体语句，执行后再判断条件是否成立，如果不成立则跳出循环。 例：使用while结构计算1+2+3+…+100 clear sum=0; k=0; while k100 k=k+1; sum=sum+k; end sum 运行结果： sum = 5050 例3-11 从键盘输入若干个数，当输入0时结束输入，求这些数的平均值和它们之和。 程序如下： sum=0; cnt=0; val=input(Enter a number (end in 0):); while (val~=0) sum=sum+val; cnt=cnt+1; val=input(Enter a number (end in 0):); end if (cnt 0) sum mean=sum/cnt end 3.2.3 循环结构(续) 3．break语句和continue语句 与循环结构相关的语句还有break语句和continue语句。它们一般与if语句配合使用。 break语句用于终止循环的执行。当在循环体内执行到该语句时，程序将跳出循环，继续执行循环语句的下一语句。 continue语句控制跳过循环体中的某些语句。当在循环体内执行到该语句时，程序将跳过循环体中所有剩下的语句，继续下一次循环。 例3-12 求[100，200]之间第一个能被21整除的整数。 程序如下： for n=100:200 if rem(n,21)~=0 continue end break end n 3.2.3 循环结构(续) 4．循环的嵌套 如果一个循环结构的循环体又包括一个循环结构，就称为循环的嵌套，或称为多重循环结构。 例3-13 若一个数等于它的各个真因子之和，则称该数为完数，如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 3.3 函数文件 3.3.1 函数文件的基本结构 函数文件由function语句引导，其基本结构为： function 输出形参表=函数名(输入形参表) 注释说明部分 函数体语句 其中以function开头的一行为引导行，表示该M文件是一个函数文件。函数名的命名规则与变量名相同。输入形参为函数