필터 근사 — 버터워스·베셀·체비셰프

한 줄 요약 — 같은 차수로 평탄함·급준함·파형 보존 셋을 다 가질 수는 없습니다. 무엇을 버릴지 정하는 것이 필터 근사를 고르는 일입니다.


왜 이름이 여럿인가

2차 저역통과의 전달함수는 하나입니다.

H(s)=ωn2s2+2ζωns+ωn2H(s) = \frac{\omega_n^2}{s^2 + 2\zeta\omega_n s + \omega_n^2}
기호
n비트 수 또는 개수number of bits or count

ζ\zeta 하나만 다르게 잡은 것들에 이름이 붙어 있습니다.

이름ζ\zeta (2차)Q무엇을 최적화
베셀0.8660.577군지연이 평탄 (파형 보존)
버터워스0.7070.707통과대역 크기가 평탄
체비셰프 0.5 dB0.5390.928급준함 (리플 허용)
체비셰프 1 dB0.4891.02더 급준
임계 감쇠1.00.5오버슈트 0

같은 회로에 저항·커패시터 값만 다릅니다. 부품이 늘어나는 게 아닙니다.

💡 ζ=0.707\zeta = 0.707 이 4-07에서 계속 나온 이유가 이것입니다 — 버터워스입니다.


셋을 나란히

크기 응답 (통과대역)

∣H(jω)∣=11+(ωωc)2n(버터워스)|H(j\omega)| = \frac{1}{\sqrt{1 + \left(\dfrac{\omega}{\omega_c}\right)^{2n}}} \quad \text{(버터워스)}
기호
omega각주파수angular frequencyrad/s
j허수 단위imaginary unit
n비트 수 또는 개수number of bits or count

버터워스는 ω=0\omega=0 에서 미분이 2n−12n-1 차까지 0 입니다 — 그래서 “최대 평탄(maximally flat)” 이라고 합니다.

차수차단 위 1옥타브(2배)에서
2차−12 dB
4차−24 dB
8차−48 dB

차수당 −6 dB/옥타브 (−20 dB/decade) 는 셋 다 같습니다. 차이는 차단 근처의 모양입니다.

차단 근처 (1.5 fcf_c 에서, 2차 기준)

감쇠
베셀−6.4 dB (가장 무름)
버터워스−7.4 dB
체비셰프 1 dB−9.8 dB (가장 급준)

스텝 응답 — ★ 여기가 갈리는 곳

오버슈트 (2차)정착 시간 (1%)
베셀0.4 %가장 빠름
버터워스4.3 %중간
체비셰프 1 dB10.8 %가장 느림

급준할수록 오버슈트가 큽니다. 예외가 없습니다 — 주파수 영역에서 얻은 것을 시간 영역에서 갚습니다.


★ 군지연 — 파형이 왜 일그러지나

필터가 모든 주파수를 같은 시간만큼 미루면 파형은 통째로 지연될 뿐 모양이 안 변합니다.

τg(ω)=−dϕdω\tau_g(\omega) = -\frac{d\phi}{d\omega}
기호
omega각주파수angular frequencyrad/s
phi위상phase°
d거리 또는 두께distance or thicknessm
g컨덕턴스 (소신호)small-signal conductanceS

τg\tau_g 가 주파수에 따라 다르면 성분마다 다른 시간에 도착해 파형이 일그러집니다.

통과대역 군지연 변화
베셀±1 % 이내
버터워스±20 %
체비셰프 1 dB±50 % 이상

베셀은 군지연이 평탄하도록 설계된 것이고, 그 대가로 차단이 무릅니다.

언제 이게 중요한가

신호군지연이
정현파 하나상관없음 (성분이 하나)
펄스·계단결정적 — 여러 성분의 시간 관계가 모양
심전도의 QRS결정적 — 모양이 진단 근거
오디오논쟁이 있음 (귀가 위상에 둔감)
스펙트럼 분석상관없음 (크기만 봄)

⚠️ 측정 회로에서 파형의 모양을 본다면 베셀입니다. 이 자료가 다루는 신호 측정 대부분이 그렇습니다.


⚖️ 트레이드오프

1. 셋 중 무엇을

무엇을 보나?파형의 모양펄스 폭, 상승 시간, QRS베셀군지연 평탄. 차단이 무른 것을 차수로 벌충통과대역의 크기주파수 응답을 재는 것버터워스가장 흔한 기본값저지대역만 중요특정 주파수를 죽이는 것체비셰프리플을 감수하고 급준하게

의심스러우면 버터워스입니다. 어느 쪽으로도 크게 나쁘지 않은 절충입니다.

2. 급준함을 차수로 살 것인가 근사로 살 것인가

같은 저지대역 감쇠를 얻는 두 가지 방법이 있습니다.

차수를 올림체비셰프로 바꿈
부품늘어남그대로
통과대역평탄 유지리플 생김
오버슈트조금 늘어남많이 늘어남
부품 오차 민감도조금크게 늘어남 (Q 가 높아짐)

Q 가 높아지면 부품 오차에 민감해집니다. 체비셰프 1 dB 2차의 Q 는 1.02 로 버터워스(0.707)보다 44 % 높고, 5 % 저항을 쓰면 리플이 설계값에서 크게 벗어납니다.

3. 안티에일리어싱에는

13-06의 요구는 “나이퀴스트에서 충분히 감쇠” 입니다. 통과대역의 평탄함이나 파형은 그 다음 문제입니다.

그런데 오버샘플링을 쓰면 전이대역이 넓어져 2차 버터워스로 충분해집니다. 급준한 필터가 필요한 상황 자체가 줄어듭니다.


🔌 실제 쓰임새

값을 어떻게 정하나

살렌-키 저역통과에서 R1=R2=RR_1=R_2=R, C1C_1, C2C_2 이면

fc=12πRC1C2,Q=C1C22C2f_c = \frac{1}{2\pi R\sqrt{C_1C_2}}, \qquad Q = \frac{\sqrt{C_1C_2}}{2C_2}
기호
C_11번 정전용량capacitance 1F
C_22번 정전용량capacitance 2F
f_c차단 주파수cutoff frequencyHz
Q전하 / 품질 계수charge / quality factor
R저항resistanceΩ

C1/C2C_1/C_2 비가 Q 를 정합니다.

C1C2=4Q2\frac{C_1}{C_2} = 4Q^2
기호
C_11번 정전용량capacitance 1F
C_22번 정전용량capacitance 2F
Q전하 / 품질 계수charge / quality factor
근사QC1/C2C_1/C_2
베셀0.5771.33
버터워스0.7072
체비셰프 1 dB1.024.2

버터워스는 C1=2C2C_1 = 2C_2 — 외우기 쉽습니다.

고차는 2차를 겹칩니다

4차 버터워스는 Q 가 다른 2차 두 개입니다.

4차 버터워스Q
1단0.541
2단1.31

같은 2차 두 개를 잇는 것이 아닙니다. Q 0.707 짜리 둘을 이으면 −6 dB 가 되어 차단 주파수가 어긋납니다.

⚠️ Q 가 높은 단을 뒤에 둡니다. 앞에 두면 그 단의 피킹이 뒷단으로 들어가 슬루레이트슬루레이트 slew rate출력이 오를 수 있는 최고 속도로, 이걸 넘기면 정현파가 삼각파로 잘려 나옵니다.fmax=SR2πVpf_{max} = \frac{\mathrm{SR}}{2\pi V_p}자세히 보기 →나 포화에 먼저 걸립니다.

부품값이 E24 에 없을 때

계산값이 표준값에 없으면 E 계열 반올림을 해야 하는데, Q 가 높을수록 그 반올림이 아프게 옵니다. 계산기의 E 계열 기능이 오차를 %로 보여 주므로, 그 값을 넣고 다시 시뮬레이션해 확인하세요.


⚠️ 이상과 실제의 차이

이상실제왜
근사를 고르면 그 응답이 나옴부품 오차로 Q 가 흔들림5 % 저항이면 Q 도 5 %
4차는 2차 둘Q 가 다른 2차 둘극점 위치가 다름
체비셰프가 항상 급준통과대역 리플을 감수한 것공짜가 아님
베셀은 성능이 나쁨파형을 지킴목적이 다름
차수를 올리면 좋아짐오버슈트·정착시간이 나빠짐주파수↔시간
OP앰프는 이상적GBW 가 부족하면 Q 가 변함7-06

⚠️ OP앰프의 GBW 가 부족하면 설계한 Q 가 안 나옵니다. 살렌-키에서 필요한 GBW 는 대략 100 Q fc100\,Q\,f_c 입니다. 10 kHz·Q 1.3 이면 1.3 MHz 이상이 필요합니다.


시뮬레이션으로 확인

베셀·버터워스·체비셰프 — ζ 하나가 다를 뿐 회로도R1447Ωζ 를 정합니다 (R = 2ζ√(L/C))L110mHV1C1100nFINMIDOUTGND
베셀·버터워스·체비셰프 — ζ 하나가 다를 뿐워크벤치에서 열기 →
# examples/04-ac-frequency/filter-approximations.yaml
meta:
  name: 베셀·버터워스·체비셰프 — ζ 하나가 다를 뿐
  lesson: docs/04-ac-frequency/08-filter-approximations.md
  note: |
    같은 RLC 2차 저역통과입니다. R1 하나로 ζ 를 바꿉니다.

        ζ = (R/2)·√(C/L),   즉  R = 2ζ·√(L/C)

    L = 10 mH, C = 100 nF 이면 √(L/C) = 316 Ω 이므로

        R = 548 Ω  → ζ = 0.866  베셀         첫 봉우리 1.004 @ 199 µs
        R = 447 Ω  → ζ = 0.707  버터워스     첫 봉우리 1.043 @ 140 µs
        R = 309 Ω  → ζ = 0.489  체비셰프 1dB  첫 봉우리 1.172 @ 114 µs
        R = 632 Ω  → ζ = 1.0    임계 감쇠     오버슈트 없음

    f_n = 1/(2π√(LC)) = 5.03 kHz.

    R1 을 밀어 보세요. 급준할수록(ζ 가 작을수록) 오버슈트가 큽니다 —
    주파수 영역에서 얻은 것을 시간 영역에서 갚습니다. 예외가 없습니다.

    기본값은 버터워스(447 Ω)입니다.
components:
  - {ref: V1, type: vpulse, params: {v1: 0, v2: 1, td: 0, tr: 1u, tf: 1u, pw: 5m, per: 10m}}
  - {ref: R1, type: r, value: 447, label: ζ 를 정합니다 (R = 2ζ√(L/C)), sweep: {min: 100, max: 900, log: false}}
  - {ref: L1, type: l, value: 10m}
  - {ref: C1, type: c, value: 100n}
nets:
  IN:  [V1.p, R1.1]
  MID: [R1.2, L1.1]
  OUT: [L1.2, C1.1]
  GND: [V1.n, C1.2]
analyses:
  - {type: tran, tstop: 1m, tstep: 200n}
expect:
  # 버터워스(ζ=0.707)의 첫 봉우리는 140 µs 에서 1.043
  - {measure: v, node: OUT, at: 140u, analysis: tran, value: 1.043, tol: 2%}
  # 정상상태는 1 V
  - {measure: v, node: OUT, at: 900u, analysis: tran, value: 1.0, tol: 2%}

돌려 보면 140 µs 에서 첫 봉우리가 1.0433 — 오버슈트 4.33 % 입니다. 이론값과 소수점 둘째 자리까지 맞습니다.

오버슈트=exp⁡(−πζ1−ζ2)=exp⁡(−π×0.7070.707)=e−π=4.32 %\text{오버슈트} = \exp\left(\frac{-\pi\zeta}{\sqrt{1-\zeta^2}}\right) = \exp\left(\frac{-\pi\times0.707}{0.707}\right) = e^{-\pi} = 4.32\ \%
기호
e잡음 전압 밀도noise voltage densityV/√Hz

ζ=1/2\zeta = 1/\sqrt2 일 때 지수가 정확히 −π-\pi 가 되는 것이 버터워스의 깔끔한 점입니다.

★ R1 을 밀어 보세요. R=2ζL/CR = 2\zeta\sqrt{L/C} 이고 L/C=316 Ω\sqrt{L/C} = 316\ \Omega 입니다.

R1ζ무엇첫 봉우리
632 Ω1.0임계 감쇠오버슈트 없음
548 Ω0.866베셀1.004
447 Ω0.707버터워스1.043
309 Ω0.489체비셰프 1 dB1.172
150 Ω0.24(설계 아님)심하게 울림

같은 부품 세 개로 셋 다 만듭니다 — 저항값 하나만 다릅니다. 필터 근사가 “다른 회로” 가 아니라 “같은 회로의 다른 값” 이라는 것이 여기서 보입니다.


연습문제

문제 1. 심전도 프론트엔드의 저역통과(150 Hz)를 설계합니다. 버터워스와 베셀 중 어느 것을 쓰겠습니까?

연습 4-08-1 — 심전도 150 Hz 저역통과, 버터워스냐 베셀이냐 회로도RA110kΩRA210kΩCA1150nFRB110kΩRB210kΩCB1122nFQ = 0.577 (베셀)V1CA275nFCB291.9nFVP5VVN-5VU1OPA320U2OPA320VCCVCCVCCVEEVEEVEEINA1A2BUTTB1B2BESSGND
연습 4-08-1 — 심전도 150 Hz 저역통과, 버터워스냐 베셀이냐워크벤치에서 열기 →
풀이

베셀입니다.

심전도에서 판단 근거는 파형의 모양입니다 — QRS 폭, ST 분절의 높이, T 파의 대칭성. 이것들은 여러 주파수 성분의 상대적 시간 관계로 정해집니다.

통과대역 군지연 변화QRS 에 미치는 영향
버터워스±20 %고주파 성분이 늦게 도착 → QRS 가 퍼지고 비대칭
베셀±1 %통째로 지연될 뿐 모양 유지

대가: 베셀은 차단이 무릅니다. 150 Hz 에서 −3 dB 를 맞추면 300 Hz 에서 버터워스보다 몇 dB 덜 감쇠합니다.

벌충: 차수를 하나 올립니다. 4차 베셀이 2차 버터워스보다 급준하면서도 군지연이 평탄합니다. 부품이 늘지만 파형을 지키는 값입니다.

실제 심전도 규격(AAMI EC11)이 저역통과의 위상 응답을 규정하는 이유가 이것입니다.

문제 2. 4차 버터워스를 살렌-키 두 단으로 만듭니다. Q 0.541 단과 Q 1.31 단 중 어느 것을 앞에 두겠습니까?

exercises/04-08-stage-order.yaml
풀이

Q 가 낮은 단(0.541)을 앞에 둡니다.

왜:

Q 1.31 단은 차단 근처에서 피킹이 있습니다.

∣H∣peak=Q1−1/4Q2=1.311−0.146=1.42(+3.0 dB)|H|_{peak} = \frac{Q}{\sqrt{1 - 1/4Q^2}} = \frac{1.31}{\sqrt{1-0.146}} = 1.42 \quad(+3.0\ \text{dB})
기호
e잡음 전압 밀도noise voltage densityV/√Hz
k결합 계수 또는 상수coupling coefficient or constant
Q전하 / 품질 계수charge / quality factor

그 단을 앞에 두면 차단 근처의 큰 신호가 3 dB 더 커진 채로 뒷단에 들어갑니다. 그러면

문제왜
포화출력 범위를 먼저 넘음
슬루레이트큰 진폭이 빠르게 변함
왜곡위 둘의 결과

Q 가 낮은 단을 먼저 두면 그 성분이 이미 줄어든 뒤 고Q 단에 들어갑니다.

반대 논리도 있습니다 — 잡음.

고Q 단이 뒤에 있으면 그 단의 잡음이 앞단 이득으로 안 나눠집니다. 하지만 필터의 이득은 대개 1 이라 잡음 지수의 캐스케이드 관점에서 큰 차이가 없고, 포화 쪽 위험이 훨씬 큽니다.

결론: 낮은 Q 먼저. 신호가 큰 경우일수록 확실합니다.

문제 3. 체비셰프 1 dB 4차 필터를 5 % 저항으로 만들었더니 통과대역 리플이 2.5 dB 나왔습니다. 왜?

exercises/04-08-chebyshev-tolerance.yaml
풀이

Q 가 부품 오차에 민감하기 때문입니다.

살렌-키에서 Q∝C1/C2Q \propto \sqrt{C_1/C_2} 이고, 저항 비도 함께 들어갑니다. 부품이 각각 5 % 틀리면 Q 는 RSS 로 약 7 % 흔들립니다.

리플은 Q 에 훨씬 민감합니다.

4차 체비셰프 1 dB 의 두 번째 단은 Q ≈ 3.56 입니다. 이 단의 피크는

∣H∣peak≈Q=3.56(+11 dB)|H|_{peak} \approx Q = 3.56 \quad(+11\ \text{dB})
기호
e잡음 전압 밀도noise voltage densityV/√Hz
k결합 계수 또는 상수coupling coefficient or constant
Q전하 / 품질 계수charge / quality factor

Q 가 7 % 올라 3.81 이 되면

20log⁡(3.81/3.56)=+0.6 dB20\log(3.81/3.56) = +0.6\ \text{dB}

두 단이 각각 이렇게 어긋나고 서로 상쇄되지 않으면 리플이 배로 늘어납니다.

대책:

방법효과
1 % 부품Q 흔들림 1.4 % → 리플 거의 설계대로
차수를 올리고 리플을 줄임 (0.1 dB 5차)Q 가 낮아져 둔감해짐
버터워스로 바꿈Q 최대 1.31 로 훨씬 둔감
디지털로 옮김계수 오차 0 (13-11)

여기서 배울 것: 급준함을 근사로 사면 부품 정밀도로 갚습니다. 체비셰프를 쓰려면 1 % 이하 부품이 사실상 필수이고, 그 값이면 차수를 하나 올린 버터워스가 대개 더 쌉니다.


이전 ← 07 2차 시스템 — 오버슈트와 정착 시간 | 다음 → 5장. 다이오드