fig3_38.m

世界经典教材《现代控制系统》第八版MATLAB仿真源码

% Chapter 3: Figure 3.38, p. 162
%
% Model Parameters
k=10;
M1=0.02; M2=0.0005;
b1=410e-03; b2=4.1e-03;
t=[0:0.001:1.5];
% State Space Model
A=[0 0 1 0;0 0 0 1;-k/M1 k/M1 -b1/M1 0; k/M2 -k/M2 0 -b2/M2];
B=[0;0;1/M1;0]; C=[0 0 0 1]; D=[0];
% Simulated step Response
y=step(A,B,C,D,1,t);
plot(t,y), grid
xlabel('Time (sec)'), ylabel('y dot (m/s)')