ex-control

pkg load control

% 전송 함수 정의 (예: H(s) = 1 / (s^2 + 3s + 2))
s = tf('s');
H = 1 / (s^2 + 3*s + 2);

% 시스템의 전달 함수 출력
disp("전달 함수 H(s):")
disp(H)

% 시스템의 뿌리(특성 방정식의 해) 구하기
disp("시스템의 뿌리:")
roots([1, 3, 2])

% 시스템의 단계 응답 시뮬레이션 (시간 범위: 0에서 10초)
t = 0:0.01:10;
[y, t_out] = step(H, t);

% 응답 출력
disp("시스템의 단계 응답:")
plot(t_out, y)
title('단계 응답')
xlabel('시간 (초)')
ylabel('출력')
grid on

% 시스템의 안정성 분석: 뿌리 위치를 기반으로 안정성 확인
disp("시스템의 안정성 분석:")
pzmap(H)