我现在是用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]