第6.2次Matlab软件应用.ppt

7.3 MATLAB在特征值、 特征向量、 二次型中的应用 ;7.3.1 MATLAB在特征值和特征向量中的应用; 例7.11 分别求数值矩阵和符号矩阵A的特征值和特征向量。 其中 A＝ 解 （1）求数值矩阵A的特征值和特征向量。 解法① A=[-2,1,1;0,2,0;-4,1,3]；P=poly(A),r=roots(P) %求A的特征值。 运行后结果显示： ;解法② A=[-2,1,1;0,2,0;-4,1,3]； [X,D]=eig(A) 运行后结果显示： X＝ D＝ ; （2）求符号矩阵A的特征值和特征向量。 [B,C]=eigensys(A) 运行后结果显示： X＝ D＝ ;7.3.2 MATLAB在二次型中的应用 ;例7.12 求向量组a1=(1,1,0,0)’,a2=(1,0,1,0)’,a3=(-1,0,0,1)’的标准正交向量组。 解 以a1,a2,a3为列向量，输入矩阵： A=[1,1,-1;1,0,0;0,1,0;0,0,1]； B=orth(A) %求矩阵A的正交矩阵。 C=sym(B) %求矩阵B的符号矩阵。 运行后结果显示： ;运行后结果显示： B= C= 请思??矩阵B与矩阵C的区别。 ;例7.13 作一个正交变换X=PY，把下列二次型化成标准形。 f(x1,x2,x3,x4)=2x1x2+2x1x3-2x1x4-2x2x3+2x2x4+2x3x4 解 此二次型的实对称矩阵 A= ;在编辑器中建立文件f.m A=[0,1,1,-1;1,0,-1,1;1,-1,0,1;-1,1,1,0]； [U,T]=schur(A)； syms y1 y2 y3 y4 y=[y1;y2;y3;y4]； X=vpa(U,2)*y， f(y1,y2,y3,y4)=y’*T*y 运行后结果显示： ; U= T= X= f(y1,y2,y3,y4)=y1^2+y2^2-3*y3^2+y4^2