전달함수 표현법에 의한 시뮬레이션

%G(s)=1/(s^2+3s+2), U=sin(t)
gs=tf([1], [1 3 2]);
%step(gs);
t=0:0.01:4*pi;
U=sin(t);
[Y, T] = lsim(gs, U, t);
plot(T, [Y, U']); grid on;

상태변수표현법에 의한 시뮬레이션

%A=[0 1; -2 -3]; b=[0;1]; c=[1 0];
A=[0 1; -2 -3];
b=[0;1];
C=[1 0];
D=0;
gs=ss(A, b, C, D);
%step(gs);
t=0:0.01:4*pi;
U=ones(1,length(t));%sin(t);
[Yy, Tt] = lsim(gs, U, t);
plot(T, [Yy, U']); grid on;

figure(2);
plot(T, Y-Yy); grid on;