루프 정형화

편집 역사 토론
최적설계 수치해석 전자공학 마지막 수정: 2026-08-14 04:28:07

1. 개요[편집]

루프 정형화
Loop Shaping
대상개루프 $L = GK$의 특이값 모양
고전적 처방저주파 고이득 · 교차 근방 기울기 $-20$ dB/dec · 고주파 롤오프
현대적 2단 절차가중함수 정형화 → 정규 호환인수 $H_\infty$ 강건 안정화
제안McFarlane & Glover (1990, 1992)
여유 지표$\varepsilon_{\max} = 1/\gamma_{\min}$, 실무 기준 $\gamma_{\min} < 4$
도구MATLAB ncfsyn · loopsyn

제어기를 설계하지 말고, 개루프의 그림을 먼저 그려라. 제어기는 그 그림을 완성하는 도구다.

루프 정형화(loop shaping)는 폐루프 성능과 강건성을 개루프 전달함수 L=GKL = GK의 주파수 응답 모양으로 번역해 놓고, 그 모양을 직접 빚는 방식으로 제어기를 설계하는 절차다. 폐루프 사양은 서로 얽혀 있어 직접 다루기 어렵지만, 개루프는 보상기를 곱하면 보드 선도 위에서 그대로 더해지므로 손으로 그릴 수 있다는 것이 발상의 전부다. 보드가 1940년대에 정립한 이 감각은 SISO 설계의 고전이었고, 1980년대에 다변수·강건성 보장이 붙으면서 현대적 절차로 부활했다.

오늘날 루프 정형화라고 하면 대개 맥팔레인-글로버의 2단 절차를 가리킨다. 먼저 가중함수 W1,W2W_1, W_2로 원하는 개루프 모양을 자유롭게 빚고, 그 다음 그 모양을 정규 호환인수(normalized coprime factor) 불확실성에 대해 HH_\infty로 강건 안정화한다. 설계자는 익숙한 주파수 감각을 그대로 쓰고, 강건성 보장은 기계가 붙여 주는 분업이다.1

2. 왜 개루프의 모양인가[편집]

폐루프 사양은 결국 감도 함수 S=(I+L)1S = (I+L)^{-1}과 상보감도 T=L(I+L)1T = L(I+L)^{-1}에 대한 요구인데, 두 함수 모두 LL의 비선형 함수라 직접 빚기가 불편하다. 그런데 이득이 극단으로 가면 관계가 단순해진다.

σˉ(S)1σ(L)    (σ(L)1),σˉ(T)σˉ(L)    (σˉ(L)1)\bar\sigma(S) \approx \frac{1}{\underline\sigma(L)} \;\; (\underline\sigma(L) \gg 1), \qquad \bar\sigma(T) \approx \bar\sigma(L) \;\; (\bar\sigma(L) \ll 1)

저주파에서는 LL의 최소 특이값이 성능을 지배하고, 고주파에서는 최대 특이값이 강건성을 지배한다. 그러면 사양은 개루프에 대한 두 개의 울타리로 번역된다 — 저주파에서 σ(L)\underline\sigma(L)을 어떤 곡선 위로, 고주파에서 σˉ(L)\bar\sigma(L)을 어떤 곡선 아래로. 다변수에서 특이값이 두 개(최대·최소) 등장하는 것이 SISO 보드 선도와의 유일한 차이이며, 그 사이의 벌어짐이 곧 채널 간 이득 불균형이다.

문제는 교차 주파수 근방이다. 여기서는 위의 근사가 둘 다 깨지고, 1+L|1 + L|이 작아지면서 SSTT가 동시에 부풀 수 있다. 보드의 이득-위상 관계(최소위상계에서 이득 기울기가 위상을 결정한다)에 따르면 기울기 20n-20n dB/dec는 대략 위상 90n-90n^\circ에 대응하므로, 교차 근방에서 기울기를 20-20 dB/dec 근처로 완만하게 유지하라는 국룰이 나온다. 기울기가 40-40 dB/dec면 위상이 180-180^\circ에 붙어 위상 여유가 사라진다. 저주파는 적분기로 세우고, 고주파는 잡음·미모형 동특성 때문에 40-40 dB/dec 이상으로 떨어뜨리되, 교차 근방만은 완만하게 — 이 세 문장이 고전적 루프 정형화의 전부다.

3. 그릴 수 없는 모양[편집]

무엇이든 그릴 수 있는 것은 아니다. 플랜트가 물리적으로 부과하는 벽이 있고, 이를 무시한 목표 모양은 어떤 제어기로도 실현되지 않는다.

  • 우반평면 영점은 대역폭에 천장을 씌운다. 실수 우반평면 영점 zz에 대해 ωBz/2\omega_B \lesssim z/2. 지연 θ\theta가 있으면 파데 근사가 만드는 우반평면 영점 때문에 사실상 ωB1/θ\omega_B \lesssim 1/\theta다.
  • 불안정 극점은 바닥을 깐다. 우반평면 극점 pp에 대해 ωB2p\omega_B \gtrsim 2p. 둘 다 있으면 대략 z>4pz > 4p가 존재 조건이 된다.
  • 워터베드. 보드 감도 적분 정리가 lnS\ln|S|의 총면적을 못 박으므로, 저주파에서 눌러 얻은 것은 교차 바로 위에서 정확히 갚아야 한다.
  • 액추에이터 대역폭과 포화. KSKS가 곧 제어 입력이므로, 개루프를 아무리 예쁘게 그려도 KS\lVert KS \rVert_\infty가 크면 실물에서 밸브가 시달린다.

이 제약들의 유도와 정확한 부등식은 감도 함수비최소위상계에서 다룬다. 루프 정형화 관점에서 중요한 것은 이것들이 “설계를 잘하면 넘을 수 있는 벽”이 아니라는 사실이다. 아래에서 볼 γ\gamma 값이 이 벽에 부딪혔음을 알려 주는 신호등 역할을 한다.

4. 정규 호환인수 불확실성[편집]

2단 절차의 뒷부분을 이해하려면 불확실성 모형이 먼저다. 플랜트를 좌측 정규 호환인수분해

G=M~1N~,N~N~+M~M~=IG = \tilde M^{-1} \tilde N, \qquad \tilde N \tilde N^* + \tilde M \tilde M^* = I

로 쓰고, 섭동을 두 인수에 동시에 준다.

GΔ=(M~+ΔM)1(N~+ΔN),[ΔNΔM]<εG_\Delta = (\tilde M + \Delta_M)^{-1}(\tilde N + \Delta_N), \qquad \left\lVert \begin{bmatrix} \Delta_N & \Delta_M \end{bmatrix} \right\rVert_\infty < \varepsilon

이 모형이 곱셈형·가법형보다 나은 점이 결정적이다. 곱셈형 불확실성 (1+WΔ)G(1 + W\Delta)G극점의 개수를 바꾸지 못한다 — 명목이 안정이면 섭동된 플랜트도 안정이어야 한다는 암묵적 가정이 깔려 있다. 반면 호환인수 섭동은 M~\tilde M이 흔들리면서 극점이 허수축을 넘어갈 수 있어, 안정한 명목 모델에서 출발해 불안정한 실제 플랜트까지 공을 덮는다. 물리적으로 이득·위상·극점 배치가 함께 변하는 상황을 자연스럽게 표현하며, 이 공은 갭 계량(gap metric)의 공과 사실상 같다.

강건 안정화 문제 — 이 공 전체를 안정화하는 KKε\varepsilon을 최대화하는 것 — 은 작은 이득 정리에 의해 다음 HH_\infty 노름의 최소화와 동치가 된다.

γ=[KI](IGK)1[GI],εmax=1γmin\gamma = \left\lVert \begin{bmatrix} K \\ I \end{bmatrix} (I - GK)^{-1} \begin{bmatrix} G & I \end{bmatrix} \right\rVert_\infty, \qquad \varepsilon_{\max} = \frac{1}{\gamma_{\min}}

괄호 안의 블록이 정확히 감도 함수에서 말하는 갱 오브 포 S,T,KS,GSS, T, KS, GS의 묶음이라는 점을 알아보는 것이 이 절차의 핵심 통찰이다. 즉 이 하나의 숫자를 줄이면 네 전달함수가 동시에 통제된다.

5. 명시 해 — γ 반복이 없다[편집]

이 문제가 특별한 이유는 닫힌 형태의 답이 있다는 것이다. 일반적인 HH_\infty 합성은 γ\gamma를 주고 리카티 조건을 검사한 뒤 이분법으로 좁혀 들어가야 하지만, 여기서는 γmin\gamma_{\min}이 바로 계산된다.

γmin=1+λmax(XZ)\gamma_{\min} = \sqrt{1 + \lambda_{\max}(XZ)}

X,ZX, Z는 각각 제어·필터 쪽 대수 리카티 방정식

AX+XAXBBX+CC=0,AZ+ZAZCCZ+BB=0A^\top X + XA - XBB^\top X + C^\top C = 0, \qquad AZ + ZA^\top - ZC^\top CZ + BB^\top = 0

의 안정화 해다. 동등하게 εmax=1[N~    M~]H2\varepsilon_{\max} = \sqrt{1 - \lVert [\tilde N \;\; \tilde M] \rVert_H^2}로도 쓰이며(H\lVert \cdot \rVert_H는 한켈 노름), 제어기 자체도 X,ZX, Z의 명시적 조합으로 떨어진다. 가중함수를 바꿀 때마다 반복 없이 즉시 γ\gamma를 받아 볼 수 있다는 것이 실무에서 이 절차가 사랑받는 이유다. 설계 루프가 “모양 수정 → γ\gamma 확인”의 몇 초짜리 사이클이 된다.

γ\gamma를 읽는 법이 이 절차의 진짜 노하우다.

  • γmin2\gamma_{\min} \lesssim 2 — 그린 모양이 플랜트와 잘 맞는다. 그대로 간다.
  • 2γmin42 \lesssim \gamma_{\min} \lesssim 4 — 쓸 만하다. 여유가 빠듯하면 모양을 조금 눕힌다.
  • γmin4\gamma_{\min} \gtrsim 4모양 자체가 플랜트와 양립하지 않는다. 제어기를 탓할 일이 아니라 목표 대역폭이 우반평면 영점을 넘겼거나, 롤오프를 불안정 극점보다 낮은 곳에 걸었거나, 가중함수가 플랜트의 우반평면 극·영점을 상쇄하려 들었다는 뜻이다.

γ\gamma는 성능 점수가 아니라 사양 진단기다. 이 해석이 있어서 앞 절의 물리적 벽들이 설계 과정에서 숫자로 되돌아온다.2

6. 2단 절차[편집]

  1. 정형화. 전보상기 W1W_1과 후보상기 W2W_2를 골라 정형화된 플랜트 Gs=W2GW1G_s = W_2 G W_1의 특이값이 원하는 모양을 갖게 한다. W1W_1에는 대개 적분 작용(1/s1/s 또는 (s+a)/s(s+a)/s)과 저역 이득을, W2W_2에는 채널 스케일링과 센서 대역폭 반영을 넣는다. 다변수에서는 특이값 정렬(alignment)도 W2W_2의 일이다 — 교차 근방에서 최대·최소 특이값이 크게 벌어져 있으면 채널마다 여유가 달라져 설계가 어려워진다.
  2. 강건 안정화. GsG_s의 정규 호환인수에 대해 ε\varepsilon을 최대화하는 KK_\infty를 위 명시 공식으로 구한다. γmin\gamma_{\min}이 크면 1번으로 돌아간다.
  3. 조립. 최종 제어기는 K=W1KW2K = W_1 K_\infty W_2다. 실제 개루프는 W2GW1KW_2 G W_1 K_\infty가 아니라 GKGK이므로, 최종 확인은 반드시 원래 플랜트에 대해 한다.

2단계가 1단계의 모양을 얼마나 훼손하는지도 정량적으로 보장된다. 글로버-맥팔레인의 결과에 따르면 저주파에서의 최소 특이값 손실과 고주파에서의 최대 특이값 초과가 γ\gamma로 상계되므로, γ\gamma가 작게 나왔다면 그린 대로 나온다고 믿어도 된다. 이것이 “고전적 직관 + 현대적 보장”이라는 홍보 문구의 실체다.

혼합 감도 HH_\infty 설계와 비교하면 차이가 분명하다. 혼합 감도는 가중함수를 조금 바꿀 때마다 γ\gamma 반복을 다시 돌려야 하고, 성능 가중에 넣은 적분기가 플랜트의 극·영점과 상쇄되면서 불안정 상쇄를 만드는 사고가 잦다. 정규 호환인수 접근은 상쇄가 구조적으로 일어나지 않고, 가중함수의 안정성 요구도 느슨하다.

7. 남는 문제들[편집]

  • 기준값 추종은 별도로 다뤄야 한다. 강건 안정화는 외란 억제 위주의 1자유도 구조라, 오버슈트가 큰 폐루프가 나올 수 있다. 표준 처방은 목표 모델 TrefT_{\text{ref}}를 추가로 사양화하는 2자유도 확장(호일-하이드-라임비어, 1991)이며, 하이드-글로버가 VSTOL 항공기 비행제어에 적용한 사례가 이 확장의 대표적 시연이다.
  • 제어기 차수. KK_\infty의 차수는 GsG_s의 차수와 같으므로 플랜트 차수 + 가중함수 차수다. 6차 플랜트에 3차 가중을 두 개 붙이면 12차 제어기가 나온다. 실장 전 축소차수모델의 평형 절단으로 깎되, 정형화 가중은 남기고 KK_\infty만 줄이는 것이 관행이다.
  • 구조적 불확실성에는 약하다. 호환인수 공은 여전히 비구조적이라, 여러 물리 파라미터가 독립적으로 변하는 경우 보수적이다. 이때는 구조적 특이값 μ\mu와 D-K 반복으로 넘어간다. 실무 워크플로는 대개 “루프 정형화로 초안 → μ\mu로 검증 → 필요하면 μ\mu 합성”이다.
  • 모양을 그리려면 플랜트를 알아야 한다. 저주파 이득과 교차 근처 위상을 모르면 W1W_1을 고를 수 없다. 앞단에 시스템 식별 실험이 붙는 것이 정상이며, 지연이 지배적인 공정이라면 계전기 되먹임으로 임계점 하나를 찍는 것만으로도 출발점이 나온다.
  • 비선형·시변계에는 그대로 못 쓴다. 운전점마다 모양을 그리고 보간하는 게인 스케줄링이 현장 관행이지만, 보간 구간의 안정성은 별도 보장이 필요하다.

8. 관련 문서[편집]

9. Footnotes[편집]

  1. 이 분업이 실제로 잘 작동하는지는 사람에 따라 평가가 갈린다. “가중함수 두 개만 고르면 된다”는 홍보와 달리, 좋은 W1W_1을 고르는 감각은 결국 고전적 보드 설계를 손으로 해 본 경험에서 온다. 도구가 바뀌어도 배워야 할 것은 그대로라는, 이 바닥에서 반복되는 결말이다.

  2. 초보자가 γ=7\gamma = 7을 받아 들고 하는 일은 대개 가중함수의 차수를 올리는 것이다. 거의 언제나 헛수고다. γ\gamma가 큰 이유는 모양이 정교하지 않아서가 아니라 모양이 물리적으로 불가능해서이므로, 목표 교차 주파수를 반으로 내려 보는 것이 30초 만에 원인을 확인하는 방법이다.