plotresult.m

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

% 这是一个多自由度受到脉冲荷载的情况
% author:chinamaker
% email:[e-mail]
% 2003年7月1号
% 动力学与控制技术论坛 http://www.dytrol.com
figure(1)
plot(tout,yout(:,1),'r',tout,yout(:,7),'b')
figure(2)
plot(tout,yout(:,2),'r',tout,yout(:,8),'b')
figure(3)
plot(tout,yout(:,3),'r',tout,yout(:,9),'b')
figure(4)
plot(tout,yout(:,4),'r',tout,yout(:,10),'b')
figure(5)
plot(tout,yout(:,5),'r',tout,yout(:,11),'b')
figure(6)
plot(tout,yout(:,6),'r',tout,yout(:,12),'b')
lpf=tf(30,[1 30])


% 将低通滤波器加入
sysdes=ss(lpf)*sys
kx=lqry(sysdes(1,1),1,1e-4)
estx=kalman(sysdes(2,:),eye(1),20)
regx=lqgreg(estx,kx)
figure(7)
bode(regx,{1 20})