積累互相關(guān)包絡(luò)對(duì)齊的8核dsp片上并行實(shí)現(xiàn)方法
【專利摘要】本發(fā)明涉及多核DSP并行方法。本發(fā)明提供一種積累互相關(guān)包絡(luò)對(duì)齊的8核DSP片上并行實(shí)現(xiàn)方法，首先，DSP核0對(duì)一幀雷達(dá)回波數(shù)據(jù)的前Z個(gè)PRF的脈沖回波作積累互相關(guān)包絡(luò)對(duì)齊處理以獲得穩(wěn)定的全局基準(zhǔn)信號(hào)refg，為核0?7后續(xù)對(duì)齊算法中共用，其中Z為8的倍數(shù)；其次，DSP核0?7根據(jù)全局基準(zhǔn)信號(hào)refg并行分別處理后續(xù)8個(gè)PRF的積累互相關(guān)包絡(luò)對(duì)齊，每個(gè)核分別得到對(duì)齊后的回波信號(hào)和局部基準(zhǔn)信號(hào)refI，I＝[0、1、2、3、4、5、6、7]，代表當(dāng)前核序號(hào)，然后，DSP核0檢測(cè)一幀所有脈沖信號(hào)是否都已經(jīng)對(duì)齊，是則輸出對(duì)齊后的回波信號(hào)，否則DSP核0利用八核輸出的八個(gè)基準(zhǔn)信號(hào)ref【0?7】更新全局基準(zhǔn)信號(hào)，重復(fù)上述步驟。適用于積累互相關(guān)包絡(luò)對(duì)齊的8核DSP片上并行。
【專利說明】
積累互相關(guān)包絡(luò)對(duì)齊的8核DSP片上并行實(shí)現(xiàn)方法
技術(shù)領(lǐng)域
[0001 ]本發(fā)明涉及多核DSP并行方法，特別涉及多核DSP并行實(shí)現(xiàn)一種大計(jì)算量、逐個(gè)脈 沖遞歸運(yùn)算的積累互相關(guān)包絡(luò)對(duì)齊的方法。
【背景技術(shù)】
[0002] 在逆合成孔徑雷達(dá)成像(ISAR)中，包絡(luò)對(duì)齊的優(yōu)劣會(huì)對(duì)最后的成像質(zhì)量產(chǎn)生重大 影響。對(duì)于相鄰兩次回波，目標(biāo)的轉(zhuǎn)角一般小于0.01〇,因此它們的實(shí)包絡(luò)十分相似，其相關(guān) 系數(shù)一般達(dá)0.95以上?；谏鲜鍪聦?shí)，國內(nèi)外許多專家學(xué)者對(duì)包絡(luò)對(duì)齊算法進(jìn)行了大量的 研究，從公開發(fā)表的文獻(xiàn)可以看出，主要有以下三種途徑：
[0003] 1、互相關(guān)函數(shù)最大準(zhǔn)則；
[0004] 2、最小熵準(zhǔn)則；
[0005] 3、最小歐式距離準(zhǔn)則。
[0006] 本發(fā)明主要涉及研究得最多、應(yīng)用得最廣泛的第一種準(zhǔn)則，并且采用了一種能有 效避免包絡(luò)漂移和突跳誤差的積累基準(zhǔn)信號(hào)互相關(guān)包絡(luò)對(duì)齊計(jì)算方法。該計(jì)算方法最大的 特點(diǎn)是其基準(zhǔn)信號(hào)不是單一的某個(gè)回波信號(hào)，而是加權(quán)積累實(shí)時(shí)更新的。這樣能很大程度 地消融異常信號(hào)在基準(zhǔn)信號(hào)中的權(quán)重，可有效消除目標(biāo)回波中起伏部分對(duì)運(yùn)動(dòng)補(bǔ)償?shù)挠?響；同時(shí)基準(zhǔn)信號(hào)中與待對(duì)齊回波越鄰近的信號(hào)占越高的權(quán)值，能獲得更好的對(duì)齊效果。缺 點(diǎn)就是此算法只能按脈沖回波順次運(yùn)算，不能并行計(jì)算提高效率。
[0007] 近年來，各種算法對(duì)DSP平臺(tái)的要求日益提高，TI公司在2010年11月推出了業(yè)界性 能最高的多核DSP處理芯片TMS320C6678，一舉超越了業(yè)界所有其他DSP內(nèi)核，成為首款獲得 最高定點(diǎn)性能和浮點(diǎn)性能平分的DSP。其擁有8個(gè)采用Keystone架構(gòu)的C66x內(nèi)核，每個(gè)內(nèi)核 最高運(yùn)行在1.25GHz的最高頻率上，提供160GFL0PS的高性能。得益于C6678強(qiáng)大的性能，將 其應(yīng)用于ISAR中的包絡(luò)對(duì)齊，能有效解決其大數(shù)據(jù)量、復(fù)雜算法的問題。

【發(fā)明內(nèi)容】

[0008] 本發(fā)明所要解決的技術(shù)問題，就是提供一種積累互相關(guān)包絡(luò)對(duì)齊的8核DSP片上并 行實(shí)現(xiàn)方法以實(shí)現(xiàn)在不影響算法性能的前提下，將順序執(zhí)行的積累基準(zhǔn)信號(hào)互相關(guān)包絡(luò)對(duì) 齊算法并行地在8核DSP上實(shí)現(xiàn)，大大提尚其運(yùn)算時(shí)間。
[0009] 本發(fā)明解決所述技術(shù)問題，采用的技術(shù)方案是，積累互相關(guān)包絡(luò)對(duì)齊的8核DSP片 上并行實(shí)現(xiàn)方法，包括以下步驟：
[0010] 步驟1、DSP核0對(duì)一幀雷達(dá)回波數(shù)據(jù)的前Z個(gè)PRF的脈沖回波作積累互相關(guān)包絡(luò)對(duì) 齊處理以獲得穩(wěn)定的全局基準(zhǔn)信號(hào)ref g，為核0-7后續(xù)對(duì)齊算法中共用，其中Z為8的倍數(shù)；
[0011] 步驟2、DSP核0-7根據(jù)全局基準(zhǔn)信號(hào)refg并行分別處理后續(xù)8個(gè)PRF的積累互相關(guān) 包絡(luò)對(duì)齊，每個(gè)核分別得到對(duì)齊后的回波信號(hào)和局部基準(zhǔn)信號(hào)^心，1 = [0、1、2、3、4、5、6、 7 ]，代表當(dāng)前核序號(hào)，進(jìn)入步驟3;
[0012] 步驟3、DSP核0檢測(cè)一幀所有脈沖信號(hào)是否都已經(jīng)對(duì)齊，是則輸出對(duì)齊后的回波信 號(hào)，否則進(jìn)入步驟4;
[0013] 步驟4，DSP核0利用八核輸出的八個(gè)基準(zhǔn)信號(hào)ref【o-7】更新全局基準(zhǔn)信號(hào)，重復(fù)步驟 2〇
[0014] 具體的，所述步驟1中，DSP核0對(duì)一幀雷達(dá)回波數(shù)據(jù)的前Z個(gè)PRF的脈沖回波作積累 互相關(guān)包絡(luò)對(duì)齊處理以獲得穩(wěn)定的全局基準(zhǔn)信號(hào)ref g，包括以下步驟：
[0015] 針對(duì)待處理的一幀N*M維的雷達(dá)回波數(shù)據(jù)，其中，N為回波脈沖個(gè)數(shù)，Μ為每個(gè)脈沖 所含的距離單元數(shù)；以第一個(gè)脈沖為初始的基準(zhǔn)信號(hào)ref，取第二個(gè)脈沖回波作為當(dāng)前回波 記為St;
[0016] A1，核0計(jì)算基準(zhǔn)信號(hào)ref與當(dāng)前回波St的互相關(guān)函數(shù)R，并獲取互相關(guān)函數(shù)峰值 Rmax及其對(duì)應(yīng)的坐標(biāo)記做Xmax，進(jìn)入步驟B1;
[0017] B1，核0分別取出峰值點(diǎn)(Xmax，Rmax)左右相鄰的互相關(guān)函數(shù)值及其對(duì)應(yīng)的坐標(biāo)分別 記為(Xi，Ri)及(X2，R2)，這三點(diǎn)可擬合出一條凹口向下的二次項(xiàng)曲線，解算出該曲線的最大 值點(diǎn)對(duì)應(yīng)的坐#
，進(jìn)入步驟C1;
[0018] C1，核0利用此時(shí)的坐標(biāo)X'構(gòu)造補(bǔ)償項(xiàng):^并與當(dāng)前回波st點(diǎn)乘(對(duì)應(yīng)距離 單元相乘）
以獲得對(duì)齊后的當(dāng)前脈沖回波st'并存儲(chǔ)，進(jìn)入步驟D1;
[0019] D1,根據(jù)對(duì)齊后的當(dāng)前脈沖回波St'計(jì)算新的基準(zhǔn)信號(hào):ref = γ · ref+st'，其中， 更新權(quán)系數(shù)〇< γ <1;至此當(dāng)前脈沖回波的包絡(luò)對(duì)齊完成，進(jìn)入步驟El;
[0020] El，核0取出下一個(gè)脈沖回波，根據(jù)新的基準(zhǔn)信號(hào)，重復(fù)步驟A，直到前Z個(gè)回波數(shù)據(jù) 全部完成對(duì)齊，獲得穩(wěn)定的全局基準(zhǔn)信號(hào)refg。
[0021] 進(jìn)一步的，在積累互相關(guān)包絡(luò)對(duì)齊算法的并行處理中，還包括利用快速傅里葉變 換(FFT)對(duì)基準(zhǔn)信號(hào)r ef與當(dāng)前回波s t進(jìn)行插值處理：
[0022]假設(shè)，包絡(luò)對(duì)齊的精度要求為1/8個(gè)距離單元，則需要對(duì)基準(zhǔn)信號(hào)ref與當(dāng)前回波 st進(jìn)行8倍插值，在DSP中可采用兩次不同點(diǎn)數(shù)的FFT來實(shí)現(xiàn)8倍插值，8卩：IFFT[FFT( st，M)， 8M]，其中，Μ為每個(gè)脈沖所含的距離單元數(shù)。
[0023]進(jìn)一步的，核0利用快速傅里葉變換(FFT)提高互相關(guān)函數(shù)計(jì)算效率；
[0024] 計(jì)算基準(zhǔn)信號(hào)ref與當(dāng)前回波st的互相關(guān)函數(shù)R時(shí)，利用3次FFT即可實(shí)現(xiàn)，即：R = IFFT[FFT(st).*FFT(ref)]〇
[0025]具體的，所述步驟2-4中，包括以下步驟：
[0026] A2,核0-7分別計(jì)算全局基準(zhǔn)信號(hào)refg與當(dāng)前回波st的互相關(guān)函數(shù)R，并獲取互相關(guān) 函數(shù)峰值Rmax及其對(duì)應(yīng)的坐標(biāo)記做Xmax，進(jìn)入步驟B2;
[0027] B2,核0-7分別取出峰值點(diǎn)(Xmax，Rmax)左右相鄰的互相關(guān)函數(shù)值及其對(duì)應(yīng)的坐標(biāo)分 別記為(Xi，Ri)及(X2，R2)，這三點(diǎn)可擬合出一條凹口向下的二次項(xiàng)曲線，解算出該曲線的最 大值點(diǎn)對(duì)應(yīng)的坐標(biāo)
丨入步驟C2;
[0028] C2,核0-7分別利用此時(shí)的坐標(biāo)X'構(gòu)造補(bǔ)償項(xiàng):e-\并與當(dāng)前回波st點(diǎn)乘(對(duì) 應(yīng)距離單元相乘）
以獲得對(duì)齊后的當(dāng)前脈沖回波St'并存儲(chǔ)，以此八 核分別得到八個(gè)對(duì)齊后的脈沖回波，并分別根據(jù)對(duì)齊后的當(dāng)前脈沖回波St'計(jì)算對(duì)應(yīng)核的 局部基準(zhǔn)信號(hào):ref I = St至此當(dāng)前脈沖回波的包絡(luò)對(duì)齊完成，進(jìn)入步驟D2;
[0029] D2,核0根據(jù)八核輸出的局部基準(zhǔn)信號(hào)更新全局基準(zhǔn)信號(hào)：
[0030] refg= γ 8refg+ γ 7ref7+ γ 6ref6+ γ 5ref5+ γ 4ref4+ γ 3ref3+ γ 2ref 2+ γ Vef 1 + ref0；0^i γ ^Ξ1 ；
[0031] 其中，refO為核0計(jì)算出的局部基準(zhǔn)信號(hào)，refl為核1計(jì)算出的局部基準(zhǔn)信號(hào)，ref2 為核2計(jì)算出的局部基準(zhǔn)信號(hào)，ref 3為核3計(jì)算出的局部基準(zhǔn)信號(hào)，ref 4為核4計(jì)算出的局部 基準(zhǔn)信號(hào)，ref5為核5計(jì)算出的局部基準(zhǔn)信號(hào)，ref6為核6計(jì)算出的局部基準(zhǔn)信號(hào)，ref7為核 7計(jì)算出的局部基準(zhǔn)信號(hào)；
[0032] 取出下一個(gè)脈沖回波，根據(jù)新的全局基準(zhǔn)信號(hào)，重復(fù)步驟A2,直到一幀回波數(shù)據(jù)全 部完成對(duì)齊，輸出對(duì)齊后的回波信號(hào)。
[0033]進(jìn)一步的，在積累互相關(guān)包絡(luò)對(duì)齊算法的并行處理中，還包括利用傅里葉變化對(duì) 全局基準(zhǔn)信號(hào)refg與當(dāng)前回波st進(jìn)行插值處理：
[0034]假設(shè)，包絡(luò)對(duì)齊的精度要求為1/8個(gè)距離單元，則需要對(duì)全局基準(zhǔn)信號(hào)refg與當(dāng)前 回波st進(jìn)行8倍插值，在DSP中可采用兩次不同點(diǎn)數(shù)的FFT來實(shí)現(xiàn)8倍插值，8卩：IFFT[FFT( st， Μ)，8M]，其中，Μ為每個(gè)脈沖所含的距離單元數(shù)。
[0035]進(jìn)一步的，核0利用FFT提尚互相關(guān)函數(shù)計(jì)算效率；
[0036]計(jì)算全局基準(zhǔn)信號(hào)refg與當(dāng)前回波st的互相關(guān)函數(shù)R時(shí)，利用3次FFT即可實(shí)現(xiàn)，即： R = IFFT[FFT(st). *FFT(refg)]。
[0037]具體的，DSP的核0作為主核，收發(fā)雷達(dá)回波數(shù)據(jù)、控制核0-7運(yùn)行及實(shí)現(xiàn)積累互相 關(guān)包絡(luò)對(duì)齊計(jì)算。
[0038]進(jìn)一步的，核0-7之間的通信及協(xié)同工作是通過改變?cè)O(shè)在DSP的共享存儲(chǔ)區(qū)的一個(gè) 狀態(tài)標(biāo)志變量f lag[8]來實(shí)現(xiàn)的，該變量是int型的8維數(shù)組，分別對(duì)應(yīng)每個(gè)核的狀態(tài)標(biāo)志， 且每個(gè)核都能訪問和修改該變量，初始值都為〇。
[0039] 進(jìn)一步的，DSP核0獲得穩(wěn)定的全局基準(zhǔn)信號(hào)后，核0修改所有核的狀態(tài)標(biāo)志，使其 遞增1;當(dāng)每個(gè)分核檢測(cè)到本核的狀態(tài)標(biāo)志遞增1時(shí)，則利用全局基準(zhǔn)信號(hào)并行分別處理后 續(xù)8個(gè)PRF的積累互相關(guān)包絡(luò)對(duì)齊，并同時(shí)分別得到8個(gè)對(duì)齊后的回波和8個(gè)局部基準(zhǔn)信號(hào)， 然后使本核的狀態(tài)標(biāo)志遞增1;當(dāng)核〇檢測(cè)到8個(gè)核的狀態(tài)標(biāo)志都已增加1，則利用8個(gè)局部基 準(zhǔn)信號(hào)更新全局基準(zhǔn)信號(hào)，更新完后核〇再次將所有核的狀態(tài)標(biāo)志遞增1，重復(fù)步驟2。
[0040] 本發(fā)明的有益效果是:本發(fā)明充分利用8核DSP的高性能特點(diǎn)，將積累互相關(guān)包絡(luò) 對(duì)齊算法并行地在8核DSP上實(shí)現(xiàn)，達(dá)到提尚算法計(jì)算效率目的；同時(shí)，本發(fā)明將一幀待對(duì)齊 的雷達(dá)回波數(shù)據(jù)的前Z個(gè)脈沖回波先順序運(yùn)算獲得相對(duì)穩(wěn)定的基準(zhǔn)信號(hào)，這樣能保證整幀 回波的對(duì)齊效果;另一方面，其基準(zhǔn)信號(hào)不是單一的某個(gè)回波信號(hào)，而是加權(quán)積累實(shí)時(shí)更新 的，這樣能很大程度地消融異常信號(hào)在基準(zhǔn)信號(hào)中的權(quán)重，可有效消除目標(biāo)回波中起伏部 分對(duì)運(yùn)動(dòng)補(bǔ)償?shù)挠绊?；同時(shí)基準(zhǔn)信號(hào)中與待對(duì)齊回波越鄰近的信號(hào)占越高的權(quán)值，能獲得 更好的對(duì)齊效果。
【附圖說明】
[0041 ]圖1所示為本發(fā)名積累互相關(guān)包絡(luò)對(duì)齊的8核DSP片上并行實(shí)現(xiàn)方法實(shí)施例中8核 DSP(TMS320C6678)片上并行實(shí)現(xiàn)方法實(shí)現(xiàn)框圖；
[0042]圖2所示為本發(fā)名積累互相關(guān)包絡(luò)對(duì)齊的8核DSP片上并行實(shí)現(xiàn)方法實(shí)施例中DSP 上分配流程圖；
[0043]圖3所示為本發(fā)名積累互相關(guān)包絡(luò)對(duì)齊的8核DSP片上并行實(shí)現(xiàn)方法實(shí)施例中互相 關(guān)函數(shù)最大值二項(xiàng)式內(nèi)插法選取最大值點(diǎn)示意圖；
[0044]圖4為本發(fā)名積累互相關(guān)包絡(luò)對(duì)齊的8核DSP片上并行實(shí)現(xiàn)方法實(shí)施例中第一個(gè)脈 沖回波處理流程圖；
[0045]圖5所示為本發(fā)名積累互相關(guān)包絡(luò)對(duì)齊的8核DSP片上并行實(shí)現(xiàn)方法實(shí)施例中前16 個(gè)脈沖處理流程圖。
[0046] 以下結(jié)合實(shí)施例的【具體實(shí)施方式】，對(duì)本發(fā)明的上述內(nèi)容再作進(jìn)一步的詳細(xì)說明。 但不應(yīng)將此理解為本發(fā)明上述主題的范圍僅限于以下的實(shí)例。在不脫離本發(fā)明上述技術(shù)思 想情況下，根據(jù)本領(lǐng)域普通技術(shù)知識(shí)和慣用手段做出的各種替換或變更，均應(yīng)包括在本發(fā) 明的范圍內(nèi)。
【具體實(shí)施方式】
[0047] 下面結(jié)合附圖及實(shí)施例詳細(xì)描述本發(fā)明的技術(shù)方案：
[0048]本發(fā)明針對(duì)現(xiàn)有技術(shù)中能有效避免包絡(luò)漂移和突跳誤差的積累基準(zhǔn)信號(hào)互相關(guān) 包絡(luò)對(duì)齊的計(jì)算方法只能按脈沖回波順次運(yùn)算，不能并行計(jì)算提高效率的問題，提供一種 積累互相關(guān)包絡(luò)對(duì)齊的8核DSP片上并行實(shí)現(xiàn)方法，首先，DSP核0對(duì)一幀雷達(dá)回波數(shù)據(jù)的前Z 個(gè)PRF的脈沖回波作積累互相關(guān)包絡(luò)對(duì)齊處理以獲得穩(wěn)定的全局基準(zhǔn)信號(hào)ref g，為核0-7后 續(xù)對(duì)齊算法中共用，其中Z為8的倍數(shù);其次，DSP核0-7根據(jù)全局基準(zhǔn)信號(hào)refg并行分別處理 后續(xù)8個(gè)PRF的積累互相關(guān)包絡(luò)對(duì)齊，每個(gè)核分別得到對(duì)齊后的回波信號(hào)和局部基準(zhǔn)信號(hào) ^6，1 = [0、1、2、3、4、5、6、7]，代表當(dāng)前核序號(hào)，然后05?核0檢測(cè)一幀所有脈沖信號(hào)是否都 已經(jīng)對(duì)齊，是則輸出對(duì)齊后的回波信號(hào)，否則SP核0利用八核輸出的八個(gè)基準(zhǔn)信號(hào)ref【〇-7】更 新全局基準(zhǔn)信號(hào)，重復(fù)上述步驟。本發(fā)明充分利用8核DSP的高性能特點(diǎn)，將積累互相關(guān)包絡(luò) 對(duì)齊算法并行地在8核DSP上實(shí)現(xiàn)，達(dá)到提尚算法計(jì)算效率目的；同時(shí)，本發(fā)明將一幀待對(duì)齊 的雷達(dá)回波數(shù)據(jù)的前Z個(gè)脈沖回波先順序運(yùn)算獲得相對(duì)穩(wěn)定的基準(zhǔn)信號(hào)，這樣能保證整幀 回波的對(duì)齊效果;另一方面，其基準(zhǔn)信號(hào)不是單一的某個(gè)回波信號(hào)，而是加權(quán)積累實(shí)時(shí)更新 的，這樣能很大程度地消融異常信號(hào)在基準(zhǔn)信號(hào)中的權(quán)重，可有效消除目標(biāo)回波中起伏部 分對(duì)運(yùn)動(dòng)補(bǔ)償?shù)挠绊?；同時(shí)基準(zhǔn)信號(hào)中與待對(duì)齊回波越鄰近的信號(hào)占越高的權(quán)值，能獲得 更好的對(duì)齊效果。另外，在積累互相關(guān)包絡(luò)對(duì)齊中的插值處理和計(jì)算互相關(guān)函數(shù)本發(fā)明利 用快速傅里葉變換(FFT)提高了計(jì)算效率。
[0049] 實(shí)施例
[0050] 本例以TI公司在2010年11月推出了業(yè)界性能最高的多核DSP處理芯片 TMS320C6678為例，具體闡述本發(fā)明的
【發(fā)明內(nèi)容】
，為了方便描述本發(fā)明的內(nèi)容，首先作以下 技術(shù)術(shù)語定義：
[0051]定義1:逆合成孔徑雷達(dá)(ISAR)成像；
[0052] ISAR是依靠雷達(dá)與目標(biāo)間的相對(duì)運(yùn)動(dòng)(一般是雷達(dá)不動(dòng))提高橫向分辨率，利用大 帶寬提高距離分辨率，以此對(duì)一般雷達(dá)目標(biāo)，如飛機(jī)、艦船、導(dǎo)彈等進(jìn)行二維成像。
[0053]定義2:-幀數(shù)據(jù)；
[0054] 一幀數(shù)據(jù)是指成一幅ISAR像所需的若干脈沖。
[0055] 定義3:包絡(luò)對(duì)齊；
[0056]待成像目標(biāo)相對(duì)于雷達(dá)存在平動(dòng)分量和轉(zhuǎn)動(dòng)分量，平動(dòng)分量為目標(biāo)相對(duì)于雷達(dá)射 線的姿態(tài)保持不變，作平移運(yùn)動(dòng)。平動(dòng)分量會(huì)使目標(biāo)回波包絡(luò)平移，這對(duì)最后成像沒有貢 獻(xiàn)，需要將其效果補(bǔ)償?shù)?。平?dòng)補(bǔ)償又包括包絡(luò)對(duì)齊和初相校正，而包絡(luò)對(duì)齊就是將一幀數(shù) 據(jù)的所有脈沖的目標(biāo)回波實(shí)包絡(luò)逐個(gè)對(duì)齊。
[0057]定義4:共享存儲(chǔ)區(qū)；
[0058]本發(fā)明中的共享存儲(chǔ)區(qū)是C6678多核共享存儲(chǔ)區(qū)域，是指地址為OxOCOOOOOO到 0x0C3FFFFF共4M的存儲(chǔ)區(qū)域，可供8核訪問。
[0059] 本例在不影響算法性能的前提下，將順序執(zhí)行的積累基準(zhǔn)信號(hào)互相關(guān)包絡(luò)對(duì)齊算 法并行地在C6678上實(shí)現(xiàn)，大大提高其運(yùn)算時(shí)間，如圖1所示，C6678核0先對(duì)一幀雷達(dá)回波數(shù) 據(jù)的前前Z個(gè)，其中Z為8的倍數(shù);本例為16個(gè)(第0-15個(gè))PRF的脈沖回波作積累互相關(guān)包絡(luò) 對(duì)齊處理以獲得相對(duì)穩(wěn)定的基準(zhǔn)信號(hào)，然后核0修改所有核的狀態(tài)標(biāo)志，使其遞增1;第二 步，DSP的核0-7-旦檢測(cè)到本核標(biāo)志位已改變成預(yù)期值，則并行分別處理后續(xù)8個(gè)（如當(dāng) flag[0-8] = 1時(shí)，核0-7分別處理第16-23個(gè)脈沖回波)PRF的積累互相關(guān)包絡(luò)對(duì)齊，計(jì)算完 成后會(huì)在程序已設(shè)定地址輸出對(duì)齊后的回波信號(hào)和基準(zhǔn)信號(hào)，同時(shí)將相應(yīng)標(biāo)志位再遞增1; 第三步，核0檢測(cè)所有核的標(biāo)志位狀態(tài)，一旦都已增加1，則8核都已完成本次運(yùn)算，此時(shí)利用 這8個(gè)核輸出的基準(zhǔn)信號(hào)計(jì)算新的基準(zhǔn)信號(hào);循環(huán)執(zhí)行第二和第三步直至這一幀的所有脈 沖信號(hào)都對(duì)齊。最后，核〇將所有核的狀態(tài)標(biāo)志都置為〇,即flag[0_7]=0,使整個(gè)DSP可執(zhí)行 下一幀數(shù)據(jù)的包絡(luò)對(duì)齊。其中，核〇如圖2所示，將多核DSP的核0作為主核，執(zhí)行收發(fā)數(shù)據(jù)、控 制其它核以及計(jì)算的操作，而多核之間的通信及協(xié)同工作是通過改變?cè)O(shè)在DSP的共享存儲(chǔ) 區(qū)的一個(gè)狀態(tài)標(biāo)志變量flag[8]來實(shí)現(xiàn)的，該變量是int型的8維數(shù)組，分別對(duì)應(yīng)每個(gè)核的狀 態(tài)標(biāo)志，且每個(gè)核都能訪問和修改該變量，初始值都為0。充分利用8核DSP的高性能特點(diǎn)，將 積累互相關(guān)包絡(luò)對(duì)齊算法并行地在8核DSP上實(shí)現(xiàn)，達(dá)到提高算法計(jì)算效率目的；同時(shí)，本發(fā) 明將一幀待對(duì)齊的雷達(dá)回波數(shù)據(jù)的前16個(gè)脈沖回波先順序運(yùn)算獲得相對(duì)穩(wěn)定的基準(zhǔn)信號(hào)， 這樣能保證整幀回波的對(duì)齊效果;另一方面，其基準(zhǔn)信號(hào)不是單一的某個(gè)回波信號(hào)，而是加 權(quán)積累實(shí)時(shí)更新的，這樣能很大程度地消融異常信號(hào)在基準(zhǔn)信號(hào)中的權(quán)重，可有效消除目 標(biāo)回波中起伏部分對(duì)運(yùn)動(dòng)補(bǔ)償?shù)挠绊?；同時(shí)基準(zhǔn)信號(hào)中與待對(duì)齊回波越鄰近的信號(hào)占越高 的權(quán)值，能獲得更好的對(duì)齊效果。
[0060] 下面具體以一幀256個(gè)脈沖回波作為待對(duì)齊回波數(shù)據(jù)，具體說明本發(fā)明；
[0061] 首先，初始化參數(shù)及開辟共享存儲(chǔ)區(qū)變量。初始化參數(shù)有：一幀待處理數(shù)據(jù)有256 個(gè)脈沖回波，每個(gè)脈沖回波包含1024個(gè)采樣點(diǎn)。在C6678共享存儲(chǔ)區(qū)設(shè)定狀態(tài)標(biāo)志位數(shù)組， 并賦予初始值'〇'，即flag[0_7]=0。同時(shí)開辟9個(gè)長度相同的基準(zhǔn)信號(hào)空間，分別是包絡(luò)對(duì) 齊中用到得基準(zhǔn)信號(hào)和核0-7每次運(yùn)算后得到的基準(zhǔn)信號(hào)。另外，待包絡(luò)對(duì)齊的數(shù)據(jù)和包絡(luò) 對(duì)齊后的數(shù)據(jù)都放在C6678的外部存儲(chǔ)空間。
[0062]其次，對(duì)于一幀待處理數(shù)據(jù)的前16個(gè)脈沖回波的包絡(luò)對(duì)齊，只核0執(zhí)行;核0接收一 幀待對(duì)齊回波數(shù)據(jù)，并將接收到的第〇個(gè)脈沖回波作為第一基準(zhǔn)信號(hào)，從第1個(gè)脈沖回波開 始順序執(zhí)行積累互相關(guān)包絡(luò)對(duì)齊算法，直到前16個(gè)脈沖都對(duì)齊完，然后將所有核標(biāo)志位置 ' 1 '。其中，具體的積累互相關(guān)包絡(luò)對(duì)齊處理具體包括：
[0063]針對(duì)待處理的一幀N*M維的雷達(dá)回波數(shù)據(jù)，其中，N為回波脈沖個(gè)數(shù)，Μ為每個(gè)脈沖 所含的距離單元數(shù)；以第一個(gè)脈沖為初始的基準(zhǔn)信號(hào)ref，取第二個(gè)脈沖回波作為當(dāng)前回波 記為st，具體如圖3所示；
[0064] A1，核0計(jì)算基準(zhǔn)信號(hào)ref與當(dāng)前回波st的互相關(guān)函數(shù)R，并獲取互相關(guān)函數(shù)峰值 Rmax及其對(duì)應(yīng)的坐標(biāo)記做Xmax，進(jìn)入步驟B1;
[0065] B1，核0分別取出峰值點(diǎn)(Xmax，RMax)左右相鄰的互相關(guān)函數(shù)值及其對(duì)應(yīng)的坐標(biāo)分別 記為(Xi，Ri)及(X2，R2)，這三點(diǎn)可擬合出一條凹口向下的二次項(xiàng)曲線，解算出該曲線的最大 值點(diǎn)對(duì)應(yīng)的坐楊
進(jìn)入步驟C1;
[0066] C1，核0利用此時(shí)的坐標(biāo)X'構(gòu)造補(bǔ)償項(xiàng)并與當(dāng)前回波st點(diǎn)乘(對(duì)應(yīng)距離 單元相乘)e | e φ以獲得對(duì)齊后的當(dāng)前脈沖回波st '并存儲(chǔ)，進(jìn)入步驟D1;
[0067] D1，根據(jù)對(duì)齊后的當(dāng)前脈沖回波st'計(jì)算新的基準(zhǔn)信號(hào):ref = γ · ref+st'，其中， 更新權(quán)系數(shù)〇< γ <1;至此當(dāng)前脈沖回波的包絡(luò)對(duì)齊完成，進(jìn)入步驟El;
[0068] El，核0取出下一個(gè)脈沖回波，根據(jù)新的基準(zhǔn)信號(hào)，重復(fù)步驟A，直到前Z個(gè)回波數(shù)據(jù) 全部完成對(duì)齊，獲得穩(wěn)定的全局基準(zhǔn)信號(hào)refg。
[0069] 如圖4及圖5所示，在上述積累互相關(guān)包絡(luò)對(duì)齊算法的并行處理中，還包括利用快 速傅里葉變換(FFT)對(duì)基準(zhǔn)信號(hào)ref與當(dāng)前回波s t進(jìn)行插值處理：
[0070] 假設(shè)，包絡(luò)對(duì)齊的精度要求為1/8個(gè)距離單元，則需要對(duì)基準(zhǔn)信號(hào)ref與當(dāng)前回波 st進(jìn)行8倍插值，在DSP中可采用兩次不同點(diǎn)數(shù)的FFT來實(shí)現(xiàn)8倍插值，8卩：IFFT[FFT( st，M)， 8M]，其中，Μ為每個(gè)脈沖所含的距離單元數(shù)。同時(shí)，核0利用快速傅里葉變換(FFT)提高互相 關(guān)函數(shù)計(jì)算效率；
[0071] 計(jì)算基準(zhǔn)信號(hào)ref與當(dāng)前回波St的互相關(guān)函數(shù)R時(shí)，利用3次FFT即可實(shí)現(xiàn)，即：R = IFFT[FFT(st).*FFT(ref)]〇
[0072] 然后，建立循環(huán)體。循環(huán)體內(nèi)，設(shè)立循環(huán)次數(shù)為η(循環(huán)條件為:f〇r(n = 0;n〈30;n+ + ))，核〇-核7同時(shí)執(zhí)行：當(dāng)每個(gè)核（當(dāng)前核序號(hào)為CoreNum)檢測(cè)到當(dāng)前核的狀態(tài)標(biāo)志位為 flag[CoreNum] =2n+l時(shí)，就取出相應(yīng)的第（16+8n+CoreNum)個(gè)脈沖回波做積累互相關(guān)包絡(luò) 對(duì)齊，然后給出相應(yīng)的對(duì)齊后回波和當(dāng)前核基準(zhǔn)信號(hào)，同時(shí)更新其狀態(tài)標(biāo)志位：flag [CoreNum] =f lag[CoreNum] + l;其中，具體的積累互相關(guān)包絡(luò)對(duì)齊處理與前16個(gè)脈沖的計(jì) 算方式一致，具體包括以下步驟：
[0073] A2,核0-7分別計(jì)算全局基準(zhǔn)信號(hào)refg與當(dāng)前回波st的互相關(guān)函數(shù)R，并獲取互相關(guān) 函數(shù)峰值Rmax及其對(duì)應(yīng)的坐標(biāo)記做Xmax，進(jìn)入步驟B2;
[0074] B2,核0-7分別取出峰值點(diǎn)(Xmax，Rmax)左右相鄰的互相關(guān)函數(shù)值及其對(duì)應(yīng)的坐標(biāo)分 別記為(Xi，Ri)及(X2，R2)，這三點(diǎn)可擬合出一條凹口向下的二次項(xiàng)曲線，解算出該曲線的最 大值點(diǎn)對(duì)應(yīng)的坐相 進(jìn)入步驟C2;
[0075] C2,核0-7方別不」用見h、」tfJ鈀你Λ恂?仲償項(xiàng)：e-如如M1]，并與當(dāng)前回波st點(diǎn)乘（對(duì) 應(yīng)距離單元相乘）? 以獲得對(duì)齊后的當(dāng)前脈沖回波St'并存儲(chǔ)，以此八 ~hf ^ ^ 核分別得到八個(gè)對(duì)齊后的脈沖回波，并分別根據(jù)對(duì)齊后的當(dāng)前脈沖回波St'計(jì)算對(duì)應(yīng)核的 局部基準(zhǔn)信號(hào):ref I = St'；至此當(dāng)前脈沖回波的包絡(luò)對(duì)齊完成，進(jìn)入步驟D2;
[0076] D2,核0根據(jù)八核輸出的局部基準(zhǔn)信號(hào)更新全局基準(zhǔn)信號(hào)：
[0077] refg= γ 8refg+ γ 7ref7+ γ 6ref6+ γ 5ref5+ γ 4ref4+ γ 3ref3+ γ 2ref 2+ γ Vef 1 + ref0；0^i γ ^Ξ1 ；
[0078] 其中，refO為核0計(jì)算出的局部基準(zhǔn)信號(hào)，refl為核1計(jì)算出的局部基準(zhǔn)信號(hào)，ref2 為核2計(jì)算出的局部基準(zhǔn)信號(hào)，ref 3為核3計(jì)算出的局部基準(zhǔn)信號(hào)，ref 4為核4計(jì)算出的局部 基準(zhǔn)信號(hào)，ref5為核5計(jì)算出的局部基準(zhǔn)信號(hào)，ref6為核6計(jì)算出的局部基準(zhǔn)信號(hào)，ref7為核 7計(jì)算出的局部基準(zhǔn)信號(hào)；
[0079] 取出下一個(gè)脈沖回波，根據(jù)新的全局基準(zhǔn)信號(hào)，重復(fù)步驟A2,直到一幀回波數(shù)據(jù)全 部完成對(duì)齊，輸出對(duì)齊后的回波信號(hào)。
[0080] 而在積累互相關(guān)包絡(luò)對(duì)齊算法的并行處理中，還包括利用傅里葉變化對(duì)全局基準(zhǔn) 信號(hào)refg與當(dāng)前回波St進(jìn)行插值處理：
[0081] 假設(shè)，包絡(luò)對(duì)齊的精度要求為1/8個(gè)距離單元，則需要對(duì)全局基準(zhǔn)信號(hào)refg與當(dāng)前 回波st進(jìn)行8倍插值，在DSP中可采用兩次不同點(diǎn)數(shù)的FFT來實(shí)現(xiàn)8倍插值，8卩：IFFT[FFT( st， Μ)，8M]，其中，Μ為每個(gè)脈沖所含的距離單元數(shù)。首先，包絡(luò)對(duì)齊精度要求達(dá)到1/8個(gè)距離單 元，所以需要對(duì)兩個(gè)待對(duì)齊的信號(hào)進(jìn)行8倍插值，這在DSP上實(shí)現(xiàn)時(shí)是將長度為Ν的待插值信 號(hào)Ν點(diǎn)FFT變換到頻域，再在其兩端分別補(bǔ)3.5Ν個(gè)零，然后8Ν點(diǎn)IFFT變回原來時(shí)域。其次，利 用時(shí)域卷積相當(dāng)于頻域點(diǎn)乘的理論，待對(duì)齊脈沖信號(hào)與基準(zhǔn)信號(hào)的時(shí)域互相關(guān)是通過將兩 信號(hào)FFT變換到頻域，再將兩者點(diǎn)乘結(jié)果IFFT變換回時(shí)域?qū)崿F(xiàn)。
[0082]同時(shí)，核0-7也利用快速傅里葉變換(FFT)提高互相關(guān)函數(shù)計(jì)算效率；
[0083]計(jì)算全局基準(zhǔn)信號(hào)refg與當(dāng)前回波st的互相關(guān)函數(shù)R時(shí)，利用3次FFT即可實(shí)現(xiàn)，即： R = IFFT[FFT(st). *FFT(refg)]。
[0084] 在上述循環(huán)體內(nèi)，只核0執(zhí)行。核0判斷若8個(gè)核的標(biāo)志位值都為2n+2時(shí)，則根據(jù)上 述計(jì)算公式更新當(dāng)前基準(zhǔn)信號(hào)。然后將所有核的標(biāo)志位增加'1'。
[0085] 最后，當(dāng)核0檢測(cè)到一幀數(shù)據(jù)全部已經(jīng)對(duì)齊后，輸出對(duì)其后的回波信號(hào)，跳出循環(huán) 體，只核〇執(zhí)行:核〇發(fā)送數(shù)據(jù)并將狀態(tài)標(biāo)志位還原成初始值' 〇 '以備下一幀數(shù)據(jù)的計(jì)算。 [0086] 經(jīng)試驗(yàn)證實(shí)，積累互相關(guān)包絡(luò)對(duì)齊算法在TS201上根本無法實(shí)現(xiàn)，若在C6678(主頻 為1GHz時(shí))上順序計(jì)算256個(gè)脈沖回波(每個(gè)脈沖有1024個(gè)采樣點(diǎn)）的積累互相關(guān)包絡(luò)對(duì)齊 算法則需要1268. lms，而按照本發(fā)明的方法則只需要61.4ms，大大提高了計(jì)算效率。
【主權(quán)項(xiàng)】
1. 積累互相關(guān)包絡(luò)對(duì)齊的8核DSP片上并行實(shí)現(xiàn)方法，其特征在于，包括w下步驟： 步驟1、DSP核0對(duì)一帖雷達(dá)回波數(shù)據(jù)的前Z個(gè)PRF的脈沖回波作積累互相關(guān)包絡(luò)對(duì)齊處 理W獲得穩(wěn)定的全局基準(zhǔn)信號(hào)refg，為核0-7后續(xù)對(duì)齊算法中共用，其中Z為8的倍數(shù)； 步驟2、DSP核0-7根據(jù)全局基準(zhǔn)信號(hào)refg并行分別處理后續(xù)8個(gè)PRF的積累互相關(guān)包絡(luò)對(duì) 齊，每個(gè)核分別得到對(duì)齊后的回波信號(hào)和局部基準(zhǔn)信號(hào)refI，I = [0、l、2、3、4、5、6、7]，代表 當(dāng)前核序號(hào)，進(jìn)入步驟3; 步驟3、DSP核0檢測(cè)一帖所有脈沖信號(hào)是否都已經(jīng)對(duì)齊，是則輸出對(duì)齊后的回波信號(hào)， 否則進(jìn)入步驟4; 步驟4，DSP核0利用八核輸出的八個(gè)基準(zhǔn)信號(hào)ref【〇-7】更新全局基準(zhǔn)信號(hào)，重復(fù)步驟2。2. 根據(jù)權(quán)利要求1所述的積累互相關(guān)包絡(luò)對(duì)齊的8核DSP片上并行實(shí)現(xiàn)方法，其特征在 于:所述步驟1中，DSP核0對(duì)一帖雷達(dá)回波數(shù)據(jù)的前Z個(gè)PRF的脈沖回波作積累互相關(guān)包絡(luò)對(duì) 齊處理W獲得穩(wěn)定的全局基準(zhǔn)信號(hào)refg,包括W下步驟： 針對(duì)待處理的一帖N*M維的雷達(dá)回波數(shù)據(jù)，其中，N為回波脈沖個(gè)數(shù)，Μ為每個(gè)脈沖所含 的距離單元數(shù)；W第一個(gè)脈沖為初始的基準(zhǔn)信號(hào)ref,取第二個(gè)脈沖回波作為當(dāng)前回波記為 St； A1，核0計(jì)算基準(zhǔn)信號(hào)ref與當(dāng)前回波St的互相關(guān)函數(shù)R，并獲取互相關(guān)函數(shù)峰值Rmax及其 對(duì)應(yīng)的坐標(biāo)記做Xmax，進(jìn)入步驟B1; B1，核0分別取出峰值點(diǎn)(Xmax，Rmax)左右相鄰的互相關(guān)函數(shù)值及其對(duì)應(yīng)的坐標(biāo)分別記為 (Xi，Ri)及(X2，R2)，運(yùn)Ξ點(diǎn)可擬合出一條凹口向下的二次項(xiàng)曲線，解算出該曲線的最大值點(diǎn) 對(duì)應(yīng)的坐標(biāo)進(jìn)入步驟C1; C1，核0利用此時(shí)的坐標(biāo)X'構(gòu)造補(bǔ)償項(xiàng)：?？趬?。并與當(dāng)前回波St點(diǎn)乘(對(duì)應(yīng)距離單元 相乘：^獲得對(duì)齊后的當(dāng)前脈沖回波St'并存儲(chǔ)，進(jìn)入步驟D1; D1，根據(jù)對(duì)齊后的當(dāng)前脈沖回波St'計(jì)算新的基準(zhǔn)信號(hào)：ref= 丫 -ref+st'，其中，更新 權(quán)系數(shù)0《丫《1;至此當(dāng)前脈沖回波的包絡(luò)對(duì)齊完成，進(jìn)入步驟E1; E1，核0取出下一個(gè)脈沖回波，根據(jù)新的基準(zhǔn)信號(hào)，重復(fù)步驟A，直到前Z個(gè)回波數(shù)據(jù)全部 完成對(duì)齊，獲得穩(wěn)定的全局基準(zhǔn)信號(hào)refg。3. 根據(jù)權(quán)利要求2所述的積累互相關(guān)包絡(luò)對(duì)齊的8核DSP片上并行實(shí)現(xiàn)方法，其特征在 于，在積累互相關(guān)包絡(luò)對(duì)齊算法的并行處理中，還包括利用快速傅里葉變換(FFT)對(duì)基準(zhǔn)信 號(hào)ref與當(dāng)前回波St進(jìn)行插值處理： 假設(shè)，包絡(luò)對(duì)齊的精度要求為1/8個(gè)距離單元，則需要對(duì)基準(zhǔn)信號(hào)ref與當(dāng)前回波St進(jìn)行 8倍插值，在DSP中可采用兩次不同點(diǎn)數(shù)的FFT來實(shí)現(xiàn)8倍插值，即：IFFT[FFT(St，M)，8M]，其 中，Μ為每個(gè)脈沖所含的距離單元數(shù)。4. 根據(jù)權(quán)利要求2所述的積累互相關(guān)包絡(luò)對(duì)齊的8核DSP片上并行實(shí)現(xiàn)方法，其特征在 于，核0利用快速傅里葉變換(FFT)提高互相關(guān)函數(shù)計(jì)算效率； 計(jì)算基準(zhǔn)信號(hào)ref與當(dāng)前回波St的互相關(guān)函數(shù)則寸，利用3次FFT即可實(shí)現(xiàn)，即：R=IFFT [FFT(st).*FFT(ref)]〇5. 根據(jù)權(quán)利要求1所述的積累互相關(guān)包絡(luò)對(duì)齊的8核DSP片上并行實(shí)現(xiàn)方法，其特征在 于:所述步驟2-4中，包括W下步驟： A2,核0-7分別計(jì)算全局基準(zhǔn)信號(hào)refg與當(dāng)前回波St的互相關(guān)函數(shù)R，并獲取互相關(guān)函數(shù) 峰值Rmax及其對(duì)應(yīng)的坐標(biāo)記做Xmax，進(jìn)入步驟B2; B2，核0-7分別取出峰值點(diǎn)(Xmax，Rmax)左右相鄰的互相關(guān)函數(shù)值及其對(duì)應(yīng)的坐標(biāo)分別記 為(Xi，Ri)及(X2，R2)，運(yùn)Ξ點(diǎn)可擬合出一條凹口向下的二次項(xiàng)曲線，解算出該曲線的最大值 點(diǎn)對(duì)應(yīng)的坐標(biāo)進(jìn)入步驟C2; C2,核0-7分別利用此時(shí)的坐標(biāo)X'構(gòu)造補(bǔ)償項(xiàng):。-口請(qǐng)。并與當(dāng)前回波St點(diǎn)乘(對(duì)應(yīng)距 離單元相乘）W獲得對(duì)齊后的當(dāng)前脈沖回波St'并存儲(chǔ)，W此八核分 別得到八個(gè)對(duì)齊后的脈沖回波，并分別根據(jù)對(duì)齊后的當(dāng)前脈沖回波St'計(jì)算對(duì)應(yīng)核的局部 基準(zhǔn)信號(hào):ref 1 = St'；至此當(dāng)前脈沖回波的包絡(luò)對(duì)齊完成，進(jìn)入步驟D2; D2，核0根據(jù)八核輸出的局部基準(zhǔn)信號(hào)更新全局基準(zhǔn)信號(hào)： refg= 丫 8refg+ 丫 7ref7+ 丫 6ref6+ 丫 5fe巧+ 丫 4ref4+ 丫 3ref3+ 丫 2ref化丫 Vef 1+ref0; 0 《丫《1; 其中，refO為核0計(jì)算出的局部基準(zhǔn)信號(hào)，ref!為核1計(jì)算出的局部基準(zhǔn)信號(hào)，ref2為核 2計(jì)算出的局部基準(zhǔn)信號(hào)，ref3為核3計(jì)算出的局部基準(zhǔn)信號(hào)，ref4為核4計(jì)算出的局部基準(zhǔn) 信號(hào)，re巧為核5計(jì)算出的局部基準(zhǔn)信號(hào)，ref6為核6計(jì)算出的局部基準(zhǔn)信號(hào)，ref7為核7計(jì) 算出的局部基準(zhǔn)信號(hào)； 取出下一個(gè)脈沖回波，根據(jù)新的全局基準(zhǔn)信號(hào)，重復(fù)步驟A2,直到一帖回波數(shù)據(jù)全部完 成對(duì)齊，輸出對(duì)齊后的回波信號(hào)。6. 根據(jù)權(quán)利要求5所述的積累互相關(guān)包絡(luò)對(duì)齊的8核DSP片上并行實(shí)現(xiàn)方法，其特征在 于，在積累互相關(guān)包絡(luò)對(duì)齊算法的并行處理中，還包括利用傅里葉變化對(duì)全局基準(zhǔn)信號(hào) refg與當(dāng)前回波St進(jìn)行插值處理： 假設(shè)，包絡(luò)對(duì)齊的精度要求為1/8個(gè)距離單元，則需要對(duì)全局基準(zhǔn)信號(hào)refg與當(dāng)前回波 St進(jìn)行8倍插值，在DSP中可采用兩次不同點(diǎn)數(shù)的FFT來實(shí)現(xiàn)8倍插值，即：IFFT[FFT( St，M)， 8M]，其中，Μ為每個(gè)脈沖所含的距離單元數(shù)。7. 根據(jù)權(quán)利要求5所述的積累互相關(guān)包絡(luò)對(duì)齊的8核DSP片上并行實(shí)現(xiàn)方法，其特征在 于，核0利用FFT提高互相關(guān)函數(shù)計(jì)算效率； 計(jì)算全局基準(zhǔn)信號(hào)refg與當(dāng)前回波St的互相關(guān)函數(shù)則寸，利用3次FFT即可實(shí)現(xiàn)，即：R = IFFT[FFT(st).沖FT(refg)]。8. 根據(jù)權(quán)利要求1所述的積累互相關(guān)包絡(luò)對(duì)齊的8核DSP片上并行實(shí)現(xiàn)方法，其特征在 于，DSP的核0作為主核，收發(fā)雷達(dá)回波數(shù)據(jù)、控制核0-7運(yùn)行及實(shí)現(xiàn)積累互相關(guān)包絡(luò)對(duì)齊計(jì) 算。9. 根據(jù)權(quán)利要求8所述的積累互相關(guān)包絡(luò)對(duì)齊的8核DSP片上并行實(shí)現(xiàn)方法，其特征在 于，核0-7之間的通信及協(xié)同工作是通過改變?cè)O(shè)在DSP的共享存儲(chǔ)區(qū)的一個(gè)狀態(tài)標(biāo)志變量 flag[引來實(shí)現(xiàn)的，該變量是int型的8維數(shù)組，分別對(duì)應(yīng)每個(gè)核的狀態(tài)標(biāo)志，且每個(gè)核都能 訪問和修改該變量，初始值都為ο。10.根據(jù)權(quán)利要求9所述的積累互相關(guān)包絡(luò)對(duì)齊的8核DSP片上并行實(shí)現(xiàn)方法，其特征在 于，DSP核0獲得穩(wěn)定的全局基準(zhǔn)信號(hào)后，核0修改所有核的狀態(tài)標(biāo)志，使其遞增1;當(dāng)每個(gè)分 核檢測(cè)到本核的狀態(tài)標(biāo)志遞增1時(shí)，則利用全局基準(zhǔn)信號(hào)并行分別處理后續(xù)8個(gè)PRF的積累 互相關(guān)包絡(luò)對(duì)齊，并同時(shí)分別得到8個(gè)對(duì)齊后的回波和8個(gè)局部基準(zhǔn)信號(hào)，然后使本核的狀 態(tài)標(biāo)志遞增1;當(dāng)核0檢測(cè)到8個(gè)核的狀態(tài)標(biāo)志都已增加1，則利用8個(gè)局部基準(zhǔn)信號(hào)更新全局 基準(zhǔn)信號(hào)，更新完后核0再次將所有核的狀態(tài)標(biāo)志遞增1，重復(fù)步驟2。
【文檔編號(hào)】G01S13/90GK105974416SQ201610597721
【公開日】2016年9月28日
【申請(qǐng)日】2016年7月26日
【發(fā)明人】帥曉飛, 朱玉軍, 張桂梅, 郭俊良, 金力, 徐建, 張之峰, 劉永強(qiáng)
【申請(qǐng)人】四川電子軍工集團(tuán)裝備技術(shù)有限公司