積累互相關(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ù)有限公司