TWI843240B - 主動降噪濾波器的生成方法、存儲介質及耳機 - Google Patents
主動降噪濾波器的生成方法、存儲介質及耳機 Download PDFInfo
- Publication number
- TWI843240B TWI843240B TW111139684A TW111139684A TWI843240B TW I843240 B TWI843240 B TW I843240B TW 111139684 A TW111139684 A TW 111139684A TW 111139684 A TW111139684 A TW 111139684A TW I843240 B TWI843240 B TW I843240B
- Authority
- TW
- Taiwan
- Prior art keywords
- noise reduction
- signal
- filter
- parameters
- active noise
- Prior art date
Links
- 230000009467 reduction Effects 0.000 title claims abstract description 202
- 238000000034 method Methods 0.000 title claims abstract description 41
- 230000006870 function Effects 0.000 claims abstract description 72
- 238000012546 transfer Methods 0.000 claims abstract description 61
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 42
- 230000003044 adaptive effect Effects 0.000 claims abstract description 29
- 238000001914 filtration Methods 0.000 claims abstract description 28
- 238000004364 calculation method Methods 0.000 claims abstract description 4
- 230000004044 response Effects 0.000 claims description 50
- 230000015654 memory Effects 0.000 claims description 15
- 238000004891 communication Methods 0.000 claims description 10
- 230000000694 effects Effects 0.000 abstract description 10
- 238000010586 diagram Methods 0.000 description 10
- 238000013459 approach Methods 0.000 description 3
- 238000004590 computer program Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 230000005236 sound signal Effects 0.000 description 3
- 230000009471 action Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
Images
Abstract
本發明涉及主動降噪技術領域,公開一種主動降噪濾波器的生成方法、存儲介質及耳機。方法包括:獲取物理降噪信號及混合信號,根據混合信號與物理降噪信號,計算輸入信號,根據自適應濾波演算法,對輸入信號與物理降噪信號作自適應濾波,得到傳遞函數,根據傳遞函數,生成主動降噪濾波器。一方面,本實施例無需精確地測算物理聲學路徑或播放聲學路徑的傳遞函數,根據自適應演算法便可生成主動降噪濾波器,以便能夠主動降噪,從而提高降噪效率。另一方面,由於無需測算物理聲學路徑或播放聲學路徑的傳遞函數,便也沒有引入上述兩者的測算誤差,從而提高降噪精度和降噪效果。
Description
本發明涉及主動降噪技術領域,具體涉及一種主動降噪濾波器的生成方法、存儲介質及耳機。
主動降噪耳機能夠產生與雜訊信號的幅值相等和相位相反的主動降噪信號,使得主動降噪信號能夠抵消雜訊信號,從而得到降噪目的。
請參閱圖1,主動降噪耳機包括兩個重要的聲學路徑,即從雜訊源到人耳(戴著耳機時)的物理聲學路徑P(z),以及從耳機揚聲器到達人耳的播放聲學路徑G(z),此處為了實現主動降噪,主動降噪濾波器的降雜訊學路徑為H(z)。假設聲學路徑P(z)的傳遞函數為H1,降雜訊學路徑H(z) 的傳遞函數為H2,聲學路徑G(z)的傳遞函數為H3,由圖1可知:三個傳遞函數的關係是:
。
假設雜訊信號為
,那麼戴上耳機後聽到的經過物理降噪後的雜訊大小為:
,其中,
是
的時域運算式。
打開主動降噪耳機的降噪功能,人耳聽到的殘留雜訊信號e(n)如下式所示,該殘留雜訊在理想情況下應該為0,即:
其中,
是
的時域運算式。
為了得到降雜訊學路徑H(z)的傳遞函數H2,通常的做法是分別求取傳遞函數H1和傳遞函數H3,再根據傳遞函數H1和傳遞函數H3得到傳遞函數H2,此種作法需要先獲得傳遞函數H1和傳遞函數H3,但是要獲取準確的P(z)和G(z)相對來說是比較困難的,因為P(z)和方向性有關,與播放音源的音箱位置以及播放的聲源類型有關;G(z)也因人而異,每個人佩戴耳機的方式、耳廓形狀等外界因素各不相同,G(z)也會不同。因此,要獲得準確的P(z)和G(z)相對而言是很困難的,並且在測量P(z)和G(z)的傳輸路徑時會引入兩個誤差,最終在估計傳遞函數H2時,由於P(z)和G(z)誤差的疊加效應,會導致估計傳遞函數H2的誤差大大增加。
本發明實施例的一個目的旨在提供一種主動降噪濾波器的生成方法、存儲介質及耳機,用於解決上述現有技術存在的技術缺陷。
在第一方面,本發明實施例提供一種主動降噪濾波器的生成方法,包括:
獲取物理降噪信號,所述物理降噪信號為雜訊信號透過耳機後被後饋麥克風接收的信號;
獲取混合信號,所述混合信號為當播放相同所述雜訊信號且所述耳機在直通狀態下播放直通信號時,所述後饋麥克風所接收的信號;
根據所述混合信號與所述物理降噪信號,計算輸入信號;
根據自適應濾波演算法,對所述輸入信號與所述物理降噪信號作自適應濾波,得到傳遞函數;
根據所述傳遞函數,生成主動降噪濾波器。
可選地,所述根據所述傳遞函數,生成主動降噪濾波器包括:
根據所述傳遞函數,計算FIR濾波器的頻率回應參數及其頻響曲線;
根據所述FIR濾波器的頻率回應參數,生成n階IIR濾波器的參數;
根據所述n階IIR濾波器的參數,生成主動降噪濾波器,其中,n為正整數。
可選地,所述根據所述FIR濾波器的頻率回應參數,生成n階IIR濾波器的參數包括:
根據所述FIR濾波器的頻率回應參數與設定的濾波器階數n,結合離散傅立葉逆變換演算法,生成n階IIR濾波器的參數。
可選地,所述根據所述n階IIR濾波器的參數,生成主動降噪濾波器包括:
根據所述n階IIR濾波器的參數,結合離散傅立葉變換演算法,計算所述n階IIR濾波器的頻率回應參數;
根據所述n階IIR濾波器的頻率回應參數,生成n階IIR濾波器的頻響曲線;
根據所述n階IIR濾波器的頻響曲線,將所述n階IIR濾波器的參數降階至m階IIR濾波器的參數,其中,m為正整數且2<m<n;
根據所述m階IIR濾波器的參數,生成主動降噪濾波器。
可選地,所述根據所述m階IIR濾波器參數,生成主動降噪濾波器包括:
將所述m階IIR濾波器的參數對應的傳遞函數轉換成多個二階IIR濾波器級聯的二次分式模型;
根據多個所述二階IIR濾波器的參數,生成主動降噪濾波器。
可選地,所述直通信號為前饋麥克風採樣的雜訊信號。
可選地,所述根據所述混合信號與所述物理降噪信號,計算輸入信號包括:
將所述混合信號減去所述物理降噪信號,得到輸入信號。
可選地,所述自適應濾波演算法包括歸一化最小均方演算法,所述根據自適應濾波演算法,對所述輸入信號與物理降噪信號作自適應濾波,得到傳遞函數包括:
根據歸一化最小均方演算法,對所述輸入信號與物理降噪信號作自適應濾波,得到傳遞函數。
可選地,所述根據歸一化最小均方演算法,對所述輸入信號與物理降噪信號作自適應濾波,得到傳遞函數包括:
根據所述輸入信號,結合自適應濾波演算法,得到參考輸出;
將所述參考輸出減去所述物理降噪信號,得到誤差;
將所述誤差回饋至自適應濾波演算法模組,以使所述自適應濾波演算法模組調整主動降雜訊學路徑的傳遞函數,直至誤差接近或等於0,記錄最終的傳遞函數。
在第二方面,本發明實施例提供一種存儲介質,存儲有電腦可執行指令,所述電腦可執行指令用於使電子設備執行上述的主動降噪濾波器的生成方法。
在第三方面,本發明實施例提供一種耳機,包括:
至少一個處理器;以及,
與所述至少一個處理器通信連接的記憶體;其中,
所述記憶體存儲有可被所述至少一個處理器執行的指令,所述指令被所述至少一個處理器執行,以使所述至少一個處理器能夠執行上述的主動降噪濾波器的生成方法。
有益效果如下:在本發明實施例提供的主動降噪濾波器的生成方法中,獲取物理降噪信號,物理降噪信號為雜訊信號透過耳機後被後饋麥克風接收的信號,獲取混合信號,混合信號為當播放相同雜訊信號且耳機在直通狀態下播放直通信號時後饋麥克風所接收的信號,根據混合信號與物理降噪信號,計算輸入信號,根據自適應演算法,對輸入信號與物理降噪信號作自適應濾波,得到傳遞函數,根據傳遞函數,生成主動降噪濾波器,因此,一方面,本實施例無需精確地測算物理聲學路徑或播放聲學路徑的傳遞函數,根據自適應演算法便可生成主動降噪濾波器,以便能夠主動降噪,從而提高降噪效率。另一方面,由於無需測算物理聲學路徑或播放聲學路徑的傳遞函數,便也沒有引入上述兩者的測算誤差,從而提高降噪精度和降噪效果。
為了使本發明的目的、技術方案及優點更加清楚明白,以下結合附圖及實施例,對本發明進行進一步詳細說明。應當理解,此處所描述的具體實施例僅用以解釋本發明,並不用於限定本發明。基於本發明中的實施例,本領域普通技術人員在沒有作出創造性勞動前提下所獲得的所有其他實施例,都屬於本發明保護的範圍。
需要說明的是,如果不衝突,本發明實施例中的各個特徵可以相互結合,均在本發明的保護範圍之內。另外,雖然在裝置示意圖中進行了功能模組劃分,在流程圖中示出了邏輯順序,但是在某些情況下,可以以不同於裝置中的模組劃分,或流程圖中的順序執行所示出或描述的步驟。再者,本發明所採用的“第一”、“第二”、“第三”等字樣並不對數據和執行次序進行限定,僅是對功能和作用基本相同的相同項或相似項進行區分。
本發明實施例提供一種耳機。請參閱圖2,耳機100包括外殼10及安裝於外殼10上的前饋麥克風11、主動降噪濾波器12、揚聲器13、後饋麥克風14及控制器15。
前饋麥克風11安裝於耳機100的殼體外部,用於採樣雜訊源16播放的雜訊信號,雜訊信號透過外殼10被後饋麥克風14接收到。其中,雜訊源16可以為任意形式的雜訊源,例如音箱等。雜訊信號可以為任意合適形式的雜訊,諸如掃頻雜訊或粉紅雜訊等。
可以理解的是,雜訊信號透過外殼10後的信號被後饋麥克風14所採樣,此信號可被視為經過物理降噪後的信號,亦即為物理降噪信號。
主動降噪濾波器12受控制器15的控制,通過揚聲器13產生主動降噪信號,其中,控制器15根據前饋麥克風11採樣的雜訊信號,控制主動降噪濾波器12的工作狀態,例如,當前饋麥克風11採樣到雜訊信號時,控制器15可以啟動主動降噪濾波器12的主動降噪功能,當前饋麥克風11未採樣到雜訊信號時,控制器15可以關閉主動降噪濾波器12的主動降噪功能。
揚聲器13用於播放主動降噪信號,在理想情況下,主動降噪信號與雜訊信號的幅值相等和相位相反。當主動降噪濾波器12的主動降噪功能被控制器15關閉時,直通信號未經過主動降噪濾波器12的主動降噪處理,便可直接通過揚聲器13傳輸到外部環境。當主動降噪濾波器12的主動降噪功能被控制器15啟動時,直通信號經過主動降噪濾波器12的主動降噪處理,得到主動降噪信號,主動降噪信號通過揚聲器13傳輸到外部環境。
可以理解的是,揚聲器13可引入傳遞函數,如前所述,為了實現主動降噪目的,現有技術是需要測算出此傳遞函數的。
後饋麥克風14用於採樣雜訊信號透過外殼10後的物理降噪信號,和/或,採樣揚聲器13播放的主動降噪信號。
控制器15集成有自適應降噪演算法模組,可根據前饋麥克風11和/或後饋麥克風14採樣的信號,控制主動降噪濾波器12作自適應濾波和降噪,舉例而言,首先,雜訊源16的雜訊信號透過外殼10後被後饋麥克風14接收到,亦即後饋麥克風14接收到物理降噪信號,控制器15記錄物理降噪信號。
接著,控制器15關閉主動降噪濾波器12的主動降噪功能,亦即將主動降噪濾波器12的工作狀態設置為直通狀態,在直通狀態下,控制器15產生的直通信號不經過任何處理,直接通過主動降噪濾波器12傳輸至揚聲器13,再通過揚聲器13傳輸到外部。
再接著,雜訊源16播放同一雜訊信號,當前饋麥克風11將採樣的雜訊信號傳輸至控制器15時,控制器15同步控制主動降噪濾波器12在直通狀態下向揚聲器13發送直通信號,揚聲器13播放直通信號,於是,後饋麥克風14接收到由直通信號與物理降噪信號混合成的混合信號,亦即,混合信號為直通信號與物理降噪信號疊加得到,直通信號可以為前饋麥克風11採集的雜訊信號,也可以為隨意選擇的音頻信號。
請參閱圖3,A(z)為雜訊信號u(n)透過外殼10後傳輸到後饋麥克風的聲學路徑,其中,A(z)可為未知系統的聲學路徑。G(z)為主動降噪濾波器12處於直通狀態下,雜訊信號u(n)被前饋麥克風11採集後依次傳輸至主動降噪濾波器12和揚聲器13的聲學路徑,雜訊信號u(n)可作為直通信號。
在t1時刻,控制器15暫停主動降噪濾波器12的工作,雜訊信號u(n)經過聲學路徑A(z)後,可得到物理降噪信號d(n),所述物理降噪信號d(n)可被後饋麥克風14接收到。
在t2時刻,控制器15將主動降噪濾波器12的工作狀態設置為直通狀態,並且控制雜訊源16播放與t1時刻相同的雜訊信號u(n)。前饋麥克風11採樣雜訊信號u(n),並將採樣的雜訊信號u(n)經過聲學路徑G(z)後,可得到混合信號f(n),亦即採樣的雜訊信號u(n)作為直通信號輸入到主動降噪濾波器12,主動降噪濾波器12在直通狀態下,將直通信號傳輸至揚聲器13,以使揚聲器13輸出所述直通信號,後饋麥克風14同時採樣到由直通信號與物理降噪信號疊加形成的混合信號f(n),並將混合信號f(n)傳輸給控制器15。控制器15將混合信號f(n)減去物理降噪信號d(n),得到輸入信號x(n)。
控制器15對主動降噪濾波器12構建一條主動降雜訊學路徑B(z),其中,主動降雜訊學路徑B(z)的傳遞函數H0是由控制器15根據輸入信號x(n)與物理降噪信號d(n),結合自適應濾波演算法計算得到,並且,控制器15根據傳遞函數H0,計算數字濾波器參數,並將數字濾波器參數填入主動降噪濾波器12,例如,控制器15根據輸入信號x(n),結合自適應濾波演算法,得到參考輸出y(n),參考輸出y(n)減去物理降噪信號d(n),得到誤差e(n)。誤差e(n)回饋至控制器15的自適應濾波演算法模組,自適應濾波演算法模組再次調整主動降雜訊學路徑B(z)的傳遞函數H0,直至誤差e(n)接近或等於0,記錄此時的傳遞函數H0。後續,在應用過程中,如遇到雜訊信號,主動降噪濾波器12便可對雜訊信號進行有效地主動降噪,避免雜訊信號對用戶的干擾。
由上述描述可知,相對於現有技術,本實施例無需測算兩條聲學路徑的傳遞函數H1和傳遞函數H3,只需測算或調整一條聲學路徑即可,亦即只需測算主動降噪濾波器的傳遞函數H0,都可實現主動降噪目的。
作為本發明實施例另一方面,本發明實施例提供一種主動降噪濾波器的生成方法。請參閱圖4,主動降噪濾波器的生成方法S300包括:
S31、獲取物理降噪信號,物理降噪信號為雜訊信號透過耳機後被後饋麥克風接收的信號;
S32、獲取混合信號,混合信號為當播放相同雜訊信號且耳機在直通狀態下播放直通信號時,後饋麥克風所接收的信號;
作為示例但非限定的是,直通信號為耳機處於直通狀態下播放的音頻信號,其中,直通信號可以為前饋麥克風採樣的雜訊信號,由於前饋麥克風採樣的雜訊信號與物理降噪信號的傳輸環境或介質相同,有利於後續步驟能夠快速高效地計算出最優的傳遞函數。可以理解的是,直通信號還可以為任意選擇的音頻信號。
S33、根據混合信號與物理降噪信號,計算輸入信號;
在一些實施例中,耳機可以對混合信號與物理降噪信號作任意合適的處理,以得到輸入信號。在一些實施例中,耳機將混合信號減去物理降噪信號,得到輸入信號。
S34、根據自適應濾波演算法,對輸入信號與物理降噪信號作自適應濾波,得到傳遞函數;
在一些實施例中,自適應濾波演算法包括歸一化最小均方演算法,S34包括:根據歸一化最小均方演算法,對輸入信號與物理降噪信號作自適應濾波,得到傳遞函數,舉例而言:
歸一化最小均方演算法權重更新的向量形式為:
其中,
是第n次迭代時的權重向量,
是在
的基礎上更新出的權重參數向量,
是第n次迭代時的輸入向量,
是第n次迭代時主動降噪濾波器輸出的參考輸出與物理降噪信號之間的誤差。
是步長因數,
的取值影響著主動降噪濾波器的收斂速度和誤差,歸一化最小均方演算法中的步長為一個隨時間變化的量,且定義如下:
其中,
是n時刻估計出的信號的功率,
,
是修正的步長常量,
,
,是一個很小的常數,其目的是為了避免
時分母為零的情況,同時也避免輸入信號功率過小時產生較大的步長,此處取值為
。
請結合圖3,參考輸出y(n)表示為:
其中,
為主動降噪濾波器的權重係數,主動降噪濾波器輸出的參考輸出與物理降噪信號之間的誤差信號為:
在濾波器優化設計中,採用某種最小代價函數或者某個性能指標來衡量濾波器的好壞,而最常用的指標就是均方誤差,也把這種衡量濾波器好壞的方法叫做均方誤差準則。用公式表示如下:
其中,
是均方誤差,
表示物理降噪信號,
表示輸入信號
經過主動降噪濾波器處理後的信號,
表示輸入
時,濾波器的參考輸出和物理降噪信號之間的誤差。根據上式,耳機尋找最優濾波器權重係數
,使得
信號無限接近於
,誤差信號
無限接近於0,此時均方誤差達到最小值。
對權向量求導得到均方誤差的梯度
:
為了使基於歸一化最小均方誤差準則設計的濾波器性能達到最優,需要找到誤差性能曲面上的最小值,由此即可得到最優濾波器參數。在誤差性能曲面上沿著曲面的切線方向,也即負梯度方向進行尋找,沿著
負梯度方向調整濾波器權重係數
。設第n次迭代得到的濾波器抽頭權值向量為
,並設該次迭代得到的均方誤差是
,那麼第n + 1次迭代得到的濾波器係數可由下式求出:
其中
是該次迭代時的梯度向量,
就是該次迭代的方向向量,
是第n次迭代時所用的步長,又將其稱作收斂因數。精確計算梯度
是十分困難的,一種粗略的但十分有效的計算
方法是直接取誤差二次方
作為均方誤差
的估計值,即:
其中,
為:
則均方誤差估計值為:
於是得到主動降噪濾波器抽頭參數的更新方式為:
當均方誤差最小時,得到最佳濾波器權重係數向量
,則即輸出信號向量
,傳遞函數
,因為採用的是FIR濾波器,得到濾波器參數b=W,a=1,因此,可得到由濾波器參數a和b形成的傳遞函數。
S35、根據傳遞函數,生成主動降噪濾波器。
綜上所述,一方面,本實施例無需精確地測算物理聲學路徑或播放聲學路徑的傳遞函數,根據自適應演算法便可生成主動降噪濾波器,以便能夠主動降噪,從而提高降噪效率。另一方面,由於無需測算物理聲學路徑或播放聲學路徑的傳遞函數,便也沒有引入上述兩者的測算誤差,從而提高降噪精度和降噪效果。
在一些實施例中,請參閱圖5,S35包括:
S351、根據傳遞函數,計算FIR濾波器的頻率回應參數及其頻響曲線;
S352、根據FIR濾波器的頻率回應參數,生成n階IIR濾波器的參數;
S353、根據n階IIR濾波器的參數,生成主動降噪濾波器,其中,n為正整數。
在步驟S351中,在一些實施例中,耳機可以根據傳遞函數,結合離散傅立葉變換演算法,計算FIR濾波器的頻率回應參數,根據FIR濾波器的頻率回應參數,繪製FIR濾波器的頻響曲線。
在步驟S352中,在一些實施例中,耳機根據FIR濾波器的頻率回應參數與設定的濾波器階數n,結合離散傅立葉逆變換演算法,生成n階IIR濾波器的參數。其中,頻率回應參數包括h1參數與w1參數,h1包含了離散系統對應區間內N個頻率等分點的頻率回應,N為正整數,w1為N個頻率等分點的值,然後根據h1和w1,設定想要的IIR濾波器階數n,進行離散傅立葉逆變換,得到新的IIR濾波器參數b_new,a_new,亦即,IIR濾波器的參數b_new和a_new可形成n階IIR濾波器。
在步驟S353中,在一些實施例中,耳機根據n階IIR濾波器的參數,結合離散傅立葉變換演算法,計算n階IIR濾波器的頻率回應參數,根據n階IIR濾波器的頻率回應參數,生成n階IIR濾波器的頻響曲線,根據n階IIR濾波器的頻響曲線,將n階IIR濾波器的參數降階至m階IIR濾波器的參數,其中,m為正整數且2<m<n,根據m階IIR濾波器的參數,生成主動降噪濾波器。
舉例而言,在步驟S352中可得到n階IIR濾波器的參數b_new和a_new,進行離散傅立葉變換求取離散系統的頻率回應參數,頻率回應參數包括h21和w21,並繪製頻響曲線。設計小於n的IIR濾波器階數i,根據頻率回應參數與設定的濾波器階數i,離散傅立葉逆變換得到i階IIR濾波器參數,得到的i階IIR濾波器的參數為b_new21和a_new21。根據i階IIR濾波器的參數b_new21和a_new21,進行離散傅立葉變換求其頻率回應參數h31和w31,繪製頻響曲線。
將i階IIR濾波器頻響曲線與n階IIR濾波器頻響曲線作相似比對,若相似程度小於或等於預設相似閾值,則提高i值以得到新i值,然後再使用上述方法,得到新i階IIR濾波器頻響曲線。
將新i階IIR濾波器頻響曲線與n階IIR濾波器頻響曲線作相似比對,以此類推,直至將IIR濾波器的階數降低至m。若相似程度大於預設相似閾值,則降階成功,判斷i是否等於m,若等於m,則記錄i=m。若不等於m,賦值n=i,將i設置為用戶期望的階數,再使用上述方法,可得到新i階IIR濾波器頻響曲線。
舉例而言,耳機根據FIR濾波器頻響曲線,生成512階IIR濾波器頻響曲線。耳機先設置i=64,根據上述方法,得到64階IIR濾波器頻響曲線。
耳機比對512階IIR濾波器頻響曲線與64階IIR濾波器頻響曲線的相似程度,若相似程度小於或等於預設相似閾值,則降階失敗,失敗的原因是64階IIR濾波器和512階IIR濾波器階數差距太大,無法用64階IIR濾波器來表示512階IIR濾波器,因此增加擬合的濾波器階數。
於是,將64階增加到128階,使用128階IIR濾波器頻響曲線擬合512階IIR濾波器頻響曲線,若相似程度大於預設相似閾值,則降階成功。
由於對於實際應用來說,128階IIR濾波器消耗較大,實現起來依然較為困難,因此,還需要重複上面的步驟,再次把剛剛得到的128階IIR濾波器降階轉換為64階,64階再轉換為16階等,依次類推。
在一些實施例中,m=16,由於16階IIR濾波器在工程上實現依然存在難度,因此,還需要將16階IIR濾波器轉換為級聯形式的多個二階IIR濾波器,在一些實施例中,耳機將m階IIR濾波器的參數對應的傳遞函數轉換成多個二階IIR濾波器級聯的二次分式模型,根據多個二階IIR濾波器的參數,生成主動降噪濾波器,舉例而言,耳機根據函數tf2sos,將m階IIR濾波器轉換成級聯形式的多個二階IIR濾波器,根據多個二階IIR濾波器的參數,生成主動降噪濾波器。
為了表現採用本實施例提供的主動降噪濾波器的降噪效果,本文結合圖6與圖7分別作出說明,其中,在圖6中,第一曲線51表示經過物理降噪人工耳接收到的外界雜訊,第二曲線52表示在主動降噪濾波器的作用下後饋麥克風接收到的殘留雜訊。在圖7中,第三曲線61表示經過物理降噪人工耳接收到的外界雜訊,第四曲線62表示在主動降噪濾波器的作用下後饋麥克風接收到的殘留雜訊。
主動降噪耳機的性能指標一般有降噪頻寬和降噪深度,降噪頻寬,指的是耳機可以處理雜訊頻率的範圍。由於不同種類的聲音頻率不一樣,因此,降噪頻寬越大,覆蓋的頻率越多,耳機能降噪的聲音種類也越多。降噪深度,指的是針對某個頻率的雜訊,可以降低音量的多少,數值越大表示降噪效果越好。一般用降噪深度最大值,作為整個耳機的降噪深度。降噪頻寬表示可以處理的聲音種類,實際處理後降噪的效果好壞,是由該頻率下的降噪深度決定的。由圖6與圖7可知,降噪頻寬分別在50hz-5khz、50hz-10khz之間,降噪深度大致在20-35dB範圍內。該主動降噪演算法的降噪頻寬和降噪深度都挺可觀,具有一定的實用價值。
需要說明的是,在上述各個實施方式中,上述各步驟之間並不必然存在一定的先後順序,本領域普通技術人員,根據本發明實施方式的描述可以理解,不同實施方式中,上述各步驟可以有不同的執行順序,亦即,可以並行執行,亦可以交換執行等等。
請參閱圖8,圖8為本發明實施例提供的一種耳機的電路結構示意圖。如圖8所示,耳機700包括一個或多個處理器71以及記憶體72。其中,圖8中以一個處理器71為例。
處理器71和記憶體72可以通過匯流排或者其他方式連接,圖8中以通過匯流排連接為例。
記憶體72作為一種非易失性電腦可讀存儲介質,可用於存儲非易失性軟體程式、非易失性電腦可執行程式以及模組,如本發明實施例中的主動降噪濾波器的生成方法對應的程式指令/模組。處理器71通過運行存儲在記憶體72中的非易失性軟體程式、指令以及模組,從而實現上述方法實施例提供的主動降噪濾波器的生成方法的功能。
記憶體72可以包括高速隨機存取記憶體,還可以包括非易失性記憶體,例如至少一個磁片記憶體件、閃存器件、或其他非易失性固態記憶體件。在一些實施例中,記憶體72可選包括相對於處理器71遠程設置的記憶體,這些遠程記憶體可以通過網路連接至處理器71。上述網路的實例包括但不限於互聯網、企業內部網、區域網、移動通信網及其組合。
所述程式指令/模組存儲在所述記憶體72中,當被所述一個或者多個處理器71執行時,執行上述任意方法實施例中的主動降噪濾波器的生成方法。
本發明實施例還提供了一種存儲介質,所述存儲介質存儲有電腦可執行指令,該電腦可執行指令被一個或多個處理器執行,例如圖8中的一個處理器71,可使得上述一個或多個處理器可執行上述任意方法實施例中的主動降噪濾波器的生成方法。
本發明實施例還提供了一種電腦程式產品,所述電腦程式產品包括存儲在非易失性電腦可讀存儲介質上的電腦程式,所述電腦程式包括程式指令,當所述程式指令被耳機執行時,使所述耳機執行任一項所述的主動降噪濾波器的生成方法。
以上所描述的裝置或設備實施例僅僅是示意性的,其中所述作為分離部件說明的單元模組可以是或者也可以不是物理上分開的,作為模組單元顯示的部件可以是或者也可以不是物理單元,即可以位於一個地方,或者也可以分佈到多個網路模組單元上。可以根據實際的需要選擇其中的部分或者全部模組來實現本實施例方案的目的。
通過以上的實施方式的描述,本領域的技術人員可以清楚地瞭解到各實施方式可借助軟體加通用硬體平臺的方式來實現,當然也可以通過硬體。基於這樣的理解,上述技術方案本質上或者說對相關技術做出貢獻的部分可以以軟體產品的形式體現出來,該電腦軟體產品可以存儲在電腦可讀存儲介質中,如ROM/RAM、磁碟、光碟等,包括若干指令用以使得一臺電腦設備(可以是個人電腦,伺服器,或者網路設備等)執行各個實施例或者實施例的某些部分所述的方法。
最後應說明的是:以上實施例僅用以說明本發明的技術方案,而非對其限制;在本發明的思路下,以上實施例或者不同實施例中的技術特徵之間也可以進行組合,步驟可以以任意順序實現,並存在如上所述的本發明的不同方面的許多其他變化,為了簡明,它們沒有在細節中提供;儘管參照前述實施例對本發明進行了詳細的說明,本領域的普通技術人員應當理解:其依然可以對前述各實施例所記載的技術方案進行修改,或者對其中部分技術特徵進行等同替換;而這些修改或者替換,並不使相應技術方案的本質脫離本發明各實施例技術方案的範圍。
10:外殼
11:前饋麥克風
12:主動降噪濾波器
13:揚聲器
14:後饋麥克風
15:控制器
16:雜訊源
700:耳機
71:處理器
72:記憶體
S300:主動降噪濾波器的生成方法
S31、S32、S33、S34、S35、S351、S352、S353:步驟
圖1為現有技術提供的聲學路徑的架構圖;
圖2為本發明實施例提供的一種耳機的結構示意圖;
圖3為圖2所示耳機的聲學路徑的架構圖;
圖4為本發明實施例提供的一種主動降噪濾波器的生成方法的流程示意圖;
圖5為圖4所示的S35的流程示意圖;
圖6為採用圖4所示的主動降噪濾波器進行降噪的第一種效果示意圖;
圖7為採用圖4所示的主動降噪濾波器進行降噪的第二種效果示意圖;
圖8為本發明實施例提供的一種耳機的電路結構示意圖。
S300:主動降噪濾波器的生成方法
S31、S32、S33、S34、S35:步驟
Claims (9)
- 一種主動降噪濾波器的生成方法,其中,包括:獲取物理降噪信號,所述物理降噪信號為雜訊信號透過耳機後被後饋麥克風接收的信號;獲取混合信號,所述混合信號為當播放相同所述雜訊信號且所述耳機在直通狀態下播放直通信號時,所述後饋麥克風所接收的信號;根據所述混合信號與所述物理降噪信號,計算輸入信號;根據自適應濾波演算法,對所述輸入信號與所述物理降噪信號作自適應濾波,得到傳遞函數;根據所述傳遞函數,生成主動降噪濾波器,包括:根據所述傳遞函數,計算FIR濾波器的頻率回應參數及其頻響曲線;根據所述FIR濾波器的頻率回應參數,生成n階IIR濾波器的參數;根據所述n階IIR濾波器的參數,生成主動降噪濾波器,其中,n為正整數。
- 如請求項第1項所述的方法,其中,所述根據所述FIR濾波器的頻率回應參數,生成n階IIR濾波器的參數包括:根據所述FIR濾波器的頻率回應參數與設定的濾波器階數n,結合離散傅立葉逆變換演算法,生成n階IIR濾波器的參數。
- 如請求項第1項所述的方法,其中,所述根據所述n階IIR濾波器的參數,生成主動降噪濾波器包括:根據所述n階IIR濾波器的參數,結合離散傅立葉變換演算法,計算所述n階IIR濾波器的頻率回應參數;根據所述n階IIR濾波器的頻率回應參數,生成n階IIR濾波器的頻響曲線;根據所述n階IIR濾波器的頻響曲線,將所述n階IIR濾波器的參數降階至m階IIR濾波器的參數,其中,m為正整數且2<m<n;根據所述m階IIR濾波器的參數,生成主動降噪濾波器。
- 如請求項第3項所述的方法,其中,所述根據所述m階IIR濾波器參數,生成主動降噪濾波器包括:將所述m階IIR濾波器的參數對應的傳遞函數轉換成多個二階IIR濾波器級聯的二次分式模型;根據多個所述二階IIR濾波器的參數,生成主動降噪濾波器。6、如請求項第1項所述的方法,其中,所述直通信號為前饋麥克風採樣的雜訊信號。
- 如請求項第1至4任一項所述的方法,其中,所述根據所述混合信號與所述物理降噪信號,計算輸入信號包括:將所述混合信號減去所述物理降噪信號,得到輸入信號。
- 如請求項第1至4任一項所述的方法,其中,所述自適應濾波演算法包括歸一化最小均方演算法,所述根據自適應濾波演算法,對所述輸入信號與物理降噪信號作自適應濾波,得到傳遞函數包括: 根據歸一化最小均方演算法,對所述輸入信號與物理降噪信號作自適應濾波,得到傳遞函數。
- 如請求項第6項所述的方法,其中,所述根據歸一化最小均方演算法,對所述輸入信號與物理降噪信號作自適應濾波,得到傳遞函數包括:根據所述輸入信號,結合自適應濾波演算法,得到參考輸出;將所述參考輸出減去所述物理降噪信號,得到誤差;將所述誤差回饋至自適應濾波演算法模組,以使所述自適應濾波演算法模組調整主動降雜訊學路徑的傳遞函數,直至誤差接近或等於0,記錄最終的傳遞函數。
- 一種存儲介質,存儲有電腦可執行指令,所述電腦可執行指令用於使電子設備執行如請求項第1至7任一項所述的主動降噪濾波器的生成方法。
- 一種耳機,包括外殼、前饋麥克風、主動降噪濾波器、揚聲器、後饋麥克風及控制器,前饋麥克風、主動降噪濾波器、揚聲器、後饋麥克風及控制器都安裝於所述外殼上,其中,所述控制器包括:至少一個處理器;以及,與所述至少一個處理器通信連接的記憶體;其中,所述記憶體存儲有可被所述至少一個處理器執行的指令,所述指令被所述至少一個處理器執行,以使所述至少一個處理器能夠執行如請求項第1至7任一項所述的主動降噪濾波器的生成方法。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111225489.9A CN114339513A (zh) | 2021-10-21 | 2021-10-21 | 主动降噪滤波器的生成方法、存储介质及耳机 |
CN202111225489.9 | 2021-10-21 |
Publications (2)
Publication Number | Publication Date |
---|---|
TW202341747A TW202341747A (zh) | 2023-10-16 |
TWI843240B true TWI843240B (zh) | 2024-05-21 |
Family
ID=
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113490098A (zh) | 2021-07-07 | 2021-10-08 | 东莞市逸音电子科技有限公司 | 一种anc耳机主动降噪滤波器主动优化算法 |
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113490098A (zh) | 2021-07-07 | 2021-10-08 | 东莞市逸音电子科技有限公司 | 一种anc耳机主动降噪滤波器主动优化算法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8693699B2 (en) | Method for adaptive control and equalization of electroacoustic channels | |
TWI611704B (zh) | 自調式主動噪聲消除方法、系統及耳機裝置 | |
US20160300562A1 (en) | Adaptive feedback control for earbuds, headphones, and handsets | |
JP4469898B2 (ja) | 外耳道共鳴補正装置 | |
CN113574593B (zh) | 调谐方法、制造方法、计算机可读存储介质和调谐系统 | |
WO2020248164A1 (zh) | 主动降噪方法、装置、芯片、主动控制系统和存储介质 | |
CN108540895A (zh) | 智能均衡器设计方法及具有智能均衡器的降噪耳机 | |
WO2020051769A1 (zh) | 主动降噪方法及耳机 | |
CN112562624B (zh) | 主动降噪滤波器设计方法、降噪方法、系统及电子设备 | |
CN114787911A (zh) | 耳戴式播放设备的噪声消除系统和信号处理方法 | |
Fabry et al. | Acoustic equalization for headphones using a fixed feed-forward filter | |
CN113409755A (zh) | 主动降噪方法、装置及主动降噪耳机 | |
US20230131827A1 (en) | Method for generating active noise reduction filter, storage medium and earphone | |
TWI843240B (zh) | 主動降噪濾波器的生成方法、存儲介質及耳機 | |
CN107666637B (zh) | 自调式主动噪声消除方法、系统及耳机装置 | |
TW201626814A (zh) | 揚聲器之頻率響應補償系統 | |
CN111866666B (zh) | 数字降噪滤波器生成方法、相关设备及可读存储介质 | |
CN113488016B (zh) | 系数确定方法及装置 | |
US11206004B1 (en) | Automatic equalization for consistent headphone playback | |
TWI746001B (zh) | 頭戴式裝置及其立體音效控制方法 | |
CN113345401A (zh) | 可穿戴设备的主动降噪系统的校准方法及装置、存储介质、终端 | |
Hilgemann et al. | Design of IIR filters for active noise control by constrained optimization | |
TWI837867B (zh) | 聲音補償方法及頭戴式裝置 | |
US11790882B2 (en) | Active noise cancellation filter adaptation with ear cavity frequency response compensation | |
JP7457149B2 (ja) | 低周波数フィードバックループ変調を有するウェアラブルアクティブノイズ低減(anr)デバイス |