五Matlab编程入门二.ppt

1、数学实验,Matlab 编程 (二),function 输出形参列表=函数名(输入形参列表) % 注释说明部分(可选) 函数体语句(必须),第一行为引导行，表示该 M文件是函数文件 函数名的命名规则与变量名相同 ( 必须以字母开头 ) 当输出行参多于一个时，用方括号括起来 函数必须是一个单独的 M文件 函数文件名必须与函数名一致 以百分号开始的语句为注释语句,函数文件,函数文件由 function 语句引导,function y = myfliplr(x) %myfliplr Flip matrix in left/right direction. % returns X with row p

2、reserved and columns % flipped in the left/right direction % % X = 1 2 3 becomes 3 2 1 % 4 5 6 6 5 4 % num_col=size(x,2); % 取 x 的列数 y = x(:,num_col:-1:1);,函数文件举例,myfliplr.m,function printyh(n) % 打印杨辉三角形, 本函数没有输出参数 yh=1; disp(yh); if n=1, return; end yh=1,1; disp(yh); for k=3:n yh_old=yh; k2=ceil(k/2

3、); for l=2:k2 yh(l)=yh_old(l-1)+yh_old(l); end yh(k2+1:k)=yh(k-k2:-1:1); disp(yh); end,函数文件举例,printyh.m,例：编写函数文件求半径为r的圆的面积和周长。 函数文件如下： function s,p=fcircle(r) %CIRCLE calculate the area and perimeter of a circle of radii r %r 圆半径 %s 圆面积 %p 圆周长 s=pi*r*r; p=2*pi*r;,函数申明行,帮助文本区,函数体,输入参数,输出参数,程序2： 建立一个函

4、数，将华氏温度转变为摄氏温度，已知公式为： c=5*(f-32)/9 首先建立函数文件f2c.m。 function c=f2c(f) 函数定义行 % compute the Fahrenheit temperature to Celsius temperature c=5*(f-32)/9 函数主体 然后在MATLAB的命令窗口调用该函数文件。 clear; y=input(Input Fahrenheit temperature：); x=f2c(y) 输出情况为： Input Fahrenheit temperature：70 c = 21.1111 x = 21.1111,利用函数文件

5、，实现直角坐标(x,y)与极坐标(,)之间的转换。【数学基础=sqrt(x*x+y*y),=atan(y/x)】 函数文件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,函数调用的一般格式,输出实参列表=函数名(输入实参列表),函数调用时，实参的顺序应与函数定义时的形参的顺序一致。 函数可以嵌套

6、调用，即一个函数可以被其它函数调用，甚至可以被它自身调用，此时称为递归调用。 函数所传递的参数具有可调性，Matlab 用两个永久变量 nargin 和 nargout 分别记录调用该函数时的输入实参和输出实参的个数。,函数调用,% 函数文件 factor.m function f=factor(n) if (n=1) f=1; else f=n*factor(n-1); end,% main.m % 计算 s=1!+2!+3!+4!+5! % s=0; for i=1:5 s=s+factor(i); end fprintf( s=%g n,s),例：利用函数的递归调用计算 n!,递归函数举

7、例,% ex4argin.m function out=ex4argin(a,b) if (nargin=1) out=a; elseif (nargin=2) out=a*b; end,参数的可调性举例,% ex4argout.m function p,q=ex4argout(a,b) if (nargout=1) p=a+b; elseif (nargout=2) p=a+b; q=a-b; end,函数通过输入和输出参数与其它 M 文件进行数据传递。,当函数调用完毕后，该函数文件中定义的所有局部变量都将被释放，即全部被清除。,如果在若干个 M文件中，都把某个变量定义为全局变量，则这些函数

8、将公共使用这一变量。所有函数都可以对它进行存取和修改操作。,函数文件中的变量都是局部的，即一个函数文件中定义的变量不能被另一个函数文件或其它 M 文件使用。,全局变量与局部变量,定义全局变量是 M文件间传递信息的一种手段。,变量名列表中的各个变量用空格隔开，不能用逗号！,全局变量给函数间的数据传递带来了方便，但却破坏了函数对变量的封装，降低了程序的可读性，因而在结构化程序设计中，全局变量是不受欢迎的。特别是当程序较大，子程序较多时，全局变量将个程序调试和维护带来不便，故不提倡使用全局变量。,global 变量名列表,全局变量的定义,程序示例,clear; global a b a=1; b=3

9、; y=mysquaresum(a,b); fprintf( a=%g, b=%g n,a,b); z=myproduct(a,b); fprintf( a=%g, b=%gn,a,b);,function square_sum=mysquaresum(x,y) square_sum=x2 + y2;,function product=myproduct(x,y) global a product=x*y; a=x+y;,子函数,function avg, med = ex4subfun(x) % 主函数 n = length(x); avg = mean(x, n); med = media

10、n(x, n); function a = mean(x, n) % 子函数 a = sum(x)/n; function m = median(x, n) % 子函数 x = sort(x); if rem(n, 2) = 1 m = x(n+1)/2); else m = (x(n/2)+x(n/2+1)/2; end,子函数,函数文件中可以含有一个或多个子函数 子函数由 function 语句引导 主函数必须位于最前面，子函数出现的次序任意 子函数只能被主函数和位于同一个函数文件中的其它子函数调用 除了用 global 定义的全局变量外，子函数中的变量都是局部变量，子函数与主函数及其它子

11、函数之间通过输入、输出参数进行数据传递,调用一个函数时，Matlab 会首先检查该函数是否为一个子函数。,课堂练习,例题. 将A矩阵用for循环写出来,A = 1 2 3 4 5 2 3 4 5 6 3 4 5 6 7 4 5 6 7 8,例：编写m函数，计算函数值 x (x=10),例题：从键盘输入若干个数，当输入0时结束输入，求这些数的平均值和它们之和。 程序如下： sum=0; cnt=0; val=input(Enter a number (end in 0):); while (val=0) sum=sum+val; cnt=cnt+1; val=input(Enter a number (end in 0):); end if (cnt 0) sum mean=sum/cnt end,例题：求100，500之间前5个能被21整除的整数。,clc,clear, s=1; n=100:500; for i=1:length(n) if rem(n(i),21)=0 m(s)=n(i); s=s+1; continue, end if s=6 break