MATLAB程序设计第三讲.ppt

《MATLAB程序设计第三讲.ppt》由会员分享，可在线阅读，更多相关《MATLAB程序设计第三讲.ppt（58页珍藏版）》请在得力文库 - 分享文档赚钱的网站上搜索。

1、MATLAB程序设计程序设计 杨凯杨凯2主主 要要 内内 容容文件设计文件设计一、一、预备知识：变量及运算符预备知识：变量及运算符二、二、M文件介绍文件介绍三、程序控制结构三、程序控制结构四、函数文件四、函数文件五、程序调试五、程序调试3一、变量及运算符一、变量及运算符1 1 变量变量n1.1 变量的作用范围变量的作用范围n 局部变量：每个局部变量：每个MATLAB函数都有自己的局部变函数都有自己的局部变量。局部变量的作用范围仅限于本函数，一旦运量。局部变量的作用范围仅限于本函数，一旦运行超过本函数，变量的值将不再保留。行超过本函数，变量的值将不再保留。n 全局变量：用全局变量：用global

2、关键字进行声明，作用范围关键字进行声明，作用范围整个整个M文件。文件。4【例例】全局变量应用示例。全局变量应用示例。先建立函数文件先建立函数文件wadd.m，该函数将输入的参数加，该函数将输入的参数加权相加。权相加。function f=wadd(x,y)global ALPHA BETAf=ALPHA*x+BETA*y;在命令窗口中输入：在命令窗口中输入：global ALPHA BETAALPHA=1;BETA=2;s=wadd(1,2)5将变量设置成全局变量的两种方式：将变量设置成全局变量的两种方式：（1）将该变量作为函数参数进行传递；）将该变量作为函数参数进行传递；（2）将该变量声明为

3、全局变量；）将该变量声明为全局变量；推荐：采用函数方式。推荐：采用函数方式。1.2变量的存活期：变量的存活期：Persistent关键词声明：静态变量。固定变量关键词声明：静态变量。固定变量（1）只能用在函数中；（）只能用在函数中；（2）其他函数不能获取）其他函数不能获取（3）函数存在时）函数存在时MATLAB不将它们从内存清除不将它们从内存清除（4）清除函数或编辑该）清除函数或编辑该M文件，文件，MATLAB清除所清除所有固定变量。有固定变量。6n1.3 表达式表达式n 运算符：运算符：（1）算术运算符：）算术运算符：运算符说明+加-减.*乘./右除.左除7n算术运算符算术运算符x=ab是方

4、程是方程a*x=b的解；的解；x=b/a是方程是方程x*a=b的解。的解。.幂.转置复数共轭转置*矩阵相乘/矩阵右除矩阵左除矩阵的幂8n比较运算符比较运算符比较运算符比较两个数的大小比较运算符比较两个数的大小运算符运算符描述描述运算符运算符描述描述=大于或等于大于或等于大于大于=不等于不等于9n逻辑运算符逻辑运算符逻辑运算符判断对象或对象之间某种逻辑关系逻辑运算符判断对象或对象之间某种逻辑关系针对数组：针对数组：&and与与|or或或not非非举例：举例：A=0 1 1 0 1；B=1 1 0 0 1;xor两个数组中相同位置的元素只有一个非两个数组中相同位置的元素只有一个非0时时 返回返回1

5、，否则返回，否则返回0；（异或）；（异或）any()如果矢量中有一个元素非如果矢量中有一个元素非0，返回，返回1；否；否 则返回则返回0；10n举例：举例：C=0 1 2;0-3 8;0 5 0;all()：如果矢量中所有元素非：如果矢量中所有元素非0，返回，返回1，否则返回，否则返回0表达式：表达式：&符号两端表达式为真，符号两端表达式为真，true|符号两端表达式有一个为真，符号两端表达式有一个为真，true 1.4运算符的优先级运算符的优先级小括号小括号转置、幂、复数共轭转置、矩阵的幂转置、幂、复数共轭转置、矩阵的幂否否乘除等乘除等加减加减冒号操作符冒号操作符比比较运算符较运算符ANDO

6、R&|保险：不确定的加括号保证保险：不确定的加括号保证11n处理字符串表达式处理字符串表达式n eval:处理包含处理包含MATLAB表达式、语句或函数调用表达式、语句或函数调用的字符串的字符串neval(string);n【例例】：t=1/(m+n-1)for m=1:k for n=1:k a(m,n)=eval(t);end end12二、文件介绍二、文件介绍n用用MATLAB语言编写的程序，称为语言编写的程序，称为M文件。文件。M文文件可以根据调用方式的不同分为两类：命令文件件可以根据调用方式的不同分为两类：命令文件(Script File)和函数文件和函数文件(Function Fi

7、le)。n命令文件：没有输入参数，也不返回输出参数。命令文件：没有输入参数，也不返回输出参数。n函数文件：可以输入参数，也可返回输出参数。函数文件：可以输入参数，也可返回输出参数。13nM M文件的建立与编辑文件的建立与编辑 M文件是一个文本文件，它可以用任何编辑文件是一个文本文件，它可以用任何编辑程序来建立和编辑，而一般常用且最为方便的是程序来建立和编辑，而一般常用且最为方便的是使用使用MATLAB提供的文本编辑器。提供的文本编辑器。n建立新的建立新的M M文件：文件：（1）菜单操作）菜单操作FileOpen （2）命令操作）命令操作 edit （3）命令按钮操作）命令按钮操作n编辑已有的编

8、辑已有的M M文件：文件：同上同上14【例例】分别建立命令文件和函数文件，将华氏温度分别建立命令文件和函数文件，将华氏温度f转换为摄氏温度转换为摄氏温度c。程序程序1：首先建立命令文件并以文件名首先建立命令文件并以文件名f2c.m存盘。存盘。clear;%清除工作空间中的变量清除工作空间中的变量f=input(Input Fahrenheit temperature：);c=5*(f-32)/9然后在然后在MATLAB的命令窗口中输入的命令窗口中输入f2c，将会执行该，将会执行该命令文件，执行情况为：命令文件，执行情况为：Input Fahrenheit temperature：73c=22.

9、777815程序程序2：首先建立函数文件：首先建立函数文件f2c.m。function c=f2c(f)c=5*(f-32)/9然后在然后在MATLAB的命令窗口调用该函数文件。的命令窗口调用该函数文件。y=input(Input Fahrenheit temperature：);x=f2c(y)输出情况为：输出情况为：Input Fahrenheit temperature：70c=21.1111x=21.1111 16三、程序控制结构三、程序控制结构n 1数据的输入数据的输入 从键盘输入数据，则可以使用从键盘输入数据，则可以使用input函数来进行，函数来进行，该函数的调用格式为：该函数的

10、调用格式为：A=input(提示信息，选项提示信息，选项)；其中提示信息为一个字符串，用于提示用户输其中提示信息为一个字符串，用于提示用户输入什么样的数据。入什么样的数据。如果在如果在input函数调用时采用函数调用时采用s选项，则允许选项，则允许用户输入一个字符串。例如，想输入一个人的姓用户输入一个字符串。例如，想输入一个人的姓名，可采用命令：名，可采用命令：xm=input(Whats your name?,s);172、数据的输出、数据的输出 MATLAB提供的命令窗口输出函数主要有提供的命令窗口输出函数主要有disp函函数，其调用格式为数，其调用格式为 disp(输出项输出项)其中输出

11、项既可以为字符串，也可以为矩阵。其中输出项既可以为字符串，也可以为矩阵。例如例如 A=Hello,MATLAB;disp(A)输出为：输出为：Hello,MATLAB18【例例】输入输入x,y的值，并将它们的值互换后输出。的值，并将它们的值互换后输出。程序如下：程序如下：x=input(Input x please.);y=input(Input y please.);z=x;x=y;y=z;disp(x);disp(y);19【例例】求一元二次方程求一元二次方程ax2+bx+c=0的根。的根。程序如下：程序如下：a=input(a=?);b=input(b=?);c=input(c=?);d

12、=b*b-4*a*c;x=(-b+sqrt(d)/(2*a),(-b-sqrt(d)/(2*a);disp(x1=,num2str(x(1),x2=,num2str(x(2);203程序的暂停程序的暂停 暂停程序的执行可以使用暂停程序的执行可以使用pause函数，其调用函数，其调用格式为：格式为：pause(延迟秒数延迟秒数)如果省略延迟时间，直接使用如果省略延迟时间，直接使用pause，则将暂，则将暂停程序，直到用户按任一键后程序继续执行。停程序，直到用户按任一键后程序继续执行。若要强行中止程序的运行可使用若要强行中止程序的运行可使用Ctrl+BREAK命令。命令。214、流程控制、流程控制

13、4.1选择结构选择结构if语句语句 在在MATLAB中，中，if语句有语句有3种格式。种格式。(1)单分支单分支if语句：语句：if 条件条件 语句组语句组 end当条件成立时，则执行语句组，执行完之后继续执当条件成立时，则执行语句组，执行完之后继续执行行if语句的后继语句，若条件不成立，则直接执行语句的后继语句，若条件不成立，则直接执行if语句的后继语句。语句的后继语句。22(2)双分支双分支if语句：语句：if 条件条件 语句组语句组1 else 语句组语句组2 end 当条件成立时，执行语句组当条件成立时，执行语句组1，否则执行语句组，否则执行语句组2，语句组，语句组1或语句组或语句组2

14、执行后，再执行执行后，再执行if语句的后继语句的后继语句。语句。23【例例】计算分段函数的值。计算分段函数的值。程序如下：程序如下：x=input(请输入请输入x的值的值:);if xA(3)&A(1)+A(3)A(2)&A(2)+A(3)A(1)p=(A(1)+A(2)+A(3)/2;s=sqrt(p*(p-A(1)*(p-A(2)*(p-A(3);disp(s);else disp(不能构成一个三角形。不能构成一个三角形。)end 25(3)多分支多分支if语句：语句：if 条件条件1 语句组语句组1 elseif 条件条件2 语句组语句组2 elseif 条件条件m 语句组语句组m el

15、se 语句组语句组n end语句用于实现多分支选择结构。语句用于实现多分支选择结构。26【例例】输入一个字符，若为大写字母，则输出其对应的小写输入一个字符，若为大写字母，则输出其对应的小写字母；若为小写字母，则输出其对应的大写字母；若为数字母；若为小写字母，则输出其对应的大写字母；若为数字字符则输出其对应的数值，若为其他字符则原样输出。字字符则输出其对应的数值，若为其他字符则原样输出。c=input(请输入一个字符请输入一个字符,s);if c=A&c=a&c=0&c=9 disp(abs(c)-abs(0);else disp(c);end27nSwitch语句语句根据表达式的值不同，执行不

16、同的语句，其格式为：根据表达式的值不同，执行不同的语句，其格式为：switch 表达式表达式 case 表达式表达式1 语句组语句组1 case 表达式表达式2 语句组语句组2 case 表达式表达式m 语句组语句组m otherwise 语句组语句组n end28【例例】根据变量根据变量 num 的值来决定显示的内容。的值来决定显示的内容。n num=input(请输入一个数请输入一个数);n switch numn case-1n disp(I am a teacher.);n case 0n disp(I am a student.);n case 1n disp(You are a t

17、eacher.);n otherwisen disp(You are a student.);n end29【例】某商场对顾客所购买的商品实行打折销售，某商场对顾客所购买的商品实行打折销售，标准如下标准如下(商品价格用商品价格用price来表示来表示)：price200 没有折扣没有折扣 200price500 3%折扣折扣 500price1000 5%折扣折扣 1000price2500 8%折扣折扣 2500price 0)sum mean=sum/cntend40n循环的嵌套：n如果一个循环结构的循环体又包括一个循环结构，就称为循环的嵌套，或称为多重循环结构。n多重循环的嵌套层数可以是

18、任意的。可以按照嵌套层数，分别叫做二重循环、三重循环等。处于内部的循环叫作内循环，处于外部的循环叫作外循环。41【例例】求求100,1000以内的全部素数。以内的全部素数。n=0;for m=100:1000 flag=1;j=m-1;i=2;while i=j&flag if rem(m,i)=0 flag=0;end i=i+1;end if flag n=n+1;prime(n)=m;end end prime%变量prime存放素数42【例例】若一个数等于它的各个真因子之和，则称该数若一个数等于它的各个真因子之和，则称该数为完数，如为完数，如6=1+2+3，所以，所以6是完数。求是完数

19、。求1,500之之间的全部完数。间的全部完数。for m=1:500s=0;for k=1:m/2if rem(m,k)=0s=s+k;endendif m=s disp(m);endend43break语句和语句和continue语句语句 与循环结构相关的语句还有与循环结构相关的语句还有break语句和语句和continue语句。它们一般与语句。它们一般与if语句配合使用。语句配合使用。break语句用于终止循环的执行。当在循环体内执行语句用于终止循环的执行。当在循环体内执行到该语句时，程序将跳出循环，继续执行循环语到该语句时，程序将跳出循环，继续执行循环语句的下一语句。句的下一语句。con

20、tinue语句控制跳过循环体中的某些语句。当在循语句控制跳过循环体中的某些语句。当在循环体内执行到该语句时，程序将跳过循环体中所环体内执行到该语句时，程序将跳过循环体中所有剩下的语句，继续下一次循环。有剩下的语句，继续下一次循环。44【例例】求求100，200之间第一个能被之间第一个能被21整除的整数。整除的整数。程序如下：程序如下：for n=100:200if rem(n,21)=0 continueendbreakendn 45四、函数文件四、函数文件n函数文件是另一种形式的函数文件是另一种形式的M文件，每一个函数文件文件，每一个函数文件都定义一个函数。事实上，都定义一个函数。事实上，M

21、ATLAB提供的标准提供的标准函数大部分都是由函数文件定义的。函数大部分都是由函数文件定义的。1、函数文件的基本结构函数文件的基本结构 函数文件由函数文件由function语句引导，其基本结构为：语句引导，其基本结构为：function 输出形参表输出形参表=函数名函数名(输入形参表输入形参表)注释说明部分注释说明部分 函数体语句函数体语句46n其中以其中以function开头的一行为引导行，表示该开头的一行为引导行，表示该M文文件是一个函数文件。函数名的命名规则与变量名件是一个函数文件。函数名的命名规则与变量名相同。输入形参为函数的输入参数，输出形参为相同。输入形参为函数的输入参数，输出形参

22、为函数的输出参数。当输出形参多于一个时，则应函数的输出参数。当输出形参多于一个时，则应该用方括号括起来。该用方括号括起来。n【例】编写函数文件求小于任意自然数n的Fibonacci数列各项。47function f=ffib(n)%用于求Fibonacci数列的函数文件%f=ffib(n)f=1,1;i=1;while f(i)+f(i+1)n f(i+2)=f(i)+f(i+1);i=i+1;end 48【例】编写函数文件求半径为编写函数文件求半径为r的圆的面积和周长。的圆的面积和周长。函数文件如下：函数文件如下：function s,p=fcircle(r)%CIRCLE calculat

23、e the area and perimeter of a circle of radii r%r 圆半径圆半径%s 圆面积圆面积%p 圆周长圆周长%2004年年7月月30日编日编s=pi*r*r;p=2*pi*r;492、函数调用、函数调用函数调用的一般格式是：函数调用的一般格式是：输出实参表输出实参表=函数名函数名(输入实参表输入实参表)要注意的是，函数调用时各实参出现的顺序、个要注意的是，函数调用时各实参出现的顺序、个数，应与函数定义时形参的顺序、个数一致，否数，应与函数定义时形参的顺序、个数一致，否则会出错。函数调用时，先将实参传递给相应的则会出错。函数调用时，先将实参传递给相应的形参

24、，从而实现参数传递，然后再执行函数的功形参，从而实现参数传递，然后再执行函数的功能。能。50【例】利用函数文件，实现直角坐标利用函数文件，实现直角坐标(x,y)与极坐标与极坐标(,)之间之间的转换。的转换。函数文件函数文件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);在在MATLAB中，函数可以嵌套调用，即一个

25、函数可以调用别中，函数可以嵌套调用，即一个函数可以调用别的函数，甚至调用它自身。一个函数调用它自身称为函数的函数，甚至调用它自身。一个函数调用它自身称为函数的递归调用。的递归调用。51【例】利用函数的递归调用，求n！。n function f=factor(n)n if n=1n f=1;n elsen f=factor(n-1)*n;n endn return;%返回n 在命令文件main2.m中调用函数文件factor.m：n for i=1:10n fac(i)=factor(i);n endn fac程序运行结果是：fac=Columns 1 through 6 1 2 6 24 12

26、0 720 Columns 7 through 10 504 40320 362880 362880052n3 函数所传递参数的可调性函数所传递参数的可调性n MATLAB在函数调用上有一个与众不同之处：函数所传递参数数目的可调性。凭借这一点，一个函数可完成多种功能。n 在调用函数时，MATLAB用两个永久变量nargin和nargout分别记录调用该函数时的输入实参和输出实参的个数。只要在函数文件中包含这两个变量，就可以准确地知道该函数文件被调用时的输入输出参数个数，从而决定函数如何进行处理。53【例】nargin用法示例n函数文件examp.m：n function fout=charra

27、y(a,b,c)n if nargin=1n fout=a;n elseif nargin=2n fout=a+b;n elseif nargin=3n fout=(a*b*c)/2;n end命令文件mydemo.m：x=1:3;y=1;2;3;examp(x)examp(x,y)examp(x,y,3)执行mydemo.m后的输出是：ans=1 2 3 ans=2 4 6 ans=2154调试器调试器1Debug菜单项菜单项该菜单项用于程序调试，需要与该菜单项用于程序调试，需要与Breakpoints菜单项菜单项配合使用。配合使用。2Breakpoints菜单项菜单项 该菜单项共有该菜单项

28、共有6个菜单命令，前两个是用于在程个菜单命令，前两个是用于在程序中设置和清除断点的，后序中设置和清除断点的，后4个是设置停止条件的，个是设置停止条件的，用于临时停止用于临时停止M文件的执行，并给用户一个检查文件的执行，并给用户一个检查局部变量的机会。局部变量的机会。55习习 题题1、猜数游戏。首先由计算机产生、猜数游戏。首先由计算机产生1,100之间之间的随机整数，然后由用户猜测所产生的随机的随机整数，然后由用户猜测所产生的随机数。根据用户猜测的情况给出不同提示，如数。根据用户猜测的情况给出不同提示，如猜测的数大于产生的数，则显示猜测的数大于产生的数，则显示“High”，小于则显示小于则显示“

29、Low”，等于则显示，等于则显示“You won”，同时退出游戏。用户最多可以猜，同时退出游戏。用户最多可以猜7次。次。56习习 题题1、猜数游戏。首先由计算机产生、猜数游戏。首先由计算机产生1,100之间之间的随机整数，然后由用户猜测所产生的随机的随机整数，然后由用户猜测所产生的随机数。根据用户猜测的情况给出不同提示，如数。根据用户猜测的情况给出不同提示，如猜测的数大于产生的数，则显示猜测的数大于产生的数，则显示“High”，小于则显示小于则显示“Low”，等于则显示，等于则显示“You won”，同时退出游戏。用户最多可以猜，同时退出游戏。用户最多可以猜7次。次。572、从键盘、从键盘输入

30、5个数，求其中最大数和最小数。并在屏幕上提示输入“请输入您要显示数字的序号”，输入后，回显“您要显示数字为（此处为实际需要显示序号对应的数字）”3.求Fibonacci数列 (1)大于4000的最小项。(2)5000之内的项数。58 作业格式要求作业格式要求1、交电子文档和纸质文档（打印）；、交电子文档和纸质文档（打印）；2、电子文档可以、电子文档可以U盘提交，也可以发到电子邮箱盘提交，也可以发到电子邮箱:yangkai_3、电子文档文件名命名规则、电子文档文件名命名规则：学号：学号+姓名姓名如如20050001杨凯杨凯4、两周内提交有效。过期无效。、两周内提交有效。过期无效。5、作业参与期末成绩、作业参与期末成绩