5-MATLAB程序设计内容资料.ppt

5.3.2 函数调用 函数调用的一般格式是： [输出实参表]=函数名(输入实参表) 注：函数调用时各实参出现的顺序、个数，应与函数定义时形参的顺序、个数一致，否则会出错。函数调用时，先将实参传递给相应的形参，从而实现参数传递，然后再执行函数的功能。 例5-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： x=input(Please input x=:); y=input(Please input y=:); [rho,the]=tran(x,y); rho the 在MATLAB中，函数可以嵌套调用，即一个函数可 以调用别的函数，甚至调用它自身。一个函数调用 它自身称为函数的递归调用。 例5-16 利用函数的递归调用，求n！。 递归调用函数文件factor.m如下： function f=factor(n) if n=1 f=1; else f=factor(n-1)*n; %递归调用求(n-1)! end 5.3.3 函数参数的可调性 两个永久变量： nargin：记录调用该函数时的输入实参个数。 nargout：记录调用该函数时的输出实参个数。 只要在函数文件中包含这两个变量，就可以准确地 知道该函数文件被调用时的输入输出参数个数，从 而决定函数如何进行处理。 例5-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 命令文件mydemo.m： x=[1:3]; y=[1;2;3]; examp(x) examp(x,y) examp(x,y,3) 5.3.4 全局变量与局部变量 全局变量用global命令定义，格式：global 变量 例5-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) 5.4 程序举例 例5-19 猜数游戏。 首先由计算机产生[1,100]之间的随机整数，然 后由用户猜测所产生的随机数。根据用户猜测 的情况给出不同提示，如猜测的数大于产生的 数，则显示“High”，小于则显示 “Low”，等 于则显示“You won”，同时退出游戏。用户最 多可以猜7次。 例5-20 用筛选法求某自然数范围内的全部素数。 素数是大于1，且除了1和它本身以外，不能被其他任何整数所整除的整数。用筛选法求素数的基本思想是：要找出2~m之间的全部素数，首先在2~m中划去2的倍数(不包括2)，然后划去3的倍数(不包括3)，由于4已被划去，再找5的倍数 (不包括5)，…，直到再划去不超过的数的倍数，剩下的数都是素数。 n=0; for m=100:1000 flag=1; j=m-1; i=2; while i=j flag if rem(m,i)==0 flag=0; end i=i+1; end if flag n=n+1; prime(n)=m; end end prime 【例】 求[100,1000]以内的全部素数。 %变量prime存放素数 例5-21 Fibonacci数列定义如下： f1=1 f2=1 fn=fn-1+fn-2 (n2) 编写函数文件求小于任意自然数n的Fibonacci数列各项。 function f=ffib(n) %用于求Fibonacci数列的函数文件 %f=ffib(n) %2008年5月20日编 f=[1,1]; i=1; while f(i)+f(i+1)n f(i