開門見山,動態範圍的定義就是訊號的最大幅值和最小幅值比值的對數(單位dB),
動態範圍會受到系統中各個環節的影響。例如同樣是這段音樂,在一個40dB背景噪聲的環境中播放,那麼由於掩蔽效應等因素的影響,最終實際聽到的音樂動態範圍實際上是50dB。
動態範圍會受到系統中各個環節的影響。例如同樣是這段音樂,在一個40dB背景噪聲的環境中播放,那麼由於掩蔽效應等因素的影響,最終實際聽到的音樂動態範圍實際上是50dB。如果對響度和動態範圍感興趣的可以移步看鬼斧神工的文章《詳解音訊中的響度與動態範圍。(推力、耳放、響度戰爭、高解析度等)》
動態範圍控制(Dynamic Range Control,DRC)將輸入音訊訊號的動態範圍對映到指定的動態範圍,可以使聲音聽起來更柔和或更大聲,即一種訊號幅度調節方式。通常對映後的動態範圍小於對映前的動態範圍,因此稱之為動態範圍壓縮。音訊訊號可以進行整體的動態範圍控制,也可以劃分為若干子帶分別進行動態範圍控制。DRC廣泛應用於音訊訊號處理領域,例如助聽器中最常見的寬動態範圍壓縮方法(Wide Dynamic Range Compression,WDRC)、音訊訊號處理中最常用的自動增益控制(Automatic Gain Control,AGC)方法等。DRC側重於從數位域上對增益進行有選擇、有目的地調整,以保證最大訊號有合適的headroom,中間訊號得到合適的amplify,小訊號或者底噪被cutoff掉。
作用
動態範圍控制的型別
DRC的增益調整基本上是三段:對底噪,中等幅度訊號,幅度較大的訊號的三段不同的處理,比如:
整體而言演演算法較為簡單,其中最關鍵的地方在於曲線的設計。DRC的曲線一般具備下圖這樣的結構,該圖表示的是來濾除小訊號的曲線圖,如果想要動態拉伸的功能,那麼整體曲線需要高於$y=x$這條直線。
靜態曲線(圖片來自《Digital Audio Signal Processing-Udo Zölzer》第七章)
名詞解釋
- NT:噪聲門限(Noise Gate),低於NT的輸入訊號將會被濾除
- LT:限制器閾值(Limiter threshold),高於CT且低於LT的訊號也會被壓縮,而高於LT的訊號將會被直接限幅
- CT:壓縮器閾值(Compressor threshold),高於ET且低於CT的訊號將會保持不變
- ET:擴充套件器閾值(Expander threshold),高於NT且低於ET的訊號將會被抑制
MALTAB的Audio Toolbox中實現了 動態範圍壓縮器、動態範圍限幅器、動態範圍擴充套件器和噪聲門。本教學還提供了動態範圍限制器各個階段的動態範圍限制範例。
DRC的原理其實很簡單,就是通過設計一條曲線,將輸入語音幅度$x(n)$通過曲線進行對映得到另一語音幅度值$y(n)$。然後計算兩者之間的差值得到增益值$g_c(n)$,然後根據設定的attack time和release time進行增益平滑以及計算make-up增益,最後再應用到語音得到處理後的語音。基本流程如下
通用的動態範圍控制系統 (來自Matlab官網)
從線性值轉換為dB值
$$x_{dB}(n) = 20log_{10}(x(n))$$
通過將$x_{dB}(n) $訊號傳遞到靜態特徵方程(DRC靜態曲線),假設設計的曲線函數關係為$f(x)$
$$x_{sc}(x_{db})=f(x_{db}(n))$$
計算差值:$g_c(n)=x_{sc}(n)-x_{dB}(n)$
不同型別的DRC都具有不同的靜態特徵和不同的可調屬性:
在下面這些靜態特性圖中,expander、limiter和compressor 各有一個 knee width 的拐點寬度。
增益平滑減少了增益的急劇變化,避免偽像和不自然的聲音。expander和noiseGate具有相同的平滑方程(都是擴充套件器),limiter和compressor具有相同的平滑方程(都是壓縮器)。
增益平滑的型別由啟動時間(attack time)、釋放時間(release time)和保持時間(hold time)組合指定。啟動時間和釋放時間對應於增益訊號從其最終值的10%上升到90%所需的時間。保持時間是應用增益之前的延遲時間。
expander 和 noiseGate 的平滑方程
$$g_s[n]=\left\{\begin{array}{cl}
\alpha_A g_s[n-1]+\left(1-\alpha_A\right) g_c[n] & \text { if }\left(C_A>k\right) \&\left(g_c[n] \leq g_s[n-1]\right) \\
g_S[n-1] & \text { if } C_A \leq k \\
\alpha_{R}g_s[n-1]+\left(1-\alpha_R\right) g_c[n] & \text { if }\left(C_R>k\right) \&\left(g_c[n]>g_s[n-1]\right) \\
g_S[n-1] & \text { if } C_R \leq k
\end{array}\right.$$
$\alpha_A$和$\alpha_R$由取樣率和指定的啟動和釋放時間確認:$\alpha_A=exp(\frac{-log(9)}{Fs*T_A})$,$\alpha_R=exp(\frac{-log(9)}{Fs*T_R})$
$k$是以樣本為單位的指定保持時間。$C_A$和$C_R$分別是啟動和釋放hold counters
compressor 和 limiter 的平滑方程
$$g_s(n)=\left\{\begin{matrix}
a_Ag_s(n-1)+(1-\alpha_A)g_c(n),\ \ g_c\le g_s(n-1) \\
a_Rg_s(n-1)+(1-\alpha_R)g_c(n),\ \ g_c > g_s(n-1)
\end{matrix}\right.$$
$\alpha_A$和$\alpha_R$由取樣率和指定的啟動和釋放時間確認:$\alpha_A=exp(\frac{-log(9)}{Fs*T_A})$,$\alpha_R=exp(\frac{-log(9)}{Fs*T_R})$
增益平滑範例
輸入訊號兩步動態範圍壓縮的一個簡單案例。在此範例中,壓縮器的閾值為 –10 dB,壓縮比為 5,拐點較硬。
下圖為了顯示了增益平滑的幾種變化。
補償增益適用於compressors 和 limiters,其中訊號的較高 dB 部分被衰減或阻塞。dB 降低會顯著降低總訊號功率。在這些情況下,增益平滑後應用補償增益以增加訊號功率。在MATLAB的 Audio Toolbox 中,可以指定一個設定的補償量或將補償模式指定為'auto'。
'auto'補償增益確保 0 dB 輸入產生 0 dB 輸出。例如,假設具有軟拐點的compressors 的靜態特性:
$$x_{s c}\left(x_{d B}\right)=\left\{\begin{array}{cc}
x_{d B} & x_{d B}<\left(T-\frac{W}{2}\right) \\
x_{d B}+\frac{\left(\frac{1}{R}-1\right)\left(x_{d B}-T+\frac{W}{2}\right)^2}{2 W} & \left(T-\frac{W}{2}\right) \leq x_{d B} \leq\left(T+\frac{W}{2}\right) \\
T+\frac{\left(x_{d B}-T\right)}{R} & x_{d B}>\left(T+\frac{W}{2}\right)
\end{array}\right.$$
$T$是threshold,$W$是knee width,$R$是compression ratio。計算出的自動補償增益是在 0 dB 時評估的靜態特性方程的負值:
$$\text { MAKE-UP GAIN }=-x_{s c}(0)=\left\{\begin{array}{cl}
0 & \frac{W}{2}<T \\
-\frac{\left(\frac{1}{R}-1\right)\left(T-\frac{W}{2}\right)^2}{2 W} & -\frac{W}{2} \leq T \leq \frac{W}{2} \\
-T+\frac{T}{R} & -\frac{W}{2}>T
\end{array}\right.$$
將增益從dB域變換到線性域:$g_{lin}(n)=10^{\frac{g_m(n)}{20}}$
$$y(n)=x(n)*g_{lin}(n)$$
範例:動態範圍限制器
本例中描述的音訊訊號是 0.5 秒間隔的鼓音。limiter屬性是:
此範例提供了動態範圍限制器系統各個階段的視覺化演練。
線性到 dB 轉換
增益計算
靜態特性將 dB 訊號限制在 –15 dB。為了確定導致此限制的 dB 增益,增益計算從靜態特性處理的 dB 訊號中減去原始 dB 訊號。
增益平滑
當所應用的增益突然增加時,相對較短的啟動時間規格會導致陡峭的曲線。相對較長的釋放時間會導致應用增益逐漸減小。
補償增益
假設限制器具有 1 dB 補償增益值。補償增益被新增到平滑增益訊號。
dB 到線性轉換
以 dB 為單位的增益逐個元素轉換為線性標度。
應用計算增益
原始訊號乘以線性增益
【書籍】Digital Audio Signal Processing 2nd Edition
【音訊探險記】DRC--動態範圍控制
【MATLAB】Dynamic Range Control
【MATLAB】Multiband Dynamic Range Compression
【知乎】詳解音訊中的響度與動態範圍。(推力、耳放、響度戰爭、高解析度等)
【論文】Zolzer, Udo. "Dynamic Range Control." Digital Audio Signal Processing. 2nd ed. Chichester, UK: Wiley, 2008.
【論文】Giannoulis, Dimitrios, Michael Massberg, and Joshua D. Reiss. "Digital Dynamic Range Compressor Design –– A Tutorial And Analysis." Journal of Audio Engineering Society. Vol. 60, Issue 6, 2012, pp. 399–408.