第10章--MATLAB7-程序设计.pptVIP

例1，利用函数的递归调用，求n! n!本身就是以递归的形式定义的： 显然，求n!需要求(n-1)!，可以采用递归调用。 第10章--MATLAB7-程序设计全文共70页，当前为第64页。 递归调用函数文件factor.m: function f=factor(n) if n=1 f=1; else f=factor(n-1)*n; end 第10章--MATLAB7-程序设计全文共70页，当前为第65页。 例2，利用函数factor.m，求1!+2!+3!+4!+5!。 在命令文件main1.m： s=0; for i=1:5 s=s+factor(i); end s 第10章--MATLAB7-程序设计全文共70页，当前为第66页。 10.5.2、程序的优化 p337页 强调矩阵整体运算，尽可能避免循环运算 例 1 for i=1:20000 %用循环 x(i)=0.1*pi*i; y(i)=cos(i); end 可简化为： x=0.1*pi:0.1*pi:2000*pi;%不用循环 y=cos(x); 例2 x=0:0.02:6; y=1./((x-0.3).^2+0.01)+1./((x-0.9).^2+0.04)-6 注意：表达式中点．不能省去，否则出错。 第10章--MATLAB7-程序设计全文共70页，当前为第67页。 例3，计算下列式子的值。 t=0:3; sum(1./2.^t) 第10章--MATLAB7-程序设计全文共70页，当前为第68页。 例4, 编制三个函数,要求任意输入两个数值后,用两个函数分别求出它们的和与它们绝对值之和,再用第三个函数将这两个和相乘. 函数zha.m内容: function ch=zha(x,y) ch=zhb(x,y)*zhc(x,y); 函数zhb.m内容: function ch=zhb(x,y) ch=abs(x)+abs(y); 函数zhc.m内容： function ch=zhc(x,y) ch=x+y 缺点:函数名太多. 第10章--MATLAB7-程序设计全文共70页，当前为第69页。 子函数 P318页 例10-26, 编制一函数,要求任意输入两个数值后,用两个子函数分别求出它们的和与它们绝对值之和,再将这两个和相乘. 函数zha.m: function ch=zha(x,y) ch=zhb(x,y)*zhc(x,y); function ch=zhb(x,y) %子函数zhb ch=abs(x)+abs(y); function ch=zhc(x,y) %子函数zhc ch=x+y 第10章--MATLAB7-程序设计全文共70页，当前为第70页。 三、循环结构 循环是指按照给定的条件,重复执行指定的语句。 1、for语句 格式： for 变量=初值:增值:终值 循环体语句 end 功能：当变量的值为初值至终值时，执行循环体语句；否则for语句的后继语句。 第10章--MATLAB7-程序设计全文共70页，当前为第32页。 例1，一个3位整数各位数字的立方和等于该数本身则称该数为水仙数。输出100~999之间的全部水仙数。 设m的百位数字是m1，十位数字是m2，个位数字是m3。其中m=100,101,…,999。 数学模型： 若m=m13+m23+m33 ，则输出m。 第10章--MATLAB7-程序设计全文共70页，当前为第33页。 设m的百位数字是m1，十位数字是m2，个位数字是m3。其中m=100,101,…,999。 数学模型：若m=m13+m23+m33 ，则输出m。 命令文件TT.m的内容： for m=100:999 m1=fix(m/100);%求m的百位数字 m2=rem(fix(m/10),10);%求m的十位数字 m3=rem(m,10);%求m的个位数字 if m==m1^3+m2^3+m3^3 disp(m) end end 第10章--MATLAB7-程序设计全文共70页，当前为第34页。 命令文件TT.m的内容： for m=100:999 %第1行 m1=fix(m/100);%求m的百位数字 %第2行 m2=rem(fix(m/10),10);%求m的十位数字 %第3行 m3=rem(m,10);%求m的个位数字