CN102299737B - 一种多路快跳频信号处理方法 - Google Patents
一种多路快跳频信号处理方法 Download PDFInfo
- Publication number
- CN102299737B CN102299737B CN201110244405.6A CN201110244405A CN102299737B CN 102299737 B CN102299737 B CN 102299737B CN 201110244405 A CN201110244405 A CN 201110244405A CN 102299737 B CN102299737 B CN 102299737B
- Authority
- CN
- China
- Prior art keywords
- sequence
- road
- value
- demodulation
- branch road
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Landscapes
- Digital Transmission Methods That Use Modulated Carrier Waves (AREA)
Abstract
本发明公开了一种多路快跳频信号处理方法,包括对以4N/Tb的采样率获得的采样序列进行预处理,形成两个4N点序列并分别进行FFT变换的步骤;对预处理后序列进行解调和参数估计的步骤。采用本发明在保证算法精度的同时大大的减少了运算规模,降低了星上处理复杂度。
Description
技术领域
本发明公开了一种多路快跳频信号处理方法。
背景技术
新一代通信卫星通过采用跳频通信技术,提供抗干扰移动通信业务,具有抗干扰能力强,支持的用户多的特点。
该***是一个支持多用户的卫星跳频通信***,由于地面用户数很多,为了节省星上资源,星上解跳时多个用户公用一个解跳频本振,所以卫星上行链路的多用户采用先频分复用FDMA再跳频的体制。在FSK调制方式中,频率的一次变换传输一个码元符号,FDMA信道间隔4倍的符号速率,星上使用模拟解跳,将解跳以后的固定中频多路信号使用一个AD统一采样,然后进行数字处理,需要进行数字分路,解调,同步误差估计,频偏估计信道功率估计等处理,按照常规的方法完成这些处理需要非常大的星上处理资源。
发明内容
本发明的技术解决问题是:克服现有技术的不足,提供了一种多路快跳频信号处理方法。采用本发明所述方法在保证算法精度的同时大大的减少了运算规模,降低了星上处理复杂度。
本发明的技术解决方案是:
利用FSK频率调制的特点,将FSK解调、参数估计与FDMA信号数字下变频分路相结合,利用数字分路运算直接产生解调和参数估计算法中需要用到的相关卷积的结果。具体方法步骤如下:
(1)对输入的多路FDMA信号以4N/Tb的采样率进行AD进行采样,获得采样共4N点的采样序列。
(2)采样序列处理:
将采样序列分为前2N点采样序列和后2N点采样序列。
分别在前2N点采样序列后和后2N点采样序列前补充2N点个0值,构成两个4N点的序列。
(3)对两个4N点的序列进行FFT变换,分别获得对应于前半符号2N点采样后序列的累加值sfi和对应于后半符号2N点采样后序列的累加值sbi,0≤i≤4N-1。
按频点将sfi和sbi相加,获得序列si;若对信息序列进行解调,则转入步骤(4),若对信息序列进行参数估计处理,则转入步骤(5)。
(4)信息序列解调:
根据步骤(3)中获得的序列si,对第i路信息序列进行判决获得对输入信息序列的解调结果。
(5)参数估计处理:
根据步骤(4)中获得的解调结果,选取每路中的f0或者f1支路的前半符号和后半符号累加值,进行参数估计,分别获得第i路信息序列的功率值、频率偏移和跳频同步误差。
本发明与现有技术相比具有如下优点:
(1)采用本发明所述方法可以使得处理性能达到理论值,常规的数字分路算法用来处理多路快跳频FDMA信号性能比较差,原因是跳频同步误差会被分路的低通滤波器滤掉,无法提取出同步误差信号,从而导致整个***同步不准,传输性能下降。
(2)本发明方法中数字处理算法运算量小,原因是将N路信号的解调以及参数估计中运算量大的运算都集中到数字分路中的FFT运算中,从而大大的节省了数字处理的硬件规模,需要的数字电路的逻辑门数仅有常规数字分路解调的三分之一左右。
附图说明
图1为数字分路原理图;
图2为FSK解调原理图;
图3为跳频同步误差提取原理图;
图4为功率估计原理图;
图5为鉴频原理图;
图6为本发明流程图。
具体实施方式
下面就结合附图对本发明做进一步介绍。
对于本发明相关的技术进行初步的介绍。
数字分路技术
一般数字分路算法采用FFT分路算法,FFT分路算法可以等效为图1的处理过程。图1所示的方法等效为对输入的一路FDMA信号只使用一个AD统一采样,然后分别利用多个不同的载波进行下变频,得到一组分离的零中频基带信号,再对获得的零中频基带信号进行低通滤波后便可得到对FDMA信道的分路输出。
解调技术
在调频***中每跳传输一个码元,由于跳频的影响,在码元转换时刻会产生了一个随机相位,因此必须采用非相干解调算法见图2。提取位于f0和f1两个频点的码元序列,并分别进行累加,然后再利用累加后的序列进行判断。
跳频同步误差估计方法
利用一个跳频驻留时间内,前半跳时间中的信号能量的累加平均值与后半跳时间中的信号能量的累加平均值比较,给出跳频超前或滞后信息。如果前半码元能量大于后半码元能量,则地面跳频相位超前星上跳频相位;如果前半码元能量小于后半码元能量,则地面跳频相位滞后星上跳频相位,见图3。
功率估计方法
分别对f0和f1支路进行匹配相关然后取模值,取出模值大的一路作为功率值,见图4。
频偏估计方法
在1跳传输1比特的快跳***中的载波频率在不停的跳变,由于技术实现上的困难和复杂性,无法保证载波跳变过程中相位的连续性。因此将鉴频放在每一跳中,即每次利用一跳中的采样点计算得到一次频偏。根据解调判决的结果,选取f0或者f1支路的变频结果,然后按照图5的方法求出频偏。
本发明技术方案
本发明中对多路快跳频信号分路解调参数估计联合处理方法的思想是:将图1中的低通滤波器采用矩形窗函数,同时将分路间隔设置为符号速率,正好可以将FFT算法等效为FDMA的FSK信号的非相干解调。从同步误差估计和频偏估计算法可以看出这两个算法都需要前半码元的累加值和后半码元的累加值,所以进一步取矩形窗函数的宽度为Tb/2,其中Tb为一个调制符号的持续时间,这样正好产生前半码元的累加值和后半码元的累加值,可以供跳时差估计、频偏估计、解调使用。
FFT变换的数学表达式为:
其中,
WN=e-j2π/N (2)
则,
由此可得:
对于任何一个给定的k值,则可以认为,X(k)是序列x(n)下变频的结果,下变频的频率f为:
根据采样定理可知模拟的采样频率fs对应着采样后数字域的2π,这样可以得出:
对一般FDMA的信号,预先规定每一路信号带宽相同,等间隔的在频域排列,信道间隔B可以从表达式(7)中得出:
这样就可以使用FFT来实现等间隔的数字下变频,下变频以后再进行低通滤波,滤掉其它频点的信号,得到N路零中频基带信号。对于图1中数字分路里的低通滤波器组通过线性运算等效到FFT运算之前,可以用一个低通滤波器来实现,提高处理效率。对于这一经典的数字分路算法在“多抽样率数字信号处理”,美国人RE.克劳切的这本书里面有详细的介绍,在这里就不详细介绍。
这里我们要处理的FDMA信号是N路FSK信号,排列方式是:信道间隔为4个基本带宽fB,其中fB=1/Tb,在四个基本带宽中的fB,2fB频点,对应FSK调制的f0,f1,FSK信号调制谱的范围从0到3fB,从3fB到4fB这一段为相邻信道的保护带,N路信号的总带为4N/Tb。
利用数字分路的原理,将信道间隔假设为1/Tb,信号总带宽还是4N/Tb,进行4N点的FFT运算,此时下变频频率间隔为1/Tb,下变频载波为0≤k≤4N-1,FFT输出的4k+1,4k+2处的值对应着第k路FSK信号f0,f1支路的相关值。
由于图3所示的跳频同步误差估计以及图5所示的跳频信道频偏估计都需要用到前半符号采样点累加值和后半符号采样点累加值,这时考虑将数字分路中的低通滤波器设为:
滤波以后的信号y(n)=x(n)*h(n),
将(9)代入(10)可得:
采用了公式(9)所示的滤波器对信号滤波,可以等效为对信号半个周期采样点进行累加。根据数字分路的原理,将滤波器放在FFT之前实现,可以提高运算效率,具体做法是采样序列x(n)与h(n)对应位置相乘,然后进行旋转变换,最后进行FFT运算,FFT结果就是各路信号的前半符号和后半符号的采样点累加值。
基于上述原理,本发明所述方法的步骤如下:
(1)对输入的多路FDMA信号使用一个AD进行采样,采样率为4N/Tb,可获得4N点的采样后序列。
(2)对获得的4N点的采样序列进行处理:
(21)4N点的采样序列分为前2N点的采用序列和后2N点的采样序列。
(22)分别在前2N点采样序列后和后2N点采样序列前补充2N点个0值,构成两个4N点的序列;
(3)基于上述的原理,对两个4N点序列分别进行快速傅里叶FFT变换,对应于包含前2N点的序列,获得变换后的序列sfi;对应于包含后2N点的序列,获得变换后的序列sbi,0≤i≤4N-1。
将sfi和sbi中的点按照FSK调制中的频点f0和f1对应相加,获得序列si。si即为叠加后分别位于f0和f1两个频点的码元序列的叠加。
若需要对信息序列进行解调结果进行判决,则转入步骤(4);若需要进行参数估计处理,则转入步骤(5)。
(4)对于第k路信号的解调,需要用第k路信号f0,f1支路的能量值进行比较判决,第k路的f0支路对应着si的4*(k-1)+1,f1支路对应着si的4*(k-1)+2,取出si序列4*(k-1)+1点的模值和4*(k-2)+2点的模值比较,若|s4(k-1)+1|≥|s4(k-1)+2|,则第k路解调结果为0,否则解调结果为1,完成第k路的解调,k取值范围为1≤k≤N;
(5)参数估计处理
功率值
根据第k路f0支路、f1支路模值比较结果,取模值大的值为第k路功率值;
频率偏移
根据第k路f0支路f1支路模值比较结果,取出对应位置的sfi、sbi分别求相位,然后求出前半符号和后半符号的相位差,再除以半个符号的时间Tb/2,就得到频偏。
跳频同步误差
根据第k路f0支路f1支路模值比较结果,取出对应位置的sfi、sbi分别取模值,前半符号和后半符号的模值差就是跳频同步误差。
本发明未详细说明部分属本领域技术人员公知常识。
Claims (1)
1.一种多路快跳频信号处理方法,用于对经过解跳频处理之后获得的包含信息序列的多路FDMA信号进行解调和参数估计处理,其特征在于包括以下步骤:
(1)对输入的多路FDMA信号以4N/Tb的采样率进行AD采样,获得采样共4N点的采样序列;
(2)采样序列处理:
将采样序列分为前2N点采样序列和后2N点采样序列;
分别在前2N点采样序列后和后2N点采样序列前补充2N点个0值,构成两个4N点的序列;
(3)对两个4N点的序列进行FFT变换,分别获得对应于前半符号2N点采样后序列的累加值sfi和对应于后半符号2N点采样后序列的累加值sbi,0≤i≤4N-1;
按频点将sfi和sbi相加,获得序列si;若对信息序列进行解调,则转入步骤(4),若对信息序列进行参数估计处理,则转入步骤(5);
(4)信息序列解调:
根据步骤(3)中获得的序列si,对第i路信息序列进行判决获得对输入信息序列的解调结果;对于第k路信号的解调,需要用第k路信号f0,f1支路的能量值进行比较判决,第k路的f0支路对应着si的4*(k-1)+1,f1支路对应着si的4*(k-1)+2,取出si序列4*(k-1)+1点的模值和4*(k-2)+2点的模值比较,若|s4(k-1)+1|≥|s4(k-1)+2|,则第k路解调结果为0,否则解调结果为1,完成第k路的解调,k取值范围为1≤k≤N,f0=1/Tb,f1=2/Tb;
(5)参数估计处理:
根据步骤(4)中获得的解调结果,选取每路中的f0或者f1支路的前半符 号和后半符号累加值,进行参数估计,分别获得第i路信息序列的功率值、频率偏移和跳频同步误差,根据第k路f0支路、f1支路模值比较结果,取模值大的值为第k路功率值,根据第k路f0支路f1支路模值比较结果,取出对应位置的 sfi、 sbi分别求相位,然后求出前半符号和后半符号的相位差,再除以半个符号的时间Tb/2,就得到频偏,根据第k路f0支路f1支路模值比较结果,取出对应位置的sfi 、sbi 分别取模值,前半符号和后半符号的模值差就是跳频同步误差。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110244405.6A CN102299737B (zh) | 2011-08-23 | 2011-08-23 | 一种多路快跳频信号处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110244405.6A CN102299737B (zh) | 2011-08-23 | 2011-08-23 | 一种多路快跳频信号处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102299737A CN102299737A (zh) | 2011-12-28 |
CN102299737B true CN102299737B (zh) | 2014-04-02 |
Family
ID=45359930
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110244405.6A Active CN102299737B (zh) | 2011-08-23 | 2011-08-23 | 一种多路快跳频信号处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102299737B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102801663B (zh) * | 2012-06-28 | 2014-12-17 | 西安空间无线电技术研究所 | 一种用于深空探测器的多档码速率估计方法 |
CN105141338A (zh) * | 2015-06-11 | 2015-12-09 | 西安佳中电子技术有限公司 | 一种超短波信号数字化处理电路 |
CN105024670B (zh) * | 2015-07-17 | 2017-09-22 | 西安空间无线电技术研究所 | 跳频背景下的多调制方式多速率信号的快速数字滤波方法 |
CN110501728B (zh) * | 2018-05-16 | 2022-03-29 | 清华大学 | 定位基站跳时信号的鉴频方法及鉴频装置 |
CN110095791B (zh) * | 2019-05-28 | 2021-04-02 | 中国电子科技集团公司第五十四研究所 | 一种多载波fsk调制信号解调方法 |
CN110519192A (zh) * | 2019-08-23 | 2019-11-29 | 苏州浪潮智能科技有限公司 | 一种fsk信号解调方法、装置、嵌入式***及存储介质 |
CN110784864B (zh) * | 2019-11-01 | 2022-03-18 | 中国电子科技集团公司第三十研究所 | 一种星地链路设备射频指纹识别及安全接入认证方法 |
CN113595588B (zh) * | 2021-06-11 | 2022-05-17 | 杭州电子科技大学 | 一种基于时频谱熵的跳频信号感知方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101594167A (zh) * | 2008-05-27 | 2009-12-02 | 北京威讯紫晶科技有限公司 | 一种无线网络中联合信号检测和频率估计的方法 |
CN101667989A (zh) * | 2009-09-16 | 2010-03-10 | 中兴通讯股份有限公司 | 信号载波频率及相位估计方法及装置 |
CN102065043A (zh) * | 2010-12-09 | 2011-05-18 | 北京理工大学 | 一种高速通信***中的频域并行解调方法 |
-
2011
- 2011-08-23 CN CN201110244405.6A patent/CN102299737B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101594167A (zh) * | 2008-05-27 | 2009-12-02 | 北京威讯紫晶科技有限公司 | 一种无线网络中联合信号检测和频率估计的方法 |
CN101667989A (zh) * | 2009-09-16 | 2010-03-10 | 中兴通讯股份有限公司 | 信号载波频率及相位估计方法及装置 |
CN102065043A (zh) * | 2010-12-09 | 2011-05-18 | 北京理工大学 | 一种高速通信***中的频域并行解调方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102299737A (zh) | 2011-12-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102299737B (zh) | 一种多路快跳频信号处理方法 | |
CN102122997B (zh) | Lte主同步信号检测的方法、装置及终端 | |
CN103733516B (zh) | 接收机、发射机以及用于数字多子频带处理的方法 | |
CN102065043B (zh) | 一种高速通信***中的频域并行解调方法 | |
CN101567701B (zh) | 一种高效多路数字下变频器*** | |
US20140079164A1 (en) | Full-Range Pilot-Assisted Frequency Offset Estimation for OFDM Communication Systems | |
CN101827052B (zh) | 一种用于lte***时间同步和频率同步的方法及装置 | |
CN103268308B (zh) | 支持混合基dft的计算装置及方法 | |
WO2014099224A1 (en) | Narrow-band preamble for orthogonal frequency-division multiplexing system | |
CN106027179A (zh) | 一种基于综合互素分析的宽带频谱感知方法及其装置 | |
CN101667818B (zh) | 对实信号进行信道化处理的方法及装置 | |
CN102148679B (zh) | 低复杂度的宽带信号数字选频方法 | |
CN104730478A (zh) | 正交频分复用的磁共振相控阵接收方法 | |
CN102035764A (zh) | 非理想同步单载波频分多址***的信道估计装置及方法 | |
CN101395880A (zh) | 用于产生正交采样时钟信号的方法和设备 | |
CN101502018B (zh) | 同步信号的产生和检测 | |
RU122818U1 (ru) | Демодулятор фазоманипулированных сигналов | |
CN103916199A (zh) | 一种天线信号的时延和相位调整装置与方法 | |
CN106160969A (zh) | 一种lte下行同步数据发射配置与接收方法 | |
CN110113273B (zh) | 一种基于频谱对称的时频重叠信号盲分离方法 | |
CN105791197B (zh) | 一种偏移正交多载波***的脉冲成形滤波器 | |
JP3717038B2 (ja) | Ofdm信号のダイバーシティ受信装置 | |
CN106027454A (zh) | 基于频率补偿的qam信号的频偏估计方法 | |
CN102769444A (zh) | 用于解调多路abpsk信号的数字滤波器组 | |
CN102281096A (zh) | 一种资源复用数字信道化实现*** |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |