ME! 发表于 2013-2-28 16:57

我现在是用matlab画振型图,请教eig命令和null命令

=eig(kk,mm);      
V2=V(:,4);




对应公式det(kk-(lambda)^2*mm)=0,解这个特征方程,就可以得到特征值,再将每一个特征值代入det(kk-(lambda)^2*mm)=0求出对应的特征向量,然后画振型图
=sort(diag(D));                      % 给特征值和特征向量排序
A1=kk-lambda(4)*mm;% 组成齐次方程组,求i阶特征值lambda(i)的特征向量
V1=null(A1);   % 求齐次方程组的解,i阶特征值lambda(i)的特征向量

这两种办法结果不能对应
分别画的振型图是

页: [1]
查看完整版本: 我现在是用matlab画振型图,请教eig命令和null命令