MATLAB简介.ppt

《MATLAB简介.ppt》由会员分享，可在线阅读，更多相关《MATLAB简介.ppt（39页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、Matlab简介简介西南交通大学峨眉校区基础课部数学教研室西南交通大学峨眉校区基础课部数学教研室1MATLABMATLAB的发展史的发展史的发展史的发展史 MATLAB名字由名字由MATrix和和 LABoratory 两词的前三个字母两词的前三个字母组合而成。那是组合而成。那是20世纪七十年代，时任美国新墨西哥大学计世纪七十年代，时任美国新墨西哥大学计算机科学系主任的算机科学系主任的Cleve Moler出于减轻学生编程负担的动机，出于减轻学生编程负担的动机，为学生设计了一组调用为学生设计了一组调用LINPACK和和EISPACK矩阵软件工具矩阵软件工具包库程序的的包库程序的的“通俗易用通俗

2、易用”的接口，此即用的接口，此即用FORTRAN编写编写的萌芽状态的的萌芽状态的MATLAB。1984年由年由Little、Moler、Steve Bangert合作成立合作成立MathWorks公司，并把公司，并把MATLAB正式推向市场。从这时起，正式推向市场。从这时起，MATLAB的的内核采用内核采用C语言编写，而且除原有的数值计算能力外，还新增语言编写，而且除原有的数值计算能力外，还新增了数据图视功能。了数据图视功能。1997年仲春，年仲春，MATLAB5.0版问世，紧接着是版问世，紧接着是5.1、5.2，以及，以及和和1999年春的年春的5.3版。现今的版。现今的MATLAB拥有更丰

3、富的数据类型拥有更丰富的数据类型和结构、更友善的面向对象、更加快速精良的图形可视、更广和结构、更友善的面向对象、更加快速精良的图形可视、更广博的数学和数据分析资源、更多的应用开发工具。博的数学和数据分析资源、更多的应用开发工具。前一页前一页后一页后一页返回返回2MATLABMATLAB语言的主要特点语言的主要特点语言的主要特点语言的主要特点（1）具有丰富的数学功能）具有丰富的数学功能 包括矩阵各种运算。如：正交变换、三角分解、特征值、包括矩阵各种运算。如：正交变换、三角分解、特征值、常见的特殊矩阵等。常见的特殊矩阵等。包括各种特殊函数。如：贝塞尔函数、勒让德函数、伽码包括各种特殊函数。如：贝塞

4、尔函数、勒让德函数、伽码函数、贝塔函数、椭圆函数等。函数、贝塔函数、椭圆函数等。包括各种数学运算功能。如：数值微分、数值积分、插值、包括各种数学运算功能。如：数值微分、数值积分、插值、求极值、方程求根、求极值、方程求根、FFT、常微分方程的数值解等。常微分方程的数值解等。（2）具有很好的图视系统）具有很好的图视系统 可方便地画出两维和三维图形。可方便地画出两维和三维图形。图形用户界面图形用户界面GUI制作工具，可以制作用户菜单和控件。制作工具，可以制作用户菜单和控件。使用者可以根据自己的需求编写出满意的图形界面。使用者可以根据自己的需求编写出满意的图形界面。高级图形处理。如：色彩控制、句柄图形

5、、动画等。高级图形处理。如：色彩控制、句柄图形、动画等。前一页前一页后一页后一页3MATLABMATLAB语言的主要特点语言的主要特点语言的主要特点语言的主要特点（续）续）续）续）（3）可以直接处理声言和图形文件）可以直接处理声言和图形文件。（4）具有若干功能强大的应用工具箱。）具有若干功能强大的应用工具箱。（5）使用方便，具有很好的扩张功能。）使用方便，具有很好的扩张功能。声言文件。如：声言文件。如：WAV文件（例：文件（例：wavreadwavread，soundsound等）。等）。图形文件。如：图形文件。如：bmp、gif、pcx、tif 、jpeg等文件。等文件。如：如：SIMULI

6、NK、COMM、DSP、SIGNAL等等16种工具箱。种工具箱。（6）具有很好的帮助功能）具有很好的帮助功能 可以可以M文件转变为独立于平台的文件转变为独立于平台的EXE可执行文件。可执行文件。使用使用MATLAB语言编写的程序可以直接运行，无需编译。语言编写的程序可以直接运行，无需编译。提供十分详细的帮助文件（提供十分详细的帮助文件（PDF、HTML、demo文件）。文件）。联机查询指令：联机查询指令：help指令指令（例：（例：help elfun，help exp，help simulink），），lookfor关键词（例：关键词（例：lookfor fourier）。）。MATLAB的

7、应用接口程序的应用接口程序API是是MATLAB提供的十分重要的提供的十分重要的组件组件，由，由 一系列接口一系列接口指令指令组成组成。用户就可在。用户就可在FORTRAN或或C中中，把把MATLAB当作计算引擎使用当作计算引擎使用。前一页前一页后一页后一页4MATLAB简介简介1 MATLAB的计算功能的计算功能（1）MATLAB的数据类型的数据类型 现有四种基本数据类型：双精度数组、字符串数组、元胞现有四种基本数据类型：双精度数组、字符串数组、元胞数组、构架数组。数组、构架数组。元胞数组（元胞数组（Cell Array）如同银行里的保险箱库一样。如同银行里的保险箱库一样。该数组的基本组分是

8、元胞（该数组的基本组分是元胞（Cell），），以下标来区分。以下标来区分。元胞可以存放任何类型、任何大小的数组。元胞可以存放任何类型、任何大小的数组。构架数组（构架数组（Structure Array）也能存放各类数据。也能存放各类数据。该数组的基本组分是构架（该数组的基本组分是构架（Structure），），以下标来区分。以下标来区分。构架必须在划分构架必须在划分“域域”后才能使用。后才能使用。数据不能存放于构架，只能存放在域中。数据不能存放于构架，只能存放在域中。构架的域可以存放任何类型、任何大小的数组。构架的域可以存放任何类型、任何大小的数组。不同构架的同名域中存放的内容可不同。（不同构

9、架的同名域中存放的内容可不同。（例例 1_2_2）同一个元胞数组中各元胞的内容可以不同。（同一个元胞数组中各元胞的内容可以不同。（例例 1_2_1）前一页前一页后一页后一页5一.矩阵：1.矩阵的建立与表示法：在命令窗口中输入：A=1,2,3;4,5,6;7,8,9 可以得到：A=1 2 3 4 5 6 7 8 9若要显示整行或整列，则可以用(:)冒号来表示。冒(:)代表矩阵中行(ROWS)或列(COLUMNS)的全部。例如执行命令：A(:,2),就会显示第2列的全部，结果为：ans=2 5 8前一页前一页后一页后一页6l其他特殊矩阵的生成方法：1）、eye(m,n)或eye(m)产生m*n 或

10、 m*m的单位矩阵。例如：eye(3,4)与eye(3)分别产生如下矩阵：1 0 0 0 1 0 0 0 1 0 0 0 1 0 0 0 1 0 0 0 12）、zeros(m,n)或 zeros(m)产生m*n 或m*m 的零矩阵。例如：zeros(3,4)与zeros(3)分别产生如下矩阵：0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0前一页前一页后一页后一页73）、ones(m,n)或ones(m)产生m*n或m*m的全部元素为1的矩 阵。例如：ones(3,4)与ones(3)分别产生如下矩阵：1 1 1 1 1 1 1 1 1 1 1 1 1 1

11、1 1 1 1 1 1 12.常用矩阵函数：1）、d=eig(A)返回矩阵A的特征值所组成的列向量；v,d=eig(A)返回分别由矩阵A的特征向量和特征值（以其为主对角线元素，其余元素为零）的两个矩阵。例如执行命令：v,d=eig(A)结果为：前一页前一页后一页后一页8 v=d=0.2320 0.7858 0.4082 16.1168 0 0 0.5253 0.0868 -0.8165 0 -1.1168 0 0.8187 -0.6123 0.4082 0 0 -0.0000 其中v(:,i)为d(i,i)所对应的特征向量。2）、det(A)计算行列式A的值。例如：det(A)结果为：ans=

12、0前一页前一页后一页后一页9 3）、expm(A)对矩阵A求幂。例如：expm(A)结果为：ans=1.0e+006*1.1189 1.3748 1.6307 2.5339 3.1134 3.6929 3.9489 4.8520 5.75524）、inv(A)求矩阵A的逆。例如：inv(A)结果为：Warning:Matrix is close to singular or badly scaled.Results may be inaccurate.RCOND=2.055969e-018.ans=1.0e+016*-0.4504 0.9007 -0.4504 0.9007 -1.8014 0

13、.9007 -0.4504 0.9007 -0.4504前一页前一页后一页后一页105）、orth(A)返回对应于A的正交化矩阵。例如：orth(A)结果为：ans=0.2148 0.8872 0.5206 0.2496 0.8263 -0.38796）、poly(A)若A为一矩阵，则返回A的特征多项式。例如：poly(A)结果为：ans=1.0000 -15.0000 -18.0000 -0.0000 若A为一向量，则返回以A的元素为根的特征多项式。例如：r=1,2,3;p=poly(r)结果为：p=1 -6 11 -67）、rank(A)计算矩阵A的秩。例如：r=rank(A)结果为：r=

14、2 前一页前一页后一页后一页113.矩阵的四则运算符号：加 “+”减 “”乘 “*”除 “/”共轭转置 “”非共轭转置“.”例如：b=1+2i;3+4ib=1.0000+2.0000i 3.0000+4.0000ibans=1.0000-2.0000i 3.0000-4.0000ib.ans=1.0000+2.0000i 3.0000+4.0000i前一页前一页后一页后一页124.矩阵分解：1）、q,r=qr(A)将矩阵A做正交化分解，使得A=q*r。q为单位矩阵(unitary matrix),其范数（norm）为1。r为对角化的上三角矩阵。例如：q,r=qr(A)q=-0.1231 0.9

15、045 0.4082 -0.4924 0.3015 -0.8165 -0.8616 -0.3015 0.4082 r=-8.1240 -9.6011 -11.0782 0 0.9045 1.8091 0 0 -0.0000 norm(q)ans=1.0000前一页前一页后一页后一页132）、L,U=lu(A)将矩阵A做对角线分解，使得A=L*U，L为下三角矩阵（lower triangular matrix），U为上三角矩阵(upper triangular matrix)。例如：L,U=lu(A)L=0.1429 1.0000 0 0.5714 0.5000 1.0000 1.0000 0

16、0 U=7.0000 8.0000 9.0000 0 0.8571 1.7143 0 0 0.0000前一页前一页后一页后一页14二、多项式：多项式是用向量形式来表示，从最右边算起，第一个为0阶系数，第二个为1阶系数，依次类推。例如一个一元三次多项式：4x3+3x2+2x+1 用4 3 2 1来表示。1.多项式的运算：1）、乘：conv指令执行多项式的相乘运算，指令格式为：z=conv(x,y)例如：x=1 3 5;y=2 4 6;z=conv(x,y)z=2 10 28 38 30前一页前一页后一页后一页15如果要对两个以上的多项式进行相乘，可以重复使用conv指令，例如：（x,y同上）co

17、nv(conv(x,y),x)ans=2 16 68 172 284 280 1502）、分解：与1）相反，用deconv指令，其指令格式为：z,r=deconv(x,y)表示x除以 y商为z，余数为r。例如：z,r=deconv(z,x)z=2 4 6 r=0 0 0 0 0前一页前一页后一页后一页163）、求根：roots指令用于求多项式的根。例如：fx=1 3 2;rootoffx=roots(fx)rootoffx=-2 -1 4）、polyval(p,x)计算多项式p在x出的值，其中x可以是点或向量或矩阵。例如：p=1 -6 11 -6；x=1;p1=polyval(p,x)结果为：

18、p1=0 x=1,2,3;p2=polyval(p,x)结果为：p2=0 0 0前一页前一页后一页后一页17 x=A;p3=polyval(p,x)结果为：p3=0 0 0 6 24 60 120 210 336 5）、polyder(p)求p的微分多项式。例如：p=1-6 11-6;dp=polyder(p)dp=3 -12 11前一页前一页后一页后一页18 6）、r,p,k=residue(x,y)求x/y的部分因式分解。若多项式x,y都没有重根,则可把x/y的比值表示为x/y=r1/(s-p1)+r2/(s-p2)+.+rn/(s-pn)+ks例如 用residue指令求x/(x2+3x

19、+2)的部分因式分解：x=1 0;y=1 3 2;r,p,k=residue(x,y)r=2 -1 p=-2 -1 k=前一页前一页后一页后一页19l当输入三个参数 r,p,k 时，该函数又会生成原来的函数。例如：?x,y=residue(r,p,k)x=1 0 y=1 3 2前一页前一页后一页后一页20三 符号变量、符号表达式、抽象函数：函数sym用于生成符号变量和符号表达式，如：x=sym(x)a=sym(alpha)分别创建变量x,alpha f=sym(a*x2+b*x+c)创建变量表达式f，但要注意此式并没有自动创建变量a,b,c,x。可以用函数syms对多个变量同时定义，如：sym

20、s a b c x 函数sym也可以用来表示确定的函数，如 f=sym(f(x)生成函数f(x)。前一页前一页后一页后一页21四 常见符号计算：1.微分：diff是求微分最常用的函数。其输入参数既可以是函数表达式，也可以是符号矩阵。Diff(f,x,n)表示对f关于x求n阶导数。例如：1).下面程序段将生成表达式sin(ax),并分别对其中的x和a求导。?syms a x?f=sin(a*x);?df=diff(f,x)df=cos(a*x)*a?dfa=diff(f,a)dfa=cos(a*x)*x前一页前一页后一页后一页222）、若输入参数为矩阵，将对矩阵中的每个元素求导。?syms a

21、x?A=-sin(a*x),sin(a*x);cos(a*x),cos(a*x)A=-sin(a*x),sin(a*x)cos(a*x),cos(a*x)?dy=diff(A,x)dy=-cos(a*x)*a,cos(a*x)*a -sin(a*x)*a,-sin(a*x)*a前一页前一页后一页后一页233）、可用函数jacobian来计算Jacobi矩阵。?syms r l f?x=r*cos(l)*cos(f);?y=r*cos(l)*sin(f);?z=r*sin(l);?J=jacobian(x;y;z,r l f)J=cos(l)*cos(f),-r*sin(l)*cos(f),-r

22、*cos(l)*sin(f)cos(l)*sin(f),-r*sin(l)*sin(f),r*cos(l)*cos(f)sin(l),r*cos(l),0 前一页前一页后一页后一页242.积分：用函数int来求符号表达式的积分。命令格式为：int(f,r,x0,x1)其中f为所要积分的表达式，r为积分变量，若为定积分，则x0,x1为积分上下限。例：?sym x;?sym k real?f=exp(-(k*x)2)f=exp(-k2*x2)?int(f,x,-inf,inf)ans=signum(k)/k*pi(1/2)前一页前一页后一页后一页253.级数求和：函数用于对符号表达式求和。例：?s

23、yms k;?s1=symsum(1/k2,1,inf)s1=1/6*pi2 前一页前一页后一页后一页264.极限：用函数limit来求表达式的极限。函数limit的常用调用格式：l 数学表达式 命令格式 l Limit(f),或limit(f,x)l Limit(f,x,a),或 limit(f,a)l Limit(f,x,a,left)l Limit(f,x,a,right)前一页前一页后一页后一页275.化简：1）、collect(f)将表达式中相同次幂的项合并，也可以再输入一个参数指定以哪个变量的幂次合并。2）、expand(f)将表达式展开。3）、horner(f)将表达式转换为嵌套

24、格式。4）、factor(f)将表达式分解因式，并且分解后的多项式的所有系数都为有理数。5）、simplify(f)利用函数规则对表达式进行化简。前一页前一页后一页后一页28前一页前一页后一页后一页29MATLAB的的程序设计程序设计 前一页前一页后一页后一页301 脚本文件和函数文件脚本文件和函数文件11 M脚本文件脚本文件 对于一些比较简单的问题对于一些比较简单的问题，在指令窗中直接输入指令计算，在指令窗中直接输入指令计算。对于复杂计算，采用脚本文件（对于复杂计算，采用脚本文件（Script file）最为合适最为合适。MATLAB只是按文件所写的指令执行只是按文件所写的指令执行。M脚本文

25、件的特点是：脚本文件的特点是：脚本文件的构成比较简单，只是一串按用户意图排列而成的脚本文件的构成比较简单，只是一串按用户意图排列而成的（包括控制流向指令在内的）（包括控制流向指令在内的）MATLAB指令集合。指令集合。脚本文件运行后脚本文件运行后，所产生的所有变量都驻留在，所产生的所有变量都驻留在 MATLAB基基本工作空间本工作空间（Base workspace）中中。只要用户不使用清除指只要用户不使用清除指令（令（clear），），MATLAB指令窗不关闭，这些变量将一直保指令窗不关闭，这些变量将一直保存在存在基本工作空间基本工作空间中中。M文件有两种形式文件有两种形式：脚本文件（脚本文件

26、（Script File）和函数文件和函数文件（Function File）。）。这两种文件的扩展名，均为这两种文件的扩展名，均为“.m”。前一页前一页后一页后一页311 脚本文件和函数文件（脚本文件和函数文件（续续续续1 1）12 M函数文件函数文件 与脚本文件不同与脚本文件不同，函数文件犹如一个，函数文件犹如一个“黑箱黑箱”，把一些数据，把一些数据送进并经加工处理，再把结果送出来。送进并经加工处理，再把结果送出来。MATLAB提供的函数指令大部分都是由函数文件定义的。提供的函数指令大部分都是由函数文件定义的。M函数文件的特点是：函数文件的特点是：从形式上看从形式上看，与脚本文件不同，与脚本

27、文件不同，函数文件的笫一行总是以，函数文件的笫一行总是以“function”引导的引导的“函数申明行函数申明行”。从运行上看从运行上看，与脚本文件运行不同，与脚本文件运行不同，每当函数文件运行，每当函数文件运行，MATLAB就会专门为它开辟一个临时工作空间，称为就会专门为它开辟一个临时工作空间，称为函数函数工作空间工作空间（Function workspace）。当执行文件最后一条。当执行文件最后一条指令时指令时，就结束该函数文件的运行，同时该临时函数空间，就结束该函数文件的运行，同时该临时函数空间及其所有的中间变量就立即被清除。及其所有的中间变量就立即被清除。MATLAB允许使用比允许使用比

28、“标称数目标称数目”较少的输入输出宗量，较少的输入输出宗量，实现对函数的调用实现对函数的调用。前一页前一页后一页后一页321 脚本文件和函数文件（脚本文件和函数文件（续续续续2 2）13 M文件的一般结构文件的一般结构 由于从结构上看由于从结构上看，脚本文件只是比函数文件少一个，脚本文件只是比函数文件少一个“函数申函数申明行明行”，所以只须描述清楚函数文件的结构，所以只须描述清楚函数文件的结构。典型典型 M函数文件的结构如下函数文件的结构如下：函数申明行：位于函数文件的首行，以关键字函数申明行：位于函数文件的首行，以关键字 functio 开头，开头，函数名以及函数的输入输出宗量都在这一行被定

29、义。函数名以及函数的输入输出宗量都在这一行被定义。笫一注释行：紧随函数申明行之后以笫一注释行：紧随函数申明行之后以%开头笫一注释行。该开头笫一注释行。该行供行供lookfor关键词查询和关键词查询和 help在线帮助使用在线帮助使用。在线帮助文本区在线帮助文本区：笫一注释行及其之后的连续以：笫一注释行及其之后的连续以%开头的所开头的所有注释行构成整个在线帮助文本。有注释行构成整个在线帮助文本。编写和修改记录：与在线帮助文本区相隔一个编写和修改记录：与在线帮助文本区相隔一个“空空”行，也行，也以以%开头，标志编写及修改该开头，标志编写及修改该M文件的作者和日期等文件的作者和日期等。函数体：为清晰

30、起见，它与前面的注释以函数体：为清晰起见，它与前面的注释以“空空”行相隔。行相隔。前一页前一页后一页后一页332 函数调用和参数传递函数调用和参数传递21 局部变量和全局变量局部变量和全局变量 局部（局部（Local）变量：它存在于函数空间内部的中间变量，变量：它存在于函数空间内部的中间变量，产生于该函数的运行过程中，其影响范围也仅限于该函数本产生于该函数的运行过程中，其影响范围也仅限于该函数本身身。全局（全局（Global）变量：通过变量：通过 global 指令，指令，MATLAB也允许几也允许几个不同的函数空间以及基本工作空间共享同一个变量，这种被个不同的函数空间以及基本工作空间共享同一

31、个变量，这种被共享的变量称为全局变量。共享的变量称为全局变量。22 函数调用函数调用 在在MATLAB中，调用函数的常用形式是：中，调用函数的常用形式是：输出参数输出参数1,输出参数输出参数2,=函数名函数名(输入参数输入参数1,输入参数输入参数2,)函数调用可以嵌套，一个函数可以调用别的函数，甚至调函数调用可以嵌套，一个函数可以调用别的函数，甚至调用它自己用它自己（递归调用）。（递归调用）。前一页前一页后一页后一页342 函数调用和参数传递（函数调用和参数传递（续）续）续）续）23 参数传递参数传递 MATLAB在函数调用上有一个与众不同之处在函数调用上有一个与众不同之处：函数所传递的：函数

32、所传递的参数具有可调性参数具有可调性。传递参数数目的可调性来源于如下两个传递参数数目的可调性来源于如下两个MATLAB永久变量：永久变量：函数体内的函数体内的 nargin 给出调用该函数时的输入参数数目。给出调用该函数时的输入参数数目。函数体内的函数体内的 nargout 给出调用该函数时的输出参数数目。给出调用该函数时的输出参数数目。只要在函数文件中包括这两个变量，就可以知道该函数文件只要在函数文件中包括这两个变量，就可以知道该函数文件调用时的输入参数和输出参数数目。调用时的输入参数和输出参数数目。值得注意：值得注意：nargin、nargout 本身都是函数本身都是函数，不是变量，所不是

33、变量，所以用户不能赋值，也不能显示。以用户不能赋值，也不能显示。“变长度变长度”输入输出宗量：输入输出宗量：varargin、varrgout。具有接受具有接受“任意多输入任意多输入”、返回、返回“任意多输出任意多输出”的能力的能力。跨空间变量传递：跨空间变量传递：evalin。（参考（参考：circle.m，am1.m）前一页前一页后一页后一页3523 MATLAB的程序结构和控制流的程序结构和控制流231 程序结构程序结构 循环结构：循环结构：MATLAB提供两种循环方式。提供两种循环方式。顺序结构顺序结构 分支结构：分支结构：ifelseend。forend 循环和循环和while-en

34、d循环。循环。232 程序流控制程序流控制 常用指令：常用指令：return，echo，input，pause，keyboard，break。switch-case 结构。结构。try-catch 结构。结构。警示指令：警示指令：error，warning。前一页前一页后一页后一页3623 MATLAB的程序结构和控制流（的程序结构和控制流（续）续）续）续）233 图形用户界面（图形用户界面（GUI）编程编程 现代的主流应用程序已经从命令行的交互方式转变为以图现代的主流应用程序已经从命令行的交互方式转变为以图形界面为主的交互方式，这主要是由于它给用户带来了操作形界面为主的交互方式，这主要是由于

35、它给用户带来了操作和控制的方便与灵活性。和控制的方便与灵活性。（面向对象编程）（面向对象编程）MATLAB能够以比较简单的方式实现一系列的图形界面功能够以比较简单的方式实现一系列的图形界面功能。通过对控件、菜单属性的设置和能。通过对控件、菜单属性的设置和 Callback 的编写，就能的编写，就能够满足大多数用户的需求。够满足大多数用户的需求。控件的控件的 Callback 属性属性：Callback 属性的取值是字符串，可属性的取值是字符串，可以是某个以是某个M文件名或一小段文件名或一小段MATLAB语句语句。当用户激活控件当用户激活控件对象（对象（例例如如：在控件对象图标上单击鼠标左键：在

36、控件对象图标上单击鼠标左键）时，应用程）时，应用程序就运行序就运行 Callback 属性定义的子程序。属性定义的子程序。菜单的菜单的 Callback 属性属性：Callback 属性的取值是字符串，可属性的取值是字符串，可以是某个以是某个M文件名或一小段文件名或一小段MATLAB语句语句。当用户激活菜单当用户激活菜单对象时，若没有子菜单就运行对象时，若没有子菜单就运行 Callback 属性定义的子程序。属性定义的子程序。若有，先运行若有，先运行 Callback 属性定义的子程序，再显示子菜单。属性定义的子程序，再显示子菜单。前一页前一页后一页后一页3724 M文件的调试文件的调试 编写

37、编写 M文件时，错误（文件时，错误（Bug）在所难免。错误有两种：语法在所难免。错误有两种：语法（Syntax）错误和运行（错误和运行（Run-time）错误。错误。语法错误是指变量名、函数名的误写，标点符号的缺、漏等。语法错误是指变量名、函数名的误写，标点符号的缺、漏等。对于这类错误，通常能在运行时发现，终止执行，并给出相应对于这类错误，通常能在运行时发现，终止执行，并给出相应的错误原因以及所在行号。的错误原因以及所在行号。运行错误是算法本身引起的，发生在运行过程中。相对语法运行错误是算法本身引起的，发生在运行过程中。相对语法错误而言，运行错误较难处理错误而言，运行错误较难处理。尤其是。尤其

38、是M函数文件，它一旦运函数文件，它一旦运行停止，其中间变量被删除一空，错误很难查找。行停止，其中间变量被删除一空，错误很难查找。有两种调试方法：直接调试法和工具调试法。有两种调试方法：直接调试法和工具调试法。前一页前一页后一页后一页3824 M文件的调试（文件的调试（续续续续1 1）直接调试法：可以用下面方法发现某些运行错误。直接调试法：可以用下面方法发现某些运行错误。在在M文件中，将某些语句后面的分号去掉，文件中，将某些语句后面的分号去掉，迫使迫使M文件输文件输出一些中间计算结果，以便发现可能的错误。出一些中间计算结果，以便发现可能的错误。在适当的位置，添加显示某些关键变量值的语句（包括使在

39、适当的位置，添加显示某些关键变量值的语句（包括使用用 disp 在内）。在内）。利用利用 echo 指令，使运行时在屏幕上逐行显示文件内容。指令，使运行时在屏幕上逐行显示文件内容。echo on 能显示能显示M脚本文件；脚本文件；echo FunNsme on 能显示名为能显示名为FunNsme 的的M函数文件。函数文件。在原在原M脚本或函数文件的适当位置，增添指令脚本或函数文件的适当位置，增添指令 keyboard。keyboard 语句可以设置程序的断点语句可以设置程序的断点。通过将原通过将原M函数文件的函数申明行注释掉，可使一个中间函数文件的函数申明行注释掉，可使一个中间变量难于观察的变量难于观察的M函数文件变为一个所有变量都保留在基函数文件变为一个所有变量都保留在基本工作空间中的本工作空间中的M脚本文件。脚本文件。前一页前一页后一页后一页39