基于传递函数的汽车瞬态响应特性仿真

k11=[-23147,-33124]; k22=[-38318,-69972]; u=20; m=3018; a=1.84; b=1.88; z=10437; for i=[1,2] k1=k11(i); k2=k22(i); a11=(k1+k2)/(m*u); a12=(a*k1-b*k2)/(m*u^2)-1; a21=(a*k1-b*k2)/z; a22=(a^2*k1+b^2*k2)/(u*z); b11=-k1/(m*u); % b12=-k2/(m*u); b21=-a*k1/z; % b22=b*k2/z; aa=[a11 a12;a21 a22]; bb=[b11;b21]; cc=[1 0; 0 1]; dd=[0 ; 0]; delta=pi/18; g=ss(aa,bb,cc,dd);%状态空间方程 opt = stepDataOptions('InputOffset',0,'StepAmplitude',delta); [y,t] = step(g,10,opt); plot(t,y(:,2),'Linewidth',1); hold on; end legend('k1=-23147N/rad,k2=-38318N/rad','k1=-33124N/rad,k2=-69972N/rad'); xlabel('仿真时间（s）'); ylabel('横摆角速度（rad/s）');