CN103634255B - 快速估算msk信号载波的方法 - Google Patents

快速估算msk信号载波的方法 Download PDF

Info

Publication number
CN103634255B
CN103634255B CN201210296783.3A CN201210296783A CN103634255B CN 103634255 B CN103634255 B CN 103634255B CN 201210296783 A CN201210296783 A CN 201210296783A CN 103634255 B CN103634255 B CN 103634255B
Authority
CN
China
Prior art keywords
signal
msk
frequency
carrier
nco
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.)
Expired - Fee Related
Application number
CN201210296783.3A
Other languages
English (en)
Other versions
CN103634255A (zh
Inventor
张玲
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
CETC 10 Research Institute
Original Assignee
CETC 10 Research Institute
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by CETC 10 Research Institute filed Critical CETC 10 Research Institute
Priority to CN201210296783.3A priority Critical patent/CN103634255B/zh
Publication of CN103634255A publication Critical patent/CN103634255A/zh
Application granted granted Critical
Publication of CN103634255B publication Critical patent/CN103634255B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Digital Transmission Methods That Use Modulated Carrier Waves (AREA)

Abstract

本发明提出一种快速估算MSK信号载波的方法，旨在提供一种能在测控系统中更准确快速的对MSK信号载波频率进行估算的方法。本发明通过下述技术方案予以实现：首先对来自测控系统的MSK测控信号，经模/数转换器A/D采样后转换为数字信号输入FPGA。FPGA将该数字信号与本地数字控制振荡器NCO的输出信号相乘，再经过数字低通滤波器滤除高频分量，对残留信号进行平方处理和快速傅立叶变换FFT，查找出幅度值最大的两个频率点，再经过最大值有效验证算法，对提取出的两个频率点进行是否有效的判断，有效则保留，无效则剔除，最后根据保留下的有效频率点计算出载波频率的估算值送入MSK接收机模块完成MSK信号的解调。

Description

快速估算MSK信号载波的方法
技术领域
本发明涉及一种测控系统中MSK信号快速载波估算的方法。
背景技术
最小频移键控（MSK)调制是一种二进制连续相位的恒包络相位调制,MSK信号是一种相位连续、包络恒定,占用带宽最小的二进制正交2FSK信号,因其频带利用率高,带外辐射小,被广泛地利用在数字通信系统中。
在目前通常采用MSK信号体制工作的一些测控系统中，MSK信号的解调一般采用差分数字解调算法，在接收机处理过程中当每个符号的采样点数大于5时，该算法解调性能已经接近MSK的最佳解调性能，并且此算法对频偏和相差有一定的容忍度。所以当载波多普勒变化较小时，采用差分数字解调算法可不必对载波频率进行估算。
然而，在无线通信系统中，由于多普勒频移和振荡器的不稳定等因素，使得接收信号的载波和本地载波有时存在较大偏差，此偏差的存在会严重影响MSK解调器的性能。在应用过程中发现，由于载波多普勒太大，变化率太快，会直接导致接收机接收解调性能下降，系统误码率增加。在应用过程中还发现，相对较大的载波多普勒频率，当调制的信息速率相对低速时，此时如果不对载波频率进行估算，接收机采用同样的接收解调算法，接收解调性能较最佳解调性能有所下降。所以为了使MSK接收机实现最佳MSK接收解调性能，对MSK信号载波进行估算是必要的。
发明内容
本发明的任务是针对现有技术测控系统中存在的上述不足之处，提供一种能够正确快速估算MSK信号载波的方法，尤其是能够准确快速的对测控系统中MSK信号进行载波估算的方法，以解决在测控系统中由于载波多普勒太大，变化率太快，导致接收机接收解调性能下降，系统误码率增加的问题。
本发明解决其技术问题通过以下措施予以实现：首先对来自测控系统的MSK测控信号，经模/数转换器A/D采样后转换为数字信号输入FPGA，FPGA将该数字信号与本地数字控制振荡器NCO的输出信号相乘，再经过数字低通滤波器滤除高频分量后，对残留信号进行平方处理，然后对平方后的数字信号S'(t)2进行快速傅立叶FFT变换，在FFT输出信号中查找出幅度值最大的两个频率点fO1、fO2，加入最大值有效验证算法，利用MSK信号特点，根据公式推导得出MSK调制信息频率fs=|fO1-fO2|，再对提取出的两个频率点进行是否有效的判断，有效则保留，无效则剔除，最后根据保留下的有效频率点计算出载波频率的估算值，再把估算的载波频率值送入MSK接收机模块完成MSK信号的解调。
本发明相比于现有技术具有如下有益效果：
本发明利用MSK的信号特点对载波信号进行准确快速的估算。相比未进行载波估算的接收机，在相同的信噪比下，加入了载波估算模块的MSK接收机其性能能够更加接近MSK最佳接收解调性能。在应用过程中，在FPGA内进行信号处理时可以改变快速傅里叶变换(FFT)计算的采样率及分辨率以适应不同的最大多普勒频率及动态变化率。
本发明采用的载波估算方法相比一般扩频信号载波频率估算方法，针对MSK信号的频谱特点提出了幅度值最大频率点是否有效的判决条件，提高了载波频率估算的可靠性，避免出现错误的载波频率估算值被预置到接收模块中，导致接收机无法正常工作的情况。解决了测控系统载波多普勒太大，变化率太快，易导致接收机接收解调性能下降，系统误码率增加的问题。
附图说明
下面结合附图和实施例对本发明进一步说明。
图1是本发明MSK信号载波估算模块的电路原理框图。
具体实施方式
参阅图1。在以下描述的实施例中，来自测控系统的MSK测控信号通过模/数转换器A/D采样电路转换为数字信号后输入FPGA，输入FPGA的数字信号可表示为其中dk为调制信息，f0为接收信号载波频率，fs为调制信息频率。
在FPGA中，输入的数字信号S(t)与本地数字控制振荡器NCO输出的本地载波信号Fnco相乘后，经过数字低通滤波器滤除高频分量，滤波后的数字信号可以表示为其中Δf＝f0-fnco。接下来在FPGA内部对滤波后的数字信号S'(t)进行平方处理，平方后的数字信号表示为其中Δf＝f0-fnco，Δf即为所求的多普勒频率值。从平方后的数字信号S'(t)2的公式分析可以得出，平方去除了因调制信息引起的频率变化，而平方后的数字信号S'(t)2只包含了2种单一正弦信号分量。然后对平方后的数字信号S'(t)2信号进行FFT变换，注意做FFT变换时根据系统中最大多普勒频率及动态变化率选取合适的采样率及分辨率。FFT变换后从FFT的输出信号中提取出幅度值最大的两个频率点，分别为
利用fO1、fO2计算可得出接收多普勒频率值Δf＝(fO1+fO2)/4，从而计算出接收载波频率f0＝Δf+fnco。
为了验证提取出的两个频率点fO1、fO2的正确性，加入最大值有效验证算法。FFT变换后提取出的幅度值最大的两个频率点fO1和fO2，根据公式推导可得出|fO1-fO2|＝fs，其中fs为MSK调制信息频率，为已知条件。总结验证频率点正确与否的判决条件如下：对两个最大频率点fO1，fO2做减法运算，求模后，其差值应该等于调制信息频率fs。如果符合上诉判决条件，则保留频率点值，不符合，则剔除。
最后保留下来的正确的频率点值fO1，fO2，根据计算公式两式相加可得出Δf＝(fO1+fO2)/4，Δf即为多普勒频率值，然后利用公式f0＝Δf+fnco可计算出接收信号的载波频率，最后把估算的载波频率值f0送入MSK接收机模块完成MSK信号的解调。

Claims (5)

1.一种快速估算MSK信号载波的方法，具有如下技术特征：首先对来自测控系统的最小频移键控MSK测控信号，经模/数转换器A/D采样后转换为数字信号输入FPGA，FPGA将该数字信号与本地数字控制振荡器NCO的输出信号相乘，再经过数字低通滤波器滤除高频分量后，对残留信号进行平方处理，然后对平方后的数字信号S'(t)2进行快速傅立叶FFT变换，在FFT输出信号中查找出幅度值最大的两个频率点fO1、fO2，加入最大值有效验证算法，根据公式推导得出MSK调制信息频率fs＝|fO1-fO2|，再对提取出的两个频率点进行是否有效的判断，有效则保留，无效则剔除；验证频率点正确与否的判决条件如下：对两个最大频率点fO1，fO2做减法运算，求模后，其差值等于MSK调制信息频率fs，如果符合上述判决条件，则保留频率点值，不符合，则剔除，其中，fsMSK调制信息频率为已知条件；最后保留下来的正确频率点值fO1和fO2，根据计算公式两式相加得出接收机的接收信号多普勒频率值Δf＝(fO1+fO2)/4，然后利用公式f0＝Δf+fnco计算出接收信号载波频率f0，根据保留下的有效频率点计算出接收信号载波频率的估算值，再把估算的接收信号载波频率值送入MSK接收机模块完成MSK信号的解调，其中，fnco为本地NCO的本地载波信号频率。
2.按权利要求1所述的快速估算MSK信号载波的方法，其特征在于，MSK测控信号经A/D采样后转换为数字信号的表达式为
式中dk为调制信息，f0为接收信号载波频率。
3.按权利要求2所述的快速估算MSK信号载波的方法，其特征在于，滤波后的数字信号
其中，Δf为接收信号多普勒频率值。
4.按权利要求3所述的快速估算MSK信号载波的方法，其特征在于，本地NCO产生的本地载波信号表达式为Fnco＝cos(2πfncot)，由A/D模/数转换器转换输入FPGA的数字信号S(t)与本地NCO输出信号Fnco相乘后，再经过数字低通滤波器滤除高频分量，余下的残留信号可表达为
5.按权利要求4所述的快速估算MSK信号载波的方法，其特征在于，平方后的数字信号为 即平方后的数字信号包含了2种单一正弦信号分量，且平方后的数字信号S'(t)2与调制信息dk无关。
CN201210296783.3A 2012-08-20 2012-08-20 快速估算msk信号载波的方法 Expired - Fee Related CN103634255B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210296783.3A CN103634255B (zh) 2012-08-20 2012-08-20 快速估算msk信号载波的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210296783.3A CN103634255B (zh) 2012-08-20 2012-08-20 快速估算msk信号载波的方法

Publications (2)

Publication Number Publication Date
CN103634255A CN103634255A (zh) 2014-03-12
CN103634255B true CN103634255B (zh) 2017-04-05

Family

ID=50214889

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210296783.3A Expired - Fee Related CN103634255B (zh) 2012-08-20 2012-08-20 快速估算msk信号载波的方法

Country Status (1)

Country Link
CN (1) CN103634255B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103973622B (zh) * 2014-05-26 2017-03-01 许昌学院 一种快捷的msk通信检测解调方法
CN105629060B (zh) * 2015-12-24 2018-05-29 电子科技大学 基于最优基带滤波的电网频率测量方法和装置
CN105704081B (zh) * 2016-01-12 2018-10-23 中国电子科技集团公司第十研究所 多个偏移载波调幅信号的优选和解调方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1175171A (zh) * 1997-07-24 1998-03-04 北京信威通信技术有限公司 扩频通信系统中载波恢复和补偿的方法及其装置
US6631170B1 (en) * 1999-12-03 2003-10-07 Nokia Corporation Radio frequency receiver

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1175171A (zh) * 1997-07-24 1998-03-04 北京信威通信技术有限公司 扩频通信系统中载波恢复和补偿的方法及其装置
US6631170B1 (en) * 1999-12-03 2003-10-07 Nokia Corporation Radio frequency receiver

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
《MSK数字调制解调及其实现技术研究》;唐良伟;《中国优秀硕士学位论文全文数据库信息科技辑》;20071215(第06期);正文第36-58页 *
《一种卫星信号载波频率精确估计算法》;黄富彪等;《现代电子技术》;20120201;第35卷(第3期);第544-547页 *
《正弦信号的直接FFT参数估计与相位差分法对比研究》;李辉等;《电子与信息学报》;20100331;第32卷(第3期);第129-131，134页 *

Also Published As

Publication number Publication date
CN103634255A (zh) 2014-03-12

Similar Documents

Publication Publication Date Title
CN101938445B (zh) Iq不平衡估计及补偿方法
CN108055224B (zh) 16qam载波同步环路同步锁定检测方法
CN106603454B (zh) 一种gfsk数字双通道解调方法
CN101479948B (zh) 相邻信道干扰检测装置及方法
CN103634255B (zh) 快速估算msk信号载波的方法
CN108512791A (zh) 基于定时频偏补偿的星载ais解调方法
CN103547344A (zh) 一种周跳的检测方法、装置及接收机
CN109889461B (zh) 一种低复杂度并行的载波恢复系统及其方法
CN107483078A (zh) 一种船舶vdes系统asm体制接收频偏估计实现方法
CN106508104B (zh) 一种扩展遥测相干接收机频偏估计范围的方法
CN105704081A (zh) 多个偏移载波调幅信号的优选和解调方法
CN111343116B (zh) 密集信号环境下基于星座轨迹图的调制识别方法
CN105516041A (zh) 一种低信噪比下自适应数字解调系统
CN102647381A (zh) 一种mpsk相干光通信系统中频率偏移估计方法及装置
CN108650203B (zh) 一种基于侦察接收机的调制模式识别方法
CN116016072B (zh) 零中频结构低复杂度msk正交解调装置及其解调方法
CN105391539B (zh) Qpsk载波同步锁定检测器
CN1269339C (zh) 数字qam接收器的载波恢复装置
CN103024888B (zh) 接收机数字自动增益控制方法与装置
CN113765838B (zh) 一种dpsk信号的解调方法、装置、设备及存储介质
CN1193562C (zh) 用于角调制信号的接收装置
CN105811590A (zh) 无线充电设备中实现调幅信号解调与解码的方法
CN104242928B (zh) 锁频环中一种基于叉积算法的频率锁定检测器处理方法
CN104580047B (zh) 一种fsk解调的软件实现方法
CN115580515B (zh) 基于多路判决的msk解调方法、装置、设备和存储介质

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20170405

Termination date: 20190820