本發(fā)明涉及一種數(shù)字通信技術(shù)，特別是一種基于壓縮感知的直擴(kuò)MSK信號(hào)二維聯(lián)合捕獲方法。
背景技術(shù)：
：目前，擴(kuò)頻系統(tǒng)多采用BPSK、QPSK調(diào)制方式，但這兩種調(diào)制方式無(wú)法適用于存在嚴(yán)重非線性失真、多普勒頻移與多徑衰落的應(yīng)用領(lǐng)域。直擴(kuò)MSK通信系統(tǒng)同時(shí)具有擴(kuò)頻系統(tǒng)的低截獲性、多用戶隨機(jī)選址能力、抗干擾性能強(qiáng)等優(yōu)點(diǎn)和MSK信號(hào)的包絡(luò)恒定、頻譜利用率高、能量集中、旁瓣衰減快、帶外輻射功率低、對(duì)非線性失真不敏感等優(yōu)點(diǎn)，在戰(zhàn)術(shù)數(shù)據(jù)鏈、民用航空地空數(shù)據(jù)鏈、導(dǎo)彈制導(dǎo)指令傳輸、衛(wèi)星通信等領(lǐng)域得到了廣泛應(yīng)用。所以，直擴(kuò)MSK信號(hào)在直擴(kuò)BPSK/QPSK信號(hào)無(wú)法適用的領(lǐng)域仍然具有很好的應(yīng)用前景。從接收直擴(kuò)MSK信號(hào)中恢復(fù)發(fā)送數(shù)據(jù)符號(hào)的前提條件是信號(hào)的同步，包括擴(kuò)頻碼同步、載波的頻率和相位同步。信號(hào)同步分為粗同步和精同步，本發(fā)明側(cè)重于信號(hào)粗同步，即偽碼和多普勒頻偏捕獲。經(jīng)典的直擴(kuò)MSK信號(hào)捕獲方法有兩類(lèi)，一類(lèi)是基于滑動(dòng)相關(guān)的碼捕獲，這種方法實(shí)現(xiàn)簡(jiǎn)單但捕獲時(shí)間很長(zhǎng)；另一類(lèi)是基于匹配濾波的碼捕獲。這兩種方法的捕獲相關(guān)峰極易受多普勒頻偏的影響，顯然不適用于高動(dòng)態(tài)下的捕獲。針對(duì)存在大多普勒頻偏的高動(dòng)態(tài)環(huán)境，有很多學(xué)者提出了解決方法，如延遲相關(guān)法估計(jì)偽碼相差與多普勒聯(lián)合估計(jì)、部分匹配濾波與FFT相結(jié)合(PMF-FFT)等。然而這些方法在多普勒估計(jì)范圍與估計(jì)誤差之間大多較難取得平衡。所以，高動(dòng)態(tài)、低信噪比下直擴(kuò)MSK通信系統(tǒng)的信號(hào)同步成為技術(shù)難點(diǎn)。壓縮感知是一種新型的采樣理論，由E.J.Candes、J.Romberg、T.Tao和D.L.Donoho等科學(xué)家于2004年提出，其利用信號(hào)的稀疏特性，能夠在低于Nyquist采樣率的條件下通過(guò)隨機(jī)采樣獲得離散的信號(hào)值，接著采用非線性重建算法完整的重構(gòu)期望信號(hào)。本發(fā)明將壓縮感知理論運(yùn)用于直擴(kuò)MSK信號(hào)的快速捕獲中，設(shè)計(jì)出一種基于壓縮感知的直擴(kuò)MSK信號(hào)二維聯(lián)合捕獲方法，可以在高動(dòng)態(tài)、低信噪比環(huán)境下準(zhǔn)確、快速地完成擴(kuò)頻信號(hào)的偽碼相位和多普勒頻偏聯(lián)合捕獲。技術(shù)實(shí)現(xiàn)要素：本發(fā)明的目的在于提供一種基于壓縮感知的直擴(kuò)MSK信號(hào)二維聯(lián)合捕獲方法，實(shí)現(xiàn)高動(dòng)態(tài)、低信噪比下直擴(kuò)MSK信號(hào)的快速捕獲，同時(shí)完成偽碼相位和多普勒頻偏的聯(lián)合估計(jì)。包括以下步驟：對(duì)中頻直擴(kuò)MSK信號(hào)進(jìn)行處理獲得近似直擴(kuò)BPSK信號(hào)；構(gòu)建壓縮感知方程；取近似直擴(kuò)BPSK信號(hào)的一個(gè)擴(kuò)頻碼周期的樣點(diǎn)序列作為觀測(cè)序列，與根據(jù)本地偽碼構(gòu)造的第一冗余字典用正交匹配追蹤算法，捕獲碼相位誤差與多普勒頻偏的粗估計(jì)范圍；縮小多普勒估計(jì)范圍，觀測(cè)序列與構(gòu)造好的第二冗余字典做正交匹配追蹤重構(gòu)，得到多普勒頻偏精確的估計(jì)值，進(jìn)而完成碼捕獲環(huán)節(jié)的工作。本發(fā)明與現(xiàn)有技術(shù)相比，具有以下優(yōu)點(diǎn)：(1)冗余字典的重構(gòu)只需要一個(gè)擴(kuò)頻碼周期的采樣數(shù)據(jù)，所以該算法的捕獲時(shí)間短，大大提高了捕獲速度；(2)重構(gòu)結(jié)果即可得到所需捕獲的碼相位誤差和多普勒頻偏，不需要再用匹配濾波等方式對(duì)信號(hào)再做處理，在一定程度上減少了計(jì)算量；(3)用二次捕獲的方式，進(jìn)一步提高了捕獲精度，同時(shí)也擴(kuò)大了多普勒搜索范圍，適用于高動(dòng)態(tài)環(huán)境；(4)正交匹配追蹤算法重構(gòu)過(guò)程中引入了擴(kuò)頻增益，且后續(xù)操作無(wú)信噪比損失，可用于低信噪比的環(huán)境；(5)對(duì)直擴(kuò)MSK信號(hào)進(jìn)行中頻匹配處理，將其轉(zhuǎn)換為直擴(kuò)BPSK信號(hào)，解決了MSK信號(hào)形式復(fù)雜的問(wèn)題，可將適用于直擴(kuò)BPSK信號(hào)的快速捕獲方法應(yīng)用于直擴(kuò)MSK信號(hào)。下面結(jié)合說(shuō)明書(shū)附圖對(duì)本發(fā)明作進(jìn)一步描述。附圖說(shuō)明圖1是本發(fā)明方法流程圖。圖2是本發(fā)明實(shí)施例的正交匹配追蹤算法重構(gòu)算法流程圖。圖3是本發(fā)明實(shí)施例的檢測(cè)概率隨輸入信噪比的變化曲線示意圖。具體實(shí)施方式結(jié)合圖1，本發(fā)明的技術(shù)方案為：對(duì)中頻直擴(kuò)MSK信號(hào)進(jìn)行處理獲得近似直擴(kuò)BPSK信號(hào)；構(gòu)建壓縮感知方程；取近似直擴(kuò)BPSK信號(hào)的一個(gè)擴(kuò)頻碼周期的樣點(diǎn)序列作為觀測(cè)序列，與根據(jù)本地偽碼構(gòu)造的第一冗余字典用正交匹配追蹤算法，捕獲碼相位誤差與多普勒頻偏的粗估計(jì)范圍；縮小多普勒估計(jì)范圍，觀測(cè)序列與構(gòu)造好的第二冗余字典做正交匹配追蹤重構(gòu)，得到多普勒頻偏精確的估計(jì)值，進(jìn)而完成碼捕獲環(huán)節(jié)的工作。當(dāng)多普勒頻偏取值范圍很大時(shí)，冗余字典A會(huì)構(gòu)造得很大，相應(yīng)的算法復(fù)雜度也很大，所以本發(fā)明采用二次捕獲的思想對(duì)算法進(jìn)行改進(jìn)。即一次粗捕獲確定碼相位差、多普勒頻率大致范圍，第二次在小的搜索區(qū)間搜索獲得精確的估計(jì)值。發(fā)送信號(hào)采用“導(dǎo)頻碼+發(fā)送數(shù)據(jù)”的形式，在發(fā)送數(shù)據(jù)的前面插入M個(gè)全1數(shù)據(jù)符號(hào)，用以接收機(jī)端的擴(kuò)頻碼和載波同步。導(dǎo)頻碼將多普勒頻偏、碼相位等補(bǔ)償過(guò)后，數(shù)據(jù)信號(hào)部分才可進(jìn)行精確的解擴(kuò)解調(diào)。直擴(kuò)MSK信號(hào)可采用串行方式產(chǎn)生，即擴(kuò)頻后的信號(hào)與載波cos(2πf1t)進(jìn)行BPSK調(diào)制得到直擴(kuò)BPSK信號(hào)，再經(jīng)轉(zhuǎn)換濾波器g(t)即可產(chǎn)生直擴(kuò)MSK信號(hào)。轉(zhuǎn)換濾波器的沖激響應(yīng)為g(t)=sin2πf2t0≤t≤T0others---(1)]]>其中，fc為載波頻率，T為擴(kuò)頻碼碼片寬度。針對(duì)直擴(kuò)MSK信號(hào)的這種串行產(chǎn)生方式，在接收端采用與轉(zhuǎn)換濾波器g(t)相匹配的濾波器h(t)進(jìn)行接收，對(duì)應(yīng)的接收轉(zhuǎn)換濾波器的沖激響應(yīng)為h(t)=g(T-t)=-sin2πf2t0≤t≤T0others---(2)]]>轉(zhuǎn)換濾波器的輸出信號(hào)再經(jīng)f1下變頻、低通濾波，再進(jìn)行K倍抽取，得到近似直擴(kuò)BPSK基帶信號(hào)。此抽取倍數(shù)與過(guò)采樣倍數(shù)一致，此時(shí)輸出基帶信號(hào)可表示為其中，τ為碼相位誤差；為初始相位，為了分析的方便，設(shè)為0，Rc為擴(kuò)頻碼速率，fd為多普勒頻偏，ωn為復(fù)高斯白噪聲，d為傳輸符號(hào)，c為擴(kuò)頻碼，n為樣點(diǎn)序號(hào)由于導(dǎo)頻碼部分?jǐn)?shù)據(jù)為1，式(3)可改寫(xiě)為r(n)=c(n-τ)ej2πfdnRc+ωn---(4)]]>令則接收信號(hào)的一般表示形式為r‾(n,τm,υm)=c(n-τm)ej2πυmn+ωn---(5)]]>其中，上標(biāo)m表示時(shí)延τ及多普勒υ為某一定值。則用矩陣表示接收信號(hào)為r(τm,υm)=[r‾(1,τm,υm)r‾(2,τm,υm)...r‾(MN,τm,υm)]T---(6)]]>其中，N為擴(kuò)頻碼周期。首先確定模型的壓縮感知(CS)方程：y＝A·s(7)其中，y為取自接收信號(hào)r的長(zhǎng)度為一個(gè)擴(kuò)頻碼周期的觀測(cè)值；A為冗余字典；s為信號(hào)在冗余字典A上的映射?；趬嚎s感知的碼捕獲算法最重要的就是冗余字典A的構(gòu)造。以時(shí)延τ、多普勒υ、信號(hào)強(qiáng)度r建立冗余坐標(biāo)系。將時(shí)延-多普勒域劃分成Nτ×Nf個(gè)區(qū)域。p，q的取值分別代表時(shí)延τ與多普勒值υ，則p＝1,2,…,Nτ，q＝0,1,…,Nf。對(duì)于第一冗余字典A1，其粗捕獲的多普勒分辨率為Δυ1，時(shí)延分辨率為Δτ；則多普勒估計(jì)范圍為[0,Δυ1Nf]，時(shí)延估計(jì)范圍為[0,ΔτNτ]。此時(shí)第一冗余字典A1可表示為：A1(n,q,p)=c(n-Δτ·p)ej2πΔυ1·qn=1,2,...,N;q=1,2,...,Nf1;p=1,2,...,Nτ---(8)]]>A1可以看成由Nτ個(gè)N×Nf二維數(shù)組組成的三維矩陣。A1中的每N×Nf二維數(shù)組可以看成時(shí)延固定的，而每列的多普勒頻偏是固定的。第一冗余字典A1上的映射s為K稀疏，此處K＝1：s(q,p)=1(τm,υm)=(τp,υq)p=1,2,...,Nτq=1,2,...,Nf0others---(9)]]>然后采用正交匹配追蹤(OMP)重構(gòu)算法，重構(gòu)出信號(hào)r在變換域的稀疏表示中不為0的索引即為A中p和q的估計(jì)值，即可求得fd與τ的估計(jì)值：τ=p^·Δτ---(10)]]>fd=q^×Δυ1·Rc---(11)]]>粗捕獲只能將多普勒頻偏確定在某個(gè)范圍：[(q^-0.5)×Δυ,(q^+0.5)×Δυ]---(12)]]>精捕獲的多普勒捕獲范圍由粗捕獲的捕獲結(jié)果確定。滿足條件：Δυ1＝Δυ2×Nf2(13)則精捕獲的冗余字典為：A2(n,m)=c(n-Δτ·p^)ej2π[Δυ1·(q^-0.5)+Δυ2·m]n=1,2,...,N;m=1,2,...,Nf2---(14)]]>與粗捕獲過(guò)程相同，冗余字典A2上的映射s為K稀疏，此處K＝1：s(i)=1υm=υi,i=1,2,...,Nf20others---(15)]]>其中，A2為精捕獲的冗余字典；Δυ2為精捕獲的多普勒分辨率；Nf2為精捕獲的多普勒取值個(gè)數(shù)。在粗捕獲的基礎(chǔ)上，精捕獲過(guò)程的估計(jì)值為則本發(fā)明最終的多普勒估計(jì)值為：fd=(q^-0.5)·Δυ1+i^·Δυ2---(16)]]>OMP算法步驟如下：Step1：初始化r0＝y(tǒng)，t＝1；Step2：找到索引λt，使得：Step3：令Λt＝Λt-1∪{λt}，Step4：求y＝Atst得最小二乘解；Step5：更新殘差Step6：t＝t+1，如果t≤K，則返回Step2，否則停止迭代進(jìn)入Step7；Step7：重構(gòu)所得在Λt處有非零值，其值分別為最后一次迭代所得rt為殘差，t為迭代次數(shù)，K為迭代總數(shù)，Λt為t次迭代的索引集合，λt為t次迭代的索引，At為按索引Λt選出的冗余字典的列集合，αn為At的第n列，為重構(gòu)出信號(hào)在變換域的稀疏表示。實(shí)施例系統(tǒng)采樣頻率fs＝245.52MHz，中頻頻率fc＝76.725MHz，過(guò)采樣倍數(shù)K＝12，擴(kuò)頻碼速率為Rc＝20.46Mchip/s，數(shù)據(jù)速率為20kbps，擴(kuò)頻碼采用Gold序列，碼長(zhǎng)N＝1023，導(dǎo)頻碼符號(hào)個(gè)數(shù)M＝20。接收中頻直擴(kuò)MSK信號(hào)經(jīng)沖激響應(yīng)為h(t)的接收轉(zhuǎn)換濾波器，得到近似直擴(kuò)BPSK信號(hào)。接收轉(zhuǎn)換濾波器的沖激響應(yīng)為h(t)=-sin2πf2t0≤t≤T0others---(18)]]>其中，為擴(kuò)頻碼周期。轉(zhuǎn)換濾波器的頻率響應(yīng)為D(f)=-sinπTc(f-f2)πTc(f-f2)---(19)]]>本發(fā)明濾波器的設(shè)計(jì)采用凸優(yōu)化技術(shù)，首先需要將濾波器設(shè)計(jì)問(wèn)題轉(zhuǎn)化為凸優(yōu)化問(wèn)題，可以將轉(zhuǎn)換濾波器的切比雪夫逼近建立為凸優(yōu)化模型：minimizesupω∈[0,π]|H(ω)-D(ω)|---(20)]]>其中，sup為下確界；ω＝2πf為角頻率；D(ω)為給定的頻率響應(yīng)函數(shù)；H(ω)為所設(shè)計(jì)的濾波器的頻率響應(yīng)，h(n)為濾波器系數(shù)，N0為濾波器階數(shù)。實(shí)際應(yīng)用的濾波器系數(shù)h(n)通過(guò)Matlab軟件的cvx工具箱求得。經(jīng)接收轉(zhuǎn)換濾波器后的信號(hào)再經(jīng)f1下變頻，其中下變頻后低通濾波，再K倍抽取，此抽取倍數(shù)與過(guò)采樣倍數(shù)一致。由于導(dǎo)頻碼部分?jǐn)?shù)據(jù)為1，此時(shí)得到基帶信號(hào)為其中，τ為碼相位誤差；為初始相位；Rc為擴(kuò)頻碼速率；fd為多普勒頻偏；ωn為高斯復(fù)白噪聲；為了方便分析，設(shè)定初始相位令則接收信號(hào)的一般表示形式為：r‾(n,τm,υm)=c(n-τm)ej2πυmn+ωn---(22)]]>其中，上標(biāo)m表示時(shí)延τ及多普勒υ為某一定值。則用矩陣表示接收信號(hào)為：r(τm,υm)=[r‾(1,τm,υm)r‾(2,τm,υm)...r‾(M,τm,υm)]M×1T---(23)]]>首先確定模型的CS方程為：y＝A·s(24)其中，y為取自接收信號(hào)r的長(zhǎng)度為一個(gè)擴(kuò)頻碼周期的觀測(cè)值；A為冗余字典；s為信號(hào)在冗余字典A上的映射?；趬嚎s感知的碼捕獲算法最重要的就是冗余字典A的構(gòu)造。本發(fā)明采用二次捕獲的方法，即一次粗捕獲確定碼相位誤差和多普勒頻偏范圍，第二次精捕獲確定多普勒頻偏的精確估計(jì)值。粗捕獲過(guò)程將時(shí)延-多普勒域劃分成Nω＝Nτ·Nf個(gè)區(qū)域。p，q的取值分別代表碼相位誤差τ與多普勒υ，則p＝0,1,…,Nτ，其中多普勒分辨率為Δυ1，時(shí)延分辨率為Δτ；則多普勒估計(jì)范圍為時(shí)延估計(jì)范圍為[0,ΔτNτ]。則粗捕獲過(guò)程的冗余字典為：A1(n,q,p)=c(n-Δτ·p)ej2πΔυ1·qn=1,2,...,N;q=1,2,...,Nf1;p=1,2,...,Nτ---(25)]]>A1可以看成由Nτ個(gè)二維數(shù)組組成的三維矩陣。A1中的每二維數(shù)組可以看成時(shí)延固定的，而每列的多普勒頻偏是固定的。冗余字典A1上的映射s為K稀疏，此處K＝1：s(q,p)=1(τm,υm)=(τp,υq)p=1,2,...,Nτq=1,2,...,Nf10others---(26)]]>最后采用正交匹配追蹤(OMP)重構(gòu)算法，重構(gòu)出信號(hào)r在變換域的稀疏表示中不為0的索引即為A1中p和q的估計(jì)值。如圖2所示，OMP算法步驟如下：Step1：初始化r0＝y(tǒng)，t＝1；Step2：找到索引λt，使得：Step3：令Λt＝Λt-1∪{λt}，Step4：求y＝Atst得最小二乘解；Step5：更新殘差Step6：t＝t+1，如果t≤K，則返回第(2)步，否則停止迭代進(jìn)入第(7)步；Step7：重構(gòu)所得在Λt處有非零值，其值分別為最后一次迭代所得OMP算法估計(jì)出粗捕獲過(guò)程的碼相位誤差和多普勒的索引估計(jì)值后，即可求得時(shí)延估計(jì)值為：τ=p^·Δτ---(27)]]>粗捕獲確定的多普勒頻偏的范圍為：[(q^-0.5)×Δυ1,(q^+0.5)×Δυ1]---(28)]]>精捕獲過(guò)程在粗捕獲的基礎(chǔ)上，將多普勒進(jìn)一步精確地估計(jì)出來(lái)，則精捕獲的多普勒捕獲范圍由粗捕獲的捕獲結(jié)果確定。捕獲的分辨率和范圍滿足條件：Δυ1＝Δυ2×Nf2(29)則精捕獲的冗余字典為：A2(n,m)=c(n-Δτ·p^)ej2π[Δυ1·(q^-0.5)+Δυ2·m]n=1,2,...,N;m=1,2,...,Nf2---(30)]]>與粗捕獲過(guò)程相同，冗余字典A2上的映射s為K稀疏，此處K＝1：s(i)=1υm=υi,i=1,2,...,Nf20others---(31)]]>最后仍然采用正交匹配追蹤(OMP)重構(gòu)算法，求出多普勒估計(jì)值在粗捕獲的基礎(chǔ)上，精捕獲的過(guò)程可以求出本發(fā)明的多普勒頻偏估計(jì)值為：fd^=(q^-0.5)·Δυ1+m^·Δυ2---(32)]]>本發(fā)明的實(shí)例的參數(shù)設(shè)置為：Δτ＝1，Nτ＝1023，Δυ1＝1500，Δυ2＝10，計(jì)算式(28)可得本發(fā)明多普勒估計(jì)范圍：[0,450]kHz；碼相位誤差估計(jì)范圍：[0,1023]。圖3為輸入信噪比SNR＝[-25dB,-5dB]，τ＝512.5chip，fd＝116kHz時(shí)，得到的檢測(cè)概率和虛警概率隨信噪比的變化曲線。由圖可以看出當(dāng)信噪比達(dá)到-20dB時(shí)，檢測(cè)概率達(dá)到0.9，可見(jiàn)本發(fā)明設(shè)計(jì)的偽碼-多普勒聯(lián)合捕獲方法能夠在高動(dòng)態(tài)環(huán)境下對(duì)偽碼相位和載波頻偏進(jìn)行準(zhǔn)確的捕獲。當(dāng)前第1頁(yè)1 2 3