第四讲应用MATLAB解决高等代数问题幻灯片.ppt

得简化的行阶梯形矩阵为 ans = 1 0 -11/9 0 0 1 5/9 0 0 0 0 1 0 0 0 0 最简矩阵中的有三个不全为零的行向量，所以向量组的秩为3，显然第一列、第二列、第四列线性无关，所以对应于原向量一个极大无关组为a1,a2 a4，最简矩阵中第三列向量有两个非零元素-11/9，5/9，它们是方程组 的解（x1=-11/9,x2=5/9)，也是方程组 的解，所以 线性表出 被最大无关组 2）输入向量及命令如下： a1=[1;-1;2;4];a2=[0;3;1;2];a3=[3;0;7;14]; a4=[1;-1;2;0];a5=[2;1;5;6]; A=[a1 a2 a3 a4 a5] rref(A) 得最简行阶梯形矩阵 ans = 1 0 3 0 1 0 1 1 0 1 0 0 0 1 1 0 0 0 0 0 由此可知向量组的秩为3，第1列，第2列，第4列的向量是线性无关的，所以a1,a2,a4是极大无关组。最简矩阵中第三列向量有两个非零元素3,1，它们是方程组 的解（x1=3,x2=1)，也是方程组 的解，所以 线性表出 被最大无关组 最简矩阵中第五列向量有三个非零元素1,1,1，它们是方程组 的解(x1=1,x2=1,x4=1)，也是方程组 的解，所以 线性表出 被最大无关组 注意：两个例子中输入的向量和命令有所不同，请同学们思考为什么？ 五.矩阵的特征值和特征向量 1 矩阵的特征值和特征向量 设A是n阶方阵，k是一个数，如果存在一非零的列向量X使得AX=kX成立，则称数k为A的征值，非零列向量X称为方阵A的属于特征值K的一个特征向量。 用MATLAB的命令 eig可以求出矩阵A的特征值和特征向量的方法有两种 法一）只求A的特征值命令为eig(A) 法二）同时求出特征值和特征向量用命令[p d]=eig(A) 例　求方阵 特征值和特征向量。 解：先输入矩阵的数据，然后用eig的两种使用方法求解，命令如下 A=[3 0 4;0 6 0;4 0 3]; eig(A) [p d]=eig(A) 第一个命令eig(A)的结果为 ans = -1 6 7 p = 0.7071 0 0.7071 0 -1.0000 0 -0.7071 0 0.7071 d = -1 0 0 0 6 0 0 0 7 命令［p d]=eig(A)计算结果为 北京大学《高等代数》Ｐ２８４ 求矩阵 应用eig(A)得 ans = -1 -1 5 应用[p d]=eig(A)结果为 p= 2131/3543 709/1284 780/1351 408/2299 -369/463 780/1351 -747/959 294/1201 780/1351 d = -1 0 0 0 -1 0 0 0 5 ２ 矩阵的相似对角化 设三阶矩阵Ａ有三个线性无关的特征向量a1,a2,a3，对应的特征值为k1,k2,k3,现定义两个矩阵 Ｐ＝［a1 a2 a3] Aa1=k1a1,Aa2=k2a2,Aa3=k3a3 矩阵形式为AP=P∧或Ａ＝Ｐ∧Ｐ-1 说明矩阵A与对角矩阵相似。利用特征矩阵向量和特征值的方法可以求矩阵A的相似对角矩阵。矩阵的相似对角化方法可用计算一矩阵的