第5讲MATLAB程序设计讲述.pptVIP

例13 编写函数文件求半径为r的圆的面积和周长。 函数文件如下： function [s,p]=fcircle(r) %CIRCLE calculate the area and perimeter of a circle of radii r %r 圆半径 %s 圆面积 %p 圆周长 ? % auther:aaa s=pi*r*r; p=2*pi*r; 将以上函数文件以文件名fcircle.m存入c:\matlab??\work下，然后在MATLAB命令窗口调用该函数： [s,p]=fcircle(10) 6.5.2 函数调用 函数调用的一般格式是： [输出实参表]=函数名(输入实参表) 例14 利用函数文件，实现直角坐标(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 例15 利用函数的递归调用，求n！。 递归调用函数文件factor.m： function f=factor(n) if n=1 f=1; else f=factor(n-1)*n; end 6.5.3 函数所传递参数的可调性 在调用函数时，MATLAB用两个永久变量nargin和nargout分别记录调用该函数时的输入实参和输出实参的个数。只要在函数文件中包含这两个变量，就可以准确地知道该函数文件被调用时的输入输出参数个数，从而决定函数如何进行处理。 用于功能较复杂的函数。 例16 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) 6.6 全局变量和局部变量 全局变量用命令global定义。函数文件的内部变量是局部的，不加载到工作空间中，与其他函数文件及MATLAB工作空间相互隔离； 定义全局变量时，变量之间用空格分隔。 例17 全局变量应用示例 先建立函数文件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) 6.7 文件操作 6.7.1 文件的打开与关闭 ? 1. 打开文件： fopen函数的调用格式为： fid= fopen(‘文件名’，‘打开方式’，‘读写数据格式’) fid：文件句柄。一般是正整数，fid=-1表示文件打开不成功。 2. 关闭文件 调用格式为： sta=fclose(fid) 打开方式： 默认情况：二进制文件格式，可添加‘b’。 如果是文本文件，加上‘t’。如‘rt’。 数据格式： ‘native’ or ‘n’??????系统默认 ieee-le or l???? IEEE 浮点数（ little-endian byte ordering） ieee-be or b???? IEEE 浮点数（ big-endian byte ordering） 6.7.2 二进制文件读写操作 1. 读二进制文件 fread 函数可以读取二进制文件的