Matlab-程序设计和M文件课件

第5章程序设计和M文件

5.1程序控制5.2M文件结构5.3函数的使用5.4程序的调试5.5匿名函数、inline对象和函数句柄5.6函数绘图5.7数值分析1精选2021版课件5.1程序控制三种常用的结构：顺序结构、分支结构和循环结构。4种程序流程控制语句：分支控制语句、循环控制语句、错误控制语句和流程控制语句。2精选2021版课件5.1.1分支控制语句——

实现满足一定条件就执行相应分支的功能。1.if结构

if条件1

语句段1

elseif条件2

语句段2 ......

else

语句段n

end3精选2021版课件例5－1x=-5:0.1:5;ifx>=1y=x.^2-1;elseif-1<x<1y=0;elsey=-x.^2-1;endplot(x,y);4精选2021版课件2.switch结构——

常用于各种条件的列举。switch表达式

case值1

语句段1

case值2

语句段2 ...

otherwise

语句段n

end说明:(1)表达式只能是标量或字符串；(2)case后面的值可以是标量、字符串或元胞数组;(3)switch和end必须配对使用。5精选2021版课件例5-2

使用switch结构判断学生成绩的等级，90分以上为优，80～90为良，70～80为中，60～70为及格，60分以下为不及格。score=98;s1=fix(score/10); %取十位数switchs1

case{9,10}s='优'

case8s='良'

case7s='中'

case6s='及格'

otherwises='不及格'ends=优6精选2021版课件5.1.2循环控制语句——

实现将某段程序重复执行。1.for循环

for循环变量=array 循环体

end

说明：array可以是向量也可以是矩阵，循环执行的次数就是array的列数，每次循环中循环变量依次取array的各列并执行循环体，直到array所有列取完。练习7精选2021版课件forn=1:5 %循环5次forn=-1:0.1:1 %循环21次forn=linspace(-2*pi,2*pi,5) %循环5次a=eye(2,3);forn=a %循环3次，n为列向量例5－3x=eye(5);len=length(x)forn=1:lenx(n,n)=nendlen=5x=1000001000001000001000001x=1000002000001000001000001…x=10000020000030000040000058精选2021版课件or:x=eye(5);forn=1:length(x)x(n,:)=x(n,:)*nend9精选2021版课件2.while循环——常用于预先知道循环条件或循环结束条件的情况。

while条件表达式 循环体

end例5-4使用while循环计算f(x)=sin(x)，x<5且每次x递增ln(x)。x=0.1;whilex<5

f=sin(x)x=x+log(x)endf=-4.0651e+002-1.5921e+003ix=5.5640+9.2717i10精选2021版课件3.break和continue语句——控制循环的流程。（1）break语句break语句使包含break的最内层for或while循环强制终止，并立即跳出该循环结构，执行end后面的命令，break一般与if语句结合使用。（2）continue语句continue语句与break不同的是continue只结束本次for或while循环，而继续进行下次循环，continue一般也与if语句结合使用。11精选2021版课件eg.for循环中的break语句forii=1:5ifii==3;break;endfprintf('ii=%d\n',ii);enddisp('Endofloop!');结果:ii=1ii=2Endofloop!12精选2021版课件eg.for中的continue语句forii=1:5ifii==3;continue;end

fprintf('ii=%d\n',ii);enddisp('Endofloop!');结果:ii=1ii=2ii=4ii=5Endofloop!13精选2021版课件【例5－5】计算f(x)=ex，使用break语句当f(x)>1000时终止计算。x=0:10;forn=xf=exp(n);

iff>1000break

endendn=7f=1.0966e+00314精选2021版课件例5-6使用for循环将字符串中的数值取出，遇到非数值则跳过。str='Theresultis100.';len=length(str);s=[];forn=1:len

ifstr(n)>'9'|str(n)<'0' %非数值时continue

ends=[sstr(n)]ends=1s=10s=10015精选2021版课件相当于:str='Theresultis100.';len=length(str);s=[]forn=1:len

ifstr(n)>='0‘&str(n)<='9's=[sstr(n)]endend16精选2021版课件5.1.3错误控制语句——

当程序可能出现运行错误时，可以使用错误控制结构来捕获和处理错误，避免程序出错而不能继续运行。

try

语句段1

catch

语句段2

end执行过程:先试探地执行语句段1，如果出现错误则将错误信息赋给保留的lasterr变量，并放弃语句段1转而执行语句段2中的语句；如果语句段2正确则结束；当语句段1和语句段2都错误，则程序出错。如果语句段1正确则不执行语句段2就结束；17精选2021版课件例5－7>>a=[135];>>b=[123];>>tryc=a.*b'catchc=a.*b;end>>cc=1615>>[lastmsg,lastid]=lasterrlastmsg=Errorusing==>timesMatrixdimensionsmustagree.lastid=MATLAB:dimagree18精选2021版课件例5-7-2a=[123]b=[456]tryc=a*b;disp('tryisexcute')catchc=a.*b;disp('catchisexcute')endc19精选2021版课件5.1.4流程控制命令1.return命令return命令:(1)用于结束程序的执行，并立即返回到上一级调用函数(2)终止键盘方式注意：当程序进入死(无限)循环时，则按Ctrl+break键来终止程序的运行。2.keyboard命令keyboard命令用来使程序暂停运行，等待键盘命令，命令窗口出现“K>>”提示符，当键盘输入“return”后，程序才继续运行。keyboard命令可以用来在程序调试或程序执行时修改变量。20精选2021版课件a=[135];keyboardtryc=a.*b';catchc=a.*b;endc21精选2021版课件3.input命令a=input(‘Placeenteradata:’);

r=input(‘str’,‘s’)%从键盘中输入字符串保存到变量r4.disp命令disp命令是较常用的显示命令，常用来显示字符串型的信息提示。5.pause命令pause命令用来使程序暂停运行，当用户按任意键才继续执行。常用于程序调试或查看中间结果，也可以用来控制执行的速度。pause的命令格式如下：pause(n)

%暂停n秒练习22精选2021版课件5.2M文件结构MATLAB7.3的程序如果要保存则使用扩展名是“.m”的M文件，M文件有两种即M脚本文件（ScriptFile）和M函数文件（FunctionFile）。23精选2021版课件5.2.1M文件的一般结构1.函数声明行函数声明行是在M函数文件的第一行，只有M函数文件必须有，以“function”引导并指定函数名、输入和输出参数，M脚本文件没有函数声明行。2.H1行H1行是帮助文字的第一行。3.帮助文本4.程序代码

程序代码由MATLAB语句和注释语句构成。24精选2021版课件5.2.2M脚本文件和M函数文件1.M脚本文件格式:通常为一连串指令，没有函数声明行，无输入和输出参数.M脚本文件的说明：（1）运行脚本文件时，只是简单地按顺序从文件中读取一条条命令，送到MATLAB命令窗口中去执行；（2）M脚本文件运行产生的变量都驻留在MATLAB的基本工作空间中，在命令窗口中运行的命令可以使用这些变量；（3）脚本文件的命令可以访问工作空间的所有数据，因此要注意避免工作空间和脚本文件中的同名变量相互覆盖，一般在M脚本文件的开头使用“clear”命令清除工作空间的变量。25精选2021版课件例5-8

在M文件编辑／调试器窗口中编写M脚本文件得出10~50范围内的所有素数。%EX5_8计算10～50内的素数cleary1=primes(50); %获得50以内的所有素数n=length(primes(10)); %获得10以内素数的个数y=y1(n+1:end)在命令窗口中输入：>>ex5_8y=1113171923293137414347在工作空间中就可以查看到变量y1、n和y，并可以修改和使用这些变量。练习

PRIMES(N)isarowvectoroftheprimenumberslessthanorequaltoN.Aprimenumberisonethathasnofactorsotherthan1anditself.26精选2021版课件2.M函数文件格式：有函数声明行，可以输入参数和输出参数。M函数文件的说明如下：（1）M函数文件在运行过程中产生的变量都存放在函数本身的工作空间中。函数的工作空间随具体的M函数文件调用而产生，并随调用结束而删除，如果运行多个函数则产生多个临时的函数工作空间；（2）当文件执行完最后一条命令或遇到“return”命令时就结束函数文件的运行，同时函数工作空间的变量被清除；（3）一个M函数文件至少要定义一个函数。27精选2021版课件函数声明行的格式：function[输出参数列表]=函数名(输入参数列表)

说明：函数名是函数的名称，保存时最好函数名与文件名一致，当不一致时，MATLAB以文件名为准；输入参数列表是函数接收的输入参数，多个参数间用“,”分隔；输出参数列表是函数运算的结果，多个参数间用“,”分隔。28精选2021版课件函数定义的形式举例：定义多输入输出参数的函数function[x,y,z]=sphere(theta,phi,rho)无输入参数的函数functiony=sphere(orsphere())无输出参数的函数functionsphere(x)(orfunction[]=sphere(x))29精选2021版课件例5-9

将例5-7的计算行向量乘积的运算使用M函数文件保存。functionc=ex5_9(a,b)%EX5_9(a,b)returnstheproductofaandb%cistheproduct%aisthemultiplier%bisthemultiplicand%copyright2007-10-01

tryc=a.*b';catchc=a.*b;end调用该函数：>>z=ex5_9([123],[456])z=4101830精选2021版课件例：functiondistance=dist2(x1,y1,x2,y2)%DIST2Calculatethedistancebetweentwopoints%FunctionDIST2calculatesthedistancebetween%twopoints(x1,y1)and(x2,y2)inaCartesiancoordinatesystem.%%Callingsequence:%res=dist2(xl,y1,x2,y2)%Definevariables:%x1--x-positionofpoint1%y1--y-positionofpoint1%x2--x-positionofpoint2%y2--y-positionofpoint2%distance--Distancebetweenpoints%Recordofrevisions:%DateProgrammerDescriptionofchange%12/15/98S.J.ChapmanOriginalcode%Calculatedistance.distance=sqrt((x2-x1).^2+(y2-y1).^2);31精选2021版课件使用该函数的脚本文件如下：%Scriptfile:test_dist2.m%Purpose:%Thisprogramtestsfunctiondist2.%Recordofrevisions:%DateProgrammerDescriptionofchange%12/15/98S.J.Chapmanoriginalcode%Definevariables:%ax--x-positionofpointa%ay--y-positionofpointa%bx--x-positionofpointb%by--y-positionofpointb%result--Distancebetweenthepoints32精选2021版课件%Getinputdata.disp('Calculatethedistancebetweentwopoints:');ax=input('Enterxvalueofpointa:');ay=input('Enteryvalueofpointa:');bx=input('Enterxvalueofpointb:');by=input('Enteryvalueofpointb:');%Evaluatefunctionresult=dist2(ax,ay,bx,by);%Writeoutresult.fprintf('Thedistancebetweenpointsaandbis%f\n',result);33精选2021版课件直接在命令窗口调用该函数：>>test_dist2Calculatethedistancebetweentwopoints:Enterxvalueofpointa:1Enteryvalueofpointa:1Enterxvalueofpointb:4Enteryvalueofpointb:5Thedistancebetweenpointsaandbis5.00000034精选2021版课件>>helpdist2DIST2CalculatethedistancebetweentwopointsFunctionDIST2calculatesthedistancebetweentwopoints(xl,yl)and(x2,y2)inaCartesiancoordinatesystem.Callingsequence:res=dist2(xl,y1,x2,y2)……35精选2021版课件5.3函数的使用模块化程序设计:将一个大任务分成若干小的子任务,用一个一个的函数完成这些子任务,通过函数间的调用完成复杂任务.优点:易于维护和修改程序.注意:主函数可以调用子函数,子函数可以调用其他子函数或本身,但子函数不能调用主函数.36精选2021版课件5.3.1主函数和子函数一个M函数文件中可以包含一个或多个函数.1.主函数主函数是出现在文件最上方的函数，即第一行声明的函数.一个M文件只能有一个主函数，通常主函数名与M函数文件名相同。2.子函数（1）子函数的次序无任何限制；（2）子函数只能被同一文件中的函数（主函数或子函数）调用，不能被其它文件的函数调用；（3）同一文件的主函数和子函数运行时的工作空间是相互独立的。eg.ex5_1037精选2021版课件5.3.2函数的输入输出参数Matlab函数调用过程实际上也是参数的传递过程.函数A调用了函数B,称A为“主调函数”,而B为“被调用函数”.函数通过输入参数接收数据,经函数运算后由输出参数输出结果.38精选2021版课件1.参数的传递函数的参数传递是将主调函数中的变量值传给被调函数的输入参数。（1）函数参数传递的是数值；（2）参数的存储空间：被调函数的输入参数存放在函数的工作空间中，与MATLAB的工作空间是独立的，当调用结束时函数的工作空间被清除，输入参数就消失了。eg.ex5_1039精选2021版课件2.输入输出参数的个数Matlab中调用函数时的参数个数和函数声明时的不一定相同，可用专用函数确定实际调用时输入输出参数的个数。(1)nargin和nargout函数

nargin('fun')

%获取函数fun的输入参数个数

nargout('fun')

%获取函数fun的输出参数个数说明：fun是函数名，当nargin和nargout函数在函数体内时fun可省略，在函数外时fun不省略。40精选2021版课件补充:可选参数MATLAB的许多函数有可选参数,它们如何判断用到了哪些参数以及进行了什么操作呢?可以用8个特殊函数来获得可选参数的信息,及报告参数中的错误.41精选2021版课件▲nargin该函数返回调用函数时实际用的输入参数的个数.▲nargout该函数返回调用函数时实际用的输出参数的个数.▲nargchk当调用函数时用了太少或太多输入参数时,该函数返回一个标准的错误信息.▲error当参数有致命错误时,显示错误信息并异常中断产生错误的函数.▲warning当参数无致命错误时,显示警告信息且函数继续执行.▲inputname该函数根据参数个数返回实际的变量名.42精选2021版课件>>helpnargchkNARGCHKValidate(验证)numberofinputarguments.

MSG=NARGCHK(LOW,HIGH,N)returnsanappropriateerrormessagestringifNisnotbetweenLOWandHIGH.Ifitis,NARGCHKreturnsanemptymatrix.Example:

error(nargchk(1,3,nargin,'struct'))

43精选2021版课件eg.使用可选参数function[mag,angle]=polar_value(x,y)%POLAR_VALUEConverts(x,y)to(r,theta)%FunctionPOLAR_VALUEconvertsaninput(x,y)%valueinto(r,theta),withthetaindegrees.%Itillustratestheuseofoptionalarguments,%Definevariables:%angle--Angleindegrees%msg--Errormessage%mag--Magnitude44精选2021版课件%x--Inputxvalue%y--Inputyvalue(optional)%Recordofrevisions:%DateProgrammerDescriptionofchange%12/16/98S.J..Chapmanoriginalcode%Checkforalegalnumberofinputarguments.msg=nargchk(1,2,nargin);error(msg);45精选2021版课件%Iftheyargumentismissing,setitto0.ifnargin<2y=0;end%Checkfor(0,0)inputarguments,andprintout%awarningmessage.ifx==0&y==0msg='Bothxandyarezero:angleismeaningless!';warning(msg);end46精选2021版课件%Nowcalculatethemagnitude.mag=sqrt(x.^2+y.^2);%Ifthesecondoutputargumentispresent,calculate%angleindegrees.ifnargout==2angle=atan(y,x)*180/pi;end47精选2021版课件执行函数:>>[magangle]=polar_value???Errorusing==>polar-valueNotenoughinputarguments.>>[magangle]=polar_value(1,-1,1)???Errorusing==>polarvalueToomanyinputarguments.>>[magangle]=polar_value(1)mag=1angle=048精选2021版课件>>[magangle]=polar_value(1,-1)mag=1.4142angle=-45>>mag=polar_value(1,-1)mag=1.4142>>[magangle]=polar_value(0,0)Warning:Bothxandyarezero:angleismeaningless!Ind:\book\matlab\chap5\polar_value.matline32mag=0angle=049精选2021版课件(2)varargin和varargout函数用于处理复杂输入输出参数，varargin和varargout函数需要与分支结构结合使用，将函数调用时实际传递的参数构成元胞数组，通过访问元胞数组中各元素内容来获得输入输出变量。格式:

functiony=fun(varargin)%输入参数为varargin的函数funfunctionvarargout=fun(x)%输出参数为varargout的函数funfunctionvarargout=fun(varargin)50精选2021版课件5.3.3局部变量、全局变量和静态变量变量的作用域和生存期：作用域是变量的作用范围，生存期是变量的生存时间。根据变量作用域的不同,可将其分为局部变量和全局变量。根据变量生存期的不同，可将其分为自动变量和静态变量。注意：函数名与同一作用域的变量名不能重复,否则函数无法被调用.51精选2021版课件1.局部变量局部变量（LocalVariables）的作用范围是定义该变量的函数内部。声明局部变量：不需要特别声明变量的存储：在独立的函数工作空间中变量结束：变量所在的函数执行完毕时。52精选2021版课件2.全局变量全局变量具有全局作用域，可在不同函数和Matlab工作空间中共享。作用：减少参数传递过程，提高程序执行效率。缺点：全局变量在任何定义过的函数中都可以修改，使用时应十分小心。53精选2021版课件（1）声明全局变量全局变量在使用前必须用“global”声明，而且每个要共享全局变量的函数和工作空间，都必须逐个用“global”对该变量加以声明。（2）清除全局变量使用clear命令，命令格式如下：clearglobal变量名 %清除某个全局变量clearglobal %清除所有的全局变量

eg.li5_1354精选2021版课件3.静态变量格式：persistentx;注意:它只能在function里声明,且只有这个function才能认识它.声明后需初始化，注意不能直接赋值;而是先用isempty(x)判断x是否已经赋值，若没有，则可以赋值给x。如：functiony=test()

persistenta;

ifisempty(a)

a=0;

end

a=a+1;

y=a;之后第一次调用y=test;结果y=1，第二次调用y=test;结果y=2。就是说a记录了每次调用function后的结果。如果a为自动变量，则每次调用y=test后结果都是1.55精选2021版课件5.3.4嵌套函数、私有函数和重载函数1.嵌套函数在MATLAB中一个函数的内部还可以定义一个或多个函数，这种定义在其他函数内部的函数就称为嵌套函数。此时，每个函数必须用end结束。调用原则：(1)外层函数可以调用内一层函数，但不能隔层调用；(2)同层嵌套函数可以互相调用(3)内层函数可以调用外一层函数及和外层函数同层的其他函数。56精选2021版课件functionA(x)functionB(x,y)functionC(z)endendfunctionD(x)functionE(z,y)endendendA可以调用：B，DB可以调用：C，D，AC可以调用：B，D，ED可以调用：E，B，AE可以调用：D，B，C57精选2021版课件2.私有函数私有函数是限制访问权限的函数，私有函数存放在“private”子目录中，只能被其直接父目录的M函数文件所调用。3.重载函数重载函数是指两个函数使用相同的名称，处理的功能相似，但参数类型或个数不同，重载函数通常放在不同的文件夹下，文件夹名称以“@”开头后面跟一个数据类型名。58精选2021版课件5.3.5函数的工作过程和P码文件1.函数的搜索过程当在MATLAB中输入一个标识符时，首先确认是不是变量名，若不是，则做如下搜索：检查是否是本M函数文件内部的子函数；检查是否是“private”目录下的私有函数；检查是否在当前路径中；检查是否在搜索路径中。59精选2021版课件2.P码文件P码就是伪代码（Pseudocode），一个M文件第一次被调用时，MATLAB就将其进行编译并生成P码文件存放在内存中，生成的P码文件与原M文件名相同，其扩展名为“.p”，P码文件的保密性好。pcodeFile1.m,File2.m......-inplace %生成File1.p,File2.p……文件60精选2021版课件3.函数的工作空间每一个M函数运行时都有一个内存区，称为函数的工作空间。清除函数的工作空间：clearfunctions

%清除所有编译过的M函数文件和MEX文件工作空间clearfunctionfunname%清除某个编译过的函数工作空间61精选2021版课件5.4程序的调试

5.4.1直接检测（1）对于需要检测的变量可以通过删除语句行末尾的分号，或在程序的适当位置加显示变量值的语句，将结果显示在命令窗口中；（2）以脚本文件的方式来执行该函数。调试函数时，可以将该函数的第一句函数声明行前加“%”，并给输入参数赋值；（3）在程序的适当位置添加“keyboard”语句，当程序运行至此句会暂停运行，并在命令窗口显示“k>>”提示符，这时就可以在命令窗口查看和修改各变量的内容。62精选2021版课件5.4.2M文件编辑／调试器窗口1.“Debug”菜单（1）Step（F10）（2）Stepin（F11）和Stepout（shift+F11）（3）Run/Continue（F5）（4）Set/ClearBreakpoints（F12）（5）Set/ModifyConditionalBreakpoint…（6）StopifErrors/Warnings…（7）ExitDebugMode63精选2021版课件5.5匿名函数、inline对象和函数句柄

5.5.1匿名函数匿名函数是面向命令行代码的函数形式，通常只有一句很简单的语句。fhandle=@(arg1,arg2,......)(expr)

%创建匿名函数说明：fhandle是函数句柄；arg1,arg2,......是参数列表，也可以省略；expr是函数表达式。例：5-16>>fhnd1=@(x)(1+exp(-x)); %创建匿名函数>>rf1=fhnd1(2) %调用匿名函数rf1=1.135364精选2021版课件5.5.2inline对象创建inline对象就是使用inline函数将字符串转换成inline对象。inline_fun=inline(‘string’,arg1,arg2,…)

%创建inline对象例：5-17>>f=inline('sin(x)*exp(-z*x)','x','z') %创建inline对象f>>class(f) %查看f的类型ans=inline>>y1=f(1,2) %调用inline对象fy1=0.1139>>y2=feval(f,1,2)y2=0.113965精选2021版课件5.5.3函数句柄函数句柄（Function_Handle）包含了函数的路径、函数名、类型以及可能存在的重载方法，即函数是否为内部函数、M或P文件、子函数、私有函数等。1.创建函数句柄fhandle=@fun %创建函数句柄2.函数句柄的调用[y1,y2,…]=fhandle(arg1,arg2…)%调用函数句柄fhandle[y1,y2,…]=feval(fhandle,arg1,arg2…)[y1,y2,…]=feval('fun',arg1,arg2…)66精选2021版课件5.6函数绘图

5.6.1fplot命令fplot命令可以绘制函数的曲线。fplot(fun,limits,tol,Linespec)

%绘制函数fun的曲线fplot(fun,limits,n)

说明：fun是函数句柄或函数名；limits是自变量的取值范围[xminxmax]或[xminxmaxyminymax]；tol是相对误差度，默认为2e-3，可省略；Linespec是线型，与plot命令的设置相同，可省略；n是绘制的点数，当n≥1时至少绘制n+1个点。例：>>fhnd=@(x)(exp(-x).*sin(x));>>fplot(fhnd,[020],'b-.*')67精选2021版课件5.6.2ezplot命令ezplot命令也是用于绘制函数在某一自变量区域内的图形，ezplot命令格式如下：ezplot(fun,[min,max])

%绘制函数fun的曲线ezplot(x,y,[tmin,tmax])

说明：fun是函数句柄或函数名；[min,max]是自变量的范围，可省略，省略时默认为[-2*pi,2*pi]，当fun是二元函数，则范围是[xminxmaxyminymax]；x和y是指x(t)和y(t)；[tmin,tmax]是t的范围，可省略，省略时范围是[0,2*pi]。例：>>ezplot('x^2+y^2-1',[-3,3,-3,3])68精选2021版课件5.7数值分析

5.7.1求最小值和过零点1.一元函数的最小值fminbnd函数可以获得一元函数在给定区间内的最小值。x=fminbnd(fun,x1,x2)%寻找最小值2.多元函数的最小值fminsearch函数可以获得多元函数的最小值。[x,fval]=fminsearch(fun,x0) %寻找最小值3.一元函数的过零点x=fzero