知道刚度矩阵和质量矩阵后用matlab怎么求固有频率和固有振型
问题描述：
知道刚度矩阵和质量矩阵后用matlab怎么求固有频率和固有振型
答
转载
特征向量都求出来了,用哪一种归一化还不跟玩儿一样吗?okok已经回答你了,这里再贴一下,主要因为已经消失几天了,出来透口气：
复制内容到剪贴板代码: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
其实完全可以不用循环,这是上学初学MATLAB一个月左右时写的程序,目的是帮助同班美女从繁重的结构动力学作业中解脱出来并有充分的时间买化妆品and keep us amused...,当时还没人学MATLAB这玩意儿,全拿计算器算动力学矩阵的问题,显然即使最多只有4×4的K矩阵也足够让人受折磨了...% 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
转载