fig3_26.m

来自「包含了控制理论的各种仿真程序」· M 代码 · 共 26 行

%  Figure 3.26     Feedback Control of Dynamic Systems, 4e
%                        Franklin, Powell, Emami
%
% fig3_26.m 
clf;
num=[2 1];    % form transfer function
den=[1 2 5];
sys=tf(num,den); % form system
t=0:.1:6;        % define time vector
y=impulse(sys,t); %compute impulse response
plot(t,y)
hold on
% compute 2e^(-t)
ne=2;
de=[1 1];
syse=tf(ne,de);
ye=impulse(syse,t);
plot(t,ye,':r')
plot(t,-ye,':r'), grid
axis([0 6 -2 2])
hold off
title('Figure 3.26')
xlabel('Time (sec)')
ylabel('h(t)')