特征向量都求出来了,用哪一种归一化还不跟玩儿一样吗?okok已经回答你了,这里再贴一下,主要因为已经消失几天了,出来透口气:复制内容到剪贴板代码:functionziyouzhendong1(k,m)
%k=600*[1-10;-13-2;0-25];m=diag([11.52]);
clc
[C,B]=eig(inv(m)*k);
n=size(m,2);
fori=1:n
Ct=C(:,i);
zhenxing(:,i)=Ct/Ct(1);
end
w=sqrt(diag(B));
B1=diag(B)';
Eigen=inv(m)*k;
fori=1:n
t1(:,i)=B1(i)*zhenxing(:,i);
t2(:,i)=Eigen*zhenxing(:,i);
end
B,zhenxing,t1,t2,w
其实完全可以不用循环,这是上学初学MATLAB一个月左右时写的程序,目的是帮助同班美女从繁重的结构动力学作业中解脱出来并有充分的时间买化妆品andkeepusamused...,当时还没人学MATLAB这玩意儿,全拿计算器算动力学矩阵的问题,显然即使最多只有4×4的K矩阵也足够让人受折磨了...
转载