如何出一道计算机仿真题,计算机仿真试题

(2)仿真框图：

(3)程序：

m1=40;

m2=400;

k1=158;

k2=23;

k1=1000*k1;k2=1000*k2; %给参数赋值

p(1)=600;

p(2)=1000;

p(3)=2000;

p(4)=4000; %阻尼系数的值

for n=1:4

b=p(n);fangz; %调用状态方程

A=[0 0 1 0;0 0 0 1;-(k1+k2)/m1 k2/m1 -b/m1 b/m1;k2/m2 -k2/m2 b/m2 -b/m2]; B=[0; 0; k1/m1; 0];

C=zeros(1,4);C(2)=1;

D=0;

sys(n)=ss(A,B,C,D); %建立系统状态空间模型

end

figure('name','悬架系统的阶跃响应','numbertitle','off');

t=0:0.01:25;

step(sys(1),sys(2),sys(3),sys(4),t); %系统单位阶跃响应曲线

grid;title('簧上质量位移响应曲线'); %添加曲线标题

xlabel('响应时间');ylabel('响应值'); %添加X、Y坐标轴名称

figure('name','悬架系统的bode图','numbertitle','off');

w=logspace(0.3,3.0);

bode(sys(1),sys(2),sys(3),sys(4),w); %绘制系统伯德图

grid;title('簧上质量的bode图');

disp('display:matrix A ,B');

disp('matrix A='),disp(A), %输出状态空间变量系数矩阵A

disp('matrix B='),disp(B), %输出状态空间输入系数矩阵B