MATLAB8.X程序的设计及典型应用第三章 M文件初步.pptVIP

【例3-14】global指令定义全局变量示例。编写计算圆周周长和面积的函数文件ffun_SC: function [C,S]=ffun_SC %FFUN_SC Calculate the area and the the perimeter of a circle global R C=2*pi*R; S=pi*R^2; end 编写脚本文件并命名为exm3_14: clear global R R=4; [C,S]=ffun_SC 在指令窗中执行指令： exm3_14 C = 25.1327 S = 50.2655 3.4 MATLAB函数类别和句柄函数 函数文件是扩展名为.m的M文件中的一种，而函数文件又可以被细分为主函数、子函数、私用函数、匿名函数和内联函数等。本文只介绍主函数、子函数和内联函数，并简要介绍函数句柄。 3.4.1 主函数 MATLAB允许一个函数文件包含多个函数。其中第一个出现的函数被称为主函数(Primary function)，该文件中的其他函数都被称为子函数(Subfunction)。保存时所用函数文件名与主函数定义名相同。外部程序只能对主函数进行调用。 主函数的特点为： 一般为“与保存文件名相同”的那个函数。 在当前目录、搜索路径上，列出文件名的函数。 在指令窗中或其他函数中，可以直接被调用的函数。 M函数文件中，有第一个function引导的函数。 采用help functionname可获取函数所携带的帮助信息。 3.2.4 for循环结构 for语句是一种基本的实现循环结构的语句，能够以确定的次数执行某一段程序。for语句的格式如下： 　　for 循环变量=表达式 　　 循环体指令语句组 　　end 【说明】 “表达式”可以是MATLAB指令产生的数组，也可以是任意给定的一个数组。循环变量从“表达式”中的第一个数值（或第一列数组）一直循环到“表达式”的最后一个数值（或最后一列数组）。 【例3-7】运用for…end循环结构计算 的值，（其中m=100） 运用for…end循环结构计算，编写文件名为exm3_7的脚本文件： clear y=0; for n=0:100 y=y+1/(2^n+1); end Y 在指令窗中执行文件exm3_7.m，运行结果为： y = 1.2645 a=0+eps;%利用极小数代替0，避免计算时出错 b=3;n=1000;h=(b-a)/n; x=a:h:b;y=0; f=log(x)./(2*x.^2+1); for i=1:n s(i)=(f(i)+f(i+1))*h/2; %计算细分后梯形的面积 y=y+s(i); %迭代法 end y 【例3-8】运用循环结构计算积分 在指令窗中执行该文件，运行结果为： y = -0.7730 【例3-8】一个三位整数各位数字的立方和等于该数本身，则称该数为水仙花数。运用for…end循环结构输出全部水仙花数。 clear for m=[100:999] m1=fix(m/100); %求m的百位数字,函数fix()功能为截断取整 m2=rem(fix(m/10),10); %求m的十位数字,rem()为求余函数 m3=rem(m,10); %求m的个位数字 if m==m1*m1*m1+m2*m2*m2+m3*m3*m3 disp(m) end end 运行结果为： 153 370 371 407 3.2.5 while循环结构 　　while结构是MATLAB语言实现循环结构的另一种基本方式，以不定的次数重复执行某一段程序。 　　while循环结构的一般形式如下： 　　while 逻辑表达式 　　　　指令语句组 　　end 执行时，只要逻辑表达式为“逻辑真”（非0），就执行指令语句组，执行后再返回到while引导的逻辑表达式处，继续判断；如果逻辑表达式为“逻辑假”（0），则跳出循环。 通常，逻辑表达式的值为一个标量，但数组也同样有效。如果逻辑表达式的值为数组，要求数组的所有元素都为“逻辑真”（非0），MATLAB才会执行循环体中的指令语句组。如果逻辑表达式为“空”，MATLAB则认为表达式的值为逻辑假，不执行循环体语句指令组。 clear val=input(Enter a number (end in 0):); sum=0;cnt=0; while (val~=0) sum=sum+val; cnt=cnt+1; val=input(Enter a number (end in 0):)