北京科技大学matlab和数学实验.ppt

矩阵的基本运算 解线性方程组 用数值方法计算定积分 矩阵特征值、特征向量 矩阵的基本运算 注意 k 是一个数， A 是一个矩阵 k * A A\B AX=B, X=A -1 B, A 必须是方阵 数 乘 矩阵的左除 矩阵的右除 A/B XB=A,X=AB -1 , B 必须是方阵 矩阵的行列式 det(A) A 必须为方阵 矩阵的逆 inv(A) A 必须为方阵， |A| ? 0 矩阵的乘幂 A^n A 必须为方阵， n 是正整数 矩阵行变换化简 rref(A) 求 A 阶梯形的行最简形式 矩阵的特征值、特征向量、特征多项式 [V,D]=eig(A) 例 1 A=[1,-1;2,4]; [V,D]=eig(A) ans V= -985/1393 1292/2889 985/1393 -2584/2889 方阵 A 的特 征向量矩阵 D= 2 0 0 3 方阵 A 的特 征值矩阵 矩阵的特征值、特征向量、特征多项式 p=poly(A) 若 A 为矩阵，则 p 为 A 的特征多项式系数； 若 A 为行向量，则 p 为以 A 为根的特征多项式系数。 例 1 A=[1,-1;2,4]; p=poly(A) poly2str(p,x) poly2str(p,x) 得到多项式的习惯形式 ans p=[1 -5 6] x^2-5x+6 解 线 性 方 程 组 1 、逆矩阵法（求逆法） X = 1.4000 0.4000 解： 例 1 ： 求方程组的解 2 3 4 1 x y x y ? ? ? ? ? ? ? A=[2,3;1,-1]; b=[4;1] X=inv(A)*b 相当于 2 3 4 1 1 1 x y ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ans 方程的解是： x=1.4, y=0.4 P94 A=[2,3;1,-1]; b=[4;1] X=A\b 逆矩阵法（左除与 右 除法） 例 1 ： 求方程组的解 2 3 4 1 x y x y ? ? ? ? ? ? ? 解 线 性 方 程 组 解： ans X = 1.4000 0.4000 方程的解是： x=1.4, y=0.4 相当于 AX=b,X=A\b 2 、初等变换法 解 线 性 方 程 组 在线性代数中用消元法求线性方程组的通解的过程为： 1 、 用初等变换化线性方程组为阶梯形方程组，把最 后的恒等式“0=0”去掉； 2 、如果剩下的方程当中最后的一个等式是零等于非 零的数，那么方程无解。否则有解； 3 、在有解的情况下： ? 如果阶梯形方程组中方程的个数 r 等于未知量 的个数，那么方程组有唯一的解； ? 如果阶梯形方程组中方程的个数 r 小于未知量 的个数，那么方程组有无穷多个解。 例 8 求齐次线性方程组的通解 1 2 3 4 1 2 3 4 1 2 3 4 8 10 2 0 2 4 5 0 3 8 6 2 0 x x x x x x x x x x x x ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 解： Matlab 命令为 1 0 4 0 0 1 -3/4 -1/4 0 0 0 0 ans= A=[1 -8 10 2;2 4 5 -1;3 8 6 -2]; 系数矩阵 rref(A) 行的最简形式 解 线 性 方 程 组 分析： 将 0=0 的一行去掉，则原方程组等价于 1 3 2 3 4 4 3 1 4 4 x x x x x ? ? ? ? ? ? ? ? ? 方程的个数 未知量个数 有无穷多个解 取 3 4 1 3 x x ? ? ? ? ? ? 3 4 0 4 x x ? ? ? ? ? 得 1 2 4 0 x x ? ? ? ? ? ? 取 得 1 2 0 1 x x ? ? ? ? ? 基础解系为 1 4 0 1 3 ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ， 2 0 1 0 4 ? ? ? ? ? ? ? ? ? ? ? ? ? ? 所以方程的通解为 1 2 1 2 3 4 4 0 0 1 1 0 3 4 x x k k x x ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 其中 k 1 , k 2 是任意实数 解 线