1、第五章 MATLAB 程序设计,Matlab语言与其他计算机语言一样，可以用来编程。充分利用 Matlab数据结构的特点，可以使程序结构简单，编程效率高。本章主要介绍有关Matlab程序控制结构以及数学上常用到的程序设 计的基本方法。,教学过程,M文件介绍 程序控制结构 函数变量及其变量作用域 子函数与局部函数 程序设计的辅助函数 程序的调试和优化,5.1 M文件介绍,5.1.1 M文件概述用MATLAB语言编写的程序，称为M文件。M文件是由若干Matlab命 令构成的，它可以完成某些操作，也可以实现某种算法。实际上，Matlab提供的内部函数以及各种工具箱，都是利用Matlab命 令开发的M
2、文件。用户也可以结合自己的工作需要，开发具体的程序或 工具箱。,M文件可以根据调用方式的不同分为两类：命令文件(Script File)和函 数文件(Function File)。它们的扩展名均为m，主要区别在于：（1）命令文件没有输入参数，也不返回输出参数，而函数文件可以带 输入参数，也可以返回输出参数。 （2）命令文件对Matlab工作空间中的变量进行操作文件中所有命令的 执行结果也完全返回工作空间中，而函数文件中定义的变量为局部变 量，当函数文件执行完毕时，这些变量被清除。 （3）命令文件可以直接运行，在Matlab命令窗口输入命令文件的名 字，就会顺序执行命令文件中的命令，而函数文件不
3、能直接运行，要以 函数调用的方式来调用。,函数m文件的格式：function 返回变量=函数名（输入变量）注释说明语句段程序语句段 特定规则： （1）函数m文件第一行必须以单词function作为引导词，必 须遵循如下形式：function =() （2）m文件的文件名必须是 .m。 （3）程序中的变量均为局部变量，不保存在工作空间中。其变量只在函数运行期间有效。,例5-1 分别建立命令文件和函数文件，将华氏温度f转换为摄氏温度c。 程序1： 首先建立命令文件并以文件名f2c.m存盘。 clear; %清除工作空间中的变量 f=input(Input Fahrenheit temperatur
4、e：); c=5*(f-32)/9 然后在MATLAB的命令窗口中输入f2c，将会执行该命令文件，执行情况为： Input Fahrenheit temperature：73 c =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：70 c =21.1111 x =21.1111,例5-2 建立一个命
5、令文件将变量a,b的值互换，然后运行该命 令文件。首先建立命令文件并以文件名exch.m存盘：clear;a=1:10; b=11,12,13,14;15,16,17,18;c=a;a=b;b=c;ab然后在MATLAB的命令窗口中输入exch，将会执行该命令 文件。,程序2：建立一个函数文件将变量a,b的值互换，然后在命令窗口 调用该函数文件。首先建立函数文件fexch.m：function a,b=exch(a,b)c=a;a=b;b=c;然后在MATLAB的命令窗口调用该函数文件：clear;x=1:10; y=11,12,13,14;15,16,17,18;x,y=fexch(x,y)
6、,5.1.2 M文件的建立与打开M文件是一个文本文件，它可以用任何编辑程序来建立和 编辑，而一般常用且最为方便的是使用MATLAB提供的文 本编辑器meditor。 1建立新的M文件为建立新的M文件，启动MATLAB文本编辑器有3种方法：(1) 菜单操作。从MATLAB主窗口的File菜单中选择New菜单项，再选择M-file命令，屏幕上将出现MATLAB 文本编辑器窗口。 (2) 命令操作。在MATLAB命令窗口输入命令edit，启动MATLAB文本编辑器后，输入M文件的内容并存盘。 (3) 命令按钮操作。单击MATLAB主窗口工具栏上的New M-File命令按钮，启动MATLAB文本编辑
7、器后，输入M文件的内容并存盘。,2打开已有的M文件打开已有的M文件，也有3种方法： （1）菜单操作。从MATLAB主窗口的File菜单中选择Open 命令，则屏幕出现Open对话框，在Open对话框中选中所需 打开的M文件。在文档窗口可以对打开的M文件进行编辑修 改，编辑完成后，将M文件存盘。(2) 命令操作。在MATLAB命令窗口输入命令：edit 文件 名，则打开指定的M文件。(3) 命令按钮操作。单击MATLAB主窗口工具栏上的Open File命令按钮，再从弹出的对话框中选择所需打开的M文件。,5.2 程序控制结构,Matlab语言的程序结构与其它高级语言是一致的，分为顺 序结构，选择
8、结构，循环结构。5.2.1 顺序结构 1数据的输入 从键盘输入数据，则可以使用input函数来进行，该函数的 调用格式为：A=input(提示信息，选项)； 其中提示信息为一个字符串，用于提示用户输入什么样的数 据。如果在input函数调用时采用s选项，则允许用户输入一个 字符串。例如，想输入一个人的姓名，可采用命令： xm=input(Whats your name?,s);,2数据的输出MATLAB提供的命令窗口输出函数主要有disp函、 fprintf函数。 disp调用格式为disp(输出项) 其中输出项既可以为字符串，也可以为矩阵。,例5-3 求一元二次方程ax2 +bx+c=0的根
9、。 程序如下： 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); disp(x1=,num2str(x(1),x2=,num2str(x(2);,fprintf函数最常见的使用方式用以下例子说明.若键入命令fprintf（圆周率pi=%10.9f，pi） 则会按浮点型输出含9位小数，1位整数的圆周率近似值，其 输出结果为 圆周率pi=3.141592654若键入命令n=23;fprintf(n=%d,n) 则会按整型数输出n值，其输出结果为n=23,若键入命
10、令n=23;fprintf(n=%f,n) n=23.0000003程序的暂停暂停程序的执行可以使用pause函数，其调用格式为：pause(延迟秒数)如果省略延迟时间，直接使用pause，则将暂停程序，直 到用户按任一键后程序继续执行。若要强行中止程序的运行可使用Ctrl+C命令。,5.2.2 选择结构 1if语句在MATLAB中，if语句有3种格式。(1) 单分支if语句： if 条件语句组end当条件成立时，则执行语句组，执行完之后继续执行if语 句的后继语句，若条件不成立，则直接执行if语句的后继语 句。,(2) 双分支if语句： if 条件语句组1else语句组2end当条件成立时，
11、执行语句组1，否则执行语句组2，语句组 1或语句组2执行后，再执行if语句的后继语句。,例5-4 计算分段函数的值。程序如下： x=input(请输入x的值:); if x=0y= (x+sqrt(pi)/exp(2); elsey=log(x+sqrt(1+x*x)/2; end y,练习：计算以下分段函数的值,(3) 多分支if语句： if 条件1语句组1elseif 条件2语句组2elseif 条件m语句组melse语句组nend 语句用于实现多分支选择结构。,例5-5 输入一个字符，若为大写字母，则输出其对应的小写字母；若为小写字母，则输出其对应的大写字母；若为数字字符则输出其对应的数
12、值，若为其他字符则原样输出。 c=input(请输入一个字符,s); if c=A end,2switch语句switch语句根据表达式的取值不同，分别执行不同的语 句，其语句格式为：switch 表达式case 表达式1语句组1case 表达式2语句组2case 表达式m语句组motherwise语句组nend,当表达式的值等于表达式1的值时，执行语句组1，当表达 式的值等于表达式2的值时，执行语句组2，当表达式的 值等于表达式m的值时，执行语句组m，当表达式的值不等 于case所列的表达式的值时，执行语句组n。当任意一个分 支的语句执行完后，直接执行switch语句的下一句。,例5-6 某
13、商场对顾客所购买的商品实行打折销售，标准如下 (商品价格用price来表示)： price200 没有折扣200price500 3%折扣500price1000 5%折扣1000price2500 8%折扣2500price5000 10%折扣 5000price 14%折扣 输入所售商品的价格，求其实际销售价格。,程序如下： price=input(请输入商品价格); switch fix(price/100) case 0,1 %价格小于200rate=0;case 2,3,4 %价格大于等于200但小于500rate=3/100;case num2cell(5:9) %价格大于等于50
14、0但小于1000rate=5/100;case num2cell(10:24) %价格大于等于1000但小于2500rate=8/100;case num2cell(25:49) %价格大于等于2500但小于5000rate=10/100;otherwise %价格大于等于5000rate=14/100; end price=price*(1-rate) %输出商品实际销售价格,3try语句 语句格式为： try语句组1 catch语句组2 end try语句先试探性执行语句组1，如果语句组1在执行过程中 出现错误，则将错误信息赋给保留的lasterr变量，并转去执 行语句组2。,例5-7 矩
15、阵乘法运算要求两矩阵的维数相容，否则会出错。 先求两矩阵的乘积，若出错，则自动转去求两矩阵的点乘。 程序如下： A=1,2,3;4,5,6; B=7,8,9;10,11,12; tryC=A*B; catchC=A.*B; end C lasterr %显示出错原因,5.2.3 循环结构循环是指按照给定的条件，重复执行指定的语句，这是 一种十分重要的程序结构。1for语句 for语句的格式为： for 循环变量=表达式1:表达式2:表达式3循环体语句end 其中表达式1的值为循环变量的初值，表达式2的值为步长， 表达式3的值为循环变量的终值。步长为1时，表达式2可以 省略。,例5-8 一个三位
16、整数各位数字的立方和等于该数本身则称该 数为水仙花数。输出全部水仙花数。程序如下： for m=100:999 m1=fix(m/100); %求m的百位数字 m2=rem(fix(m/10),10); %求m的十位数字 m3=rem(m,10); %求m的个位数字 if m=m1*m1*m1+m2*m2*m2+m3*m3*m3 disp(m) end end,例5-9 已知 ，当n=100时，求y的值。程序如下： y=0; n=100; for i=1:ny=y+1/(2*i-1); end y,在实际MATLAB编程中，采用循环语句会降低其执行速 度，所以前面的程序通常由下面的程序来代替：
17、 n=100; i=1:2:2*n-1; y=sum(1./i); y,for语句更一般的格式为：for 循环变量=矩阵表达式循环体语句end执行过程是依次将矩阵的各列元素赋给循环变量，然后执 行循环体语句，直至各列元素处理完毕。,例5-10 写出下列程序的执行结果。s=0;a=12,13,14;15,16,17;18,19,20;21,22,23;for k=as=s+k;enddisp(s);,练习：编程输入范德蒙型的矩阵,2while语句while语句的一般格式为：while (条件)循环体语句end其执行过程为：若条件成立，则执行循环体语句，执行后 再判断条件是否成立，如果不成立则跳出
18、循环。,例5-11 从键盘输入若干个数，当输入0时结束输入，求这些 数的平均值和它们之和。 程序如下： sum=0; cnt=0; val=input(Enter a number (end in 0):); while (val=0)sum=sum+val;cnt=cnt+1;val=input(Enter a number (end in 0):); end if (cnt 0)summean=sum/cnt end,3break语句和continue语句与循环结构相关的语句还有break语句和continue语句。它 们一般与if语句配合使用。break语句用于终止循环的执行。当在循环体
19、内执行到该 语句时，程序将跳出循环，继续执行循环语句的下一语句。continue语句控制跳过循环体中的某些语句。当在循环体 内执行到该语句时，程序将跳过循环体中所有剩下的语句， 继续下一次循环。,例5-12 求100，200之间第一个能被21整除的整数。 程序如下： for n=100:200 if rem(n,21)=0continue end break end n,例5-13：问题：从1到n的任何一个自然数, 只要对n反复进行 下列两种运算: 1）如果n是偶数, 就除以2 ; 2）如果n是奇数, 就乘以3加1， 最后的结果总是1。 这个问题大约是在二十世纪五十年代被提出来的。在西方它 常
20、被称为西拉古斯Syracuse)猜想，因为据说这个问题首先是 在美国的西拉古斯大学被研究的；而在东方，这个问题由将 它带到日本的日本数学家角谷静夫的名字命名，被称作角谷 猜想。 角谷静夫曾用计算机验算到71011 ，并未出现反例。 1992年李文斯(G.T.Leavens)和孚门南(M.Vermeulen)也以计算 机对小于5.61013的正整数进行验证，也未发现反例。,n=input(请输入一个大于1的正整数 n=); if n1 if rem(n,2)=0 n=n/2elsen=n*3+1 end end,4循环的嵌套如果一个循环结构的循环体又包括一个循环结构，就称为 循环的嵌套，或称为多
21、重循环结构。例5-13 若一个数等于它的各个真因子之和，则称该数为完 数，如6=1+2+3，所以6是完数。求1,500之间的全部完数。,程序如下： for m=1:500 s=0; for k=1:m/2 if rem(m,k)=0 s=s+k; end end if m=sdisp(m); end end,5.2.5 程序举例 1、 猜数游戏。首先由计算机产生1,100之间的随机整数， 然后由用户猜测所产生的随机数。根据用户猜测的情况给出 不同提示，如猜测的数大于产生的数，则显示“您的数字较 大!”，小于则显示“您的数字较小!”，等于则显示“YOU WIN! ”，同时退出游戏。用户最多可以猜
22、6次。,y=round(10+89*rand(); for k=1:6; x=input(第,num2str(k),次输入一个两位数(输完请按回车):); if(xy) 您的数字较小! else if(x=y) msgbox(YOU WIN!); return; else 您的数字较大! end end if(k=6) msgbox(YOU LOSE! GAME OVER!) end end,2、用筛选法求某自然数范围内的全部素数。 素数是大于1，且除了1和它本身以外，不能被其他任何 整数所整除的整数。用筛选法求素数的基本思想是：要找出 2m之间的全部素数，首先在2m中划去2的倍数(不包括2)
23、， 然后划去3的倍数(不包括3)，由于4已被划去，再找5的倍数 (不包括5)，直到再划去不超过的数的倍数，剩下的数都 是素数。,m=input(m=); p=2:m;for i=2:sqrt(m)n=find(rem(p,i)=0endp,例：有五个水手带了一只猴子来到南太平洋的一个荒岛上，发现那里有 一大堆椰子。由于旅途的颠簸，大家都很疲倦，很快就入睡了。第一个 水手醒来后，把椰子平分成五堆，并将多余的一个椰子给了猴子，他私 藏了一堆后便又去睡了。第二、第三、第四、第五个水手也陆续起来， 和第个水手一样，把椰子平分成五堆后，将恰好多的一个也给猴子，并 私藏了一堆，再去入睡。天亮以后，大家把余
24、下的椰子重新等分成五 堆，每人分一堆，正好余一个再给猴子。试问原先共有几个椰子？,求解这一问题可以用递推算法。首先分析椰子数目的 变化规律，设最初的椰子数为p0，即第一个水手所处理 之前的椰子数，用p1、p2、p3、p4、p5 分别表示五个 水手对椰子动了手脚以后剩余的椰子数目，则根据问题 有 p(k+1)=4/5*(p(k)-1)再用x表示最后每个水手平分得到的椰子数，于是有 x=1/5*(p5-1) 所以 p5 = 5x +1 利用逆向递推的方法，有,n=input(input n:); for x=1:n p=5*x+1; for k=1:5 p=5*p/4+1; end if p=fi
25、x(p)breakend end disp(x,p),例、猴子吃桃问题：猴子第一天摘下若干个桃子，当即吃了一半，还不 过瘾，又多吃了一个，第二天又将剩下的桃子吃了一半，又多吃了一个。 以后每天早上都吃了前一天剩下的一半零一个。到到第十天早上想再 吃，只剩一个桃子。求第一天共摘多少桃子？,number=1; x(10)=number; fprintf(计算结果如下：n); for i=9:-1:1x(i)=(x(i+1)+1)*2;fprintf( 第%d天有%d个桃子n,i,x(i); end,例、统计一个字符串里有几个字母？程序如下：,function k=f(s) m,n=size(s);
26、 x=isletter(s); if x(1)+x(n)=2 %若该字符串首尾均是字母k=1; elseif x(1)+x(n)=1 %若该字符串首或尾有一个且仅有一个是字母k=0;elseif x(1)+x(n)=0 %若该字符串首尾均不是字母（比如标点符号或空格）k=-1; end for i=1:nif x(i)=1k=k+1; end end,5.3 函数变量及其变量作用域,在MATLAB语言中，变量可以分为输入变量、输出变量和函数内使用 的变量。输入变量相当于函数的入口数据，也是一个函数操作的主要对象，从 某种意义上说，函数的功能在于对输入变量进行一定的操作从而实现一 定的功能。函数
27、的输入变量为局部变量，函数对输入变量的一切操作和 修改如果不依靠输出变量的话，将不会影响工作区间中该变量的值。,5.3.1 变量的输入和输出规则MATLAB可以有任意数量的输入和输出变量。这些参数的特性和规 则如下： 函数式M文件可以没有输入和输出变量。 函数可以用比M文件中的函数定义行所规定的输入输出 变量更少的变量进行调用。但是不能用比规定的输入输出变 量更多的变量进行调用。 在一次调用中所用到的输入和输出变量的个数可以通过 分别调用函数nargin和nargout来确定。因为nargin和nargout是 函数而不是变量，所以用户不能用诸如nargin=nargin+pi之类 的语句对它
28、们进行重新赋值。 当一个函数被调用时，输入变量并没有被复制到函数的 工作区间中，但是它们的值在这个函数是可读的。应当注意 的是，如果输入变量的任何值被改变了，这个输入变量组就 被复制到了函数的工作区。,如果一个函数声明了一个或者多个输出变量，但是用户在使用的时 候又不想要输出参数，则只要不把输出变量赋值给任何变量就可以了； 或者在函数结束之前用函数clear删除这些变量。函数可以通过在函数声明中将varargin作为最后的输入参数，接受 可变的任意个数的输入参数。Varargin是一个预先定义的单元数组，这个 单元数组的第i个单元就是varargin程序的位置算起的第i个变量。通过函数声明行中
29、将varargout作为最后的输出变量，函数可以接受 任意个数的变量形式的输出参数。Varargout也是一个预先定义的单元数 组，这个单元数组的第i个单元就是从varargout的出现位置算起的第i个变 量。函数nargchk和nargoutch分别提供了对有效地输入和输出变量个数的 简单错误校验，因为如果函数调用的输入或者输出变量的个数多于函数 定义中出现的个数，函数都自动地返回一个错误，因此虽然这些函数的 作用有限，但是在一个函数定义声明了任意数目的输入变量和输出变量 的时候却是非常有用的。,例、nargin函数的初级使用方法。 解：本程序实现如下功能，当调用过程时小于或等于一个变量时，
30、系统 提示错误的输入，当有两个变量时，程序将两个数相加，当有3个变量 时，将前两个数相加并减去第3个。程序如下： function d=nargintest(a,b,c) if nargin=1 error(Not enough input arguments) elseif nargin=2 d=a+b; elseif nargin=3 d=a+b+c end,此外，使用nargin可以查找函数输入变量的个数，例如，想要查找上 述nargintest函数的参数个数，可以使用命令 nargin(nargintest)5.3.2 局部变量 局部变量是在函数内部使用的变量，其影响范围只能在本函数内
31、， 每个函数在运行时，都占有独立的函数工作空间，此工作空间和 MATLAB的工作空间时相互独立的，局部变量仅存在于函数的工作空间 内。当函数执行完毕之后，该变量即自行消失。,5.3.3 全局变量在Matlab中，函数内部定义的变量都是局部变量，它们不被加载到 工作区间中。有时，用户需要使用全局变量，这时要使用global函数来进 行定义，而且在任何使用该全局变量的函数中都应加以定义，即使是在 命令窗口也不例外。,例 全局变量应用示例。先建立函数文件wadd.m，该函数将输入的参数加权相加。function f=wadd(x,y)global ALPHA BETAf=ALPHA*x+BETA*y
32、;在命令窗口中输入：global ALPHA BETAALPHA=1;BETA=2;s=wadd(1,2),5.3.4 永久变量 除了通过全局变量共享数据外，函数式M文件还可以通过声明一个变 量persistent来对函数中重复使用和递归调用的变量的访问进行限制，使 用格式形如persistent (X Y Z)。永久变量与全局变量类似，但是它的范围 被限制在声明这些变量的函数内部，不允许在其他的函数中对它们进行 改变。只要M文件还在MATLAB 7的内存中，永久变量就存在。,5.4 函数的分类,5.4.1 主函数 M文件中的第一个函数就叫做主函数，前边章节中所引用的函数事实 上都是主函数，主
33、函数之后可以是任意数量的子函数， 它们可以作为主 程序的子程序。一般来说，在命令窗口或是其他的M文件只能调用主函 数，调用的时候就是直接调用其函数名。 比如，函数average的M文件average.m如下： function y = average(x) % AVERAGE Mean of vector elements. y = sum(x)/length(x); % Actual computation,5.4.2 匿名函数 匿名函数提供了一种创建简单程序的方法，使用它用户可以不必每次 都编写M文件。用户可以在MATLAB 7的命令窗口或是其他任意M文件 和脚本文件中使用匿名函数。 匿名
34、函数的格式如下所示。 fhandle = (arglist) expr,例、匿名函数举例 本例列举3个匿名函数的例子，分别有0个，1个和两个输入参数。 t=() datestr(now); t() 28-Oct-2010 11:21:34sqrt=(x) x.2 sqrt(12)sumAxBy=(x,y)(14*x+41*y) sumAxBy(3,7),5.4.3 嵌套式函数 在MATLAB 7中，可以在一个函数的内部定义一个或多个其他的函 数，这些在内部定义的函数被称作嵌套式函数，应当注意的是，在嵌套 式函数的内部也可以定义嵌套式函数。定义嵌套式函数时，只需在另一个M文件的内部定义该函数即可
35、，同 其他M文件一样， 嵌套式函数包含有M文件的所有基本部分。,5.4.4 子函数与其他的高级语言一样，在MATLAB 7语言中也可以很方便地定义子 函数，用来扩充函数的功能。在函数文件中题头定义的函数为主函数， 而在函数体内定义的其他函数都被视为子函数。子函数只能为主函数或 同一主函数下的其他子函数所使用。,例、编制一函数，要求任意输入两个数值后，用两个子函数分别求出它 们的和与它们的绝对值的和，在将这两个和相乘。 function ch=zihanshu(x,y) ch=zihanshu1(x,y)*zihanshu2(x,y); function ch=zihanshu1(x,y) ch
36、=abs(x)+abs(y); function ch=zihanshu2(x,y) ch=x+y;,5.4.5局部函数MATLAB 7语言中把放置在目录private下的函数称为局部函数，这些 函数只有private目录的父目录中的函数才可以调用，其他目录下的函数 不能调用。局部函数与子函数所不同的是，局部函数可以被其父目录下的所有函 数所调用，而子函数则只能被其所在的M文件的主函数所调用。所以， 局部函数在可用的范围上大于子函数；在函数编辑的结构上，局部函数 与一般的函数文件的编辑相同，而子函数只能在主函数文件中编辑。,5.4.6 函数句柄 函数句柄提供了一种间接访问函数的手段，用户可以很
37、方便地调用其 他函数；提供函数调用过程中的可靠性；减少程序设计中的冗余； 同时 可以在使用函数的过程中保存函数相关的信息，尤其是关于函数执行的 信息。函数句柄的创建定义函数句柄的定义只需在提示符后添加相应函数的函数名，使用 格式形如： fhandle=functionname 其中fhandle为所创建的函数句柄，functionname为所创建的函数 。,trigFun = sin, cos, tan; plot(trigFun2(-pi:0.01:pi),函数句柄的调用function x = plotFHandle(fhandle, data) plot(data, fhandle(da
38、ta) 在命令窗口中调用程序如下: plotFHandle(sin, -pi:0.01:pi),函数句柄的操作,5.4 程序设计的辅助函数,5.4.1 执行函数 执行函数及其功能,容错函数程序设计的好坏在很大程度上取决于其容错能力的大小。MATLAB 7 语言提供了相应的报错及警告函数error、warning、lasterr、lastwarn以 及errortrap on/off等，可以很方便地实现这方面的功能。,时间运算函数 常用的时间控制函数及其功能,5.5 程序的调试和优化,程序的调试 程序的错误种类 错误的识别和程序调试 程序的优化 程序优劣的分析 程序优化的技巧 有效地使用内存,5
39、.5.1 程序的错误种类,在MATLAB 7的表达式中可能存在两种类型的错误，即语法错误和运 行错误 ： 语法错误 语法错误发生在M文件程序代码的生成过程中，一般是由函数 参数输入类型有误或是矩阵运算阶数不符等引起。 运行错误 运行错误一般指在程序运行过程中，出现溢出或是死循环等异 常现象。, min(x,y) ? Undefined function or variable x. x=0;y=-9; y=minmun(x,y) ? Undefined command/function minmun. A=ones(3); B=magic(4); C=B*b ? Undefined funct
40、ion or variable b. D=A*B ? Error using = mtimes Inner matrix dimensions must agree.,5.5.2 错误的识别和程序调试,对于简单的问题，使用下面的一种或几种方法可以方便地求解：将函数中输出关键值的行的分号(;)去掉，这样，这些运算的中间结 果将在命令窗口中予以显示，用户可以据此来检查中间结果的正确性。在函数中添加一些语句，用来显示用户认为很重要的变量的值。使用keyboard命令中断程序，该命令实现函数工作区间和命令窗 口工作区间的交互，从而获得用户所需要的信息，使用该命令后，程 序将处于调试状态，此时命令窗口的
41、提示符由“”变为“K”，用户可 以进行相应的操作。在函数头前加“%”，这样就将函数式M文件变为脚本式M文件，而 脚本式M文件运行时，其工作区间就是MATLAB 7的工作区间，这样 在出现错误的时候就可以查询这个工作区间。,错误识别函数及其功能,程序断点的设置,程序优劣的分析,在MATLAB 7语言中，使用profile函数以及计时函数 tic 和toc来分析程序中各个部分的耗时情况，从而帮助 用户找出程序中需要改进的地方。其中profile在计算相 对耗时以及查找文件执行过程中瓶颈问题时更为有效， 而tic和toc函数在计算绝对耗时时更为有效。,程序优化的技巧,程序的向量化操作 循环运算是MATLAB中的最大弱点，在程序设计中，应当尽 量避免使用循环运算。用户可以通过将M文件向量化来优化 M文件，所谓向量化就是使用向量和矩阵运算来代替for循环 和while循环。 数据的预定义 使用for循环和while循环来增加数据结构的大小时，将影响 系统和内存的使用。提高效率的可行办法是进行预定义，,
