线代实验4(二次型).ppt

2、教材P2224、5题。 * * 数学实验之二次型 空间解析几何与线性代数实验之四 一、特征值与特征向量 格式1：T=eig (A) T为特征值列向量 　 格式2： [D,T] = eig (A) 其中 T为由特征值构成的对角阵 D为由特征向量构成的矩阵 [例5] 求方阵 的特征值与特征向量 创建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 特征值 的全部特征向量为： 其中 k1 , k2 不能同时为零. 特征值 的全部特征向量为： ，k3 不能为零 。 二、化二次型为标准型 格式 ： [P, T]=schur(A) 其中A为方阵，T为A的特征值所构成 的对角形矩阵， P为T对应的正交变换的 正交矩阵 。 [ 例6 ] 求一个正交变换，将二次型 化为标准型 二次型所对应的矩阵 运行结果： 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 创建schur01.m文件 A=[1 1 0 -1; 1 1 -1 0; 0 -1 1 1;-1 0 1 1]; [P, T]=schur(A) 所作的正交变换为： X=PY 二次型的标准型为： T = -1.0000 0 0 0 0 1.0000 0 0 0 0 1.0000 0 0 0 0 3.0000 1. 顺序主子式判断法 ⑴ 求二次型 F=X’AX的矩阵A的各阶顺序主子式 Di (i=1,2,3…..); ⑵ 判断 Di 是否大于0 . 三、 正定二次型的判定 建立函数文件shxu.m function [C,M]=shxu(A) % C为A的各阶顺序主子式组成的向量 % M为判定向量: if C(i)0, then M(i)=1; % others M(i)=0 n=size(A); C=[ ];M=[ ]; for i=1:n(1) A1=A([1:i],[1:i]); D=det(A1); C=[C D]; if D0 m=1; else m=0; end M=[M,m]; end 2、特征值判别法 ⑴ 求二次型 f =X’AX的矩阵A的全部特 征值 （i=1,2,……）； ⑵ 判断 是否大于0 . 建立函数文件tezh.m function [T,M]=tezh(A) n=size(A); T=(eig(A)); M=[ ]; for i=1:n(1) if T(i)0 m=1; else m=0; end M=[M,m]; end