[工程科技]MATLAB软件第4章.ppt

例11 求1+2+3+……+100的和 i=0; s=0; for i=1:100 s=s+i; end s 例12 一个三位整数各位数字的立方和等于该数本身则称该数为水仙花数。输出全部水仙花数。 程序如下： 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 例13 已知 ， 当n=100时，求y的值。 程序如下： y=0; n=100; for i=1:n y=y+1/(2*i-1); end y 在实际MATLAB编程中，采用循环语句会降低其执行速度，所以前面的程序通常由下面的程序来代替： n=100; i=1:2:2*n-1; y=sum(1./i); y for语句更一般的格式为： for 循环变量=矩阵表达式 循环体语句 end 执行过程是依次将矩阵的各列元素赋给循环变量，然后执行循环体语句，直至各列元素处理完毕。 例14 写出下列程序的执行结果。 s=0; a=[12,13,14;15,16,17;18,19,20;21,22,23]; for k=a s=s+k; end disp(s); 3．break语句和continue语句 与循环结构相关的语句还有break语句和continue语句。它们一般与if语句配合使用。 break语句用于终止循环的执行。当在循环体内执行到该语句时，程序将跳出循环，继续执行循环语句的下一语句。 continue语句控制跳过循环体中的某些语句。当在循环体内执行到该语句时，程序将跳过循环体中所有剩下的语句，继续下一次循环。 例15 求[100，200]之间第一个能被21整除的整数。 程序如下： for n=100:200 if rem(n,21)~=0 continue end break end n 4．循环的嵌套 如果一个循环结构的循环体又包括一个循环结构，就称为循环的嵌套，或称为多重循环结构。 例16 有一数: 11+12+13…+110+21+22+23…+210+31+32+33…+310 s=0; for a=1:3 for b=1:10 s=s+a^b; end end s 编写程序：猜数游戏 首先由计算机产生[1,100]之间的随机整数，然后由用户猜测所产生的随机数。根据用户猜测的情况给出不同提示，如猜测的数大于产生的数，则显示“高了”，小于则显示“低了”，等于则显示“你赢了”，同时退出游戏。用户最多可以猜7次。若 7次猜不中，则显示“你输了” 4.3 函数文件 4.3.1 函数文件的基本结构 函数文件由function语句引导，其基本结构为： function 输出形参表=函数名(输入形参表) 注释说明部分 函数体语句 其中以function开头的一行为引导行，表示该M文件是一个函数文件。函数名的命名规则与变量名相同。输入形参为函数的输入参数，输出形参为函数的输出参数。当输出形参多于一个时，则应该用方括号括起来。 1.函数申明行 ：位于函数文件的首行，以关键字 function开 头，函数名以及函数的输入输出参量都在这一行定义。函数名 和文件名必须相同。 2.注释行：紧随函数申明行之后以 %开头的第一行注释行，包括大写体的函数文件名；运用关键词要描述的函数功能。 3.函数体：这部分内容由实现该 M函数文件功能的Matlab指令 组成。它接收输入参量，进行程序流控制，得到输出参量。 例17 编写函数文件求半径为r的圆的面积和周长。 函数文件如下： function [s,p]=fcircle(r) %CIRCLE calculate the area and perimeter of a circle of radii r %r 圆半径 %s 圆面积 %p 圆周长 s=pi*r*r; p=2*pi*r; 4.3.2 函数调用 函数调用的一般格式是： [输出实参表]=函数名(输入实参表) 要注意的是，函数调用时各实