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） elsestatements（语句组） end,if 条件语句,switch expression （表达式）case value1 （表达式1）
3、statement1（语句组1）case value2 （表达式2）statement2（语句组2）. .case valuem （表达式m）statementm（语句组m）otherwisestatement （语句组） end,switch 语句,for variable=expressionstatement（循环体） 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 direction % % X = 1 2
5、 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:nyh_old=yh; k2=ceil(k/2);for l=2:k2yh(l)=yh_old(l-1)+yh_old(l);endyh(k2+1:k)=yh(k-k2:-1:1); di
6、sp(yh); end,函数文件举例,printyh.m,函数调用的一般格式,输出实参列表=函数名(输入实参列表),函数调用时，实参的顺序应与函数定义时的形参的顺序一致。 函数可以嵌套调用，即一个函数可以被其它函数调用，甚至可以被它自身调用，此时称为递归调用。函数所传递的参数具有可调性，Matlab 用两个永久变量 nargin 和 nargout 分别记录调用该函数时的输入实参和输出实参的个数。,函数调用,% 函数文件 factor.m function f=factor(n) if (n=1)f=1; elsef=n*factor(n-1); end,% main.m % 计算 s=1!+
7、2!+3!+4!+5! % s=0; for i=1:5s=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=2)p=a+b; q=a-b; end,函数通过输入和输出参数与
8、其它 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) global a product=x*y; a=x+y;,子函数
10、,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) = 1m = x(n+1)/2); elsem = (x(n/2)+x(n/2+1)/2; end,子函数,函数文件中可以含有一个或多个子函数子函数由 function 语句引导主函数必须位于最前面，子函数出现的次序任意子函数只能被主函数和位于同一个函数文件中的其它子函数调用除了用 global 定义的全局变量外，子函数中的变量都是局部变量，子函数与主函数及其它子函数之间通过输入、输出参数进行数据传递,调用一个函数时，Matlab 会首先检查该函数是否为一个子函数。,
