Audio Processing38 오디오 필터(22)- 바이쿼드 필터의 주파수 응답 바이쿼드 필터의 주파수 응답 지난 글에서 바이쿼드 필터에 대해 설명한적이 있습니다. 지난 글에서는 바이쿼드 필터의 주파수 특성을 GNU Octave의 함수를 이용해서 검토했습니다. 이번에는 FFT나 Octave 함수를 사용하지 않고 바이쿼드 필터의 주파수 특성을 구하는 내용을 다루도록 하겠습니다. 바이쿼드 필터의 전달함수 바이쿼드 필터의 전달함수는 다음과 같이 정리됩니다. $$H(z)=\frac{b_0+b_1z^{-1}+b_2z^{-2}}{1+a_1z^{-1}+a_2z^{-2}}$$ z를 $j\omega$로 치환하면 디지털 필터 전달함수를 복소수 공간의 전달함수 로 변환이 됩니다. $$H(e^{j\omega})=\frac{b_0+b_1e^{-j\omega}+b_2e^{-j2\omega}}{1+a_1e^{.. 2021. 11. 11. 오디오 필터(21)-FIR 카이저 윈도우(Kaiser window) 카이저 윈도우로 FIR 필터 구하기 FIR 필터를 구하는 방법 중에 윈도우 방법에서 카이저 윈도우(Kaiser window) 함수를 만들어서 기존의 윈도우 함수 대신 사용하는 것입니다. 카이저 윈도우는 필터를 만들기 위한 설계요건에 따라 윈도우 형태가 변경되는 특징이 있습니다. 이는 마치 IIR 필터의 통과 필터를 설계하는 방법을 FIR 필터로 옮긴것 같습니다. 카이저 윈도우 구하기 카이저 윈도우는 구현되는 필터의 조건에 따라서 윈도우를 형성하는 파라미터를 구해줘야 합니다. 윈도우의 통과 대역 주파수 $f_{pass}$와 저지 대역 주파수 $f_{stop}$과 통과 대역의 리플(ripple) $\delta_{pass}$와 저지 대역의 리플 $\delta_{stop}$, 샘플 주파수 $f_s$를 설정합니다.. 2021. 11. 1. 오디오 필터(20)-콘보루션을 이용한 FIR 필터의 실행 콘보루션과 FIR 필터 FIR 필터의 계수를 구하는 방법으로 윈도우 방식과 주파수 샘플링 방식에 대해 설명한 적이 있습니다. IIR 필터를 실행한 것과 같이 이번에도 오디오 신호를 입력하고 FIR 필터를 통과해서 출력하는 방법을 알아보겠습니다. 먼저 FIR 필터 실행에 기반이 되는 콘보루션에 대해서도 알아보겠습니다. 콘보루션 콘보루션(convolution) 또는 합성곱은 하나의 신호 함수에 가중 함수를 곱해서 그 결과는 얻는 방법입니다. 기본적인 방법은 하나의 함수를 시간 반전하고 $\tau$ 만큼 시간 이동해서 이를 적분하는 것입니다. 이때 함수의 순서는 중요하지 않습니다. 어느 함수를 반전하던지 관계없습니다. $$(f\ast g)(t)=\int_{\tau = -\infty}^{\infty}f(\tau.. 2021. 10. 29. 오디오 필터(19)-주파수 샘플 방법으로 FIR 필터 구하기 주파수 샘플링 주파수 샘플링 방법으로 FIR 필터를 구하는 방법은 참으로 기막힌 기법입니다. 시간축에서 계속 변화하는 아날로그 신호를 샘플링으로 디지털 신호로 변화하는 방법이 있었습니다. 이 방법을 주파수 축에서 적용하는 해서 필터를 구하는 방법이 주파수 샘플링 방법입니다. 이 글에서는 주파수 샘플링 방법으로 원하는 특성의 필터를 쉽게 구하는 방법에 대해서 알아보겠습니다. 사실 이 방법이 실제로 필터를 구하기 더 쉽습니다. 먼저 주파수 축에서 구현될 필터 모양을 결정해야 합니다. 수식으로 해도 되고 실제 원하는 필터의 모양에서 필터의 크기에 관한 정보를 수집해야 합니다. 그리고 주파수 축에서 각주파수로 [0, 2$\pi$] 사이에 필터의 모양에 따라 N개의 샘플을 등간격으로 채취합니다. 이를 H(k)로 .. 2021. 10. 29. 오디오 필터(18) - 윈도우 방법을 이용한 FIR 필터 구하기 윈도우 방법을 이용한 FIR 필터 구하기 지난 글에서 이상적인 필터의 모양을 기반으로 필터 $h(n)$을 구하는 방법을 알아봤습니다. 그런데 $h(n)$은 시간이 -t이고 주파수가 $-\omega $인 조건까지 사용해야 하기 때문에 공학적으로 사용할 수 없습니다. 이를 비 인과성(non-causalty)라고 합니다. 비 인과성은 현재, 과거, 미래의 모든 신호가 필요합니다. -t와 $-\omega $는 미래의 정보입니다. 그런데 우리는 미래의 시간을 만들 수 없습니다. 그러니 인과성(causalty)만 고려해야 합니다. 즉 현재와 과거의 신호만으로 공학적으로 의미 있는 신호처리를 해야 합니다. 이 글에서는 지난 글의 비인과성 조건의 필터를 인과성 조건으로 변경해서 저역 통과 필터를 만들어 보겠습니다. S.. 2021. 10. 29. 오디오 필터(17)-FIR 필터 기반의 통과 필터 구하기 FIR 필터 구하기 FIR 필터는 다른글에서 다룬 신호의 궤환이 있는 IIR 필터와 다르게 궤환이 없이 유한한 필터 계수만으로 구현되는 필터로 선형 위상을 가지는 특징이 있습니다. 다만 IIR 필터 보다 연산량이 많은 단점이 있습니다. 이 글과 다른 글에서 FIR을 이용한 필터를 만들고 실행하는 방법을 알아보겠습니다. 이상적인 필터의 특징을 이용한 필터 계수 구하기 윈도우 방식을 이용한 FIR 필터 구하기 주파수 샘플링 방식을 이용한 FIR 필터 구하기 콘보루션으로 FIR 필터 실행하기 이상적인 필터 지난 글에서 이상적인 필터와 실제 구현 가능한 필터 형태에 대해 설명했습니다. FIR(finity impulse response) 필터는 이상적인 필터를 출발점으로 해서 구현 가능한 형태로 유도합니다. 유도.. 2021. 10. 29. 오디오 필터(16)-IIR 필터 실행, 오디오 볼륨 처리 IIR 필터 다루기 이전의 글에서 디지털 필터의 전달 함수를 위한 계수들을 구하는 방법을 알아봤습니다. 이때 디지털 필터 전달 함수의 형태에 따라 IIR 또는 FIR가 있으며 형태가 조금 다릅니다. 이에 관한 내용은 다른 글에서 다루겠습니다. 이글에서는 GNU Octave에서 오디오를 다루는 방법, 그리고 이를 이용해서 오디오의 볼륨을 조정하는 방법, 마지막으로 디지털 필터를 적용해서 오디오 신호를 처리하고 출력하는 방법을 알아보겠습니다. 유투브를 보시면 오디오 파일을 6초 정도로 짧게 만들어 실험에 사용하는 방법과 GNU Octave로 오디오 파일을 출력 하는 방법 까지 설명해 놨습니다. IIR 필터 실행 방법 IIR 필터는 다음과 같이 처리됩니다. 입력 신호를 받아서 필터 계수를 적용한 IIR 필터에.. 2021. 10. 29. 오디오 필터 실험(15)-원시 데이터 얻기 오디오 필터 실험 여러 글에서 디지털 필터의 계수를 구하는 방법을 알아봤습니다. 이 중에서 피킹 필터와 저역 통과 필터를 실제 오디오 신호에 적용해서 기록하는 과정을 몇 개의 글에 나눠서 설명하겠습니다. 그중에서 첫 번째 내용으로 오디오 파일을 읽고, 다른 형식으로 데이터를 추출하고, 오디오 파일을 재생하는 프로그램의 사용 방법을 먼저 다뤄보겠습니다. OcenAudio의 설치 오디오 재생 프로그램으로 유명한 프로그램으로 알송, zaudio, fooba2000등이 있고, 오디오 편집과 효과 추가가 가능한 프로그램으로 audition, goldwave등과 프로용으로는 cubase, potool, logic 등이 있습니다. 편집이 가능한 오디오용 프로그램은 대부분 유료입니다. ocenaudio는 무료이면서, .. 2021. 10. 29. 오디오 필터(14)-IIR 버터워스 고역 통과 필터(Butterworth High Pass Filter) 설계 고역 통과 필터 설계 다른 글에서 사용했던 유사한 방법으로 Butterworth HPF를 설계하겠습니다. 이번엔 필터의 Q를 이용하지 않고, 필터의 차단 주파수와 감쇄량을 이용하는 방법을 알아보겠습니다. 다음과 같은 조건으로 고역 통과 필터를 만들겠습니다. 통과대역의 특성으로 차단 주파수와 통과대역 주파수$f_{c}=F_{p}=1000Hz$로 하겠습니다. 그러면 $F_{p}$ 에서 감쇄량은 -3dB가 됩니다. 그리고 저지 대역의 주파수 $F_{s}=80Hz$이며 감쇄량이 -80dB입니다. 이때 샘플링률 $f_{s}=44100Hz$, 샘플 구간 $T_{s}=1/f_{s}$입니다. 필터의 전달 함수 구하기 먼저 필터의 조건으로 다음과 같이 필터를 계산하기 위한 사전 조치를 취합니다. $$A_{p}=10^{\.. 2021. 10. 29. 오디오 필터(13)-IIR 버터워스 저역 통과 필터(Butterworth Low Pass Filter) 설계 저역통과 필터 오디오 필터로 가장 유용한 형태는 버터워스 필터입니다. 통과대역의 특성이 가장 평탄하기 때문입니다. 물론 이를 발전시킨 필터도 있습니다만, 사용성이나 실용적으로 구현하기에도 버터워스 필터는 오디오에 적합합니다. 다른 글에서는 필터의 Q를 이용해서 저역 통과 필터의 계수를 구했습니다. https://medialink.tistory.com/74?category=958130 이글에서는 필터의 일반적인 규격 조건을 통해 버터워스의 s-영역(s-domain)(라플라스 영역)에서 구하고 이를 z-영역(z-domain)으로 변환해서 필요한 디지털 계수를 구하는 방법을 알아보겠습니다. 그리고 필터의 특성을 GNU Octave에서 알아보도록 하겠습니다. 버터워스 필터를 구하는 방식 아날로그 신호 영역인 s.. 2021. 10. 29. 오디오 필터(12)-IIR 노치 필터(Notch Filter, NF) 노치 필터 노치 필터(notch filter)는 대역 저지 필터(band rejection filter)의 하나로 매우 좁은 주파수 대역의 신호를 차단할 목적으로 사용합니다. 이 필터는 오디오 신호처리에서는 적응 필터에서 하울링 억제를 위해 사용됩니다. 또는 특정 주파수의 노이즈를 제거할 때도 사용합니다. 노치 필터를 만드는 방법도 여러 가지가 있습니다. 특히 피킹 필터를 대신 사용해도 됩니다. 필터의 계수는 다음과 같이 생성하면 됩니다. 여기서 Q에 따라 노치 필터가 형성하는 대역이 달라집니다. 이 부분에 대한 내용은 다음을 참조해 주세요 https://medialink.tistory.com/73?category=958130 https://medialink.tistory.com/77?category=9.. 2021. 10. 29. 오디오 필터(11)-IIR 전대역 통과필터(All Pass Filter, APF) 전대역 통과 필터 전대역 통과 필터(all pass filter, APF)는 LPF등과 다르게 모든 신호를 그냥 통과시킵니다. 그러나 필터 주파수를 기준으로 위상을 변경시킵니다. 그래서 위상 특성의 수정을 위해서 주로 사용합니다. 이를 이용해서 인공 잔향기를 만들기도 합니다. APF의 필터 계수 APF의 필터 계수는 다음과 같이 구합니다. 이때 Q에 따라서 위상의 변화 정도가 결정됩니다. 이부분에 대햔 내용은 아래를 참조해 주세요 https://medialink.tistory.com/73?category=958130 https://medialink.tistory.com/77?category=958130 $$h=\frac{sin(\omega_0)}{2Q}$$ $$b_{0}=a_{2}=1-h$$ $$b_{1.. 2021. 10. 29. 오디오 필터(10)-IIR 대역통과 필터(Band Pass Filter, BPF) 대역 통과 필터 대역통과 필터(band pass filter, BPF)는 특정 주파수 사이의 신호만 통과하는 필터입니다. 대역통과 필터는 오디오 신호의 대역을 조정해서 불필요한 잡음의 영향을 줄이거나, 음향적인 효과를주기 위해 사용합니다. 특히 스피커를 위한 크로스오버에서도 활용됩니다. 이부분은 다시 다루기로 하겠습니다. 필터 합성 이 필터를 만드는 방법도 여러 가지가 있지만, 다음과 같이 저역통과 필터와 고역통과 필터를 조합해서 구현하기도 합니다. 그러니 두 필터를 형성하는 필터 계수를 각각 구하고 조합하는 방법을 사용합니다. 이때 필터의 순서는 대역 통과 필터를 구현하는데 아무런 관계가 없습니다. 두 필터를 합성하는 방법은 두필터의 전달 함수를 곱하면됩니다. 그런데 다른 글에서 전달함수의 결과로 필터.. 2021. 10. 29. 오디오 필터(9)-IIR 이퀄라이저(Equalizer), 하이 쉘빙 필터(High Shelving Filter, HSF) 하이 쉘빙 필터 하이 쉘빙 필터(hight shelving filterm, HSF)는 다른 글에서 설명한 로 쉘빙 필터(low shelving filter, LSF)와 반대로 고음 영역의 이득을 조정합니다. 하아 쉐빙 필터도 조정하고자 하는 주파수, 이득 드리고 필터의 기울기를 인자로 합니다. 필터를 구하는 내용은 두 필터가 동일하며 필터 계수는 다음과 같습니다. 이 부분에 대한 내용은 아래를 참조해 주세요 https://medialink.tistory.com/73?category=958130 https://medialink.tistory.com/77?category=958130 $$\omega_{0}=\frac{2\pi f}{f_{s}}$$ $$A=10^{\frac{G}{40}}$$ $$c_{0}=co.. 2021. 10. 29. 오디오 필터(8)-IIR 이퀄라이저(Equalizer), 로 쉘빙 필터(Low Shelving Filter, LSF) 로 쉘빙 필터 이퀄라이저에서 사용되는 로 쉘빙 필터(low shelving filter, LSF)는 저역 부분의 이득을 조정할 목적으로 사용되는 이퀄라이저용 필터입니다. 필터 주파수를 기준으로 낮은 주파수 부분의 신호를 증폭하거나 감소합니다. 로 쉘빙 필터는 조정하고자 하는 주파수와 조정하는 이득 그리고 필터의 기울기를 인자로 합니다. 피킹 필터와 다르게 Q 대신 slope을 이용해서 조정합니다. 그런데 두 인자는 변환이 가능합니다. 관련 내용은 다른 글을 참조해 주세요. 필터를 구하는 기본 내용은 피킹 필터와 같으니 피킹 필터 부분으로 참조하면 됩니다. 필터 계수는 다음과 같습니다. 여기서 $f$는 필터의 주파수, $f_{s}$는 샘플링 주파수, G는 dB로 표시되는 필터의 이득, slope는 필터의 .. 2021. 10. 29. 오디오 필터(7)-IIR 이퀄라이저(Equalizer), 피킹 필터(Peaking Filter, PF) 피킹 필터 오디오 시스템에서 사용 용도가 많은 필터가 아마도 피킹 필터(peaking filter)입니다. 이 필터는 중심 주파수, 이득, 그리고 품질 팩터 Q로 조정이 가능합니다. 이 필터는 특정 주파수 근처의 신호를 증가 또는 감소하여 음향적인 특징을 조정할 때 사용합니다. 조정의 결과도 매우 효과적이어서 피킹 필터의 활용도는 매우 높습니다. 이 필터는 필터 중심 주파수를 기준으로 근처의 주파수의 이득에 영향을 미칩니다. 이는 필터 중심 주파수 부근에 다른 주파수에도 영향을 줍니다. 이 필터의 품질 팩터 Q에 따라서 영향을 주는 주파수 범위가 달라집니다. Q가 크다는 것은 매우 좁은 범위의 주변 주파수 영역에만 영향을 미친다는 의미 합니다. 필터 계수는 다음과 같습니다. 이부분에 대한 내용은 아래를 .. 2021. 10. 29. 오디오 필터(6)-이퀄라저(Equalizer)의 품질팩터(Quality Factor, Q) 품질 팩터 이퀄라이저에 사용되는 필터는 품질 팩터(Q)를 필터의 인자로 사용합니다. Q는 필터가 어느 정도 범위까지 영향을 미치는지 결정합니다. 일반적으로 Q 작을수록 필터가 주파수 축에서 영향을 미치는 범위가 큽니다. Q는 영향을 미치는 주파수에 범위 의해 다음과 같이 정의됩니다. 아래는 피킹 필터의 Q를 설명하는 그림입니다. $$Q=\frac{f_{c}}{BW}=\frac{f_{c}}{f_{1}-f_{0}}$$ $$f_{0,1}=f_{c}\left (\sqrt{1+\frac{1}{4Q^{2}}}+\frac{1}{2Q} \right )$$ Slope 반면 쉘빙 필터는 Q보다는 slope이라는 용어를 사용합니다. 이는 쉘빙 필터의 모양에서 필터가 기울어진 정도를 표시하려고 만든 개념입니다. Q와 slop.. 2021. 10. 29. 오디오 필터(5)-이퀄라이저와 그래픽 이퀄라이저 오디오의 이퀄라이저 이퀄라이저(equalizer)는 공학적 기능으로 볼 때 평활기입니다. 즉 어떤 장치의 주파수 축 특성이 평탄하지 않을 때 이를 보정(compensation)해서 평활화하기 위해 사용합니다. 이를 오디오 장치에 개념을 확장해서 활용하고 있습니다. 오디오에서 이퀄라이저는 크게 두 가지 종류가 있는데 파라메트릭 이퀄라이저(Parameteric Equalizer, PEQ)와 그래픽 이퀄라이저(Graphic Equalizer, GEQ)입니다. 파라메트릭 이퀄라이저 일반적으로 오디오에서 이퀄라이저라고 하면, 파라메트릭 이퀄라이저를 의미합니다. 이퀄라이저는 입력되는 오디오 신호의 음색을 조정할 목적으로 사용됩니다. 예를 들어 사람마다 음성 특성이 다른데 각 특성이 작 두드러지도록 이퀄라이저로 조정.. 2021. 10. 29. 오디오 필터(4)-IIR 고역 통과 필터(Hi Pass Filter, HPF) 고역 통과 필터 다른 글에서 다룬 저역 통과 필터와 같은 방법으로 필터 계수를 구할 수 있습니다. https://medialink.tistory.com/74?category=958130 고역 통과 필터는 아래 그림과 같이 차단 주파수를 기준으로 차단 주파수보다 높은 주파수는 그대로 출력하고 다른 주파수 부분은 출력을 억제하는 특성을 가지고 있습니다. 이 부분에 대한 내용은 아래를 참조해 주세요 https://medialink.tistory.com/73?category=958130 https://medialink.tistory.com/77?category=958130 $$h=\frac{sin(2\pi f)}{2Q}$$ $$b_{0}=b_{1}=\frac{1+cos(2\pi f)}{2}$$ $$b_{1}= .. 2021. 10. 29. 오디오 필터(3)-IIR 저역 통과 필터 (Low Pass Filter, LPF) 저역 통과 필터 저역통과 필터는 필터의 차단 주파수를 기준으로 이보다 낮은 주파수는 입력과 출력이 같도록 통과 시키고, 나머지 부분은의 주파수는 통과를 차단하는 필터입니다. 이 필터는 오디오용 필터의 기본이 되는 중요한 필터입니다. 이중에서 음성확성에 사용되어 좋은 음질을 제공하기도 합니다. 이 글에서는 저역 통과 필터의 필터 계수를 구하는 방법을 알아 보겠습니다. 필터 계수는 신호 처리 과정에서 IIR 필터에 의해 실제 입력되는 신호를 처리과정에 사용되기 때문입니다. 다른 글에서 IIR을 다루도록 하겠습니다. Q를 이용한 저역 통과 필터의 설계 디지털 저역 통과 필터(low pass filter, LPF)를 구현하는 방법도 여러 가지가 있습니다. 여기서는 아날로그 필터의 Q값을 이용해서 필터를 구현하는.. 2021. 10. 29. 이전 1 2 다음