시간 영역에서 복소평면으로
Z-변환은 수열 x_n을 복소 변수의 함수 X(z)로 매핑합니다:
X(z) = Σ_{n=−∞}^{∞} x_n · z^{−n}
변수 z는 복소평면을 매개변수화합니다. 이 평면의 서로 다른 영역은 필터의 서로 다른 정성적 동작에 해당합니다.
기하학적 영역
| 영역 | z | 동작 | ||
|---|---|---|---|---|
| 단위원 내부 | < 1 | 안정적 극점: 감쇠 응답 | ||
| 단위원 | = 1 | 주파수 축: z = e^{i2πf} | ||
| 단위원 외부 | > 1 | 불안정 극점: 성장 응답 |
단위원은 이산시간 안정성에서 연속시간(Laplace) 안정성의 허수축과 같은 역할을 합니다.
Laplace 변환과의 관계
연속시간 시스템의 경우 Laplace 변환은 변수 s를 사용합니다. 허수축 s = iω는 주파수 응답이 존재하는 곳입니다. 안정성: 극점은 Re(s) < 0(좌반평면)을 가져야 합니다.
쌍선형 변환은 s → z로 매핑합니다: z = (1 + Ts/2) / (1 − Ts/2). 이것은 좌반평면을 단위원 내부로 매핑합니다 — '좌반평면 안정'의 기하학적 번역입니다.
쌍선형 변환을 등각 사상으로
쌍선형 변환 z = (1 + Ts/2) / (1 − Ts/2)는 뫼비우스 변환 — 복소평면의 등각(각도 보존) 사상입니다.
주요 기하학적 성질:
- s = iω(허수축)을 |z| = 1(단위원)로 매핑
- Re(s) < 0(좌반평면)을 |z| < 1(단위원 내부)로 매핑
- Re(s) > 0(우반평면)을 |z| > 1(단위원 외부)로 매핑
- 주파수 왜곡: 매핑 ω → f는 비선형 — ω_analog = (2/T)·tan(πf_digital)
이 왜곡은 높은 주파수를 Nyquist 지점 쪽으로 압축합니다. 설계자는 쌍선형 변환을 적용하기 전에 아날로그 명세를 사전 왜곡하여 이를 보정합니다.
Butterworth 극점: 원 궤적
Butterworth 필터는 통과대역에서 최대 평탄도를 달성하기 위해 아날로그 극점을 s-평면에 있는 반지름 ω_c인 원 위에 배치합니다.
N차 Butterworth 필터의 경우 극점은 다음 위치에 있습니다:
s_k = ω_c · e^{iπ(2k+N−1)/(2N)} k = 0, 1, …, N−1의 경우
이것은 반지름 ω_c인 원의 좌반에 균등하게 배치합니다. (우반의 극점은 불안정할 것입니다; 좌반평면의 극점만 유지됩니다.)
원 궤적이 왜 최대 평탄 통과대역을 만듭니까?
Butterworth 다항식 |B_N(jω)|² = 1 + (ω/ω_c)^{2N}은 모든 극점이 |s| = ω_c에 있습니다. 동일 반지름 제약은 모든 극점이 ω = ω_c에서 크기 응답에 동등하게 기여함을 의미합니다. 최대 평탄도 정리: 이 원 위의 극점을 가진 모든 N차 다항식 중에서 Butterworth 다항식은 ω = 0에서 0과 같은 도함수를 가장 많이 가집니다.
Chebyshev 극점: 타원 궤적
Chebyshev 극점은 s-평면의 타원 위에 있습니다(원이 아님). 타원은 리플 매개변수 ε에 의해 결정되는 장축과 단축을 가집니다. 동일 리플 통과대역은 Chebyshev 다항식의 등리플 성질에서 나옵니다.
Elliptic 극점: 타원 함수 궤적
Elliptic(Cauer) 필터 극점도 타원 위에 있습니다 — 하지만 극점과 영점 모두 주파수 응답에 기여합니다. 영점은 허수축 위에 있습니다(통지대의 유한 감쇠 극점). 타원 함수 매핑은 양쪽 대역에서 동일 리플을 달성하기 위해 영점을 최적으로 배포합니다.
Butterworth 극점 위치 계산
4차 Butterworth 필터의 경우 ω_c = 1(정규화), 극점은 다음 위치에 있습니다:
s_k = e^{iπ(2k+3)/8} k = 0, 1, 2, 3의 경우
k=0: s₀ = e^{i3π/8} (좌반평면)
k=1: s₁ = e^{i5π/8} (좌반평면)
k=2: s₂ = e^{i7π/8} (좌반평면)
k=3: s₃ = e^{i9π/8} (좌반평면)
이 네 개의 극점은 단위원 위의 동일 각도 간격에 있으며, 모두 음의 실수부를 가집니다(좌반평면).
단위원으로부터 극점까지의 거리
이론적 안정성은 |p| < 1을 필요로 합니다. 실제로는 두 가지 추가 우려가 발생합니다.
안정성 여유
IIR 필터의 안정성 여유는 모든 극점에서 단위원까지의 최소 거리입니다: min_k (1 − |p_k|).
|p| = 0.99의 극점은 이론적으로 안정적이지만 1%의 여유만 남깁니다. 유한 정밀도 산술(계수 표현의 반올림 & 반올림 오차의 누적)은 효과적으로 극점을 움직일 수 있습니다. 계수 양자화가 극점을 0.99에서 1.001로 이동하면 필터가 불안정해집니다.
기하학적 결과
단위원에 매우 가까운 극점은 매우 뾰족한 주파수 응답 피크를 생성합니다 — 좁은 대역폭 공진기. 하지만 좁은 공진기는 높은 정밀도를 필요로 합니다: 작은 계수 오류는 피크 주파수를 크게 이동합니다.
기하학적 트레이드오프: 피크 날카로움 ∝ 1 / (1 − |p|). |p| → 1일수록 날카로움 → ∞이지만 안정성 여유 → 0이고 계수 오류에 대한 민감도 → ∞.
2차 섹션
단일 다항식으로 구현된 고차 IIR 필터는 수치적으로 민감합니다 — 단일 계수를 반올림하면 많은 극점을 이동할 수 있습니다. 표준 솔루션: 2차 섹션(biquads) 캐스케이드로 구현하며, 각각 하나의 켤레 극점 쌍과 하나의 켤레 영점 쌍만 가집니다. 한 섹션의 오류는 다른 섹션의 극점을 교란할 수 없습니다.