CN111711586B - 一种应用于通信调制方式识别的时频协同滤波方法及系统 - Google Patents

一种应用于通信调制方式识别的时频协同滤波方法及系统 Download PDF

Info

Publication number
CN111711586B
CN111711586B CN202010540939.2A CN202010540939A CN111711586B CN 111711586 B CN111711586 B CN 111711586B CN 202010540939 A CN202010540939 A CN 202010540939A CN 111711586 B CN111711586 B CN 111711586B
Authority
CN
China
Prior art keywords
frequency
analytic signal
time
instantaneous
signal
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
Application number
CN202010540939.2A
Other languages
English (en)
Other versions
CN111711586A (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.)
Jilin University
Original Assignee
Jilin 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 Jilin University filed Critical Jilin University
Priority to CN202010540939.2A priority Critical patent/CN111711586B/zh
Publication of CN111711586A publication Critical patent/CN111711586A/zh
Application granted granted Critical
Publication of CN111711586B publication Critical patent/CN111711586B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00—Modulated-carrier systems
    • H04L27/0012—Modulated-carrier systems arrangements for identifying the type of modulation
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00—Baseband systems
    • H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/03—Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
    • H04L25/03006—Arrangements for removing intersymbol interference
    • H04L25/03159—Arrangements for removing intersymbol interference operating in the frequency domain

Landscapes

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

Abstract

本发明涉及一种应用于通信调制方式识别的时频协同滤波方法及系统，所述滤波方法包括：确定解析信号；对所述解析信号进行瞬时频率估计，得到所述解析信号的瞬时频率；将所述解析信号输入至下变频混频器，混频器频率采用信号的瞬时频率，得到下变频混频器输出数据序列；将所述变换后的下变频混频器输出数据序列输入至低通滤波器，得到滤波后的数据序列；将滤波后的数据序列输入至上变频混频器，得到协同滤波后的高信噪比解析信号。本发明中的上述方法能够有效滤除干扰噪声，提升信噪比。

Description

一种应用于通信调制方式识别的时频协同滤波方法及系统
技术领域
本发明涉及通信调制领域，特别是涉及一种应用于通信调制方式识别的时频协同滤波方法及系统。
背景技术
通信调制方式识别是在未知调制信息内容的前提下，判断出通信信号的调制方式，并估计出相应的调制参数。在民用方面，相关职能部门可以通过分析接收到的无线电信号识别出未注册的发射机；在军事方面，通过对截获信号调制方式的识别与分类，有针对性地进行信号解调，以便实现对敌方通信内容的侦听，或对敌方通信实施欺骗式干扰。
通信调制方式识别系统的构成主要有信号接收机、数据预处理、特征提取、分类识别和显示。预处理的目的是削弱噪声的影响，加强有用的信息特征；特征提取是通过对预处理数据的分析与变换，得到最能反映不同调制类型本质的特征；分类识别是利用基于特征空间的信息对待识别对象分类进行。
当通信信号信噪比较低时，通信调制方式识别率会受到严重影响。解决该问题的可行策略是对接收原始数据进行数字滤波处理，减小噪声对特征提取的影响。
频率调制信号广泛应用于无线电通讯领域，典型应用如线性调频通信，非线性调频通信(跳频通信和正弦调频)及各种频率调制通信信号。频率调制信号属于典型非平稳信号，其频率随时间变化。基于平稳信号处理的滤波方法不能有效滤除干扰噪声。时频分析是分析和处理非平稳信号的有力工具，它将信号表示为时间和频率的联合函数，清楚描述信号频率随时间变化的关系。平稳信号的滤波方法，通常采用时域或频域进行滤波，滤波器设有起始频率和截止频率，因此滤波器带宽固定。由于频率调制信号的频率是变化的，若带通滤波器通带带宽设置过小，虽然能够大大减小白噪声噪声影响，但可能滤除有用信号；若带通滤波器通带带宽过大，则会使更多噪声通过滤波器，不利于信噪比提升。
发明内容
本发明的目的是提供一种应用于通信调制方式识别的时频协同滤波方法及系统，有效滤除干扰噪声，并提升信噪比。
为实现上述目的，本发明提供了如下方案：
一种应用于通信调制方式识别的时频协同滤波方法，所述滤波方法包括：
确定解析信号；所述解析信号为含噪声的解析信号；
对所述解析信号进行迭代的瞬时频率估计，得到所述解析信号的准确的瞬时频率；所述瞬时频率用于后续上变频混频器和下变频混频器的工作频率；
将所述解析信号输入至下变频混频器，得到下变频混频器输出数据序列；所述下变频混频器的工作频率为所述解析信号的准确的瞬时频率；
将所述下变频混频器输出数据序列输入至低通滤波器，得到滤波后的数据序列；
将滤波后的数据序列输入至上变频混频器，得到协同滤波后的高信噪比解析信号；所述上变频混频器的工作频率为所述解析信号的准确的瞬时频率。
可选的，所述对所述解析信号进行迭代的瞬时频率估计，得到所述解析信号的准确的瞬时频率具体包括：
步骤1：对所述解析信号进行Wigner-Ville分布表示，并确定时频谱；
步骤2：对每个时间点按频率分量取解析信号的能量最大点对应频率值作为该时刻频率估计值，得到解析信号瞬时频率；
步骤3：基于所述估计的瞬时频率确定待估解析信号的瞬时相位；
步骤4：对所述瞬时相位进行重构，得到估计的解析信号；
步骤5：对所述估计的解析信号进行Wigner-Ville分布表示，得到新的时频谱图；
步骤6：基于所述新的时频谱图，采用步骤2方法求取其瞬时频率，判断本次瞬时频率和上次估计瞬时频率之差在固定间隔内的变化量的和是否达到最小值；
步骤7：若达到最小值，则停止运算，输出此时的瞬时频率，否则返回所述步骤3。
可选的，所述解析信号具体采用以下公式表示：
Figure BDA0002538922350000031
其中，s(t)为频率调制信号，a(t)为信号时变幅度，
Figure BDA0002538922350000032
为信号相位，n(t)为高斯白噪声。
可选的，所述对所述解析信号进行Wigner-Ville分布表示，并确定时频谱具体采用以下公式：
Figure BDA0002538922350000033
其中，Wn(t,f)表示时间t为横轴，频率为纵轴的Wigner-Ville时频分布；
Figure BDA0002538922350000034
表示时刻t时延为τ/2的解析信号；
Figure BDA0002538922350000035
表示时刻t，时延为-τ/2解析信号的复共轭，f表示瞬时频率，t表示时间，n表示迭代运算时第n次Wigner-Ville时频分布；e-j2πft表示核函数，τ表示时间延迟。
可选的，所述对每个时间点按频率分量取解析信号的能量最大点对应频率值作为该时刻频率估计值，得到解析信号瞬时频率具体采用以下公式：
Figure BDA0002538922350000036
其中，
Figure BDA0002538922350000037
表示解析信号的瞬时频率估计，f表示瞬时频率，Wn(t,f)表示第n次Wigner-Ville时频分布，n表示迭代运算时第n次Wigner-Ville时频分布。
所述基于所述估计的瞬时频率确定待估解析信号的瞬时相位具体采用以下公式：
Figure BDA0002538922350000038
其中，
Figure BDA0002538922350000039
表示解析信号的瞬时相位，
Figure BDA00025389223500000310
表示时变信号的瞬时频率估计值。
可选的，对所述瞬时相位进行重构，得到估计的解析信号具体采用以下公式：
Figure BDA00025389223500000311
其中，
Figure BDA00025389223500000312
表示估计的解析信号，a(t)表示瞬时幅度，
Figure BDA00025389223500000313
表示单位解析信号的复数域表示，
Figure BDA00025389223500000314
表示瞬时相位，j表示虚部。
可选的，对所述估计的解析信号进行Wignal-Ville分布表示，得到新的时频谱图具体采用以下公式：
Figure BDA0002538922350000041
其中，Wn+1(t,f)表示第n+1次Wigner-Ville分布，
Figure BDA0002538922350000042
表示第n+1次解析信号估计，时间为t，时延为τ/2，
Figure BDA0002538922350000043
表示第n+1次解析信号估计值的复共轭，e-j2πft表示核函数。
可选的，所述本次瞬时频率和上次估计瞬时频率之差在固定间隔内的变化量的和具体采用以下公式：
Figure BDA0002538922350000044
其中，
Figure BDA0002538922350000045
表示第n次瞬时频率估计，
Figure BDA0002538922350000046
表示第n+1次瞬时频率估计。
本发明另外提供一种应用于通信调制方式识别的时频协同滤波系统，所述系统包括：
解析信号确定模块，用于确定解析信号；所述解析信号为含噪声的解析信号；
瞬时频率估计模块，用于对所述解析信号进行迭代的瞬时频率估计，得到所述解析信号的准确的瞬时频率；所述瞬时频率用于后续上变频混频器和下变频混频器的工作频率；
第一输入模块，用于将所述解析信号输入至下变频混频器，得到下变频混频器输出数据序列；所述下变频混频器的工作频率为所述解析信号的准确的瞬时频率；
滤波模块，用于将所述下变频混频器输出数据序列输入至低通滤波器，得到滤波后的数据序列；
第二输入模块，用于将滤波后的数据序列输入至上变频混频器，得到协同滤波后的高信噪比解析信号；所述上变频混频器的工作频率为所述解析信号的准确的瞬时频率。
根据本发明提供的具体实施例，本发明公开了以下技术效果：
本发明中的上述方法通过对解析信号进行瞬时频率估计，得到所述解析信号的瞬时频率，并作为混频器工作频率；将解析信号输入至下变频混频器，得到下变频混频器输出数据序列；将变换后的下变频混频器输出数据序列输入至低通滤波器，得到滤波后的数据序列；将滤波后的数据序列输入至上变频混频器，得到协同滤波后的数据序列，可有效滤除干扰噪声，保留有用信号，提升信噪比。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例中心频率可变的带通滤波器原理示意图；
图2为本发明实施例应用于通信调制方式识别的时频协同滤波方法流程图；
图3为本发明实施例时变带通滤波器结构原理图；
图4为本发明实施例x(t)频谱图；
图5为本发明实施例y(t)频谱图；
图6为本发明实施例w(t)频谱图；
图7为本发明实施例s(t)频谱图；
图8为本发明实施例带通滤波器示意图；
图9为本发明实施例用于通信调制方式识别的时频协同滤波系统结构示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
本发明的目的是提供一种应用于通信调制方式识别的时频协同滤波方法及系统，有效滤除干扰噪声，并提升信噪比。
为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。
本发明拟采用一种中心频率可调的自适应带通滤波器，主要考虑接收信号为非平稳信号，即信号频率随时间变化，而干扰噪声为高斯白噪声情况。虽然调频信号在长时间内占用较宽频带，但在确定时刻或一小段时间内该信号仍属窄带信号，而噪声一直是宽谱噪声，为此设计一个根据信号瞬时频率自动调节中心频率的带通滤波器，滤除瞬时信号频率以外的宽谱噪声，随着中心频率的改变依次形成一系列的带通滤波器，其原理示意图如图1所示。以线性调频信号为例，线形调频信号持续时间t1—t3内频率由f1变为f3，若时间t1到t2较短，可以在t1-t2时间内设计一个带通滤波器，其通带频率为f1-f2，这样在该段时间内仅有f1到f2间较小通带内的噪声保留，其他频段噪声滤除；接着在t2到t3时间内设计通带为f2到f3的带通滤波器，实际系统可能时间分段更多。这样的时变带通滤波器要优于固定带宽(f1到f3通带)滤波器。
图2为本发明实施例应用于通信调制方式识别的时频协同滤波方法流程图，图3为本发明实施例时变带通滤波器结构原理图，如图2和图3所示，所述方法包括：
步骤101：确定解析信号。
其中，解析信号即为待滤波非平稳信号，具体表达式如下：
Figure BDA0002538922350000061
其中，s(t)为频率调制信号，a(t)为信号时变幅度，
Figure BDA0002538922350000062
为信号相位，n(t)为高斯白噪声。
若输入的信号为实信号x(t)，则需要经过希尔伯特变换，形成实信号x(t)的虚部j*HT[x(t)]，信号实部和虚部相加构成复信号输入，
Figure BDA0002538922350000064
步骤102：对所述解析信号进行迭代的瞬时频率估计，得到所述解析信号的准确的瞬时频率；所述瞬时频率用于后续上变频混频器和下变频混频器的工作频率。
具体包括以下步骤：
步骤1：对所述解析信号进行Wigner-Ville分布表示，并确定时频谱。
具体采用以下公式：
Figure BDA0002538922350000063
步骤2：对每个时间点按频率分量取解析信号的时频谱的能量最大值点对应的频率值作为当前时刻频率估计值，得到解析信号瞬时频率估计。
具体采用以下公式：
Figure BDA0002538922350000071
步骤3：基于所述估计的瞬时频率确定估计解析信号的瞬时相位。
具体采用以下公式：
Figure BDA0002538922350000072
步骤4：对所述瞬时相位进行重构，得到估计的解调信号。
具体采用以下公式：
Figure BDA0002538922350000073
步骤5：对所述估计的解析信号进行Wigner-Ville分布表示，得到新的时频谱图，并利用步骤2估计其瞬时频率
Figure BDA0002538922350000074
具体采用以下公式：
Figure BDA0002538922350000075
步骤6：基于所述新的时频谱图，采用步骤2方法求取其瞬时频率，判断本次瞬时频率和上次估计瞬时频率之差在固定间隔内的变化量的和是否达到最小值。
所述瞬时频率差在固定间隔内的变化量的和具体采用以下公式：
Figure BDA0002538922350000076
其中
Figure BDA0002538922350000077
为瞬时频率在离散时间点t处的估计频率值，N为分析数据离散化后长度。
步骤7：若达到最小值，则停止运算，输出此段数据序列的瞬时频率，否则返回所述步骤2。
实际程序运行中，若Var值出现极小值点或采用两次Var的相对变化率小于5％
Figure BDA0002538922350000078
即可终止迭代。
步骤103：将所述解析信号输入至下变频混频器，得到下变频混频器输出数据序列y(t)；所述下变频混频器的工作频率为所述解析信号的准确的瞬时频率。
其中，用公式表示为：
Figure BDA0002538922350000079
混频器角频率为
Figure BDA00025389223500000710
步骤104：将所述变换后的下变频混频器输出数据序列输入至低通滤波器，得到滤波后的数据序列，滤波后的数据序列频谱w(t)表示。
步骤105：将滤波后的数据序列输入至上变频混频器，得到协同滤波后的信噪比高的解析信号；所述上变频混频器的工作频率为所述解析信号的准确的瞬时频率。
带通滤波工作原理：下变频混频器输出频谱Y(jω)(如图5所示)就等于x(t)的频谱(如图4所示)向左移动ωc，X(ω)中ω＝ωc附近的频谱(信号的主要能量)就移进低通滤波器(截止频率为ω0)的通带内，而宽谱噪声的多数频谱能量被低通滤波器去除。低通滤波器的输出w(t)(如图6所示)进一步与复指数正弦
Figure BDA0002538922350000081
相乘(相当于上变频混频)，其输出信号
Figure BDA0002538922350000082
的傅立叶变换为
Figure BDA0002538922350000083
由傅立叶变换频移特性可知
Figure BDA0002538922350000084
(如图7所示)，就是W(jω)向右移ωc。至此原始信号的中心频率又移到原来的位置，这相当于x(t)经过带宽为2ω0的理想带通滤波器得到调频信号s(t)的近似信号
Figure BDA0002538922350000085
而噪声大大减弱。这样仅通过设计合理的低通滤波器(其截止频率ω0可设定为信号瞬时频率值的5％，即ω0＝5％ωc)就能构建与带通滤波器(通带频率为[ωc-ω0，ωc+ω0]，如图8所示同样效果。
图9为本发明实施例用于通信调制方式识别的时频协同滤波系统结构示意图，如图9所示，所述系统包括：解析信号确定模块201、瞬时频率估计模块202、第一输入模块203、滤波模块204以及第二输入模块205。
解析信号确定模块201，用于确定解析信号；所述解析信号为含噪声的解析信号。
瞬时频率估计模块202，用于对所述解析信号进行迭代的瞬时频率估计，得到所述解析信号的准确的瞬时频率；所述瞬时频率用于后续上变频混频器和下变频混频器的工作频率。
第一输入模块203，用于将所述解析信号输入至下变频混频器，得到下变频混频器输出数据序列；所述下变频混频器的工作频率为所述解析信号的准确的瞬时频率。
滤波模块204，用于将所述变换后的下变频混频器输出数据序列输入至低通滤波器，得到滤波后的数据序列。
第二输入模块205，用于将滤波后的数据序列输入至上变频混频器，得到协同滤波后的高信噪比解析信号；所述上变频混频器的工作频率为所述解析信号的准确的瞬时频率。
本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处。综上所述，本说明书内容不应理解为对本发明的限制。

Claims (9)

1.一种应用于通信调制方式识别的时频协同滤波方法，其特征在于，所述滤波方法包括：
确定解析信号；所述解析信号为含噪声的解析信号；
对所述解析信号进行迭代的瞬时频率估计，得到所述解析信号的准确的瞬时频率；所述瞬时频率用于后续上变频混频器和下变频混频器的工作频率；所述对所述解析信号进行迭代的瞬时频率估计，得到所述解析信号的准确的瞬时频率具体包括：
步骤1：对所述解析信号进行Wigner-Ville分布表示，并确定时频谱；
步骤2：对每个时间点按频率分量取解析信号的能量最大点对应频率值作为该时刻频率估计值，得到解析信号瞬时频率；
步骤3：基于所述估计的瞬时频率确定待估解析信号的瞬时相位；
步骤4：对所述瞬时相位进行重构，得到估计的解析信号；
步骤5：对所述估计的解析信号进行Wigner-Ville分布表示，得到新的时频谱图；
步骤6：基于所述新的时频谱图，采用步骤2方法求取其瞬时频率，判断本次瞬时频率和上次估计瞬时频率之差在固定间隔内的变化量的和是否达到最小值；
步骤7：若达到最小值，则停止运算，输出此时的瞬时频率，否则返回所述步骤3；
将所述解析信号输入至下变频混频器，得到下变频混频器输出数据序列；所述下变频混频器的工作频率为所述解析信号的准确的瞬时频率；
将所述下变频混频器输出数据序列输入至低通滤波器，得到滤波后的数据序列；
将滤波后的数据序列输入至上变频混频器，得到协同滤波后的高信噪比解析信号；所述上变频混频器的工作频率为所述解析信号的准确的瞬时频率。
2.根据权利要求1所述的应用于通信调制方式识别的时频协同滤波方法，其特征在于，所述解析信号具体采用以下公式表示：
Figure FDA0003046203970000011
其中，s(t)为频率调制信号，a(t)为信号时变幅度，
Figure FDA0003046203970000012
为信号相位，n(t)为高斯白噪声。
3.根据权利要求1所述的应用于通信调制方式识别的时频协同滤波方法，其特征在于，所述对所述解析信号进行Wigner-Ville分布表示，并确定时频谱具体采用以下公式：
Figure FDA0003046203970000021
其中Wn(t,f)表示时间t为横轴，频率为纵轴的Wigner-Ville时频分布；
Figure FDA0003046203970000022
表示时刻t时延为τ/2的解析信号；
Figure FDA0003046203970000023
表示时刻t，时延为-τ/2解析信号的复共轭，f表示瞬时频率，t表示时间，n表示迭代运算时第n次Wigner-Ville时频分布；e-j2πft表示核函数，τ表示时间延迟。
4.根据权利要求1所述的应用于通信调制方式识别的时频协同滤波方法，其特征在于，所述对每个时间点按频率分量取解析信号的能量最大点对应频率值作为该时刻频率估计值，得到解析信号瞬时频率具体采用以下公式：
Figure FDA0003046203970000024
其中，
Figure FDA0003046203970000025
表示解析信号的瞬时频率估计，f表示瞬时频率，Wn(t,f)表示第n次Wigner-Ville时频分布，n表示迭代运算时第n次Wigner-Ville时频分布。
5.根据权利要求1所述的应用于通信调制方式识别的时频协同滤波方法，其特征在于，所述基于所述估计的瞬时频率确定待估解析信号的瞬时相位具体采用以下公式：
Figure FDA0003046203970000026
其中，
Figure FDA0003046203970000027
表示解析信号的瞬时相位，
Figure FDA0003046203970000028
表示时变信号的瞬时频率估计值，τ表示时间延迟。
6.根据权利要求1所述的应用于通信调制方式识别的时频协同滤波方法，其特征在于，对所述瞬时相位进行重构，得到估计的解析信号具体采用以下公式：
Figure FDA0003046203970000029
其中，
Figure FDA00030462039700000210
表示估计的解析信号，a(t)表示瞬时幅度，
Figure FDA00030462039700000211
表示单位解析信号的复数域表示，
Figure FDA00030462039700000212
表示瞬时相位，j表示虚部。
7.根据权利要求1所述的应用于通信调制方式识别的时频协同滤波方法，其特征在于，对所述估计的解析信号进行Wigner-Ville分布表示，得到新的时频谱图具体采用以下公式：
Figure FDA00030462039700000213
其中，Wn+1(t,f)表示第n+1次Wigner-Ville分布，
Figure FDA00030462039700000214
表示第n+1次解析信号估计，时间为t，时延为τ/2，
Figure FDA00030462039700000215
表示第n+1次解析信号估计值的复共轭，e-j2πft表示核函数。
8.根据权利要求1所述的应用于通信调制方式识别的时频协同滤波方法，其特征在于，所述本次瞬时频率和上次估计瞬时频率之差在固定间隔内的变化量的和具体采用以下公式：
Figure FDA0003046203970000031
其中，
Figure FDA0003046203970000032
表示第n次瞬时频率估计，
Figure FDA0003046203970000033
表示第n+1次瞬时频率估计。
9.一种应用于通信调制方式识别的时频协同滤波系统，其特征在于，所述系统包括：
解析信号确定模块，用于确定解析信号；所述解析信号为含噪声的解析信号；
瞬时频率估计模块，用于对所述解析信号进行迭代的瞬时频率估计，得到所述解析信号的准确的瞬时频率；所述瞬时频率用于后续上变频混频器和下变频混频器的工作频率；所述对所述解析信号进行迭代的瞬时频率估计，得到所述解析信号的准确的瞬时频率具体包括：
步骤1：对所述解析信号进行Wigner-Ville分布表示，并确定时频谱；
步骤2：对每个时间点按频率分量取解析信号的能量最大点对应频率值作为该时刻频率估计值，得到解析信号瞬时频率；
步骤3：基于所述估计的瞬时频率确定待估解析信号的瞬时相位；
步骤4：对所述瞬时相位进行重构，得到估计的解析信号；
步骤5：对所述估计的解析信号进行Wigner-Ville分布表示，得到新的时频谱图；
步骤6：基于所述新的时频谱图，采用步骤2方法求取其瞬时频率，判断本次瞬时频率和上次估计瞬时频率之差在固定间隔内的变化量的和是否达到最小值；
步骤7：若达到最小值，则停止运算，输出此时的瞬时频率，否则返回所述步骤3；
第一输入模块，用于将所述解析信号输入至下变频混频器，得到下变频混频器输出数据序列；所述下变频混频器的工作频率为所述解析信号的准确的瞬时频率；
滤波模块，用于将所述下变频混频器输出数据序列输入至低通滤波器，得到滤波后的数据序列；
第二输入模块，用于将滤波后的数据序列输入至上变频混频器，得到协同滤波后的高信噪比解析信号；所述上变频混频器的工作频率为所述解析信号的准确的瞬时频率。
CN202010540939.2A 2020-06-15 2020-06-15 一种应用于通信调制方式识别的时频协同滤波方法及系统 Active CN111711586B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010540939.2A CN111711586B (zh) 2020-06-15 2020-06-15 一种应用于通信调制方式识别的时频协同滤波方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010540939.2A CN111711586B (zh) 2020-06-15 2020-06-15 一种应用于通信调制方式识别的时频协同滤波方法及系统

Publications (2)

Publication Number Publication Date
CN111711586A CN111711586A (zh) 2020-09-25
CN111711586B true CN111711586B (zh) 2021-07-13

Family

ID=72540050

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010540939.2A Active CN111711586B (zh) 2020-06-15 2020-06-15 一种应用于通信调制方式识别的时频协同滤波方法及系统

Country Status (1)

Country Link
CN (1) CN111711586B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114488105B (zh) * 2022-04-15 2022-08-23 四川锐明智通科技有限公司 一种基于运动特征及方向模板滤波的雷达目标检测方法
CN116626568B (zh) * 2023-07-19 2023-10-03 中国科学院精密测量科学与技术创新研究院 一种基于混频迭代的磁共振信号频率测量方法

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1866754A (zh) * 2005-05-17 2006-11-22 松下电器产业株式会社 具有镜像抑制混频器和带通滤波器的接收if系统
CN101610076A (zh) * 2009-07-21 2009-12-23 中国电子科技集团公司第四十一研究所 智能跟踪的可调谐带通滤波装置
CN103048488A (zh) * 2012-09-03 2013-04-17 中山大学 一种汽车加速度信号的消噪方法
CN108234370A (zh) * 2017-12-22 2018-06-29 西安电子科技大学 基于卷积神经网络的通信信号调制方式识别方法
US10432249B1 (en) * 2018-06-14 2019-10-01 The Boeing Company Universal nonlinear variable delay filter
CN110444843A (zh) * 2019-09-04 2019-11-12 上海海事大学 基于频率可调结构可切换的带通滤波器
CN111046791A (zh) * 2019-12-11 2020-04-21 河海大学 一种基于含可变因子广义s变换的电流信号滤波去噪方法
CN111092283A (zh) * 2020-01-03 2020-05-01 西安电子科技大学 一种具有可调陷波的超宽带带通滤波器及应用

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8040265B2 (en) * 2007-11-29 2011-10-18 Hrl Laboratories, Llc Down-converter and up-converter for time-encoded signals
US8090332B2 (en) * 2007-12-12 2012-01-03 Qualcomm, Incorporated Tracking filter for a receiver
CN102196111B (zh) * 2011-04-14 2014-05-07 西安烽火电子科技有限责任公司 一种短波两路调制解调器

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1866754A (zh) * 2005-05-17 2006-11-22 松下电器产业株式会社 具有镜像抑制混频器和带通滤波器的接收if系统
CN101610076A (zh) * 2009-07-21 2009-12-23 中国电子科技集团公司第四十一研究所 智能跟踪的可调谐带通滤波装置
CN103048488A (zh) * 2012-09-03 2013-04-17 中山大学 一种汽车加速度信号的消噪方法
CN108234370A (zh) * 2017-12-22 2018-06-29 西安电子科技大学 基于卷积神经网络的通信信号调制方式识别方法
US10432249B1 (en) * 2018-06-14 2019-10-01 The Boeing Company Universal nonlinear variable delay filter
CN110444843A (zh) * 2019-09-04 2019-11-12 上海海事大学 基于频率可调结构可切换的带通滤波器
CN111046791A (zh) * 2019-12-11 2020-04-21 河海大学 一种基于含可变因子广义s变换的电流信号滤波去噪方法
CN111092283A (zh) * 2020-01-03 2020-05-01 西安电子科技大学 一种具有可调陷波的超宽带带通滤波器及应用

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Synthesis and design of asymmetrical dual-band bandpass filters based on equivalent network simplification;Lenoir;《IEEE Transactions on Microwave Theory and Techniques》;20061231;全文 *
基于数字下变频的低通滤波器设计;周遐;《昆明冶金高等专科学校学报》;20101231;全文 *
用PSCF 设计自适应带通滤波器;马胜前;《兰州大学学报》;19991231;全文 *

Also Published As

Publication number Publication date
CN111711586A (zh) 2020-09-25

Similar Documents

Publication Publication Date Title
CN108469560B (zh) 一种基于快速s变换时频空间模型的电磁干扰客观复杂度评估方法
CN111711586B (zh) 一种应用于通信调制方式识别的时频协同滤波方法及系统
KR102169078B1 (ko) 광대역 측정 신호들을 측정하기 위한 방법 및 측정 장치
CN102546499B (zh) 一种实线性调频信号的分数阶信道化接收方法
US20200119759A1 (en) Noise suppression device, noise suppression method, and reception device and reception method using same
CN105158740A (zh) 基于高精度频率估计的噪声调幅干扰抑制方法
Kanaa et al. A robust parameter estimation of FHSS signals using time–frequency analysis in a non-cooperative environment
CN105163340A (zh) 一种适用于tdd-lte系统外干扰的分析方法
CN112003803B (zh) 一种vhf、uhf频段航空电台信号的侦收设备
Adjemov et al. Methods and algorithms of broadband HF signals dispersion distortion compensation
CN108900445B (zh) 一种信号符号率估计的方法及装置
CN111562438A (zh) 一种基于fft和相位差的正弦信号频率估计方法和装置
Sisi et al. Carrier frequency and symbol rate estimation based on cyclic spectrum
Liu et al. A novel signal separation algorithm for wideband spectrum sensing in cognitive networks
CN109104215A (zh) 一种基于小波变换的跳频信号码元速率盲估计方法
CN109076038B (zh) 用于对频带中包含的信号的参数进行估计的方法
CN110113273B (zh) 一种基于频谱对称的时频重叠信号盲分离方法
CN102790650B (zh) 一种宽带lfm信号的分数域奇偶交替信道化接收方法
CN109460614B (zh) 基于瞬时带宽的信号时间-频率分解方法
Chee et al. IF estimation of FSK signals using adaptive smoothed windowed cross Wigner–Ville distribution
CN104320360B (zh) 一种基于分数阶傅里叶变换的线性调频信号时延估计方法
Chee et al. Adaptive windowed cross Wigner–Ville distribution as an optimum phase estimator for PSK signals
KR101943734B1 (ko) 도약 주파수 추정 장치 및 방법
Djukanović et al. A parametric method for non-stationary interference suppression in direct sequence spread-spectrum systems
CN113206714B (zh) 一种基于误码率的接收机脉冲信号干扰评估方法

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