control.m

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

% 这是一个多自由度受到脉冲荷载的情况
% author:chinamaker
% email:[e-mail]
% 2003年7月1号
% 动力学与控制技术论坛 http://www.dytrol.com
load('kobe.mat')                           % 导入地震波
M=[2500 0 0;0 2500 0;0 0 2500];                                      % 结构质量矩阵
K=[2500 -1250 0;-1250 2500 -1250;0 -1250 1250]*1000;     % 结构刚度矩阵
C=[10000 -5000 0;-5000 10000 -5000;0 -5000 5000];        % 结构阻尼矩阵
ling=zeros(3);                 % 定义零矩阵
danwei=eye(3);                 % 定义单位矩阵
% 转换到状态空间
A=[ling danwei;-K*inv(M) -C*inv(M)];
B=[ling;danwei];
C=[eye(6)];
D=[ling;ling];
% 检查系统的可控性
disp('系统的可控性'),M=ctrb(A,B)
disp('可控性矩阵的阶次'),rank(M)
sys=ss(A,B,C,D);
set(sys,'inputn',{'u1','u2','u3'})    % 定义输入矩阵
set(sys,'outputn',{'x1','x2','x3','x4','x5','x6'})  % 定义输出状态矩阵
kx=lqry(sys,eye(6),5e-2*danwei)                  
cl=feedback(sys,kx)