CN116073853A - 一种射频自干扰消除系统及方法 - Google Patents

一种射频自干扰消除系统及方法 Download PDF

Info

Publication number
CN116073853A
CN116073853A CN202211728329.0A CN202211728329A CN116073853A CN 116073853 A CN116073853 A CN 116073853A CN 202211728329 A CN202211728329 A CN 202211728329A CN 116073853 A CN116073853 A CN 116073853A
Authority
CN
China
Prior art keywords
amplitude
phase control
signal
radio frequency
residual 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.)
Pending
Application number
CN202211728329.0A
Other languages
English (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.)
Beijing Jiaotong University
Original Assignee
Beijing Jiaotong 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 Beijing Jiaotong University filed Critical Beijing Jiaotong University
Priority to CN202211728329.0A priority Critical patent/CN116073853A/zh
Publication of CN116073853A publication Critical patent/CN116073853A/zh
Pending legal-status Critical Current

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/38—Transceivers, i.e. devices in which transmitter and receiver form a structural unit and in which at least one part is used for functions of transmitting and receiving
    • H04B1/40—Circuits
    • H04B1/50—Circuits using different frequencies for the two directions of communication
    • H04B1/52—Hybrid arrangements, i.e. arrangements for transition from single-path two-direction transmission to single-direction transmission on each of two paths or vice versa
    • H04B1/525—Hybrid arrangements, i.e. arrangements for transition from single-path two-direction transmission to single-direction transmission on each of two paths or vice versa with means for reducing leakage of transmitter signal into the receiver
    • 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/38—Transceivers, i.e. devices in which transmitter and receiver form a structural unit and in which at least one part is used for functions of transmitting and receiving
    • H04B1/40—Circuits
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D30/00—Reducing energy consumption in communication networks
    • Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks

Abstract

本发明涉及一种射频自干扰消除系统及方法。利用初始随机产生的若干幅相控制信号对接收信号进行射频自干扰消除，得到若干组残余信号及其功率；选取当前残余信号功率较大的部分幅相控制信号执行莱维随机游走及二元化处理，产生修正幅相控制信号；该组修正的幅相控制信号分别产生修正残余信号；比较当前残余信号功率与修正残余信号功率，保留每组比较结果的较小值，用于更新残余信号功率和幅相控制信号；将更新结果与当前残余信号功率最小值和对应的幅相控制信号分别构成下一轮迭代的残余信号功率及幅相控制信号，算法迭代；当算法满足停止准则时，系统选取残余信号功率最小值对应的幅相位控制信号作为最优幅相控制信号，完成系统射频自干扰消除。

Description

一种射频自干扰消除系统及方法
技术领域
本发明涉及无线通信技术领域，特别是涉及一种射频自干扰消除系统及方法。
背景技术
相较于时分双工与频分双工技术，同时同频全双工技术理论上可在同一时隙、同一频段实现数据的双向传输，可实现频谱利用率的倍增，具有广阔的应用前景。然而，同时同频全双工技术面临强自干扰问题。全双工设备接收信号时会受到自身发射信号的干扰，由于通信节点的发射天线与接收天线之间的距离很近，自干扰信号强度远大于有用信号强度，接收机难以检测到有用信号，严重影响接收性能。因此，自干扰消除是同时同频全双工系统的关键技术之一
现有的自干扰消除方案包括被动消除和数字消除、射频消除三种方式。被动消除技术通常依赖于硬件设施，例如采用环形器进行收发隔离、增加全双工设备的收发天线间距等，该方式硬件成本较高。数字消除技术主要通过信道估计重建自干扰信号，通常可实现的干扰抑制性能较低。射频消除技术通过射频电路调整信号幅度、相位，重构出与自干扰信号等幅、反向的信号，实现干扰抵消，该方案研究广泛、自干扰消除效果显著。
发明内容
本发明的目的在于提供一种射频自干扰消除系统及方法，通过算法处理单元及幅相控制单元重构出与自干扰信号等幅反向的对消信号，生成的对消信号与含有用信号的混合信号合路相加，从而实现自干扰消除。对消信号重构过程中，引入莱维随机游走全局搜索最优的幅相控制信号，提高自干扰对消信号的精度。
为达到以上目的，本发明提供了一种射频自干扰消除系统，包括发送信号耦合单元、射频隔离单元、信号合成单元、接收信号耦合单元、功率检测单元、算法处理单元、数模转换单元及幅相控制单元；
所述发送信号耦合单元用于将发送端射频信号按比例分成两路，分别输出至所述射频隔离单元及幅相控制单元；
所述射频隔离单元用于隔离发送端和接收端之间的射频信号；
所述信号合成单元用于将两路输入信号合并为一路输出信号；
所述接收信号耦合单元用于将所述信号合成单元输出的射频信号按比例分成两路，分别输出至接收端及所述功率测量单元；
所述功率测量单元用于测量并获取所输入的合并信号功率值；
所述算法处理单元用于执行所设计的射频自干扰消除方法；
所述数模转换单元用于将所述算法处理单元输出的数字幅相控制信号转换为模拟幅相控制信号；
所述幅相控制单元用于接收所述数模转换单元输出的模拟幅相控制信号，并调整信号幅度和相位，输出射频对消信号至所述信号合成单元。
在上述方案的基础上，所述发送信号耦合单元将发送端射频信号分为两路，其中一路输出至射频隔离单元第一端口，经射频隔离单元第二端口输出至天线，另一路输出至所述幅相控制单元；
射频隔离单元第三端口的输出信号为射频隔离单元第二端口输入的接收射频信号及由射频隔离单元第一端口泄露至射频隔离单元第三端口的射频自干扰信号的叠加；
所述信号合成单元将所述幅相控制单元输出的射频对消信号与射频隔离单元第三端口输出的信号合并为一路残余信号，输出至所述接收信号耦合单元；
所述接收信号耦合单元将信号合成单元的输出信号分为两路，其中一路输出至接收端，另一路输出至所述功率测量单元；
所述功率测量单元测量接收信号耦合单元输出的残余信号功率，并将功率值输出至所述算法处理单元；
所述算法处理单元执行射频自干扰消除算法，产生数字幅相控制信号，并输出至所述数模转换单元；
所述数模转换单元将算法处理单元输出的数字幅相控制信号转换为模拟幅相控制信号，并输出至所述幅相控制单元；
所述幅相控制单元根据模拟幅相控制信号，调节信号幅度和相位并输出射频对消信号至所述信号合成单元。
基于上述方案的基础上，提供一种射频自干扰消除方法，具体包括如下步骤：
步骤一：初始化标志位F＝1，迭代计数器t＝0，迭代计数器溢出上限tmax，N个幅相控制信号以及莱维随机游走参数集合；
步骤二：N个幅相控制信号构成当前幅相控制信号集合；依次利用N个幅相控制信号产生N个残余信号，得到当前残余信号功率集合；
步骤三：判断标志位F是否为1，若是，转至步骤四；否则，转至步骤九；
步骤四：从残余信号功率集合中选取最小值，记录其幅相控制信号；剩余N-1个残余信号功率及其对应的幅相控制信号分别构成集合；
步骤五：对剩余N-1个幅相控制信号分别执行莱维随机游走，产生N-1个修正幅相控制信号，采用二元化处理后构成集合；
步骤六：依次利用N-1个修正幅相控制信号产生N-1个修正残余信号，得到修正残余信号功率集合；
步骤七：依序比较N-1组当前残余信号功率与修正残余信号功率，每组比较结果取二者较小值，构成残余信号功率集合及对应的幅相控制信号集合；
步骤八：将比较所得残余信号功率集合与步骤四中残余信号功率最小值构成下一轮迭代的当前残余信号功率集合，集合中残余信号功率对应的幅相控制信号构成下一轮迭代的当前幅相控制信号集合，置标志位F＝2，t＝t+1，并转至步骤三；
步骤九：判断算法是否满足停止准则，若是，转至步骤十；否则，转至步骤十一；
步骤十：从当前残余信号功率集合中选取功率最小值，将其对应的幅相位控制信号作为最优幅相控制信号，并进行输出；
步骤十一：置标志位F＝1，转至步骤三。
在上述方案的基础上，步骤一中所述幅相控制信号由长度为2M的二进制数组xn，1≤n≤N，表示，其中，M为数模转换单元采样位数，二进制数组xn的前M位表示幅度控制信号，后M位表示相位控制信号；
所述莱维随机游走参数集由二元组<α,β>表示，其中，α为步长缩放因子，β为莱维指数。
在上述方案的基础上，步骤三中所述当前幅相控制信号集合表示为
Figure BDA0004028157890000051
当前残余信号功率集合表示为
Figure BDA0004028157890000052
Figure BDA0004028157890000053
其中，Pn为第n个残余信号的功率。
在上述方案的基础上，步骤四中所述从残余信号功率集合中选取最小值
Figure BDA0004028157890000054
记录其幅相控制信号
Figure BDA0004028157890000055
其中，n*为残余信号功率最小值下标，由公式(1)计算：
Figure BDA0004028157890000056
其中，argmin(·)(·)为取集合中最小元素下标运算符，
Figure BDA0004028157890000057
为当前残余信号功率集合；
所述剩余N-1个残余信号功率集合可表示为
Figure BDA0004028157890000058
\为集合差集运算符；其对应的幅相控制信号集合可表示为
Figure BDA0004028157890000059
Figure BDA00040281578900000510
在上述方案的基础上，步骤五中所述对剩余N-1个幅相控制信号分别执行莱维随机游走，产生N-1个修正幅相控制信号，可由公式(2)表示：
Figure BDA0004028157890000061
其中，xi为集合
Figure BDA0004028157890000062
中的幅相控制信号，
Figure BDA0004028157890000063
为修正幅相控制信号，α为步长缩放因子，
Figure BDA0004028157890000064
为修正幅相控制信号的随机步长向量，si服从莱维分布，由公式(3)产生：
Figure BDA0004028157890000065
其中，⊙为哈达玛积运算符，|·|为绝对值运算符，
Figure BDA0004028157890000066
为均值为零向量协方差矩阵为σ2I2M的高斯随机向量，
Figure BDA0004028157890000067
IQ为Q阶单位阵，β为莱维指数，Γ(·)为伽马函数，sin(·)为正弦函数，π为圆周率，
Figure BDA0004028157890000068
是均值为零向量，协方差矩阵为I2M的高斯随机向量；
所述对N-1个修正幅相控制信号采用二元化处理，由公式(4)计算：
Figure BDA0004028157890000069
其中，
Figure BDA00040281578900000610
为
Figure BDA00040281578900000611
的第m个元素，∈是区间为(0,1)的均匀分布随机变量，M为数模转换单元采样位数；
所述二元化处理后修正幅相控制信号集合可表示为
Figure BDA00040281578900000612
Figure BDA00040281578900000613
Figure BDA00040281578900000614
为修正幅相控制信号。
在上述方案的基础上，步骤六中所述修正残余信号功率集合可表示为
Figure BDA0004028157890000071
Figure BDA0004028157890000072
为修正残余信号功率。
在上述方案的基础上，步骤七中所述依序比较N-1组当前残余信号功率与修正残余信号功率，每组比较结果取二者较小值，构成残余信号功率集合及对应的幅相控制信号集合，分别由公式(5)和(6)表示：
Figure BDA0004028157890000073
Figure BDA0004028157890000074
其中，
Figure BDA0004028157890000075
为比较后的残余信号功率集合，N为幅相控制信号个数，n*为步骤四中残余信号功率最小值下标，Pi为集合
Figure BDA0004028157890000076
的元素，可由公式(7)表示：
Figure BDA0004028157890000077
其中，min{·,·}为取最小值运算符，
Figure BDA0004028157890000078
为修正残余信号功率，
Figure BDA0004028157890000079
为比较后的幅相控制信号集合，xi为集合
Figure BDA00040281578900000710
的元素，可由公式(8)表示：
Figure BDA00040281578900000711
其中，
Figure BDA00040281578900000715
为指示函数，当约束条件C为真时，
Figure BDA00040281578900000716
否则，
Figure BDA00040281578900000717
在上述方案的基础上，步骤八中所述将比较所得N-1个残余信号功率集合与步骤四中残余信号功率最小值构成下一轮迭代的当前残余信号功率集合，可由公式(9)表示：
Figure BDA00040281578900000712
其中，
Figure BDA00040281578900000713
为当前残余信号功率集合，
Figure BDA00040281578900000714
为比较后的残余信号功率集合，
Figure BDA0004028157890000081
为步骤四中残余信号功率最小值；所述集合中残余信号功率对应的幅相控制信号构成下一轮迭代的当前幅相控制信号集合，可由公式(10)表示：
Figure BDA0004028157890000082
其中，
Figure BDA0004028157890000083
为当前幅相控制信号集合，
Figure BDA0004028157890000084
为比较后的幅相控制信号集合，
Figure BDA0004028157890000085
为步骤四中残余信号功率最小值对应的幅相控制信号。
在上述方案的基础上，步骤九中所述算法停止准则，包括：若t>tmax，则迭代计数器溢出，算法停止；否则，迭代计数器未溢出，算法继续迭代；
在上述方案的基础上，步骤十中所述从当前残余信号功率集合中选取功率最小值，将其对应的幅相位控制信号作为最优幅相控制信号，由公式(11)和(12)表示：
Figure BDA0004028157890000086
Figure BDA0004028157890000087
其中，
Figure BDA0004028157890000088
为当前残余信号功率集合，nopt为当前残余信号功率集合中选取功率最小值对应下标；xopt为输出的最优幅相控制信号。本发明的有益效果：
本系统以及方法可满足无线通信射频自干扰消除需求，在干扰信号重建的过程中，通过莱维随机游走进行全局搜索，能够快速、高效地搜索到全局最优的幅相控制信号，以重构与自干扰信号等幅且反向的对消信号，实施简单，通用性强，鲁棒性高，对提高无线通信性能具有重要意义。
附图说明
本发明有如下附图：
图1一种射频自干扰消除系统结构框图；
图2一种射频自干扰消除方法获得的残余信号功率与迭代轮数曲线仿真图；
图3一种射频自干扰消除方法流程框图。
具体实施方式
下面结合图1～2和具体实施例，进一步阐明本发明，应理解这些实施例仅用于说明本发明而不用于限制本发明的范围，在阅读了本发明之后，本领域技术人员对本发明的各种等价形式的修改均落于本申请所附权利要求所限定的范围。
一种射频自干扰消除系统如图1所示包括：发送信号耦合单元将发送端射频信号分为两路，其中一路输出至射频隔离单元1端口，经2端口输出至天线，另一路输出至幅相控制单元；射频隔离单元3端口的输出信号为2端口输入的接收射频信号及由1端口泄露至3端口的射频自干扰信号的叠加；信号合成单元将幅相控制单元输出的射频对消信号与射频隔离单元3端口输出的信号合并为一路残余信号，输出至接收信号耦合单元；接收信号耦合单元将信号合成单元的输出信号分为两路，其中一路输出至接收端，另一路输出至功率测量单元；功率测量单元测量接收信号耦合单元输出的残余信号功率，并将功率值输出至算法处理单元；算法处理单元执行射频自干扰消除算法，产生数字幅相控制信号，并输出至数模转换单元；数模转换单元将算法处理单元输出的数字幅相控制信号转换为模拟幅相控制信号，并输出至幅相控制单元；幅相控制单元根据模拟幅相控制信号，调节信号幅度和相位并输出射频对消信号至信号合成单元。
采用上述方案的射频自干扰消除方法，包括以下步骤：
步骤一：初始化标志位F＝1，迭代计数器t＝0，迭代计数器溢出上限tmax，N个幅相控制信号以及莱维随机游走参数集合；
步骤二：N个幅相控制信号构成当前幅相控制信号集合；依次利用N个幅相控制信号产生N个残余信号，得到当前残余信号功率集合；
步骤三：判断标志位F是否为1，若是，转至步骤四；否则，转至步骤九；
步骤四：从残余信号功率集合中选取最小值，记录其幅相控制信号；剩余N-1个残余信号功率及其对应的幅相控制信号分别构成集合；
步骤五：对剩余N-1个幅相控制信号分别执行莱维随机游走，产生N-1个修正幅相控制信号，采用二元化处理后构成集合；
步骤六：依次利用N-1个修正幅相控制信号产生N-1个修正残余信号，得到修正残余信号功率集合；
步骤七：依序比较N-1组当前残余信号功率与修正残余信号功率，每组比较结果取二者较小值，构成残余信号功率集合及对应的幅相控制信号集合；
步骤八：将比较所得残余信号功率集合与步骤四中残余信号功率最小值构成下一轮迭代的当前残余信号功率集合，集合中残余信号功率对应的幅相控制信号构成下一轮迭代的当前幅相控制信号集合，置标志位F＝2，t＝t+1，并转至步骤三；
步骤九：判断算法是否满足停止准则，若是，转至步骤十；否则，转至步骤十一；
步骤十：从当前残余信号功率集合中选取功率最小值，将其对应的幅相位控制信号作为最优幅相控制信号，并进行输出；
步骤十一：置标志位F＝1，转至步骤三。
下面以幅相控制信号数目N＝4，数模转换单元采样位数M＝2为例，对无线通信系统进行射频自干扰消除，包括：
所述步骤一中分别初始化4个数组长度为4的幅相控制信号为x1＝[1,1,0,1]，x2＝[0,0,1,1]，x3＝[1,1,0,0]和x4＝[0,1,0,0]，初始化莱维随机游走参数α＝1和β＝1.5，初始化迭代计数器数值t＝0，假设由4个幅相控制信号分别产生的残余信号功率为P1＝0.21，P2＝0.12，P3＝0.33和P4＝0.25。
所述步骤三中当前幅相控制信号集合可表示为
Figure BDA0004028157890000111
Figure BDA0004028157890000112
当前残余信号功率集合可表示为
Figure BDA0004028157890000113
Figure BDA0004028157890000114
所述步骤四中从残余信号功率集合中选取最小值P2，记录其幅相控制信号x2，其中，残余信号功率最小值下标由公式
Figure BDA0004028157890000121
Figure BDA0004028157890000122
计算。
所述步骤五中假设产生修正幅相控制信号的随机步长向量为s1＝[2.06,-0.13,-0.45,-0.65]，s2＝[-0.69,-0.64,1.00,0.91]，s3＝[0.48,0.15,1.87,-0.85]和s4＝[-0.55,-1.91,-0.74,0.31]；剩余3个残余信号功率集合可表示为
Figure BDA0004028157890000123
其对应幅相控制信号集合可表示为
Figure BDA0004028157890000124
利用公式(2)对集合
Figure BDA0004028157890000125
中的3个幅相控制信号分别执行莱维随机游走，产生3个修正幅相控制信号：
Figure BDA0004028157890000126
Figure BDA0004028157890000127
Figure BDA0004028157890000128
和
Figure BDA0004028157890000129
假设∈＝0.5，利用公式(4)对上述3个修正幅相控制信号分别采用二元化处理，得到
Figure BDA00040281578900001210
Figure BDA00040281578900001211
Figure BDA00040281578900001212
对应的残余信号功率假设为
Figure BDA00040281578900001213
和
Figure BDA00040281578900001214
Figure BDA00040281578900001215
二元化处理后修正幅相控制信号可表示为
Figure BDA00040281578900001216
Figure BDA00040281578900001217
所述步骤六中修正残余信号功率集合可表示为
Figure BDA00040281578900001218
Figure BDA00040281578900001219
所述步骤七中依序比较3组当前残余信号功率与修正残余信号功率，每组比较结果取二者较小值，构成残余信号功率集合
Figure BDA00040281578900001220
Figure BDA00040281578900001221
其中，Pi由公式(7)计算：P1＝P1，
Figure BDA00040281578900001222
及
Figure BDA00040281578900001223
Figure BDA00040281578900001224
对应的幅相控制信号集合
Figure BDA00040281578900001225
其中，xi由公式(8)计算：x1＝x1，
Figure BDA00040281578900001226
和
Figure BDA00040281578900001227
所述步骤八中将比较所得残余信号功率集合与步骤四中残余信号功率最小值构成下一轮迭代的当前残余信号功率集合，即
Figure BDA0004028157890000131
Figure BDA0004028157890000132
集合中残余信号功率对应的幅相控制信号构成下一轮迭代的当前幅相控制信号集合，即
Figure BDA0004028157890000133
所述步骤十中假设算法已满足停止准则，从当前残余信号功率集合中选取功率最小值，即
Figure BDA0004028157890000134
将其对应的幅相位控制信号作为最优幅相控制信号，即xopt＝x4，完成射频自干扰消除。
采用本发明获得残余信号功率与算法迭代轮数曲线仿真图如图2所示：
从附图2中可看出：1)残余信号功率均随着算法迭代轮数的增加而显著下降；2)在迭代至约400轮时，射频自干扰消除后的残余信号功率趋向于零；基于上述两点可认为本方法有效解决无线通信系统射频自干扰消除问题。
本说明书中未作详细描述的内容属于本领域专业技术人员公知的现有技术。

Claims (10)

1.一种射频自干扰消除系统，其特征在于，包括发送信号耦合单元、射频隔离单元、信号合成单元、接收信号耦合单元、功率检测单元、算法处理单元、数模转换单元及幅相控制单元；
所述发送信号耦合单元用于将发送端射频信号按比例分成两路，分别输出至所述射频隔离单元及幅相控制单元；
所述射频隔离单元用于隔离发送端和接收端之间的射频信号；
所述信号合成单元用于将两路输入信号合并为一路输出信号；
所述接收信号耦合单元用于将所述信号合成单元输出的射频信号按比例分成两路，分别输出至接收端及所述功率测量单元；
所述功率测量单元用于测量并获取所输入的合并信号功率值；
所述算法处理单元用于执行所设计的射频自干扰消除方法；
所述数模转换单元用于将所述算法处理单元输出的数字幅相控制信号转换为模拟幅相控制信号；
所述幅相控制单元用于接收所述数模转换单元输出的模拟幅相控制信号，并调整信号幅度和相位，输出射频对消信号至所述信号合成单元。
2.根据权利要求1所述的一种射频自干扰消除系统，其特征在于，
所述发送信号耦合单元将发送端射频信号分为两路，其中一路输出至射频隔离单元第一端口，经射频隔离单元第二端口输出至天线，另一路输出至所述幅相控制单元；
射频隔离单元第三端口的输出信号为射频隔离单元第二端口输入的接收射频信号及由射频隔离单元第一端口泄露至射频隔离单元第三端口的射频自干扰信号的叠加；
所述信号合成单元将所述幅相控制单元输出的射频对消信号与射频隔离单元第三端口输出的信号合并为一路残余信号，输出至所述接收信号耦合单元；
所述接收信号耦合单元将信号合成单元的输出信号分为两路，其中一路输出至接收端，另一路输出至所述功率测量单元；
所述功率测量单元测量接收信号耦合单元输出的残余信号功率，并将功率值输出至所述算法处理单元；
所述算法处理单元执行射频自干扰消除算法，产生数字幅相控制信号，并输出至所述数模转换单元；
所述数模转换单元将算法处理单元输出的数字幅相控制信号转换为模拟幅相控制信号，并输出至所述幅相控制单元；
所述幅相控制单元根据模拟幅相控制信号，调节信号幅度和相位并输出射频对消信号至所述信号合成单元。
3.一种适用于射频自干扰消除系统的射频自干扰消除方法，具体包括如下步骤：
步骤一：初始化标志位F＝1，迭代计数器t＝0，迭代计数器溢出上限tmax，N个幅相控制信号以及莱维随机游走参数集合；
步骤二：N个幅相控制信号构成当前幅相控制信号集合；依次利用N个幅相控制信号产生N个残余信号，得到当前残余信号功率集合；
步骤三：判断标志位F是否为1，若是，转至步骤四；否则，转至步骤九；
步骤四：从残余信号功率集合中选取最小值，记录其幅相控制信号；剩余N-1个残余信号功率及其对应的幅相控制信号分别构成集合；
步骤五：对剩余N-1个幅相控制信号分别执行莱维随机游走，产生N-1个修正幅相控制信号，采用二元化处理后构成集合；
步骤六：依次利用N-1个修正幅相控制信号产生N-1个修正残余信号，得到修正残余信号功率集合；
步骤七：依序比较N-1组当前残余信号功率与修正残余信号功率，每组比较结果取二者较小值，构成残余信号功率集合及对应的幅相控制信号集合；
步骤八：将比较所得残余信号功率集合与步骤四中残余信号功率最小值构成下一轮迭代的当前残余信号功率集合，集合中残余信号功率对应的幅相控制信号构成下一轮迭代的当前幅相控制信号集合，置标志位F＝2，t＝t+1，并转至步骤三；
步骤九：判断算法是否满足停止准则，若是，转至步骤十；否则，转至步骤十一；
步骤十：从当前残余信号功率集合中选取功率最小值，将其对应的幅相位控制信号作为最优幅相控制信号，并进行输出；
步骤十一：置标志位F＝1，转至步骤三。
4.根据权利要求3所述的一种适用于射频自干扰消除系统的射频自干扰消除方法，其特征在于，步骤一中所述幅相控制信号由长度为2M的二进制数组xn，1≤n≤N，表示，其中，M为数模转换单元采样位数，二进制数组xn的前M位表示幅度控制信号，后M位表示相位控制信号；
所述莱维随机游走参数集由二元组<α,β>表示，其中，α为步长缩放因子，β为莱维指数。
5.根据权利要求3所述的一种适用于射频自干扰消除系统的射频自干扰消除方法，其特征在于，步骤三中所述当前幅相控制信号集合表示为
Figure FDA0004028157880000041
当前残余信号功率集合表示为
Figure FDA0004028157880000042
其中，Pn为第n个残余信号的功率。
6.根据权利要求3所述的一种适用于射频自干扰消除系统的射频自干扰消除方法，其特征在于，步骤四中所述从残余信号功率集合中选取最小值
Figure FDA0004028157880000043
记录其幅相控制信号
Figure FDA0004028157880000044
其中，n*为残余信号功率最小值下标，由公式(1)计算：
Figure FDA0004028157880000045
其中，argmin(·)(·)为取集合中最小元素下标运算符，
Figure FDA0004028157880000046
为当前残余信号功率集合；
所述剩余N-1个残余信号功率集合可表示为
Figure FDA0004028157880000047
\为集合差集运算符；其对应的幅相控制信号集合可表示为
Figure FDA0004028157880000048
Figure FDA0004028157880000049
7.根据权利要求3所述的一种适用于射频自干扰消除系统的射频自干扰消除方法，其特征在于，步骤五中所述对剩余N-1个幅相控制信号分别执行莱维随机游走，产生N-1个修正幅相控制信号，可由公式(2)表示：
Figure FDA0004028157880000051
其中，xi为集合
Figure FDA0004028157880000052
中的幅相控制信号，
Figure FDA0004028157880000053
为修正幅相控制信号，α为步长缩放因子，
Figure FDA0004028157880000054
为修正幅相控制信号的随机步长向量，si服从莱维分布，由公式(3)产生：
Figure FDA0004028157880000055
其中，⊙为哈达玛积运算符，|·|为绝对值运算符，
Figure FDA0004028157880000056
为均值为零向量协方差矩阵为σ2I2M的高斯随机向量，
Figure FDA0004028157880000057
IQ为Q阶单位阵，β为莱维指数，Γ(·)为伽马函数，sin(·)为正弦函数，π为圆周率，
Figure FDA0004028157880000058
是均值为零向量，协方差矩阵为I2M的高斯随机向量；
对N-1个修正幅相控制信号采用二元化处理，由公式(4)计算：
Figure FDA0004028157880000059
其中，
Figure FDA00040281578800000510
为
Figure FDA00040281578800000511
的第m个元素，∈是区间为(0,1)的均匀分布随机变量，M为数模转换单元采样位数；
所述二元化处理后修正幅相控制信号集合可表示为
Figure FDA00040281578800000512
Figure FDA00040281578800000513
Figure FDA00040281578800000514
为修正幅相控制信号。
8.根据权利要求3所述的一种适用于射频自干扰消除系统的射频自干扰消除方法，其特征在于，步骤六中所述修正残余信号功率集合可表示为
Figure FDA0004028157880000061
Pi levy为修正残余信号功率。
9.根据权利要求3所述的一种适用于射频自干扰消除系统的射频自干扰消除方法，其特征在于，步骤七中所述依序比较N-1组当前残余信号功率与修正残余信号功率，每组比较结果取二者较小值，构成残余信号功率集合及对应的幅相控制信号集合，分别由公式(5)和(6)表示：
Figure FDA0004028157880000062
Figure FDA0004028157880000063
其中，
Figure FDA0004028157880000064
为比较后的残余信号功率集合，N为幅相控制信号个数，n*为步骤四中残余信号功率最小值下标，Pi为集合
Figure FDA0004028157880000065
的元素，可由公式(7)表示：
Pi＝min{Pi,Pi levy}                             (7)，
其中，min{·,·}为取最小值运算符，Pi levy为修正残余信号功率，
Figure FDA0004028157880000066
为比较后的幅相控制信号集合，xi为集合
Figure FDA0004028157880000067
的元素，可由公式(8)表示：
Figure FDA0004028157880000068
其中，
Figure FDA0004028157880000069
为指示函数，当约束条件C为真时，
Figure FDA00040281578800000610
否则，
Figure FDA00040281578800000611
步骤八中所述将比较所得N-1个残余信号功率集合与步骤四中残余信号功率最小值构成下一轮迭代的当前残余信号功率集合，可由公式(9)表示：
Figure FDA0004028157880000071
其中，
Figure FDA0004028157880000072
为当前残余信号功率集合，
Figure FDA0004028157880000073
为比较后的残余信号功率集合，
Figure FDA0004028157880000074
为步骤四中残余信号功率最小值；所述集合中残余信号功率对应的幅相控制信号构成下一轮迭代的当前幅相控制信号集合，可由公式(10)表示：
Figure FDA0004028157880000075
其中，
Figure FDA0004028157880000076
为当前幅相控制信号集合，
Figure FDA0004028157880000077
为比较后的幅相控制信号集合，
Figure FDA0004028157880000078
为步骤四中残余信号功率最小值对应的幅相控制信号。
10.根据权利要求3所述的一种适用于射频自干扰消除系统的射频自干扰消除方法，其特征在于，步骤九中所述算法停止准则，包括：若t>tmax，则迭代计数器溢出，算法停止；否则，迭代计数器未溢出，算法继续迭代；
步骤十中所述从当前残余信号功率集合中选取功率最小值，将其对应的幅相位控制信号作为最优幅相控制信号，由公式(11)和(12)表示：
Figure FDA0004028157880000079
Figure FDA00040281578800000710
其中，
Figure FDA00040281578800000711
为当前残余信号功率集合，nopt为当前残余信号功率集合中选取功率最小值对应下标；xopt为输出的最优幅相控制信号。
CN202211728329.0A 2022-12-30 2022-12-30 一种射频自干扰消除系统及方法 Pending CN116073853A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211728329.0A CN116073853A (zh) 2022-12-30 2022-12-30 一种射频自干扰消除系统及方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211728329.0A CN116073853A (zh) 2022-12-30 2022-12-30 一种射频自干扰消除系统及方法

Publications (1)

Publication Number Publication Date
CN116073853A true CN116073853A (zh) 2023-05-05

Family

ID=86177999

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211728329.0A Pending CN116073853A (zh) 2022-12-30 2022-12-30 一种射频自干扰消除系统及方法

Country Status (1)

Country Link
CN (1) CN116073853A (zh)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103916148A (zh) * 2013-01-05 2014-07-09 华为技术有限公司 一种自适应射频干扰抵消装置、方法以及接收机
CN104935540A (zh) * 2015-05-21 2015-09-23 电子科技大学 一种同时同频全双工极限自干扰抵消方法
US20160057641A1 (en) * 2013-03-22 2016-02-25 Telefonaktiebolaget L M Ericsson (Publ) Interference Estimation and Devices Therefore
CN105812083A (zh) * 2016-04-13 2016-07-27 电子科技大学 一种同时同频全双工系统中的射频快速自干扰抵消方法
CN206389365U (zh) * 2017-01-25 2017-08-08 北京科技大学 一种多抽头非等功分全双工模拟自干扰消除装置
CN107846236A (zh) * 2017-11-07 2018-03-27 成都信息工程大学 消除全双工自干扰抑制过程中射频泄露信号影响的方法
US20180357997A1 (en) * 2017-06-09 2018-12-13 Foster-Miller, Inc. Acoustic and rf cancellation systems and methods
CN115441898A (zh) * 2022-09-02 2022-12-06 北京理工大学 一种多抽头结构的全双工模拟域自干扰消除系统

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103916148A (zh) * 2013-01-05 2014-07-09 华为技术有限公司 一种自适应射频干扰抵消装置、方法以及接收机
US20160057641A1 (en) * 2013-03-22 2016-02-25 Telefonaktiebolaget L M Ericsson (Publ) Interference Estimation and Devices Therefore
CN104935540A (zh) * 2015-05-21 2015-09-23 电子科技大学 一种同时同频全双工极限自干扰抵消方法
CN105812083A (zh) * 2016-04-13 2016-07-27 电子科技大学 一种同时同频全双工系统中的射频快速自干扰抵消方法
CN206389365U (zh) * 2017-01-25 2017-08-08 北京科技大学 一种多抽头非等功分全双工模拟自干扰消除装置
US20180357997A1 (en) * 2017-06-09 2018-12-13 Foster-Miller, Inc. Acoustic and rf cancellation systems and methods
CN107846236A (zh) * 2017-11-07 2018-03-27 成都信息工程大学 消除全双工自干扰抑制过程中射频泄露信号影响的方法
CN115441898A (zh) * 2022-09-02 2022-12-06 北京理工大学 一种多抽头结构的全双工模拟域自干扰消除系统

Similar Documents

Publication Publication Date Title
CN1732634B (zh) 一种用于自相关滤波的系统和方法
CN110113279B (zh) 一种移动跳频水声通信多普勒因子估计方法
JP2013500662A (ja) 特にマルチスタンダードなソフトウェア無線、および／またはコグニティブ無線の使用のための並列アナログ−デジタル変換器中のアナログ欠陥の訂正方法
CN109061296B (zh) 一种射频脉冲信号的高精度载频估计方法
Li et al. Distributed MIMO radar based on sparse sensing: Analysis and efficient implementation
CN112180320A (zh) 一种无人机无源定位系统及方法
US5793797A (en) Data transmisson system with a low peak-to-average power ratio based on distorting small amplitude signals
JP2021166362A (ja) 信号推定装置、信号推定方法、及び、コンピュータプログラム
CN103701515A (zh) 一种数字多波束形成方法
CN116073853A (zh) 一种射频自干扰消除系统及方法
US7830949B2 (en) Cross correlation circuits and methods
Kong A novel and efficient mixed-signal compressed sensing for wide-band cognitive radio
CN111934651A (zh) 一种空时频级联自适应滤波处理方法、装置及设备
Yang et al. Minimax and WLS designs of digital FIR filters using SOCP for aliasing errors reduction in BI-DAC
Anderson et al. Real-number convolutional codes for speech-like quasi-stationary sources (Corresp.)
CN217181205U (zh) 一种宽带dbf侦察通道校准系统
KR100506487B1 (ko) 단일 주파수, 다중 송신 통신망의 코드벡터 검출 장치
Khakpour et al. Adaptive noise cancellation based architecture for correction of gain and offset mismatch errors in time-interleaved ADC
CN113612571A (zh) 一种基于广义多分数傅里叶变换的多分量安全传输方法
Beydoun et al. Optimal digital reconstruction and calibration for multichannel Time Interleaved ΣΔ ADC based on Comb-filters
CN106341102B (zh) 一种用于宽带系统的数字式相位非线性校正方法及装置
KR100769955B1 (ko) 무선 통신 시스템 수신단의 검파 장치
CN114826347B (zh) 无线通信系统的波束成形方法、装置及存储介质
CN113300728B (zh) 一种基于非正交多址上行链路的自适应消除干扰方法
CN117607916B (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