本發明屬于高(gao)速(su)數(shu)傳，具體(ti)涉及一種關于低信噪比下(xia)非數(shu)據(ju)輔(fu)助的(de)高(gao)精度頻(pin)偏估計(ji)補償方(fang)法，用于寬帶衛星通信系(xi)統(tong)以及衛星高(gao)速(su)數(shu)傳系(xi)統(tong)。

背景技術：

1、進入21世紀以(yi)來(lai)，隨著寬帶互聯(lian)網(wang)大(da)面(mian)積普及(ji)、通(tong)信(xin)技術的(de)不(bu)斷革新換代以(yi)及(ji)衛(wei)星(xing)有效載(zai)荷(he)技術、多(duo)媒體技術等(deng)發(fa)展趨于成熟(shu)，寬帶衛(wei)星(xing)通(tong)信(xin)系統(tong)(tong)及(ji)高速數傳(chuan)(chuan)系統(tong)(tong)已(yi)成為如今一個重(zhong)要(yao)的(de)研究發(fa)展方向(xiang)。歐洲電信(xin)標(biao)(biao)準(zhun)化組織最新提出的(de)dvb-s2x標(biao)(biao)準(zhun)較之前的(de)dvb-s2標(biao)(biao)準(zhun)新增了(le)多(duo)種使用場景(jing)，如以(yi)vlsnr為特征的(de)甚小孔徑(jing)終端(vsat)、以(yi)超幀跳波束為特征的(de)ip語音傳(chuan)(chuan)輸(shu)、蜂窩回傳(chuan)(chuan)網(wang)絡(luo)等(deng)。此外，面(mian)對應(ying)用場景(jing)的(de)復雜多(duo)樣(yang)化，dvb-s2x可以(yi)結(jie)合自適應(ying)調(diao)制編(bian)碼技術(acm)自動選擇并(bing)提供最優(you)化的(de)調(diao)制解調(diao)方案，進一步優(you)化衛(wei)星(xing)鏈路，獲得(de)較大(da)的(de)效率提升，盡可能逼近理(li)論(lun)香(xiang)農限。

2、全(quan)球(qiu)衛星通(tong)信(xin)(xin)業(ye)務需求質量的(de)不斷提升以及高通(tong)量、高速(su)(su)率(lv)的(de)信(xin)(xin)息傳輸(shu)(shu)需求，使(shi)得可用(yong)空間頻(pin)譜資(zi)源(yuan)已經日(ri)益緊(jin)缺。當(dang)前(qian)針對較(jiao)高的(de)數(shu)據傳輸(shu)(shu)速(su)(su)率(lv)需求，更(geng)高頻(pin)譜利用(yong)率(lv)的(de)高階調(diao)制方式如(ru)16apsk、32apsk、最高為256apsk等已投入使(shi)用(yong)，更(geng)加高效(xiao)的(de)編譯碼方式如(ru)bch碼、ldpc碼也被采(cai)用(yong)，這使(shi)得解調(diao)信(xin)(xin)噪比門限要求得到進(jin)一步提高。這些因素使(shi)得如(ru)何在惡劣信(xin)(xin)噪比環境下對信(xin)(xin)號的(de)頻(pin)率(lv)偏移(yi)以及相位偏移(yi)進(jin)行高精度估計補(bu)償成為一種具有較(jiao)高價(jia)值的(de)研究(jiu)方向，否則將會(hui)成為信(xin)(xin)號接收解調(diao)的(de)瓶頸，嚴重(zhong)影(ying)響解調(diao)系統后續其它(ta)模塊的(de)正(zheng)常運行。

3、對于載波(bo)恢復算法(fa)的(de)研究，按照(zhao)是否使(shi)用(yong)導(dao)頻數(shu)據(ju)(ju)以(yi)及編譯碼(ma)(ma)信(xin)息(xi)，還可將(jiang)其分為以(yi)下三種(zhong)算法(fa)：一種(zhong)為利用(yong)當前數(shu)據(ju)(ju)幀中的(de)導(dao)頻信(xin)息(xi)進(jin)行(xing)相關(guan)(guan)計(ji)算進(jin)而進(jin)行(xing)載波(bo)同步(bu)的(de)數(shu)據(ju)(ju)輔助(zhu)方(fang)法(fa)(data?aided,da)，另一種(zhong)則為未直接(jie)使(shi)用(yong)到已知先驗數(shu)據(ju)(ju)的(de)非(fei)數(shu)據(ju)(ju)輔助(zhu)方(fang)法(fa)(nondata?aided,nda)，第三種(zhong)方(fang)法(fa)一般(ban)通過使(shi)用(yong)譯碼(ma)(ma)軟信(xin)息(xi)或者利用(yong)相關(guan)(guan)調制信(xin)息(xi)進(jin)行(xing)相關(guan)(guan)參數(shu)估計(ji)的(de)編碼(ma)(ma)輔助(zhu)算法(fa)(code?aided，ca)。

4、數據輔(fu)助(zhu)(zhu)(da)類算(suan)(suan)(suan)法(fa)(fa)主(zhu)要使(shi)用(yong)(yong)到(dao)數據幀(zhen)(zhen)中導頻塊等(deng)(deng)已(yi)知信(xin)息(xi)(xi)，如l&r算(suan)(suan)(suan)法(fa)(fa)、m&m算(suan)(suan)(suan)法(fa)(fa)、fitz算(suan)(suan)(suan)法(fa)(fa)、kay算(suan)(suan)(suan)法(fa)(fa)、l&w算(suan)(suan)(suan)法(fa)(fa)等(deng)(deng)，其(qi)中l&r算(suan)(suan)(suan)法(fa)(fa)與fitz算(suan)(suan)(suan)法(fa)(fa)雖然具(ju)有較(jiao)好(hao)的(de)(de)(de)(de)抗噪性(xing)(xing)能(neng)(neng)，但估(gu)偏范(fan)圍(wei)小，估(gu)計精度較(jiao)差，而(er)kay、l&w算(suan)(suan)(suan)法(fa)(fa)也僅(jin)是擴(kuo)大(da)了估(gu)計范(fan)圍(wei)，對(dui)于抗噪聲能(neng)(neng)力及估(gu)計精度方面還是較(jiao)差。并且該類算(suan)(suan)(suan)法(fa)(fa)需要使(shi)用(yong)(yong)到(dao)大(da)量導頻塊數據，會(hui)(hui)導致數據幀(zhen)(zhen)傳輸(shu)容量產(chan)生一定程度的(de)(de)(de)(de)損失。編碼輔(fu)助(zhu)(zhu)類(ca)算(suan)(suan)(suan)法(fa)(fa)最早(zao)由wangrok.oh等(deng)(deng)人于2001年提出，其(qi)本質上是將譯碼器輸(shu)出的(de)(de)(de)(de)譯碼軟信(xin)息(xi)(xi)與迭代性(xing)(xing)能(neng)(neng)和(he)同步相結合，來對(dui)載波相偏頻偏估(gu)計進(jin)行(xing)輔(fu)助(zhu)(zhu)。該算(suan)(suan)(suan)法(fa)(fa)在低信(xin)噪比環境(jing)下具(ju)有較(jiao)好(hao)的(de)(de)(de)(de)性(xing)(xing)能(neng)(neng)，但其(qi)存(cun)在著(zhu)明顯的(de)(de)(de)(de)弊端，其(qi)需要借助(zhu)(zhu)多次譯碼迭代，進(jin)而(er)帶來一定時延(yan)處理問(wen)題，影響整個系統的(de)(de)(de)(de)運行(xing)效(xiao)率并且計算(suan)(suan)(suan)復雜(za)度較(jiao)大(da)，而(er)繁(fan)瑣的(de)(de)(de)(de)迭代結構會(hui)(hui)進(jin)一步增加(jia)fpga的(de)(de)(de)(de)實現復雜(za)度。

5、非數(shu)據輔助方法是在(zai)(zai)基于最大似然理論的(de)基礎上(shang)，利用離散時間觀測(ce)量去估計單音信(xin)號的(de)頻(pin)(pin)率信(xin)息(xi)。dd算法作為(wei)(wei)一(yi)種較為(wei)(wei)經典的(de)算法，其(qi)以鎖(suo)相環為(wei)(wei)核(he)心，通過(guo)判(pan)決(jue)、鑒相以及環路濾波等步驟相結合對(dui)數(shu)據頻(pin)(pin)偏(pian)進(jin)(jin)行(xing)估計補償，但其(qi)對(dui)部分(fen)(fen)高(gao)(gao)階調制(zhi)方式，鑒相會(hui)(hui)較為(wei)(wei)困(kun)難，并且整(zheng)個環路需(xu)要一(yi)定(ding)(ding)數(shu)據去訓練收斂，較低的(de)信(xin)噪(zao)比(bi)(bi)會(hui)(hui)導(dao)致判(pan)決(jue)錯(cuo)誤，進(jin)(jin)而導(dao)致失鎖(suo)。在(zai)(zai)低信(xin)噪(zao)比(bi)(bi)場景(jing)下(xia)(xia)，傳統(tong)的(de)頻(pin)(pin)偏(pian)估計算法很難應用，尤其(qi)是面對(dui)無先驗數(shu)據以及多(duo)種調制(zhi)的(de)情況下(xia)(xia)。本發(fa)明在(zai)(zai)場景(jing)下(xia)(xia)具有一(yi)定(ding)(ding)優勢，可在(zai)(zai)esn0分(fen)(fen)別為(wei)(wei)-3db與6db的(de)情況下(xia)(xia)對(dui)opsk數(shu)據、8psk/8apsk數(shu)據進(jin)(jin)行(xing)高(gao)(gao)精(jing)度且穩定(ding)(ding)的(de)頻(pin)(pin)偏(pian)估計，并對(dui)其(qi)進(jin)(jin)行(xing)補償，具備較低的(de)計算復雜度，易(yi)于實現。

技術實現思路

1、本發明(ming)所要解決的技術問(wen)題是：

2、在(zai)(zai)低信噪(zao)比條件下(xia)，為(wei)了克(ke)服隨機(ji)噪(zao)聲的影(ying)響(xiang)，傳統(tong)的頻偏(pian)(pian)估(gu)計算(suan)(suan)法一般利用(yong)數(shu)據輔助類(da)算(suan)(suan)法通(tong)過導(dao)頻數(shu)據提取載(zai)波相(xiang)關信息，但該(gai)類算(suan)(suan)法估(gu)計精度易(yi)受信噪(zao)比影(ying)響(xiang)，且在(zai)(zai)衛星通(tong)信系(xi)統(tong)以及高(gao)速數(shu)傳系(xi)統(tong)中存在(zai)(zai)無導(dao)頻模式的數(shu)據幀。在(zai)(zai)這(zhe)些綜合場景下(xia)，傳統(tong)的頻偏(pian)(pian)估(gu)計算(suan)(suan)法難以直接應用(yong)，

3、針對在低信噪比(bi)條件下(xia)上訴頻(pin)偏估(gu)計(ji)方法存在的問題(ti)，本發明提供(gong)一種低信噪比(bi)下(xia)非數(shu)據(ju)輔助的高精(jing)度頻(pin)偏估(gu)計(ji)補(bu)償(chang)方法。在無需先驗數(shu)據(ju)的前提下(xia)，可(ke)通過粗(cu)估(gu)計(ji)與(yu)細(xi)估(gu)計(ji)兩級頻(pin)偏估(gu)計(ji)模塊獲得精(jing)度較高的估(gu)計(ji)值(zhi)，并對數(shu)據(ju)進行補(bu)償(chang)。

4、為了解決上述(shu)技術問(wen)題，本(ben)發明采用的技術方(fang)案為：

5、一種低信(xin)噪比(bi)下(xia)非數據輔助的高精度(du)頻偏估計補償方法，其特征在(zai)于，包括：

6、s1：對輸入的信(xin)號進(jin)行降階操(cao)作獲得四相星座點信(xin)號；

7、s2：將(jiang)s1得到(dao)的信(xin)號由直(zhi)角(jiao)坐(zuo)標系轉(zhuan)到(dao)極坐(zuo)標系，并(bing)將(jiang)其相位(wei)進行4倍放(fang)大，隨后重新映射(she)到(dao)直(zhi)角(jiao)坐(zuo)標系，得到(dao)去(qu)除調制信(xin)息(xi)的信(xin)號；

8、s3：對(dui)s2得(de)到(dao)的信號(hao)進行(xing)分段fft計算(suan)，得(de)到(dao)信號(hao)疊加頻譜并(bing)進行(xing)峰(feng)值搜索(suo)，獲得(de)周期圖峰(feng)值索(suo)引(yin)；

9、s4：根(gen)據s3得(de)到的頻譜(pu)峰值索(suo)引(yin)以及fft點數、調制階數獲得(de)頻偏粗(cu)估計(ji)值，完(wan)成(cheng)粗(cu)頻偏補(bu)償；

10、s5：對s4得到的(de)(de)補償(chang)后的(de)(de)信(xin)號(hao)進(jin)行s1～s2的(de)(de)降階去調制操作，并(bing)進(jin)行抽(chou)取因子為s的(de)(de)重(zhong)采樣，得到采樣率(lv)為fs/s的(de)(de)信(xin)號(hao)；

11、s6：對s5得到的重(zhong)采樣信(xin)號進(jin)行(xing)分段fft計(ji)算，依據信(xin)號疊(die)加(jia)頻譜獲得其峰(feng)值索引，之后進(jin)行(xing)插值計(ji)算獲得索引補償因子，逼近峰(feng)值實際位置；

12、s7：根據補償后(hou)的(de)峰值索引、抽取因子s以及fft點數(shu)參數(shu)進行細頻(pin)偏估計并完(wan)成對信號的(de)補償。

13、本發明(ming)進一步(bu)的技術方案(an)：所述s1之(zhi)前(qian)還包(bao)括：

14、s0：將接(jie)收(shou)數(shu)據符(fu)號(hao)經過定(ding)時(shi)恢復(fu)以及幀同步處理，消除符(fu)號(hao)定(ding)時(shi)偏得(de)到(dao)待糾偏數(shu)據，同時(shi)獲得(de)相應的調制(zhi)方(fang)式(shi)信息，所述調制(zhi)方(fang)式(shi)包括(kuo)非(fei)qpsk。

15、本發(fa)明進(jin)一步的技(ji)術方案：若調制(zhi)方式非qpsk，所述s1具(ju)體(ti)為：

16、對輸入信號的外環星座(zuo)點做(zuo)m次方降階操作，獲得與qpsk調制一致(zhi)的四相星座(zuo)點。

17、本發明進一步的(de)技術方案：對輸入信號的(de)外環星座(zuo)點做(zuo)m次(ci)方降(jiang)階(jie)操(cao)作(zuo)，獲(huo)得與qpsk調(diao)制一致的(de)四相星座(zuo)點，包括：

18、若(ruo)待糾偏數據(ju)調制類型為8psk或8apsk，則需對(dui)其符號對(dui)應外環(huan)星(xing)座(zuo)(zuo)點(dian)(dian)(dian)(dian)進(jin)行2次方降階操作(zuo)；若(ruo)轉換后的星(xing)座(zuo)(zuo)點(dian)(dian)(dian)(dian)與(yu)qpsk標準星(xing)座(zuo)(zuo)點(dian)(dian)(dian)(dian)相(xiang)(xiang)比(bi)存(cun)在π/4的相(xiang)(xiang)位差，則需對(dui)該星(xing)座(zuo)(zuo)點(dian)(dian)(dian)(dian)進(jin)行相(xiang)(xiang)位補償，得到與(yu)qpsk相(xiang)(xiang)位一(yi)致的四相(xiang)(xiang)星(xing)座(zuo)(zuo)點(dian)(dian)(dian)(dian)；

19、如待糾偏數據調制(zhi)類型為(wei)16apsk，則降階算法(fa)中的m值表示為(wei)：

20、m＝y/4

21、其中，y為(wei)數(shu)據(ju)星(xing)(xing)座圖對(dui)應的(de)外環數(shu)據(ju)點(dian)數(shu)；之后需(xu)劃定內外圈(quan)功率(lv)門限，對(dui)大(da)于門限的(de)外圈(quan)信號(hao)做m次降階(jie)，對(dui)小(xiao)于功率(lv)門限的(de)內圈(quan)信號(hao)置0；若轉(zhuan)換后的(de)星(xing)(xing)座點(dian)與qpsk標準(zhun)星(xing)(xing)座點(dian)相(xiang)比存在π/4的(de)相(xiang)位(wei)差，則需(xu)對(dui)該星(xing)(xing)座點(dian)進行相(xiang)位(wei)補償，得到與qpsk調制一致的(de)四相(xiang)星(xing)(xing)座點(dian)。

22、本發明進一步的技術方案：所述(shu)s2包括：

23、s21：對信號做能(neng)量歸(gui)一化，獲得信號s1(i)；

24、s22：將(jiang)s1(i)由直(zhi)角極(ji)坐標系轉化到極(ji)坐標系，得到信號的相位(wei)(wei)信息(xi)θi＝2πfdits+φ0+φz，其中，fd為(wei)調制頻率信息(xi)，φ0為(wei)初始相位(wei)(wei)，φz為(wei)相位(wei)(wei)域高斯白噪(zao)聲分量；

25、對θi進行4倍(bei)放(fang)大去除信號調制(zhi)信息得(de)到：

26、4θi＝8πfdits+4φ0+4φz

27、將其映射到直(zhi)角坐標系有：

28、s2(i)＝exp(j(8πfdits+4φ0+4φz)。

29、本(ben)發明(ming)進一步的技術(shu)方案：所(suo)述s3包(bao)括：

30、以fft計算點(dian)數n為(wei)長度，對信號s2(i)進(jin)(jin)行分段，并對每段信號進(jin)(jin)行fft變換，用(yong)公式表示為(wei)：

31、

32、其(qi)中，hn為分段信(xin)號，k＝0,1,…,n-1；

33、利用(yong)公式|fk|＝fkfk*/n求(qiu)得頻(pin)譜幅值，對所(suo)有(you)數(shu)據段(duan)的頻(pin)譜幅值進行(xing)疊(die)加求(qiu)和，在求(qiu)和頻(pin)譜結果中搜索出最大峰值索引kmax1。

34、本發明(ming)進一步的技術方案：所述s4包括：

35、由最大似然估計理論得到頻偏(pian)的估計表達(da)式如(ru)下：

36、

37、式中(zhong)，ts為(wei)(wei)采樣間隔、kmax1為(wei)(wei)頻譜峰值索引、n為(wei)(wei)fft點數、m為(wei)(wei)qpsk調制指(zhi)數，頻偏粗(cu)估計(ji)值

38、本發明進(jin)一(yi)步的技術方案：所述s6包括：

39、s61：對(dui)重采(cai)樣后的(de)信(xin)號x(n)′以n為長度進行(xing)分段，按照s3中的(de)步(bu)驟獲得信(xin)號疊(die)加頻(pin)譜(pu)的(de)峰值索引kmax2；

40、s62：細化峰值(zhi)兩(liang)側索引(yin)之間(jian)的(de)(de)頻(pin)譜；具體的(de)(de)，先計算每(mei)段信(xin)號對(dui)應(ying)的(de)(de)頻(pin)譜幅值(zhi)點(dian)|xn(kmax2±0.5)|對(dui)該(gai)段信(xin)號[kmax2-1,kmax2+1]之間(jian)的(de)(de)頻(pin)譜進(jin)(jin)行細化；通過對(dui)每(mei)段信(xin)號后綴填補長度為(wei)n的(de)(de)全(quan)0數據，由公式(shi)計算得到(dao)每(mei)段信(xin)號的(de)(de)|xn(kmax2±0.5)|，最后進(jin)(jin)行求和平均得到(dao)

41、s63：計(ji)算(suan)補(bu)償(chang)因子δ，并對峰(feng)值索引kmax2進行補(bu)償(chang)，通過小數(shu)部分的(de)補(bu)償(chang)使其(qi)更逼近實(shi)際(ji)頻(pin)(pin)譜峰(feng)值位置；具體的(de)，結合fft計(ji)算(suan)表達(da)(da)式以及(ji)頻(pin)(pin)偏估計(ji)表達(da)(da)式，獲得(de)fft表達(da)(da)式的(de)另一(yi)種形式：

42、

43、其中(zhong)，kmax為(wei)(wei)峰值(zhi)索(suo)引(yin)，δf為(wei)(wei)頻(pin)偏(pian)值(zhi)，ts為(wei)(wei)采樣間隔，n為(wei)(wei)fft計算點數，θ為(wei)(wei)數據初始相位；則補償因子δ的值(zhi)在-0.5到0.5之間，得(de)出下列關系：

44、

45、依據(ju)三(san)角(jiao)函數性質(zhi)對(dui)上(shang)式進行(xing)簡化(hua)有：

46、

47、計算補(bu)償(chang)因子δ的表達式為：

48、

49、最后利用δ補償后的峰值索引表(biao)示為(wei)：kmax＝kmax2-δ。

50、本發明進一步的(de)技術方(fang)案：所述s7包括：

51、根據最(zui)大似然估計理論，頻偏(pian)細估計值(zhi)與補償后的峰值(zhi)索引關系式如(ru)下(xia)所示：

52、

53、其中，fs為(wei)(wei)采(cai)樣間隔、kmax為(wei)(wei)補(bu)償后的(de)頻譜峰值索引、n為(wei)(wei)fft變(bian)換(huan)點數、m為(wei)(wei)qpsk調制指(zhi)數、s為(wei)(wei)重采(cai)樣抽取因(yin)子，頻偏細估計值為(wei)(wei)

54、一(yi)種低(di)信噪比下非數據輔助的(de)高精(jing)度頻偏估計補償方法的(de)應用(yong)，其特征在(zai)于應用(yong)于寬(kuan)帶(dai)衛星通信物理層解調系(xi)統載(zai)波恢復。

55、本發明的有(you)益(yi)效果在(zai)于：

56、本發明涉及(ji)一(yi)種低信(xin)(xin)(xin)噪(zao)比下非數(shu)據(ju)(ju)(ju)輔助的(de)(de)(de)(de)(de)(de)(de)高精度頻(pin)偏(pian)估計補(bu)(bu)償(chang)方法，其無需利(li)用(yong)已知的(de)(de)(de)(de)(de)(de)(de)先驗數(shu)據(ju)(ju)(ju)，可(ke)應用(yong)于(yu)寬帶(dai)衛(wei)(wei)星物理層解(jie)調(diao)系(xi)統中的(de)(de)(de)(de)(de)(de)(de)載(zai)波恢(hui)復，包括：對(dui)(dui)經過符號(hao)(hao)定時(shi)恢(hui)復以及(ji)幀同步(bu)的(de)(de)(de)(de)(de)(de)(de)信(xin)(xin)(xin)號(hao)(hao)進(jin)行(xing)(xing)m次乘方降階獲得(de)四(si)相(xiang)星座(zuo)點(dian)(dian)信(xin)(xin)(xin)號(hao)(hao)；對(dui)(dui)信(xin)(xin)(xin)號(hao)(hao)進(jin)行(xing)(xing)能(neng)量歸一(yi)化操(cao)(cao)作，將(jiang)其由(you)直角(jiao)坐(zuo)標系(xi)轉(zhuan)為(wei)極坐(zuo)標形式(shi)，并將(jiang)其相(xiang)位數(shu)據(ju)(ju)(ju)擴大4倍，重(zhong)新映射到直角(jiao)坐(zuo)標系(xi)，得(de)到去除調(diao)制信(xin)(xin)(xin)息的(de)(de)(de)(de)(de)(de)(de)信(xin)(xin)(xin)號(hao)(hao)；根(gen)據(ju)(ju)(ju)fft點(dian)(dian)數(shu)對(dui)(dui)信(xin)(xin)(xin)號(hao)(hao)進(jin)行(xing)(xing)分段，對(dui)(dui)每段信(xin)(xin)(xin)號(hao)(hao)完(wan)成fft變(bian)換(huan)，并對(dui)(dui)其頻(pin)譜進(jin)行(xing)(xing)累(lei)加(jia)，依據(ju)(ju)(ju)信(xin)(xin)(xin)號(hao)(hao)的(de)(de)(de)(de)(de)(de)(de)累(lei)加(jia)頻(pin)譜進(jin)行(xing)(xing)峰(feng)(feng)值(zhi)搜(sou)(sou)索，獲得(de)周期圖峰(feng)(feng)值(zhi)索引(yin)(yin)；根(gen)據(ju)(ju)(ju)頻(pin)譜峰(feng)(feng)值(zhi)索引(yin)(yin)以及(ji)fft點(dian)(dian)數(shu)、調(diao)制階數(shu)等信(xin)(xin)(xin)息獲得(de)頻(pin)偏(pian)粗估計值(zhi)，完(wan)成對(dui)(dui)信(xin)(xin)(xin)號(hao)(hao)的(de)(de)(de)(de)(de)(de)(de)第一(yi)級粗補(bu)(bu)償(chang)；對(dui)(dui)粗偏(pian)補(bu)(bu)償(chang)后的(de)(de)(de)(de)(de)(de)(de)信(xin)(xin)(xin)號(hao)(hao)進(jin)行(xing)(xing)相(xiang)同的(de)(de)(de)(de)(de)(de)(de)降階去調(diao)制操(cao)(cao)作，并進(jin)行(xing)(xing)抽取(qu)因子為(wei)s的(de)(de)(de)(de)(de)(de)(de)重(zhong)采樣(yang)，得(de)到采樣(yang)率為(wei)fs/s的(de)(de)(de)(de)(de)(de)(de)信(xin)(xin)(xin)號(hao)(hao)；將(jiang)重(zhong)采樣(yang)的(de)(de)(de)(de)(de)(de)(de)信(xin)(xin)(xin)號(hao)(hao)進(jin)行(xing)(xing)分段fft計算并對(dui)(dui)其頻(pin)譜進(jin)行(xing)(xing)累(lei)加(jia)，對(dui)(dui)累(lei)加(jia)頻(pin)譜進(jin)行(xing)(xing)搜(sou)(sou)峰(feng)(feng)操(cao)(cao)作獲得(de)峰(feng)(feng)值(zhi)索引(yin)(yin)；進(jin)行(xing)(xing)插值(zhi)計算獲得(de)索引(yin)(yin)補(bu)(bu)償(chang)因子，逼(bi)近峰(feng)(feng)值(zhi)實際位置；根(gen)據(ju)(ju)(ju)補(bu)(bu)償(chang)后的(de)(de)(de)(de)(de)(de)(de)峰(feng)(feng)值(zhi)索引(yin)(yin)、抽取(qu)因子s以及(ji)fft點(dian)(dian)數(shu)等參數(shu)進(jin)行(xing)(xing)細(xi)頻(pin)偏(pian)估計并完(wan)成對(dui)(dui)信(xin)(xin)(xin)號(hao)(hao)的(de)(de)(de)(de)(de)(de)(de)補(bu)(bu)償(chang)。本發明可(ke)在低信(xin)(xin)(xin)噪(zao)比環(huan)境下進(jin)行(xing)(xing)高精度頻(pin)偏(pian)估計，有效提升解(jie)調(diao)性能(neng)，適用(yong)于(yu)寬帶(dai)衛(wei)(wei)星通信(xin)(xin)(xin)系(xi)統以及(ji)衛(wei)(wei)星高速(su)數(shu)傳系(xi)統。