基于聯(lián)合稀疏表示和同時(shí)稀疏近似的并行坐標(biāo)下降去噪算法

何選森 徐 麗

1(廣州商學(xué)院信息技術(shù)與工程學(xué)院 廣東 廣州 511363)2(湖南大學(xué)信息科學(xué)與工程學(xué)院 湖南 長(zhǎng)沙 410082)

0 引 言

信號(hào)的稀疏表示(Sparse representation)[1]是將信號(hào)分解為超完備字典(Over-complete dictionary)[2]少量原子的線性組合。由于觀測(cè)信號(hào)存在噪聲,信號(hào)處理的任務(wù)之一是去噪(Denoising)。而稀疏表示具有潛在的降噪能力[3],聯(lián)合稀疏表示(Joint Sparse Representation,JSR)[4]和同時(shí)稀疏近似(Simultaneous Sparse Approximation,SSA)[5]已被廣泛地應(yīng)用于音頻和圖像去噪。收縮(Shrinkage)則是基于稀疏表示的著名降(去)噪技術(shù)[6],如經(jīng)典的小波收縮是利用正交變換實(shí)現(xiàn)降噪[7],而新的趨勢(shì)是采用超完備變換替代正交變換。另外,基本追蹤(Basis Pursuit,BP)[8]是通過優(yōu)化將信號(hào)分解為字典原子的最佳疊加,從而使表示系數(shù)的l1范數(shù)(或l0范數(shù))最小化[9],而基本追蹤去噪(BP Denoising,BPDN)[10]對(duì)于消除高斯白噪聲是最佳的。隨著超完備字典的發(fā)展,迭代收縮(Iterative Shrinkage,IS)算法[11]已成為一種高效的數(shù)值技術(shù)。Elad[12]提出了四種IS算法以解決BPDN問題,其中的并行坐標(biāo)下降(Parallel Coordinate Descent,PCD)算法[12]具有最快的搜索速度。文獻(xiàn)[13]對(duì)PCD進(jìn)行了收斂性分析,并引入了一種正則化函數(shù)形式,使得PCD能對(duì)坐標(biāo)優(yōu)化提供解析的解。

對(duì)于音頻信號(hào),用PCD分別處理每一幀能夠獲得很好的去噪性能。然而,當(dāng)信號(hào)很長(zhǎng)時(shí),分割的音頻幀數(shù)量會(huì)很大,用PCD對(duì)每幀進(jìn)行降噪使得計(jì)算量急劇增加,造成極大的運(yùn)行時(shí)間成本。為此,本文通過構(gòu)建以信號(hào)(矩陣)為操作對(duì)象的時(shí)域處理框架,提出一種啟發(fā)式矩陣型PCD算法(Joint-PCD)。在新的時(shí)域框架中,把信號(hào)中的每個(gè)語音幀看作一個(gè)列向量生成信號(hào)矩陣,Joint-PCD算法是對(duì)一個(gè)信號(hào)(矩陣)實(shí)施降噪,而不僅僅只對(duì)一幀處理,因而可以降低算法的運(yùn)行時(shí)間成本。

1 稀疏信號(hào)的降噪

若在同一時(shí)刻,信號(hào)向量x∈RN中只有很少分量是非零的,其他絕大多數(shù)分量都為零值,則稱x服從稀疏分布,即x可由幾個(gè)單位能量的原子dk的線性組合來近似[14]。由所有原子構(gòu)成的集合{dk,k=1,2，…,K}稱為字典D∈RN×K。設(shè)向量z∈RK是稀疏的,則信號(hào)x可以表示為：

x=Dz

(1)

式中:z是由x的K個(gè)重構(gòu)系數(shù)組成的列向量。當(dāng)字典D是正交的,K=N,則表示是唯一的。當(dāng)字典D是冗余(超完備)的,K>N,則重構(gòu)每個(gè)信號(hào),就有無限多個(gè)可能的系數(shù)集。正是這種非唯一性提供了表示方法的可選擇性。一種更好的表示法是選擇具有最少的非零系數(shù)的解[15]。

s.t.x=Dz

式(2)是尋找欠定(K>N)線性方程組x=Dz的最稀疏解[15]。類似地,式(2)也可以表示為[15]：

s.t.x=Dz

式中:l1范數(shù)是求向量各元素絕對(duì)值之和。由于l1范數(shù)是l0范數(shù)的凸化[15],因此,式(4)也可看作是式(2)的凸化。

對(duì)于含噪聲的觀測(cè)信號(hào)向量:

y=x+v

(5)

式中:x∈RN是不含噪聲的干凈信號(hào)；v∈RN是有限能量(表示為δ)的高斯零均值白噪聲向量。對(duì)含噪聲信號(hào)y進(jìn)行去噪處理的優(yōu)化問題可表示為[16]:

式(6)還可以表示為拉格朗日乘數(shù)形式[9]:

式中:λ是正的參數(shù),它用于平衡近似誤差與系數(shù)向量稀疏性二者之間的關(guān)系。本質(zhì)上,式(7)就是基于含噪聲信號(hào)向量y和超完備字典D的PCD算法的目標(biāo)函數(shù)[12]。

從y中去除高斯白噪聲是對(duì)以下函數(shù)最小化[12]:

式中:等號(hào)右邊的第一項(xiàng)為干凈信號(hào)x與觀測(cè)信號(hào)y之間的對(duì)數(shù)似然；第二項(xiàng)的ρ(x)表示在x上的先驗(yàn)值。f(x)的展開式與ρ(x)有關(guān)。將式(1)代入式(8)中得到[12]:

式中:1K表示元素為1的K×1列向量。

IS就是解決式(9)所對(duì)應(yīng)的混合l2-lp(p≤1)優(yōu)化問題的算法。比較式(9)和式(7)可知,在PCD算法中,先驗(yàn)函數(shù)ρ(z)是l1范數(shù),而且PCD的操作對(duì)象是向量。

2 新的時(shí)域處理框架

用超完備字典原子的線性組合來表示信號(hào)向量,就稱為簡(jiǎn)單稀疏近似(Simple sparse approximation),在PCD算法中使用的就是簡(jiǎn)單稀疏近似。對(duì)于很長(zhǎng)的音頻信號(hào),其分割的幀數(shù)就很多,利用PCD分別處理每一幀將造成算法運(yùn)行時(shí)間成本的劇增。為此,需要將以向量為操作對(duì)象的模式拓展成以信號(hào)(矩陣)為操作對(duì)象的模式。

對(duì)多個(gè)向量同時(shí)稀疏近似理論的分析和研究[17],產(chǎn)生出求解SSA問題的策略[18]。在通信系統(tǒng)中,對(duì)于被同一信道噪聲污染的音頻信號(hào),其分割的(幀)向量也都具有相同的分布特性。對(duì)不同的幀向量,可以使用同一個(gè)字典的不同原子的線性組合來近似計(jì)算,即以聯(lián)合稀疏的模式促進(jìn)它們非零分量的耦合。因此,為了同時(shí)(或并行)地對(duì)一個(gè)信號(hào)的所有幀(向量)進(jìn)行降噪操作,需要將傳統(tǒng)的幀處理模式轉(zhuǎn)變成如下新的時(shí)域處理框架。

根據(jù)分割規(guī)則,將音頻信號(hào)分割成一系列的幀,每個(gè)幀作為一個(gè)列向量生成矩陣,使得信號(hào)與矩陣唯一的對(duì)應(yīng)。由于一個(gè)矩陣中所有的列向量對(duì)應(yīng)于同一個(gè)音頻信號(hào)的不同幀,因此這些列向量具有相同的稀疏性。另外,由于每個(gè)列向量都受到同一個(gè)信道噪聲的污染,因而不同噪聲分量之間也不需要統(tǒng)計(jì)獨(dú)立?；谶@樣的信號(hào)矩陣,可以同時(shí)處理多個(gè)列向量以節(jié)省降噪的時(shí)間開銷。

3 Joint-PCD啟發(fā)式算法

在新的時(shí)域處理框架之下,將簡(jiǎn)單稀疏近似的概念推廣到同時(shí)稀疏近似(SSA),并以此為基礎(chǔ),將PCD擴(kuò)展為矩陣型的Joint-PCD算法,其基本原理如下。

令y1=x1+v1、y2=x2+v2分別為干凈信號(hào)x1,x2∈RN的含噪聲觀測(cè)量,其中v1,v2∈RN為加性噪聲。假設(shè)信號(hào)x1、x2可以由超完備字典D∈RN×K近似表示為x1=Dz1、x2=Dz2,其中z1,z2∈RK為x1、x2的重建系數(shù)。由式(7)可得:

聯(lián)合考慮以上這兩個(gè)優(yōu)化問題,則有：

(11)

再令X=[x1,x2]、Y=[y1,y2]、Z=[z1,z2],則將式(11)簡(jiǎn)化為如下標(biāo)準(zhǔn)的稀疏編碼問題:

對(duì)Z和D二者來說,式(12)不是一個(gè)凸優(yōu)化,但當(dāng)固定其中一個(gè)時(shí),式(12)就成為凸優(yōu)化問題[4]。本文中,字典D是給定的,而通過優(yōu)化系數(shù)矩陣Z來構(gòu)造Joint-PCD算法。

式(12)是對(duì)兩個(gè)信號(hào)y1、y2同時(shí)去噪的優(yōu)化問題。顯然,優(yōu)化式(12)可以推廣到多個(gè)信號(hào)向量的情況。考慮矩陣形式的觀測(cè)信號(hào)模型為[19]：

Y=X+V

(13)

式中:X∈RN×P是包含P個(gè)向量xi∈RN的干凈信號(hào)矩陣；V∈RN×P是P個(gè)向量vi∈RN組成的噪聲矩陣；Y∈RN×P是含噪聲觀測(cè)信號(hào)矩陣?；诼?lián)合稀疏表示(JSR)[4],信號(hào)X可以由字典D表示為X=DZ,于是觀察信號(hào)矩陣Y為[19]：

Y=DZ+V

(14)

式中:ρ(·)是任意標(biāo)量函數(shù),而ρ(Z)是矩陣Z稀疏性的懲罰函數(shù)。式(16)就是Joint-PCD算法的目標(biāo)函數(shù)。

由式(15)和式(16)可歸納出構(gòu)造Joint-PCD算法的基本思想。為了限制對(duì)觀測(cè)信號(hào)Y進(jìn)行逼近的字典原子的總數(shù),需要使Z中非零行向量的數(shù)量最小化。為此,從兩方面來考慮。首先,參與對(duì)Y近似的字典原子總數(shù)越少越好;其次,每個(gè)參與近似的原子要對(duì)盡可能多的Y的列向量做出貢獻(xiàn)。也就是說,系數(shù)矩陣Z是稀疏的,即它的大多數(shù)行向量都是零向量,而Z的每個(gè)非零行向量都應(yīng)該有盡可能多的非零元素。基于此,將稀疏性懲罰函數(shù)ρ(Z)設(shè)置為行l(wèi)0準(zhǔn)范數(shù)(Row-l0quasi-norm)的松弛版:

把式(17)代入到式(16)中,則觀測(cè)信號(hào)矩陣Y的降噪任務(wù)變成為混合范數(shù)的優(yōu)化問題:

同樣地,由式(12)可知,懲罰函數(shù)ρ(Z)是矩陣Z的l1范數(shù)。在本文中,對(duì)ρ(Z)的l1范數(shù)使用偶對(duì)稱函數(shù):

式中:zij是矩陣Z的第(i,j)個(gè)元素。

這里需要提示,在式(12)中Z和Y都是由兩個(gè)向量構(gòu)成的矩陣,然而,在下文中,將兩個(gè)向量情況推廣到多個(gè)向量的情況,即式(12)中的Z和Y是由P(P≤2)個(gè)列向量組成的矩陣Z∈RK×P和Y∈RN×P,而不僅僅是由2個(gè)列向量組成的矩陣Z∈RK×2和Y∈RN×2。

U(D,Y,Zk)=diag-1{DTD}DT(Y-DZk)+Zk

(20)

一般地,超完備字典D的列是歸一化的,即diag-1(

DTD)=I,其中I是單位矩陣,則可將式(20)簡(jiǎn)化為:

U(D,Y,Zk)=DT(Y-DZk)+Zk

(21)

矩陣Z的更新迭代值Zk+1為:

Zk+1=Zk+μQk

(22)

式中：μ是迭代步長(zhǎng)參數(shù)。

Qk=S[U(D,Y,Zk),λ]-Zk

(23)

式中：λ是收縮函數(shù)S(·)的閾值。從四種IS算法的推導(dǎo)過程[12]可知,式(23)中的收縮函數(shù)S(·)是一個(gè)軟閾值函數(shù),其定義為：

式中:θ為軟收縮函數(shù)S(·)的閾值。

另外,Joint-PCD算法的性能會(huì)受到式(22)中步長(zhǎng)參數(shù)μ影響。對(duì)足夠小的μ>0,步長(zhǎng)要保證懲罰函數(shù)ρ(·)的可行下降。本文使用線搜索[11]來尋找最佳的步長(zhǎng)μ值。線搜索作為一維優(yōu)化,它用于解決以下的問題:

式中:ρ(·)可以是矩陣的rx范數(shù),也可以是矩陣的l1范數(shù)。

若信號(hào)矩陣Y退化成一個(gè)列向量y,系數(shù)矩陣Z也退化成一個(gè)列向量z,則Joint-PCD算法退化為PCD算法。因此,將Joint-PCD稱為啟發(fā)式算法。

在音頻信號(hào)處理中,經(jīng)常使用離散余弦變換(DCT)字典和Gabor字典。一個(gè)窗口式DCT字典為DDCT=[d1,d2,…,dK],它的第k個(gè)原子定義為[20]：

一個(gè)窗口式Gabor字典為DG=[d(k,φ)](k,φ)∈Ψ,其原子以連續(xù)的參數(shù)集Ψ=[1,K]×[0,2π]為索引,其定義為[20]：

式中:1≤k≤K,0≤t≤N-1,N是音頻幀的長(zhǎng)度,K是字典的尺寸；wd是加權(quán)窗口；φ是相位。Gabor字典連續(xù)索引原子d(k,φ)的分解可以用離散字典中的原子對(duì)來表示[20]。

式中：原子對(duì)是相同頻率且相位為零的余弦和正弦對(duì)。為簡(jiǎn)化計(jì)算,詞典的列都做歸一化處理,如Gabor字典的正弦和余弦原子都使用對(duì)應(yīng)的單位范數(shù)版本。

綜上所述,對(duì)啟發(fā)式Joint-PCD算法描述見算法1。

輸入:觀測(cè)矩陣Y，字典D，閾值λ，函數(shù)ρ(·)，迭代次數(shù)M。

Begin:設(shè)置k=1,迭代開始。

ifk

誤差:計(jì)算E=Y-DZk;

映射:計(jì)算ET=DTE;

收縮:計(jì)算ES=S(ET+Zk,λ);

線收縮:尋找μ0以獲得Zk+μ(ES-Zk)的下降;

放松:更新迭代Zk+1=Zk+μ0(ES-Zk);

Next:k←k+1;

End

4 仿真實(shí)驗(yàn)及結(jié)果分析

4.1 仿真環(huán)境和性能指標(biāo)

為了驗(yàn)證Joint-PCD算法降噪的有效性,對(duì)Joint-PCD和PCD在音頻信號(hào)降噪性能上進(jìn)行對(duì)比分析。不同算法的仿真環(huán)境是完全一樣的,以便得到公平的結(jié)論。

仿真平臺(tái)為筆記本電腦,其CPU為Intel(R) Celeron(R) 1007U CPU-1.50 GHz,內(nèi)存4 GB,操作系統(tǒng)為Windows 10。仿真軟件平臺(tái)為MATLAB 9(R2016a)。

降噪的性能用信噪比(signal to noise ratio,SNR)和均方根誤差(root mean square error,RMSE)來衡量:

測(cè)試用的五個(gè)音頻音樂信源[21]的類型為.wav,信號(hào)采樣率為44.1 MHz。每個(gè)信源的樣本數(shù)為L(zhǎng)=216=65 536。五個(gè)信源均服從超高斯分布,它們之間的區(qū)別體現(xiàn)在具有不同的峭(峰)度(kurtosis)值。信源X=[x1,x2,x3,x4,x5]T的歸一化峰度值分別為2.400 8、4.776 6、1.052 8、4.690 2、3.133 7,其對(duì)應(yīng)的時(shí)域波形如圖1所示。

圖1 五路音樂音頻源信號(hào)的波形

降噪的仿真過程如下。音頻信源與高斯噪聲相加生成含噪聲信號(hào)。首先,每個(gè)含噪聲時(shí)域信號(hào)乘以長(zhǎng)度為64的漢明窗生成語音幀,連續(xù)兩幀之間的重疊樣本數(shù)為32,即每幀的前后各有16個(gè)樣本與相鄰幀的樣本重疊。分幀過程是通過使用Voicebox[22]中的MATLAB函數(shù)enframe來實(shí)現(xiàn)。然后,將PCD和Joint-PCD算法分別應(yīng)用于分割后的信號(hào)以進(jìn)行降噪。每次使用PCD是對(duì)音頻信號(hào)的一幀進(jìn)行處理;而使用Joint-PCD算法每次是對(duì)一個(gè)信號(hào)(而不僅是一幀)進(jìn)行處理。最后,重建降噪后的信號(hào),即每個(gè)降噪的幀乘以逆漢明窗,并且在相鄰兩幀重疊的樣本中,分別在兩邊去除50%的樣本,僅保留每一幀的中間部分的樣本,將經(jīng)過處理的各個(gè)片段簡(jiǎn)單串聯(lián)起來即可得到重新合成的音頻信號(hào)。

4.2 算法參數(shù)設(shè)置

在PCD算法中,需要給定軟收縮函數(shù)S(·)的閾值λ。常用的閾值規(guī)則為全局閾值(universal threshold)[7]和最小最大閾值(minimax threshold)[7]。本文采用全局閾值的設(shè)置,滿足以下條件:

式中:L是信號(hào)的樣本數(shù)量；σ是高斯噪聲的標(biāo)準(zhǔn)差。顯然,要確定閾值λ,則必須先確定噪聲的強(qiáng)度σ。

為了確定噪聲強(qiáng)度對(duì)降噪性能的影響,采用經(jīng)典的小波降噪方法進(jìn)行仿真。對(duì)5個(gè)信源X加噪聲形成觀測(cè)信號(hào)Y,當(dāng)Y的信噪比分別為10、15、20時(shí),利用小波重復(fù)進(jìn)行100次降噪實(shí)驗(yàn),并計(jì)算在每次實(shí)驗(yàn)中5個(gè)信號(hào)的平均SNR值和RMSE值,其結(jié)果如圖2所示。

圖2 5個(gè)信號(hào)的平均SNR和RMSE值

從圖2可知,隨著輸入SNR的降低(噪聲強(qiáng)度的增大),小波的降噪性能變差。當(dāng)輸入SNR為10時(shí),經(jīng)過降噪后5個(gè)信號(hào)的平均SNR小于9.8,即降噪后的信噪比小于輸入的信噪比。對(duì)音頻信號(hào)來說,這是一個(gè)非常惡劣的環(huán)境。

為了便于分析,對(duì)不同的輸入SNR,計(jì)算出每個(gè)信源被施加高斯噪聲的強(qiáng)度σ。其結(jié)果如表1所示。

表1 輸入SNR與噪聲強(qiáng)度σ的對(duì)應(yīng)關(guān)系

從表1中的數(shù)據(jù)可知,由于每個(gè)信源的幅度(功率)不同,對(duì)同一個(gè)輸入SNR值,其對(duì)應(yīng)于不同信源的噪聲強(qiáng)度也是不同的。對(duì)于輸入SNR為10的情況,信號(hào)x2、x3、x4、x5上施加的噪聲強(qiáng)度接近0.1,但信號(hào)x1上施加的噪聲強(qiáng)度為0.2。為獲得一個(gè)統(tǒng)一的噪聲強(qiáng)度,在下面的仿真中,5個(gè)信源都被施加σ=0.1的高斯噪聲形成觀測(cè)信號(hào)。

給定噪聲強(qiáng)度后,再來確定軟閾值λ。由于信號(hào)采樣點(diǎn)數(shù)為L(zhǎng),則(2log2L)1/2=4.709 6,由文獻(xiàn)[7]可知,在給定σ=1的情況下,其參考的閾值為λ*=3.31。將真實(shí)的σ=0.1代入可得λ≤σ(3.31=0.331,這也只是λ的取值范圍,還需要通過仿真來尋找PCD算法最佳的閾值λ。仿真條件為：每個(gè)含噪信號(hào)都被分割成長(zhǎng)度為64的幀,PCD的迭代次數(shù)設(shè)置為M=5,閾值λ從0.050到0.331變化,更新步長(zhǎng)為0.01,每個(gè)λ值做一次實(shí)驗(yàn),超完備字典是DCT和Gabor。各個(gè)信號(hào)的性能指標(biāo)SNR和RMSE如圖3所示。

(a) DCT

(b) Gabor圖3 每個(gè)信號(hào)對(duì)不同λ的SNR和RMSE值

從圖3可知,無論是DCT還是Gabor詞典,當(dāng)λ=0.1時(shí),每個(gè)信號(hào)的SNR達(dá)到最大值,而RMSE也達(dá)到對(duì)應(yīng)的最小值。因此在PCD算法中,取λ=0.1作為最佳的閾值。另外,從圖3還可看出,采用不同的超完備字典,PCD算法的降噪性能是有差異的。為了比較DCT和Gabor的性能,將迭代次數(shù)設(shè)置為M=15,并分別計(jì)算PCD對(duì)每個(gè)信號(hào)在每次迭代中的SNR和RMSE值。其結(jié)果如圖4所示。

圖4 每個(gè)信號(hào)對(duì)不同字典的SNR和RMSE值

為便于分析,僅觀察信號(hào)1(signal1)的情況。從圖4可知,使用Gabor字典的SNR和RMSE值都優(yōu)于DCT字典對(duì)應(yīng)的指標(biāo)。其他信號(hào)的SNR和RMSE值也有相似的關(guān)系。為此,在PCD和Joint-PCD算法中,使用Gabor字典。

另外,在算法的參數(shù)中,還需要確定懲罰函數(shù)ρ(·)及迭代步長(zhǎng)μ的值。在下面的仿真中,對(duì)于PCD算法,采用向量的l1范數(shù)求解式(7);對(duì)于Joint-PCD算法,采用矩陣的l1范數(shù)求解式(12),采用矩陣的rx范數(shù)求解式(18)。為了提高算法的性能,各算法的步長(zhǎng)μ(0<μ<1)采用MATLAB函數(shù)fminbnd通過在線搜索最優(yōu)的μ值作為其迭代步長(zhǎng)。這樣就完成了PCD和Joint-PCD算法中所有參數(shù)的設(shè)置。

4.3 算法性能比較

在對(duì)算法去噪性能的比較過程中,為方便起見,把rx范數(shù)和l1范數(shù)的Joint-PCD分別表示為Joint-PCD withrxnorm和Joint-PCD withl1norm。

關(guān)于算法的性能,從兩個(gè)方面對(duì)PCD和Joint-PCD進(jìn)行比較。(1) 在相同條件下算法所消耗的運(yùn)行時(shí)間成本;(2) 去噪的性能指標(biāo)值SNR和RMSE。在這個(gè)仿真實(shí)驗(yàn)中,將算法的最大迭代次數(shù)M設(shè)置為50以進(jìn)行全面的考察。其運(yùn)行時(shí)間如圖5所示。

圖5 算法的運(yùn)行時(shí)間比較

可以看出,算法的運(yùn)行時(shí)間隨迭代次數(shù)的增加基本上呈線性增長(zhǎng)。顯然,rx范數(shù)和l1范數(shù)的Joint-PCD算法的運(yùn)行時(shí)間增長(zhǎng)率幾乎是相同的;而PCD算法的運(yùn)行時(shí)間增長(zhǎng)率卻遠(yuǎn)遠(yuǎn)大于Joint-PCD的增長(zhǎng)率。

圖6分別給出了PCD算法、Joint-PCD withl1norm和Joint-PCD withrxnorm算法在該仿真中去噪性能指標(biāo)值SNR和RMSE的波形圖。

(a) PCD

(b) Joint-PCD with l1 norm

(c) Joint-PCD with rx norm圖6 PCD和Joint-PCD算法對(duì)每個(gè)信號(hào)的SNR和RMSE值

可以看出,PCD與Joint-PCD算法的降噪性能幾乎是完全相同的,而且經(jīng)過大約25次迭代,三種算法的降噪指標(biāo)SNR和RMSE值都逐漸收斂到其穩(wěn)定值。

為了比較算法的平均降噪性能,將每種算法對(duì)5個(gè)信號(hào)的平均SNR和平均RMSE值繪制在圖7中。

(a)

(b)圖7 5個(gè)信號(hào)的平均SNR和RMSE指標(biāo)

可以看出,PCD和Joint-PCD算法的平均SNR和RMSE值在算法的初始階段隨著迭代次數(shù)的增加具有一定的波動(dòng),這是因?yàn)槊總€(gè)信號(hào)的SNR和RMSE值具有不同的波動(dòng)范圍。對(duì)于PCD和Joint-PCD withrxnorm算法,平均降噪指標(biāo)大約在經(jīng)過20次迭代后結(jié)束波動(dòng),趨于平穩(wěn);而Joint-PCD withl1norm算法大約在經(jīng)25次迭代后其降噪指標(biāo)趨于平穩(wěn)。也就是說,PCD和Joint-PCD的平均降噪性都能在25次迭代后收斂到穩(wěn)定的指標(biāo)值。

為了從數(shù)值上直觀地比較PCD和Joint-PCD在降噪性能上的差異,在以下仿真中,將算法的最大迭代次數(shù)M設(shè)置為25。表2給出PCD和Joint-PCD的運(yùn)行時(shí)間。

表2 PCD和Joint-PCD算法的運(yùn)行時(shí)間 單位:s

可以看出,l1范數(shù)的Joint-PCD運(yùn)行時(shí)間是最少的,而rx范數(shù)的Joint-PCD運(yùn)行時(shí)間也僅僅是多了0.379 2 s,因此,兩種Joint-PCD算法的運(yùn)行時(shí)間幾乎是相同的。相反地,PCD算法的運(yùn)行時(shí)間大約是Joint-PCD的4.6倍,這說明在音頻信號(hào)的降噪過程中,Joint-PCD算法降低了計(jì)算負(fù)擔(dān)。顯然,與PCD算法相比較,Joint-PCD算法的收斂速度得到了有效的提高。

對(duì)于PCD和Joint-PCD算法的降噪性能,采用每個(gè)信號(hào)在整個(gè)迭代過程中SNR和RMSE的平均值來度量。表3給出了每個(gè)信號(hào)的平均降噪性能指標(biāo)。

表3 每個(gè)信號(hào)在整個(gè)迭代過程的平均SNR和RMSE值

可以看出,Joint-PCD withl1norm算法對(duì)于第2、第4和第5個(gè)信號(hào)的平均SNR和RMSE值是最優(yōu)的;而PCD算法僅對(duì)于第1和第3個(gè)信號(hào)的平均指標(biāo)值是最優(yōu)的。另外,PCD算法和Joint-PCD算法的平均SNR值的差別僅體現(xiàn)在小數(shù)點(diǎn)后面第二位;而算法的平均RMSE值的差別僅發(fā)生在小數(shù)點(diǎn)后面第四位上。因此,從去噪性能指標(biāo)的差異性來看,利用Joint-PCD算法代替PCD算法對(duì)音頻信號(hào)進(jìn)行降噪處理,其降噪的效果幾乎完全相同。

5 結(jié) 語

為了解決PCD算法在實(shí)際應(yīng)用中的計(jì)算成本問題,本文定義一種新的音頻信號(hào)時(shí)域處理框架,將每個(gè)被分割的音頻幀作為一個(gè)列向量生成信號(hào)矩陣。在聯(lián)合稀疏表示(JSR)和同時(shí)稀疏近似(SSA)的基礎(chǔ)上,利用新的處理框架,提出以矩陣為操作對(duì)象的Joint-PCD算法。仿真的結(jié)果表明,Joint-PCD算法不僅能夠提供與PCD算法幾乎完全相同的降噪性能,而且Joint-PCD算法提供了更高的計(jì)算效率,節(jié)約了算法的運(yùn)行時(shí)間成本，從而加快了算法的收斂速度。更為重要的是,與PCD算法相比較,Joint-PCD算法在大規(guī)模數(shù)據(jù)處理和實(shí)時(shí)信號(hào)處理應(yīng)用方面具有巨大的潛在優(yōu)勢(shì)。