多自由度.m

文件多自由度.m为单独的文件

% 这是一个多自由度受到脉冲荷载的情况
% author:chinamaker
% email:[e-mail]
% 2003年7月1号
% 动力学与控制技术论坛 http://www.dytrol.com

% 数据文件
load('EI.dat')                           % 导入地震波
M=[2923.38 0 0;0 2923.38 0;0 0 2923.38];                                      % 结构质量矩阵
K=[2782.12 -1391.06 0;-1391.06 2782.12 -1391.06;0 -1391.06 1391.06]*1000;     % 结构刚度矩阵
C=[12747.48 -6373.74 0;-6373.74 12747.48 -6373.74;0 -6373.74 6373.74];        % 结构阻尼矩阵
ling=zeros(3);                 % 定义零矩阵
danwei=eye(3);
A=[ling danwei;-K*inv(M) -C*inv(M)];
B=[0 0 0 0 0 1]';
C=[danwei,ling];
D=[0 0 0]';
% 检查系统的可控性
disp('系统的可控性'),M=ctrb(A,B)
disp('可控性矩阵的阶次'),rank(M)
sys=ss(A,B,C,D);
t=0:0.02:20
y=impulse(sys,t);
plot(t,y(:,1),'r',t,y(:,2),'g',t,y(:,3),'b')