摘要：根據(jù)BPSK調(diào)制信號調(diào)制機理和平方倍頻法原理，在FPGA平臺上設(shè)計實現(xiàn)了BPSK調(diào)制信號載波頻率估計單元。利用ModelSim仿真環(huán)境對載頻估計功能進行仿真，驗證了平方倍頻法對BPSK信號進行載波信號估計的有效性。仿真表明基于FPGA的BPSK信號載頻估計單元，有較高的估計精度，且實現(xiàn)原理簡單，有一定的實際應(yīng)用價值。

0 引言

BPSK即二進制相移鍵控，是直擴信號中經(jīng)常使用的一種調(diào)制方式，利用載波的相位變化傳遞數(shù)字信息，信號的振幅、頻率保持恒定。BP SK調(diào)制方式具有較高傳輸效率、誤碼率低，不易受信道特性變化影響等特點，而且調(diào)制電路簡單易行，頻譜密度低，處理增益高，具有良好的低截獲概率可能，廣泛應(yīng)用于雷達、保密通信和導(dǎo)航定位等領(lǐng)域。對BPSK信號的載頻估計為后續(xù)的跟蹤捕獲等處理提供載頻參數(shù)，具有重大意義。隨著信號處理技術(shù)和檢測技術(shù)的飛速發(fā)展，涌現(xiàn)出了很多估計載波頻率的方法，如平方倍頻法、小波相關(guān)法等。

本文根據(jù)BPSK調(diào)制信號調(diào)制機理和平方倍頻法原理，在FPGA平臺下對基于平方倍頻法的BPSK調(diào)制信號載頻估計單元進行設(shè)計，并在Mode lsim6.5b環(huán)境下進行仿真驗證和結(jié)果分析。

1 平方倍頻法頻率估計原理

BPSK調(diào)制信號用初始相位0和π分別表示二進制“1”和“0”，BPSK信號的時域數(shù)學(xué)表達式可以表示為：
Sbpsk(t)=A·D(t)cos(2πft+φ) (1)

式中，A代表振幅，D(t)代表二進制信息，將D(t)與載波相乘，因D(t)只有兩種值，即“+1”和“-1”，分別代表“0”和“1”，使得BPSK調(diào)制信號只有兩種相位，則BPSK調(diào)制信號的生成原理圖如圖1所示。

根據(jù)BPSK的調(diào)制原理，利用二進制信息對載波信號進行相位調(diào)制，使載波信號相位突變，即BPSK信號同時含有載波信息和二進制信息。因此對BPSK調(diào)制信號的載頻估計應(yīng)該首先將二進制信息造成的相位突變消除，只留下載波或與載波有關(guān)的成分，再進行載頻估計。因為BPSK調(diào)制信號的二進制信息是±1構(gòu)成的序列，可以用平方處理消除二進制信息的影響，提取其中僅與載波有關(guān)的成分進行載波頻率估計。

BPSK調(diào)制信號的數(shù)學(xué)表達式如式(1)，對其平方，結(jié)果如式(2)。

由式(3)可知，平方后信號中包含BPSK調(diào)制信號載頻信號的平方項和直流分量。對平方結(jié)果進行傅里葉變換，求其頻譜，搜索頻譜峰值，并將譜峰位置輸出，則可以得到2fc的估計值，最后除以2即可得到BPSK調(diào)制信號的載波頻率。

2 基于FPGA的載頻估計單元設(shè)計

載波頻率估計單元首先實現(xiàn)對BPSK調(diào)制信號進行平方處理，然后將平方后的信號進行快速傅里葉變換(FFT)，對頻域進行二倍頻的頻域采樣點輸出，最后通過FFT變換的頻率分辨率與輸出采樣點的比例關(guān)系完成載頻的估計。

載頻估計單元實現(xiàn)邏輯框圖如圖2所示。

其中，F(xiàn)FT運算點數(shù)為cal_index，采樣頻率為fsample，而峰值所在位置為xk_index，而峰值所在位置是載波二倍頻所在位置，則載波頻率fcarry的計算公式為：

根據(jù)以上分析，對載頻估計模塊的設(shè)計主要分為四部分：乘法器、FFT單元、平方求和單元、判決單元。根據(jù)平方倍頻法原理以及邏輯框圖2，在FPGA平臺上設(shè)計的載波頻率估計單元FPGA結(jié)構(gòu)如圖3所示。

表1為圖3中的載頻估計單元的輸入輸出接口，以及各個接口實現(xiàn)的功能。

圖3中所設(shè)計的載頻估計模塊采用的算法是平方倍頻法，根據(jù)平方倍頻法原理，對輸入信號首先要進行平方處理。本單元使用XilinX公司提供的Multiplier IP核，版本為4.0。Multiplier IP核的兩個輸入信號為8位的有符號定點數(shù)，輸出信號是16位的有符號定點數(shù)。將平方后的信號進行傅里葉變換之前，需要對信號進行預(yù)處理。圖3中所示的載頻估計單元的關(guān)鍵模塊是傅里葉變換模塊。本單元使用的傅里葉模塊是由Xilinx公司提供的Fast Fourier Transform IP核，版本號為7.1。該IP核要求輸入數(shù)據(jù)為復(fù)數(shù)形式，因為通過乘法器計算后的數(shù)據(jù)是實數(shù)，因此對數(shù)據(jù)的預(yù)處理是加上一個為0的虛部，同時為了減少傅里葉變換的計算量，減少計算時間，這里將輸入數(shù)據(jù)進行截短，只留數(shù)據(jù)的前8位，然后傳送給Fast Fourier transform IP核進行計算。