自动控制原理仿真MATLAB源程序，例子比较容易懂，含注释

% G = tf(10,conv([0.1 1],[0.5 1])) % Gc = feedback(G,1) % step(Gc) % ess = 1 - dcgain(Gc) %稳态误差 % Kp = dcgain(G) % ess1 = 1/(1+Kp) %稳态误差 % [num,den] = tfdata(Gc, 'v') % Kv = dcgain([num 0],den) % essv = 1/Kv % Ka = dcgain([num 0 0],den) % essa = 1/Ka % inf 无穷大 % figure(2) % impulse(Gc) %脉冲响应 figure(3) t = 0:0.1:100; T = [ 1 5 10]; for n = 1:3 G1 = tf([1], [T(n) 1]) y(:,n) = impulse(G1, t); % 返回的是列向量 y1 = impulse(G1, t); plot(t,y) title('脉冲响应曲线') xlabel('t/s') ylabel('amplitude') hold on end legend('T = 1','T = 5','T = 10') figure(4) plot(t,y) title('脉冲响应曲线') xlabel('t/s') ylabel('amplitude') legend('T = 1','T = 5','T = 10')