CN110113278A - A kind of Modulation Mode Recognition method based on all-digital receiver - Google Patents

A kind of Modulation Mode Recognition method based on all-digital receiver Download PDF

Info

Publication number
CN110113278A
CN110113278A CN201910362765.2A CN201910362765A CN110113278A CN 110113278 A CN110113278 A CN 110113278A CN 201910362765 A CN201910362765 A CN 201910362765A CN 110113278 A CN110113278 A CN 110113278A
Authority
CN
China
Prior art keywords
signal
frequency
mtone
modulation
phase
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.)
Granted
Application number
CN201910362765.2A
Other languages
Chinese (zh)
Other versions
CN110113278B (en
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.)
Hangzhou Tianzhi Rongtong Technology Co ltd
Original Assignee
Hangzhou Dianzi University
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 Hangzhou Dianzi University filed Critical Hangzhou Dianzi University
Priority to CN201910362765.2A priority Critical patent/CN110113278B/en
Publication of CN110113278A publication Critical patent/CN110113278A/en
Application granted granted Critical
Publication of CN110113278B publication Critical patent/CN110113278B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/0012Modulated-carrier systems arrangements for identifying the type of modulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Digital Transmission Methods That Use Modulated Carrier Waves (AREA)

Abstract

The Modulation Mode Recognition method based on all-digital receiver that the invention discloses a kind of.The present invention estimates carrier frequency and three dB bandwidth the following steps are included: step 1, do FFT transform to the measurement and control signal that accepting device receives;Step 2 completes down coversion, filtering and resampling using the carrier frequency that estimation obtains;Step 3 passes through phase-locked loop tracking frequency offset to the low frequency signal after Orthogonal Decomposition, distinguishes fsk signal and AM_FM, MTONE_FM signal using frequency-tracking curve;Step 4 identifies AM_FM, MTONE_FM signal by doing FFT transform to the frequency-tracking curve for including modulation intelligence.The present invention has operand simple, easy to accomplish on hardware.The present invention can accurately be judged in all-digital receiver the narrowband FM signal of fsk signal and complex modulated that recognition correct rate is better than based on traditional signal transient frequency spectrum discerning method.

Description

一种基于全数字接收机的调制方式识别方法A Modulation Method Identification Method Based on All-Digital Receiver

技术领域technical field

本发明主要针对信号检测与频率估计,主要涉及一种基于全数字接收机的调制方式识别方法。The invention is mainly aimed at signal detection and frequency estimation, and mainly relates to a modulation mode identification method based on an all-digital receiver.

背景技术Background technique

在实际电子侦察环境中,发送的信号往往会伴有宽带数字调制信号和窄带信号混合出现的情况,这就需要对接收到的通信信号实现调制方式识别研究。并且在信号各种参数(信号带宽、调制方式)未知的情况下,使得信号侦察识别难度变大。In the actual electronic reconnaissance environment, the transmitted signal is often accompanied by a mixture of wideband digital modulation signal and narrowband signal, which requires research on the modulation method identification of the received communication signal. And in the case that the various parameters of the signal (signal bandwidth, modulation mode) are unknown, it makes the signal reconnaissance and identification more difficult.

目前采用到的调制方式识别方法主要有基于信号循环谱的识别方法研究,基于信号高阶累积量的识别方法研究,基于信号星座图识别方法,基于瞬时频率谱识别方法。基于瞬时频率谱识别方法虽然操作起来简单,但是不能对复合调制的窄带FM信号实现调制方式识别,并且容易受噪声等因素影响,在识别时往往要求较高的信噪比。基于信号星座图的调制方式识别方法,是一种PSK/QAM信号调制识别,要想实现调制方式识别,必须完成接收信号载波频率的准确估计以及信号完全同步。基于信号循环谱分析的通信信号调制识别研究,通过提取信号循环谱包络幅度的方差特征来分类识别以及使用高阶统计量特征这一方法,通过各阶累积量参数实现调制方式识别。这两种方法计算量大,参数提取计算过程难度高,信号的实时在线分析差,难以在实际工程项目中应用。The currently used modulation method identification methods mainly include the identification method research based on signal cyclic spectrum, the identification method research based on signal high-order cumulant, the identification method based on signal constellation diagram, and the identification method based on instantaneous frequency spectrum. Although the identification method based on the instantaneous frequency spectrum is simple to operate, it cannot realize the modulation mode identification of the composite modulated narrow-band FM signal, and is easily affected by factors such as noise, and often requires a higher signal-to-noise ratio for identification. The modulation identification method based on the signal constellation diagram is a kind of PSK/QAM signal modulation identification. In order to realize the modulation identification, the accurate estimation of the carrier frequency of the received signal and the complete synchronization of the signal must be completed. Research on communication signal modulation recognition based on signal cyclic spectrum analysis, by extracting the variance feature of the signal cyclic spectrum envelope amplitude to classify and identify, and using the method of high-order statistical features, the modulation mode identification is realized through the cumulant parameters of each order. These two methods have a large amount of calculation, the calculation process of parameter extraction is difficult, and the real-time online analysis of signals is poor, so it is difficult to apply in actual engineering projects.

基于以上分析可知,目前的调制识别方法主要针对PSK/QAM信号,对于窄宽带模数混合信号的识别,目前的识别方法识别性能快速下降。并且目前的识别算法计算量大,在硬件上实现困难。Based on the above analysis, it can be seen that the current modulation recognition methods are mainly aimed at PSK/QAM signals, and for the recognition of narrow-bandwidth analog-digital mixed signals, the recognition performance of the current recognition methods drops rapidly. Moreover, the current recognition algorithm has a large amount of calculation and is difficult to implement on hardware.

发明内容Contents of the invention

本发明的目的是针对现有技术的不足,实现FSK、AM_FM、MTONE_FM的信号侦察。本发明提出了一种基于全数字接收机的调制方式识别方法。本发明算法具有更好的识别性能,并且能够区分出复合调制的窄带FM信号,且算法的复杂度不高,能够在硬件上实现。The purpose of the present invention is to realize the signal detection of FSK, AM_FM and MTONE_FM aiming at the deficiencies of the prior art. The invention proposes a modulation mode identification method based on an all-digital receiver. The algorithm of the invention has better identification performance, and can distinguish composite modulated narrow-band FM signals, and the complexity of the algorithm is not high, and can be realized on hardware.

本发明解决其技术问题所采用的技术方案包括以下步骤:The technical solution adopted by the present invention to solve its technical problems comprises the following steps:

步骤1、对接受设备收到的测控信号做FFT变换,估计载频以及3dB带宽;Step 1. Perform FFT transformation on the measurement and control signal received by the receiving device, and estimate the carrier frequency and 3dB bandwidth;

步骤2、利用估计得到的载频完成下变频、滤波和重采样;Step 2, using the estimated carrier frequency to complete down-conversion, filtering and re-sampling;

步骤3、对正交分解后的低频信号通过锁相环路跟踪频偏,利用频率跟踪曲线区分出FSK信号和AM_FM、MTONE_FM信号;Step 3, track the frequency offset through the phase-locked loop for the low-frequency signal after the orthogonal decomposition, and use the frequency tracking curve to distinguish the FSK signal from the AM_FM and MTONE_FM signals;

步骤4、通过对包含有调制信息的频率跟踪曲线做FFT变换识别出AM_FM、MTONE_FM信号。Step 4. Identify the AM_FM and MTONE_FM signals by performing FFT transformation on the frequency tracking curve containing the modulation information.

步骤1具体实现如下:Step 1 is specifically implemented as follows:

1-1.对接收到的信号做1024点FFT变换,找出幅度最大值P对应的横坐标值,计算出信号频率f01-1. Perform 1024-point FFT transformation on the received signal, find out the abscissa value corresponding to the maximum amplitude P, and calculate the signal frequency f 0 ;

1-2.找到最大值点处所对应频率的底噪大小SNRest,由此计算出3dB带宽的门限值,高于门限值部分对应的频率范围即为信号的带宽threshold:1-2. Find the noise floor SNR est of the frequency corresponding to the maximum point, and calculate the threshold value of the 3dB bandwidth. The frequency range corresponding to the part higher than the threshold value is the bandwidth threshold of the signal:

步骤2具体实现如下:The specific implementation of step 2 is as follows:

2-1.根据计算得到的信号频率f0,对信号进行下变频、低通滤波、CIC抽取以降低数据速率;2-1. According to the calculated signal frequency f 0 , the signal is down-converted, low-pass filtered, and CIC extracted to reduce the data rate;

此过程表示为:This process is expressed as:

根据接收信号频率初估计的结果对信号进行混频,低通滤波滤除和频分量,得到FM复合调制信号同相信号和正交信号可表示为:According to the result of the initial estimation of the frequency of the received signal, the signal is mixed, and the sum frequency component is filtered out by low-pass filtering, and the in-phase signal and quadrature signal of the FM composite modulation signal are obtained, which can be expressed as:

2-2.根据3dB带宽,区别出窄带(AM_FM、MTONE_FM)信号和宽带FSK信号。若是窄带信号需要对采样速率进行CIC500倍抽取来降低采样速率,若是宽带信号则需要对采样速率进行CIC2倍抽取。2-2. According to the 3dB bandwidth, distinguish between narrowband (AM_FM, MTONE_FM) signals and wideband FSK signals. If it is a narrowband signal, it needs to decimate the sampling rate by 500 times of CIC to reduce the sampling rate, and if it is a broadband signal, it needs to decimate the sampling rate by 2 times of CIC.

步骤3具体实现如下:Step 3 is specifically implemented as follows:

3-1.正交分解后的信号和本地压控振荡器的初始频率相乘,低通滤波,得到同相分量和正交分量。此时信号I路和Q路可表示为:3-1. The quadrature decomposed signal is multiplied by the initial frequency of the local voltage-controlled oscillator, and low-pass filtered to obtain the in-phase component and the quadrature component. At this time, the signal I and Q channels can be expressed as:

式中,TS为信号的采样间隔,Kf为调制频偏,Δf为经过下变频后的信号和压控振荡器初始频差,m(t)为调制信号;In the formula, T S is the sampling interval of the signal, K f is the modulation frequency deviation, Δf is the initial frequency difference between the down-converted signal and the VCO, and m(t) is the modulation signal;

fc表示信号的载波频率,n表示采样点个数;f c represents the carrier frequency of the signal, and n represents the number of sampling points;

将这两路正交的信号作为鉴相器的输入信号。鉴相方程可表示为:Take these two quadrature signals as the input signals of the phase detector. The phase discrimination equation can be expressed as:

因为反正切与除法运算在硬件上难以实现,可以使用如下鉴相方程来代替:Because arctangent and division operations are difficult to implement on hardware, the following phase identification equation can be used instead:

3-2.鉴相方程得到相位误差信息,鉴相误差结果经过二阶环路滤波滤除和频分量,差频分量作为压控振荡器的输入电压,环路稳定后误差趋近于零,实现频率跟踪。3-2. The phase error information is obtained by the phase detection equation. The phase detection error result is filtered by the second-order loop to filter out the sum frequency component, and the difference frequency component is used as the input voltage of the voltage-controlled oscillator. After the loop is stable, the error approaches zero. Implement frequency tracking.

步骤4具体实现如下:Step 4 is specifically implemented as follows:

4-1.基于全数字锁相环的载波同步环路跟踪的是频率信息,AM、MTONE信号调制在频率上,因此频率跟踪曲线反映的是调制信息,通过对锁相环跟踪的含频率信息的信号做FFT变换,根据调制信息频谱特点不一样,可以实现调制方式识别。而FSK信号跟踪的频率会出现跳变,通过跟踪频率结合信号的频谱图可以区分出窄带(AM_FM、MTONE_FM)信号和宽带FSK信号。4-1. The carrier synchronization loop based on the all-digital phase-locked loop tracks the frequency information, and the AM and MTONE signals are modulated on the frequency, so the frequency tracking curve reflects the modulation information, and the frequency information contained in the tracking of the phase-locked loop The signal is transformed by FFT, and the modulation mode identification can be realized according to the different spectrum characteristics of the modulation information. The frequency tracked by the FSK signal will jump, and the narrowband (AM_FM, MTONE_FM) signal and the wideband FSK signal can be distinguished by tracking the frequency combined with the signal's spectrogram.

4-2.AM_FM信号经过载波环路跟踪得到频率信息,而AM信号调制在频率上,通过对频率信息做FFT运算,可以得到一个明显的谱峰值。而通过对MTONE_FM信号的频率信息做FFT运算,因为MTONE信号调制在频率上,从仿真结果可以看到有多个谱峰值。显然,可以通过搜索谱峰值的个数来区分出AM_FM信号和MTONE_FM信号。具体方法是计算出对应信号频谱的平均值,找出比平均值大的多的波峰值,从而确定频率冲激分量的个数。4-2. The frequency information of the AM_FM signal is tracked by the carrier loop, and the AM signal is modulated on the frequency. By performing FFT operation on the frequency information, an obvious spectral peak can be obtained. By performing FFT operation on the frequency information of the MTONE_FM signal, because the MTONE signal is modulated on the frequency, it can be seen from the simulation results that there are multiple spectral peaks. Obviously, the AM_FM signal and the MTONE_FM signal can be distinguished by searching for the number of spectral peaks. The specific method is to calculate the average value of the corresponding signal spectrum, and find out the peak value which is much larger than the average value, so as to determine the number of frequency impulse components.

本发明有益效果如下:The beneficial effects of the present invention are as follows:

1.本发明针对侦察信号参数未知的情况下,如载波频率,调制方式等都是未知的。可以实现调制方式识别。1. The present invention is aimed at the situation where the reconnaissance signal parameters are unknown, such as carrier frequency, modulation mode, etc. are all unknown. Modulation identification can be realized.

2.本发明通过锁相环跟踪的频率信息含有调制信息,实现复合调制信号(AM_FM、MTONE_FM)的调制方式识别,识别率高于常用的基于信号瞬时频谱的识别方法。2. The present invention contains modulation information through the frequency information tracked by the phase-locked loop, and realizes the modulation mode identification of the composite modulation signal (AM_FM, MTONE_FM), and the identification rate is higher than the commonly used identification method based on the instantaneous frequency spectrum of the signal.

3.本发明计算量简单,易于硬件实现,在实际的工程项目中可以很好的使用,有较高的现实意义。3. The present invention has a simple calculation amount, is easy to realize by hardware, can be well used in actual engineering projects, and has high practical significance.

附图说明Description of drawings

图1本发明中基于载波同步环路提取用户信息的调制识别方法框图;Fig. 1 is a block diagram of a modulation identification method based on a carrier synchronization loop to extract user information in the present invention;

图2 3dB带宽估计流程图;Figure 2 3dB bandwidth estimation flow chart;

图3 FSK信号载波同步频偏跟踪曲线;Figure 3 FSK signal carrier synchronization frequency offset tracking curve;

图4 AM_FM信号载波同步跟踪的频率信息;Figure 4 Frequency information of AM_FM signal carrier synchronous tracking;

图5 MTONE_FM信号载波同步跟踪的频率信息;Fig. 5 Frequency information of MTONE_FM signal carrier synchronous tracking;

图6 AM_FM信号跟踪的包含有调制信息的频率信号FFT频谱图;Fig. 6 FFT spectrum diagram of frequency signal containing modulation information tracked by AM_FM signal;

图7 MTONE_FM信号跟踪的包含有调制信息的频率信号FFT频谱图;Fig. 7 FFT spectrum diagram of frequency signal containing modulation information tracked by MTONE_FM signal;

具体实施方式Detailed ways

下面结合附图对本发明的具体实施例做进一步的说明。Specific embodiments of the present invention will be further described below in conjunction with the accompanying drawings.

如图1所示,一种基于全数字接收机的调制模式识别方法,具体实现过程如下:As shown in Figure 1, a modulation mode recognition method based on an all-digital receiver, the specific implementation process is as follows:

步骤1、对接受设备收到的测控信号做FFT变换,估计载频以及3dB带宽,具体实现流程图如图2所示;Step 1. Perform FFT transformation on the measurement and control signal received by the receiving device, and estimate the carrier frequency and 3dB bandwidth. The specific implementation flow chart is shown in Figure 2;

步骤2、利用估计得到的频率完成下变频、滤波和重采样等;Step 2, using the estimated frequency to complete down-conversion, filtering and resampling, etc.;

步骤3、对正交分解后的低频信号通过锁相环路跟踪频偏,利用频率跟踪曲线区分出FSK信号和AM_FM、MTONE_FM信号分别如图2,图3,图4所示;Step 3, track the frequency offset through the phase-locked loop for the low-frequency signal after the orthogonal decomposition, and use the frequency tracking curve to distinguish the FSK signal from the AM_FM and MTONE_FM signals as shown in Figure 2, Figure 3, and Figure 4 respectively;

步骤4、通过对包含有调制信息的频率跟踪曲线做FFT变换识别出AM_FM、MTONE_FM信号,如图5、6所示。Step 4. Identify the AM_FM and MTONE_FM signals by performing FFT transformation on the frequency tracking curve containing modulation information, as shown in Figures 5 and 6.

步骤1具体实现如下:Step 1 is specifically implemented as follows:

1-1.对接收到的信号做1024点FFT变换,找出幅度最大值P对应的横坐标值,计算出载频fc1-1. Perform 1024-point FFT transformation on the received signal, find out the abscissa value corresponding to the maximum amplitude value P, and calculate the carrier frequency f c ;

1-2.找到最大值点处对应频率的底噪大小SNRest,3dB带宽估计流程图如图2所示。由此可以计算出3dB带宽的门限值,高于门限值部分对应的频率范围即为信号的带宽;1-2. Find the noise floor SNR est of the corresponding frequency at the maximum point, and the flow chart of 3dB bandwidth estimation is shown in Figure 2. From this, the threshold value of the 3dB bandwidth can be calculated, and the frequency range corresponding to the part higher than the threshold value is the bandwidth of the signal;

步骤2具体实现如下:The specific implementation of step 2 is as follows:

2-1.根据计算得到的载频,对信号做正交分解变为I、Q两路,通过低通滤波滤除混频后的信号,再经过CIC抽取以降低数据速率;2-1. According to the calculated carrier frequency, the signal is orthogonally decomposed into two channels of I and Q, and the mixed signal is filtered out by low-pass filtering, and then extracted by CIC to reduce the data rate;

此过程可以表示为:This process can be expressed as:

根据接收信号频率初估计的结果对信号进行混频后通过低通滤波器滤除产生的和频分量,得到FM复合调制信号同相和正交信号可表示为:According to the result of the initial estimation of the frequency of the received signal, the signal is mixed and the sum frequency component is filtered out by a low-pass filter, and the in-phase and quadrature signals of the FM composite modulation signal can be expressed as:

2-2.根据3dB带宽,区别出窄带(AM_FM、MTONE_FM)信号和宽带FSK信号。若是窄带信号需要对采样速率进行CIC500倍抽取来降低采样速率,若是宽带信号则需要对采样速率进行CIC2倍抽取。2-2. According to the 3dB bandwidth, distinguish between narrowband (AM_FM, MTONE_FM) signals and wideband FSK signals. If it is a narrowband signal, it needs to decimate the sampling rate by 500 times of CIC to reduce the sampling rate, and if it is a broadband signal, it needs to decimate the sampling rate by 2 times of CIC.

步骤3具体实现如下:The specific implementation of step 3 is as follows:

3-1.正交分解后的信号和本地压控振荡器的初始频率相乘后得到同相分量和正交分量。此时信号I路和Q路可表示为:3-1. The quadrature-decomposed signal is multiplied by the initial frequency of the local voltage-controlled oscillator to obtain an in-phase component and a quadrature component. At this time, the signal I channel and Q channel can be expressed as:

式中,TS为信号的采样间隔,Kf为调制频偏,Δf表示经过下变频后的信号和压控振荡器初始频差,m(t)为调制信号。In the formula, T S is the sampling interval of the signal, K f is the modulation frequency deviation, Δf is the initial frequency difference between the down-converted signal and the VCO, and m(t) is the modulation signal.

将这两路正交的信号作为鉴相器的输入信号。鉴相方程可表示为:Take these two quadrature signals as the input signals of the phase detector. The phase discrimination equation can be expressed as:

因为反正切与除法运算在硬件上难以实现,可以使用如下鉴相方程来代替:Because arctangent and division operations are difficult to implement on hardware, the following phase identification equation can be used instead:

3-2.鉴相方程得到相位误差信息,鉴相误差结果经过二阶环路滤波滤除和频分量,差频分量作为压控振荡器的输入电压,环路稳定后误差趋近于零,实现频率跟踪。3-2. The phase error information is obtained by the phase detection equation. The phase detection error result is filtered by the second-order loop to filter out the sum frequency component, and the difference frequency component is used as the input voltage of the voltage-controlled oscillator. After the loop is stable, the error approaches zero. Implement frequency tracking.

步骤4具体实现如下:Step 4 is specifically implemented as follows:

4-1.基于全数字锁相环的载波同步环路跟踪的是频率信息,如图3,图4,图5所示。AM、MTONE信号调制在频率上,因此频率跟踪曲线反映的是调制信息,通过对锁相环跟踪的含频率信息的信号做FFT变换,根据调制信息频谱特点不一样,可以实现调制方式识别。而FSK信号跟踪的频率会出现跳变,通过跟踪频率结合信号的频谱图可以区分出窄带(AM_FM、MTONE_FM)信号和宽带FSK信号。4-1. The carrier synchronization loop based on the all-digital phase-locked loop tracks frequency information, as shown in Figure 3, Figure 4, and Figure 5. AM and MTONE signals are modulated on frequency, so the frequency tracking curve reflects the modulation information. By performing FFT transformation on the signal containing frequency information tracked by the phase-locked loop, the modulation mode identification can be realized according to the different spectrum characteristics of the modulation information. The frequency tracked by the FSK signal will jump, and the narrowband (AM_FM, MTONE_FM) signal and the wideband FSK signal can be distinguished by tracking the frequency combined with the signal's spectrogram.

4-2.AM_FM信号经过载波环路跟踪得到频率信息,而AM信号调制在频率上,通过对频率信息做FFT运算,可以得到一个明显的谱峰值,如图6所示。而通过对MTONE_FM信号的频率信息做FFT运算,因为MTONE信号调制在频率上,从仿真结果可以看到有多个谱峰值,如图7所示。显然,可以通过搜索谱峰值的个数来区分出AM_FM信号和MTONE_FM信号。具体方法是计算出对应信号频谱的平均值,找出比平均值大的多的波峰值,从而确定频率冲激分量的个数。4-2. The frequency information of the AM_FM signal is tracked by the carrier loop, and the AM signal is modulated on the frequency. By performing FFT operation on the frequency information, an obvious spectral peak can be obtained, as shown in Figure 6. By performing FFT operation on the frequency information of the MTONE_FM signal, because the MTONE signal is modulated on the frequency, it can be seen from the simulation results that there are multiple spectral peaks, as shown in Figure 7. Obviously, the AM_FM signal and the MTONE_FM signal can be distinguished by searching for the number of spectral peaks. The specific method is to calculate the average value of the corresponding signal spectrum, and find out the peak value which is much larger than the average value, so as to determine the number of frequency impulse components.

实施例:Example:

信号为中国电子科技集团公司利用Agilent E4438C发送已调信号,并用采样频率为200M的AD采样板采集AM_FM、MTONE_FM、FSK信号,其中FSK信号的带宽为6.4MHz,AM_FM、MTONE_FM信号的带宽为200kHz。The signals are modulated signals sent by China Electronics Technology Group Corporation using Agilent E4438C, and AM_FM, MTONE_FM, and FSK signals are collected by an AD sampling board with a sampling frequency of 200M. The bandwidth of the FSK signal is 6.4MHz, and the bandwidth of the AM_FM and MTONE_FM signals is 200kHz.

先将输入信号做FFT变换,得到信号的频谱图,估计出信号的载频并找出3dB带宽,区别出窄带和宽带信号。接下来将估计的载频作为下变频频率,完成下变频,CIC抽取和低通滤波。若是窄带信号,需要对采样频率进行CIC500倍抽取来降低采样速率,送入到载波同步环路中。若是宽带信号需要进行CIC2倍抽取,送入到载波同步环路中跟踪频率信息。在调制方式识别中,可以依据带宽和跟踪的频率信息是否跳变来区分出FSK信号和AM_FM、MTONE_FM信号。因为跟踪的频率曲线反映的是调制信号的信息,可以通过对包含有调制信息的频率信息做FFT变换,来区分出复合调制的窄带FM信号。通过测试结果得出,本发明的方法测试结果在信噪比为2dB、4dB和6dB时,识别率分别对应95.5%、98.5%和100。而目前常用的瞬时频谱调制方式识别方法很难将这两种信号区分出来。Firstly, the input signal is transformed by FFT to obtain the spectrum diagram of the signal, the carrier frequency of the signal is estimated and the 3dB bandwidth is found, and the narrowband and broadband signals are distinguished. Next, the estimated carrier frequency is used as the down-conversion frequency to complete the down-conversion, CIC extraction and low-pass filtering. If it is a narrowband signal, the sampling frequency needs to be decimated by CIC500 times to reduce the sampling rate and sent to the carrier synchronization loop. If the broadband signal needs to be extracted twice by CIC, it is sent to the carrier synchronization loop to track the frequency information. In the identification of the modulation mode, the FSK signal and the AM_FM and MTONE_FM signals can be distinguished according to whether the bandwidth and the tracked frequency information jump. Because the tracked frequency curve reflects the information of the modulated signal, the narrowband FM signal of complex modulation can be distinguished by performing FFT transformation on the frequency information containing the modulated information. According to the test results, the test results of the method of the present invention correspond to 95.5%, 98.5% and 100 respectively when the signal-to-noise ratio is 2dB, 4dB and 6dB. However, it is difficult to distinguish these two signals by the commonly used instantaneous spectrum modulation identification method at present.

实施例结果表明,本专利的基于全数字接收机的窄宽带信号调制方式识别算法具有运算量简单,易于硬件实现,仿真运行速度快。本专利算法相较于瞬时谱特征识别方法,可以识别出复合调制的窄带FM信号,具有比较高的现实意义。The result of the embodiment shows that the identification algorithm of the narrow-bandwidth signal modulation mode based on the all-digital receiver of this patent has the advantages of simple computation, easy hardware implementation, and fast simulation operation speed. Compared with the instantaneous spectrum feature recognition method, the patented algorithm can recognize composite modulated narrow-band FM signals, and has relatively high practical significance.

最后需要注意的是,公布实施例的目的在于帮助进一步理解本发明,但是本领域的技术人员可以理解:在不脱离本发明及所附的权利要求的精神和范围内,各种替换和修改都是可能的。因此,本发明不应局限于实施例所公开的内容,本发明要求保护的范围以权利要求书界定的范围为准。Finally, it should be noted that the purpose of the disclosed embodiments is to help further understand the present invention, but those skilled in the art can understand that various replacements and modifications can be made without departing from the spirit and scope of the present invention and the appended claims. It is possible. Therefore, the present invention should not be limited to the content disclosed in the embodiments, and the protection scope of the present invention is subject to the scope defined in the claims.

Claims (1)

1.一种基于全数字接收机的调制方式识别方法,其特征在于包括以下步骤:1. A modulation method identification method based on an all-digital receiver, characterized in that it may further comprise the steps: 步骤1、对接受设备收到的测控信号做FFT变换,估计载频以及3dB带宽;Step 1. Perform FFT transformation on the measurement and control signal received by the receiving device, and estimate the carrier frequency and 3dB bandwidth; 步骤2、利用估计得到的载频完成下变频、滤波和重采样;Step 2, using the estimated carrier frequency to complete down-conversion, filtering and re-sampling; 步骤3、对正交分解后的低频信号通过锁相环路跟踪频偏,利用频率跟踪曲线区分出FSK信号和AM_FM、MTONE_FM信号;Step 3, track the frequency offset through the phase-locked loop for the low-frequency signal after the orthogonal decomposition, and use the frequency tracking curve to distinguish the FSK signal from the AM_FM and MTONE_FM signals; 步骤4、通过对包含有调制信息的频率跟踪曲线做FFT变换识别出AM_FM、MTONE_FM信号;Step 4, identify the AM_FM and MTONE_FM signals by performing FFT transformation on the frequency tracking curve containing the modulation information; 步骤1具体实现如下:Step 1 is specifically implemented as follows: 1-1.对接收到的信号做1024点FFT变换,找出幅度最大值P对应的横坐标值,估计出信号频率f01-1. Perform 1024-point FFT transformation on the received signal, find out the abscissa value corresponding to the maximum amplitude P, and estimate the signal frequency f 0 ; 1-2.找到最大值点处所对应频率的底噪大小SNRest,由此计算出3dB带宽的门限值,高于门限值部分对应的频率范围即为信号的带宽threshold:1-2. Find the noise floor SNR est of the frequency corresponding to the maximum point, and calculate the threshold value of the 3dB bandwidth. The frequency range corresponding to the part higher than the threshold value is the bandwidth threshold of the signal: 步骤2具体实现如下:Step 2 is specifically implemented as follows: 2-1.根据计算得到的信号频率f0,对信号进行下变频、低通滤波、CIC抽取以降低数据速率;2-1. According to the calculated signal frequency f 0 , the signal is down-converted, low-pass filtered, and CIC extracted to reduce the data rate; 此过程表示为:This process is expressed as: 根据接收信号频率初估计的结果对信号进行混频,低通滤波滤除和频分量,得到FM复合调制信号同相信号和正交信号可表示为:According to the result of the initial estimation of the frequency of the received signal, the signal is mixed, and the sum frequency component is filtered out by low-pass filtering, and the in-phase signal and quadrature signal of the FM composite modulation signal are obtained, which can be expressed as: 2-2.根据3dB带宽,区别出窄带(AM_FM、MTONE_FM)信号和宽带FSK信号;若是窄带信号需要对采样速率进行CIC500倍抽取来降低采样速率,若是宽带信号则需要对采样速率进行CIC2倍抽取;2-2. According to the 3dB bandwidth, distinguish the narrowband (AM_FM, MTONE_FM) signal and the wideband FSK signal; if the narrowband signal needs to decimate the sampling rate by 500 times of CIC to reduce the sampling rate, if it is a broadband signal, it needs to decimate the sampling rate by 2 times of CIC ; 步骤3具体实现如下:Step 3 is specifically implemented as follows: 3-1.正交分解后的信号和本地压控振荡器的初始频率相乘,低通滤波,得到同相分量和正交分量;此时信号I路和Q路可表示为:3-1. The signal after quadrature decomposition is multiplied by the initial frequency of the local voltage-controlled oscillator, and low-pass filtered to obtain the in-phase component and quadrature component; at this time, the signal I and Q can be expressed as: 式中,TS为信号的采样间隔,Kf为调制频偏,Δf为经过下变频后的信号和压控振荡器初始频差,m(t)为调制信号;In the formula, T S is the sampling interval of the signal, K f is the modulation frequency deviation, Δf is the initial frequency difference between the down-converted signal and the VCO, and m(t) is the modulation signal; fc表示信号的载波频率,n表示采样点个数;f c represents the carrier frequency of the signal, and n represents the number of sampling points; 将这两路正交的信号作为鉴相器的输入信号;鉴相方程可表示为:The two orthogonal signals are used as the input signal of the phase detector; the phase detection equation can be expressed as: 因为反正切与除法运算在硬件上难以实现,可以使用如下鉴相方程来代替:Because arctangent and division operations are difficult to implement on hardware, the following phase identification equation can be used instead: 3-2.鉴相方程得到相位误差信息,鉴相误差结果经过二阶环路滤波滤除和频分量,差频分量作为压控振荡器的输入电压,环路稳定后误差趋近于零,实现频率跟踪;3-2. The phase error information is obtained by the phase detection equation. The phase detection error result is filtered by the second-order loop to filter out the sum frequency component, and the difference frequency component is used as the input voltage of the voltage-controlled oscillator. After the loop is stable, the error approaches zero. Realize frequency tracking; 步骤4具体实现如下:Step 4 is specifically implemented as follows: 4-1.基于全数字锁相环的载波同步环路跟踪的是频率信息,AM、MTONE信号调制在频率上,因此频率跟踪曲线反映的是调制信息,通过对锁相环跟踪的含频率信息的信号做FFT变换,根据调制信息频谱特点不一样,可以实现调制方式识别;而FSK信号跟踪的频率会出现跳变,通过跟踪频率结合信号的频谱图可以区分出窄带(AM_FM、MTONE_FM)信号和宽带FSK信号;4-1. The carrier synchronization loop based on the all-digital phase-locked loop tracks the frequency information, and the AM and MTONE signals are modulated on the frequency, so the frequency tracking curve reflects the modulation information, and the frequency information contained in the tracking of the phase-locked loop According to the different characteristics of the modulation information spectrum, the modulation mode identification can be realized; while the frequency tracked by the FSK signal will jump, and the narrowband (AM_FM, MTONE_FM) signal can be distinguished by tracking the frequency combined with the signal spectrum diagram. Broadband FSK signal; 4-2.AM_FM信号经过载波环路跟踪得到频率信息,而AM信号调制在频率上,通过对频率信息做FFT运算,可以得到一个明显的谱峰值;而通过对MTONE_FM信号的频率信息做FFT运算,因为MTONE信号调制在频率上,从仿真结果可以看到有多个谱峰值;显然,可以通过搜索谱峰值的个数来区分出AM_FM信号和MTONE_FM信号;具体方法是计算出对应信号频谱的平均值,找出比平均值大的多的波峰值,从而确定频率冲激分量的个数。4-2. The frequency information of the AM_FM signal is tracked by the carrier loop, and the AM signal is modulated on the frequency. By performing FFT operation on the frequency information, an obvious spectral peak can be obtained; and by performing FFT operation on the frequency information of the MTONE_FM signal , because the MTONE signal is modulated on the frequency, it can be seen from the simulation results that there are multiple spectral peaks; obviously, the AM_FM signal and the MTONE_FM signal can be distinguished by searching the number of spectral peaks; the specific method is to calculate the average of the corresponding signal spectrum value, find the peak value that is much larger than the average value, so as to determine the number of frequency impulse components.
CN201910362765.2A 2019-04-30 2019-04-30 A Modulation Identification Method Based on All-Digital Receiver Active CN110113278B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910362765.2A CN110113278B (en) 2019-04-30 2019-04-30 A Modulation Identification Method Based on All-Digital Receiver

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910362765.2A CN110113278B (en) 2019-04-30 2019-04-30 A Modulation Identification Method Based on All-Digital Receiver

Publications (2)

Publication Number Publication Date
CN110113278A true CN110113278A (en) 2019-08-09
CN110113278B CN110113278B (en) 2021-06-29

Family

ID=67487897

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910362765.2A Active CN110113278B (en) 2019-04-30 2019-04-30 A Modulation Identification Method Based on All-Digital Receiver

Country Status (1)

Country Link
CN (1) CN110113278B (en)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111245756A (en) * 2020-02-19 2020-06-05 杭州电子科技大学 Modulation identification method of composite signal based on cascaded SVM and all-digital receiver
CN111343116A (en) * 2020-02-25 2020-06-26 杭州电子科技大学 Modulation identification method based on constellation trajectory map in dense signal environment
CN112162153A (en) * 2020-08-31 2021-01-01 南京亿杰明信息技术有限公司 Frequency Estimation Method of Sine Wave Signal Based on Phase Line Fitting
CN112422465A (en) * 2019-10-09 2021-02-26 上海矢元电子有限公司 Signal modulation identification equipment
CN112565127A (en) * 2020-11-26 2021-03-26 杭州电子科技大学 Measurement and control signal modulation identification method based on residual error network and frequency domain graph
CN113014523A (en) * 2021-02-25 2021-06-22 电子科技大学 Composite modulation signal blind identification method based on fusion characteristics
CN116566782A (en) * 2023-05-22 2023-08-08 杭州电子科技大学 Continuous phase signal modulation identification method based on constellation-like diagram characteristics
CN117201249A (en) * 2023-11-01 2023-12-08 中孚安全技术有限公司 Signal modulation mode identification method, system and device
CN116566782B (en) * 2023-05-22 2025-07-04 杭州电子科技大学 A method for continuous phase signal modulation recognition based on constellation-like features

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2365923C1 (en) * 2008-07-03 2009-08-27 Федеральное государственное образовательное учреждение высшего профессионального образования "Военный авиационный инженерный университет" (г. Воронеж) Министерства обороны Российской Федерации Device for recognition of radio signals
US9071497B2 (en) * 2012-09-06 2015-06-30 Electronics And Telecommunications Research Institute Apparatus and method for modulation classification in wireless communication system
CN108134753A (en) * 2016-12-01 2018-06-08 南京理工大学 The recognition methods of broadband signal modulation system
CN108206797A (en) * 2016-12-20 2018-06-26 天津大学(青岛)海洋工程研究院有限公司 A kind of modulation mode of communication signal self-adaptive identification method based on software radio
CN108650203A (en) * 2018-07-23 2018-10-12 杭州电子科技大学 A kind of signal type Identification method based on reconnaissance receiver
CN109560794A (en) * 2018-11-09 2019-04-02 南京邮电大学 Vector Modulation radio-frequency signal generator

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2365923C1 (en) * 2008-07-03 2009-08-27 Федеральное государственное образовательное учреждение высшего профессионального образования "Военный авиационный инженерный университет" (г. Воронеж) Министерства обороны Российской Федерации Device for recognition of radio signals
US9071497B2 (en) * 2012-09-06 2015-06-30 Electronics And Telecommunications Research Institute Apparatus and method for modulation classification in wireless communication system
CN108134753A (en) * 2016-12-01 2018-06-08 南京理工大学 The recognition methods of broadband signal modulation system
CN108206797A (en) * 2016-12-20 2018-06-26 天津大学(青岛)海洋工程研究院有限公司 A kind of modulation mode of communication signal self-adaptive identification method based on software radio
CN108650203A (en) * 2018-07-23 2018-10-12 杭州电子科技大学 A kind of signal type Identification method based on reconnaissance receiver
CN109560794A (en) * 2018-11-09 2019-04-02 南京邮电大学 Vector Modulation radio-frequency signal generator

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
D. BOUDREAU等: ""A fast automatic modulation recognition algorithm and its implementation in a spectrum monitoring application"", 《MILCOM 2000 PROCEEDINGS. 21ST CENTURY MILITARY COMMUNICATIONS. ARCHITECTURES AND TECHNOLOGIES FOR INFORMATION SUPERIORITY》 *
张鸣: ""软件无线电接收机中信号调制样式识别算法的研究"", 《中国优秀硕士学位论文全文数据库》 *
李庆岩: ""无线信号调制模式识别技术的研究"", 《中国优秀硕士学位论文全文数据库》 *
范小骞: ""电子对抗中通信信号数字调制自动识别的研究"", 《中国优秀硕士学位论文全文数据库》 *
陈宇杨: ""一种基于数字接收机的信号调制方式识别技术"", 《软件导刊》 *

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112422465B (en) * 2019-10-09 2021-10-22 上海矢元电子股份有限公司 Signal modulation identification equipment
CN112422465A (en) * 2019-10-09 2021-02-26 上海矢元电子有限公司 Signal modulation identification equipment
CN111245756A (en) * 2020-02-19 2020-06-05 杭州电子科技大学 Modulation identification method of composite signal based on cascaded SVM and all-digital receiver
CN111245756B (en) * 2020-02-19 2022-07-22 杭州电子科技大学 Modulation identification method of composite signal based on cascaded SVM and all-digital receiver
CN111343116B (en) * 2020-02-25 2022-06-10 杭州电子科技大学 Modulation identification method based on constellation locus diagram under dense signal environment
CN111343116A (en) * 2020-02-25 2020-06-26 杭州电子科技大学 Modulation identification method based on constellation trajectory map in dense signal environment
CN112162153A (en) * 2020-08-31 2021-01-01 南京亿杰明信息技术有限公司 Frequency Estimation Method of Sine Wave Signal Based on Phase Line Fitting
CN112565127A (en) * 2020-11-26 2021-03-26 杭州电子科技大学 Measurement and control signal modulation identification method based on residual error network and frequency domain graph
CN113014523A (en) * 2021-02-25 2021-06-22 电子科技大学 Composite modulation signal blind identification method based on fusion characteristics
CN113014523B (en) * 2021-02-25 2022-01-11 电子科技大学 A Blind Recognition Method of Composite Modulated Signal Based on Fusion Features
CN116566782A (en) * 2023-05-22 2023-08-08 杭州电子科技大学 Continuous phase signal modulation identification method based on constellation-like diagram characteristics
CN116566782B (en) * 2023-05-22 2025-07-04 杭州电子科技大学 A method for continuous phase signal modulation recognition based on constellation-like features
CN117201249A (en) * 2023-11-01 2023-12-08 中孚安全技术有限公司 Signal modulation mode identification method, system and device
CN117201249B (en) * 2023-11-01 2024-02-20 中孚安全技术有限公司 Signal modulation mode identification method, system and device

Also Published As

Publication number Publication date
CN110113278B (en) 2021-06-29

Similar Documents

Publication Publication Date Title
CN110113278A (en) A kind of Modulation Mode Recognition method based on all-digital receiver
CN109495410B (en) High dynamic PCM/FM signal carrier frequency accurate estimation method
JP3271504B2 (en) Frequency estimation circuit and AFC circuit using the same
CN104022981B (en) A kind of blind carrier frequency bias estimation of quadrature amplitude modulation signal
CN107911329A (en) A kind of signal analyzer ofdm signal demodulation method
CN103856431B (en) The lower open-loop carrier synchronized algorithm of MPSK modulation based on fast Fourier transform
CN105450380B (en) Non-cooperating communication carrier synchronization system based on FFT carrier frequency estimation and Costas ring
CN111935046B (en) Low-complexity frequency shift keying signal symbol rate estimation method
CN116232823B (en) Amplitude-phase modulation signal identification method based on combination of multi-stage blind digital receivers
CN109067680A (en) A kind of carrier frequency bias estimation and its device of baseband signal
CN102571652A (en) Method for estimating global navigation satellite system (GNSS) baseband signal
CN107786479A (en) A kind of big frequency deviation compensation system of QPSK carrier auxiliaries and method
JP2002504774A (en) Quadrature-free RF receiver that receives angle-modulated signals directly
CN114531329B (en) A multi-channel MSK signal carrier frequency estimation method, system and application
CN106526631B (en) A high-precision method for estimating carrier frequency of Beidou satellite signal at B1 frequency point
CN106027454A (en) Frequency offset estimation method of QAM (Quadrature Amplitude Modulation) signal based on frequency compensation
US8149956B1 (en) Method of automated demodulation and classification of phase-shift-keying signals using hysteretic differential zero-crossing time samples
CN111343116A (en) Modulation identification method based on constellation trajectory map in dense signal environment
CN104199063B (en) A kind of blind frequency discriminator processing method based on cross product algorithm
CN107528805B (en) A PSK signal synchronization method and device suitable for signal analyzer
JP2020010195A (en) Frequency estimation device and tracking receiver
CN106341123B (en) A filtering method and device for single-tone interference
CN111245580A (en) A system and method for calculating signal-to-noise ratio based on hardware logic circuit
CN109521269A (en) A kind of am signals digitlization frequency measuring method
CN111245756B (en) Modulation identification method of composite signal based on cascaded SVM and all-digital receiver

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right

Effective date of registration: 20210908

Address after: 310018 room 412, No. 68, Keye Road, Jianggan District, Hangzhou City, Zhejiang Province

Patentee after: Hangzhou Tianzhi RONGTONG Technology Co.,Ltd.

Address before: 310018 No. 2 street, Xiasha Higher Education Zone, Hangzhou, Zhejiang

Patentee before: HANGZHOU DIANZI University

TR01 Transfer of patent right