《1.3 MATLAB程序设计》由会员分享,可在线阅读,更多相关《1.3 MATLAB程序设计(34页珍藏版)》请在金锄头文库上搜索。
1、第3章 MATLAB程序设计,3.1 M文件3.2 程序控制结构3.3 函数文件3.4 程序调试,3.1.1 M文件概述 用MATLAB语言编写的程序,称为M文件。 M文件可以根据调用方式的不同分为两类:命令文件(Script File)和函数文件(Function File)。,3.1 M文件,例3-1 分别建立命令文件和函数文件,将华氏温度f转换为摄氏温度c。程序1:首先建立命令文件并以文件名f2c.m存盘。clear; %清除工作空间中的变量f=input(Input Fahrenheit temperature:);c=5*(f-32)/9然后在MATLAB的命令窗口中输入f2c:In
2、put Fahrenheit temperature:73c = 22.7778,程序2:首先建立函数文件f2c.m。function c=f2c(f)c=5*(f-32)/9然后在MATLAB的命令窗口调用该函数文件。clear;y=input(Input Fahrenheit temperature:);x=f2c(y)输出情况为:Input Fahrenheit temperature:70c = 21.1111x = 21.1111,3.1.2 M文件的建立与打开 1建立新的M文件 2打开已有的M文件,1数据的输入从键盘输入数据,则可以使用input函数来进行,该函数的调用格式为: A
3、=input(提示信息,选项);其中提示信息为一个字符串,用于提示用户输入什么样的数据。 如果在input函数调用时采用s选项,则允许用户输入一个字符串。例如,想输入一个人的姓名,可采用命令:xm=input(Whats your name?,s);,3.2 程序控制结构,2数据的输出 MATLAB提供的命令窗口输出函数主要有disp函数,其调用格式为 disp(输出项) 其中输出项既可以为字符串,也可以为矩阵。,例3-2 输入x,y的值,并将它们的值互换后输出。 程序如下: x=input(Input x please.); y=input(Input y please.); z=x; x=
4、y; y=z; disp(x); disp(y);,例3-3 求一元二次方程ax2 +bx+c=0的根。程序如下:a=input(a=?);b=input(b=?);c=input(c=?);d=b*b-4*a*c;x=(-b+sqrt(d)/(2*a),(-b-sqrt(d)/(2*a),3程序的暂停 暂停程序的执行可以使用pause函数,其调用格式为: pause(延迟秒数) 如果省略延迟时间,直接使用pause,则将暂停程序,直到用户按任一键后程序继续执行。 若要强行中止程序的运行可使用Ctrl+C命令。,3.2.2 选择结构1if语句 在MATLAB中,if语句有3种格式。 (1) 单
5、分支if语句: if 条件 语句组 end 当条件成立时,则执行语句组,执行完之后继续执行if语句的后继语句,若条件不成立,则直接执行if语句的后继语句。,(2) 双分支if语句: if 条件 语句组1 else 语句组2 end 当条件成立时,执行语句组1,否则执行语句组2,语句组1或语句组2执行后,再执行if语句的后继语句。,例3-4 计算分段函数的值。程序如下:x=input(请输入x的值:);if x=0 y= (x+sqrt(pi)/exp(2);else y=log(x+sqrt(1+x*x)/2;endy,(3) 多分支if语句: if 条件1 语句组1 elseif 条件2 语
6、句组2 elseif 条件m 语句组m else 语句组n end 语句用于实现多分支选择结构。,2switch语句 switch语句根据表达式的取值不同,分别执行不同的语句,其语句格式为: switch 表达式 case 表达式1 语句组1 case 表达式2 语句组2 case 表达式m 语句组m otherwise 语句组n end,当表达式的值等于表达式1的值时,执行语句组1当表达式的值等于表达式2的值时,执行语句组2,当表达式的值等于表达式m的值时,执行语句组m当表达式的值不等于case所列的表达式的值时,执行语句组n。当任意一个分支的语句执行完后,直接执行switch语句的下一句。
7、,例3-6 某商场对顾客所购买的商品实行打折销售,标准如下(商品价格用price来表示): price200 没有折扣 200price500 3%折扣 500price1000 5%折扣 1000price2500 8%折扣 2500price 0) sum1 mean=sum1/cntend,3break语句和continue语句一般与if语句配合使用。break语句用于终止循环的执行。当在循环体内执行到该语句时,程序将跳出循环,继续执行循环语句的下一语句。continue语句控制跳过循环体中的某些语句。当在循环体内执行到该语句时,程序将跳过循环体中所有剩下的语句,继续下一次循环。,例3-
8、12 求100,200之间第一个能被21整除的整数。程序如下:for n=100:200 if rem(n,21)=0 continue end breakendn,4循环的嵌套一个循环结构的循环体又包括一个循环结构 例3-13 若一个数等于它的各个真因子之和,则称该数为完数,如6=1+2+3,所以6是完数。求1,500之间的全部完数。for m=1:500s=0;for k=1:m/2if rem(m,k)=0s=s+k;endendif m=s disp(m);endend,3.3.1 函数文件的基本结构 函数文件由function语句引导,其基本结构为: function 输出形参表=函
9、数名(输入形参表) 注释说明部分 函数体语句 其中以function开头的一行为引导行,表示该M文件是一个函数文件。函数名的命名规则与变量名相同。输入形参为函数的输入参数,输出形参为函数的输出参数。当输出形参多于一个时,则应该用方括号括起来。,3.3 函数文件,例3-14 编写函数文件求半径为r的圆的面积和周长。 函数文件如下: function s,p=fcircle(r)%calculate the area and perimeter of a circle of radii r%r 圆半径%s 圆面积%p 圆周长s=pi*r*r;p=2*pi*r;,3.3.2 函数调用 函数调用的一般格式是: 输出实参表=函数名(输入实参表) 要注意的是,函数调用时各实参出现的顺序、个数,应与函数定义时形参的顺序、个数一致,否则会出错。函数调用时,先将实参传递给相应的形参,从而实现参数传递,然后再执行函数的功能。,