第6章 线性代数实验.pptVIP

例1： A=[16 2 3 13; 5 11 10 8; 9 7 6 12; 4 14 15 1]; poly(A) ％直接求取 ans = 1 -34 -80 2720 -2.5612e-012 A=sym(A); poly(A) ％运用符号工具箱 ans =x^4-34*x^3-80*x^2+2720*x poly2sym应用举例 格式： f=poly2sym(P) 或 f=poly2sym(P,x) 例： 例2：用不同形式表示多项式 P=[1 2 3 4 5 6]; % 先由系数按降幂顺序排列表示多项式 f=poly2sym(P, x) % 以 x 为算子表示多项式 f =x^5 + 2*x^4 + 3*x^3 + 4*x^2 + 5*x + 6 P=sym2poly(f) P = 1 2 3 4 5 6 rref应用举例 例3：求向量组 jordan应用举例 例5：设矩阵 ,求一可逆矩阵 ,使 为对角矩阵. 6.2线性方程组求解 我们将线性方程组的求解问题分为两类：一类是方程组求唯一解或求特解，另一类是方程组求无穷解即通解。可以通过系数矩阵的秩来判断： ◆若系数矩阵的秩 r=n(n为方程组中未知变量的个数)，则有唯一解 ◆若系数矩阵的秩 rn，则可能有无穷解。 ◆线性方程组的无穷解 = 对应齐次方程组的通解+非齐次方程组的一个特解；其特解的求法属于解的第一类问题，通解部分属第二类问题。 (二） 第二类问题求解 在 Matlab 中，函数 null用来求解齐次线性方程组A·X=0 的解空间，实际上是求出解空间的一组基(基础解系)。 格式：z = null(A,’r’) % z 的列向量是方程 AX=0 的有理基 例8 求方程组 的通解 非齐次线性方程组需要先判断方程组是否有解，若有解，再去求通解。因此，步骤为： 第一步：判断 AX=b 是否有解，若有解则进行第二步 第二步：求 AX=b 的一个特解 第三步：求 AX=0 的通解 第四步：AX=b 的通解= AX=0 的通解+AX=b 的一个特解。 例9：求解方程组 解法二：在Matlab编辑器中建立 M文件：LX07212.m A=[1 1 -3 -1;3 -1 -3 4;1 5 -9 -8]; b=[1 4 0]; B=[A b]; C=rref(B) %求增广矩阵的行最简形，得最简同解方程组。 例10：求解方程组 （一） 平面上线性变换的几何意义 例11： 设x为二维平面上第一象限中的一个单位方块，其四个顶点的数据可写成 把不同的A矩阵作用于此组数据，可以得到多种多样的结果yi=Ai*x。用程序实现变换计算，并画出x及yi图形： x?[0,1,1,0;0,0,1,1]; subplot(2,3,1), fill([x(1,:),0],[x(2,:),0],r) A1?[?1,0;0,1], y1?A1*x subplot(2,3,2), fill([y1(1,:),0],[y1(2,:),0],g) … 几种变换的行列式与特征值 看出的基本关系 可以看出，矩阵A1使原图对纵轴生成镜像，矩阵A2使原图在横轴方向膨胀，矩阵A3使原图在纵轴方向压缩，矩阵A4使原图向右方剪切变形，矩阵A5使原图沿逆时针方向旋转t?pi/6。分别计算出这五个矩阵的行列式和特征值，如上页。 对二维空间（平面），一个变换所造成的图形的面积变化，取决于该变换的行列式。A1，A4和A5的行列式绝对值都是1，所以它们不会使变换后图形的面积发生改变。而A2和A3的行列式分别为1.5和0.2， 二维矩阵特征值的几何意义 二维矩阵的特征值表示该变换在原图形的特征向量的方向上的放大量。 例如矩阵A1在第一特征向量 方向的特征值为 ，即横轴正方向的增益为?1，其结果是把原图中横轴正方向的部分变换到新图的负方向去了； A1在第二特征向量 的方向的特征值为 ，即纵轴正方向的增益为1，因而保持了新图和原图在纵轴方向尺度不变。 （二） 对称矩阵与二次型主轴 对称矩阵的特点是所有元素关于主对角线对称，即A’?A。所以对称矩阵一定是方阵。前面曾要求读者特别注意A是对称矩阵时x与Ax的对应关系，其特点就是Ax呈椭圆形状，在椭圆的两个主轴方向，Ax与x在一条直线上长度差λ倍，即Ax?λx。当Ax与x方向相同时，λ为正数；当A