2.1 LQR 설계 문제에서 PID 제어기 동조
플랜트(제어 대상)가 미분 방정식으로 표현되는 2차 시스템이라 가정하면 식 (1)과 같다.
식 (1)에서 $y(t)$는 출력(output), $\omega_{n}$은 고유 주파수(natural frequency), $\zeta$는 감쇠비(damping
ratio), $u(t)$는 제어 입력(control input)이다. 그리고 출력 $y(t)$를 최종적으로 안정된 상태 $y(\infty)=0$이
되도록 제어하는 것이 목표이다. 식 (1)의 미분 방정식에서 초기 조건은 고려하지 않는다. LQR 제어시스템의 구조는 상태 궤환(state feedback)으로 오차를 제어할 수 없는 단점이
있다. 따라서 상태방정식의 상태 변수에서 출력 변수를 추가해야 명령 추종 성능을 만족하는 제어기 설계가 가능하다. 출력 변수와 출력의 미분값이 포함된
형태로 새로운 상태 변수를 정의하면 식 (2)와 같다.
식 (1)의 미분 방정식에서 새로운 상태 변수에 대한 식 (2)를 적용하면 식 (3)과 같은 상태 공간 방정식으로 표현할 수 있다.
PID 제어기와 LQR 제어시스템의 구조를 등가화하기 위해 순방향 루프에서 식 (5)와 같이 출력 변수에 적분 변수를 추가한 상태방정식을 유도하면 식 (6)과 같다.
그리고 식 (8)과 같은 가격 함수(cost function)를 고려해야 LQR 제어문제로 변환할 수 있다.
식 (8)에서 $Q$는 양의 반한정 대칭 행렬이고 $R=\rho I$는 양의 한정 대칭 행렬로 정의된다. LQR 제어기 설계 문제는 식 (8)에서 고려된 가격함수를 최소화하기 위해 설계요소인 $Q$와 $R$의 해를 구하는 것이며, 최적 제어 입력은 식 (9)와 같이 정의된다.[9]
여기서 제어 이득 $G$는 $R^{-1}B^{T}K$이며 $K=K^{T}$이다. 그리고 식 (9)에서의 $K$는 식 (10)에서 정의된 Riccati 방정식의 해로 구할 수 있다.
적분요소가 추가된 상태방정식에서 $x(t)\in R^{3}$ 이기 때문에 식 (9)에서 정의된 제어 이득 $G$는 식 (11)과 같이 분해할 수 있다.
식 (9)에서 유도된 죄적 제어 입력 $u(t)$는 식 (6)에서 유도된 상태방정식을 이용하여 식 (12)와 같이 분해 가능하면 PID 제어기의 형태로 정의할 수 있다. 일반적인 PID 제어의 형태는 식 (13)과 같다.
여기서, $K_{p},\: K_{I},\: K_{D}$ 는 각각, 비례 이득, 적분 시간, 미분 시간으로 PID 제어기 파라미터 요소들이다. 식 (13)에서 표현된 일반적인 PID 제어기의 설계요소와 식 (11)에서 정의된 제어 이득의 관계를 이용하면 LQR 제어기의 설계 파라미터들과 PID 제어기의 설계 파라미터들의 등가 관계를 식 (14)와 같이 정의할 수 있다.
따라서, PID 제어기 설계 파라미터인 $K_{p}$(비례 이득), $K_{I}$(적분 시간), $K_{D}$(미분 시간)을 동조하는 문제가 LQR
제어기 파라미터 동조문제로 변형할 수 있다. LQR 제어의 가중치 요소 $Q$와 $R$을 선정하고 식 (11)에서 유도된 제어 이득 $G$를 구하면 식 (14)에서 유도된 제어 이득과 PID 제어 설계 파라미터를 구할 수 있다.[7] 따라서 그림 1의 LQR 제어시스템의 구조와 그림 2의 PID 제어기 시스템의 구조는 등가가 된다.
그림 1. LQR 제어시스템의 블록 선도
Fig. 1. Block diagram of LQR control system
그림 2. PID 제어시스템의 블록선도
Fig. 2. Block diagram of PID control system
2.2 루프형성기법을 이용한 LQ-PID 제어기 동조
루프 전달함수의 특이값 동조기법[11]을 이용하여 LQR 의 설계 파라미터 Q와 R을 선택한다. 그림 2에서 플랜트 출력 측의 절단점 ①에서 구한 루프 전달함수 $H(s)$와 절단점 ②에서의 루프 전달함수 $G_{L}(s)$는 식 (15)와 식 (16)과 같이 정의된다.
단일 입출력 시스템에서는 식 (15)와 식 (16)과 같이 같다. LQR 제어 문제의 주파수역 등가(Frequency domain equality)를 사용하여 설계 파라미터 $Q=N^{T}N$과 $R=\rho
I$를 결정할 수 있고, LQR의 주파수역 특성은 식 (17)과 같다.
그리고 설계 파라미터 $N$은 식 (18)과 같이 분해할 수 있다.
명령 추종과 외란 제거 성능을 위한 저주파수 영역의 성능을 위한 조건은 식 (19)와 같이 정의할 수 있고 보드 선도를 나타내면 그림 3과 같다.
그림 3. 명령 추종 및 외란 제거에 대한 루프 형상
Fig. 3. Loop Shaping for command following and disturbance rejection
센서 잡음 및 모델링 오차에 대한 저감도를 위한 고주파수 영역의 성능을 위한 조건은 식 (20)과 같이 정의할 수 있고 보드 선도를 나타내면 그림 4와 같다.
2.3 영점 배치기법을 이용한 루프 형상
플랜트가 2차 시스템으로 가정할 때 두 극점이 존재한다. 제어기의 영점이 플랜트의 최대 극점과 멀리 떨어지면 고주파수 영역에서 낮은 이득을 갖게 된다.
그러나 제어기를 이용하여 설정되는 영점이 복소 평면의 원점으로부터 멀어지면 많은 에너지가 필요한 단점이 있다. 따라서 플랜트의 최대 극점과 제어기의
두 영점이 같게 되는 경우가 가장 효율적인 경우이다. 플랜트의 극점이 실근을 갖는 경우 식 (18)의 설계 파라미터 $n_{2}$의 범위는 식 (20)과 같다.
그림 4. 센서 잡음 및 모델링 오차에 대한 루프 형상
Fig. 4. Loop Shaping for Sensor Noise and Modeling error
플랜트의 극점이 공액 복소근을 갖는 경우 설계 파라미터 $n_{2}$의 범위는 식 (21)과 같다.