麦克风阵列声源定位实现
) p3 X8 W. z3 V; x( l麥克風陣列音源定位系統是利用麥克風陣列接收音訊,然後經由適當的演算法估算出音源入射麥克風陣列的方向角,即判斷出音源的位置方向。本專題著重於此演算法的程式撰寫,並實際測試二維及三維空間中單一音源入射的方向角,期望於演算速度及準確度上能有最好的效果。此系統完成後可搭配攝影機應用於視訊會議,隨時定位出發言者;也可應用於監視系統、玩具等。 麥克風陣列音源定位系統主要是利用同一音源到麥克風陣列中每支麥克風的距離不全相同,因此同一音源的訊號傳遞到每支麥克風會有時間差TDOA ( Time Difference of Arrival ),利用求得的TDOA,代入推論出的方向角公式即可得到音源入射的方向角。而本專題著重於估算TDOA及方向角公式的演算法上,我們使用兩種演算法實作測試,其中一種演算法只適用於二維空間,而另一種適用於二維或三維空間的音源定位。
: ^/ d1 ^+ t; F有許多相關的論文是以頻域的方法分析,但若要將此系統實作在微處理器上,運算較簡單時域分析方法會有較好的效用,且對初學訊號分析的我們來說也較容易實現,所以我們主要以時域的方法來建立演算法。
0 ?5 l$ z3 _! \9 b3 Q% z) `+ v此系統主要的實作流程為:「語料取樣」、「訊號前處理」、「估算TDOA」、「計算方向角」。下圖1-1為麥克風陣列音源定位系統流程圖,詳細說明如下: ) o! c6 |8 z" v& F9 H/ h
由兩支以上麥克風組成麥克風陣列接收音訊。 + d$ B- P; x: \5 C7 X
麥克風接收的訊號經過一放大器後,再經由四輸入USB動態訊號擷取模組,依設定的取樣頻率取樣訊號,輸入電腦分析。
( A- L ]1 n4 T+ J6 D) ]. ]2 Z# G! R將每支麥克風接收的訊號音量標準化,並切割成多個音框(frame)。
) x! ]( j& x1 I* U& k1 m$ X設定音量門檻值,去除音量小的部份。(算出一段訊號中所有音框的音量,取其中音框音量最小值的k倍為門檻值,k值由經驗得到。) " Q" `) I1 Q5 |. }6 t- j
以Cross Correlation的方法估算TDOA。 ' |9 H4 T5 R I8 U" }& q6 S
將求得的TDOA代入方向角公式,即可得到音源入射的方向角。 ) }& w. ~1 a# J) e, ]# H& N' x
file:///C:/Users/ADMINI~1/AppData/Local/Temp/msohtmlclip1/01/clip_image002.gif TDOA的計算方法: V- d+ i y, i6 S: o3 P% j
TDOA (Time Difference of Arrival) 是指同一音源的訊號傳遞到兩支麥克風的時間差由麥克風所取樣的訊號來估算。而時域上常見的計算方法有四種:AMDF(Average Magnitude Difference Function)、ratio AMDF、最小平方法、CrossCorrelation。我们採用最常見且穩定性高的Cross Correlation。 1 m O: S7 { J* J4 j
Cross Correlation % a3 i( Z! c1 H b5 ]
x1及x2為兩支麥克風所取樣到的兩段訊號,經由下式可算出所對應λ值的cross correlation值。
. L! N' L% U7 }: A* f# s0 q
1 g+ X: E5 o, yfile:///C:/Users/ADMINI~1/AppData/Local/Temp/msohtmlclip1/01/clip_image004.gif & B" r' \+ D* i6 t1 g
# a1 L* m% |* y8 y3 w+ I( O; @
file:///C:/Users/ADMINI~1/AppData/Local/Temp/msohtmlclip1/01/clip_image006.gif
; K( J; I9 J* e% a最大合理平移範圍當音源與兩支麥克風的中心成90度時(音源與兩支麥克風距離相等),音源傳遞到兩支麥克風的時間差為0;當音源與兩支麥克風三者成一直線時,音源傳遞到兩支麥克風的時間差最大。 file:///C:/Users/ADMINI~1/AppData/Local/Temp/msohtmlclip1/01/clip_image008.gif
* Q! v6 `- g) R: u8 R3 o———————————————— 版权声明:本文为CSDN博主「narytowzy」的原创文章,遵循CC 4.0 BY-SA版权协议,转载请附上原文出处链接及本声明。 原文链接:https://blog.csdn.net/narytowzy/article/details/46771215 ) n2 y( C1 P. `0 F, Y8 `! e$ `
/ X6 J. S0 I, S( O+ A! N1 b. m
: j! D+ }; S7 D2 T |