shjsbshs

% ① 比例环节 G1(s)=1
sys1 = tf(1, 1);
step(sys1);
title('比例环节 G1(s)=1 的单位阶跃响应');
figure;
% ② 惯性环节 G1(s)=1/(s + 1)
sys2 = tf(1, [lbk]1 1[rbk]);
step(sys2);
title('惯性环节 G1(s)=1/(s + 1) 的单位阶跃响应');
figure;
% ③ 积分环节 G1(s)=1/s
sys3 = tf(1, [lbk]1 0[rbk]);
step(sys3);
title('积分环节 G1(s)=1/s 的单位阶跃响应');
figure;
% ④ 微分环节 G1(s)=s
sys4 = tf([lbk]1 0[rbk], 1);
step(sys4);
title('微分环节 G1(s)=s 的单位阶跃响应');
figure;
% ⑤ 比例+微分环节（PD）G1(s)=s + 2
sys5 = tf([lbk]1 2[rbk], 1);
step(sys5);
title('比例+微分环节（PD）G1(s)=s + 2 的单位阶跃响应');
figure;
% ⑥ 比例+积分环节（PI）G1(s)=1 + 1/s
sys6 = tf([lbk]1 1[rbk], [lbk]1 0[rbk]);
step(sys6);
title('比例+积分环节（PI）G1(s)=1 + 1/s 的单位阶跃响应');
figure;
T_values = [lbk]0.1, 1, 10[rbk]; % 选取不同的T值
for i = 1:length(T_values)
T = T_values(i);
sys = tf(1, [lbk]T 1[rbk]);
step(sys);
hold on;
end
legend('T = 0.1', 'T = 1', 'T = 10');
title('G(s)=1/(T*s + 1) 随T变化的单位阶跃响应');