工程技术中的一些问题,如振动和稳定性问题,常归结为求一.doc

工程技术中的一些问题，如振动和稳定性问题，常归结为求一个方阵的特征值和特征向量。 特征值与特征向量的求法 设A为n阶方阵，如果数“”和n维列向量x使得关系式成立，则称为方阵A的特征值，非零向量x称为A对应于特征值“”的特征向量。 求矩阵的特征值和特征向量 A=[-2 1 1;0 2 0;-4 1 3]; [V,D]=eig(A) V = -0.7071 -0.2425 0.3015 0 0 0.9045 -0.7071 -0.9701 0.3015 D = -1 0 0 0 2 0 0 0 2 即：特征值-1对应特征向量(-0.7071 0 -0.7071)T 特征值2对应特征向量(-0.2425 0 -0.9701)T和(-0.3015 0.9045 -0.3015)T 提高特征值的计算精度 [T,B] = balance(A) %求相似变换矩阵T和平衡矩阵B，满足。 B = balance(A) %求平衡矩阵B 复对角矩阵转化为实对角矩阵 [V,D] = cdf2rdf (v,d) %将复对角阵d变为实对角阵D，在对角线上，用2×2实数块代替共轭复数对。 例1-91 A=[1 2 3;0 4 5;0 -5 4]; [v,d]=eig(A) v = 1.0000 -0.0191 - 0.4002i -0.0191 + 0.4002i 0 0 - 0.6479i 0 + 0.6479i 0 0.6479 0.6479 d = 1.0000 0 0 0 4.0000 + 5.0000i 0 0 0 4.0000 - 5.0000i [V,D]=cdf2rdf(v,d) V = 1.0000 -0.0191 -0.4002 0 0 -0.6479 0 0.6479 0 D = 1.0000 0 0 0 4.0000 5.0000 0 -5.0000 4.0000 正交基 B=orth(A) %将矩阵A正交规范化，B的列与A的列具有相同的空间，B的列向量是正交向量，且满足：B*B = eye(rank(A))。 将矩阵正交规范化。 A=[4 0 0; 0 3 1; 0 1 3]; B=orth(A) Q=B*B 则显示结果为 P = 1.0000 0 0 0 0.7071 -0.7071 0 0.7071 0.7071 Q = 1.0000 0 0 0 1.0000 0 0 0 1.0000 二次型 求一个正交变换X=PY，把二次型 化成标准形。 解：先写出二次型的实对称矩阵 在Matlab编辑器中建立M文件如下： A=[0 1 1 -1;1 0 -1 1;1 -1 0 1;-1 1 1 0]; [P,D]=schur(A) syms y1 y2 y3 y4 y=[y1;y2;y3;y4]； X=vpa(P,2)*y %vpa表示可变精度计算，这里取2位精度 f=[y1 y2 y3 y4]*D*y 运行后结果显示如下： P = 780/989 780/3691 1/2 -390/1351 780/3691 780/989 -1/2 390/1351 780/1351 -780/1351 -1/2 390/1351 0 0 1/2 1170/1351 D = 1