计算机模拟第10课.pptVIP

4.3 函数文件 4.3.1 函数文件的基本结构 函数文件由function语句引导，其基本结构为： function [输出形参表]=函数名(输入形参表) 注释说明部分 函数体语句 其中以function开头的一行为函数声明行，表示该M文件是一个函数文件。函数名的命名规则与变量名相同。输入形参为函数的输入参数，输出形参为函数的输出参数。当输出形参多于一个时，则应该用方括号括起来。 说明： （1）关于函数文件名。函数文件名必须与函数名相同（若不同，调用时使用函数文件名）。 （2）关于注释说明部分。（包括3部分内容） 以%开头的第一注释行。这行一般括大写的函数名和函数功能简要描述。 第一注释行之后连续的注释行。通常包括输入输出参数的含义、调用格式说明等信息，构成全部在线帮助文本。 与在线帮助文本相隔一空行的注释行。包括函数文件编写和修改的信息。 例4-14 编写函数文件求半径为r的圆的面积和周长。 函数文件如下： function [s,p]=fcircle(r) %CIRCLE calculate the area and perimeter of a circle of radii r %r 圆半径 %s 圆面积 %p 圆周长 %2009年10月15日编 s=pi*r*r; p=2*pi*r; 4.3.2 函数调用 函数调用的一般格式是： [输出实参表]=函数名(输入实参表) 注意：函数调用时各实参出现的顺序、个数，应与函数定义时形参的顺序、个数一致，否则会出错。 函数调用时，先将实参传递给相应的形参，从而实现参数传递，然后再执行函数的功能。 例4-15 利用函数文件，实现直角坐标(x,y)与极坐标(ρ,θ)之间的转换。 函数文件tran.m： function [rho,theta]=tran(x,y) rho=sqrt(x*x+y*y); theta=atan(y/x); 调用tran.m的命令文件main1.m： x1=input(Please input x=:); y1=input(Please input y=:); [rho,the]=tran(x1,y1); rho the MATLAB中，函数可以嵌套调用，即一个函数可以调用别的函数，甚至调用它自身。一个函数调用它自身称为函数的递归调用。 例4-16 利用函数的递归调用，求n！。 递归调用函数文件factor.m如下： function f=factor(n) if n=1 f=1; else f=factor(n-1)*n; %递归调用求(n-1)! end 4.3.3 函数参数的可调性 在调用函数时，MATLAB用两个永久变量 nargin 和 nargout 分别记录调用该函数时的输入实参和输出实参的个数。只要在函数文件中包含这两个变量，就可以准确地知道该函数文件被调用时的输入输出参数个数，从而决定函数如何进行处理。 例4-17 nargin用法示例。 函数文件examp.m： function fout=charray(a,b,c) if nargin==1 fout=a; elseif nargin==2 fout=a+b; elseif nargin==3 fout=(a*b*c)/2; end 命令文件mdt.m： x=[1:3]; y=[1;2;3]; examp(x) examp(x,y) examp(x,y,3) 4.3.4 全局变量与局部变量 函数文件中的变量是局部的，与其他函数文件及matlab基本工作空间相互隔离。 如果在若干函数中，都把某一变量定义为全局变量，那么这些函数将公用这一变量。 全局变量用global命令定义，格式为： global 变量名表 例4-18 全局变量应用示例。 先建立函数文件wadd.m，该函数将输入的参数加权相加。 function f=wadd(x,y) global ALPHA BETA f=ALPHA*x+BETA*y; 在命令窗口中输入： global ALPHA BETA ALPHA=1; BETA=2; s=wadd(1,2) 4.3.5 主函数与子函数 M文件可以包含任意多个函数，主函数是指M文件中的 第一个函数。子函数是除了主函数之外的函数。 一个M文件可以有任意多个子函数，但只能有一个主函数。 例4-19 主