matlab编程二.ppt

1、数学实验,Matlab 编程 (二),用 Matlab 语言编写的程序称为 M 文件。,M 文件,M 文件以 .m 为扩展名。,M 文件的建立与打开,关系运算,关系操作符能用来比较两个同样大小的数组，或用来比较一个数组和一个标量，在后一种情况，标量和数组中的每一个元素相比较，比较结果与数组大小一样。,关系运算符,逻辑运算,逻辑运算符,all(x)如果向量 X 中所有元素都非零，则返回 1，否则返回 0,any(x) 如果向量 X 中存在非零元素，则返回 1，否则返回 0,若 x 为矩阵，则 any 和 all 按列运算，返回一个 0-1 向量,运算优先级,高,低,disp(X),输出变量 X

2、的值，可以是数值矩阵或字符串,相关命令,pause 或 pause(n),fprintf(fid,format,variables),按指定的格式将变量的值输出到指定的文件,暂停程序的运行，以秒为单位,if expression1 （条件1） statements1（语句组1） elseif expression2 （条件2） statements2（语句组2） . . elseif expressionm （条件m） statementsm（语句组m） else statements（语句组） end,if 条件语句,switch expression （表达式） case value1 （

3、表达式1） statement1（语句组1） case value2 （表达式2） statement2（语句组2） . . case valuem （表达式m） statementm（语句组m） otherwise statement （语句组） end,switch 语句,for variable=expression statement（循环体） end,循环,while expression （条件） statement（循环体） end,break 、 continue、return,break 和 continue 一般与 if 语句配合使用,return 语句用于退出正在运行的脚

4、本或函数，通常用在函数文件中。,function 输出形参列表=函数名(输入形参列表) % 注释说明部分(可选) 函数体语句(必须),第一行为引导行，表示该 M文件是函数文件 函数名的命名规则与变量名相同 ( 必须以字母开头 ) 当输出行参多于一个时，用方括号括起来 函数必须是一个单独的 M文件 函数文件名必须与函数名一致 以百分号开始的语句为注释语句,函数文件,函数文件由 function 语句引导,function y = myfliplr(x) % returns X with row preserved and columns % flipped in the left/right d

5、irection % % 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); for l=2:k2 yh(l)=yh_old(l-1)+yh_old(l); end yh(

6、k2+1:k)=yh(k-k2:-1:1); disp(yh); end,函数文件举例,printyh.m,函数调用的一般格式,输出实参列表=函数名(输入实参列表),函数调用时，实参的顺序应与函数定义时的形参的顺序一致。 函数可以嵌套调用，即一个函数可以被其它函数调用，甚至可以被它自身调用，此时称为递归调用。 函数所传递的参数具有可调性，Matlab 用两个永久变量 nargin 和 nargout 分别记录调用该函数时的输入实参和输出实参的个数。,函数调用,% 函数文件 factor.m function f=factor(n) if (n=1) f=1; else f=n*factor(n

7、-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!,递归函数举例,% 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=

8、2) p=a+b; q=a-b; end,函数通过输入和输出参数与其它 M 文件进行数据传递。,当函数调用完毕后，该函数文件中定义的所有局部变量都将被释放，即全部被清除。,如果在若干个 M文件中，都把某个变量定义为全局变量，则这些函数将公共使用这一变量。所有函数都可以对它进行存取和修改操作。,函数文件中的变量都是局部的，即一个函数文件中定义的变量不能被另一个函数文件或其它 M 文件使用。,全局变量与局部变量,定义全局变量是 M文件间传递信息的一种手段。,变量名列表中的各个变量用空格隔开，不能用逗号！,全局变量给函数间的数据传递带来了方便，但却破坏了函数对变量的封装，降低了程序的可读性，因而在结

9、构化程序设计中，全局变量是不受欢迎的。特别是当程序较大，子程序较多时，全局变量将个程序调试和维护带来不便，故不提倡使用全局变量。,global 变量名列表,全局变量的定义,程序示例,clear; global a b a=1; b=3; 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)

10、 global a product=x*y; a=x+y;,子函数,function avg, med = ex4subfun(x) % 主函数 n = length(x); avg = mean(x, n); med = median(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,子函数,函数文件中可以含有一个或多个子函数 子函数由 funct