《线控》matlab程序作业.docxVIP

A=[0 1 0 0;0 0 1 0;0 0 0 1;-680 -176 -86 -6]; B=[0;0;0;1]; C=[100 20 10 0]; D=[0]; CME14=ss(A,B,C,D); CME14tf=tf(CME14) CME14tf = 10 s^2 + 20 s + 100 ------------------------------------------------ s^4 + 6 s^3 + 86 s^2 + 176 s + 680 %CME2_4 A=[0 1 0 0;0 0 1 0;0 0 0 1;-680 -176 -86 -6]; B=[0;0;0;1]; C=[100 20 10 0]; D=[0]; CME14=ss(A,B,C,D); CME14tf=tf(CME14) y_impulse=impulse(CME14) subplot(3,1,1); plot(y_impulse);grid; y_step=step(CME14); subplot(3,1,2); plot(y_step);grid; t=[0:.01:1]; u=[zeros(size(t))]; x0=[4;3;2;1]; y=lsim(CME14,u,t,x0); subplot(3,1,3); plot(y);grid; [Tdcf