特征向量都求出來(lái)了,用哪一種歸一化還不跟玩兒一樣嗎?okok已經(jīng)回答你了,這里再貼一下,主要因?yàn)橐呀?jīng)消失幾天了,出來(lái)透口氣：復(fù)制內(nèi)容到剪貼板代碼:function ziyouzhendong1(k,m)
% k=600*[1 -1 0;-1 3 -2;0 -2 5];m=diag([1 1.5 2]);
clc
[C,B]=eig(inv(m)*k);
n=size(m,2);
for i=1:n
Ct=C(:,i);
zhenxing(:,i)=Ct/Ct(1);
end
w=sqrt(diag(B));
B1=diag(B)';
Eigen=inv(m)*k;
for i=1:n
t1(:,i)=B1(i)*zhenxing(:,i);
t2(:,i)=Eigen*zhenxing(:,i);
end
B,zhenxing,t1,t2,w其實(shí)完全可以不用循環(huán),這是上學(xué)初學(xué)MATLAB一個(gè)月左右時(shí)寫(xiě)的程序,目的是幫助同班美女從繁重的結(jié)構(gòu)動(dòng)力學(xué)作業(yè)中解脫出來(lái)并有充分的時(shí)間買(mǎi)化妝品and keep us amused...,當(dāng)時(shí)還沒(méi)人學(xué)MATLAB這玩意兒,全拿計(jì)算器算動(dòng)力學(xué)矩陣的問(wèn)題,顯然即使最多只有4×4的K矩陣也足夠讓人受折磨了...



轉(zhuǎn)載