第四章MATLAB程序设计.pptVIP

0.5 sin( ) 6 x y e x ? ? ? ? 例 4.11 求 定积分。 实例分析 a=0;b=3*pi; n=1000 ; h=(b-a)/n; x=a:h:b; f=exp(-0.5*x).*sin(x+pi/6); for i=1:n s(i)= (f(i)+f(i+1))*h/2; end s=sum(s) 可直接调用函数 trapz() ，求数值积分！ 实例分析 for 语句更一般的格式： for 循环变量 = 矩阵表达式 循环体语句 end 执行过程是依次将矩阵的 各列元素 赋给循环变量，然后 执行循环体语句，直至各列元素处理完毕。实际上， “表达式 1: 表达式 2: 表达式 3” 是一个仅为一行的矩阵 ( 行 向量 ) ，因而 列向量 是单个数据。 第三节 循环结构 例 4.12 已知 5 个学生 4 门功课的成绩，求每名学生的总成绩。 s=0; a=[65,76,56,78;98,83,74,85; ... 76,67,78,79;98,58,42,73; ... 67,89,76,87]; for k=a s=s+k; end disp(s); 实例分析 3.2 while 语句 while 语句的一般格式为： while ( 条件 ) 循环体语句 end 其执行过程为：若条件成立，则执行循环体语句，执行后 再判断条件是否成立，如果不成立则跳出循环体，向下继 续执行。 第三节 循环结构 while 循环和 for 循环的区别： while 循环结构的循环体被执行的 次数不是确定的，而 for 结构中循环体的执行次数是确定的。 例 4.13 求 a=1+2+…+100 ， b=1/i (i=1~10) sum=0; i=0; while(i=100) sum=sum+i; i=i+1; end sum sum = 5050 for I=1:10 A(I)=1/(I+1) ; end A 实例分析 其它与流程控制有关的语句 break 语句 循环结构 continue 语句 循环结构 第三节 循环结构 第三节 循环结构 break 语句 ：导致包含该指令的循环强制终止，即跳 出所在的一级循环。该语句一般与 if 语句配合使用 continue 语句 ：跳过位于其后的其它指令，执行循环 的下一次迭代循环 y=0; i=1; while 1 % 逻辑真 f=1/i/i; y=y+f; if i==100 break; end i=i+1; end y 例 4.14 用 while 语句实现求 a=1+1/2 2 +…+1/100 2 向量运算来代替循环操作 n=100; i=1:n; f=1./i.^2; y=sum(f) %sum() 求和函数 实例分析 例 4.15 求 [100, 200] 之间第一个能被 21 整除的整数。 for n=100:200 if rem(n,21)~=0 continue end break end n 实例分析 3.3 循环的嵌套 ? 如果一个循环结构的循环体又包括一个循环结构， 称为 循环的嵌套 （多重循环结构） ? 处于内部的循环叫作 内循环 ； 处于外部的循环叫作 外循环 。 ? 在设计多重循环时，要特别注意内、外循环之间的 关系，以及各语句放置的位置。 第三节 循环结构 例 4.16 若一个数等于它的各个真因子之和，则称该数为完数， 如 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 实例分析 m=input(m=); p=2:m; for i=2:sqrt(m) n=find(rem(p,i)==0p~=i); p(n)=[]; end p 实例分析 例 4.17 用筛选法求某自然数范围内的全部 素数 。 一个数，如果只 有 1 和它本身两 个因数，这样的 数叫做 质数 ，又 称 素数 第四节 函数文件 function 输出形参表 = 函数名 ( 输入形参表 ) 注释说明部分 函数体语句 4.1 函数文件的基本结构 注意：当形参表多于 1 个时，需要中括号 [] 如： [m,n]=size(A) 说 明： (1) 函数文件名。当函数文件名与函数名不同时， MATLAB 将忽略函数名而确认函数文件名，因此调用时使用函数文 件名； (2) 注释说明部分。注释说明包括三部分内容： ①函数文件引导行之后以 % 开头的第一注释行； ②第一注释行及之后连续的注释行； ③与在线帮助文本相隔一空行的注释行。 (3) r