线性方程组与二次型.pptVIP

四、特征值与特征向量 格式1：T=eig (A) 其中T为特征值列向量 　 格式2： [D,T] = eig (A) 其中T为由特征值构成的对角阵 D为由特征向量构成的矩阵 [例6] 求方阵 的特征值与特征向量 创建eig01.m文件 A=[2 2 -2; 2 5 -4; -2 -4 5]; C=eig(A) [D,T]=eig(A) 运行结果： C = 1.0000 1.0000 10.0000 D = -0.2981 0.8944 0.3333 -0.5963 -0.4472 0.6667 -0.7454 0 -0.6667 T = 1.0000 0 0 0 1.0000 0 0 0 10.0000 特征值 = = 1的全部特征向量为： 其中 k1 , k2 不能同时为零. 特征值 = 10 的全部特征向量为： ，k3 不能为零 。 五、 二次型化标准型 格式 ： [P, T]=schur(A) 其中A为方阵，T为A的特征值所构成 的对角形矩阵， P为T对应的正交变换的 正交矩阵 。 [ 例7 ] 求一个正交变换，将二次型 化为标准型 二次型所对应的矩阵 创建schur01.m文件 A=[1 1 0 -1; 1 1 -1 0; 0 -1 1 1;-1 0 1 1]; [P, T]=schur(A) 运行结果： P = -0.5000 0.7071 0.0000 0.5000 0.5000 -0.0000 0.7071 0.5000 0.5000 0.7071 0.0000 -0.5000 -0.5000 0.0000 0.7071 -0.5000 T = -1.0000 0 0 0 0 1.0000 0 0 0 0 1.0000 0 0 0 0 3.0000 所作的正交变换为： X=PY 二次型的标准型为： 六、 正定二次型的判定 1. 顺序主子式判断法 ⑴ 求二次型 的矩阵A的各阶顺序主子式 Di (i=1,2,3…..); ⑵ 判断 Di 是否大于0 . 建立函数文件shxu.m function C=shxu(A) % C为A的各阶顺序主子式组成的向量 n=size(A); C=[ ]; for i=1:n(1) A1=A([1:i],[1:i]); D=det(A1); C=[C D]; end 2、特征值判别法 ⑴ 求二次型 的矩阵A的全部特征值 （i=1,2,……）； ⑵ 判断 是否大于0 . [例8] 判定下列二次型是否正定 二次型矩阵 方法一： 顺序主子式法 创建shxu01.m文件 A=[1 -1 2 1 -1 3 0 -3 2 0 9 -6 1 -3 -6 19 ]; C=shxu(A) 运行结果： C = 1 2 6 24 方法二 特征值法 创建tezh01.m文件 A=[1 -1 2 1 -1 3 0 -3 2 0 9 -6 1 -3 -6 19]; T=eig(A) 运行结果： T = 0.0643 2.2421 7.4945 22.1991 上机实验题 1、已知向量组A