CN109962714B - 一种电磁频谱伞罩数字域自干扰抑制方法及系统 - Google Patents

一种电磁频谱伞罩数字域自干扰抑制方法及系统 Download PDF

Info

Publication number
CN109962714B
CN109962714B CN201910272810.5A CN201910272810A CN109962714B CN 109962714 B CN109962714 B CN 109962714B CN 201910272810 A CN201910272810 A CN 201910272810A CN 109962714 B CN109962714 B CN 109962714B
Authority
CN
China
Prior art keywords
signal
interference
self
umbrella cover
baseband
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
CN201910272810.5A
Other languages
English (en)
Other versions
CN109962714A (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.)
University of Electronic Science and Technology of China
Original Assignee
University of Electronic Science and Technology of China
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 University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN201910272810.5A priority Critical patent/CN109962714B/zh
Publication of CN109962714A publication Critical patent/CN109962714A/zh
Application granted granted Critical
Publication of CN109962714B publication Critical patent/CN109962714B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B1/00—Details of transmission systems, not covered by a single one of groups H04B3/00 - H04B13/00; Details of transmission systems not characterised by the medium used for transmission
    • H04B1/02—Transmitters
    • H04B1/03—Constructional details, e.g. casings, housings
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B1/00—Details of transmission systems, not covered by a single one of groups H04B3/00 - H04B13/00; Details of transmission systems not characterised by the medium used for transmission
    • H04B1/02—Transmitters
    • H04B1/04—Circuits
    • H04B1/0475—Circuits with means for limiting noise, interference or distortion
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B1/00—Details of transmission systems, not covered by a single one of groups H04B3/00 - H04B13/00; Details of transmission systems not characterised by the medium used for transmission
    • H04B1/06—Receivers
    • H04B1/10—Means associated with receiver for limiting or suppressing noise or interference
    • H04B1/1027—Means associated with receiver for limiting or suppressing noise or interference assessing signal quality or detecting noise/interference for the received signal
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B1/00—Details of transmission systems, not covered by a single one of groups H04B3/00 - H04B13/00; Details of transmission systems not characterised by the medium used for transmission
    • H04B1/06—Receivers
    • H04B1/10—Means associated with receiver for limiting or suppressing noise or interference
    • H04B1/12—Neutralising, balancing, or compensation arrangements
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00—Modulated-carrier systems
    • H04L27/0014—Carrier regulation
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00—Modulated-carrier systems
    • H04L27/0014—Carrier regulation
    • H04L2027/0024—Carrier regulation at the receiver end
    • H04L2027/0026—Correction of carrier offset

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Noise Elimination (AREA)

Abstract

本发明公开了一种电磁频谱伞罩数字域自干扰抑制方法及系统,所述包括以下步骤：S1.伞罩内的通信发射机产生通信信号；S2.伞罩干扰机根据伞罩干扰基带序列产生自干扰信号；S3.伞罩自干扰抑制接收机对接收到的信号进行预处理，得到基带信号；S4.对基带信号进行频偏估计与补偿，并以频偏补偿后的信号为基础，从时延、幅度、相位三个维度，对本地参考干扰信号进行估计与补偿，完成自干扰信号重建；S5.根据本地干扰重建信号进行伞罩数字域自干扰抑制；S6.对抑制得到的信号进行解调处理，恢复出通信发射机发送的数据。本发明实现了对伞罩自干扰信号频偏、时延、幅度、相位四个参数的高精度跟踪与补偿，有效解决了干扰信号的时频幅相四维对齐的问题。

Description

一种电磁频谱伞罩数字域自干扰抑制方法及系统
技术领域
本发明涉及电磁频谱伞罩系统，特别是涉及一种电磁频谱伞罩数字域自干扰抑制方法及系统。
背景技术
在特定区域，同时存在通信、探测、干扰、指控、导航等多种无线电收发设备。当大功率电子系统开机工作时，产生的同时同频干扰信号的辐射功率远超过我方电台接收机的灵敏度，严重时直接阻塞接收机的设备前端，造成这些设备接收中断。
同时，物理层安全技术不断发展，其中一种基本研究思路就是以各种窃听信道模型下的保密容量分析为基础，如果主信道比窃听信道的可靠性更好，就可以确保在一定码率上的可靠通信。本电磁频谱伞罩数字域自干扰抑制方法就是基于此思路，通过发送人为的干扰信号，使其干扰非法用户的正常通信，并通过自干扰抑制方法降低干扰信号对合法用户的影响，从而提升系统安全性。
电磁频谱伞罩数字域自干扰抑制方法，充分利用已知的干扰特征信息，在本地重建和消除干扰信号的影响，保证己方的正常通信。其中，精确地重建和跟踪伞罩自干扰信号是实现干扰抑制的关键所在。
由于收发信机之间的相对运动和时钟抖动的影响，伞罩自干扰信号的传播时延、幅度、相位、收发频偏均随时间变化，时频幅相的精准对齐是实现射频域高干扰抑制比的前提，因此实现电磁频谱伞罩数字域自干扰抑制的关键是干扰信号的时频幅相四维对齐。现有的干扰信号四维对齐研究主要集中在同时同频全双工，其研究对象是同一通信设备的发射信号对接收信号的干扰，即本机自干扰。
分析同时同频全双工具有以下特点：收发信道相对固定；本机自干扰信号的传播时延、幅度、相位均为定值；不存在收发频差。因此同时同频全双工通信技术不能解决干扰信号传播时延、幅度、相位、收发频偏的时变问题。
发明内容
本发明的目的在于克服现有技术的不足，提供一种电磁频谱伞罩数字域自干扰抑制方法及系统，实现了对伞罩自干扰信号频偏、时延、幅度、相位四个参数的高精度跟踪与补偿，有效解决了干扰信号的时频幅相四维对齐的问题。
本发明的目的是通过以下技术方案来实现的：一种电磁频谱伞罩数字域自干扰抑制方法，包括以下步骤：
S1.伞罩内的通信发射机产生通信信号s(t)并进行发射；
S2.伞罩干扰机根据伞罩干扰基带序列yb(n)产生自干扰信号y(t)，并向特定区域进行发射，提高该区域电磁波干扰强度，形成电磁频谱伞罩；
S3.伞罩自干扰抑制接收机对接收到的信号r(t)进行预处理，得到基带信号；
S4.伞罩自干扰抑制接收机在数字基带，对接收到的信号r(t)所对应的基带信号进行频偏估计与补偿，并将本地已知的伞罩干扰基带序列yb(n)作为本地参考干扰信号，以频偏补偿后的信号为基础，从时延、幅度、相位三个维度，对本地参考干扰信号进行估计与补偿，完成自干扰信号重建，得到本地干扰重建信号；
S5.伞罩自干扰抑制接收机根据时频幅相补偿得到的本地干扰重建信号进行伞罩数字域自干扰抑制；
S6.伞罩自干扰抑制接收机对抑制得到的信号进行解调处理，恢复出通信发射机发送的数据。
进一步地，所述步骤S2包括：
S201.将伞罩干扰基带序列yb(n)进行上变频，得到：
Figure BDA0002018954960000021
其中，fc为发送信号的载波频率，Ts为信号采样率；
S202.对信号y(n)进行数模转换，得到连续的自干扰信号y(t)；
S203.将产生自干扰信号y(t)向特定区域进行发送，提高该区域电磁波干扰强度，形成电磁频谱伞罩。
进一步地，所述步骤S3包括：
伞罩自干扰抑制接收机对接收到的信号r(t)依次进行模数转换、下变频和成型滤波，得到对应的基带信号。
进一步地，所述步骤S4包括以下子步骤：
S401.经过无线信道后，设接收到的干扰信号部分为ry(t)，接收到的有用信号部分为rs(t)，瞬时加性高斯白噪声为n(t)，将伞罩接收机接收到的信号对应的基带信号表示为：
Figure BDA0002018954960000022
式中，A1为幅度衰减系数，τ为相应的时延，fr为伞罩中接收干扰信号的频率，θ为相位偏移，
Figure BDA0002018954960000031
为自干扰信号y(t)经过无线传输信道后的信道衰减系数，fΔ＝fr-fc为收发两端信号的频率偏差，rs+n(t)为接收信号中有用信号和噪声部分对应的基带信号；
S402.在数字基带，对接收到的信号r(t)进行频偏估计与补偿：
将伞罩自干扰抑制接收机接收到信号r(t)对应的复基带信号采样序列表示为r(n)：
Figure BDA0002018954960000032
其中，
Figure BDA0002018954960000033
为自干扰信号经无线传输信道后的信道衰减系数；rs+n(n)为接收信号中有用信号和噪声部分对应的复基带信号采样序列；
Figure BDA0002018954960000034
为伞罩干扰发射信号和干扰接收信号之间的时间差τ归一化分数时延；
对r(n)进行处理，推导出接收信号中干扰部分的频偏估计表达式为：
Figure BDA0002018954960000035
其中，m为每一帧中选取的导引序列长度，L为进行共轭相乘的两复信号项之间的间隔帧长；
根据频偏估计值进行频偏补偿，设
Figure BDA0002018954960000036
则接收信号的频偏已被完全补偿，补偿后的信号为：
Figure BDA0002018954960000037
其中，α·yb(n-D)为接收信号中频偏补偿后的干扰部分，
Figure BDA0002018954960000038
为接收信号中频偏补偿后的有用部分和噪声部分之和；
S403.对频偏补偿后的信号
Figure BDA0002018954960000039
进行符号级同步：将信号
Figure BDA00020189549600000310
与干扰发射信号的采样序列作滑动相关，通过相关峰的位置确定
Figure BDA00020189549600000311
序列中干扰部分的起始位置；
S404.同步至单个符号精度后，对
Figure BDA00020189549600000312
进行分数级时延估计，并根据时延估计值对干扰参考信号进行时延补偿：
若接收信号的频偏已被完全补偿，补偿后的信号为：
Figure BDA0002018954960000041
将频偏补偿后的接收信号
Figure BDA0002018954960000042
与本地参考干扰信号yref(n)＝yb(n)作FFT变换，从时域变换到频域，进行共轭相乘和幅度归一化处理，有：
Figure BDA0002018954960000043
其中，Yref(ω)为本地参考干扰信号yref(n)的FFT变换，
Figure BDA0002018954960000044
为接收信号
Figure BDA0002018954960000045
的FFT变换；
对
Figure BDA0002018954960000046
取相位，得到信号
Figure BDA0002018954960000047
所在频段ω的相位函数
Figure BDA0002018954960000048
该相位函数是一条相位直线方程，为：
Figure BDA0002018954960000049
其中，ωd＜ω＜ωu，其中ωd是信号的低端频率，ωu是信号的高端频率，该直线的斜率直接对应了两信号之间的时间差
Figure BDA00020189549600000410
即估计出信号yb(n)和信号
Figure BDA00020189549600000411
之间的时间差，为：
Figure BDA00020189549600000412
根据估计值，采用Farrow滤波器分数时延对接收信号的时延进行补偿，假设滤波器的频率响应为：
Figure BDA00020189549600000413
其中，
Figure BDA00020189549600000414
采用
Figure BDA00020189549600000415
的多项式来逼近，即：
Figure BDA00020189549600000416
使H(ejω)逼近
Figure BDA00020189549600000417
即设计一组滤波器系数Cn,m，使其满足：
Figure BDA00020189549600000418
其中，ω∈[ω0,ω1]表示数字频率上线性时延带宽，
Figure BDA00020189549600000419
表示分数时延参数的范围；
利用上述Farrow滤波器，对本地参考干扰信号进行时延补偿后，对应的基带信号为：
Figure BDA00020189549600000420
S405.对基带信号
Figure BDA00020189549600000421
到信号
Figure BDA00020189549600000422
之间的幅相衰减进行估计与补偿：
基于时延补偿后得到的本地参考干扰对应的基带信号为
Figure BDA0002018954960000051
则由最大似然估计，得出：
Figure BDA0002018954960000052
计算幅度衰减的估计值为：
Figure BDA0002018954960000053
计算相位衰减的估计值为：
Figure BDA0002018954960000054
对
Figure BDA0002018954960000055
进行幅相补偿，得到本地干扰重建信号为：
Figure BDA0002018954960000056
一种电磁频谱伞罩数字域自干扰抑制系统，包括通信发射机、伞罩干扰机和伞罩自干扰抑制接收机；
所述通信发射机，用于产生通信信号并进行发射；
所述伞罩干扰机，用于根据伞罩干扰基带序列产生自干扰信号，并向特定区域进行发射，提高该区域电磁波干扰强度，形成电磁频谱伞罩；
所述伞罩自干扰抑制接收机，用于接收电磁频谱伞罩中的信号变换到基带，在数字基带对基带信号进行频偏估计与补偿，以频偏补偿后的信号为基础，从时延、幅度、相位三个维度，对本地参考干扰信号进行估计与补偿，完成自干扰信号重建，得到本地干扰重建信号；根据本地干扰重建信号进行伞罩数字域自干扰抑制；然后对抑制得到的信号进行解调处理，恢复出通信发射机发送的数据。
其中，所述伞罩自干扰抑制接收机包括：信号接收转换单元，用于接收电磁频谱伞罩中的信号，对接收到的信号进行预处理，得到基带信号；自干扰重建单元，用于在数字基带对基带信号进行频偏估计与补偿，以频偏补偿后的信号为基础，从时延、幅度、相位三个维度，对本地参考干扰信号进行估计与补偿，完成自干扰信号重建，得到本地干扰重建信号；自干扰抑制单元，用于根据时频幅相补偿后的本地干扰重建信号进行伞罩数字域自干扰抑制；解调恢复单元，用于对抑制得到的信号进行解调处理，恢复出通信发射机发送的数据。
本发明的有益效果是：本发明实现了对伞罩自干扰信号频偏、时延、幅度、相位四个参数的高精度跟踪与补偿，有效解决了干扰信号的时频幅相四维对齐的问题，伞罩自干扰的抑制不受通信信号的类型制约，可以与任何已有的通信链路相结合，兼容性较强。
附图说明
图1为本发明的方法流程图；
图2为电磁频谱伞罩数字域自干扰抑制原理；
图3为电磁频谱伞罩数字域自干扰抑制实验验证环境；
图4为电磁频谱伞罩数字域自干扰抑制接收机频率环工作方法；
图5为电磁频谱伞罩数字域自干扰抑制接收机时间环工作方法；
图6为电磁频谱伞罩数字域自干扰抑制接收机Farrow滤波器结构；
图7为电磁频谱伞罩数字域自干扰抑制接收机幅相环工作方法；
图8为自干扰抑制前和抑制后我方接收机信号频谱；
图9为自干扰抑制前和抑制后我方接收机信号频谱。
具体实施方式
下面结合附图进一步详细描述本发明的技术方案，但本发明的保护范围不局限于以下所述。
如图1所示，一种电磁频谱伞罩数字域自干扰抑制方法，包括以下步骤：
S1.伞罩内的通信发射机产生通信信号s(t)并进行发射；
S2.伞罩干扰机根据伞罩干扰基带序列yb(n)产生自干扰信号y(t)，并向特定区域进行发射，提高该区域电磁波干扰强度，形成电磁频谱伞罩；
S3.伞罩自干扰抑制接收机对接收到的信号r(t)进行预处理，得到基带信号；
S4.伞罩自干扰抑制接收机在数字基带，对接收到的信号r(t)所对应的基带信号进行频偏估计与补偿，并将本地已知的伞罩干扰基带序列yb(n)作为本地参考干扰信号，以频偏补偿后的信号为基础，从时延、幅度、相位三个维度，对本地参考干扰信号进行估计与补偿，完成自干扰信号重建，得到本地干扰重建信号；
S5.伞罩自干扰抑制接收机根据时频幅相补偿得到的本地干扰重建信号进行伞罩数字域自干扰抑制；
S6.伞罩自干扰抑制接收机对抑制得到的信号进行解调处理，恢复出通信发射机发送的数据。
进一步地，所述步骤S2包括：
S201.将伞罩干扰基带序列yb(n)进行上变频，得到：
Figure BDA0002018954960000071
其中，fc为发送信号的载波频率，Ts为信号采样率；
S202.对信号y(n)进行数模转换，得到连续的自干扰信号y(t)；
S203.将产生自干扰信号y(t)向特定区域进行发送，提高该区域电磁波干扰强度，形成电磁频谱伞罩。
进一步地，所述步骤S3包括：
伞罩自干扰抑制接收机对接收到的信号r(t)依次进行模数转换、下变频和成型滤波，得到对应的基带信号。
进一步地，所述步骤S4包括以下子步骤：
S401.经过无线信道后，设接收到的干扰信号部分为ry(t)，接收到的有用信号部分为rs(t)，瞬时加性高斯白噪声为n(t)，将伞罩接收机接收到的信号对应的基带信号表示为：
Figure BDA0002018954960000072
式中，A1为幅度衰减系数，τ为相应的时延，fr为伞罩中接收干扰信号的频率，θ为相位偏移，
Figure BDA0002018954960000073
为自干扰信号y(t)经过无线传输信道后的信道衰减系数，fΔ＝fr-fc为收发两端信号的频率偏差，rs+n(t)为接收信号中有用信号和噪声部分对应的基带信号；
S402.在数字基带，对接收到的信号r(t)进行频偏估计与补偿：
将伞罩自干扰抑制接收机接收到信号r(t)对应的复基带信号采样序列表示为r(n)：
Figure BDA0002018954960000074
其中，
Figure BDA0002018954960000075
为自干扰信号经无线传输信道后的信道衰减系数；rs+n(n)为接收信号中有用信号和噪声部分对应的复基带信号采样序列；
Figure BDA0002018954960000076
为伞罩干扰发射信号和干扰接收信号之间的时间差τ归一化分数时延；
对r(n)进行处理，推导出接收信号中干扰部分的频偏估计表达式为：
Figure BDA0002018954960000081
其中，m为每一帧中选取的导引序列长度，L为进行共轭相乘的两复信号项之间的间隔帧长；
根据频偏估计值进行频偏补偿，设
Figure BDA0002018954960000082
则接收信号的频偏已被完全补偿，补偿后的信号为：
Figure BDA0002018954960000083
其中，α·yb(n-D)为接收信号中频偏补偿后的干扰部分，
Figure BDA0002018954960000084
为接收信号中频偏补偿后的有用部分和噪声部分之和；
S403.对频偏补偿后的信号
Figure BDA0002018954960000085
进行符号级同步：将信号
Figure BDA0002018954960000086
与干扰发射信号的采样序列作滑动相关，通过相关峰的位置确定
Figure BDA0002018954960000087
序列中干扰部分的起始位置；
S404.同步至单个符号精度后，对
Figure BDA0002018954960000088
进行分数级时延估计，并根据时延估计值对干扰参考信号进行时延补偿：
若接收信号的频偏已被完全补偿，补偿后的信号为：
Figure BDA0002018954960000089
将频偏补偿后的接收信号
Figure BDA00020189549600000810
与本地参考干扰信号yref(n)＝yb(n)作FFT变换，从时域变换到频域，进行共轭相乘和幅度归一化处理，有：
Figure BDA00020189549600000811
其中，Yref(ω)为本地参考干扰信号yref(n)的FFT变换，
Figure BDA00020189549600000812
为接收信号
Figure BDA00020189549600000813
的FFT变换；
对
Figure BDA00020189549600000814
取相位，得到信号
Figure BDA00020189549600000815
所在频段ω的相位函数
Figure BDA00020189549600000816
该相位函数是一条相位直线方程，为：
Figure BDA00020189549600000817
其中，ωd＜ω＜ωu，其中ωd是信号的低端频率，ωu是信号的高端频率，该直线的斜率直接对应了两信号之间的时间差
Figure BDA00020189549600000819
，即估计出信号yb(n)和信号
Figure BDA00020189549600000818
之间的时间差，为：
Figure BDA0002018954960000091
根据估计值，采用Farrow滤波器分数时延对接收信号的时延进行补偿，假设滤波器的频率响应为：
Figure BDA0002018954960000092
其中，
Figure BDA0002018954960000093
采用
Figure BDA0002018954960000094
的多项式来逼近，即：
Figure BDA0002018954960000095
使H(ejω)逼近
Figure BDA0002018954960000096
即设计一组滤波器系数Cn,m，使其满足：
Figure BDA0002018954960000097
其中，ω∈[ω0,ω1]表示数字频率上线性时延带宽，
Figure BDA0002018954960000098
表示分数时延参数的范围；
利用上述Farrow滤波器，对本地参考干扰信号进行时延补偿后，对应的基带信号为：
Figure BDA0002018954960000099
S405.对基带信号
Figure BDA00020189549600000910
到信号
Figure BDA00020189549600000911
之间的幅相衰减进行估计与补偿：
基于时延补偿后得到的本地参考干扰对应的基带信号为
Figure BDA00020189549600000912
则由最大似然估计，得出：
Figure BDA00020189549600000913
计算幅度衰减的估计值为：
Figure BDA00020189549600000914
计算相位衰减的估计值为：
Figure BDA00020189549600000915
对
Figure BDA00020189549600000916
进行幅相补偿，得到本地干扰重建信号为：
Figure BDA00020189549600000917
一种电磁频谱伞罩数字域自干扰抑制系统，包括通信发射机、伞罩干扰机和伞罩自干扰抑制接收机；
所述通信发射机，用于产生通信信号并进行发射；
所述伞罩干扰机，用于根据伞罩干扰基带序列产生自干扰信号，并向特定区域进行发射，提高该区域电磁波干扰强度，形成电磁频谱伞罩；
所述伞罩自干扰抑制接收机，用于接收电磁频谱伞罩中的信号变换到基带，在数字基带对基带信号进行频偏估计与补偿，以频偏补偿后的信号为基础，从时延、幅度、相位三个维度，对本地参考干扰信号进行估计与补偿，完成自干扰信号重建，得到本地干扰重建信号；根据本地干扰重建信号进行伞罩数字域自干扰抑制；然后对抑制得到的信号进行解调处理，恢复出通信发射机发送的数据。
其中，所述伞罩自干扰抑制接收机包括：信号接收转换单元，用于接收电磁频谱伞罩中的信号，对接收到的信号进行预处理，得到基带信号；自干扰重建单元，用于在数字基带对基带信号进行频偏估计与补偿，以频偏补偿后的信号为基础，从时延、幅度、相位三个维度，对本地参考干扰信号进行估计与补偿，完成自干扰信号重建，得到本地干扰重建信号；自干扰抑制单元，用于根据时频幅相补偿后的本地干扰重建信号进行伞罩数字域自干扰抑制；解调恢复单元，用于对抑制得到的信号进行解调处理，恢复出通信发射机发送的数据。
在本申请的实施例中，以FM信号为通信发射机、M序列扩频信号为干扰发射机为例，本发明的工作原理如图2所示。伞罩干扰机发射机、FM通信发射机通过天线分别辐射出去信号1、2，接收天线接收到合信号3，送至伞罩自干扰抑制接收机。自干扰抑制接收机充分利用已知的干扰特征信息，在本地重建和跟踪干扰信号，进行数字域干扰抵消，实现己方的正常通信，伞罩自干扰抑制实验验证环境如图3所示。
以图4所示的伞罩自干扰抑制接收机频率环工作方法为例，频率环包括频率估计和补偿，将接收信号4进行缓存，从每帧中选取导引序列开始与结束所处位置，设置对应间隔帧数，使这两帧导引序列对应共轭相乘，对乘积累加、取相位，得到频偏估计值，完成频偏对齐，得到频偏校正后的接收信号5，实现信号载波同步。
以图5所示的伞罩自干扰抑制接收机时间环工作方法为例，时间环包括时间估计和补偿，接收信号5和本地参考信号6通过相关完成符号同步，再作FFT变换，从时域变换到频域，进行共轭相乘和幅度归一化处理，得到一条相位直线方程，直线的斜率即表示干扰发射信号和接收信号之间的时间差。再通过Farrow滤波器分数延时完成时间对齐，得到重建信号7。
以图6所示的伞罩自干扰抑制接收机Farrow滤波器结构为例，Farrow结构有M组N阶FIR滤波器构成，其中虚线框图中具体画出了第M组N阶FIR滤波器(CM-1)的内部结构，其余M-1组N阶FIR滤波器的实现方式与CM-1相同。Farrow结构所需要的滤波器阶数远大于前两种滤波器结构，运算量很大，它的优点是当时延变化时，仅改变时延参数D，就可以获得不同的分数时延，不用重新加载系数，节省了存储空间，降低了硬件实现的复杂度。
以图7所示的伞罩自干扰抑制接收机幅相环工作方法为例，幅相环包括幅度估计和补偿，接收信号5与参考信号7进行共轭相乘和归一化处理，得到幅度估计值，完成幅度对齐，最后得到时频幅相四维对齐后的重建信号8。
为了对本发明提出的一种电磁频谱伞罩数字域自干扰抑制方法的性能进行验证，对电磁频谱伞罩实验设备参数设置如下：
Figure BDA0002018954960000111
图8给出了通信发射机保持静默时，伞罩自干扰抑制前和抑制后的我方接收机信号功率谱。当我方接收机的干噪比为30dB时，自干扰抑制度为28.5dB。自干扰抑制度基本不受自干扰发射功率的影响。
图9给出了伞罩内的通信发射机正常发射FM调制语音信号时，自干扰抑制前和抑制后的我方接收机信号功率谱。我方天线输出口的干信比为10dB。自干扰被抑制前，期望的FM调制信号被淹没在干扰信号频谱中，无法正常解调和播放语音信号；自干扰抑制后，能够看到期望的FM调制信号频谱轮廓，可以正常解调和播放语音信号。FM收音机一直无法正常解调和播放语音信号。
综上所述，考虑无线通信系统中伞罩数字域自干扰抑制的传输场景，针对加性白高斯噪声无线传播信道，在加性白高斯噪声无线自干扰信道中，本文提出的一种电磁频谱伞罩数字域自干扰抑制方法在指定区域可以干扰非法通信，同时根据干扰特征在本地重建并抵消干扰，以保证己方正常通信的集干扰与通信为一体。从仿真结果中容易看出：自干扰抑制前，通信发射机与我方接收机的通信被阻塞；自干扰抑制后，通信发射机与我方接收机可以正常进行通信。非法接收机一直不能够正常接收通信信号。该结果初步验证了本方法的可行性。
这里已经通过具体的实施例子对本发明进行了详细描述，提供上述实施例的描述为了使本领域的技术人员制造或适用本发明，这些实施例的各种修改对于本领域的技术人员来说是容易理解的。本发明并不限于这些例子，或其中的某些方面。本发明的范围通过附加的权利要求进行详细说明。
上述说明示出并描述了本发明的一个优选实施例，但如前所述，应当理解本发明并非局限于本文所披露的形式，不应看作是对其他实施例的排除，而可用于各种其他组合、修改和环境，并能够在本文所述发明构想范围内，通过上述教导或相关领域的技术或知识进行改动。而本领域人员所进行的改动和变化不脱离本发明的精神和范围，则都应在本发明所附权利要求的保护范围内。

Claims (5)

1.一种电磁频谱伞罩数字域自干扰抑制方法，其特征在于：包括以下步骤：
S1.伞罩内的通信发射机产生通信信号s(t)并进行发射；
S2.伞罩干扰机根据伞罩干扰基带序列yb(n)产生自干扰信号y(t)，并向特定区域进行发射，提高该区域电磁波干扰强度，形成电磁频谱伞罩；
S3.伞罩自干扰抑制接收机对接收到的信号r(t)进行预处理，得到基带信号；
S4.伞罩自干扰抑制接收机在数字基带，对接收到的信号r(t)所对应的基带信号进行频偏估计与补偿，并将本地已知的伞罩干扰基带序列yb(n)作为本地参考干扰信号，以频偏补偿后的信号为基础，从时延、幅度、相位三个维度，对本地参考干扰信号进行估计与补偿，完成自干扰信号重建，得到本地干扰重建信号；
所述步骤S4包括以下子步骤：
S401.经过无线信道后，设接收到的干扰信号部分为ry(t)，接收到的有用信号部分为rs(t)，瞬时加性高斯白噪声为n(t)，将伞罩接收机接收到的信号对应的基带信号表示为：
Figure FDA0002463451670000011
式中，A1为幅度衰减系数，τ为相应的时延，fr为伞罩中接收干扰信号的频率，θ为相位偏移，
Figure FDA0002463451670000015
为自干扰信号y(t)经过无线传输信道后的信道衰减系数，fΔ＝fr-fc为收发两端信号的频率偏差，fc为发送信号的载波频率，rs+n(t)为接收信号中有用信号和噪声部分对应的基带信号；
S402.在数字基带，对接收到的信号r(t)进行频偏估计与补偿：
将伞罩自干扰抑制接收机接收到信号r(t)对应的复基带信号采样序列表示为r(n)：
Figure FDA0002463451670000012
其中，
Figure FDA0002463451670000013
为自干扰信号经无线传输信道后的信道衰减系数；rs+n(n)为接收信号中有用信号和噪声部分对应的复基带信号采样序列；
Figure FDA0002463451670000014
为伞罩干扰发射信号和干扰接收信号之间的时间差τ归一化分数时延，Ts为信号采样率；
对r(n)进行处理，推导出接收信号中干扰部分的频偏估计表达式为：
Figure FDA0002463451670000021
其中，m为每一帧中选取的导引序列长度，L为进行共轭相乘的两复信号项之间的间隔帧长；
根据频偏估计值进行频偏补偿，设
Figure FDA0002463451670000022
则接收信号的频偏已被完全补偿，补偿后的信号为：
Figure FDA0002463451670000023
其中，α·yb(n-D)为接收信号中频偏补偿后的干扰部分，
Figure FDA0002463451670000024
为接收信号中频偏补偿后的有用部分和噪声部分之和；
S403.对频偏补偿后的信号
Figure FDA0002463451670000025
进行符号级同步：将信号
Figure FDA0002463451670000026
与干扰发射信号的采样序列作滑动相关，通过相关峰的位置确定
Figure FDA0002463451670000027
序列中干扰部分的起始位置；
S404.同步至单个符号精度后，对
Figure FDA0002463451670000028
进行分数级时延估计，并根据时延估计值对干扰参考信号进行时延补偿：
若接收信号的频偏已被完全补偿，补偿后的信号为：
Figure FDA0002463451670000029
将频偏补偿后的接收信号
Figure FDA00024634516700000210
与本地参考干扰信号yref(n)＝yb(n)作FFT变换，从时域变换到频域，进行共轭相乘和幅度归一化处理，有：
Figure FDA00024634516700000211
其中，Yref(ω)为本地参考干扰信号yref(n)的FFT变换，
Figure FDA00024634516700000212
为接收信号
Figure FDA00024634516700000213
的FFT变换；
对
Figure FDA00024634516700000214
取相位，得到信号
Figure FDA00024634516700000215
所在频段ω的相位函数
Figure FDA00024634516700000216
该相位函数是一条相位直线方程，为：
Figure FDA00024634516700000217
其中，ωd＜ω＜ωu，其中ωd是信号的低端频率，ωu是信号的高端频率，该直线的斜率直接对应了两信号之间的时间差
Figure FDA0002463451670000031
即估计出信号yb(n)和信号
Figure FDA0002463451670000032
之间的时间差，为：
Figure FDA0002463451670000033
根据估计值，采用Farrow滤波器分数时延对接收信号的时延进行补偿，假设Farrow滤波器由M组N阶FIR滤波器构成，滤波器的频率响应为：
Figure FDA0002463451670000034
其中，
Figure FDA0002463451670000035
采用
Figure FDA0002463451670000036
的多项式来逼近，即：
Figure FDA0002463451670000037
使H(ejω)逼近
Figure FDA0002463451670000038
即设计一组滤波器系数Cn,m，使其满足：
Figure FDA0002463451670000039
其中，ω∈[ω0,ω1]表示数字频率上线性时延带宽，
Figure FDA00024634516700000310
表示分数时延参数的范围；
利用上述Farrow滤波器，对本地参考干扰信号进行时延补偿后，对应的基带信号为：
Figure FDA00024634516700000311
S405.对基带信号
Figure FDA00024634516700000312
到信号
Figure FDA00024634516700000313
之间的幅相衰减进行估计与补偿：
基于时延补偿后得到的本地参考干扰对应的基带信号为
Figure FDA00024634516700000314
则由最大似然估计，得出：
Figure FDA00024634516700000315
计算幅度衰减的估计值为：
Figure FDA00024634516700000316
计算相位衰减的估计值为：
Figure FDA00024634516700000317
对
Figure FDA00024634516700000318
进行幅相补偿，得到本地干扰重建信号为：
Figure FDA00024634516700000319
S5.伞罩自干扰抑制接收机根据时频幅相补偿得到的本地干扰重建信号进行伞罩数字域自干扰抑制；
S6.伞罩自干扰抑制接收机对抑制得到的信号进行解调处理，恢复出通信发射机发送的数据。
2.根据权利要求1所述的一种电磁频谱伞罩数字域自干扰抑制方法，其特征在于：所述步骤S2包括：
S201.将伞罩干扰基带序列yb(n)进行上变频，得到：
Figure FDA0002463451670000041
其中，fc为发送信号的载波频率，Ts为信号采样率；
S202.对信号y(n)进行数模转换，得到连续的自干扰信号y(t)；
S203.将产生自干扰信号y(t)向特定区域进行发送，提高该区域电磁波干扰强度，形成电磁频谱伞罩。
3.根据权利要求1所述的一种电磁频谱伞罩数字域自干扰抑制方法，其特征在于：所述步骤S3包括：
伞罩自干扰抑制接收机对接收到的信号r(t)依次进行模数转换、下变频和成型滤波，得到对应的基带信号。
4.一种电磁频谱伞罩数字域自干扰抑制系统，采用如权利要求1～3中任意一项所述的电磁频谱伞罩数字域自干扰抑制方法，其特征在于：包括通信发射机、伞罩干扰机和伞罩自干扰抑制接收机；
所述通信发射机，用于产生通信信号并进行发射；
所述伞罩干扰机，用于根据伞罩干扰基带序列产生自干扰信号，并向特定区域进行发射，提高该区域电磁波干扰强度，形成电磁频谱伞罩；
所述伞罩自干扰抑制接收机，用于接收电磁频谱伞罩中的信号变换到基带，在数字基带对基带信号进行频偏估计与补偿，以频偏补偿后的信号为基础，从时延、幅度、相位三个维度，对本地参考干扰信号进行估计与补偿，完成自干扰信号重建，得到本地干扰重建信号；根据本地干扰重建信号进行伞罩数字域自干扰抑制；然后对抑制得到的信号进行解调处理，恢复出通信发射机发送的数据。
5.根据权利要求4所述的一种电磁频谱伞罩数字域自干扰抑制系统，其特征在于：所述伞罩自干扰抑制接收机包括：
信号接收转换单元，用于接收电磁频谱伞罩中的信号，对接收到的信号进行预处理，得到基带信号；
自干扰重建单元，用于在数字基带对基带信号进行频偏估计与补偿，以频偏补偿后的信号为基础，从时延、幅度、相位三个维度，对本地参考干扰信号进行估计与补偿，完成自干扰信号重建，得到本地干扰重建信号；
自干扰抑制单元，用于根据时频幅相补偿后的本地干扰重建信号进行伞罩数字域自干扰抑制；
解调恢复单元，用于对抑制得到的信号进行解调处理，恢复出通信发射机发送的数据。
CN201910272810.5A 2019-04-04 2019-04-04 一种电磁频谱伞罩数字域自干扰抑制方法及系统 Active CN109962714B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910272810.5A CN109962714B (zh) 2019-04-04 2019-04-04 一种电磁频谱伞罩数字域自干扰抑制方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910272810.5A CN109962714B (zh) 2019-04-04 2019-04-04 一种电磁频谱伞罩数字域自干扰抑制方法及系统

Publications (2)

Publication Number Publication Date
CN109962714A CN109962714A (zh) 2019-07-02
CN109962714B true CN109962714B (zh) 2020-07-03

Family

ID=67025851

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910272810.5A Active CN109962714B (zh) 2019-04-04 2019-04-04 一种电磁频谱伞罩数字域自干扰抑制方法及系统

Country Status (1)

Country Link
CN (1) CN109962714B (zh)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111654458B (zh) * 2020-05-15 2021-10-08 电子科技大学 一种电磁频谱伞罩频域快速自干扰抑制方法及系统
CN112162159A (zh) * 2020-09-08 2021-01-01 中电科仪器仪表有限公司 一种基于信号识别的电磁干扰现场测试方法
US11456901B2 (en) 2020-09-25 2022-09-27 Cypress Semiconductor Corporation Transmit spur detection and mitigation for wireless communications devices
CN113708880B (zh) * 2021-07-30 2022-11-08 电子科技大学 一种基于共生干扰抑制的频谱共生实现方法
CN113691267B (zh) * 2021-09-06 2023-03-10 中国工程物理研究院电子工程研究所 一种同时双频带内载波聚合接收端自干扰抑制方法
CN113872898B (zh) * 2021-09-18 2024-03-29 中山大学 一种基于有界成分分析的数字域自干扰抑制的方法和系统
CN113765617B (zh) * 2021-09-30 2023-09-22 电子科技大学 一种基于反射放大面对抗同频干扰的方法
CN114204946A (zh) * 2022-01-27 2022-03-18 南通大学 一种收发同路的空间激光通信终端后向杂散光抑制方法
CN116112324B (zh) * 2023-02-02 2024-10-15 电子科技大学 一种同时同频自干扰收发联合抑制方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106508102B (zh) * 2010-03-05 2014-02-12 西南电子电信技术研究所 一种自干扰通信方法及系统
CN107276713A (zh) * 2017-06-27 2017-10-20 电子科技大学 一种三维区域电磁频谱伞罩系统及其组建方法
CN107359896A (zh) * 2017-06-27 2017-11-17 电子科技大学 一种三维区域内电磁波设备的干扰抑制系统及方法
CN107465481A (zh) * 2017-08-28 2017-12-12 电子科技大学 一种自干扰功率排序的电磁频谱伞罩干扰抑制系统与方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106508102B (zh) * 2010-03-05 2014-02-12 西南电子电信技术研究所 一种自干扰通信方法及系统
CN107276713A (zh) * 2017-06-27 2017-10-20 电子科技大学 一种三维区域电磁频谱伞罩系统及其组建方法
CN107359896A (zh) * 2017-06-27 2017-11-17 电子科技大学 一种三维区域内电磁波设备的干扰抑制系统及方法
CN107465481A (zh) * 2017-08-28 2017-12-12 电子科技大学 一种自干扰功率排序的电磁频谱伞罩干扰抑制系统与方法

Also Published As

Publication number Publication date
CN109962714A (zh) 2019-07-02

Similar Documents

Publication Publication Date Title
CN109962714B (zh) 一种电磁频谱伞罩数字域自干扰抑制方法及系统
DE68927658T2 (de) Hochfrequenzendgerät für Entstörungskommunikationssystem
US11424777B2 (en) Methods and systems for mitigation of interference signals for a wireless network receiver
US20080170560A1 (en) Training sequence for wireless communication system
CN111585690B (zh) 一种多制式同频干扰抵消方法
Brandes et al. Compensation of the impact of interference mitigation by pulse blanking in OFDM systems
Shaboyan et al. Frequency and timing synchronization for in-band full-duplex OFDM system
EP2135415B1 (en) Adaptive channel estimation
US6249518B1 (en) TDMA single antenna co-channel interference cancellation
Gutierrez et al. Joint radar-communications system implementation using software defined radios: Feasibility and results
EP2235967B1 (en) Common wave and sideband mitigation communication systems and methods for increasing communication speeds,spectral efficiency and enabling other benefits
CN113114407B (zh) 一种保密通信方法及系统
US7978796B2 (en) Recovering symbols in a communication receiver
CN111654458B (zh) 一种电磁频谱伞罩频域快速自干扰抑制方法及系统
Bora et al. An all-digital wideband ofdm-based frequency-hopping system using rf sampling data converters
Guan et al. Iterative phase noise suppression for full-duplex communication systems
EP2745406A1 (en) Method for controlling performance in a radio base station arranged for communication in tdd mode, and radio base station
JP3779554B2 (ja) 周波数オフセット推定方法及び周波数オフセット推定器
WO2005002067A2 (en) Multi-carrier spread spectrum using non-linear modification of sub-carrier bands
Alaghbari et al. Design and validation of the physical layer functions of FBMC/OQAM transceiver with improved residual phase error correction
CN116112324B (zh) 一种同时同频自干扰收发联合抑制方法
US8255000B2 (en) Process for computing a frequency offset for a UMTS communication system based on the CPICH pilot signals
WO2010138442A1 (en) Geometric detector for communicating through constant modulus (cm) interferers
CN102907007B (zh) 用于通过路径选择来改善干扰消除的方法和系统
Okamura et al. Channel compensation with virtual pilot signal and time-frequency interferometry for MIMO-OFDM

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