第4章M文件和流程控制结构

1、第第4 4章章 M M文件和流程控制结构文件和流程控制结构n用户在进行MATLAB操作时，经常需要输入较多的命令语句，而且有些还需要重复地输入，如果仍然按照上述在MATLAB命令窗口（Command Window）进行直接交互式命令操作方式，就显得比较烦琐和费时。为了改进这种情况，MATLAB提供文件编程操作方式。nMATLAB中所有的文件和函数都分别对应一个M文件（以.m为后缀的文本文件），它可以在任何文本编辑器中进行显示和编辑，最方便的是在MATLAB程序编辑器窗口中进行（图1-2）。M文件的代码是纯ASCII码字符，可以在MATLAB命令窗口（Command Window）键入命令：ty

2、pe function_name来查看。4.1 M4.1 M文件的分类文件的分类n4.1.1 4.1.1 命令文件（命令文件（Script FileScript File）n命令文件也称为脚本文件，通常只是一系列命令语句的组合，没有输入参数，也不返回参数，运行时只需要在命令窗口键入文件名即可。n例例4-1 4-1 试编制一命令文件，计算圆的面积与周长n% 计算圆的面积与周长的命令文件（rsc.m）nr=input (请输入圆的半径（mm）r=); % 使用键盘输入圆的半径rns=pi*r2; % 计算圆的面积nc=2*pi*r; % 计算圆的周长nfprintf(圆的半径 r = % 3.4f

3、 mm n,r)nfprintf(圆的面积 s = % 3.4f mm2 n,s)nfprintf(圆的周长 c = % 3.4f mm n,c)n rscn请输入圆的半径（mm）r=35n圆的半径 r = 35.0000 mm n圆的面积 s = 3848.4510 mm2 n圆的周长 c = 219.9115 mm n rscn请输入圆的半径（mm）r=58n圆的半径 r = 58.0000 mm n圆的面积 s = 10568.3177 mm2 n圆的周长 c = 364.4247 mmn命令文件在运行过程中可以调用MATLAB系统工作域内的所有数据，所产生的变量也是全局变量，并且一直保

4、存在工作域内，直至quit时为止。 n4.1.2 4.1.2 函数文件函数文件nMATLAB函数或函数文件是系统中已经设计好的为完成某一种特定的运算或是实现某一种特定的功能的子程序， MATLAB系统几十种工具箱中提供了丰富的函数库，它们是MATLAB系统的重要组成部分。在使用这些函数或函数文件时是作为命令来对待的，所以函数又称为函数命令。n函数文件常用于需要反复调用和不断改变参数的场合，它可以接受参数，也可以返回参数，它的运行必须通过语句调用。n1、函数文件的第一个可执行必须是函数文件定义语句：function输出参数=函数名（输入参数）n2、函数文件可以带有多个输出参数和输入参数，也可以没

5、有输出参数。当输入参数或是输出参数多于一个时，参数之间一般用逗号隔开。n3、函数文件的命名规则与变量名的规则相同。两个函数的函数名即使相同，但是输入参数或输出参数的数目不同，它们的含义也不同。n4、函数文件运行时，自动创建其独自的函数工作空间（function workspace），过程变量都存储这个工作空间中。函数文件执行之后只保留最后结果，不保留之间任何过程，所定义的变量也只要在函数内部起作用，并且随着该函数文件调用的结束而自动被清除。n5、应当指出，以function开头的函数文件是一个独立M文件，必须单独储存在内存中MATLAB系统所在的工作目录：D:MATLAB7work中。n将函数

6、文件存盘时，缺省状态下系统自动储存为“函数名.m”的函数文件（即函数文件名与函数名相同）。n如果函数文件名与函数名不同，MATLAB在调用时将以函数文件名为准（忽略函数名）。n例例4-2 4-2 试使用调用函数文件方法，计算渐开线齿轮的齿廓函数值。n% 1-建立计算渐开线函数的函数文件（jkxhs.m）nfunction theta,rk=jkxhs(alpha,rb)nhd=pi/180;ntheta=tan(hd.*alpha)-hd.*alpha; % 渐开线展角nrk=rb./cos(hd.*alpha); % 渐开线向径n将该文件以文件名jkxhs.m保存在MATLAB工作区路径中。

7、tancosbkkrrn% 2-建立调用命令文件nm=input( 输入齿轮模数（mm） m = );nz=input( 输入齿轮齿数 z = );nha=input( 输入齿高系数 ha = );nhd=pi/180;nalpha_f=20; % 分度圆压力角n% 计算齿轮的几何尺寸nr=m*z/2; % 分度圆半径nrb=r*cos(alpha_f*hd); % 基圆半径nra=r+ha*m; % 齿顶圆半径nalpha_a=acos(rb/ra); % 齿顶圆压力角n% 调用计算渐开线函数的函数文件nalpha=0:alpha_a/hd/5:alpha_a/hd;ntheta,rk=jk

8、xhs(alpha,rb);thetak=theta/hd; n% 输出计算结果nfprintf( 齿轮模数 m = % 6.4f mm n,m)nfprintf( 齿轮齿数 z = % 3.0f n,z)nfprintf( 齿高系数 ha = % 3.4f n,ha)nfprintf( 基圆半径 rb = % 6.4f mm n,rb)nfprintf( 分度圆半径 r = % 6.4f mm n,r)nfprintf( 齿顶圆半径ra = % 6.4f mm n,ra)nfprintf( 齿顶圆压力角 alpha_a = % 6.4f n,alpha_a/hd)ndisp 渐开线展角 th

9、etak(),thetakndisp 渐开线向径 rk(mm),rkn输出的计算结果如下：n 齿轮模数 m = 5.0000 mm n 齿轮齿数 z = 28 n 齿高系数 ha = 1.0000 n 基圆半径 rb = 65.7785 mm n 分度圆半径 r = 70.0000 mm n 齿顶圆半径 ra = 75.0000 mm n 齿顶圆压力角 alpha_a = 28.7119 n 渐开线展角thetak()nthetak = 0 0.0193 0.1563 0.5386 1.3151 2.6721n 渐开线向径rk(mm)nrk = 65.7785 66.1102 67.1224

10、68.8680 71.4430 75.0000n有时为了在命令文件中直接描述数学函数，可采用函数inline()，它的功能是建立一个内联函数，其调用格式是：n fun=inline(函数表达式，自变量列表)n其中自变量列表中每个自变量之间均需要用单引号括起来，相互之间用逗号分隔。n例如，渐开线函数方程可以用函数inline()语句描述为ntheta=inline(tan(alpha*pi/180)-pi/180.*alpha,alpha);n其中，alpha（表示压力角）是函数式中自变量。n它只能返回1个变量，不支持较复杂的语句结构。 4.2 4.2 流程控制结构流程控制结构n4.2.1 4.

11、2.1 顺序结构顺序结构n顺序结构是根据语句排列的顺序自上而下执行。主要包括数据的输入、处理、结果的输出等内容n1 1、数据的输入、数据的输入（input函数）n从键盘输入数据，可以使用input函数来进行，该函数的调用格式为：n z=input(屏幕上显示的提示信息) n当执行该函数时，系统等待从键盘输入数据后按回车键，输入的数据就存入变量z中。n如果输入的是字符，则格式是n z=input(屏幕上显示的提示信息,s)n其中s表示键盘输入字符。n2 2、数据的输出、数据的输出（disp和 fprintf函数）n1、函数disp是自由格式输出数据，用于直接输出字符或是数值，例如n a=2,4,

12、6;n disp anan disp(a)n2 4 6n注意：以上两例自由格式输出中，“a”表示字符，“(a)”表示向量a的值。n2、函数fprintf在屏幕上输出格式化数据，一般调用格式是： fprintf（屏幕上显示的提示信息 %format n,v）n其中，format用来指定数据输出时采用的格式：d表示整数；e表示实数（科学计算法形式）；f 表示浮点实数（小数形式）；g表示由系统自动选取上述两种格式之一；s表示输出字符串；nn是换行符号，以避免下一个输出数据或是提示符号也挤在同一行；nV是输出变量或是矩阵。n例如：n fprintf(The pi is %4.8f n, pi)nThe

13、 pi is 3.14159265n可见，在2个单引号之间包括输出的字串“The pi is”、输出数据的格式“%4.8f”（表示按照浮点数格式，8位小数，输出变量pi（）的数据。n fprintf(34的5次方 = %2.6e n, 345)n34的5次方 = 4.543542e+007n可见，输出“345”值为指数格式的数据，6位小数。n4.2.2 4.2.2 循环结构循环结构n循环结构是指按照给定的条件重复执行指定的语句。MATLAB提供了两种循环结构：for循环结构和while循环结构。n1 1、forfor循环结构循环结构nfor循环结构的格式是：n for i=c0:cb:cnn

14、循环体n endn其中，c0是循环变量i的初值；cb是循环变量i的步长，当cb取单位步长时可以省略；cn是循环变量i的终值。n例例4-3 4-3 用筛选法求自然数2至m范围内的全部素数。素数是大于1，且除了1和它本身以外，不能被其他任何整数所整除的整数。nm=input( 输入素数的最大范围 m = );nx=2:m;nfor k=2:sqrt(m)n n=find(rem(x,k)=0 & x=k);n x(n)=; % 去掉非素数nendnxn运算结果：n 输入素数的最大范围 m = 30nx = 2 3 5 7 11 13 17 19 23 29n2 2、whilewhile循环

15、结构循环结构nwhile循环结构的格式是：n while(条件)n 循环体n endn其中的“条件”是一个关系表达式，当关系表达式的值为真时自动执行循环体，直到关系表达式的值条件为假，转向循环体结束标志end后面语句。n应当指出，while循环执行过程中会改变while后面关系表达式的值，否则这个while循环是一个死循环，因此while循环执行次数是不定的。n在for循环结构和while循环结构中，使用break命令可终止循环结构的运行。n例例4-5 4-5 水仙花数是指一个3位数的各位数字的立方和等于该数本身，用while循环结构检索水仙花数nm=100;nwhile (m=999);n

16、m3=fix(m/100); % 保存百位数字n m2=rem(fix(m/10),10); % 保存十位数字n m1=rem(m,10); % 保存个位数字n if m=m33+m23+m13; n disp(m)n endn m=m+1;nendn运算结果（写成一行）：n 153 370 371 407 n4.2.3 4.2.3 选择结构选择结构n选择结构是指按照给定的条件是否成立，分别执行不同的语句。MATLAB提供了两种选择结构：if选择结构和switch选择结构。n1 1、ifif选择结构选择结构n if (条件1)n 语句组1 % 满足条件1时执行该语句组n elseif (条件2

17、)n 语句组2 % 满足条件2时执行该语句组n n elsen 语句组n+1 % 上面的条件均不满足时n endn例例4-6 4-6 用选择结构建立样条函数的判断函数文件nfunction f=pdc(x) % 该函数文件没有输出参数nif x0n f=0;nelseif x1n f=x;nelseif x pdc(-3.5)nans =n0n pdc(0.8)nans =n 0.8000n pdc(2.8)nans =n 0 n2 2、switchswitch选择结构选择结构nswitch 开关表达式n case 常量表达式1n 语句组1n case 常量表达式2n 语句组2n n otherwisen 语句组nn endn当开关表达式的值与对应常量表达式的值匹配时，执行相应的语句组，然