극점과 영점 편집하기 (부분)

== MATLAB 예제 ==
인공위성이 자세 제어를 위해서 엔진에서 추력을 분출해야 한다. 이때 <code>input1</code>은 한 방향으로 주는 것이고 <code>input2</code>는 한 방향으로 주었다가 반대 방향으로도 동일한 크기로 주었을 때 어떠한 시스템의 응답을 보이는지에 대해 표현한 예제이다.

<syntaxhighlight lang='matlab'>
%추가 코드 (극점, 영점 계산하기)
num = [0 0 0.001];
den = [1 0.05 0];
[z,p,k] = tf2zp(num,den); 

% input 1에 대해서 sysG (전달함수)의 output
s = tf('s');
sysG = 0.0002/s^2;
t = 0:0.01:10;
u1 = [zeros(1,500) 25*ones(1,10) zeros(1,491)];
[y1] = lsim(sysG, u1, t);
y1 = y1 * 180/pi;

% input 2에 대해서 sysG (전달함수)의 output
u2 = [zeros(1,500) 25*ones(1,10) zeros(1,100) -25*ones(1,10) zeros(1,381)];
[y2] = lsim(sysG, u2, t);
y2 = y2 * 180/pi;

figure()
subplot(2,2,1)
plot(t,u1)
title('input 1')
xlabel('Time (sec)')
ylabel('Thrust, Fc (input 1)')
grid on

subplot(2,2,2)
plot(t,y1)
title('output (input 1)')
xlabel('Time (sec)')
ylabel('Theta (radian)')
grid on

subplot(2,2,3)
plot(t,u2)
title('input 2')
xlabel('Time (sec)')
ylabel('Thrust, Fc (input 2)')
grid on

subplot(2,2,4)
plot(t,y2)
title('output 2 (input 2)')
xlabel('Time (sec)')
ylabel('theta (radian)')
grid on
</syntaxhighlight>

[https://domybestinlife.tistory.com/313 CC BY로 가져옴]