[현대제어]현대제어 디지털제어 시뮬레이션 Matched pole zero method등

소개글

현대제어공학의 마지막단계. 학부생 기준으로 최고수준 단계인 디지털제어 시뮬레이션입니다.
Matched Pole Zero method, Tustin`s method 등을 사용해 디지털제어로 변환하여 시뮬레이션 했습니다.
A,B continuous space를 Pi, Gamma Discrete space로 변환해주는 함수도 만들었습니다.

목차

8.18-(b) Simulation
- 소스코드 -
-결과-

8.18-(d) Simulation
8.21-(a),(b) 시뮬레이션
-소스코드-

본문내용

- 소스코드 -
%p818.m
clear
A = [0 1;-.1 0];
B = [0;0.1];
C = [1 0];
D = 0;
[P,G]=c2dmatrix(A,B,1);
poles = [-0.5+sqrt(3/2)*i -0.5-sqrt(3/2)*i];
K = place(P,G,poles);
Pnew = P - G*K;
SYS = ss(Pnew,B,C,D);
SYSD=c2d(SYS,1,`zoh`);
step(SYSD)

% c2dmatrix.m
% 이 함수는 Continous State Equation을
% Discrete State Equation으로 바꾸는 함수이다.
% 입력 : F, G , T , 출력 : Pi, Gamma

function [Pi, Gamma] = c2dmatrix(F,G,T)
Pi = eye(size(F));
Temp2 = zeros(size(F));
i=1;
while norm(Pi-Temp2) > 0.001,
Temp2 = Pi;
Pi = Pi + F^i*T^i/factorial(i);
i = i+1;
end
Gamma = zeros(size(G));
Temp1 = ones(size(G));
i=0;
while norm(Gamma-Temp1) > 0.001,
Temp1 = Gamma;
Temp2 = zeros(size(F));
for j = 0 : 1 :i,
Temp2 = Temp2 + (F^j*T^(j+1))/factorial(j+1);
end
Gamma = Temp2 * G;
i = i + 1;
end

-결과-

참고 자료

없음