一種單通道混沌信號盲源分離法
【技術(shù)領(lǐng)域】
[0001] 本發(fā)明屬于電子信息技術(shù)領(lǐng)域,具體涉及一種單通道混沌信號盲源分離法。
【背景技術(shù)】
[0002] 混沌信號的波形是非常不規(guī)則的,表面上看來就像噪聲,但實際上它卻是由確定 性的規(guī)則所產(chǎn)生的,這種規(guī)則有時是很簡單的。正是這種簡單的規(guī)則產(chǎn)生出復(fù)雜的波形激 發(fā)了人們對它極大的興趣。自從20世紀(jì)70年代以來,混沌學(xué)已經(jīng)成為一門新型的科學(xué)分 支。
[0003] 對于多路混沌信號的傳輸,有些采用多接口和連線,但花費成本較高;或者利用復(fù) 用和解復(fù)用技術(shù),但這些方法的缺點是采用設(shè)備復(fù)雜度較高。所以現(xiàn)實生活中,由于條件及 成本限制,往往采用單通道盲源分離法,即在多路源信號混合后,憑借單通道混合信號再恢 復(fù)出多路源信號。
[0004] 目前對于單通道盲源分離方法主要有三種類型:基于模型法、虛擬多通道法和針 對已調(diào)信號特征的方法。由于虛擬多通道法在算法的復(fù)雜度和收斂速度等方面表現(xiàn)出較好 的特性,其無需依賴信號的先驗概率特征或者使用的調(diào)制方式,直接將單通道混合信號投 射到高維空間,再與經(jīng)典的獨立成分分析(Independent Component Analysis, ICA)算法結(jié) 合恢復(fù)出源信號,所以該方法針對單通道盲源分離這種極端欠定的情況具有很大的應(yīng)用價 值。
[0005] 總體經(jīng)驗?zāi)B(tài)分解(Ensemble Empirical Mode Decomposition, EEMD)是將單通 道混合信號分解為多個本征模態(tài)函數(shù)(MFs),構(gòu)建出虛擬多通道,此方法相對于空時法、 小波分解法等分離效果得到明顯提高,對于頻譜有重疊的信號也能夠得到較好分離恢復(fù)效 果,但在實際實驗中發(fā)現(xiàn),當(dāng)直接用ICA處理MFs分量時,經(jīng)常遇到算法迭代次數(shù)過高、收 斂速度慢等問題。因此利用EEMD_ICA方法處理過程中需要人憑借經(jīng)驗進(jìn)行信號的挑選,但 智能性不高。
【發(fā)明內(nèi)容】
[0006] 本發(fā)明目的是提供一種單通道混沌信號盲源分離法,可以有效地克服現(xiàn)有技術(shù)的 缺點。
[0007] 本發(fā)明是這樣實現(xiàn)的,其特征在于采用的模塊包括有總體經(jīng)驗?zāi)B(tài)分解模塊、相 關(guān)性分析模塊、主成分分析模塊、獨立成分分析模塊,具體實施步驟是:
[0008] A、將多路混沌信號通過屏蔽導(dǎo)線送到混沌信號預(yù)處理模塊直接相加,得到預(yù)處理 單通道混沌信號x(t);
[0009] B、將混合得到的預(yù)處理單通道混沌信號x(t)送到混沌信號盲源分離模塊,進(jìn)行 總體經(jīng)驗?zāi)B(tài)分解,即EEMD分解,基于相關(guān)系數(shù)的主成分分析PCA降維和獨立成分分析,即 ICA分析,實現(xiàn)多路混沌信號通過一個輸入口采集,多路輸出口輸出;
[0010] B. 1、總體經(jīng)驗?zāi)B(tài)分解EEMD處理得到本征模態(tài)函數(shù)MF ;
[0011] a、將混合得到的預(yù)處理單通道混沌信號x(t)分別與i = 1,2, ···,N路不同的白噪 聲叫(0疊加,加入的白噪聲njt)的均值為零、標(biāo)準(zhǔn)差為常數(shù),得到加入白噪聲之后的信號 xi⑴,其中i = 1,2,…,N :
[0012]
【主權(quán)項】
1. 一種單通道混沌信號盲源分離法,其特征在于采用的模塊包括有總體經(jīng)驗?zāi)B(tài)分解 模塊、相關(guān)性分析模塊、主成分分析模塊、獨立成分分析模塊,實施步驟是: A、 將多路混沌信號通過屏蔽導(dǎo)線送到混沌信號預(yù)處理模塊直接相加,得到預(yù)處理單通 道混沌信號x(t); B、 將混合得到的預(yù)處理單通道混沌信號x(t)送到混沌信號盲源分離模塊,進(jìn)行總體 經(jīng)驗?zāi)B(tài)分解,即EEMD分解,基于相關(guān)系數(shù)的主成分分析PCA降維和獨立成分分析,即ICA 分析,實現(xiàn)多路混沌信號通過一個輸入口采集,多路輸出口輸出; B. 1、總體經(jīng)驗?zāi)B(tài)分解EEMD處理得到本征模態(tài)函數(shù)頂F ; a、 將混合得到的預(yù)處理單通道混沌信號x(t)分別與i = 1,2,···,Ν路不同的白噪聲 Hi (t)疊加,加入的白噪聲Ili (t)的均值為零、標(biāo)準(zhǔn)差為常數(shù),得到加入白噪聲之后的信號 Xi(t),其中 i = 1,2,…,N :
b、 對所得到的信號Xi⑴分別進(jìn)行經(jīng)驗?zāi)B(tài)分解EMD,得到各自的MF記為aij⑴和余 項巧⑴,其中 aiJ(t)表示加入白噪聲叫⑴后分解得到的第j個頂F:
c、 對經(jīng)EMD分解后所得到的MF集合{aij (t),…,aij (t),…aNj (t)}進(jìn)行總體平均運 \ N 算,得到的IMF為α _ (〇 =--Σ (0,即為對原信號進(jìn)行EEMD分解后所得到的第j個 IMF,j = 1,2, aj(t)的序列長度為n,由此得到IMF的mXn分量矩陣A = {ajt),… ,a^t),其中a」(t) e A,m為單通道混沌信號經(jīng)EEMD分解后的MF個數(shù),上標(biāo)T 為轉(zhuǎn)置運算; B. 2、對得到的IMF分量進(jìn)行降維: a、對得到的IMF分量根據(jù)相關(guān)度進(jìn)行第一次降維: 對得到的m個MF分量分別與混合得到的預(yù)處理單通道混沌信號x(t)求相關(guān)度'
,其中D (X)為X (t)的方差,D (a』)為aj(t)的方差,cov (X)為X (t) 的協(xié)方差,cot (ap為\_(〇的協(xié)方差,取相關(guān)度大于閾值的IMF分量,閾值為常數(shù),進(jìn) 行第一次降維,得到IMF的kXn分量矩陣B= {ajt), …七⑴}1,其中aj(t) eB, k為第一次降維后的頂F個數(shù),上標(biāo)T為轉(zhuǎn)置運算; b、對第一次降維所得到的矩陣B進(jìn)行PCA降維: R = E(BBT),RV = VA,其中B = {ai(t),…七⑴,…以⑴}1為第一次降維后的MF 分量k X η矩陣,R為k個變量IMF的自相關(guān)矩陣,V為R的k X k階特征向量矩陣,其列向量 是R的正交歸一化的特征向量,Λ為R的特征對角矩陣,λ i (i = 1,2,…,k)為Λ第i個 對角線上的元素,構(gòu)造 k個不相關(guān)的新變量Y = VTB,其中Y= Iy^y2,…"!^,對AiQ = 1,2,…,k)按降序排列后,取前面p個較大特征值所對應(yīng)的特征向量,得到kXk階特征向 量矩陣V降維后的pXn階向量矩陣C,其中p彡2 ; B. 3、將PCA降維所得到矩陣C進(jìn)行ICA處理,采用約束獨立成分分析(CICA)算法進(jìn)行 處理: a、 對C進(jìn)行白化處理,得到白化數(shù)據(jù),即觀測矩陣X,其中X = MS (n),M為信號的混合 矩陣,S (η)為源信號; b、 隨機選取權(quán)向量作為解混矩陣W的權(quán)值,即W的初始值,W為解混矩陣,同時也是混 合矩陣M的虛擬反矩陣,因為W = 所以S (n) = IT1X = WX ; c、 為了求出最終收斂的W值,根據(jù)牛頓迭代法推導(dǎo),得到迭代式 W -E(Xg(WTX))-E(g' (WTX))W,其中 X 為 C 白化后的觀測矩陣,g(u) = tanh(alU)為非二 次函數(shù)的導(dǎo)數(shù),1彡S1S 2 ; d、 循環(huán)迭代,直到收斂為止,最后得到pXn階向量矩陣Y(n),通過觀察選取分離后的 混沌信號。
【專利摘要】一種單通道混沌信號盲源分離法,屬于電子信息技術(shù)領(lǐng)域,其特征在于采用的模塊包括有總體經(jīng)驗?zāi)B(tài)分解模塊、相關(guān)性分析模塊、主成分分析模塊、獨立成分分析模塊。實施步驟是將多路混沌信號混合成的單路信號,通過經(jīng)驗?zāi)B(tài)分解模塊分解為多路本征模態(tài)函數(shù)分量,利用相關(guān)性分析對多路本征模態(tài)函數(shù)的主元進(jìn)行一次降維,再用主成分分析法對多路本征模態(tài)函數(shù)的主元進(jìn)行二次降維,最后采用獨立成分分析技術(shù)完成混合混沌信號的盲源分離,恢復(fù)出源信號。優(yōu)點是能夠有效地減少傳輸通道數(shù)量,達(dá)到降低硬件復(fù)雜度,節(jié)約設(shè)備成本的目的,并能快速有效地恢復(fù)出源信號。
【IPC分類】G06F17-16
【公開號】CN104636314
【申請?zhí)枴緾N201510049797
【發(fā)明人】郭一娜, 莫曉敏, 王曉梅, 杜雅梅, 田文艷, 卓東風(fēng)
【申請人】太原科技大學(xué)
【公開日】2015年5月20日
【申請日】2015年1月30日