第2章 数值计算功能.pptVIP

MATLAB应用相当广泛的一个领域是解决数值计算问题。这是由于MATLAB 对处理数组具有非凡的能力，使其成为许多科学与工程应用中的一个有力的工具。本章主要介绍矩阵的创建、运算以及如何运用MATLAB解决线性代数中的相关问题。 【学习目标】 能够创建矩阵，并熟练掌握矩阵的各种运算法则 掌握稀疏矩阵的基本操作。 学会运用MATLAB解决线性代数的相关问题。 掌握如何运用MATLAB实现数据的插值、曲线的拟合和零极点的求解。 2.1 矩阵及其运算 MATLAB最基本、最重要的功能就是进行矩阵运算。所有的数值功能都是以矩阵为基本单位实现的。 2.1.1 矩阵生成 1.直接生成 矩阵的生成可以通过在方括号[]中输入元素实现，不同行之间要用分号隔开，不同列用空格分开。 【例2-1】生成一个2×3矩阵。 a=[1 2 3;4 5 6] a = 1 2 3 4 5 6 可以用冒号快速建立矩阵，步长在默认的情况下为1，也可以自定义步长。 a=[1:3;4:6] a = 1 2 3 4 5 6 2.函数生成 MATLAB中提供了特殊函数可以生成特殊矩阵。 【例2-2】生成一个2×3全零矩阵，一个3×3单位矩阵。 a=zeros(2，3) %生成2×3全零阵 a = 0 0 0 0 0 0 a=eye(3，3) %生成3×3单位阵 a = 1 0 0 0 1 0 0 0 1 2.1.2 矩阵连接 矩阵连接，是联合一个或多个矩阵形成一个新矩阵。方括号[]就可以作为矩阵连接操作符。表达式A=[B C]就是水平连接矩阵B和C；A=[B；C]就是垂直连接矩阵B和C连接矩阵也可以用函数形成一个新矩阵，如cat函数。cat函数调用格式： 1) C=cat(dim，A，B) Dim指链接方向。Dim的可选值：1，垂直方向；2，水平方向；3，生成三维数组。 2) C=cat(dim，A1，A2，…) 可以实现连接多个矩阵。 【例2-3】把矩阵A，B分别在垂直方向，水平方向连接，A=[1 2 3]，B=[4 5 6]。 C1=cat(1，A，B) %垂直方向连接 C1 = 1 2 3 4 5 6 C2=cat(2，A，B) %水平方向连接 C2 = 1 2 3 4 5 6 2.1.3 矩阵元素的下标引用及操作 1.矩阵下标 在MATLAB中，矩阵元素的引用可以通过下标表示，通常有如下规则： 1) 用两个下标来表示：在矩阵A中，第i行j列的元素用A(i，j)来表示。 2) 用一个下标来表示：对于矩阵，由于MATLAB的基本运算都是对列操作的，矩阵认为是按列优先排列的一个长的列向量，从而可用单下标引用。 【例2-4】提取矩阵的元素。 A = [1 2 3; 4 5 6; 7 8 9] A = 1 2 3 4 5 6 7 8 9 可以用A(m,n)选出第m行n列的元素。 A(2,3) ans = 6 也可以用一个下标选出元素。 A(8) ans = 6 2.矩阵的抽取 用下标可以从大矩阵里面抽取某些元素组成小的矩阵，并使用()运算符。 【例2-5】提取矩阵A的前2行为一个新矩阵B。 B=A(1:2,:) % B=A([1 2],:) B = 1 2 3 4 5 6 3.删除矩阵的行和列 我们把空的中括号[]称为空矩阵，即矩阵大小为0×0。可以通过把矩阵的行或列赋值为空矩阵来删除某行或某列。 【例2-6】删除上例中A的第1，2列。 A(:,[1 2])=[] A = 3 6 9 在MATLAB里不能用双下标来删除矩阵的某个元素，例如X(1，2)=[]将给出错误信息。但可以用单下标来删除矩阵的某个或某些元素，MATLAB删除这些元素后，将剩下的元素按列编成一个长的行向量。 4.reshape函数实现改变矩阵大小 reshape函数可以实现矩阵行列大小的变化，但变化前矩阵行列值的乘积和变化后的行列值的乘积不变。 reshape(A,M,N)命令将矩阵A的所有元素分配到一个M×N的新的矩阵，当矩阵A的元