카테고리
산업계측제어기술사

Nyquist Plot

Open Loop Transfer Function Pole이 원점을 지나는 경우

Gain Margin 구하는 방법 예제

“`html

Nyquist Plot

Nyquist Contour · Frequency Response · Origin Pole · Gain Margin · Phase Margin

1. 예제: G(s) = 1/(s+1)²

G(s) = 1/(s+1)2

Frequency response를 구하기 위해

s = jω

를 대입한다.

G(jω) = 1/(jω+1)2

분모를 전개하면

G(jω) = 1/(1-ω2+j2ω)

켤레복소수를 곱하여 실수부와 허수부로 정리하면

G(jω) = [(1-ω2) – j2ω] / [(1-ω2)2 + 4ω2]

2. Magnitude와 Phase

|G(jω)| = 1 / √{ (1-ω2)2 + 4ω2 }
∠G(jω) = tan-1 [ -2ω/(1-ω2) ]

ω=0에서

|G(j0)| = 1
∠G(j0) = 0°

따라서 Nyquist plot은 실수축의

1 + j0

에서 시작한다.

3. ω → ∞

ω가 매우 커지면

|G(jω)| ≈ 1/ω2 → 0

이므로 Nyquist plot은 원점으로 접근한다.

두 pole의 영향으로 총 위상은 최종적으로 -180° 방향으로 접근한다.

G(s)=1/(s+1)²의 Nyquist plot은 1+j0에서 시작하여 하반평면을 지나 원점으로 접근하고, 음의 주파수 영역은 실수축에 대한 대칭으로 나타난다.

4. ω = 1에서의 값

ω=1을 대입하면

G(j) = 1/(j+1)2
= 1/(j2)
= -j/2

따라서

|G(j)| = 1/2
∠G(j) = -90°
Nyquist plot은 ω=1에서 -j0.5 지점을 통과한다.

5. 원점에 Pole이 있는 경우

원문 예제에서는 open-loop transfer function의 pole이 원점에 존재하는 경우를 다룬다.

G(s) = (s+2)/s2

이 경우 Nyquist contour가 원점의 pole을 직접 지나갈 수 없다.

따라서 s-plane에서 원점을 작은 반원으로 우회한다.

s = εejθ

여기서

ε → 0

이다.

Nyquist contour 위에 pole이 존재하면 그 pole을 작은 반원으로 우회해야 한다.

6. 원점 Pole 우회 경로의 Mapping

원점 부근에서는

G(s) = (s+2)/s2

에

s = εejθ

를 대입한다.

G(s) = [εejθ+2] / [ε²ej2θ]

ε→0이면 분자의 2가 지배적이므로

G(s) ≈ 2 / [ε²ej2θ]
= (2/ε²)e-j2θ

따라서

|G(s)| → ∞
∠G(s) = -2θ
s-plane에서 원점을 우회하는 매우 작은 반원은 GH-plane에서는 무한히 큰 반원으로 mapping된다.

7. Nyquist Stability Criterion

Nyquist 판정에서는 다음 세 값을 사용한다.

기호 의미
P Open-loop transfer function의 RHP pole 개수
N Nyquist plot이 -1점을 둘러싸는 횟수
Z Closed-loop characteristic equation의 RHP zero 개수
Z = P + N

여기서 N의 부호는 사용하는 encirclement 방향 정의에 따라 달라질 수 있다.

Closed-loop 안정조건은 Z=0이 되도록 -1점의 encirclement 수가 결정되는 것이다.

8. Gain Margin 예제

원문 예제의 open-loop transfer function은

G(s) = K / [s(s+3)(s+5)]

이다.

먼저 K=1로 두고 Nyquist plot을 그린다.

G(jω) = 1 / [jω(jω+3)(jω+5)]

정리하면 원문에서는

G(jω) = [-8ω2 – j(15ω-ω3)] / [64ω4 + ω2(15-ω2)2]

로 나타낸다.

9. Negative Real Axis 교점

Nyquist plot이 음의 실수축과 만나는 지점에서는 imaginary part가 0이다.

15ω – ω3 = 0
ω(15-ω2) = 0

ω≠0인 교점은

ω = √15

이다.

이때 real part는 약

Re[G(jω)] ≈ -0.0083

이다.

K=1에서 Nyquist plot은 ω=√15 rad/s에서 약 -0.0083+j0 지점을 통과한다.

10. Gain Margin 계산

안정 한계에서는 해당 교점이 critical point -1+j0에 도달한다.

K × 0.0083 = 1

따라서

K ≈ 1/0.0083
K ≈ 120.5

즉 K=1에서 시작하면 gain을 약 120.5배까지 증가시킬 수 있다.

GM = 120.5

dB로 나타내면

GMdB = 20log10(120.5)
≈ 41.6 dB
K=120.5에서 Nyquist plot이 -1+j0에 도달하여 시스템은 marginally stable이 된다.

이때 진동 주파수는

ω = √15 rad/s

이다.

11. Gain Margin의 의미

Gain Margin은 Nyquist plot이 -180° 위상에서 critical point -1까지 얼마나 더 gain을 증가시킬 수 있는지를 나타낸다.

GM = 1 / |G(jωpc)H(jωpc)|

여기서

∠G(jωpc)H(jωpc) = -180°

이다.

12. Gain Margin 필기 예제

원문의 다른 예제에서는

G(s)H(s) = K / [(s²+2s+2)(s+2)]

로 두고 K=6인 경우 Gain Margin과 Phase Margin을 계산한다.

s=jω를 대입하면 분모는

(-ω²+2+j2ω)(2+jω)

이고 정리하면

-4(ω²-1) + jω(6-ω²)

형태가 된다.

13. Phase Crossover Frequency

Gain Margin을 구할 때는 phase가 -180°가 되는 주파수를 찾는다.

imaginary part가 0이므로

ω(6-ω²) = 0

ω≠0에서

ωpc = √6

이다.

원문 필기에서는 이 주파수에서 Nyquist plot의 실수축 교점을 약

-0.3

으로 계산한다.

14. Gain Margin

음의 실수축 교점의 크기가 0.3이면

GM = 1/0.3
≈ 3.33

dB로 나타내면

GMdB = 20log10(3.33)
≈ 10.45 dB

15. Phase Margin

Phase Margin은 Nyquist plot이 unit circle,

|G(jω)H(jω)| = 1

과 만나는 Gain Crossover Frequency에서 계산한다.

원문 필기에서는

ωgc ≈ 1.233 rad/s

를 얻고, 이때 위상을 약

∠G(jωgc)H(jωgc) ≈ -112.3°

로 계산한다.

따라서 Phase Margin은

PM = 180° + ∠G(jωgc)H(jωgc)
PM ≈ 180° – 112.3°
PM ≈ 67.7°
Gain Margin은 -180° 지점에서의 크기로 구하고, Phase Margin은 |GH|=1 지점에서의 위상으로 구한다.

16. Gain Margin / Phase Margin 정리

항목 조건 계산
Gain Margin ∠GH = -180° GM = 1/|GH|
Gain Margin [dB] Phase Crossover Frequency 20log₁₀(GM)
Phase Margin |GH| = 1 PM = 180° + ∠GH

17. 주요 관계식

항목 내용
Frequency Response s=jω 대입
Nyquist Critical Point -1+j0
Closed-loop 안정 Z=0
Nyquist 관계 Z=P+N
Contour 위 Pole 작은 반원으로 우회
Gain Margin 1/|GH(jωpc)|
Phase Margin 180° + ∠GH(jωgc)

18. 핵심 기억 문장

  1. Nyquist plot은 open-loop transfer function에 s=jω를 대입하여 복소평면에 그린 frequency response이다.
  2. Nyquist 안정도 판정의 기준점은 -1+j0이다.
  3. Open-loop pole이 Nyquist contour 위에 있으면 작은 반원으로 우회한다.
  4. s-plane의 작은 우회 경로는 GH-plane에서 큰 반원으로 mapping될 수 있다.
  5. Nyquist criterion은 open-loop RHP pole P와 -1점 encirclement N을 이용하여 closed-loop RHP pole Z를 판정한다.
  6. Gain Margin은 phase가 -180°가 되는 주파수에서 구한다.
  7. Phase Margin은 |GH|=1이 되는 주파수에서 구한다.
Nyquist Plot, Nyquist Criterion, Frequency Response, Nyquist Contour, Critical Point, -1+j0, Encirclement, RHP Pole, Origin Pole, Gain Margin, Phase Margin, Phase Crossover Frequency, Gain Crossover Frequency
“`

답글 남기기

이메일 주소는 공개되지 않습니다. 필수 필드는 *로 표시됩니다

이 사이트는 Akismet을 사용하여 스팸을 줄입니다. 댓글 데이터가 어떻게 처리되는지 알아보세요.