專利名稱：基于低密度奇偶校驗(yàn)譯碼軟判決信息的迭代定時(shí)同步方法
技術(shù)領(lǐng)域：
本發(fā)明屬于數(shù)字通信技術(shù)領(lǐng)域，特別涉及到較低信噪比(SNR)情況下的基于低密度奇 偶校驗(yàn)(LDPC)譯碼軟判決信息(SDM)的迭代定時(shí)同步方法與裝置。
背景技術(shù)：
在數(shù)字通信領(lǐng)域中，接收機(jī)因系統(tǒng)熱噪聲、傳輸時(shí)延、收發(fā)端頻率源偏差、通信接收 機(jī)晶振頻率漂移、載波同步后的殘留偏差等因素影響，不能實(shí)現(xiàn)接收碼元符號(hào)的最佳定時(shí) 采樣判決。這將導(dǎo)致系統(tǒng)通信性能降低甚至失效。針對(duì)該問題，接收機(jī)通?？刹捎枚〞r(shí)同 步方法來(lái)解決。目前，實(shí)用的通信接收機(jī)定時(shí)同步方法主要有基于最小均方誤差的有數(shù)據(jù)輔助定時(shí)同 步，利用信號(hào)自身統(tǒng)計(jì)特性的無(wú)數(shù)據(jù)輔助定時(shí)同步及結(jié)合信道碼的迭代定時(shí)同步等三類方 '法。其中，采用有數(shù)據(jù)輔助定時(shí)同步已有實(shí)用方法，如Miller&Miiller定時(shí)誤差檢測(cè)方法， 見文獻(xiàn)K. Miller and M. Miiller， "Timing Recovery for Digital Synchronous Data Receivers", IEEE Trans, on Com., 1976, 24: 516-531。但該方法有以下不足使用訓(xùn)練數(shù)據(jù)多占用傳輸帶 寬，不適于單點(diǎn)對(duì)多點(diǎn)通信的廣播通信等缺陷。所以，也有人提出利用信號(hào)自身統(tǒng)計(jì)特性 的無(wú)數(shù)據(jù)輔助定時(shí)同步方法，如Gardner早遲門定時(shí)同步方法，見文獻(xiàn)F. M. Gardner, "A BPSK/QPSK Timing Error Detector for Sampled Receivers", IEEE Trans, on Com. 1986, 34(5): 423-429。但該方法在較低SNR等情況下性能較差。另外，以上兩類傳統(tǒng)定時(shí)同步方法的綜 合性描述著作，可見文獻(xiàn)Umberto Mengali et.al., "Synchronization techniques for digital receivers", hew York: Plenum Press, 1997。這兩類方法都忽視了數(shù)字通信系統(tǒng)一般都采用信 道碼的事實(shí)。同時(shí)，因估計(jì)誤差受參數(shù)估計(jì)理論中Cramer-Rao下界的約束，在較低SNR情 ^L下定時(shí)估計(jì)的可靠性都較差。而且，這些方法捕獲過程都較慢，影響了它們?cè)谳^低SNR 情況下的實(shí)用性。第三類結(jié)合信道碼的迭代定時(shí)同步方法主要是采用信道譯碼的中間結(jié)果 來(lái)涵行定時(shí)誤差估計(jì)的定時(shí)同步方法。根據(jù)最新研究，可由所采用的信道編碼的不同，分 為格型譯碼輔助迭代定時(shí)同步和分組譯碼輔助迭代定時(shí)同步兩類方法格型譯碼輔助迭代，al同歩方法采用搜索最可能滿足該信道編碼格型譯碼約束條件的幸存路徑，來(lái)判斷最佳l]i時(shí)同步采樣點(diǎn)。該方法可見文獻(xiàn)John R. Barry, et.al. "Iterative Timing Recovery",正EE Sigjal Processing Magazine, 2004, 89-102。該方法主要存在計(jì)算復(fù)雜度過高和存儲(chǔ)量過大等 l奇^^而分組譯碼輔助迭代定時(shí)同步方法主要利用分組碼(如LDPC碼)的迭代譯碼過程中軟 判決譯碼中間結(jié)果滿足LDPC校驗(yàn)方程約束數(shù)量的多少程度來(lái)判決最佳采樣點(diǎn)，即對(duì)接收 信號(hào)高倍插值采樣，并將譯碼結(jié)果滿足校驗(yàn)方程約束數(shù)量最多的采樣點(diǎn)作為最佳定時(shí)采樣點(diǎn)，來(lái)實(shí)現(xiàn)定時(shí)同步。該方法可見文獻(xiàn)Dong-U Lee, et.al., "Joint LDPC Decoding and Timing Recovery Using Code Constraint Feedback", IEEE Com. Letters, 2006, 10(3): 189-191 。該方法
的復(fù)雜性有所降低，但其迭代計(jì)算次數(shù)還是較大，性能較差，限制了它的實(shí)用性?？傊?， 在較低SNR情況下，現(xiàn)有迭代定時(shí)同步方法性能較好。但因其算法復(fù)雜及計(jì)算量大等缺陷， 目前尚未見有效的實(shí)用方法。
另外，在很多情況下，定時(shí)同步處理至少以兩種階段或操作模式進(jìn)行工作:(a)捕獲階 段；(b)跟蹤階段。在定時(shí)同步的捕獲階段，定時(shí)同步模塊采用搜索、最大似然估計(jì)等方法 來(lái)粗估計(jì)定時(shí)誤差并進(jìn)行初步矯正。在定時(shí)捕獲之后，定時(shí)同步處理進(jìn)入跟蹤階段，用最 大似然反饋環(huán)等方法使定時(shí)采樣點(diǎn)穩(wěn)定的保持于準(zhǔn)確定時(shí)采樣點(diǎn)的某個(gè)小范圍內(nèi)，以減少 定時(shí)同步誤差。
在本發(fā)明中，需要用到插值、匹配濾波、相關(guān).、LDPC迭代譯碼、單純形極小值搜索、 Miller&Miiller定時(shí)誤差檢測(cè)和數(shù)字環(huán)路濾波等現(xiàn)有技術(shù)，這些技術(shù)分別描述如下
插值方法
插值方法是數(shù)值計(jì)算中常見的方法，它能根據(jù)己知若干樣本點(diǎn)及其函數(shù)值，'來(lái)估計(jì)未 知樣本點(diǎn)的函數(shù)值。插值方法主要包括線性分段插值，Lagrange插值，牛頓插值，樣條插 值等方法。其中，線性分段插值方法性能較差，但計(jì)算量少，適用于硬件資源有限的應(yīng)用 場(chǎng)合。其他插值方法性能較好，但計(jì)算量較大，更適用于理論分析的應(yīng)用場(chǎng)合。
匹配濾波方法
在數(shù)字通信中，為了使采樣信號(hào)滿足奈奎斯特采樣定理而不產(chǎn)生碼間干擾，在通信發(fā) 送,和接收端需要分別采用對(duì)稱的均方根升余弦(SRRC)成型及匹配濾波處理。實(shí)現(xiàn)該方法 的ri波器參數(shù)主要是滾降系數(shù)。另外，在具體實(shí)施過程中，可通過數(shù)字化處理，均勻取以 中心為基準(zhǔn)對(duì)稱的若干離散點(diǎn)，并忽略離中心點(diǎn)較遠(yuǎn)處函數(shù)值較小點(diǎn)，作為該濾波器的抽 頭系數(shù)。
相關(guān)主、要是通過特定的運(yùn)算來(lái)求解未知數(shù)據(jù)與已知數(shù)據(jù)的相似程度，如兩者越相似， f尋到的相關(guān)結(jié)果越大。相關(guān)的特定計(jì)算過程如下如已知數(shù)據(jù)長(zhǎng)為L(zhǎng)(L為自然數(shù))，每次逐 赤號(hào)的從未知的數(shù)據(jù)中取連續(xù)長(zhǎng)度為L(zhǎng)的數(shù)據(jù)，且按順序?qū)⒃摂?shù)據(jù)中每個(gè)數(shù)據(jù)分別與已知 余#中對(duì)應(yīng)序號(hào)的數(shù)據(jù)相乘，然后將這L個(gè)相乘結(jié)果相加，得到相關(guān)結(jié)果。
LDPC的迭代譯碼方法
;LDPC的迭代譯碼方法可用迭代的方式以低計(jì)算復(fù)雜度實(shí)現(xiàn)LDPC碼的譯碼。該方法 ^i體實(shí)現(xiàn)最常用的是和積譯碼算法。它的詳細(xì)描述可見文獻(xiàn)MacKay D. J., "Good Error-Correcting Codes Based on Very Sparse Matrices", IEEE Trans. On Inform. Theory. 1999， 45(i5): 399-431。該算法包括以下內(nèi)容
LDPC譯碼算法先定義變量如下"c,)為變量節(jié)點(diǎn)f的先驗(yàn)信息，即對(duì)數(shù)似然比； 丄( W)表示第Jt次迭代中，從校驗(yàn)節(jié)點(diǎn)y到變量節(jié)點(diǎn)/的外信息；丄(^W)表示第ik次迭代中，從變量節(jié)點(diǎn)/到校驗(yàn)節(jié)點(diǎn)y的外信息；丄(2/")為第/fc次迭代后變量節(jié)點(diǎn)/的后驗(yàn)信息； C(/)表示和變量節(jié)點(diǎn)/相連的校驗(yàn)節(jié)點(diǎn)的集合；及(y')表示和校驗(yàn)節(jié)點(diǎn)y相連的變量節(jié)點(diǎn)的崔A 朱O o具體的LDPC譯碼算法可以分為以下三個(gè)步驟 首先，LDPC譯碼的初始化A^/o))"(c,), = 0 (2)"^")三"C三0，當(dāng)變量節(jié)點(diǎn)/與校驗(yàn)節(jié)點(diǎn)y不相連 (3) 其中，ij為自然數(shù)，表示節(jié)點(diǎn)序號(hào)；；c,為發(fā)送比特，^為接收符號(hào)；如接收系統(tǒng)采用BPSK調(diào)制，其映射為"o"4"+r，'T'i~>"-r;其他調(diào)制方式也可通過解調(diào)、解映射及并串變換的方式轉(zhuǎn)成這種映射結(jié)構(gòu)；信道為加性高斯白噪聲(AWGN)信道，AWGN ^差為CT2; log(q)為自然對(duì)數(shù)函數(shù)；戶(1,=+1|;/,)為當(dāng)前收到信號(hào)乂，則《=+1的概率，尸(、=+1|>;,) 的定義也類似。其次，LDPC碼的迭代譯碼迭代譯碼可包括兩個(gè)部分，變量節(jié)點(diǎn)的更新計(jì)算和校驗(yàn) 節(jié)點(diǎn)的更新計(jì)算，在本發(fā)明中設(shè)定固定的迭代次數(shù)K次(K為自然數(shù)，可選為20~30)，然 后判決輸出。變量節(jié)點(diǎn)的運(yùn)算(對(duì)所有的變量節(jié)點(diǎn)i):丄 ))=丄")+ S 〖4)其中，Z(q)為累加求和函數(shù)；其下標(biāo)y'ec(/)"'表示所求知的項(xiàng)的標(biāo)號(hào)y'為除去j外 的所有屬于集合cw的值。校驗(yàn)節(jié)點(diǎn)的運(yùn)算(對(duì)所有的校驗(yàn)節(jié)點(diǎn)j):丄(r/))-tan!Tr 、 J"! tanh("0/2)(5)其中，yt^l; tanh(Q)為雙曲正切函數(shù)；n 為對(duì)括號(hào)內(nèi)的項(xiàng)進(jìn)行相乘運(yùn)算的乘積函數(shù)； 為了簡(jiǎn)化處理以方便數(shù)字集成電路實(shí)現(xiàn)，可分別定義以下公式(6)(7),-log411 (8) e —l其中，w'gw(q)為符號(hào)函數(shù)；iq為絕對(duì)值函數(shù)。則以上公式(5)也可以表示為最后，LDPC譯碼后驗(yàn)信息的計(jì)算和符號(hào)硬判決計(jì)算分別如下:<formula>formula see original document page 9</formula>其中，公式(ll)中得到的A即為對(duì)應(yīng)序號(hào)i的最終譯碼后的碼字。另外，如設(shè)該LDPC 碼的校驗(yàn)矩陣為二維RxN的矩陣(N,R為自然數(shù)，分別表示為L(zhǎng)DPC碼的碼長(zhǎng)和校驗(yàn)方程的 個(gè)數(shù)，且有IKN)，這樣就有N個(gè)硬判決譯碼結(jié)果。最后，可定義該N個(gè)硬判決譯碼結(jié)果中 前(N-R)個(gè)硬判決譯碼結(jié)果(這部分對(duì)應(yīng)于編碼前的信息比特)和全部N個(gè)硬判決譯碼結(jié)果 (這部分包括編碼前的信息碼元和相應(yīng)的LDPC校驗(yàn)碼元)分別作為L(zhǎng)DPC信息碼元和LDPC 碼字。-單純形極小值搜索方法單純形極小值搜索方法是Nelder和Mead提出并經(jīng)后人改進(jìn)的極小值搜索算法。該方法 可用于求解多維函數(shù)的極值，且其具有無(wú)需求導(dǎo)運(yùn)算，算法簡(jiǎn)單魯棒，適合變?cè)獢?shù)不多及 函數(shù)值不連續(xù)的方程求極值等優(yōu)勢(shì)，適用于無(wú)導(dǎo)數(shù)可用時(shí)極小值搜索算法的求解。該方法 的具體描述及基本原理可見文獻(xiàn)Jeffrey C. L.， et.al., "Convergence Properties of the Nelder-Mead Simplex Method in Low Dimensions", SIAM Journal of Optimization, 1998， 9(1):112-147。該方法主要通過參數(shù)初始化、反射、擴(kuò)展、壓縮等步驟，每次迭代運(yùn)算時(shí)， 將初始輸入的M+1個(gè)數(shù)據(jù)點(diǎn)及函數(shù)值的一組參數(shù)進(jìn)行更新，直到更新后的數(shù)據(jù)點(diǎn)及函數(shù)值 滿足設(shè)定的誤差條件或迭代次數(shù)達(dá)到設(shè)定的次數(shù)。最后，將最后更新的數(shù)據(jù)點(diǎn)作為該方法 所求的極小值。Miller&Mimer定時(shí)誤差檢測(cè)方法Miller&Mtmer定時(shí)誤差檢測(cè)方法主要通過數(shù)據(jù)輔助來(lái)獲得定時(shí)誤差的最大似然估計(jì)， 并用該定時(shí)誤差的估計(jì)來(lái)實(shí)現(xiàn)定時(shí)同步的方法。該方法的具體描述及基本原理可見文獻(xiàn) K. Mueller and M. Mtiller, "Timing Recovery for Digital Synchronous Data Receivers", IEEE Trans, on COM.， 1976,24:516-531。該方法得到的定時(shí)誤差與接收信號(hào)及期望信號(hào)之間的關(guān) 系如下<formula>formula see original document page 9</formula>其中，k為自然數(shù)，表示采樣序號(hào)或迭代計(jì)算的次數(shù)；^是第k次迭代計(jì)算時(shí)的誤差估計(jì)； ^和《分別是序號(hào)為k的接收信號(hào)值及其期望的準(zhǔn)確參考值。<formula>formula see original document page 9</formula>環(huán)路濾波方法環(huán)路濾波方法通常用于鎖相環(huán)結(jié)構(gòu)的電路中，其作用為濾除誤差信號(hào)中的高頻成分及 噪聲，并控制環(huán)路誤差的校正速度和精度。它可用一階數(shù)字濾波器實(shí)現(xiàn)誤差信號(hào)在進(jìn)入 該數(shù)字環(huán)路濾波器后，分為兩路 一路直接乘以系數(shù)G,;另外一路經(jīng)過積分器積分后，再 乘以系數(shù)G2。之后，將這兩路信號(hào)相加，作為數(shù)字環(huán)路濾波器的輸出。其中，G,和G2為實(shí) 數(shù)；為了使該環(huán)路濾波方法穩(wěn)定，需要合理選擇參數(shù)G,和G2，使該數(shù)字濾波器傳輸函數(shù)的極點(diǎn)在單位圓內(nèi)。 發(fā)明內(nèi)容本發(fā)明的目的是為了克服已有定時(shí)同步技術(shù)的不足，提出基于LDPC譯碼軟判決信息 (SDM)的迭代定時(shí)同步方法及其裝置。本發(fā)明的方法具有能在極低SNR情況下工作，提高 接收信號(hào)的誤碼性能，實(shí)現(xiàn)復(fù)雜度較低，計(jì)算量較少等優(yōu)點(diǎn)。而且，本發(fā)明的裝置還能以 全并行或部分并行實(shí)現(xiàn)，以克服因硬件速度限制而不能實(shí)現(xiàn)高速同步的問題。所以，本發(fā) 明非常適用于較低SNR情況下采用大規(guī)模集成電路(VLSI)方式實(shí)現(xiàn)的定時(shí)同步。本發(fā)明的基于LDPC譯碼軟判決信息(SDM)的迭代定時(shí)同步方法，其特征在于，包括 以下步驟1) 將系統(tǒng)接收端收到的數(shù)據(jù)分別經(jīng)M倍符號(hào)頻率采樣、用插值方法進(jìn)行插值和均方 根升余弦(SRRC)匹配濾波得到濾波后數(shù)據(jù)；其中，M為自然數(shù)(可取為4等偶數(shù)，匹配濾 波所用的SRRC濾波器的滾降系數(shù)可根據(jù)工程實(shí)際取值為0.3-0.8);在首次執(zhí)行插值處理 時(shí)，不進(jìn)行插值(否則，用步驟4)反饋的插值位置進(jìn)行插值)；2) 對(duì)步驟l)所述濾波后數(shù)據(jù)分別進(jìn)行兩路輸出第1路數(shù)據(jù)進(jìn)行LDPC碼幀同步相關(guān) 檢測(cè)方法處理，得到LDPC碼數(shù)據(jù)幀起始位置序號(hào)和幀同步檢測(cè)成功信號(hào)；第2路數(shù)據(jù)則 以相關(guān)檢測(cè)方法得到的LDPC碼數(shù)據(jù)幀起始位置序號(hào)為濾波后新數(shù)據(jù)的起始序號(hào)，將濾波 后新數(shù)據(jù)進(jìn)行1:M的串并變換，得到M路數(shù)據(jù)，且每路各N+L個(gè)數(shù)據(jù)；其中，N、 L都 為自然數(shù)，分別表示LDPC碼的幀長(zhǎng)和LDPC碼同步頭的長(zhǎng)度；3) 在接收到所述幀同步檢測(cè)成功信號(hào)后，進(jìn)行迭代定時(shí)控制方法處理如為初次處理，則輸出定時(shí)捕獲控制信號(hào)，轉(zhuǎn)步驟4);否則轉(zhuǎn)步驟6);4) 在接收到所述定時(shí)捕獲控制信號(hào)后，對(duì)步驟2)所述M路數(shù)據(jù)用迭代定時(shí)捕獲方法 處理，得到插值位置和定時(shí)捕獲輸出數(shù)據(jù)，將得到的插值位置反饋到步驟1)，用于接收端 新接收數(shù)據(jù)的插值；5) 將所述定時(shí)捕獲輸出數(shù)據(jù)進(jìn)行K次LDPC迭代譯碼，得到N個(gè)LDPC碼字和(N-R) 個(gè)LDPC信息碼元；其中，K為自然數(shù)(可選為20 30); R為自然數(shù)，表示LDPC碼校驗(yàn)方 程的數(shù)量(即LDPC碼的校驗(yàn)矩陣為R行N列的矩陣)，且有IKN;6) 將所得的N個(gè)LDPC碼字進(jìn)行迭代定時(shí)控制方法處理，得到定時(shí)捕獲控制信號(hào)、 定時(shí)跟蹤控制信號(hào)或最終LDPC譯碼輸出控制信號(hào)；若得到定時(shí)捕獲控制信號(hào)，設(shè)置定時(shí)狀態(tài)為捕獲狀態(tài)，則轉(zhuǎn)步驟4);若得到定時(shí)跟蹤控制信號(hào)，設(shè)置定時(shí)狀態(tài)為跟蹤狀態(tài)，則轉(zhuǎn)步驟7);若得到最終LDPC譯碼輸出控制信號(hào)，將所述的LDPC信息碼元，作為系統(tǒng)最 終的LDPC譯碼輸出數(shù)據(jù)，并輸出，并轉(zhuǎn)步驟9);7) 將步驟2)所得的M路數(shù)據(jù)和N個(gè)LDPC碼字，進(jìn)行迭代定時(shí)跟蹤方法處理，得到 定時(shí)跟蹤輸出數(shù)據(jù)；8) 將步驟7)所述定時(shí)跟蹤輸出數(shù)據(jù)進(jìn)行K次LDPC迭代譯碼，得到N個(gè)LDPC碼字 和LDPC信息碼元，并轉(zhuǎn)(6);9) 接收端將接收到的新輸入數(shù)據(jù)，轉(zhuǎn)步驟1)，直到接收端接收數(shù)據(jù)完畢。 本發(fā)明采用上述方法的定時(shí)同步裝置，其特征在于，用于實(shí)現(xiàn)基于LDPC譯碼軟判決信息的迭代定時(shí)同步；該裝置具體包括以下模塊數(shù)據(jù)接收處理模塊；LDPC碼幀同步相 關(guān)檢測(cè)模塊；迭代定時(shí)捕獲模塊；迭代定時(shí)控制模塊；迭代定時(shí)跟蹤模塊；LDPC譯碼模 塊。該裝置的連接關(guān)系為數(shù)據(jù)接收處理模塊的采樣數(shù)據(jù)接收端作為整個(gè)系統(tǒng)的輸入端； 數(shù)據(jù)接收處理模塊的匹配濾波輸出端連接到LDPC碼幀同步相關(guān)檢測(cè)模塊的數(shù)據(jù)輸入端； 數(shù)據(jù)接收處理模塊中M路第一數(shù)據(jù)輸出端連接到迭代定時(shí)跟蹤模塊中的插值數(shù)桐輸入端； 數(shù)據(jù)接收處理模塊的M路第二數(shù)據(jù)輸出端連接到迭代定時(shí)捕獲模塊中的第一數(shù)據(jù)輸入端；數(shù)據(jù)接收處理模塊的M路數(shù)據(jù)第三組輸出端共(M+l)個(gè)輸出端連接到迭代定時(shí)捕獲模塊 的第二組數(shù)據(jù)輸入端的(M+l)個(gè)輸入端；LDPC碼幀同步相關(guān)檢測(cè)模塊的數(shù)據(jù)輸出端連接 到數(shù)據(jù)接收處理模塊的匹配濾波輸入端；LDPC碼幀同步相關(guān)檢測(cè)模塊的控制信號(hào)輸出端 連接到迭代定時(shí)控制模塊的控制輸入端；迭代定時(shí)捕獲模塊的數(shù)據(jù)輸出端連接到LDPC譯 碼模塊的第二數(shù)據(jù)輸入端；迭代定時(shí)捕獲模塊的插值輸出端連接到數(shù)據(jù)接收處理模塊的插 值輸入端；迭代定時(shí)跟蹤模塊的數(shù)據(jù)輸出端連接到LDPC譯碼模塊的第一數(shù)據(jù)輸入端； LDPC譯碼模塊的第一碼字輸出端連接到迭代定時(shí)跟蹤模塊的碼字?jǐn)?shù)據(jù)輸入端；LDPC譯 碼模塊的第二碼字輸出端連接到迭代定時(shí)控制模塊的數(shù)據(jù)輸入端；LDPC譯碼模塊的信息 碼元輸出端作為系統(tǒng)的輸出端；迭代定時(shí)控制模塊的定時(shí)捕獲控制輸出端連接到迭代定時(shí) 捕獲模塊的控制輸入端；迭代定時(shí)控制模塊的定時(shí)跟蹤控制輸出端連接到迭代定時(shí)跟蹤模 塊的控制輸入端；迭代定時(shí)控制模塊的最終輸出控制輸出端連接到LDPC譯碼模塊的控制 輸入端；迭代定時(shí)控制模塊的數(shù)據(jù)選擇控制輸出端連接到LDPC譯碼模塊的數(shù)據(jù)選擇控制 輸入端；其中，M為自然數(shù)，取4等偶數(shù)；'所述的LDPC碼幀同步相關(guān)檢測(cè)模塊，具體包括1:M路串并變換器；M個(gè)相關(guān)模塊; 最大值判斷器；比較器。該模塊的接口及內(nèi)部各器件的連接關(guān)系為1:M路串并變換器的 數(shù)據(jù)輸入端作為該模塊的數(shù)據(jù)輸入端；比較器的數(shù)據(jù)輸出端作為該模塊的數(shù)據(jù)輸出端；比 較器的控制信號(hào)輸出端作為該模塊的控制信號(hào)輸出端；1:M路串并變換器的M個(gè)輸出端分 別連接到M個(gè)相關(guān)模塊的輸入端；M個(gè)相關(guān)模塊的輸出端連接到最大值判斷器的M個(gè)輸 入端；最大值判斷器的輸出端連接到比較器的輸入端；其中，相關(guān)模塊用于輸入數(shù)據(jù)與本 地已知數(shù)據(jù)序列的相關(guān)計(jì)算，得到相關(guān)計(jì)算結(jié)果。所述的迭代定時(shí)捕獲模塊，具體包括第二組插值器，共包括(M+1)個(gè)插值器；(M+l) 個(gè)LDPC譯碼軟判決信息(SDM)計(jì)算器；二維聯(lián)合定時(shí)誤差和頻偏計(jì)算器；插值位置計(jì)算 器；第三插值器。該模塊的接口及內(nèi)部各器件的連接關(guān)系為該模塊有l(wèi)個(gè)控制該模塊工 作的控制輸入端；第三插值器的數(shù)據(jù)輸入端作為該迭代定時(shí)捕獲模塊的第一數(shù)據(jù)輸入端； 第二組插值器的(M+1)數(shù)據(jù)輸入端作為該模塊的第二組數(shù)據(jù)輸入端的(M+l)個(gè)輸入端；插 值位置計(jì)算器的輸出端作為該模塊的插值輸出端；第三插值器的數(shù)據(jù)輸出端作為該模塊的 數(shù)據(jù)輸出端；第二組插值器的(M+1)數(shù)據(jù)輸出端連接到(M+1)個(gè)LDPC譯碼SDM計(jì)算模塊 的輸入端；(M+l)個(gè)LDPC譯碼SDM計(jì)算模塊的輸出端連接到二維聯(lián)合定時(shí)誤差與頻偏估 計(jì)器的(M+1)數(shù)據(jù)輸入端；二維聯(lián)合定時(shí)誤差與頻偏估計(jì)器的輸出端連接到插值位置計(jì)算 器的輸入端；插值位置計(jì)算器的第一插值位置輸出端連接到第三插值器的插值位置輸入 端；插值位置計(jì)算器的插值位置第二組輸出端的(M+1)個(gè)輸出端連接到第二組插值器的 (M+l)個(gè)插值位置輸入端。所述的LDPC譯碼SDM計(jì)算模塊，包括LDPC譯碼后驗(yàn)信息計(jì)算器；取絕對(duì)值計(jì) 算器；累加器；取反器。該模塊的接口及內(nèi)部各器件的連接關(guān)系為L(zhǎng)DPC譯碼后驗(yàn)信息 計(jì)算器的輸入端作為該模塊的輸入端；取反處理器的輸出端作為該模塊的輸出端；LDPC 譯碼后驗(yàn)信息計(jì)算器的輸出端連接到取絕對(duì)值計(jì)算器的輸入端；取絕對(duì)值計(jì)算器的輸出端 連接到累加器的輸入端；累加器的輸出端連接到取反器的輸入端。所述的迭代定時(shí)控制模塊，具體包括首次接收判斷器；1:N路串并變換器；矩陣與 向量乘法模塊；存儲(chǔ)器；累加器；比較器；比較及迭代次數(shù)判斷器。該模塊的接口及內(nèi)部 各器件的連接關(guān)系為首次接收判斷器的輸入端作為該模塊的控制輸入端；1:N路串并變 換器的數(shù)據(jù)輸入端作為該模塊的數(shù)據(jù)輸入端；比較器的判決結(jié)果為"否"的輸出端及首次接 收判斷器的判決結(jié)果為"是"的輸出端作為該模塊的定時(shí)捕獲控制輸出端；比較及迭代次數(shù) 判斷器的判決結(jié)果為"否"的輸出端作為該模塊的定時(shí)跟蹤控制輸出端；比較及迭代次數(shù)判 斷器的判決結(jié)果為"是"的輸出端作為該模塊的最終輸出控制的輸出端。首次接收判斷器的 判決結(jié)果為"否"的輸出端作為該模塊的控制輸入端；1:N路串并變換器的N個(gè)輸出端連接 到矩陣與向量乘法模塊的N個(gè)向量輸入端；存儲(chǔ)器中LDPC碼校驗(yàn)矩陣H的非零數(shù)據(jù)分 別作為矩陣與向量乘法模塊的矩陣數(shù)據(jù)輸入端；矩陣與向量乘法模塊的R個(gè)向量輸出端連 接到累加器的輸入端；累加器的輸出端連接到比較器的輸入端；比較器的判決結(jié)果為"否" 的輸出端連接到比較及迭代次數(shù)判斷器的輸入端；其中，R， N為自然數(shù)，分別表示LDPC 碼校驗(yàn)方程的數(shù)量及LDPC碼的碼長(zhǎng)，LDPC碼的校驗(yàn)矩陣H為R行N列的矩陣，且有 R<N;矩陣與向量乘法模塊執(zhí)行二進(jìn)制數(shù)據(jù)模-2域的R行N列的矩陣與N維向量相乘運(yùn) 算，得到R維向量的輸出。本發(fā)明的特點(diǎn)及效果數(shù)字通信系統(tǒng)因受系統(tǒng)熱噪聲、傳輸時(shí)延和收發(fā)端頻率源偏差等因素影響，而使接收 的碼元符號(hào)速率與各種定時(shí)控制不能協(xié)同工作，就會(huì)導(dǎo)致通信系統(tǒng)性能降低甚至通信失 效。而本發(fā)明可在很大程度上糾正因以上定時(shí)不準(zhǔn)因素而引起的定時(shí)同步不準(zhǔn)確問題，大大改善數(shù)字通信系統(tǒng)中信號(hào)接收的誤碼性能。根據(jù)仿真表明在數(shù)字通信系統(tǒng)采用二進(jìn)制 移相鍵控(BPSK)調(diào)制的情況下，采用本發(fā)明的方法，LDPC編碼系統(tǒng)在較大定時(shí)誤差情況 下獲得的性能與理想同步情況下的編碼系統(tǒng)的性能偏差在0.2犯以下，滿足工程應(yīng)用中的 定時(shí)同步要求。同時(shí)，該方法還能有效的工作在使傳統(tǒng)定時(shí)同步方法失效的低SNR情況。 所以，該發(fā)明可作為較低信噪比(SNR)情況下有效可靠的定時(shí)同步方案，非常適用于深空、 衛(wèi)星通信及磁記錄系統(tǒng)等應(yīng)用場(chǎng)合。


本發(fā)明最好參考以下詳細(xì)描述和圖解本發(fā)明的附圖來(lái)理解，在這些附圖中 圖1為采用迭代定時(shí)同步的典型數(shù)字通信系統(tǒng)的基帶等效模型圖； 圖2為圖1中系統(tǒng)在無(wú)加性高斯白噪聲(AWGN)的理想情況下接收信號(hào)的波形及定時(shí) 誤差的示意圖；圖3為本發(fā)明的基于LDPC譯碼SDM的迭代定時(shí)同步裝置的實(shí)施例結(jié)構(gòu)圖；圖4為本發(fā)明的數(shù)據(jù)接收處理模塊的實(shí)施例結(jié)構(gòu)圖；圖5為本發(fā)明的LDPC碼幀同步相關(guān)檢測(cè)模塊的實(shí)施例結(jié)構(gòu)圖；圖6為本發(fā)明的迭代定時(shí)捕獲模塊的實(shí)施例結(jié)構(gòu)圖；圖7為本發(fā)明的LDPC譯碼SDM計(jì)算模塊的實(shí)施例結(jié)構(gòu)圖；圖8為本發(fā)明的迭代定時(shí)跟蹤模塊的實(shí)施例結(jié)構(gòu)圖；圖9為本發(fā)明的數(shù)字環(huán)路濾波器的實(shí)施例結(jié)構(gòu)圖；圖IO為本發(fā)明的定時(shí)誤差檢測(cè)模塊的實(shí)施例結(jié)構(gòu)圖；圖11為本發(fā)明的LDPC譯碼模塊的實(shí)施例結(jié)構(gòu)圖；圖12為本發(fā)明的迭代定時(shí)控制模塊的實(shí)施例結(jié)構(gòu)圖；圖13為本發(fā)明的實(shí)施例中采用1/2碼率LDPC碼(1944， 972)的BPSK通信系統(tǒng)進(jìn)行 基于LDPC碼SDM的迭代定時(shí)同步的仿真誤碼性能圖。
具體實(shí)施方式
本發(fā)明的方法及裝置結(jié)合附圖實(shí)施例詳細(xì)描述如下數(shù)字通信系統(tǒng)因受系統(tǒng)熱噪聲、傳輸時(shí)延和收發(fā)端頻率源偏差等因素影響，會(huì)因定時(shí) 控制不準(zhǔn)確而導(dǎo)致通信性能降低甚至失效，通?？刹捎枚〞r(shí)同步來(lái)解決該問題。為了提高定時(shí)同步的有效和可靠性，特別是考慮低SNR等情況，本發(fā)明采用基于LDPC譯碼的迭 代定時(shí)同步方法。該方法的具體實(shí)施如下1)將系統(tǒng)接收端收到的數(shù)據(jù)分別經(jīng)M倍符號(hào)頻率采樣、用插值方法進(jìn)行插值和均方 根升余弦(SRRC)匹配濾波得到濾波后數(shù)據(jù)；其中，M為自然數(shù)(可取為4等偶數(shù)，匹配濾 波所用的SRRC濾波器的滾降系數(shù)可根據(jù)工程實(shí)際取值為0.3-0.8);在首次執(zhí)行插值處理 時(shí)，不進(jìn)行插值；否則，用步驟4)反饋的插值位置進(jìn)行插值；2) 對(duì)步驟l)所述濾波后數(shù)據(jù)分別進(jìn)行兩路輸出第1路數(shù)據(jù)進(jìn)行LDPC碼幀同步相關(guān) 檢測(cè)方法處理，得到LDPC碼數(shù)據(jù)幀起始位置序號(hào)和幀同步檢測(cè)成功信號(hào)；第2路數(shù)據(jù)則 以相關(guān)檢測(cè)方法得到的LDPC碼數(shù)據(jù)幀起始位置序號(hào)為濾波后新數(shù)據(jù)的起始序號(hào)，將濾波 后新數(shù)據(jù)進(jìn)行1:M的串并變換，得到M路數(shù)據(jù)，且每路各N+L個(gè)數(shù)據(jù)；其中，N、 L都 為自然數(shù)，分別表示LDPC碼的幀長(zhǎng)和LDPC碼同步頭的長(zhǎng)度；3) 在接收到所述幀同步檢測(cè)成功信號(hào)后，進(jìn)行迭代定時(shí)控制方法處理:如為初次處理，則輸出定時(shí)捕獲控制信號(hào)，轉(zhuǎn)步驟4);否則轉(zhuǎn)步驟6);4) 在接收到所述定時(shí)捕獲控制信號(hào)后，對(duì)步驟2)所述M路數(shù)據(jù)用迭代定時(shí)捕獲方法 處理，得到插值位置和定時(shí)捕獲輸出數(shù)據(jù)，將得到的插值位置反饋到步驟1)，用于接收端 新接收數(shù)據(jù)的插值；5) 將所述定時(shí)捕獲輸出數(shù)據(jù)進(jìn)行K次LDPC迭代譯碼，得到N個(gè)LDPC碼字和(N-R) 個(gè)LDPC信息碼元；其中，K為自然數(shù)(可選為20 30); R為自然數(shù)，表示LDPC碼校驗(yàn)方 程的數(shù)量(即LDPC碼的校驗(yàn)矩陣為R行N列的矩陣)，且有R<N;6) 將0f得的N個(gè)LDPC碼字進(jìn)行迭代定時(shí)控制方法處理，得到定時(shí)捕獲控制信號(hào)、 定時(shí)跟蹤控制信號(hào)或最終LDPC譯碼輸出控制信號(hào)；若得到定時(shí)捕獲控制信號(hào)，設(shè)置定時(shí) 狀態(tài)為捕獲狀態(tài)，則轉(zhuǎn)步驟4);若得到定時(shí)跟蹤控制信號(hào)，設(shè)置定時(shí)狀態(tài)為跟蹤狀態(tài)，則 轉(zhuǎn)步驟7);若得到最終LDPC譯碼輸出控制信號(hào)，將所述的LDPC信息碼元，作為系統(tǒng)最 終的LDPC譯碼輸出數(shù)據(jù)，將其輸出，并轉(zhuǎn)步驟9);7) 將步驟2)所得的M路數(shù)據(jù)和N個(gè)LDPC碼字，進(jìn)行迭代定時(shí)跟蹤方法處理，得到 定時(shí)跟蹤輸出數(shù)據(jù)；8) 將步驟7)所述定時(shí)跟蹤輸出數(shù)據(jù)進(jìn)行K次LDPC迭代譯碼，得到N個(gè)LDPC碼字 和LDPC信息碼元，并轉(zhuǎn)(6);9) 接收端將接收到的新輸入數(shù)據(jù)，轉(zhuǎn)步驟l)，直到接收端接收數(shù)據(jù)完畢。 所述(2)中的LDPC碼幀同步相關(guān)檢測(cè)方法，用于LDPC碼幀同步相關(guān)檢測(cè)，得到LDPC碼數(shù)據(jù)幀起始位置序號(hào)和幀同步檢測(cè)成功信號(hào)；其特征在于，具體包括以下步驟(21) 將輸入的數(shù)據(jù)進(jìn)行1:M串并變換為M路數(shù)據(jù)；(22) 將(21)所得的M路數(shù)據(jù)分別與本地已知L長(zhǎng)LDPC碼幀同步頭數(shù)據(jù)(L為自然數(shù)， 表示序列長(zhǎng)度)進(jìn)行相關(guān)，得到M個(gè)相關(guān)值和對(duì)應(yīng)的數(shù)據(jù)序號(hào)；(23) 將(22)所得M個(gè)相關(guān)值進(jìn)行比較，得到最大的相關(guān)值和該相關(guān)值對(duì)應(yīng)的輸入數(shù)據(jù) 序號(hào)；(24) 將(23)所得最大相關(guān)值與給定相關(guān)門限V相比較若該最大相關(guān)值大于V，則得 到LDPC碼幀同步檢測(cè)成功信號(hào)；同時(shí)，將該最大相關(guān)值所對(duì)應(yīng)的輸入數(shù)據(jù)序號(hào)作為得到 的LDPC碼數(shù)據(jù)幀起始位置序號(hào)；其中，V為常數(shù)，取為(0.6 0.8)xL;在該方法中，主要利用了成型信號(hào)的特性。由附圖2可知，如暫不考慮AWGN的影 響，當(dāng)實(shí)際采樣點(diǎn)與最佳采樣點(diǎn)的誤差在符號(hào)周期T的12.5n/。以內(nèi)時(shí)，實(shí)際接收信號(hào)幅度與最佳接收信號(hào)的幅度差異在10%以內(nèi)，從而在相關(guān)處理時(shí)，易產(chǎn)生較大的相關(guān)結(jié)果。但 當(dāng)定時(shí)采樣的誤差在其他范圍時(shí)，實(shí)際幅度與最佳幅度的差異就比較大，就不易產(chǎn)生較大 的相關(guān)結(jié)果。另外，當(dāng)定時(shí)頻偏比較小時(shí)(如定時(shí)頻偏^2000ppm)，它對(duì)短數(shù)據(jù)序列(序列長(zhǎng) 度128SL5256)的相關(guān)檢測(cè)影響不大。所以，可利用這兩點(diǎn)，兼顧計(jì)算量和存儲(chǔ)量的要求， 采用M倍(M可選為4)上采樣數(shù)據(jù)來(lái)分別檢測(cè)當(dāng)前所接收的M路數(shù)據(jù)是否存在LDPC碼的 幀同步頭。所述(4)中的迭代定時(shí)捕獲方法，用于定時(shí)捕獲處理，得到插值位置和定時(shí)捕獲輸出數(shù) 據(jù)。其特征在于，包括以下步驟(41) 根據(jù)預(yù)設(shè)的初始定時(shí)誤差搜索范圍中的初始定時(shí)誤差和頻偏，分別計(jì)算得到M 路輸入數(shù)據(jù)的M+1組插值位置；其中，初始定時(shí)誤差搜索范圍中的初始定時(shí)誤差和頻偏的 取值如下定時(shí)誤差范圍選擇[-T/2，T/2]，且T為符號(hào)周期；定時(shí)頻偏范圍(可根據(jù)工程實(shí)際 情況)選擇[-2000ppm,2000ppm]范圍，其中，ppm表示百萬(wàn)分之一；插值位置的計(jì)算采用以 下方法第n個(gè)插值數(shù)據(jù)的插值位置-初始定時(shí)誤差+nx初始定時(shí)頻偏+n/M;(42) 將(41)所得M+l組插值位置分別對(duì)接收到的M路數(shù)據(jù)用插值方法進(jìn)行插值，得 到M+1路數(shù)據(jù)，且該數(shù)據(jù)每路各N個(gè)；(43) 將產(chǎn)生的M+1路數(shù)據(jù)，且該數(shù)據(jù)每路各N個(gè)；對(duì)每路分別進(jìn)行LDPC譯碼SDM 計(jì)算方法進(jìn)行處理，得到M+l個(gè)LDPC譯碼SDM;(44) 將(43)所得的M+l個(gè)LDPC譯碼SDM及對(duì)應(yīng)的(41)所述的初始定時(shí)誤差和頻偏， 用二維聯(lián)合定時(shí)誤差與頻偏估計(jì)方法進(jìn)行處理，得到定時(shí)誤差與頻偏；(45) 用(44)所得的定時(shí)誤差與頻偏，計(jì)算得到插值位置；其中，插值位置的計(jì)算采用 以下方法第n個(gè)插值數(shù)據(jù)的插值位置-定時(shí)誤差+nx定時(shí)頻偏+n/M;(46) 根據(jù)(45)所得的插值位置對(duì)所述的M路數(shù)據(jù)用插值方法進(jìn)行插值，得到定時(shí)捕獲 輸出數(shù)據(jù)；所述步驟(43)中的LDPC譯碼SDM計(jì)算方法，其特征在于，包括以下步驟(431) 對(duì)(43)所述的輸入長(zhǎng)為N的數(shù)據(jù)，進(jìn)行迭代譯碼處理，得到N個(gè)LDPC譯碼后 驗(yàn)信息；其中，迭代次數(shù)選3 5次；(432) 將(431)所得的N個(gè)LDPC譯碼的后驗(yàn)信息分別取絕對(duì)值，并把所得結(jié)果累加起 來(lái)，得到累加和； (433) 將(432)所得的累加和取反，得到LDPC譯碼SDM;所述(44)中的二維聯(lián)合定時(shí)誤差與頻偏估計(jì)方法，其特征在于，具體包括以下步驟(441) 將所述的M+l個(gè)LDPC譯碼SDM及對(duì)應(yīng)的定時(shí)誤差和頻偏用單純形極小值搜 索方法進(jìn)行處理，計(jì)算得到定時(shí)誤差和頻偏；(442) 根據(jù)預(yù)設(shè)的定時(shí)誤差和頻偏及設(shè)定的迭代次數(shù)(5 10次)，對(duì)所得的定時(shí)誤差和 頻偏及已進(jìn)行的迭代次數(shù)進(jìn)行判斷若所得定時(shí)誤差和頻偏滿足預(yù)設(shè)的定時(shí)誤差和頻偏要 求已進(jìn)行的迭代次數(shù)大于設(shè)定的迭代次數(shù)，則得到該方法最終輸出的定時(shí)誤差和頻偏；否其中，預(yù)設(shè)的誤差和頻偏條件如下M+l組定時(shí)誤差及頻偏序列中，最小和最大定時(shí)誤差及頻偏之間分別小于0.05T及O.Olx頻偏搜索區(qū)長(zhǎng)度，且T為符號(hào)周期；(443) 用(441)所得的定時(shí)誤差和頻偏，計(jì)算插值位置；其中，插值位置的計(jì)算采用以 下方法第n個(gè)插值數(shù)據(jù)的插值位置-定時(shí)誤差+nx定時(shí)頻偏+n/M;(444) 用(443)所得的插值位置，對(duì)所述的M路數(shù)據(jù)用插值方法進(jìn)行插值，得到插值后 數(shù)據(jù)；(445) 用(444)所得的插值后數(shù)據(jù)用LDPC譯碼SDM計(jì)算方法處理，得到LDPC譯碼 S薩；(446) 用(445)所得的LDPC譯碼SDM與(441)所述的定時(shí)誤差和頻偏用單純形極小值 搜索方法處理，得到更新后的定時(shí)誤差和頻偏；同時(shí)，將更新后的定時(shí)誤差和頻偏及對(duì)應(yīng) 的LDPC譯碼SDM的組合，根據(jù)單純形極小值搜索方法，取代原M+l個(gè)組合中性能最差 的LDPC譯碼所對(duì)應(yīng)的組合，并轉(zhuǎn)(441);其中，每個(gè)組合包括定時(shí)誤差和頻偏及相應(yīng)的 LDPC譯碼SDM。其中，該LDPC譯碼SDM的計(jì)算及最優(yōu)搜索的方法詳細(xì)描述如下 首先，本發(fā)明定義了一個(gè)LDPC譯碼SDM函數(shù)的代價(jià)函數(shù)^(r,U來(lái)判斷是否得到 最佳的定時(shí)采樣點(diǎn)。這個(gè)代價(jià)函數(shù)^(T,/p，)表示如下-z(13)其中，r為定時(shí)誤差估計(jì)，單位為秒；力，為定時(shí)頻偏估計(jì)，且單位為百萬(wàn)分之一的 赫茲；"C,k，/p，)為在給定r和力^前提下的變量節(jié)點(diǎn)Z'的先驗(yàn)信息，即對(duì)數(shù)似然比；丄(。,Wlr，/，)為在給定r和/^前提下第先次迭代中，從校驗(yàn)節(jié)點(diǎn)/到變量節(jié)點(diǎn)/的外信 息；丄(c,l2",/p，)和丄(。,Wlr,/p")的具體迭代計(jì)算可以見背景技術(shù)中LDPC的譯碼部分及 公式(1) 公式(10)。求和符號(hào)上標(biāo)Len為計(jì)算SDM所需要的少量計(jì)算次數(shù)，且計(jì)算次數(shù)可 選為3~5次。則對(duì)應(yīng)的最佳定時(shí)誤差估計(jì)r。p,與定時(shí)頻偏估計(jì)/p，，。p,可聯(lián)合求解以下優(yōu)化方 程(14)來(lái)得到[V, / , ] = Mg min {—， /鉀)} " 4、k 庁" re-r/2,7721 尸 (14) /爐，其中，argmin(/(;c》表示計(jì)算使函數(shù)/(;c)的值達(dá)到最小時(shí)自變量x的值，且/(;c)是與 x有關(guān)的函數(shù)Z定時(shí)誤差估計(jì)T的取值范圍為[-T/2,T/2]，且T為符號(hào)采樣周期，單位為秒； 定時(shí)頻偏估計(jì)/p，的頻率搜索范圍為可能存在的定時(shí)頻偏范圍， 一般可根據(jù)工程實(shí)際設(shè)置 該區(qū)間為[-2000 ppm ， 2000 ppm]。最后，對(duì)優(yōu)化方程(14)的計(jì)算可以采用背景技術(shù)中所描述的單純形極小值搜索算法來(lái) 具體實(shí)施通過以上初始'的搜索范圍定時(shí)誤差估計(jì)r和定時(shí)頻偏估計(jì)/p^，用單純形極小 值搜索算法搜索得到最佳的定時(shí)誤差估計(jì)r。p,與定時(shí)頻偏估計(jì)/，，。p,。所述步驟6)中的迭代定時(shí)控制方法，其特征在于，具體包括以下步驟(61) 若為首次接收到由所述的幀同步檢測(cè)成功信號(hào)，則輸出定時(shí)捕獲控制信號(hào)；否則， 轉(zhuǎn)(62);(62) 將所述的LDPC碼字的N個(gè)數(shù)據(jù)構(gòu)成的向量與LDPC碼的校驗(yàn)矩陣的轉(zhuǎn)置在模-2 域中相乘，得到一個(gè)R長(zhǎng)向量；其中，二進(jìn)制數(shù)據(jù)在模-2域中的運(yùn)算過程和普通運(yùn)算相同， 只是將最后的結(jié)果除以2取余數(shù)，并將該余數(shù)作為運(yùn)算結(jié)果；(63) 將(62)所得的R長(zhǎng)向量中的每個(gè)元素累加，得到累加和；(64) 將(63)所得的累加和與給定的判斷門限Q1進(jìn)行比較若累加和大于門限Q1，則 輸出定時(shí)捕獲控制信號(hào)；若累加和小于判斷門限Q1，將該累加和與給定的判斷門限Q2進(jìn) 行比較若該累加和小于判斷門限Q2或定時(shí)跟蹤方法的迭代次數(shù)等于設(shè)定的迭代次數(shù)Z, 則輸出最終LDPC譯碼輸出控制信號(hào)；否則，輸出定時(shí)跟蹤控制信號(hào)。其中，Ql， Q2為 實(shí)數(shù)，分別設(shè)為0.05R和0.01R; Z為自然數(shù)，設(shè)為5~10。所述步驟7)中的迭代定時(shí)跟蹤方法，其特征在于，具體包括以下步驟(71) 對(duì)所述(75)反饋的定時(shí)誤差，用以下公式計(jì)算插值位置第n個(gè)插值數(shù)據(jù)的插值位置=定時(shí)誤差+11/^1;其中，首次執(zhí)行該步驟時(shí)，無(wú)定時(shí)誤差反饋，可設(shè)定時(shí)誤差為0;(72) 用(71)所得的插值位置對(duì)所述的M路數(shù)據(jù)進(jìn)行插值，得到插值后的N個(gè)數(shù)據(jù)， 若(72)的執(zhí)行次數(shù)小于設(shè)定的迭代次數(shù)5~7次，則轉(zhuǎn)(73)，否則將得到插值后的N個(gè)數(shù)據(jù) 作為定時(shí)跟蹤輸出數(shù)據(jù)；(73) 將(72)所述的插值后數(shù)據(jù)和所述的N個(gè)LDPC碼字用Miller&M仙er定時(shí)誤差檢 測(cè)方法進(jìn)行計(jì)算，得到定時(shí)誤差估計(jì)；(74) 將(73)所述的定時(shí)誤差估計(jì)通過數(shù)字環(huán)路濾波處理，轉(zhuǎn)(71)。 所述各步驟中的插值方法，其特征在于，具體包括以下步驟(11) 對(duì)于所需插值的數(shù)據(jù)中的第n個(gè)數(shù)據(jù)，將該數(shù)據(jù)對(duì)應(yīng)的插值位置取整，得到一個(gè) 整數(shù)Zn;(12) 將該第n個(gè)數(shù)據(jù)的插值位置減去所述整數(shù)Zn,得到小數(shù)Qn;(13) 取所述M路數(shù)據(jù)中序號(hào)為Zn的M個(gè)并行數(shù)據(jù)，再增加一個(gè)序號(hào)為Zn+1的M 個(gè)并行數(shù)據(jù)中的第一個(gè)數(shù)據(jù)，共同組成M+1個(gè)樣本數(shù)據(jù)(對(duì)于原串行輸入數(shù)據(jù)，則為序號(hào) 從Zn*M+l開始的M+l個(gè)數(shù)據(jù))作為插值樣^:的函數(shù)值；(14) 采用線性分段插值等插值方法，以Qn作為所需插值的點(diǎn)，以樣本點(diǎn)序列 (0,1/M,…,U及(13)所得M+l個(gè)按次序的樣本數(shù)據(jù)值作為插值樣本序列，進(jìn)行插值計(jì)算， 得到插值結(jié)果。以下結(jié)合

，對(duì)本發(fā)明的裝置實(shí)施例進(jìn)行詳細(xì)說明圖l是采用迭代定時(shí)同步的典型數(shù)字通信系統(tǒng)的基帶等效模型圖。該模型的實(shí)現(xiàn)過程 如下首先，接收數(shù)據(jù)經(jīng)過現(xiàn)代信道譯碼得到譯碼后數(shù)據(jù)；其次，將該譯碼后數(shù)據(jù)進(jìn)行碼 字映射，如將{，0，、 ，1，}數(shù)據(jù)分別映射為{，1'、 ，-1，}數(shù)據(jù)，得到映射后數(shù)據(jù)；再次，將該映射后數(shù)據(jù)進(jìn)行成型濾波，得到成型濾波后數(shù)據(jù)；之后，將濾波后數(shù)據(jù)通過加性高斯白噪聲 (AWGN)信道，使發(fā)送的數(shù)據(jù)增加了噪聲；然后，將接收信號(hào)經(jīng)過匹配濾波處理，得到匹 配濾波后數(shù)據(jù)；之后，將匹配濾波后數(shù)據(jù)迸行迭代定時(shí)同步處理，得到最終的系統(tǒng)輸出數(shù) 據(jù)。其中，迭代定時(shí)同步主要是將該過程接收的數(shù)據(jù)，經(jīng)過抽樣軟判后，進(jìn)行現(xiàn)代信道譯 碼，再將譯碼結(jié)果反饋給定時(shí)同步過程，再將定時(shí)同步結(jié)果反饋給抽樣軟判處理，形成一 個(gè)迭代處理過程。該過程能充分利用定時(shí)同步與迭代譯碼的處理信息，使該系統(tǒng)在低SNR 情況下也能獲得滿意的定時(shí)同步效果；圖2是假設(shè)在無(wú)AWGN噪聲的理想情況下接收信號(hào)的波形及定時(shí)誤差示意圖；在本 發(fā)明中，由于通信系統(tǒng)發(fā)射和接收端都采用了對(duì)稱的均方根升余弦(SRRC)成型及匹配濾波 器，所以在該情況下，接收信號(hào)的波形類似升余弦波形結(jié)構(gòu)。另外，該圖也表示了因定時(shí) 不準(zhǔn)而造成的實(shí)際信號(hào)采樣點(diǎn)與最佳采樣點(diǎn)的偏差，如q。這將導(dǎo)致通信系統(tǒng)性能下降。 所以，數(shù)字通信系統(tǒng)需要高效的定時(shí)同步方法來(lái)恢復(fù)最佳定時(shí)采樣點(diǎn)。同時(shí)，定時(shí)誤差的 搜索范圍一般可取為[-T/2,T/2]。其中，T為符號(hào)周期；圖3為本發(fā)明的基于LDPC譯碼SDM的迭代定時(shí)同步裝置的實(shí)施例結(jié)構(gòu)圖。該裝置具體包括以下模塊數(shù)據(jù)接收處理模塊l; LDPC碼幀同步相關(guān)檢測(cè)模塊2;迭代定時(shí)捕獲模塊3;迭代定時(shí)控制模塊4;迭代定時(shí)跟蹤模塊5; LDPC譯碼模塊6。該裝置中各器件的連接關(guān)系及工作過程如下首先，將系統(tǒng)接收的數(shù)據(jù)在數(shù)據(jù)接收處理模塊中分別經(jīng)采樣器的M倍符號(hào)頻率采樣、數(shù)據(jù)緩存器的緩存數(shù)據(jù)，插值器1的插值和SRRC匹配濾波器的匹配濾波得到 濾波后的數(shù)據(jù)；其中，M為自然數(shù)，可取為4等偶數(shù)；在首次執(zhí)行插值處理時(shí)，不進(jìn)行插值；否則，用迭代定時(shí)捕獲模塊反饋的插值位置進(jìn)行插值。其次，對(duì)所得的濾波后數(shù)據(jù)分別進(jìn)行兩路輸出第一路數(shù)據(jù)用LDPC碼幀同步相關(guān)檢測(cè)模塊進(jìn)行LDPC碼幀同步相關(guān)檢測(cè)處 理，得到LDPC碼數(shù)據(jù)幀的起始位置序號(hào)和幀同步檢測(cè)成功信號(hào)；第二路數(shù)據(jù)則以第l路得 到的起始位置序號(hào)為濾波后新數(shù)據(jù)的起始序號(hào)，將其經(jīng)過1:M串并變換器，得到M路數(shù)據(jù)， 并保存在數(shù)據(jù)接收處理模塊中的M路數(shù)據(jù)緩存器。再次，在接收到上述幀同步檢測(cè)成功信 號(hào)后，用迭代定時(shí)控制模塊進(jìn)行迭代定時(shí)控制處理，輸出定時(shí)捕獲控制信號(hào)。之后，在接 收到定時(shí)捕獲控制信號(hào)有效后，對(duì)所述M路數(shù)據(jù)緩存器中的M路數(shù)據(jù)用迭代定時(shí)捕獲模塊 進(jìn)行迭代定時(shí)捕獲處理，得到插值位置和定時(shí)捕獲輸出數(shù)據(jù)。然后，將定時(shí)捕獲輸出數(shù)據(jù) 用LDPC譯碼模塊中的LDPC譯碼器進(jìn)行LDPC譯碼，得到LDPC碼字。之后，將得到的LDPC 碼字用迭代定時(shí)控制模塊進(jìn)行定時(shí)同步狀態(tài)的判斷與切換處理，得到以下信號(hào)定時(shí)捕獲 控制信號(hào)、定時(shí)跟蹤控制信號(hào)、數(shù)據(jù)選擇控制信號(hào)和最終LDPC譯碼輸出控制信號(hào)；然后，根據(jù)這些控制信號(hào)，分別進(jìn)行以下處理若定時(shí)捕獲控制信號(hào)有效，轉(zhuǎn)上述迭代定時(shí)捕獲 模塊進(jìn)行迭代定時(shí)捕獲處理；若定時(shí)跟蹤控制信號(hào)有效，用迭代定時(shí)跟蹤模塊進(jìn)行迭代定 時(shí)跟蹤的處理；若最終LDPC譯碼輸出控制信號(hào)有效，控制上述LDPC譯碼器，將所得的 LDPC信息碼元數(shù)據(jù)，作為系統(tǒng)最終的LDPC譯碼輸出數(shù)據(jù)并輸出。之后，繼續(xù)用迭代定時(shí) 捕獲模塊進(jìn)行定時(shí)狀態(tài)的判斷與切換和最終數(shù)據(jù)的輸出控制處理，直到數(shù)據(jù)接收結(jié)束。其中，用數(shù)據(jù)選擇控制信號(hào)控制LDPC譯碼模塊中的數(shù)據(jù)選擇器，選擇相應(yīng)的定時(shí)跟蹤輸出 數(shù)據(jù)或定時(shí)跟蹤輸出數(shù)據(jù)輸入到該數(shù)據(jù)選擇器。另外，這四個(gè)模塊可用圖4 圖12的具體圖 解來(lái)分別詳細(xì)說明。其中，在該圖中，普通線條表示串行一路數(shù)據(jù)，而標(biāo)有""^""的線 條表示多路數(shù)據(jù)，且該說明對(duì)以下各圖也適用。圖4為本發(fā)明的數(shù)據(jù)接收處理模塊的實(shí)施例結(jié)構(gòu)圖。該模塊具體包括以下器件M倍符號(hào)頻率采樣器；數(shù)據(jù)緩存器；插值器l;匹配濾波器；M路數(shù)據(jù)緩存器。該裝置中各器件 的連接關(guān)系及工作過程如下首先，將系統(tǒng)接收端收到的數(shù)據(jù)分別經(jīng)采樣器的M倍符號(hào)頻 率采樣處理，得到M倍采樣后的輸入數(shù)據(jù)；其次，將該M倍采樣后的輸入數(shù)據(jù)保存在數(shù)據(jù) 緩存器中，用于后續(xù)處理的數(shù)據(jù)緩存；再次，將該數(shù)據(jù)緩存器中的數(shù)據(jù)輸出，并用插值器 l進(jìn)行插值處理，得到插值后數(shù)據(jù)；其中，首次插值處理時(shí)，可不進(jìn)行插值處理；否則， 按外部輸入的插值位置，用現(xiàn)有的線形分段插值等方法對(duì)輸入數(shù)據(jù)進(jìn)行插值處理；之后， 將該插值后的數(shù)據(jù)用SRRC匹配濾波器進(jìn)行匹配濾波，得到濾波后數(shù)據(jù)，并將該濾波數(shù)據(jù) 輸出，用于相關(guān)檢測(cè)；最后，根據(jù)相關(guān)檢測(cè)獲得的LDPC數(shù)據(jù)幀起始數(shù)據(jù)序號(hào)，.以該序號(hào) 對(duì)應(yīng)的數(shù)據(jù)為起始數(shù)據(jù)，將該濾波后數(shù)據(jù)保存在M路數(shù)據(jù)緩存器中，用于后續(xù)處理。圖5為本發(fā)明的LDPC碼幀同步相關(guān)檢測(cè)模塊的實(shí)施例結(jié)構(gòu)圖。該模塊具體包括以下器件1:M路串并變換器；M個(gè)相關(guān)模塊；最大值判斷器；比較器。該裝置中各器件的連接關(guān)系及工作過程如下首先，將外部輸入的數(shù)據(jù)，輸入到1:M路串并變換器，輸出M 路數(shù)據(jù)。其次，將該M路數(shù)據(jù)分別輸入到M個(gè)相關(guān)模塊得到M個(gè)相關(guān)計(jì)算結(jié)果；再次， 將這M個(gè)相關(guān)計(jì)算結(jié)果用最大值判斷器判斷最大值，并得到該最大相關(guān)計(jì)算結(jié)果對(duì)應(yīng)的輸 入數(shù)據(jù)序列中第一個(gè)數(shù)據(jù)的序號(hào)n。其中，n為自然數(shù)，表示數(shù)據(jù)序號(hào)。最后，將該最大 值在比較器中判斷是否比相關(guān)門限V大若該最大值比相關(guān)門限V大，則輸出LDPC碼 的起始位置序號(hào)n和幀同步檢測(cè)成功信號(hào)。相關(guān)模塊具體包括以下器件移位寄存器；L 個(gè)乘法器；存儲(chǔ)器；累加器。該模塊主要用L長(zhǎng)的移位寄存器移位寄存輸入數(shù)據(jù)。在移位 寄存器中的L個(gè)數(shù)據(jù)分別用L個(gè)乘法器與本地?cái)?shù)據(jù)中對(duì)應(yīng)序號(hào)的數(shù)據(jù)相乘，并把該L個(gè)相 乘結(jié)果用累加器累加，得到相關(guān)結(jié)果，實(shí)現(xiàn)背景技術(shù)中的相關(guān)運(yùn)算功能。其中，L為自然 數(shù)，表示為相關(guān)運(yùn)算的長(zhǎng)度。(此處用于LDPC同步頭的相關(guān)運(yùn)算，根據(jù)同步頭長(zhǎng)度，可選 為128-256等數(shù))圖6為本發(fā)明的迭代定時(shí)捕獲模塊的實(shí)施例結(jié)構(gòu)圖。該模塊具體包括以下器件插值 器2(1) 2(M+1); LDPC譯碼軟判決信息(SDM)計(jì)算模塊1~(M+1); 二維聯(lián)合定時(shí)誤差和頻偏計(jì)算器；插值位置計(jì)算器；插值器3。該裝置中各器件的連接關(guān)系及工作過程如下該 過程如下首先，用插值位置計(jì)算器計(jì)算預(yù)設(shè)的定時(shí)誤差搜索范圍中的初始定時(shí)誤差和頻 偏，得到M路輸入數(shù)據(jù)的M+1組插值位置。其中，定時(shí)誤差和頻偏搜索范圍如下定時(shí)誤差估計(jì)范圍[-T/2,T/2]，且T為符號(hào)周期；定時(shí)頻偏估計(jì)范圍可根據(jù)實(shí)際的定時(shí)頻偏誤差 范圍選[-2000ppm,2000ppm]等范圍；初始的定時(shí)誤差和頻偏可在上述搜索范圍內(nèi)取M+l個(gè)均勻分布點(diǎn)(如可選擇M個(gè)點(diǎn)均勻分布在矩形搜索區(qū)的周長(zhǎng)上，剩下一個(gè)點(diǎn)設(shè)置在矩形搜 索區(qū)的中心位置)；第n個(gè)插值位置的計(jì)算如下定時(shí)誤差+nx定時(shí)頻偏+n/M; n為自然數(shù)， 為相應(yīng)數(shù)據(jù)的序號(hào)。其次，將這些插值位置用插值器2(1) 2(M+1)分別對(duì)保存于M路數(shù)據(jù) 緩存器的M路數(shù)據(jù)進(jìn)行插值，得到(M+1)組插值數(shù)據(jù)；再次，將該(M+1)組插值數(shù)據(jù)分別 輸入到LDPC譯碼SDM計(jì)算模塊1~(M+1)，輸出(M+1)個(gè)LDPC譯碼SDM。之后，將這 些LDPC譯碼SDM輸入到二維聯(lián)合定時(shí)誤差與頻偏估計(jì)器(LDPC譯碼SDM對(duì)應(yīng)的初始 定時(shí)誤差和頻偏已預(yù)先保存于該器件內(nèi))，經(jīng)過多次迭代后，得到定時(shí)誤差與頻偏。另外， 在該模塊中，二維聯(lián)合定時(shí)誤差與頻偏估計(jì)器可以用背景技術(shù)中單純形極小值搜索方法實(shí) 現(xiàn)，且該器件工作過程中會(huì)產(chǎn)生需要更新的定時(shí)誤差和頻偏及相應(yīng)的LDPC譯碼SDM。這 部分的實(shí)現(xiàn)可通過時(shí)分復(fù)用的方式(不同時(shí)間使用同一個(gè)器件，不同時(shí)間的使用效果相互之 間不影響)，將二維聯(lián)合定時(shí)誤差與頻偏估計(jì)器輸出的臨時(shí)的定時(shí)誤差和頻偏通過插值位置 計(jì)算器計(jì)算插值位置，并將該插值位置和上述M路輸入數(shù)據(jù)輸入到LDPC譯碼軟判決信息 (SDM)計(jì)算模塊(M+1)計(jì)算需要更新的LDPC譯碼SDM，最后將需要更新的定時(shí)誤差和頻 偏及相應(yīng)的LDPC譯碼SDM —起再用二維聯(lián)合定時(shí)誤差和頻偏計(jì)算器計(jì)算新的定時(shí)誤差 和頻偏，實(shí)現(xiàn)二維聯(lián)合定時(shí)誤差和頻偏計(jì)算器的迭代處理；其中，LDPC譯碼SDM計(jì)算模 塊可由圖7來(lái)詳細(xì)描述。圖7為本發(fā)明的LDPC譯碼SDM計(jì)算模塊的實(shí)施例結(jié)構(gòu)圖。該模塊具體包括以下器件LDPC譯碼后驗(yàn)信息計(jì)算模塊；取絕對(duì)值計(jì)算器；累加器B;取反器。該裝置中各器件的連接關(guān)系及工作過程如下首先，將輸入的數(shù)據(jù)通過LDPC譯碼后驗(yàn)信息計(jì)算模塊計(jì) 算得到N路LDPC譯碼后驗(yàn)信息數(shù)據(jù)；其次，將該N路LDPC譯碼后驗(yàn)信息數(shù)據(jù)通過取 絕對(duì)值計(jì)算器分別取絕對(duì)值，得到N路絕對(duì)值輸出數(shù)據(jù)；再次，將該N路絕對(duì)值輸出數(shù) 據(jù)通過累加器B累加，得到累加和；最后，將該累加和取反，得到LDPC譯碼SDM。其 中，LDPC譯碼后驗(yàn)信息計(jì)算模塊中的LDPC譯碼后驗(yàn)信息計(jì)算計(jì)算可由背景技術(shù)中公式 (l)-公式(ll)來(lái)描述。圖8為本發(fā)明的迭代定時(shí)跟蹤模塊的實(shí)施例結(jié)構(gòu)圖。該模塊具體包括以下器件插值器4;延遲器；定時(shí)誤差檢測(cè)器；數(shù)字環(huán)路濾波器。該裝置中各器件的連接關(guān)系及工作過 程如下首先，該模塊將該模塊自身反饋的定時(shí)誤差(首次處理設(shè)為O)，對(duì)輸入的M路數(shù) 據(jù)用插值器4進(jìn)行插值(插值位置定時(shí)誤差+n/M)，得到插值后數(shù)據(jù)；其次，接收到該插值后數(shù)據(jù)對(duì)應(yīng)的LDPC碼字(該LDPC碼字為該插值后數(shù)據(jù)進(jìn)行LDPC譯碼得到的)；再次， 將該插值后數(shù)據(jù)和LDPC碼字輸入到采用Miller&Mttller算法的定時(shí)誤差檢測(cè)器，得到定 時(shí)誤差估計(jì)。最后，將該定時(shí)誤差估計(jì)通過數(shù)字環(huán)路濾波器得到定時(shí)誤差，并將該輸出定 時(shí)誤差再次進(jìn)行插值等以上過程的迭代處理，直到迭代次數(shù)達(dá)到5~10次后，將該所得的 插值后的數(shù)據(jù)作為該模塊最后輸出的定時(shí)跟蹤輸出數(shù)據(jù)，并將其輸出。另外，在該模塊中， 數(shù)字環(huán)路濾波器可由圖9來(lái)詳細(xì)描述；定時(shí)誤差檢測(cè)器可由圖IO來(lái)詳細(xì)描述。圖9為本發(fā)明的數(shù)字環(huán)路濾波器的實(shí)施例結(jié)構(gòu)圖。該模塊具體包括以下器件積分器； 增益器G1;增益器G2;加法器。該裝置中各器件的連接關(guān)系及工作過程如下誤差信號(hào) 進(jìn)入該數(shù)字環(huán)路濾波器后，分為兩路第一路數(shù)據(jù)通過增益器Gi乘以了系數(shù)G"第二路經(jīng)過積分器積分后，再通過乘增益器G2乘以了系數(shù)G2。之后，將這兩路處理結(jié)果相加，作為數(shù)字環(huán)路濾波器的輸出。為了使數(shù)字環(huán)路濾波器穩(wěn)定，需要使該數(shù)字濾波器的傳輸函數(shù)的極點(diǎn)在單位圓內(nèi)。在本發(fā)明實(shí)施例中，可取G!和G2分別為1/4和1/64，能使該濾波 器穩(wěn)定有效的工作。圖10為本發(fā)明的定時(shí)誤差檢測(cè)模塊的實(shí)施例結(jié)構(gòu)圖。該模塊具體包括以下器件單 位延遲器；乘法器；減法器。該過程主要采用背景技術(shù)中Miller&Mttller定時(shí)誤差檢測(cè)方 法中公式(12)實(shí)現(xiàn)定時(shí)誤差檢測(cè)。該裝置中各器件的連接關(guān)系及工作過程如下首先，輸 入數(shù)據(jù)與LDPC碼字分別與對(duì)方通過單位延遲器進(jìn)行單位延遲后的數(shù)據(jù)用乘法器相乘得到 兩個(gè)數(shù)據(jù)。其次，將這兩個(gè)數(shù)據(jù)通過減法器進(jìn)行減法運(yùn)算，得到公式(12)所描述的定時(shí)誤 差估計(jì)。圖11為本發(fā)明的LDPC譯碼模塊的實(shí)施例結(jié)構(gòu)圖。該模塊具體包括以下器件數(shù)據(jù)選擇器；LDPC譯碼器。該裝置中各器件的連接關(guān)系及工作過程如下根據(jù)數(shù)據(jù)選擇控制信號(hào)控制數(shù)據(jù)選擇器中的輸入數(shù)據(jù)，并將該數(shù)據(jù)經(jīng)過LDPC譯碼器進(jìn)行譯碼，得到LDPC 碼字和LDPC信息碼元。LDPC譯碼器中執(zhí)行的LDPC譯碼過程可以由背景技術(shù)中公式(l)-公式(12)來(lái)描述。圖12為本發(fā)明的迭代定時(shí)控制模塊的實(shí)施例結(jié)構(gòu)圖。該模塊具體包括以下器件1:N 路串并變換器；存儲(chǔ)器；矩陣與向量乘法模塊；累加器；比較器；比較及迭代次數(shù)判斷器。 該裝置中各器件的連接關(guān)系及工作過程如下它在接收到LDPC碼幀同步檢測(cè)成功信號(hào)有 效后，立刻用定時(shí)捕獲控制信號(hào)控制迭代定時(shí)捕獲模塊開始工作。不然進(jìn)行以下處理首 先，將迭代定時(shí)跟蹤模塊輸出的LDPC碼字的數(shù)據(jù)用1 :N路串并變換器得到每次N個(gè)數(shù)據(jù)。 其次，將該N個(gè)數(shù)據(jù)與存儲(chǔ)器中的LDPC校驗(yàn)矩陣H的數(shù)據(jù)一起輸入到矩陣與向量乘法 模塊進(jìn)行矩陣與向量相乘的處理，得到R維向量。再次，將這R維向量的每個(gè)元素用累加 器累加，得到累加和。最后，將累加和輸入到比較器，將該累加和與比較器中設(shè)定的判斷 門限Q1(可設(shè)為0.05R)進(jìn)行比較若該累加和大于Q1，則輸出定時(shí)捕獲控制信號(hào)，控制迭 代定時(shí)捕獲模塊開始定時(shí)捕獲處理；否則，將該累加和輸入到比較及迭代次數(shù)判斷器，將 該累加和與比較及迭代次數(shù)判斷器中設(shè)定的判斷門限Q2(可設(shè)為O.OIR)進(jìn)行比較若該累 加小于Q2或定時(shí)跟蹤的迭代次數(shù)等于Z次(Z為自然數(shù)，可選為5~10),則輸出LDPC最 終數(shù)據(jù)輸出控制信號(hào)，控制LDPC譯碼器模塊，輸出整個(gè)系統(tǒng)最終的LDPC信息碼元；否 則，輸出定時(shí)跟蹤控制信號(hào)，控制迭代定時(shí)跟蹤模塊開始定時(shí)跟蹤處理。其中，R， N為 自然數(shù)，分別表示LDPC碼校驗(yàn)方程的數(shù)量及LDPC碼的碼長(zhǎng)(LDPC碼的校驗(yàn)矩陣H為R 行N列的矩陣)，且有IKN;矩陣與向量乘法模塊執(zhí)行二進(jìn)制數(shù)據(jù)模-2域的R行N列的矩 陣與N維向量相乘運(yùn)算，得到R維向量的輸出；其中的二進(jìn)制數(shù)據(jù)模-2域中的運(yùn)算和普通 運(yùn)算相同，只是將最后的結(jié)果除以2取余數(shù)，并將該余數(shù)作為運(yùn)算結(jié)果；另外，矩陣與向 量乘法模塊可由乘法器和異或器按上述運(yùn)算規(guī)則構(gòu)成，且為了減少計(jì)算量，可將值為"0" 的矩陣元素不參與運(yùn)算。圖13為本發(fā)明的實(shí)施例中采用1/2碼率LDPC碼(1944， 972)的BPSK通信系統(tǒng)進(jìn)行 基于LDPC碼SDM的迭代定時(shí)同步的仿真誤碼性能圖。在該圖中實(shí)線、虛線分別表示誤碼率和誤幀率。曲線"-"表示理想無(wú)定時(shí)誤差時(shí)的誤碼曲線；"一o—"表示在固定定時(shí)誤差為士0.5T時(shí)得到的誤碼曲線；曲線"一x—"表示在定時(shí)頻偏為士2000ppm時(shí)，采用本 發(fā)明方法仿真得到的誤碼曲線；"一 + —"表示在固定定時(shí)誤差與頻偏分別為i0.5T及 ±2000 ppm時(shí)得到的誤碼性能曲線。該仿真的系統(tǒng)模型如圖I所示。仿真的條件如下預(yù) 設(shè)定時(shí)誤差為定時(shí)誤差士0.5T;定時(shí)頻偏±2000卯m;定時(shí)誤差士0.5T及定時(shí)頻偏±2000 ppm;發(fā)送端所發(fā)送偽隨機(jī)數(shù)產(chǎn)生的BPSK調(diào)制信號(hào)；采用l/4采樣、25抽頭長(zhǎng)度和滾降 系數(shù)為0.3的SRRC成型與匹配濾波器。LDPC碼采用正EE802.11n標(biāo)準(zhǔn)草案中1/2碼率的 (1944,972)LDPC碼，且在譯碼過程中進(jìn)行20次迭代；LDPC譯碼采用置信傳播譯碼，并 用本發(fā)明戶萬(wàn)述的方法進(jìn)行定時(shí)誤差與頻偏估計(jì)。定時(shí)搜索路數(shù)采用5路搜索。仿真結(jié)束條 件為出現(xiàn)錯(cuò)誤的LDPC數(shù)據(jù)幀數(shù)量超過1萬(wàn)幀或總的仿真所用LDPC數(shù)據(jù)幀的數(shù)量超過80 萬(wàn)幀。仿真結(jié)果顯示，采用該方法，在預(yù)設(shè)定時(shí)誤差情況下，所得誤碼曲線相距理論曲線 值不超過0.2(18。且在較低SNR(Eb/N^ldB)情況下，該方法所得誤碼曲線相距理論曲線值 不超過O.ldB。所以，該方法可獲得較好的定時(shí)同步，并適用于較低SNR情況下的定時(shí)同 步等場(chǎng)合。
權(quán)利要求
1.一種基于低密度奇偶校驗(yàn)譯碼軟判決信息的迭代定時(shí)同步方法，其特征在于，該方法包括以下步驟1)將系統(tǒng)接收端收到的數(shù)據(jù)分別經(jīng)M倍符號(hào)頻率采樣、用插值方法進(jìn)行插值和均方根升余弦匹配濾波得到濾波后數(shù)據(jù)；其中，M為自然數(shù)；在首次執(zhí)行插值處理時(shí)，不進(jìn)行插值；2)對(duì)步驟1)所述濾波后數(shù)據(jù)分別進(jìn)行兩路輸出第1路數(shù)據(jù)進(jìn)行LDPC碼幀同步相關(guān)檢測(cè)方法處理，得到LDPC碼數(shù)據(jù)幀起始位置序號(hào)和幀同步檢測(cè)成功信號(hào)；第2路數(shù)據(jù)則以相關(guān)檢測(cè)方法得到的LDPC碼數(shù)據(jù)幀起始位置序號(hào)為濾波后新數(shù)據(jù)的起始序號(hào)，將濾波后新數(shù)據(jù)進(jìn)行1:M的串并變換，得到M路數(shù)據(jù)，且每路各N+L個(gè)數(shù)據(jù)；其中，N、L均為自然數(shù)，分別表示LDPC碼的幀長(zhǎng)和LDPC碼同步頭的長(zhǎng)度；3)在接收到所述幀同步檢測(cè)成功信號(hào)后，進(jìn)行迭代定時(shí)控制方法處理如為初次處理，則輸出定時(shí)捕獲控制信號(hào)，轉(zhuǎn)步驟4)；否則轉(zhuǎn)步驟6)；4)在接收到所述定時(shí)捕獲控制信號(hào)后，對(duì)步驟2)所述M路數(shù)據(jù)用迭代定時(shí)捕獲方法處理，得到插值位置和定時(shí)捕獲輸出數(shù)據(jù)，將得到的插值位置反饋到步驟1)，用于接收端新接收數(shù)據(jù)的插值；5)將所述定時(shí)捕獲輸出數(shù)據(jù)進(jìn)行K次LDPC迭代譯碼，得到N個(gè)LDPC碼字和(N-R)個(gè)LDPC信息碼元；其中，K為自然數(shù)；R為自然數(shù)，表示LDPC碼校驗(yàn)方程的數(shù)量，且有R＜N；6)將所得的N個(gè)LDPC碼字進(jìn)行迭代定時(shí)控制方法處理，得到定時(shí)捕獲控制信號(hào)、定時(shí)跟蹤控制信號(hào)或最終LDPC譯碼輸出控制信號(hào)；若得到定時(shí)捕獲控制信號(hào)，設(shè)置定時(shí)狀態(tài)為捕獲狀態(tài)，則轉(zhuǎn)步驟4)；若得到定時(shí)跟蹤控制信號(hào)，設(shè)置定時(shí)狀態(tài)為跟蹤狀態(tài)，則轉(zhuǎn)步驟7)；若得到最終LDPC譯碼輸出控制信號(hào)，將所述的LDPC信息碼元，作為系統(tǒng)最終的LDPC譯碼輸出數(shù)據(jù)，并輸出，并轉(zhuǎn)步驟9)；7)將步驟2)所得的M路數(shù)據(jù)和N個(gè)LDPC碼字，進(jìn)行迭代定時(shí)跟蹤方法處理，得到定時(shí)跟蹤輸出數(shù)據(jù)；8)將步驟7)所述定時(shí)跟蹤輸出數(shù)據(jù)進(jìn)行K次LDPC迭代譯碼，得到N個(gè)LDPC碼字和LDPC信息碼元，并轉(zhuǎn)(6)；9)接收端將接收到的新輸入數(shù)據(jù)，轉(zhuǎn)步驟1)，直到接收端接收數(shù)據(jù)完畢。
2、 如權(quán)利要求l所述方法，其特征在于，所述步驟2)中的LDPC碼幀同步相關(guān)檢測(cè) 方法，具體包括以下步驟(21) 將輸入的數(shù)據(jù)進(jìn)行1:M串并變換為M路數(shù)據(jù)；(22) 將(21)所得的M路數(shù)據(jù)分別與本地已知L長(zhǎng)LDPC碼幀同步頭數(shù)據(jù)，得到M個(gè) 相關(guān)值和對(duì)應(yīng)的數(shù)據(jù)序號(hào)，其中L為自然數(shù)，表示序列長(zhǎng)度進(jìn)行相關(guān)；(23) 將(22)所得M個(gè)相關(guān)值進(jìn)行比較，得到最大的相關(guān)值和該相關(guān)值對(duì)應(yīng)的輸入數(shù)據(jù) 序號(hào)；(24) 將(23)所得最大相關(guān)值與給定相關(guān)門限V相比較若該最大相關(guān)值大于V，貝幌 到LDPC碼幀同步檢測(cè)成功信號(hào)；同時(shí)，將該最大相關(guān)值所對(duì)應(yīng)的輸入數(shù)據(jù)序號(hào)作為得到 的LDPC碼數(shù)據(jù)幀起始位置序號(hào)；其中，V為常數(shù)。
3、 如權(quán)利要求l所述方法，其特征在于，所述步驟4)中的迭代定時(shí)捕獲方法，具體 包括以下步驟(41) 根據(jù)預(yù)設(shè)的初始定時(shí)誤差搜索范圍中的初始定時(shí)誤差和頻偏，分別計(jì)算得到M 路輸入數(shù)據(jù)的M+l組插值位置，且該每路數(shù)據(jù)為N個(gè)；其中，初始定時(shí)誤差搜索范圍中的初始定時(shí)誤差和頻偏的取值如下定時(shí)誤差范圍選擇[-T/2，T/2]，且T為符號(hào)周期；定時(shí)頻偏范圍選擇[-2000ppm,2000ppm]范圍，其中，ppm表示百萬(wàn)分之一；插值位置的計(jì)算采 用以下方法第n個(gè)插值數(shù)據(jù)的插值位置=初始定時(shí)誤差+11><初始定時(shí)頻偏+自；(42) 將(41)所得M+l組插值位置分別對(duì)接收到的M路數(shù)據(jù)用插值方法進(jìn)行插值，得 到M+1路數(shù)據(jù)；(43) 將產(chǎn)生的M+1路數(shù)據(jù)；對(duì)每路分別進(jìn)行LDPC譯碼SDM計(jì)算方法進(jìn)行處理，得 到M+l個(gè)LDPC譯碼SDM;(44) 將(43)所得的M+l個(gè)LDPC譯碼SDM及對(duì)應(yīng)的(41)所述的初始定時(shí)誤差和頻偏， 用二維聯(lián)合定時(shí)誤差與頻偏估計(jì)方法進(jìn)行處理，得到定時(shí)誤差與頻偏；(45) 用(44)所得的定時(shí)誤差與頻偏，計(jì)算得到插值位置；其中，插值位置的計(jì)算采用 以下方法第n個(gè)插值數(shù)據(jù)的插值位置-定時(shí)誤差+nx定時(shí)頻偏+n/M;(46) 根據(jù)(45)所得的插值位置對(duì)所述的M路數(shù)據(jù)用插值方法進(jìn)行插值，得到定時(shí)捕獲 輸出數(shù)據(jù)。
4、 如權(quán)利要求3所述方法，其特征在于，所述步驟(43)中的LDPC譯碼SDM計(jì)算方 法，具體包括以下步驟(431) 對(duì)(43)所述的輸入長(zhǎng)為N的數(shù)據(jù)，進(jìn)行迭代譯碼處理，得到N個(gè)LDPC譯碼后 驗(yàn)信息；(432) 將(431)所得的N個(gè)LDPC譯碼的后驗(yàn)信息分別取絕對(duì)值，并把所得結(jié)果累加起 來(lái)，得到累加和；(433) 將(432)所得的累加和取反，得到LDPC譯碼SDM。
5、 如權(quán)利要求1所述方法，其特征在于，所述步驟6)中的迭代定時(shí)控制方法，具體 包括以下步驟(61) 若為首次接收到由所述的幀同步檢測(cè)成功信號(hào)，則輸出定時(shí)捕獲控制信號(hào)；否則， 轉(zhuǎn)(62);(62) 將所述的LDPC碼字的N個(gè)數(shù)據(jù)構(gòu)成的向量與LDPC碼的校驗(yàn)矩陣的轉(zhuǎn)置在模-2 域中相乘，得到一個(gè)R長(zhǎng)向量；(63) 將(62)所得的R長(zhǎng)向量中的每個(gè)元素累加，得到累加和；(64) 將(63)所得的累加和與給定的判斷門限Q1進(jìn)行比較若累加和大于門限Q1，則 輸出定時(shí)捕獲控制信號(hào)；若累加和小于判斷門限Q1，將該累加和與給定的判斷門限Q2進(jìn) 行比較若該累加和小于判斷門限Q2或定時(shí)跟蹤方法的迭代次數(shù)等于設(shè)定的迭代次數(shù)Z， 則輸出最終LDPC譯碼輸出控制信號(hào)；否則，輸出定時(shí)跟蹤控制信號(hào)。其中，Ql， Q2為 實(shí)數(shù)；Z為自然數(shù)。
6、 采用如權(quán)利要求1所述方法的基于低密度奇偶校驗(yàn)譯碼軟判決信息的迭代定時(shí)同 步裝置，其特征在于，該裝置具體包括以下模塊數(shù)據(jù)接收處理模塊；LDPC碼幀同步相 關(guān)檢測(cè)模塊；迭代定時(shí)捕獲模塊；迭代定時(shí)控制模塊；迭代定時(shí)跟蹤模塊；LDPC譯碼模 塊。該裝置的連接關(guān)系為數(shù)據(jù)接收處理模塊的釆樣數(shù)據(jù)接收端作為整個(gè)系統(tǒng)的輸入端； 數(shù)據(jù)接收處理模塊的匹配濾波輸出端連接到LDPC碼幀同步相關(guān)檢測(cè)模塊的數(shù)據(jù)輸入端； 數(shù)據(jù)接收處理模塊中M路第一數(shù)據(jù)輸出端連接到迭代定時(shí)跟蹤模塊中的插值數(shù)據(jù)輸入端； 數(shù)據(jù)接收處理模塊的M路第二數(shù)據(jù)輸出端連接到迭代定時(shí)捕獲模塊中的第一數(shù)據(jù)輸入端； 數(shù)據(jù)接收處理模塊的M路數(shù)據(jù)第三組輸出端共(M+l)個(gè)輸出端連接到迭代定時(shí)捕獲模塊的第二組數(shù)據(jù)輸入端的(M+l)個(gè)輸入端；LDPC碼幀同步相關(guān)檢測(cè)模塊的數(shù)據(jù)輸出端連接 到數(shù)據(jù)接收處理模塊的匹配濾波輸入端；LDPC碼幀同步相關(guān)檢測(cè)模塊的控制信號(hào)輸出端 連接到迭代定時(shí)控制模塊的控制輸入端；迭代定時(shí)捕獲模塊的數(shù)據(jù)輸出端連接到LDPC譯 碼模塊的第二數(shù)據(jù)輸入端；迭代定時(shí)捕獲模塊的插值輸出端連接到數(shù)據(jù)接收處理模塊的插 值輸入端；迭代定時(shí)跟蹤模塊的數(shù)據(jù)輸出端連接到LDPC譯碼模塊的第一數(shù)據(jù)輸入端； LDPC譯碼模塊的第一碼字輸出端連接到迭代定時(shí)跟蹤模塊的碼字?jǐn)?shù)據(jù)輸入端；LDPC譯 碼模塊的第二碼字輸出端連接到迭代定時(shí)控制模塊的數(shù)據(jù)輸入端；LDPC譯碼模塊的信息 碼元輸出端作為系統(tǒng)的輸出端；迭代定時(shí)控制模塊的定時(shí)捕獲控制輸出端連接到迭代定時(shí) 捕獲模塊的控制輸入端；迭代定時(shí)控制模塊的定時(shí)跟蹤控制輸出端連接到迭代定時(shí)跟蹤模 塊的控制輸入端；迭代定時(shí)控制模塊的最終輸出控制輸出端連接到LDPC譯碼模塊的控制 輸入端；迭代定時(shí)控制模塊的數(shù)據(jù)選擇控制輸出端連接到LDPC譯碼模塊的數(shù)據(jù)選擇控制 輸入端；其中，M為自然數(shù)。
7、 如權(quán)利要求6所述的裝置，其特征在于，所述的LDPC碼幀同步相關(guān)檢測(cè)模塊， 具體包括1:M路串并變換器；M個(gè)相關(guān)模塊；最大值判斷器；比較器。該模塊的接口及 內(nèi)部各器件的連接關(guān)系為1:M路串并變換器的數(shù)據(jù)輸入端作為該模塊的數(shù)據(jù)輸入端；比 較器的數(shù)據(jù)輸出端作為該模塊的數(shù)據(jù)輸出端；比較器的控制信號(hào)輸出端作為該模塊的控制 信號(hào)輸出端；1:M路串并變換器的M個(gè)輸出端分別連接到M個(gè)相關(guān)模塊的輸入端；M個(gè) 相關(guān)模塊的輸出端連接到最大值判斷器的M個(gè)輸入端；最大值判斷器的輸出端連接到比較 器的輸入端；其中，相關(guān)模塊用于輸入數(shù)據(jù)與本地己知數(shù)據(jù)序列的相關(guān)計(jì)算，得到相關(guān)計(jì) 算結(jié)果。
8、 如權(quán)利要求6所述的裝置，其特征在于，所述的迭代定時(shí)捕獲模塊，具體包括第二組插值器，共包括(M+1)個(gè)插值器；(M+1)個(gè)LDPC譯碼軟判決信息(SDM)計(jì)算器；二 維聯(lián)合定時(shí)誤差和頻偏計(jì)算器；插值位置計(jì)算器；第三插值器。該模塊的接口及內(nèi)部各器 件的連接關(guān)系為該模塊有1個(gè)控制該模塊工作的控制輸入端；第三插值器的數(shù)據(jù)輸入端 作為該迭代定時(shí)捕獲模塊的第一數(shù)據(jù)輸入端；第二組插值器的(M+1)數(shù)據(jù)輸入端作為該模 塊的第二組數(shù)據(jù)輸入端的(M+l)個(gè)輸入端；插值位置計(jì)算器的輸出端作為該模塊的插值輸 出端；第三插值器的數(shù)據(jù)輸出端作為該模塊的數(shù)據(jù)輸出端；第二組插值器的(M+1)數(shù)據(jù)輸 出端連接到(M+1)個(gè)LDPC譯碼SDM計(jì)算模塊的輸入端；(M+l)個(gè)LDPC譯碼SDM計(jì)算 模塊的輸出端連接到二維聯(lián)合定時(shí)誤差與頻偏估計(jì)器的(M+1)數(shù)據(jù)輸入端；二維聯(lián)合定時(shí) 誤差與頻偏估計(jì)器的輸出端連接到插值位置計(jì)算器的輸入端；插值位置計(jì)算器的第一插值 位置輸出端連接到第三插值器的插值位置輸入端；插值位置計(jì)算器的插值位置第二組輸出 端的(M+1)個(gè)輸出端連接到第二組插值器的(M+1)個(gè)插值位置輸入端。
9、 如權(quán)利要求8所述的裝置，其特征在于，所述的LDPC譯碼SDM計(jì)算模i央，包括 LDPC譯碼后驗(yàn)信息計(jì)算器；取絕對(duì)值計(jì)算器；累加器；取反器。該模塊的接口及內(nèi)部各 器件的連接關(guān)系為L(zhǎng)DPC譯碼后驗(yàn)信息計(jì)算器的輸入端作為該模塊的輸入端；'取反處理 器的輸出端作為該模塊的輸出端；LDPC譯碼后驗(yàn)信息計(jì)算器的輸出端連接到取絕對(duì)值計(jì) 算器的輸入端；取絕對(duì)值計(jì)算器的輸出端連接到累加器的輸入端；累加器的輸出端連接到 取反器的輸入端。
10、 如權(quán)利要求6所述的裝置，其特征在于，所述的迭代定時(shí)控制模塊，具體包括 首次接收判斷器；1:N路串并變換器；矩陣與向量乘法模塊；存儲(chǔ)器；累加器；比較器； 比較及迭代次數(shù)判斷器。該模塊的接口及內(nèi)部各器件的連接關(guān)系為首次接收判斷器的輸 入端作為該模塊的控制輸入端；1:N路串并變換器的數(shù)據(jù)輸入端作為該模塊的數(shù)據(jù)輸入端； 比較器的判決結(jié)果為"否"的輸出端及首次接收判斷器的判決結(jié)果為"是"的輸出端作為該模 塊的定時(shí)捕獲控制輸出端；比較及迭代次數(shù)判斷器的判決結(jié)果為"否"的輸出端作為該模塊 的定時(shí)跟蹤控制輸出端；比較及迭代次數(shù)判斷器的判決結(jié)果為"是"的輸出端作為該模塊的 最終輸出控制的輸出端。首次接收判斷器的判決結(jié)果為"否"的輸出端作為該模塊的控制輸 入端；1:N路串并變換器的N個(gè)輸出端連接到矩陣與向量乘法模塊的N個(gè)向量輸入端；存儲(chǔ)器中LDPC碼校驗(yàn)矩陣H的非零數(shù)據(jù)分別作為矩陣與向量乘法模塊的矩陣數(shù)據(jù)輸入端； 矩陣與向量乘法模塊的R個(gè)向量輸出端連接到累加器的輸入端；累加器的輸出端連接到比較器的輸入端；比較器的判決結(jié)果為"否"的輸出端連接到比較及迭代次數(shù)判斷器的輸入端； 其中，R， N為自然數(shù)，分別表示LDPC碼校驗(yàn)方程的數(shù)量及LDPC碼的碼長(zhǎng)，LDPC碼 的校驗(yàn)矩陣H為R行N列的矩陣，且有IKN;矩陣與向量乘法模塊執(zhí)行二進(jìn)制數(shù)據(jù)模-2 域的R行N列的矩陣與N維向量相乘運(yùn)算，得到R維向量的輸出。
全文摘要
本發(fā)明涉及基于低密度奇偶校驗(yàn)(LDPC)譯碼軟判決信息(SDM)的迭代定時(shí)同步方法與裝置，屬于數(shù)字通信技術(shù)領(lǐng)域，其特征在于利用LDPC譯碼中的SDM來(lái)實(shí)現(xiàn)定時(shí)同步。該發(fā)明的內(nèi)容如下首先，將系統(tǒng)接收的信號(hào)經(jīng)數(shù)據(jù)接收處理模塊1的高倍采樣，插值和匹配濾波后，用LDPC碼幀同步檢測(cè)模塊獲得LDPC碼幀起始序號(hào)并觸發(fā)同步處理。其次，它用迭代定時(shí)捕獲模塊2實(shí)現(xiàn)迭代定時(shí)捕獲。之后，它用迭代定時(shí)跟蹤模塊3實(shí)現(xiàn)迭代定時(shí)跟蹤。最后，它用迭代定時(shí)控制模塊4判斷及切換定時(shí)捕獲與跟蹤，并控制LDPC譯碼模塊輸出系統(tǒng)最終的譯碼數(shù)據(jù)。該發(fā)明可應(yīng)用于在較低信噪比(SNR)及高速數(shù)字通信情況下的定時(shí)同步，且特別適用于較低SNR情況下衛(wèi)星、深空通信及高速磁記錄系統(tǒng)等應(yīng)用場(chǎng)合。
文檔編號(hào)H04L7/02GK101217355SQ20081005588
公開日2008年7月9日 申請(qǐng)日期2008年1月11日 優(yōu)先權(quán)日2008年1月11日
發(fā)明者包建榮, 詹亞鋒, 邢騰飛, 陸建華 申請(qǐng)人:清華大學(xué)