盲源分離(BSS, Blind Source Separation)

數學描述:

假設N個統計獨立的未知信號S(t)

S(t)=[s_{1}(t),s_{2}(t),s_{3}(t),...,s_{N}(t)]^{T}

經過未知信道A的傳輸後由M個傳感器檢測獲得M個觀測信號

X(t)=[x_{1}(t),x_{2}(t),x_{3}(t),...,x_{M}(t)]^{T}

整個傳輸過程的數學模型爲:

X(t)=AS(t)+n(t)

X(t)爲M維觀測矢量,S(t)爲N維未知源信號矢量,n(t)爲M維加性信道噪聲,A爲M\times N維傳遞函數矩陣。

盲源分離問題就是求一個分離矩陣W,使得觀測信號X(t)通過該矩陣,儘量的完全分離出源信號S(t)的各個組成,設Y(t)爲源信號估計矢量,則分離系統的數學描述爲:Y(t)=WX(t)

實際中,傳感器測得的信號是源信號及其延時信號的混迭,通常稱卷積混迭。

X(t)=AS(t)=\sum_{-\infty}^{\infty}A(k)S(t-k)

針對這種情況,盲反捲積方法僅僅是通過觀測信號X(t)估計信道衝激響應A(k)進而恢復源信號,這樣得到如下分離系統,也叫盲均衡系統:Y(t)=\sum_{-\infty}^{\infty}W(k)S(t-k)

對無延時的情況,Y和S相差一個信號幅度的放大倍數,對於有延時的情況,Y和S相差一個濾波器函數。

信號被“分離”指的是每個分離出來的信號y_{i}與某個源信號s_{j}有關,i和j可以不相同。信號被恢復是恢復後的信號y_{i}與某個源信號s_{j}僅僅只相差一個固定的幅度值。

進行盲分離通常作以下假設:

1、信號源的個數和傳感器的個數相同

2、源信號相互統計獨立

3、源信號各矢量均值爲0,至多有一個是高斯信號(當源信號都是高斯信號時,它們的混合信號仍是高斯信號,它們是無法進行分離的)

 

自適應濾波理論和技術是統計信息處理和非平穩隨機信號處理的主要內容,在不需要先驗知識的初始條件下,通過自學習來適應外部的自然隨機環境,因而自適應算法可以用來估計確定信號。

1、最陡下降法(Steepest Descent Method)

沿性能曲面最陡方向向下搜索曲面的最低點。曲面的最陡下降方向是曲面的負梯度方向。

Y(n)=W(n)X(n)W(n)是濾波器係數矩陣,X(n)是系統輸入矢量,Y(n)是輸出信號,期望信號爲d(n)

誤差序列e(n)=d(n)-Y(n)

按照均方誤差準則

F(e(n))=E[e^2(n)]]=E[d^2(n)-2d(n)Y(n)+Y^2(n)] =E[d^2(n)]-2E[d(n)W^T(n)X(n)]+E[W^T(n)X(n)X^T(n)W(n)]

當濾波器係數固定時,目標函數可以寫成

\xi (n)=E[d^2(n)]-2W^TP+W^TRW

其中,P=E[d(n)X(n)]是期望信號與輸入信號的互相關矢量,R=E[X(n)X^T(n)]是輸入信號的自相關矩陣

由上式可知,自適應算法的目標函數是延時線抽頭係數的二次函數,當矩陣R和矢量P已知時,可以由權矢量W直接對其求解,對上式求導可得到目標函數最小的最佳濾波係數

w_{0}=R^{-1}P

濾波係數更新W(n+1)=W(n)+1/2\mu [-\bigtriangledown (n)]

1/2表示\mu減半,《現代數字信號處理》(姚天任)一書中沒有減半,即W(n+1)=W(n)-\mu [\bigtriangledown (n)]

\bigtriangledown (n)=2RW(n)-2P

若按\mu減半的算,濾波器係數更新值爲W(n+1)=W(n)+\mu [P-RW(n)]

2、最小均方誤差(LMS)算法

3、最小二乘法(LS)

 

 

 

轉自:https://www.docin.com/p-126816392.html?docfrom=rrela

發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章