CN109962714A - 一种电磁频谱伞罩数字域自干扰抑制方法及系统 - Google Patents
一种电磁频谱伞罩数字域自干扰抑制方法及系统 Download PDFInfo
- Publication number
- CN109962714A CN109962714A CN201910272810.5A CN201910272810A CN109962714A CN 109962714 A CN109962714 A CN 109962714A CN 201910272810 A CN201910272810 A CN 201910272810A CN 109962714 A CN109962714 A CN 109962714A
- 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.)
- Granted
Links
- 238000001228 spectrum Methods 0.000 title claims abstract description 50
- 238000000034 method Methods 0.000 title claims abstract description 26
- 230000001629 suppression Effects 0.000 claims abstract description 85
- 238000004891 communication Methods 0.000 claims abstract description 48
- 230000005540 biological transmission Effects 0.000 claims abstract description 9
- 238000006243 chemical reaction Methods 0.000 claims description 9
- 230000002452 interceptive effect Effects 0.000 claims description 9
- 238000005070 sampling Methods 0.000 claims description 8
- 230000009466 transformation Effects 0.000 claims description 7
- 238000007781 pre-processing Methods 0.000 claims description 6
- 239000000654 additive Substances 0.000 claims description 5
- 230000000996 additive effect Effects 0.000 claims description 5
- 238000010606 normalization Methods 0.000 claims description 5
- 238000012545 processing Methods 0.000 claims description 5
- 238000007476 Maximum Likelihood Methods 0.000 claims description 3
- 238000001914 filtration Methods 0.000 claims description 3
- 238000011084 recovery Methods 0.000 claims description 3
- 230000004044 response Effects 0.000 claims description 3
- 230000005764 inhibitory process Effects 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000001934 delay Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011017 operating method Methods 0.000 description 1
- 230000000452 restraining effect Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
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)进行上变频，得到：
其中，fc为发送信号的载波频率，Ts为信号采样率；
S202.对信号y(n)进行数模转换，得到连续的自干扰信号y(t)；
S203.将产生自干扰信号y(t)向特定区域进行发送，提高该区域电磁波干扰强度，形成电磁频谱伞罩。
进一步地，所述步骤S3包括：
伞罩自干扰抑制接收机对接收到的信号r(t)依次进行模数转换、下变频和成型滤波，得到对应的基带信号。
进一步地，所述步骤S4包括以下子步骤：
S401.经过无线信道后，设接收到的干扰信号部分为ry(t)，接收到的有用信号部分为rs(t)，瞬时加性高斯白噪声为n(t)，将伞罩接收机接收到的信号对应的基带信号表示为：
式中，A1为幅度衰减系数，τ为相应的时延，fr为伞罩中接收干扰信号的频率，θ为相位偏移，为自干扰信号y(t)经过无线传输信道后的信道衰减系数，fΔ＝fr-fc为收发两端信号的频率偏差，rs+n(t)为接收信号中有用信号和噪声部分对应的基带信号；
S402.在数字基带，对接收到的信号r(t)进行频偏估计与补偿：
将伞罩自干扰抑制接收机接收到信号r(t)对应的复基带信号采样序列表示为r(n)：
其中，为自干扰信号经无线传输信道后的信道衰减系数；rs+n(n)为接收信号中有用信号和噪声部分对应的复基带信号采样序列；为伞罩干扰发射信号和干扰接收信号之间的时间差τ归一化分数时延；
对r(n)进行处理，推导出接收信号中干扰部分的频偏估计表达式为：
其中，m为每一帧中选取的导引序列长度，L为进行共轭相乘的两复信号项之间的间隔帧长；
根据频偏估计值进行频偏补偿，设则接收信号的频偏已被完全补偿，补偿后的信号为：
其中，α·yb(n-D)为接收信号中频偏补偿后的干扰部分，为接收信号中频偏补偿后的有用部分和噪声部分之和；
S403.对频偏补偿后的信号进行符号级同步：将信号与干扰发射信号的采样序列作滑动相关，通过相关峰的位置确定序列中干扰部分的起始位置；
S404.同步至单个符号精度后，对进行分数级时延估计，并根据时延估计值对干扰参考信号进行时延补偿：
若接收信号的频偏已被完全补偿，补偿后的信号为：
将频偏补偿后的接收信号与本地参考干扰信号yref(n)＝yb(n)作FFT变换，从时域变换到频域，进行共轭相乘和幅度归一化处理，有：
其中，Yref(ω)为本地参考干扰信号yref(n)的FFT变换，为接收信号的FFT变换；
对取相位，得到信号所在频段ω的相位函数该相位函数是一条相位直线方程，为：
其中，ωd＜ω＜ωu，其中ωd是信号的低端频率，ωu是信号的高端频率，该直线的斜率直接对应了两信号之间的时间差即估计出信号yb(n)和信号之间的时间差，为：
根据估计值，采用Farrow滤波器分数时延对接收信号的时延进行补偿，假设滤波器的频率响应为：
其中，采用的多项式来逼近，即：
使H(ejω)逼近即设计一组滤波器系数Cn,m，使其满足：
其中，ω∈[ω0,ω1]表示数字频率上线性时延带宽，表示分数时延参数的范围；
利用上述Farrow滤波器，对本地参考干扰信号进行时延补偿后，对应的基带信号为：
S405.对基带信号到信号之间的幅相衰减进行估计与补偿：
基于时延补偿后得到的本地参考干扰对应的基带信号为则由最大似然估计，得出：
计算幅度衰减的估计值为：
计算相位衰减的估计值为：
对进行幅相补偿，得到本地干扰重建信号为：
一种电磁频谱伞罩数字域自干扰抑制系统，包括通信发射机、伞罩干扰机和伞罩自干扰抑制接收机；
所述通信发射机，用于产生通信信号并进行发射；
所述伞罩干扰机，用于根据伞罩干扰基带序列产生自干扰信号，并向特定区域进行发射，提高该区域电磁波干扰强度，形成电磁频谱伞罩；
所述伞罩自干扰抑制接收机，用于接收电磁频谱伞罩中的信号变换到基带，在数字基带对基带信号进行频偏估计与补偿，以频偏补偿后的信号为基础，从时延、幅度、相位三个维度，对本地参考干扰信号进行估计与补偿，完成自干扰信号重建，得到本地干扰重建信号；根据本地干扰重建信号进行伞罩数字域自干扰抑制；然后对抑制得到的信号进行解调处理，恢复出通信发射机发送的数据。
其中，所述伞罩自干扰抑制接收机包括：信号接收转换单元，用于接收电磁频谱伞罩中的信号，对接收到的信号进行预处理，得到基带信号；自干扰重建单元，用于在数字基带对基带信号进行频偏估计与补偿，以频偏补偿后的信号为基础，从时延、幅度、相位三个维度，对本地参考干扰信号进行估计与补偿，完成自干扰信号重建，得到本地干扰重建信号；自干扰抑制单元，用于根据时频幅相补偿后的本地干扰重建信号进行伞罩数字域自干扰抑制；解调恢复单元，用于对抑制得到的信号进行解调处理，恢复出通信发射机发送的数据。
本发明的有益效果是：本发明实现了对伞罩自干扰信号频偏、时延、幅度、相位四个参数的高精度跟踪与补偿，有效解决了干扰信号的时频幅相四维对齐的问题，伞罩自干扰的抑制不受通信信号的类型制约，可以与任何已有的通信链路相结合，兼容性较强。
附图说明
图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)进行上变频，得到：
其中，fc为发送信号的载波频率，Ts为信号采样率；
S202.对信号y(n)进行数模转换，得到连续的自干扰信号y(t)；
S203.将产生自干扰信号y(t)向特定区域进行发送，提高该区域电磁波干扰强度，形成电磁频谱伞罩。
进一步地，所述步骤S3包括：
伞罩自干扰抑制接收机对接收到的信号r(t)依次进行模数转换、下变频和成型滤波，得到对应的基带信号。
进一步地，所述步骤S4包括以下子步骤：
S401.经过无线信道后，设接收到的干扰信号部分为ry(t)，接收到的有用信号部分为rs(t)，瞬时加性高斯白噪声为n(t)，将伞罩接收机接收到的信号对应的基带信号表示为：
式中，A1为幅度衰减系数，τ为相应的时延，fr为伞罩中接收干扰信号的频率，θ为相位偏移，为自干扰信号y(t)经过无线传输信道后的信道衰减系数，fΔ＝fr-fc为收发两端信号的频率偏差，rs+n(t)为接收信号中有用信号和噪声部分对应的基带信号；
S402.在数字基带，对接收到的信号r(t)进行频偏估计与补偿：
将伞罩自干扰抑制接收机接收到信号r(t)对应的复基带信号采样序列表示为r(n)：
其中，为自干扰信号经无线传输信道后的信道衰减系数；rs+n(n)为接收信号中有用信号和噪声部分对应的复基带信号采样序列；为伞罩干扰发射信号和干扰接收信号之间的时间差τ归一化分数时延；
对r(n)进行处理，推导出接收信号中干扰部分的频偏估计表达式为：
其中，m为每一帧中选取的导引序列长度，L为进行共轭相乘的两复信号项之间的间隔帧长；
根据频偏估计值进行频偏补偿，设则接收信号的频偏已被完全补偿，补偿后的信号为：
其中，α·yb(n-D)为接收信号中频偏补偿后的干扰部分，为接收信号中频偏补偿后的有用部分和噪声部分之和；
S403.对频偏补偿后的信号进行符号级同步：将信号与干扰发射信号的采样序列作滑动相关，通过相关峰的位置确定序列中干扰部分的起始位置；
S404.同步至单个符号精度后，对进行分数级时延估计，并根据时延估计值对干扰参考信号进行时延补偿：
若接收信号的频偏已被完全补偿，补偿后的信号为：
将频偏补偿后的接收信号与本地参考干扰信号yref(n)＝yb(n)作FFT变换，从时域变换到频域，进行共轭相乘和幅度归一化处理，有：
其中，Yref(ω)为本地参考干扰信号yref(n)的FFT变换，为接收信号的FFT变换；
对取相位，得到信号所在频段ω的相位函数该相位函数是一条相位直线方程，为：
其中，ωd＜ω＜ωu，其中ωd是信号的低端频率，ωu是信号的高端频率，该直线的斜率直接对应了两信号之间的时间差，即估计出信号yb(n)和信号之间的时间差，为：
根据估计值，采用Farrow滤波器分数时延对接收信号的时延进行补偿，假设滤波器的频率响应为：
其中，采用的多项式来逼近，即：
使H(ejω)逼近即设计一组滤波器系数Cn,m，使其满足：
其中，ω∈[ω0,ω1]表示数字频率上线性时延带宽，表示分数时延参数的范围；
利用上述Farrow滤波器，对本地参考干扰信号进行时延补偿后，对应的基带信号为：
S405.对基带信号到信号之间的幅相衰减进行估计与补偿：
基于时延补偿后得到的本地参考干扰对应的基带信号为则由最大似然估计，得出：
计算幅度衰减的估计值为：
计算相位衰减的估计值为：
对进行幅相补偿，得到本地干扰重建信号为：
一种电磁频谱伞罩数字域自干扰抑制系统，包括通信发射机、伞罩干扰机和伞罩自干扰抑制接收机；
所述通信发射机，用于产生通信信号并进行发射；
所述伞罩干扰机，用于根据伞罩干扰基带序列产生自干扰信号，并向特定区域进行发射，提高该区域电磁波干扰强度，形成电磁频谱伞罩；
所述伞罩自干扰抑制接收机，用于接收电磁频谱伞罩中的信号变换到基带，在数字基带对基带信号进行频偏估计与补偿，以频偏补偿后的信号为基础，从时延、幅度、相位三个维度，对本地参考干扰信号进行估计与补偿，完成自干扰信号重建，得到本地干扰重建信号；根据本地干扰重建信号进行伞罩数字域自干扰抑制；然后对抑制得到的信号进行解调处理，恢复出通信发射机发送的数据。
其中，所述伞罩自干扰抑制接收机包括：信号接收转换单元，用于接收电磁频谱伞罩中的信号，对接收到的信号进行预处理，得到基带信号；自干扰重建单元，用于在数字基带对基带信号进行频偏估计与补偿，以频偏补偿后的信号为基础，从时延、幅度、相位三个维度，对本地参考干扰信号进行估计与补偿，完成自干扰信号重建，得到本地干扰重建信号；自干扰抑制单元，用于根据时频幅相补偿后的本地干扰重建信号进行伞罩数字域自干扰抑制；解调恢复单元，用于对抑制得到的信号进行解调处理，恢复出通信发射机发送的数据。
在本申请的实施例中，以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。
为了对本发明提出的一种电磁频谱伞罩数字域自干扰抑制方法的性能进行验证，对电磁频谱伞罩实验设备参数设置如下：
图8给出了通信发射机保持静默时，伞罩自干扰抑制前和抑制后的我方接收机信号功率谱。当我方接收机的干噪比为30dB时，自干扰抑制度为28.5dB。自干扰抑制度基本不受自干扰发射功率的影响。
图9给出了伞罩内的通信发射机正常发射FM调制语音信号时，自干扰抑制前和抑制后的我方接收机信号功率谱。我方天线输出口的干信比为10dB。自干扰被抑制前，期望的FM调制信号被淹没在干扰信号频谱中，无法正常解调和播放语音信号；自干扰抑制后，能够看到期望的FM调制信号频谱轮廓，可以正常解调和播放语音信号。FM收音机一直无法正常解调和播放语音信号。
综上所述，考虑无线通信系统中伞罩数字域自干扰抑制的传输场景，针对加性白高斯噪声无线传播信道，在加性白高斯噪声无线自干扰信道中，本文提出的一种电磁频谱伞罩数字域自干扰抑制方法在指定区域可以干扰非法通信，同时根据干扰特征在本地重建并抵消干扰，以保证己方正常通信的集干扰与通信为一体。从仿真结果中容易看出：自干扰抑制前，通信发射机与我方接收机的通信被阻塞；自干扰抑制后，通信发射机与我方接收机可以正常进行通信。非法接收机一直不能够正常接收通信信号。该结果初步验证了本方法的可行性。
这里已经通过具体的实施例子对本发明进行了详细描述，提供上述实施例的描述为了使本领域的技术人员制造或适用本发明，这些实施例的各种修改对于本领域的技术人员来说是容易理解的。本发明并不限于这些例子，或其中的某些方面。本发明的范围通过附加的权利要求进行详细说明。
上述说明示出并描述了本发明的一个优选实施例，但如前所述，应当理解本发明并非局限于本文所披露的形式，不应看作是对其他实施例的排除，而可用于各种其他组合、修改和环境，并能够在本文所述发明构想范围内，通过上述教导或相关领域的技术或知识进行改动。而本领域人员所进行的改动和变化不脱离本发明的精神和范围，则都应在本发明所附权利要求的保护范围内。
Claims (6)
1.一种电磁频谱伞罩数字域自干扰抑制方法，其特征在于：包括以下步骤：
S1.伞罩内的通信发射机产生通信信号s(t)并进行发射；
S2.伞罩干扰机根据伞罩干扰基带序列yb(n)产生自干扰信号y(t)，并向特定区域进行发射，提高该区域电磁波干扰强度，形成电磁频谱伞罩；
S3.伞罩自干扰抑制接收机对接收到的信号r(t)进行预处理，得到基带信号；
S4.伞罩自干扰抑制接收机在数字基带，对接收到的信号r(t)所对应的基带信号进行频偏估计与补偿，并将本地已知的伞罩干扰基带序列yb(n)作为本地参考干扰信号，以频偏补偿后的信号为基础，从时延、幅度、相位三个维度，对本地参考干扰信号进行估计与补偿，完成自干扰信号重建，得到本地干扰重建信号；
S5.伞罩自干扰抑制接收机根据时频幅相补偿得到的本地干扰重建信号进行伞罩数字域自干扰抑制；
S6.伞罩自干扰抑制接收机对抑制得到的信号进行解调处理，恢复出通信发射机发送的数据。
2.根据权利要求1所述的一种电磁频谱伞罩数字域自干扰抑制方法，其特征在于：所述步骤S2包括：
S201.将伞罩干扰基带序列yb(n)进行上变频，得到：
其中，fc为发送信号的载波频率，Ts为信号采样率；
S202.对信号y(n)进行数模转换，得到连续的自干扰信号y(t)；
S203.将产生自干扰信号y(t)向特定区域进行发送，提高该区域电磁波干扰强度，形成电磁频谱伞罩。
3.根据权利要求1所述的一种电磁频谱伞罩数字域自干扰抑制方法，其特征在于：所述步骤S3包括：
伞罩自干扰抑制接收机对接收到的信号r(t)依次进行模数转换、下变频和成型滤波，得到对应的基带信号。
4.根据权利要求1所述的一种电磁频谱伞罩数字域自干扰抑制方法，其特征在于：所述步骤S4包括以下子步骤：
S401.经过无线信道后，设接收到的干扰信号部分为ry(t)，接收到的有用信号部分为rs(t)，瞬时加性高斯白噪声为n(t)，将伞罩接收机接收到的信号对应的基带信号表示为：
式中，A1为幅度衰减系数，τ为相应的时延，fr为伞罩中接收干扰信号的频率，θ为相位偏移，为自干扰信号y(t)经过无线传输信道后的信道衰减系数，fΔ＝fr-fc为收发两端信号的频率偏差，rs+n(t)为接收信号中有用信号和噪声部分对应的基带信号；
S402.在数字基带，对接收到的信号r(t)进行频偏估计与补偿：
将伞罩自干扰抑制接收机接收到信号r(t)对应的复基带信号采样序列表示为r(n)：
其中，为自干扰信号经无线传输信道后的信道衰减系数；rs+n(n)为接收信号中有用信号和噪声部分对应的复基带信号采样序列；为伞罩干扰发射信号和干扰接收信号之间的时间差τ归一化分数时延；
对r(n)进行处理，推导出接收信号中干扰部分的频偏估计表达式为：
其中，m为每一帧中选取的导引序列长度，L为进行共轭相乘的两复信号项之间的间隔帧长；
根据频偏估计值进行频偏补偿，设则接收信号的频偏已被完全补偿，补偿后的信号为：
其中，α·yb(n-D)为接收信号中频偏补偿后的干扰部分，为接收信号中频偏补偿后的有用部分和噪声部分之和；
S403.对频偏补偿后的信号进行符号级同步：将信号与干扰发射信号的采样序列作滑动相关，通过相关峰的位置确定序列中干扰部分的起始位置；
S404.同步至单个符号精度后，对进行分数级时延估计，并根据时延估计值对干扰参考信号进行时延补偿：
若接收信号的频偏已被完全补偿，补偿后的信号为：
将频偏补偿后的接收信号与本地参考干扰信号yref(n)＝yb(n)作FFT变换，从时域变换到频域，进行共轭相乘和幅度归一化处理，有：
其中，Yref(ω)为本地参考干扰信号yref(n)的FFT变换，为接收信号的FFT变换；
对取相位，得到信号所在频段ω的相位函数该相位函数是一条相位直线方程，为：
其中，ωd＜ω＜ωu，其中ωd是信号的低端频率，ωu是信号的高端频率，该直线的斜率直接对应了两信号之间的时间差即估计出信号yb(n)和信号之间的时间差，为：
根据估计值，采用Farrow滤波器分数时延对接收信号的时延进行补偿，假设滤波器的频率响应为：
其中，采用的多项式来逼近，即：
使H(ejω)逼近即设计一组滤波器系数Cn,m，使其满足：
其中，ω∈[ω0,ω1]表示数字频率上线性时延带宽，表示分数时延参数的范围；
利用上述Farrow滤波器，对本地参考干扰信号进行时延补偿后，对应的基带信号为：
S405.对基带信号到信号之间的幅相衰减进行估计与补偿：
基于时延补偿后得到的本地参考干扰对应的基带信号为则由最大似然估计，得出：
计算幅度衰减的估计值为：
计算相位衰减的估计值为：
对进行幅相补偿，得到本地干扰重建信号为：
5.一种电磁频谱伞罩数字域自干扰抑制系统，其特征在于：包括通信发射机、伞罩干扰机和伞罩自干扰抑制接收机；
所述通信发射机，用于产生通信信号并进行发射；
所述伞罩干扰机，用于根据伞罩干扰基带序列产生自干扰信号，并向特定区域进行发射，提高该区域电磁波干扰强度，形成电磁频谱伞罩；
所述伞罩自干扰抑制接收机，用于接收电磁频谱伞罩中的信号变换到基带，在数字基带对基带信号进行频偏估计与补偿，以频偏补偿后的信号为基础，从时延、幅度、相位三个维度，对本地参考干扰信号进行估计与补偿，完成自干扰信号重建，得到本地干扰重建信号；根据本地干扰重建信号进行伞罩数字域自干扰抑制；然后对抑制得到的信号进行解调处理，恢复出通信发射机发送的数据。
6.一种电磁频谱伞罩数字域自干扰抑制系统，其特征在于：所述伞罩自干扰抑制接收机包括：
信号接收转换单元，用于接收电磁频谱伞罩中的信号，对接收到的信号进行预处理，得到基带信号；
自干扰重建单元，用于在数字基带对基带信号进行频偏估计与补偿，以频偏补偿后的信号为基础，从时延、幅度、相位三个维度，对本地参考干扰信号进行估计与补偿，完成自干扰信号重建，得到本地干扰重建信号；
自干扰抑制单元，用于根据时频幅相补偿后的本地干扰重建信号进行伞罩数字域自干扰抑制；
解调恢复单元，用于对抑制得到的信号进行解调处理，恢复出通信发射机发送的数据。
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 true CN109962714A (zh) | 2019-07-02 |
| CN109962714B 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) |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111654458A (zh) * | 2020-05-15 | 2020-09-11 | 电子科技大学 | 一种电磁频谱伞罩频域快速自干扰抑制方法及系统 |
| CN112162159A (zh) * | 2020-09-08 | 2021-01-01 | 中电科仪器仪表有限公司 | 一种基于信号识别的电磁干扰现场测试方法 |
| CN113691267A (zh) * | 2021-09-06 | 2021-11-23 | 中国工程物理研究院电子工程研究所 | 一种同时双频带内载波聚合接收端自干扰抑制方法 |
| CN113708880A (zh) * | 2021-07-30 | 2021-11-26 | 电子科技大学 | 一种基于共生干扰抑制的频谱共生实现方法 |
| CN113765617A (zh) * | 2021-09-30 | 2021-12-07 | 电子科技大学 | 一种基于反射放大面对抗同频干扰的方法 |
| CN113872898A (zh) * | 2021-09-18 | 2021-12-31 | 中山大学 | 一种基于有界成分分析的数字域自干扰抑制的方法和系统 |
| CN114204946A (zh) * | 2022-01-27 | 2022-03-18 | 南通大学 | 一种收发同路的空间激光通信终端后向杂散光抑制方法 |
| CN114257255A (zh) * | 2020-09-25 | 2022-03-29 | 赛普拉斯半导体公司 | 针对无线通信设备的发送杂散检测和抑制 |
| CN116112324A (zh) * | 2023-02-02 | 2023-05-12 | 电子科技大学 | 一种同时同频自干扰收发联合抑制方法 |
Citations (4)
| 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 | 电子科技大学 | 一种自干扰功率排序的电磁频谱伞罩干扰抑制系统与方法 |
-
2019
- 2019-04-04 CN CN201910272810.5A patent/CN109962714B/zh active Active
Patent Citations (4)
| 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 | 电子科技大学 | 一种自干扰功率排序的电磁频谱伞罩干扰抑制系统与方法 |
Cited By (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111654458B (zh) * | 2020-05-15 | 2021-10-08 | 电子科技大学 | 一种电磁频谱伞罩频域快速自干扰抑制方法及系统 |
| CN111654458A (zh) * | 2020-05-15 | 2020-09-11 | 电子科技大学 | 一种电磁频谱伞罩频域快速自干扰抑制方法及系统 |
| CN112162159A (zh) * | 2020-09-08 | 2021-01-01 | 中电科仪器仪表有限公司 | 一种基于信号识别的电磁干扰现场测试方法 |
| CN114257255A (zh) * | 2020-09-25 | 2022-03-29 | 赛普拉斯半导体公司 | 针对无线通信设备的发送杂散检测和抑制 |
| US11804993B2 (en) | 2020-09-25 | 2023-10-31 | Cypress Semiconductor Corporation | Transmit spur detection and mitigation for wireless communications devices |
| CN114257255B (zh) * | 2020-09-25 | 2023-10-27 | 赛普拉斯半导体公司 | 针对无线通信设备的发送杂散检测和抑制 |
| CN113708880A (zh) * | 2021-07-30 | 2021-11-26 | 电子科技大学 | 一种基于共生干扰抑制的频谱共生实现方法 |
| CN113691267A (zh) * | 2021-09-06 | 2021-11-23 | 中国工程物理研究院电子工程研究所 | 一种同时双频带内载波聚合接收端自干扰抑制方法 |
| CN113872898A (zh) * | 2021-09-18 | 2021-12-31 | 中山大学 | 一种基于有界成分分析的数字域自干扰抑制的方法和系统 |
| CN113872898B (zh) * | 2021-09-18 | 2024-03-29 | 中山大学 | 一种基于有界成分分析的数字域自干扰抑制的方法和系统 |
| CN113765617B (zh) * | 2021-09-30 | 2023-09-22 | 电子科技大学 | 一种基于反射放大面对抗同频干扰的方法 |
| CN113765617A (zh) * | 2021-09-30 | 2021-12-07 | 电子科技大学 | 一种基于反射放大面对抗同频干扰的方法 |
| CN114204946A (zh) * | 2022-01-27 | 2022-03-18 | 南通大学 | 一种收发同路的空间激光通信终端后向杂散光抑制方法 |
| CN116112324A (zh) * | 2023-02-02 | 2023-05-12 | 电子科技大学 | 一种同时同频自干扰收发联合抑制方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109962714B (zh) | 2020-07-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109962714B (zh) | 一种电磁频谱伞罩数字域自干扰抑制方法及系统 | |
| DE68927658T2 (de) | Hochfrequenzendgerät für Entstörungskommunikationssystem | |
| US8483338B2 (en) | Channel estimation method and system for inter-carrier interference-limited wireless communication network | |
| US7392018B1 (en) | Channel estimation applique for wireless communications | |
| JP4343694B2 (ja) | 送信チャネル間のゲインオフセットを決定する方法 | |
| US20110105037A1 (en) | Methods and systems for interference cancellation in multi-mode coexistence modems | |
| US11424777B2 (en) | Methods and systems for mitigation of interference signals for a wireless network receiver | |
| Brandes et al. | Compensation of the impact of interference mitigation by pulse blanking in OFDM systems | |
| CN105812083A (zh) | 一种同时同频全双工系统中的射频快速自干扰抵消方法 | |
| CN104617981A (zh) | 一种基于谱估计的全双工自干扰消除方法及装置 | |
| Shaboyan et al. | Frequency and timing synchronization for in-band full-duplex OFDM system | |
| EP2135415B1 (en) | Adaptive channel estimation | |
| Gutierrez et al. | Joint radar-communications system implementation using software defined radios: Feasibility and results | |
| US6249518B1 (en) | TDMA single antenna co-channel interference cancellation | |
| KR100979151B1 (ko) | 복소 가중치 생성(ｃｗｇ) 알고리즘을 이용한 레이크 수신기에서의 채널 이득 추정 | |
| CN113114407B (zh) | 一种保密通信方法及系统 | |
| US7978796B2 (en) | Recovering symbols in a communication receiver | |
| CN111654458B (zh) | 一种电磁频谱伞罩频域快速自干扰抑制方法及系统 | |
| Nabeel et al. | Low-Complexity Linear and Non-Linear Digital Self-Interference Cancellation | |
| Bora et al. | An all-digital wideband ofdm-based frequency-hopping system using rf sampling data converters | |
| JP3779554B2 (ja) | 周波数オフセット推定方法及び周波数オフセット推定器 | |
| US20070064773A1 (en) | Multi-carrier spread spectrum using non-linear modification of sub-carrier bands | |
| Guan et al. | Iterative phase noise suppression for full-duplex communication systems | |
| US7254181B2 (en) | Transmitting apparatus and method, receiving apparatus and method, and transmitting and receiving apparatus and method | |
| Alaghbari et al. | Design and validation of the physical layer functions of FBMC/OQAM transceiver with improved residual phase error correction |
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 |