[發(fā)明專利]一種基于GPU的長碼分段重疊局部相關(guān)捕獲方法

有效

• 申請?zhí)枺?/td>: 201611006729.5
• 申請日: 2016-11-16
• 公開（公告）號: CN106772471B
• 公開（公告）日: 2019-03-22
• 發(fā)明（設(shè)計）人: 牟衛(wèi)華;唐小妹;馬春江;歐鋼;王飛雪;孫廣富;倪少杰;高毅;李壯
• 申請（專利權(quán)）人: 中國人民解放軍國防科學(xué)技術(shù)大學(xué)
• 主分類號: G01S19/29
• 分類號: G01S19/29;G01S19/30
• 代理公司: 北京中濟(jì)緯天專利代理有限公司 11429
• 代理人: 陳立新

1.一種基于GPU的長碼分段重疊局部相關(guān)捕獲方法，其特征在于，包括以下步驟：

步驟S1提取長碼生成本地參考信號序列

根據(jù)本地時間t₀，取[t₀t₀+t]時間段持續(xù)時間為t的長碼數(shù)據(jù)，根據(jù)采樣率f_s生成本地參考信號序列p(n)，p(n)為{p(0) p(1) … p(tf_s)}；

步驟S2輸入信號剝離載波

以采樣率f_s對輸入信號進(jìn)行采樣，采集(T+t)秒的輸入信號，其中T為預(yù)設(shè)的偽碼搜索不確定時間范圍，剝離載波后生成輸入信號復(fù)信號序列s(n)，s(n)為{s(0) s(1) … s[(T+t)f_s]}；

步驟S3分配數(shù)據(jù)存儲空間，傳輸數(shù)據(jù)到GPU設(shè)備端

在GPU設(shè)備端為輸入信號復(fù)信號序列s(n)和本地參考信號序列p(n)分配存儲空間，使用cudaMemcpy()函數(shù)將輸入信號復(fù)信號序列s(n)和本地參考信號序列p(n)由主機(jī)端拷貝到GPU設(shè)備端；

步驟S4輸入信號復(fù)信號序列重疊分段劃分，本地信號序列補(bǔ)零擴(kuò)展

將輸入信號復(fù)信號序列部分重疊分段劃分，得到m個序列s_i(n)i＝0,1,…,(m-1)，將本地信號序列p(n)補(bǔ)充Tf_s/m個0，得到序列p0m(n)；

步驟S5創(chuàng)建傅里葉變換計劃

調(diào)用函數(shù)cufftPlan1d()、cufftPlanMany()創(chuàng)建傅里葉變換計劃；

步驟S6批處理計算信號序列傅里葉變換

調(diào)用函數(shù)cufftExecC2C()、cufftExecR2C()，批處理計算分段劃分的m個序列s_i(n)的傅里葉變換序列S_i(k)和本地信號序列p_0m(n)的傅里葉變換序列P_0m(k)；

步驟S7頻域共軛乘運(yùn)算

在GPU中逐點(diǎn)計算S_i(k)與P_0m(k)的共軛乘積，得到序列Y_im(k)；

步驟S8逆FFT得到多段相關(guān)值序列，拼接后得到整個相關(guān)值序列

調(diào)用函數(shù)cufftExecC2C()計算Y_im(k)的逆FFT，得到m個循環(huán)相關(guān)序列y_im(n)，截取每個相關(guān)序列的前Tf_s/m個值y_i(n)，拼接后得到整個不確定時間范圍內(nèi)的相關(guān)值序列r(n)；

步驟S9計算相關(guān)值序列模值

在GPU中逐點(diǎn)計算相關(guān)值序列r(n)的模值，得到模值序列M(n)；

步驟S10信號檢測

搜索M(n)的最大值，判斷是否超過檢測門限，并確定碼相位值t_d。