Ch3-线性代数.docVIP

83-90 第3章线性代数 应用MATLAB函数，可以求解数学问题。在这一章中，我们主要介绍MATLAB在线性代数中的应用，其内容包括：行列式的求值，矩阵的基本计算，矩阵的初等变换和矩阵的秩，矩阵的分解，求解线性方程组，向量的内积和正交，方阵的特征值和特征向量，矩阵的对角化，求解二次型的标准形，以及判别二次型的正定性。 3-1 行列式的求值 在线性代数中，行列式是一个基本工具，其应用比较广泛。在MATLAB中我们只需要借助函数det 就可以求出行列式的值，其格式为： det(A) 其中A为n阶方阵。 范例3-1 求矩阵A=的行列式的值。 程序设计： clear A=[1 0 2 1;-1 2 2 3;2 3 3 1;0 1 2 1]; det(A) 运行结果： ans = 14 程序说明： clear命令的作用是清除内存中的变量。 矩阵的输入可以有两种格式，除程序中的输入方式外，还可以如下输入： A=[1，0，2，1；-1，2，2，3；2，3，3，1；0，1，2，1] 例题分析： 在线性代数中，我们可以用det(A)的值，来判断矩阵是否可逆，即矩阵是否奇异（可逆矩阵即为非奇异矩阵）。但是，在MATLAB中，由于计算机运算存在舍入误差，用det(A)=0判断A是否可逆，并不是完全有效；如用行列式的绝对值abs(det(A))tol判断时，tol的值又不好确定。所以，我们常根据矩阵的条件数cond(A)判断矩阵是否奇异（关于矩阵的条件数可以参看3-2-11节）。 范例3-2 计算行列式 程序设计： clear syms a b c d %声明变量 A=[a 1 0 0;-1 b 1 0;0 -1 c 1;0 0 -1 d]; %生成符号矩阵 DA=det(A) 运行结果： DA = a*b*c*d+a*b+a*d+c*d+1 程序说明： 函数det也可以用于计算含有变量的行列式。 3-2 矩阵的基本运算 矩阵的基本运算主要介绍矩阵的加、减，数与矩阵相乘，矩阵与矩阵相乘，矩阵相除，矩阵的幂运算，矩阵的转置和共轭，矩阵的逆和伪逆，矩阵的迹，矩阵和向量的范数以及矩阵的条件数。矩阵的基本运算命令如表3-1所示。 表3-1矩阵的基本运算命令列表 命令 意义 命令 意义 + 矩阵相加 conj 矩阵的共轭 - 矩阵相减 inv 矩阵的逆 * 矩阵相乘 pinv 矩阵的伪逆 ＼ 矩阵左除 trace 矩阵的迹 / 矩阵右除 norm 矩阵的范数 ^ 矩阵的幂 norm 向量的范数 ’ 矩阵的转置 cond 矩阵的条件数 说明：表中仅列出具有代表性的部分命令，关于该部分内容的命令的格式及使用在下面的章节中都有介绍。 3-2-1 矩阵的加、减 进行加减运算的矩阵，要求维数相同，即行数和列数都分别相同。进行加、减运算时，矩阵相应位置的元素相加、减。 范例3-3 求矩阵A=与矩阵B=的和与差。 程序设计： clear A=[1 2 3;2 1 2;3 3 1]; B=[3 2 4;2 5 3;2 3 1]; C=A+B; D=A-B; C,D 运行结果： C = 4 4 7 4 6 5 5 6 2 D = -2 0 -1 0 -4 -1 1 0 0 例题分析： 进行加减运算的矩阵必须是同型的。 在进行矩阵相加的运算时，A+B和B+A的值相同，满足加法交换率。 3-2-2数与矩阵相乘 数与矩阵相乘，是数与矩阵中的每个元素相乘。 范例3-4 求矩阵A=与5的乘积。 程序设计： clear A=[1 0 1;2 1 1;1 2 1]; B=5*A C=A*5 运行结果： B = 5 0 5 10 5 5 5 10 5 C = 5 0 5 10 5 5 5 10 5 程序说明： 5*A与A*5的值相同。 3-2-3 矩阵与矩阵相乘 设是A一个m×n的矩阵，B是一个n×s矩阵， A= B= 则A与B的乘积C（[]）是m×n的矩阵，且 在线性代数中，我们就是按照上面的公式求得矩阵与矩阵的乘积，同时，要注意：两矩阵相乘时，第一个矩阵（左矩阵）的列数必须等于第二个矩阵（右矩阵）的行数。 范例3-5 求矩阵A=与矩阵B=的乘积。 程序设计： cl