MATLAB与数学实验讲义

目录第一部分MATLAB 入门 (1)第一讲变量与函数 (1)一、变量与函数 (1)第二讲数组与矩阵 (3)一、数组 (3)二、矩阵 (4)第三讲M文件与程序设计初步 (6)一、M文件的特点与形式 (6)二. 命令式M文件 (7)三. 函数式M文件 (7)四、控制流 (8)第四讲图形处理功能 (10)一、二维图形 (11)二、三维图形 (18)第二部分Matlab与数学实验 (19)第五讲插值 (19)一、一维插值 (19)二、二维插值 (25)第六讲拟合 (27)一、线性最小二乘拟合 (27)二、非线性最小二乘拟合 (29)二、非线性最小二乘拟合 (33)第七讲数值积分与数值微分 (35)一数值积分 (35)二、数值微分 (39)第八讲常微分方程数值解 (40)一、微分方程的解析解 (40)二、常微分方程数值解 (41)三、用Matlab软件求常微分方程的数值解 (42)四、实例1：导弹追踪问题 (43)第九讲线性方程组的数值解法 (47)一、直接解法 (47)二、迭代解法 (48)第一部分MATLAB 入门第一讲变量与函数1.MATLAB作为线性系统的一种分析和仿真工具，是理工科大学生应该掌握的技术工具，它作为一种编程语言和可视化工具，可解决工程、科学计算和数学学科中许多问题. 2.MATLAB建立在向量、数组和矩阵的基础上，使用方便，人机界面直观，输出结果可视化。

3.矩阵是MA TLAB的核心4.MATLAB的进入与运行方式（两种）一、变量与函数1、变量MATLAB中变量的命名规则是：（1）变量名必须是不含空格的单个词；（2）变量名区分大小写；（3）变量名最多不超过19个字符；（4）变量名必须以字母打头，之后可以是任意字母、数字或下划线，变量名中不允许使用标点符号.特殊变量表2、数学运算符号及标点符号（1）MATLAB的每条命令后，若为逗号或无标点符号，则显示命令的结果；若命令后为分号，则禁止显示结果. （2）“%” 后面所有文字为注释.（3）“...”表示续行.3、数学函数第二讲数组与矩阵一、数组1、创建简单的数组x=[a b c d e f ]创建包含指定元素的行向量x=first：last创建从first开始，加1计数，到last结束的行向量x=first：increment：last创建从first开始，加increment计数，last结束的行向量x=linspace(first，last，n）创建从first开始，到last结束，有n个元素的行向量x=logspace(first，last，n）创建从first开始，到last结束，有n个元素的对数分隔行向量.2、数组元素的访问（1）访问一个元素：x(i)表示访问数组x的第i个元素.（2）访问一块元素：x(a ：b ：c)表示访问数组x的从第a个元素开始，以步长为b到第c个元素（但不超过c），b可以为负数，b缺损时为1.（3）直接使用元素编址序号. x([a b c d]) 表示提取数组x的第a、b、c、d个元素构成一个新的数组[x(a) x(b) x(c) x(d)].3、数组的方向前面例子中的数组都是一行数列，是行方向分布的. 称之为行向量. 数组也可以是列向量，它的数组操作和运算与行向量是一样的，唯一的区别是结果以列形式显示.产生列向量有两种方法：直接产生例c=[1；2；3；4]转置产生例b=[1 2 3 4]; c=b‟说明：以空格或逗号分隔的元素指定的是不同列的元素，而以分号分隔的元素指定了不同行的元素.4、数组的运算（1）标量-数组运算数组对标量的加、减、乘、除、乘方是数组的每个元素对该标量施加相应的加、减、乘、除、乘方运算.设：a=[a1,a2,…,an], c=标量则：a+c=[a1+c,a2+c,…,an+c]a.*c=[a1*c,a2*c,…,an*c]a./c= [a1/c,a2/c,…,an/c](右除）a.\c= [c/a1,c/a2,…,c/an] (左除）a.^c= [a1^c,a2^c,…,an^c]c.^a= [c^a1,c^a2,…,c^an]当两个数组有相同维数时，加、减、乘、除、幂运算可按元素对元素方式进行的，不同大小或维数的数组是不能进行运算的.设：a=[a1,a2,…,an], b=[b1,b2,…,bn]则：a+b= [a1+b1,a2+b2,…,an+bn]a.*b= [a1*b1,a2*b2,…,an*bn]a./b= [a1/b1,a2/b2,…,an/bn]a.\b=[b1/a1,b2/a2,…,bn/an]a.^b=[a1^b1,a2^b2,…,an^bn]二、矩阵1、矩阵的建立逗号或空格用于分隔某一行的元素，分号用于区分不同的行. 除了分号，在输入矩阵时，按Enter键也表示开始一新行. 输入矩阵时，严格要求所有行有相同的列.例m=[1 2 3 4 ；5 6 7 8；9 10 11 12]p=[1 1 1 12 2 2 23 3 3 3]特殊矩阵的建立：a=[ ] 产生一个空矩阵，当对一项操作无结果时，返回空矩阵，空矩阵的大小为零.b=zeros(m，n) 产生一个m行、n列的零矩阵c=ones(m，n) 产生一个m行、n列的元素全为1的矩阵d=eye(m，n) 产生一个m行、n列的单位矩阵2、矩阵中元素的操作（1）矩阵A的第r行：A（r，：）（2）矩阵A的第r列：A（：，r）（3）依次提取矩阵A的每一列，将A拉伸为一个列向量：A（：）（4）取矩阵A的第i1~i2行、第j1~j2列构成新矩阵:A(i1:i2, j1:j2)（5）以逆序提取矩阵A的第i1~i2行，构成新矩阵:A(i2:-1：i1，：）（6）以逆序提取矩阵A的第j1~j2列，构成新矩阵:A(:, j2:-1：j1）（7）删除A的第i1~i2行，构成新矩阵:A(i1:i2，：)=[ ]（8）删除A的第j1~j2列，构成新矩阵:A(：，j1:j2)=[ ]（9）将矩阵A和B拼接成新矩阵：[A B]；[A；B]3、矩阵的运算同标量-数组运算。

（2）矩阵-矩阵运算[1] 元素对元素的运算，同数组-数组运算。

[2]矩阵运算：矩阵加法：A+B矩阵乘法：A*B方阵的行列式：det（A）方阵的逆：inv（A）方阵的特征值与特征向量：[V，D]=eig[A]关系与逻辑运算1、关系操作符2、逻辑运算符第三讲M文件与程序设计初步一、M文件的特点与形式要说明M文件的特点，就得从MATLAB本身说起。

MATLAB实质上是一种解释性语言，就MATLAB(matlab.exe)本身来说，它并不能做任何事情，它就像DOS操作系统的一样，本身没有实现功能而只对用户发出的指令起到解释执行的作用。

像前面介绍过的命令行式的操作一样，命令先送到MATLAB系统内解释，再运行得到结果。

这样就给用户提供了最大的方便。

用户可以把所要实现的指令罗列编制成文件，再统一送入MATLAB系统中解释运行，这就是M文件。

只不过此文件必须以m为扩展名，MA TLAB 系统才能识别。

也就是说，M文件其实是一个像命令集一样的ASCII(纯文本)码文件。

因此M文件语法简单，调试容易，人机交互性强。

用户可以便用任何字处理软件对其进行编写和修改。

正是M文件的这个特点造就了MATLAB强大的可开发性和可扩展性，Mathworks 公司推出的一系列工具箱就是明证。

而正是有了这些工具箱，MATLAB才能被广泛地应用于信号处理、神经网络、鲁棒控制、系统识别、控制系统、实时工作系统、图形处理、光谱分析、模型预测、模糊逻辑、数字信号处理、定点设置、金融管理、小波分析、地图工具、交流通信、模型处理、LMI控制、概率统计、样条处理、工程规划、非线性控制设计、QFT控制设计、NAG等各个领域。

对个人用户来说，还可以利用M文件来建造和扩充属于自己的“库”。

因此，一个不了解M文件，没有掌握M文件的MATLAB使用者不能称其为一个真正的MATLAB用户。

由于商用的MA TLAB软件用C/C++语言编写而成。

因此，M文件的语法与C语言十分相似。

对广大的C语言爱好者来说，M文件的编写是相当容易的。

M文件有两种形式，命令式（Script）和函数式(Function).两者相同之处在于它们都是以m作为扩展名的文本文件，不进入命令窗口，而是由文本编辑器来创建的外部文本文件。

命令式M文件就是命令行的简单叠加，MATLAB会自动按顺序执行文件中的命令，与批处理文件类似，在MA TLAB命令窗口直接输入此文件的文件名，MA TLAB可逐一执行在此文件内的所有命令，和在命令窗口中逐行输入这些命令一样。

这样就解决了用户在命令窗中运行许多命令的麻烦;还可以避免用户做许多重复性的工作。

函数式M文件主要用以解决参数传递和函数调用的问题，它的第一句以function语句为引导。

另外，值得注意的是，命令式M文件在运行过程中可以调用MATLAB工作域内所有的数据，而且，所产生的所有变量均为全局变量。

也就是说，这些变量一旦生成，就一直保存在工作空间中，直到执行clear 或quit时为止。

在函数式M文件中的变量除特别声明外，均为局部变量。

因为MATLAB是一种解释性语言，所以即使在某个或某些函数中存在语法错误，但如果没执行到该语句时可能就不会发现该错误，这在一个成功的程序设计中是不能容许的。

要查出某目录中所有的M函数语法错误，首先应该用cd 命令进入该目录，然后运行pcode * 命令进行伪代码转换。

因为该命令会将MATLAB 函数转换成伪代码，而在转换过程中该程序将自动翻译每一条语句，所以一旦发现有语法错误，将会停止翻译，给出错误信息。

改正了该语法错误后，再重新执行pcode 命令，直到没有错误为止。

至少这样会保证目录下所有的程序不含有语法错误。

一个M文件首次被调用时，MATLAB将首先对该M文件进行语法分析。

并把生成的相应内部伪代码（P码）存放在内存中。

此后当再次调用该M文件时，将直接运行该文件在内存中的P码文件而不会对原码文件重复进行语法分析。

P码文件和原码文件具有相同的文件名，但其扩展名为“.p”,并且其运行速度要高于原码文件，但对于规模不大的文件，用户一般察觉不到这种速度上的优势。

P码文件不是只有当M文件被调用时才生成，也可被预先生成。

二. 命令式M文件由于命令式M文件的运行相当于在命令窗口中逐行输入并运行命令，因此，在编制此类文件时，只需把所要执行的命令按行编辑到指定的文件中，且变量不需要预先定义，也不存在文件名对应问题。

例：假设当前目录下有一个命令M文件：%solver.m%used to solve A*x=b%where A=[-1.5 1 2;3 -1 1;-1 3 5],b=[2.5;5;8].A=[-1.5 1 2;3 -1 1;-1 3 5];b=[2.5;5;8];x=A\b在命令窗口中执行solver命令，即可得到方程组的解。