[go: up one dir, main page]

CN116405129A - Self-interference cancellation method, system, device and medium based on improved LMS algorithm - Google Patents

Self-interference cancellation method, system, device and medium based on improved LMS algorithm Download PDF

Info

Publication number
CN116405129A
CN116405129A CN202310293398.1A CN202310293398A CN116405129A CN 116405129 A CN116405129 A CN 116405129A CN 202310293398 A CN202310293398 A CN 202310293398A CN 116405129 A CN116405129 A CN 116405129A
Authority
CN
China
Prior art keywords
signal
interference cancellation
self
error
function
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
CN202310293398.1A
Other languages
Chinese (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.)
Sun Yat Sen University
Original Assignee
Sun Yat Sen 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 Sun Yat Sen University filed Critical Sun Yat Sen University
Priority to CN202310293398.1A priority Critical patent/CN116405129A/en
Publication of CN116405129A publication Critical patent/CN116405129A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B15/00Suppression or limitation of noise or interference
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/15Correlation function computation including computation of convolution operations
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L5/00Arrangements affording multiple use of the transmission path
    • H04L5/14Two-way operation using the same type of signal, i.e. duplex
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02DCLIMATE 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/00Reducing energy consumption in communication networks
    • Y02D30/70Reducing energy consumption in communication networks in wireless communication networks

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Signal Processing (AREA)
  • Mathematical Optimization (AREA)
  • Data Mining & Analysis (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Pure & Applied Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Computing Systems (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Algebra (AREA)
  • Noise Elimination (AREA)

Abstract

本发明涉及数字信号处理技术领域,公开了一种基于改进的LMS算法的自干扰抵消方法、系统、设备和介质,其中方法包括:根据全双工ISAC系统发送端所发送的发送信号及滤波器的权重系数,计算所述滤波器的输出信号;根据所述全双工ISAC系统中基带接收信号及所述输出信号进行误差计算,得到误差信号;根据所述误差信号和改进的箕舌线函数,构造步长调整函数,并根据所述步长调整函数对所述权重系数进行迭代更新;根据更新后的权重系数对误差信号进行迭代更新,当更新后的误差信号的功率达到最小值时,输出更新后的误差信号及其对应的输出信号。本发明的改进算法使得全双工ISAC系统的数字域自干扰抵消过程可在保持较快收敛性的同时还能保持较低的计算复杂度。

Figure 202310293398

The present invention relates to the technical field of digital signal processing, and discloses a self-interference cancellation method, system, equipment and medium based on an improved LMS algorithm, wherein the method includes: according to the transmission signal and filter transmitted by the transmission end of the full-duplex ISAC system Calculate the output signal of the filter; carry out error calculation according to the baseband receiving signal and the output signal in the full-duplex ISAC system, and obtain the error signal; according to the error signal and the improved skip line function , construct a step adjustment function, and iteratively update the weight coefficient according to the step adjustment function; iteratively update the error signal according to the updated weight coefficient, when the power of the updated error signal reaches a minimum value, The updated error signal and its corresponding output signal are output. The improved algorithm of the invention enables the self-interference cancellation process in the digital domain of the full-duplex ISAC system to maintain relatively fast convergence while maintaining low computational complexity.

Figure 202310293398

Description

基于改进的LMS算法的自干扰抵消方法、系统、设备和介质Self-interference cancellation method, system, device and medium based on improved LMS algorithm

技术领域technical field

本发明涉及数字信号处理技术领域,特别是涉及一种基于改进的LMS算法的自干扰抵消方法、系统、设备和介质。The present invention relates to the technical field of digital signal processing, in particular to a self-interference cancellation method, system, device and medium based on an improved LMS algorithm.

背景技术Background technique

目前,最小均方算法——LMS算法,由于实现简单,被广泛应用于全双工ISAC系统的数字域自干扰消除中,但是对于传统固定步长算法,当步长因子选取较大值时,算法可以获取较快的收敛速度,但是稳态误差较大;当步长因子选取较小值时,算法可以获取较小的稳态误差,但收敛速度较慢。固定步长值固有缺陷无法合理协调收敛速度和稳态误差,各种变步长LMS算法陆续出现:如由覃景繁等人基于Sigmoid函数提出的变步长LMS算法——SVSLMS算法,通过在误差信号和步长之间建立新的非线性函数关系,克服了LMS算法收敛速度与稳态误差之间的矛盾,但SVSLMS算法在误差接近零处时,变步长函数变化太大,不具有缓慢变化的特性,使得SVSLMS算法在自适应稳态阶段仍具有较大的步长变化;再如利用箕舌线函数的特性提出的基于箕舌线的变步长LMS算法——TCLMS算法,该算法计算复杂度低,且在稳态误差与前述SVSLMS算法相同的前提下,具有更快的收敛速度和跟踪速度,更大程度上解决了收敛速度和稳态误差的矛盾,但由于目前应用的要求,TCLMS算法的收敛速度还需进一步提高,并且其容易受信号输入端不相关噪声的干扰而影响算法的稳定性。因此,亟需一种算法可在全双工ISAC系统的数字域自干扰消除过程中,保持较快收敛性的同时还能保持较低的计算复杂度。At present, the least mean square algorithm——LMS algorithm is widely used in digital domain self-interference cancellation of full-duplex ISAC systems due to its simple implementation. The algorithm can obtain a faster convergence speed, but the steady-state error is larger; when the step size factor is smaller, the algorithm can obtain a smaller steady-state error, but the convergence speed is slower. The inherent defect of the fixed step value cannot reasonably coordinate the convergence speed and steady-state error, and various variable-step-size LMS algorithms have emerged one after another: such as the variable-step-size LMS algorithm proposed by Qin Jingfan et al. based on the Sigmoid function——SVSLMS algorithm. Establish a new nonlinear function relationship between the LMS algorithm and the step size, which overcomes the contradiction between the convergence speed of the LMS algorithm and the steady-state error, but when the error is close to zero in the SVSLMS algorithm, the variable step size function changes too much and does not have a slow change The characteristics of the SVSLMS algorithm still have a large step size change in the adaptive steady-state stage; another example is the TCLMS algorithm based on the variable step length LMS algorithm proposed by using the characteristics of the skip line function, which calculates The complexity is low, and under the premise that the steady-state error is the same as the aforementioned SVSLMS algorithm, it has faster convergence speed and tracking speed, and solves the contradiction between convergence speed and steady-state error to a greater extent. However, due to the current application requirements, The convergence speed of the TCLMS algorithm needs to be further improved, and it is easily affected by the interference of uncorrelated noise at the signal input end, which affects the stability of the algorithm. Therefore, there is an urgent need for an algorithm that can maintain fast convergence and low computational complexity in the process of self-interference cancellation in the digital domain of a full-duplex ISAC system.

发明内容Contents of the invention

本发明提供一种基于改进的LMS算法的自干扰抵消方法、系统、设备和介质,通过对算法进行改进,使得全双工ISAC系统数字域自干扰抵消过程即使处于低信噪比状态下,也仍能保持良好的性能,并在保持较快收敛性的同时还能保持较低的计算复杂度。The present invention provides a self-interference cancellation method, system, equipment and medium based on an improved LMS algorithm. By improving the algorithm, the self-interference cancellation process in the digital domain of a full-duplex ISAC system can be performed even in a low signal-to-noise ratio state. It still maintains good performance and maintains low computational complexity while maintaining fast convergence.

为达到上述目的,本发明第一方面提供一种基于改进的LMS算法的自干扰抵消方法,包括以下步骤:In order to achieve the above object, the first aspect of the present invention provides a self-interference cancellation method based on the improved LMS algorithm, comprising the following steps:

根据全双工ISAC系统发送端所发送的发送信号及滤波器的权重系数,计算所述滤波器的输出信号;calculating the output signal of the filter according to the transmission signal sent by the transmitting end of the full-duplex ISAC system and the weight coefficient of the filter;

根据所述全双工ISAC系统中基带接收信号及所述输出信号进行误差计算,得到误差信号;Performing error calculation according to the baseband receiving signal and the output signal in the full-duplex ISAC system to obtain an error signal;

根据所述误差信号和改进的箕舌线函数,构造步长调整函数,并根据所述步长调整函数对所述权重系数进行迭代更新;Construct a step size adjustment function according to the error signal and the improved skip line function, and iteratively update the weight coefficient according to the step size adjustment function;

根据更新后的权重系数对误差信号进行迭代更新,当更新后的误差信号的功率达到最小值时,输出更新后的误差信号及其对应的输出信号。The error signal is iteratively updated according to the updated weight coefficient, and when the power of the updated error signal reaches a minimum value, the updated error signal and its corresponding output signal are output.

进一步地,在所述根据全双工ISAC系统发送端所发送的发送信号及滤波器的权重系数,计算所述滤波器的输出信号之前,包括:Further, before calculating the output signal of the filter according to the transmission signal sent by the transmitting end of the full-duplex ISAC system and the weight coefficient of the filter, the method includes:

对所述滤波器进行初始化,使所述权重系数为0。The filter is initialized so that the weight coefficient is 0.

进一步地,所述根据所述误差信号和改进的箕舌线函数,构造步长调整函数,包括:Further, the step size adjustment function is constructed according to the error signal and the improved skip line function, including:

将所述误差函数作为自变量代入改进的箕舌线函数中,根据步长调整原则,对所述误差信号添加绝对值并引入调整参数,得到所述步长调整函数。Substituting the error function as an independent variable into the improved skip line function, according to the principle of step size adjustment, adding an absolute value to the error signal and introducing adjustment parameters to obtain the step size adjustment function.

进一步地,所述步长调整函数的表达式为:Further, the expression of the step size adjustment function is:

Figure BDA0004142398750000021
Figure BDA0004142398750000021

式中,μ(n)为步长调整函数,α为第一调整参数,β为第二调整参数,m为第三调整参数,且α、β、m均大于0,e(n)为误差函数。In the formula, μ(n) is the step adjustment function, α is the first adjustment parameter, β is the second adjustment parameter, m is the third adjustment parameter, and α, β, m are all greater than 0, and e(n) is the error function.

进一步地,所述基带接收信号为所述全双工ISAC系统中经过天线域自干扰抵消、模拟域自干扰抵消及ADC量化后的信号,所述基带接收信号的表达式为:Further, the baseband received signal is a signal after antenna domain self-interference cancellation, analog domain self-interference cancellation and ADC quantization in the full-duplex ISAC system, and the expression of the baseband received signal is:

r(n)=sI(n)+d(n)+ε(n)r(n)=s I (n)+d(n)+ε(n)

式中,r(n)为基带接收信号,sI(n)为天线域自干扰抵消及模拟域自干扰抵消后的残余自干扰信号,d(n)为远端期望信号,ε(n)为加性噪声。In the formula, r(n) is the baseband received signal, s I (n) is the residual self-interference signal after antenna domain self-interference cancellation and analog domain self-interference cancellation, d(n) is the desired signal at the far end, ε(n) is additive noise.

进一步地,当干扰抵消比达到最大值时,结束迭代更新过程,并输出更新后的误差信号及其对应的输出信号;其中,所述干扰抵消比的计算公式如下:Further, when the interference cancellation ratio reaches the maximum value, the iterative update process is ended, and the updated error signal and its corresponding output signal are output; wherein, the calculation formula of the interference cancellation ratio is as follows:

Figure BDA0004142398750000031
Figure BDA0004142398750000031

式中,ICR为干扰抵消比,Pr为基带接收信号的功率,Pd为远端期望信号的功率,σε 2为加性噪声信号的功率,E{|e(∞)|2}为n趋于无穷大时的均方误差。In the formula, ICR is the interference cancellation ratio, P r is the power of baseband received signal, P d is the power of remote desired signal, σ ε 2 is the power of additive noise signal, E{|e(∞)| 2 } is The mean square error as n tends to infinity.

本发明第二方面提供一种基于改进的LMS算法的自干扰抵消系统,包括:The second aspect of the present invention provides a self-interference cancellation system based on the improved LMS algorithm, including:

输出信号获取模块,用于根据全双工ISAC系统发送端所发送的发送信号及滤波器的权重系数,计算所述滤波器的输出信号;The output signal acquisition module is used to calculate the output signal of the filter according to the transmission signal sent by the full-duplex ISAC system transmission end and the weight coefficient of the filter;

误差信号获取模块,用于根据所述全双工ISAC系统中基带接收信号及所述输出信号进行误差计算,得到误差信号;An error signal acquisition module, configured to perform error calculation according to the baseband receiving signal and the output signal in the full-duplex ISAC system, to obtain an error signal;

权重系数更新模块,用于根据所述误差信号和改进的箕舌线函数,构造步长调整函数,并根据所述步长调整函数对所述权重系数进行迭代更新;A weight coefficient updating module, configured to construct a step adjustment function according to the error signal and the improved skip line function, and iteratively update the weight coefficient according to the step adjustment function;

循环截止输出模块,用于根据更新后的权重系数对误差信号进行迭代更新,当更新后的误差信号的功率达到最小值时,输出更新后的误差信号及其对应的输出信号。The loop cutoff output module is configured to iteratively update the error signal according to the updated weight coefficient, and output the updated error signal and its corresponding output signal when the power of the updated error signal reaches a minimum value.

进一步地,所述权重系数更新模块包:Further, the weight coefficient update module package:

步长函数构造模块,用于将所述误差函数作为自变量代入改进的箕舌线函数中,根据步长调整原则,对所述误差信号添加绝对值并引入步长调整参数,得到所述步长调整函数。The step size function construction module is used for substituting the error function into the improved skip line function as an independent variable, according to the step size adjustment principle, adding an absolute value to the error signal and introducing a step size adjustment parameter to obtain the step size Long adjustment function.

本发明第三方面提供一种电子装置,包括处理器、存储器以及存储在所述存储器中且被配置为由所述处理器执行的计算机程序,所述处理器执行所述计算机程序时实现如上述第一方面中任意一项所述的基于改进的LMS算法的自干扰抵消方法。A third aspect of the present invention provides an electronic device, including a processor, a memory, and a computer program stored in the memory and configured to be executed by the processor. When the processor executes the computer program, the above-mentioned The self-interference cancellation method based on the improved LMS algorithm described in any one of the first aspect.

本发明第四方面提供一种计算机可读存储介质,所述计算机可读存储介质包括存储的计算机程序,其中,在所述计算机程序运行时控制所述计算机可读存储介质所在设备执行如上述第一方面中任意一项所述的基于改进的LMS算法的自干扰抵消方法。The fourth aspect of the present invention provides a computer-readable storage medium, the computer-readable storage medium includes a stored computer program, wherein, when the computer program is running, the device where the computer-readable storage medium is located is controlled to execute the above-mentioned first The self-interference cancellation method based on the improved LMS algorithm described in any one of the aspects.

与现有技术相比,本发明实施例的有益效果在于:Compared with the prior art, the beneficial effects of the embodiments of the present invention are:

本发明提供一种基于改进的LMS算法的自干扰抵消方法、系统、设备和介质,通过改进的箕舌线函数与误差信号建立新的非线性关系来构造步长调整函数,进而更新权重系数的步长因子,可在算法收敛过程中动态地改变步长因子的大小,使得改进的LMS算法在保持较快收敛性的同时还能保持较低的计算复杂度,且在低信噪比条件下,仍能保持良好的性能,并在全双工通感一体化数字域自干扰抵消中,在数字域的ICR高于其他比较的算法且具有较快的收敛速度。The present invention provides a self-interference cancellation method, system, device and medium based on the improved LMS algorithm, and constructs a step size adjustment function by establishing a new nonlinear relationship between the improved skip line function and the error signal, and then updates the weight coefficient The step size factor can dynamically change the size of the step size factor during the convergence process of the algorithm, so that the improved LMS algorithm can maintain a fast convergence while maintaining a low computational complexity, and it can be used under the condition of low signal-to-noise ratio , can still maintain good performance, and in the full-duplex synaesthesia integrated digital domain self-interference cancellation, the ICR in the digital domain is higher than other compared algorithms and has a faster convergence speed.

附图说明Description of drawings

为了更清楚地说明本发明的技术方案,下面将对实施方式中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施方式,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。In order to illustrate the technical solution of the present invention more clearly, the accompanying drawings used in the implementation will be briefly introduced below. Obviously, the accompanying drawings in the following description are only some implementations of the present invention. As far as the skilled person is concerned, other drawings can also be obtained based on these drawings on the premise of not paying creative work.

图1是本发明某一实施例提供的全双工ISAC系统模型图;Fig. 1 is a full-duplex ISAC system model diagram provided by a certain embodiment of the present invention;

图2是本发明某一实施例提供的基于LMS自适应滤波的自干扰抵消模型图;Fig. 2 is a self-interference cancellation model diagram based on LMS adaptive filtering provided by an embodiment of the present invention;

图3是本发明某一实施例提供的一种基于改进的LMS算法的自干扰抵消方法的流程图;Fig. 3 is a flowchart of a self-interference cancellation method based on an improved LMS algorithm provided by an embodiment of the present invention;

图4是本发明某一实施例提供的IVSSLMS、SVSLMS和TCLMS算法的步长更新变化的性能模拟结果图;Fig. 4 is the performance simulation result diagram of the step size update change of the IVSSLMS, SVSLMS and TCLMS algorithms provided by a certain embodiment of the present invention;

图5是本发明某一实施例提供的不同变步长LMS算法下的均方误差性能模拟结果图;Fig. 5 is the mean square error performance simulation result diagram under different variable step size LMS algorithms provided by a certain embodiment of the present invention;

图6是本发明某一实施例提供的信噪比为5时,不同变步长LMS算法下的均方误差性能模拟结果图;Fig. 6 is when the signal-to-noise ratio provided by a certain embodiment of the present invention is 5, the mean square error performance simulation result diagram under different variable step size LMS algorithms;

图7是本发明某一实施例提供的信噪比为1时,不同变步长LMS算法下的均方误差性能模拟结果图;Fig. 7 is the mean square error performance simulation result diagram under different variable step size LMS algorithms when the signal-to-noise ratio provided by a certain embodiment of the present invention is 1;

图8是本发明某一实施例提供的不同变步长LMS算法下的ICR对比图;Fig. 8 is a comparison diagram of ICR under different variable step size LMS algorithms provided by a certain embodiment of the present invention;

图9是本发明某一实施例提供的参数β、m不变,α分别取0.5,1,2,5时的步长因子与误差的函数模型曲线图;Fig. 9 is a function model curve diagram of the step factor and the error when the parameters β and m provided by a certain embodiment of the present invention are constant, and α is 0.5, 1, 2, and 5 respectively;

图10是本发明某一实施例提供的参数α、β不变,m分别取0.5,1,2,5时的步长因子与误差的函数模型曲线图;Fig. 10 is a function model curve diagram of the step factor and the error when the parameters α and β provided by a certain embodiment of the present invention are constant, and m is 0.5, 1, 2, and 5 respectively;

图11是本发明某一实施例提供的参数α、m不变,β分别取0.05,0.1,0.2,0.5时的步长因子与误差的函数模型曲线图;Fig. 11 is a function model curve diagram of the step factor and the error when the parameters α and m provided by a certain embodiment of the present invention are constant, and β is 0.05, 0.1, 0.2, and 0.5 respectively;

图12是一种基于改进的LMS算法的自干扰抵消系统的装置图;Fig. 12 is a device diagram of a self-interference cancellation system based on an improved LMS algorithm;

图13是本发明某一实施例提供的一种电子设备的结构图。Fig. 13 is a structural diagram of an electronic device provided by an embodiment of the present invention.

具体实施方式Detailed ways

下面结合附图和实施例,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。The technical solutions in the embodiments of the present invention will be clearly and completely described below in conjunction with the accompanying drawings and embodiments. Obviously, the described embodiments are only some of the embodiments of the present invention, not all of them. Based on the embodiments of the present invention, all other embodiments obtained by persons of ordinary skill in the art without creative efforts fall within the protection scope of the present invention.

应当理解,文中所使用的步骤编号仅是为了方便描述,不对作为对步骤执行先后顺序的限定。It should be understood that the step numbers used herein are only for convenience of description, and are not intended to limit the execution order of the steps.

应当理解,在本发明说明书中所使用的术语仅仅是出于描述特定实施例的目的而并不意在限制本发明。如在本发明说明书和所附权利要求书中所使用的那样,除非上下文清楚地指明其它情况,否则单数形式的“一”、“一个”及“该”意在包括复数形式。It should be understood that the terminology used in the description of the present invention is for the purpose of describing particular embodiments only and is not intended to limit the present invention. As used in this specification and the appended claims, the singular forms "a", "an" and "the" are intended to include plural referents unless the context clearly dictates otherwise.

术语“包括”和“包含”指示所描述特征、整体、步骤、操作、元素和/或组件的存在,但并不排除一个或多个其它特征、整体、步骤、操作、元素、组件和/或其集合的存在或添加。The terms "comprising" and "comprising" indicate the presence of described features, integers, steps, operations, elements and/or components, but do not exclude the presence of one or more other features, integers, steps, operations, elements, components and/or The presence or addition of its collection.

术语“和/或”是指相关联列出的项中的一个或多个的任何组合以及所有可能组合,并且包括这些组合。The term "and/or" refers to any combination of one or more of the associated listed items and all possible combinations, and includes such combinations.

本发明中涉及到的一些术语解释:Explanation of some terms involved in the present invention:

通感一体化:是指通信和感知两个功能融合在一起,使得通信系统同时具有通信和感知两个功能,在无线信道传输信息的同时通过主动认知并分析信道的特性,从而去感知周围环境的物理特征,从而使得通信与感知功能相互增强。Synaesthesia integration: refers to the fusion of communication and perception functions, so that the communication system has both communication and perception functions. While transmitting information through wireless channels, it actively recognizes and analyzes the characteristics of the channel to perceive the surroundings. The physical characteristics of the environment, so that communication and perception functions mutually enhance each other.

同时同频全双工:同频同时全双工技术被认为是一项有效提高频谱效率的技术,该技术是在同一个物理信道上实现两个方向信号的传输,即通过在通信双工节点的接收机处消除自身发射机信号的干扰,并在发射机信号处同时接收来自另一节点的同频信号。对比传统的时分双工(TDD)和频分双工(FDD)而言,同频同时全双工可以将频谱效率提高一倍。Simultaneous full-duplex at the same frequency: Simultaneous full-duplex at the same frequency is considered to be a technology that can effectively improve spectrum efficiency. This technology realizes the transmission of signals in two directions on the same physical channel, that is, through The receiver eliminates the interference of its own transmitter signal, and simultaneously receives the same frequency signal from another node at the transmitter signal. Compared with the traditional time division duplex (TDD) and frequency division duplex (FDD), simultaneous full duplex on the same frequency can double the spectrum efficiency.

LMS算法:最小均方算法,简称LMS算法,是一种最陡下降算法的改进算法,是在维纳滤波理论上运用速下降法后的优化延伸,最早是由Widrow和Hoff提出来的。该算法不需要已知输入信号和期望信号的统计特征,“当前时刻”的权系数是通过“上一时刻”权系数再加上一个负均方误差梯度的比例项求得。其具有计算复杂程度低、在信号为平稳信号的环境中收敛性好、其期望值无偏地收敛到维纳解和利用有限精度实现算法时的平稳性等特性,使LMS算法成为自适应算法中稳定性最好、应用最广的算法。LMS algorithm: The least mean square algorithm, referred to as the LMS algorithm, is an improved algorithm of the steepest descent algorithm. It is an optimized extension after using the rapid descent method in the Wiener filter theory. It was first proposed by Widrow and Hoff. The algorithm does not need to know the statistical characteristics of the input signal and the expected signal, and the weight coefficient of the "current moment" is obtained by adding a proportional term of the negative mean square error gradient to the weight coefficient of the "previous moment". It has the characteristics of low computational complexity, good convergence in the environment where the signal is a stationary signal, its expected value unbiasedly converges to the Wiener solution, and the stability of the algorithm when the algorithm is implemented with finite precision, making the LMS algorithm an adaptive algorithm. The most stable and widely used algorithm.

干扰抵消比:自干扰抵消后的干扰抵消比(Interference Cancellation Ratio,ICR),是判断干扰抵消能力的重要指标,ICR的值越大,表明自干扰抑制能力越好。Interference cancellation ratio: The interference cancellation ratio (Interference Cancellation Ratio, ICR) after self-interference cancellation is an important indicator for judging the interference cancellation ability. The larger the value of ICR, the better the self-interference suppression ability.

全双工ISAC系统模型如图1所示,由于全双工ISAC系统的特性,接收机在接收信号的同时,既会受到自干扰信号的影响,也会受到感知信号的影响,远端接收机亦是如此。本发明中全双工ISAC系统只考虑远端全双工通信感知节点的影响,不考虑感知目标的影响。由于全双工ISAC系统自干扰信号的功率远高于感兴趣信号的功率,因此它会使接收机饱和甚至损坏,进而使接收机无法接收所需信号,所以全双工ISAC系统最主要的挑战就是自干扰抵消。The full-duplex ISAC system model is shown in Figure 1. Due to the characteristics of the full-duplex ISAC system, the receiver will be affected by both the self-interference signal and the perception signal while receiving the signal. The remote receiver The same is true. In the present invention, the full-duplex ISAC system only considers the influence of the remote full-duplex communication sensing node, and does not consider the influence of the sensing target. The main challenge of the full-duplex ISAC system is that the power of the self-interfering signal in the full-duplex ISAC system is much higher than the power of the signal of interest, so it will saturate or even damage the receiver, which will make the receiver unable to receive the desired signal. It is self-interference cancellation.

全双工ISAC系统实现自干扰抵消分为三步:天线域、模拟域和数字域。当自干扰信号被接收端的接收天线接收时,通过天线自干扰抵消技术可以抵消一部分自干扰信号;完成天线域的自干扰抵消后,剩余的自干扰信号在进入ADC量化前,进行模拟域的自干扰抵消;完成ADC量化后的信号包括远端期望信号、部分自干扰信号和噪声信号,这些信号最终在基带解调前进行数字域自干扰抵消;至此,完成整个全双工ISAC系统的自干扰抵消过程。本发明所提出的基于改进的LMS算法的自干扰抵消方法即为全双工ISAC系统的数字域自干扰抵消,对基带接收信号中残余的自干扰信号进行消除。The self-interference cancellation of full-duplex ISAC system is divided into three steps: antenna domain, analog domain and digital domain. When the self-interference signal is received by the receiving antenna at the receiving end, a part of the self-interference signal can be canceled by the antenna self-interference cancellation technology; after the self-interference cancellation in the antenna domain is completed, the remaining self-interference signal is self-interferenced in the analog domain before being quantized by the ADC. Interference cancellation; the signal after ADC quantization includes the remote desired signal, some self-interference signals and noise signals, and these signals are finally subjected to digital domain self-interference cancellation before baseband demodulation; so far, the self-interference of the entire full-duplex ISAC system has been completed Offset process. The self-interference cancellation method based on the improved LMS algorithm proposed by the present invention is the digital domain self-interference cancellation of the full-duplex ISAC system, which eliminates the residual self-interference signal in the baseband received signal.

对于全双工ISAC系统模型的数字域自干扰抵消通常采用LMS自适应滤波法,而基于LMS自适应滤波的自干扰抵消模型如图2所示,该模型的数学描述如下:For the self-interference cancellation in the digital domain of the full-duplex ISAC system model, the LMS adaptive filtering method is usually used, and the self-interference cancellation model based on the LMS adaptive filtering is shown in Figure 2. The mathematical description of the model is as follows:

y(n)=sH(n)w(n)y(n)=s H (n)w(n)

e(n)=r(n)-y(n)e(n)=r(n)-y(n)

w(n+1)=w(n)+2μe(n)s(n)w(n+1)=w(n)+2μe(n)s(n)

式中,μ为步长因子,是一个常数;上标H表示矩阵的共轭转置;w(n)为M阶滤波器的权重系数,n为n时刻;s(n)为输入信号;y(n)为滤波器的输出信号;r(n)为基带接收信号;e(n)为误差信号,其用于向滤波器进行反馈,从而不断地调节滤波器的权重系数,在自干扰抵消模型中。但传统固定步长算法,当步长因子选取较大值时,算法可以获取较快的收敛速度,但是稳态误差较大;当步长因子选取较小值时,算法可以获取较小的稳态误差,但收敛速度较慢。固定步长值固有缺陷无法合理协调收敛速度和稳态误差,针对这个问题,各种变步长LMS算法被相继提出,其核心就是算法收敛过程中动态地改变步长因子的大小,本发明也就此问题提供了一种解决办法,通过改进的箕舌线函数与误差信号建立新的非线性关系来更新权重系数的步长因子。In the formula, μ is the step size factor, which is a constant; the superscript H represents the conjugate transposition of the matrix; w(n) is the weight coefficient of the M-order filter, n is the n moment; s(n) is the input signal; y(n) is the output signal of the filter; r(n) is the baseband receiving signal; e(n) is the error signal, which is used to feed back to the filter, so as to continuously adjust the weight coefficient of the filter. offset model. However, in the traditional fixed step size algorithm, when the step size factor is selected to a larger value, the algorithm can obtain a faster convergence speed, but the steady-state error is larger; when the step size factor is selected to a smaller value, the algorithm can obtain a smaller steady-state error. state error, but the convergence speed is slow. The inherent defect of the fixed step size value cannot reasonably coordinate the convergence speed and the steady-state error. Aiming at this problem, various variable step size LMS algorithms have been proposed one after another, the core of which is to dynamically change the size of the step size factor in the algorithm convergence process. The present invention also A solution to this problem is provided, and the step size factor of the weight coefficient is updated by establishing a new nonlinear relationship between the improved skip line function and the error signal.

在一实施例中,如图3所示,本发明第一方面提供一种基于改进的LMS算法的自干扰抵消方法,包括以下步骤:In one embodiment, as shown in FIG. 3 , the first aspect of the present invention provides a self-interference cancellation method based on an improved LMS algorithm, including the following steps:

S1、根据全双工ISAC系统发送端所发送的发送信号及滤波器的权重系数,计算滤波器的输出信号;S1. Calculate the output signal of the filter according to the transmission signal sent by the transmitting end of the full-duplex ISAC system and the weight coefficient of the filter;

具体的,滤波器输出信号的获取方法与常规LMS自适应滤波法中输出信号的获取方法相同,均是将全双工ISAC系统发送端所发送的发送信号进行共轭转置后,再与当前时刻M阶滤波器的权重系数相乘来得到。但在一具体实施例中,步骤S1之前,需对滤波器进行初始化,使权重系数为0,并为输入信号和误差信号赋初值,便于滤波器快速收敛,以更好的进行工作。Specifically, the acquisition method of the filter output signal is the same as the acquisition method of the output signal in the conventional LMS adaptive filtering method, which is to perform conjugate transposition on the transmission signal transmitted by the transmission end of the full-duplex ISAC system, and then combine it with the current It is obtained by multiplying the weight coefficients of the M-order filter at the moment. However, in a specific embodiment, before step S1, the filter needs to be initialized so that the weight coefficient is 0, and an initial value is assigned to the input signal and the error signal, so that the filter can converge quickly and work better.

S2、根据全双工ISAC系统中基带接收信号及输出信号进行误差计算,得到误差信号;S2. Perform error calculation according to the baseband receiving signal and output signal in the full-duplex ISAC system to obtain an error signal;

具体的,误差信号为全双工ISAC系统中基带接收信号与滤波器的输出信号之差,在一具体实施例中,基带接收信号为全双工ISAC系统中经过天线域自干扰抵消、模拟域自干扰抵消及ADC量化后的信号,基带接收信号的表达式为:Specifically, the error signal is the difference between the baseband received signal and the output signal of the filter in the full-duplex ISAC system. The signal after self-interference cancellation and ADC quantization, the expression of the baseband received signal is:

r(n)=sI(n)+d(n)+ε(n)r(n)=s I (n)+d(n)+ε(n)

式中,sI(n)为天线域自干扰抵消及模拟域自干扰抵消后的残余自干扰信号,d(n)为远端期望信号,ε(n)为加性噪声。where s I (n) is the residual self-interference signal after antenna domain self-interference cancellation and analog domain self-interference cancellation, d(n) is the desired signal at the far end, and ε(n) is additive noise.

根据图1及图2可知,基带接收信号为经过完成ADC量化后包括远端期望信号、部分的自干扰信号和噪声信号的r(n),该信号最终在基带解调前进行数字域的自干扰抵消。而误差信号作为LMS算法更新权系数向量的反馈信号,在自干扰抵消模型中,当自适应算法收敛后,该信号可表示为数字域经过自干扰抵消后剩余的信号,即远端期望信号;并根据本发明中改进的LMS算法(也可称为IVSSLMS算法),对滤波器中的权重系数进行更新,使最终输出的误差信号尽可能地逼近远端期望信号,从而达到自适应干扰抵消的目的。According to Figure 1 and Figure 2, it can be seen that the baseband received signal is r(n) including the remote desired signal, part of the self-interference signal and noise signal after the ADC quantization is completed, and the signal is finally self-assessed in the digital domain before baseband demodulation Interference cancellation. The error signal is used as the feedback signal for the LMS algorithm to update the weight coefficient vector. In the self-interference cancellation model, after the adaptive algorithm converges, the signal can be expressed as the remaining signal in the digital domain after self-interference cancellation, that is, the remote desired signal; And according to the improved LMS algorithm (also known as IVSSLMS algorithm) in the present invention, the weight coefficient in the filter is updated, so that the error signal of the final output is as close as possible to the far-end desired signal, thereby achieving the effect of adaptive interference cancellation Purpose.

S3、根据误差信号和改进的箕舌线函数,构造步长调整函数,并根据步长调整函数对权重系数进行迭代更新;S3. Construct a step size adjustment function according to the error signal and the improved skip line function, and iteratively update the weight coefficient according to the step size adjustment function;

在一具体实施例中,步骤S3具体包括:将误差函数作为自变量代入改进的箕舌线函数中,根据步长调整原则,对误差信号添加绝对值并引入调整参数,得到步长调整函数。本发明所提出的IVSSLMS算法通过改进的箕舌线函数与误差信号建立新的非线性关系,来构造步长调整函数,进而更新步长因子和权重系数,在复杂度和收敛性之间建立平衡。In a specific embodiment, step S3 specifically includes: substituting the error function as an argument into the improved skip line function, adding an absolute value to the error signal and introducing adjustment parameters according to the step size adjustment principle to obtain a step size adjustment function. The IVSSLMS algorithm proposed by the present invention establishes a new nonlinear relationship between the improved skip line function and the error signal to construct a step size adjustment function, and then update the step size factor and weight coefficient to establish a balance between complexity and convergence .

具体的,本申请的步长调整函数从最基础的箕舌线函数出发,对其进行改进,而改进后的箕舌线函数的表达式如下:Specifically, the step size adjustment function of the present application starts from the most basic skip line function and improves it, and the expression of the improved skip line function is as follows:

Figure BDA0004142398750000091
Figure BDA0004142398750000091

将误差函数作为自变量代入改进的箕舌线函数中,对代入改进的箕舌线函数中的误差函数引入第一调整参数,进而得到第一变步长函数,其公式表达式如下:Substituting the error function as an independent variable into the improved skip line function, introducing the first adjustment parameter to the error function substituted into the improved skip line function, and then obtaining the first variable step size function, the formula expression is as follows:

Figure BDA0004142398750000092
Figure BDA0004142398750000092

式中,μ1(n)为第一变步长函数,α为第一调整参数。In the formula, μ 1 (n) is the first variable step function, and α is the first adjustment parameter.

在改进的箕舌线函数中引入误差函数及第一调整参数而形成的第一变步长函数在LMS算法的基础上提高了收敛速度和跟踪速度。The first variable step size function formed by introducing the error function and the first adjustment parameter into the improved skip line function improves the convergence speed and tracking speed on the basis of the LMS algorithm.

在第一变步长函数中引入曲线参数来实现函数曲线可控,得到可控变步长函数,其公式表达式如下:Introduce the curve parameters in the first variable step function to realize the controllable function curve, and obtain the controllable variable step function, the formula expression is as follows:

Figure BDA0004142398750000101
Figure BDA0004142398750000101

式中,μk(n)为可控变步长函数,k为曲线参数。In the formula, μ k (n) is a controllable variable step function, and k is a curve parameter.

通过对可控变步长函数的分析可知,k值越小,曲线的最大收敛值越大,但是此时曲线下降速度越急;k值越大,曲线下降速度越缓,但是曲线的最大收敛值越小。因此,综合曲线收敛速度和稳态误差复杂度的综合分析,k值优选为2。Through the analysis of the controllable variable step size function, it can be known that the smaller the value of k, the greater the maximum convergence value of the curve, but the faster the curve descends at this time; the larger the k value, the slower the curve descends, but the maximum convergence of the curve The smaller the value. Therefore, for the comprehensive analysis of the convergence speed of the comprehensive curve and the complexity of the steady-state error, the value of k is preferably 2.

对确定k值的可控变步长函数引入控制函数取值范围的第二调整参数,并对包含误差函数的分母进行平方操作,得到第二变步长函数,其公式表达式如下:Introduce the second adjustment parameter of the value range of the control function to the controllable variable step size function that determines the value of k, and perform a square operation on the denominator including the error function to obtain the second variable step size function, and its formula expression is as follows:

Figure BDA0004142398750000102
Figure BDA0004142398750000102

式中,μ2(n)为第二变步长函数,β为第二调整参数,也是控制函数取值范围的系数。引入第二调整参数的第二变步长函数在曲线收敛速度和稳态误差中建立了平衡,并进一步降低了算法复杂度。In the formula, μ 2 (n) is the second variable step function, and β is the second adjustment parameter, which is also a coefficient controlling the value range of the function. The second variable step function introducing the second adjustment parameter establishes a balance between the curve convergence speed and the steady-state error, and further reduces the complexity of the algorithm.

为满足变步长LMS自适应滤波算法的步长调整原则,对第二变步长函数中的误差信号添加绝对值并引入第三调整参数,得到步长调整函数;其中,步长调整原则是:在初始收敛阶段或未知系统参数发生变化时,步长应比较大,以便有较快的收敛速度或对时变系统的跟踪速度;而在算法收敛后,都应保持很小的调整步长以达到很小的稳态失调。也就是说,变步长LMS算法在收敛初期有较大的步长因子,以得到较快的收敛速度;在要进入稳态阶段时,步长因子缓慢地减少至较小值,获得更好的稳态性能。因此,添加第三调整参数可动态地调整步长函数,使其满足步长调整原则,而步长因子为正数,对误差函数添加绝对值使步长调整函数成为偶函数来满足要求。In order to meet the step size adjustment principle of the variable step size LMS adaptive filtering algorithm, the absolute value is added to the error signal in the second variable step size function and the third adjustment parameter is introduced to obtain the step size adjustment function; among them, the step size adjustment principle is : In the initial convergence stage or when the unknown system parameters change, the step size should be relatively large in order to have a faster convergence speed or tracking speed of the time-varying system; after the algorithm converges, the adjustment step size should be kept small In order to achieve a small steady-state imbalance. That is to say, the variable step size LMS algorithm has a larger step size factor at the initial stage of convergence to obtain a faster convergence speed; when it is about to enter the steady state stage, the step size factor is slowly reduced to a smaller value to obtain a better steady-state performance. Therefore, adding a third adjustment parameter can dynamically adjust the step size function to meet the step size adjustment principle, and the step size factor is a positive number, and adding an absolute value to the error function makes the step size adjustment function an even function to meet the requirements.

在一具体实施例中,步长调整函数的公式表达式如下:In a specific embodiment, the formula expression of the step size adjustment function is as follows:

Figure BDA0004142398750000111
Figure BDA0004142398750000111

式中,μ(n)为步长调整函数,m为第三调整参数。In the formula, μ(n) is the step size adjustment function, and m is the third adjustment parameter.

本发明中的主要发明点即为对LMS算法进行改进,使传统的固定步长变为由改进的箕舌线函数与误差信号建立的新的非线性关系而形成的步长调整函数。本发明所提出的IVSSLMS算法在复杂度和收敛性之间建立了平衡,与经典的变步长LMS算法——SVSLMS和TCLMS算法相比,三种算法的步长更新变化的性能模拟结果如图4所示,由图4中不同函数的曲线变化可以看出,SVSLMS算法和TCLMS算法在误差因子接近0时产生突然的突变,而本发明提出的IVSSLMS算法克服了这一缺点,在误差信号接近于0时步长变化较为平缓,不仅在误差较大时提供较大的步长,而且提高了算法的收敛速度。The main invention point in the present invention is to improve the LMS algorithm, so that the traditional fixed step size becomes a step size adjustment function formed by the new nonlinear relationship established by the improved skip line function and the error signal. The IVSSLMS algorithm proposed by the present invention has established a balance between complexity and convergence. Compared with the classic variable step size LMS algorithm - SVSLMS and TCLMS algorithms, the performance simulation results of the step size update changes of the three algorithms are shown in the figure 4, as can be seen from the curve changes of different functions in Fig. 4, the SVSLMS algorithm and the TCLMS algorithm produce a sudden mutation when the error factor is close to 0, and the IVSSLMS algorithm proposed by the present invention overcomes this shortcoming, and when the error signal is close to When the step size is 0, the change is relatively gentle, which not only provides a larger step size when the error is large, but also improves the convergence speed of the algorithm.

在构建出步长调整函数后,通过以下公式来更新权重系数:After constructing the step size adjustment function, the weight coefficient is updated by the following formula:

w(n+1)=w(n)+2μ(n)e(n)s(n)w(n+1)=w(n)+2μ(n)e(n)s(n)

式中,w(n)为滤波器当前时刻的权重系数,w(n+1)为滤波器下一时刻的权重系数。In the formula, w(n) is the weight coefficient of the filter at the current moment, and w(n+1) is the weight coefficient of the filter at the next moment.

S4、根据更新后的权重系数对误差信号进行迭代更新,当更新后的误差信号的功率达到最小值时,输出更新后的误差信号及其对应的输出信号;S4. Iteratively updating the error signal according to the updated weight coefficient, and outputting the updated error signal and its corresponding output signal when the power of the updated error signal reaches a minimum value;

具体的,当权重系数更新后,进入IVSSLMS算法的循环步骤S1-S3,直至满足更新后的误差信号的功率达到最小值这一条件后,结束循环,并输出更新后的误差信号及其对应的输出信号。一般情况下,滤波器性能的评价标准通常是通过均方误差(Mean SquareError,MSE)来衡量,即均方误差越小性能越好,而均方误差公式表示如下:Specifically, after the weight coefficient is updated, enter the loop steps S1-S3 of the IVSSLMS algorithm until the power of the updated error signal reaches the minimum value, and then end the loop, and output the updated error signal and its corresponding output signal. In general, the evaluation standard of filter performance is usually measured by mean square error (Mean Square Error, MSE), that is, the smaller the mean square error, the better the performance, and the mean square error formula is expressed as follows:

MSE=E[e2(n)]MSE=E[e 2 (n)]

式中,E为求均值。LMS算法的收敛准则为使误差信号的功率最小,本发明同样使用该准则。也就是说,当自适应滤波器的权系数向量接近维纳解时,滤波器权系数向量可等效为最佳的自干扰信道估计,此时的自干扰抵消性能达到最佳。In the formula, E is the mean value. The convergence criterion of the LMS algorithm is to minimize the power of the error signal, which is also used in the present invention. That is to say, when the weight coefficient vector of the adaptive filter is close to the Wiener solution, the filter weight coefficient vector can be equivalent to the best self-interference channel estimation, and the self-interference cancellation performance at this time reaches the best.

不同变步长LMS算法下的均方误差性能模拟效果图如图5所示,在相同的输入信号和相同的仿真条件下,SVSLMS算法的收敛速度最慢,经过300次迭代开始收敛;TCLMS算法经过200次迭代开始收敛;而本发明提出的IVSSLMS算法在经过100次迭代后就开始收敛,表现出更快的收敛速度,且在第500次迭代发生时,信号发生相位变化,可见,本发明提出的改进算法相比较于其他两种算法保持较好的追踪能力。图6、7分别表示不同变步长LMS算法在低信噪比下的均方误差性能模拟效果图,从图中能看到本发明提出的IVSSLMS算法在低信噪比条件下仍可以保持良好的性能。The simulation effect diagram of the mean square error performance under different variable step size LMS algorithms is shown in Figure 5. Under the same input signal and the same simulation conditions, the convergence speed of the SVSLMS algorithm is the slowest, and it starts to converge after 300 iterations; the TCLMS algorithm Convergence begins after 200 iterations; and the IVSSLMS algorithm proposed by the present invention begins to converge after 100 iterations, showing a faster convergence speed, and when the 500th iteration occurs, the phase of the signal changes, it can be seen that the present invention Compared with the other two algorithms, the proposed improved algorithm maintains better tracking ability. Fig. 6, 7 respectively represent the mean square error performance simulation effect figure of different variable step length LMS algorithms under low signal-to-noise ratio, can see from the figure that the IVSSLMS algorithm proposed by the present invention can still keep good under low signal-to-noise ratio conditions performance.

在一具体实施例中,当干扰抵消比达到最大值时,结束迭代更新过程,并输出更新后的误差信号及其对应的输出信号;其中,干扰抵消比的计算公式如下:In a specific embodiment, when the interference cancellation ratio reaches the maximum value, the iterative update process is ended, and the updated error signal and its corresponding output signal are output; wherein, the calculation formula of the interference cancellation ratio is as follows:

Figure BDA0004142398750000121
Figure BDA0004142398750000121

式中,ICR为干扰抵消比,Pr为基带接收信号的功率,Pd为远端期望信号的功率,σε 2为加性噪声信号的功率,E{|e(∞)|2}为n趋于无穷大时的均方误差。In the formula, ICR is the interference cancellation ratio, P r is the power of baseband received signal, P d is the power of remote desired signal, σ ε 2 is the power of additive noise signal, E{|e(∞)| 2 } is The mean square error as n tends to infinity.

而数字域自干扰抵消后的干扰抵消比ICR,也是判断干扰抵消能力的重要指标,ICR的值越大,表明自干扰抵消能力越好,可以将干扰抵消比的值来作为循环过程是否结束的条件,当ICR达到最大值时,结束循环过程。The interference cancellation ratio ICR after self-interference cancellation in the digital domain is also an important indicator for judging the interference cancellation ability. The larger the ICR value, the better the self-interference cancellation ability. The value of the interference cancellation ratio can be used as the end of the cycle process. condition, when the ICR reaches the maximum value, the loop process ends.

不同变步长LMS算法下的ICR对比图如图8所示,由图8可知,将IVSSLMS算法应用到全双工ISAC数字域自干扰抵消中,该算法经过50次迭代就可以实现收敛,而TCLMS算法、SVSLMS算法则分别需要经过200次以及250次迭代才可以实现收敛过程,并且相比较于其他两种算法,该算法可以达到43dB的干扰抵消比,比其他两种算法可以提升2dB以上的抵消效果,不仅可以满足全双工ISAC系统的正常通信,还可实现较快的收敛速度且更高的干扰抵消比。The comparison chart of ICR under different variable step size LMS algorithms is shown in Figure 8. From Figure 8, it can be seen that the IVSSLMS algorithm is applied to full-duplex ISAC digital domain self-interference cancellation, and the algorithm can achieve convergence after 50 iterations, while The TCLMS algorithm and the SVSLMS algorithm need 200 and 250 iterations respectively to achieve the convergence process. Compared with the other two algorithms, this algorithm can achieve an interference cancellation ratio of 43dB, which is more than 2dB higher than the other two algorithms. The offset effect can not only meet the normal communication of the full-duplex ISAC system, but also achieve a faster convergence speed and a higher interference cancellation ratio.

因此,改进的LMS算法——IVSSLMS算法的流程总结如下表:Therefore, the process of the improved LMS algorithm - IVSSLMS algorithm is summarized in the following table:

Figure BDA0004142398750000131
Figure BDA0004142398750000131

α、β、m均大于0,且当e(n)→0时,

Figure BDA0004142398750000132
当e(n)→∞,/>
Figure BDA0004142398750000133
可以初步判断步长μ(n)的范围是(0,0.5β),β是调节μ(n)取值范围的系数,影响步长μ(n)变化的上限,α、m是调节μ(n)变化趋势的系数,影响步长的变化速率。α, β, m are all greater than 0, and when e(n)→0,
Figure BDA0004142398750000132
When e(n)→∞, />
Figure BDA0004142398750000133
It can be preliminarily judged that the range of the step size μ(n) is (0, 0.5β), β is the coefficient to adjust the value range of μ(n), and affects the upper limit of the change of the step size μ(n), α and m are the adjustments of μ( n) The coefficient of the change trend, which affects the rate of change of the step size.

本发明通过控制变量法来研究参数α、β、m对步长因子的影响,即通过固定两个参数变动另一个参数来确定其对步长因子波动的影响。当β=1,m=2,α分别取0.5,1,2,5时的步长因子与误差的函数模型曲线如图9所示,4条曲线均满足步长调整原则。当误差为4(此时误差比较大)时,图中大多数曲线都已经达到或者接近步长因子的最大值,说明此时的收敛速度非常快,有利于快速朝着误差减小的方向变化;在误差逐渐减小为零的过程中,步长因子也逐渐减小为零。但是由于稳态误差也与步长有关,在系统接近收敛时,快速变化的步长对稳态失调会产生较大的影响,这将可能导致较大的振荡。因此在系统接近收敛时,为保证算法性能,要求步长因子缓慢减小为零。所以,综合考虑以上因素,α值宜在1~4范围内选取,本发明优选α=2,此时算法在误差较大时,步长因子取值较大,保证了算法收敛速度较快,而在误差接近零时步长曲线也逐渐贴近零刻度线,并且变化幅度较为平缓,算法获得更高的稳定性。The present invention studies the influence of parameters α, β and m on the step factor by controlling the variable method, that is, by fixing two parameters and changing another parameter to determine its influence on the fluctuation of the step factor. When β=1, m=2, and α are 0.5, 1, 2, and 5 respectively, the function model curves of step factor and error are shown in Fig. 9, and all four curves satisfy the step size adjustment principle. When the error is 4 (the error is relatively large at this time), most of the curves in the figure have reached or approached the maximum value of the step size factor, indicating that the convergence speed at this time is very fast, which is conducive to quickly changing towards the direction of error reduction ; As the error gradually decreases to zero, the step factor also gradually decreases to zero. However, since the steady-state error is also related to the step size, when the system is close to convergence, the rapidly changing step size will have a greater impact on the steady-state misalignment, which may lead to larger oscillations. Therefore, when the system is close to convergence, in order to ensure the performance of the algorithm, the step size factor is required to be slowly reduced to zero. Therefore, considering the above factors, the α value should be selected within the range of 1 to 4, and the preferred α=2 of the present invention, when the algorithm has a large error, the value of the step factor is larger, which ensures that the algorithm convergence speed is faster, When the error is close to zero, the step length curve is gradually close to the zero scale line, and the change range is relatively gentle, and the algorithm obtains higher stability.

当α=2、β=1,m分别取0.5,1,2,5时的步长因子与误差的函数模型曲线如图10所示,当m值较小时,在误差e(n)接近为0时步长无法平滑下降,因此m应该大于1。为了保证步长能够快速跟踪误差同步变化,应该使m值大一点,但是当m过大时,将会使函数复杂度增加,计算量变大,同时还会使小误差范围内步长为0,因此综合考虑,m值优选为3。When α=2, β=1, and m take 0.5, 1, 2, and 5 respectively, the function model curve of the step factor and the error is shown in Figure 10. When the value of m is small, the error e(n) is close to When the step size is 0, the step size cannot decrease smoothly, so m should be greater than 1. In order to ensure that the step size can quickly track the synchronous change of the error, the value of m should be made larger, but when m is too large, the complexity of the function will increase, the amount of calculation will increase, and the step size will be 0 within the small error range. Therefore, considering all factors, the value of m is preferably 3.

当α=2,m=3,β分别取0.05,0.1,0.2,0.5时的步长因子与误差的函数模型曲线如图11所示,从迭代步长与参数β关系图可以看出,对于相同的初始误差,β取值越大,初始收敛阶段算法的收敛速度越快。若β取值过大,虽然收敛速度得到提高,但是算法收敛后e(n)对应的步长μ(n)值也较大,造成较大的稳态误差,也会导致算法发散。如果在实际应用中对收敛速度的需求较高,则可以选用较大的β值;如果对稳态误差的需求较高,可选取较小的β值。但应当注意的是,当β取值过小时,步长μ(n)迭代变化区间较小,此时变步长LMS算法也就近似退化为固定步长LMS算法,收敛速度不会有较大提高。因此,本发明优选设置α为2,m为3,β为变参数,但μ(n)与β成正比,可根据实际需求决定参数β的最优取值。When α=2, m=3, and β are respectively 0.05, 0.1, 0.2, and 0.5, the function model curves of the step size factor and the error are shown in Figure 11. It can be seen from the relationship between the iteration step size and the parameter β that for For the same initial error, the larger the value of β, the faster the convergence speed of the algorithm in the initial convergence stage. If the value of β is too large, although the convergence speed is improved, the value of the step size μ(n) corresponding to e(n) after the algorithm converges is also large, resulting in a large steady-state error and algorithm divergence. If the demand for convergence speed is high in practical applications, a larger β value can be selected; if the demand for steady-state error is high, a smaller β value can be selected. However, it should be noted that when the value of β is too small, the iterative change interval of the step size μ(n) is small, and the variable step size LMS algorithm will degenerate approximately into a fixed step size LMS algorithm at this time, and the convergence speed will not be greatly improved. improve. Therefore, the present invention preferably sets α to 2, m to 3, and β to variable parameters, but μ(n) is proportional to β, and the optimal value of parameter β can be determined according to actual needs.

本发明实施例中基于固定步长值固有缺陷无法合理协调收敛速度和稳态误差及现有改进的LMS算法中存在的问题,设计了一种基于改进的LMS算法的自干扰抵消方法,其实现了根据全双工ISAC系统发送端所发送的发送信号及滤波器的权重系数,计算滤波器的输出信号;根据全双工ISAC系统中基带接收信号及输出信号进行误差计算,得到误差信号;根据误差信号和改进的箕舌线函数,构造步长调整函数,并根据步长调整函数对权重系数进行迭代更新;根据更新后的权重系数对误差信号进行迭代更新,当更新后的误差信号的功率达到最小值时,输出更新后的误差信号及其对应的输出信号的技术方案;以使在全双工ISAC系统数字域自干扰抵消中的改进算法虽处于低信噪比状态下,仍能保持良好的性能,并在保持较快收敛性的同时还能保持较低的计算复杂度。In the embodiment of the present invention, a self-interference cancellation method based on the improved LMS algorithm is designed based on the inherent defects of the fixed step size value, which cannot reasonably coordinate the convergence speed and steady-state error, and the existing improved LMS algorithm. The output signal of the filter is calculated according to the transmission signal sent by the transmitter of the full-duplex ISAC system and the weight coefficient of the filter; the error is calculated according to the baseband receiving signal and output signal in the full-duplex ISAC system, and the error signal is obtained; according to The error signal and the improved skip line function are used to construct a step adjustment function, and the weight coefficient is iteratively updated according to the step adjustment function; the error signal is iteratively updated according to the updated weight coefficient, when the power of the updated error signal When the minimum value is reached, the technical scheme of outputting the updated error signal and its corresponding output signal; so that the improved algorithm in the self-interference cancellation in the digital domain of the full-duplex ISAC system can still maintain the low signal-to-noise ratio. Good performance and low computational complexity while maintaining fast convergence.

需要说明的是,虽然上述流程图中的各个步骤按照箭头的指示依次显示,但是这些步骤并不是必然按照箭头指示的顺序依次执行。除非本发明中有明确的说明,这些步骤的执行并没有严格的顺序限制,这些步骤可以以其它的顺序执行。It should be noted that although the various steps in the above flow chart are displayed sequentially according to the arrows, these steps are not necessarily executed sequentially in the order indicated by the arrows. Unless otherwise specified in the present invention, the execution of these steps is not strictly limited to the order, and these steps can be executed in other orders.

在另一实施例中,如图12所示,本发明第二方面提供一种基于改进的LMS算法的自干扰抵消系统,包括:In another embodiment, as shown in FIG. 12 , the second aspect of the present invention provides a self-interference cancellation system based on an improved LMS algorithm, including:

输出信号获取模块11,用于根据全双工ISAC系统发送端所发送的发送信号及滤波器的权重系数,计算滤波器的输出信号;The output signal acquisition module 11 is used to calculate the output signal of the filter according to the transmission signal sent by the full-duplex ISAC system transmission end and the weight coefficient of the filter;

误差信号获取模块12,用于根据全双工ISAC系统中基带接收信号及输出信号进行误差计算,得到误差信号;The error signal acquisition module 12 is used for performing error calculation according to the baseband receiving signal and the output signal in the full-duplex ISAC system to obtain the error signal;

权重系数更新模块13,用于根据误差信号和改进的箕舌线函数,构造步长调整函数,并根据步长调整函数对权重系数进行迭代更新;Weight coefficient update module 13, is used for according to error signal and improved skip line function, constructs step size adjustment function, and according to step size adjustment function weight coefficient is iteratively updated;

循环截止输出模块14,用于根据更新后的权重系数对误差信号进行迭代更新,当更新后的误差信号的功率达到最小值时,输出更新后的误差信号及其对应的输出信号。The loop cutoff output module 14 is configured to iteratively update the error signal according to the updated weight coefficient, and output the updated error signal and its corresponding output signal when the power of the updated error signal reaches a minimum value.

在一具体实施例中,权重系数更新模块13包括:In a specific embodiment, the weight coefficient updating module 13 includes:

步长函数构造模块,用于将误差函数作为自变量代入改进的箕舌线函数中,根据步长调整原则,对误差信号添加绝对值并引入步长调整参数,得到步长调整函数。The step size function construction module is used to substitute the error function as an independent variable into the improved skip line function, and according to the step size adjustment principle, add an absolute value to the error signal and introduce a step size adjustment parameter to obtain a step size adjustment function.

在一具体实施例中,所述系统还包括滤波器初始化模块10,用于对滤波器进行初始化,使权重系数为0。In a specific embodiment, the system further includes a filter initialization module 10, configured to initialize the filter so that the weight coefficient is 0.

需要说明的是,上述一种基于改进的LMS算法的自干扰抵消系统中的各个模块可全部或部分通过软件、硬件及其组合来实现。上述各模块可以硬件形式内嵌于或独立于计算机设备中的处理器中,也可以以软件形式存储于计算机设备中的存储器中,以便于处理器调用执行以上各个模块对应的操作。关于一种基于改进的LMS算法的自干扰抵消系统的具体限定参见上文中对于一种基于改进的LMS算法的自干扰抵消方法的限定,二者具有相同的功能和作用,在此不再赘述。It should be noted that each module in the above self-interference cancellation system based on the improved LMS algorithm can be fully or partially realized by software, hardware and a combination thereof. The above-mentioned modules can be embedded in or independent of the processor in the computer device in the form of hardware, and can also be stored in the memory of the computer device in the form of software, so that the processor can invoke and execute the corresponding operations of the above-mentioned modules. For the specific definition of a self-interference cancellation system based on the improved LMS algorithm, please refer to the above-mentioned definition of a self-interference cancellation method based on the improved LMS algorithm. The two have the same function and effect, and will not be repeated here.

本发明第三方面提供了一种电子设备,该电子设备包括:A third aspect of the present invention provides an electronic device, the electronic device comprising:

处理器、存储器和总线;Processor, memory and bus;

所述总线,用于连接所述处理器和所述存储器;The bus is used to connect the processor and the memory;

所述存储器,用于存储操作指令;The memory is used to store operation instructions;

所述处理器,用于通过调用所述操作指令,可执行指令使处理器执行如本发明的第一方面所示的一种基于改进的LMS算法的自干扰抵消方法对应的操作。The processor is configured to invoke the operation instruction, and the executable instruction causes the processor to execute the operation corresponding to the self-interference cancellation method based on the improved LMS algorithm as shown in the first aspect of the present invention.

在一个可选实施例中提供了一种电子设备,如图13所示,图13所示的电子设备5000包括:处理器5001和存储器5003。其中,处理器5001和存储器5003相连,如通过总线5002相连。可选地,电子设备5000还可以包括收发器5004。需要说明的是,实际应用中收发器5004不限于一个,该电子设备5000的结构并不构成对本发明实施例的限定。An optional embodiment provides an electronic device, as shown in FIG. 13 , the electronic device 5000 shown in FIG. 13 includes: a processor 5001 and a memory 5003 . Wherein, the processor 5001 is connected to the memory 5003 , such as through a bus 5002 . Optionally, the electronic device 5000 may further include a transceiver 5004 . It should be noted that in practical applications, the transceiver 5004 is not limited to one, and the structure of the electronic device 5000 does not limit the embodiment of the present invention.

处理器5001可以是CPU,通用处理器,DSP,ASIC,FPGA或者其他可编程逻辑器件、晶体管逻辑器件、硬件部件或者其任意组合。其可以实现或执行结合本发明公开内容所描述的各种示例性的逻辑方框,模块和电路。处理器5001也可以是实现计算功能的组合,例如包含一个或多个微处理器组合,DSP和微处理器的组合等。The processor 5001 may be a CPU, a general processor, DSP, ASIC, FPGA or other programmable logic devices, transistor logic devices, hardware components or any combination thereof. It can implement or execute the various illustrative logical blocks, modules and circuits described in connection with the present disclosure. The processor 5001 may also be a combination that implements computing functions, for example, a combination of one or more microprocessors, a combination of a DSP and a microprocessor, and the like.

总线5002可包括一通路,在上述组件之间传送信息。总线5002可以是PCI总线或EISA总线等。总线5002可以分为地址总线、数据总线、控制总线等。为便于表示,图13中仅用一条粗线表示,但并不表示仅有一根总线或一种类型的总线。Bus 5002 may include a path for communicating information between the components described above. The bus 5002 can be a PCI bus or an EISA bus, etc. The bus 5002 can be divided into address bus, data bus, control bus and so on. For ease of representation, only one thick line is used in FIG. 13 , but it does not mean that there is only one bus or one type of bus.

存储器5003可以是ROM或可存储静态信息和指令的其他类型的静态存储设备,RAM或者可存储信息和指令的其他类型的动态存储设备,也可以是EEPROM、CD-ROM或其他光盘存储、光碟存储(包括压缩光碟、激光碟、光碟、数字通用光碟、蓝光光碟等)、磁盘存储介质或者其他磁存储设备、或者能够用于携带或存储具有指令或数据结构形式的期望的程序代码并能够由计算机存取的任何其他介质,但不限于此。Memory 5003 can be ROM or other types of static storage devices that can store static information and instructions, RAM or other types of dynamic storage devices that can store information and instructions, and can also be EEPROM, CD-ROM or other optical disk storage, optical disk storage (including compact discs, laser discs, optical discs, digital versatile discs, Blu-ray discs, etc.), magnetic disk storage media or other magnetic storage devices, or can be used to carry or store desired program code in the form of instructions or data structures and can be programmed by a computer Any other medium accessed, but not limited to.

存储器5003用于存储执行本发明方案的应用程序代码,并由处理器5001来控制执行。处理器5001用于执行存储器5003中存储的应用程序代码,以实现前述任一方法实施例所示的内容。The memory 5003 is used to store application program codes for executing the solutions of the present invention, and the execution is controlled by the processor 5001 . The processor 5001 is configured to execute the application program code stored in the memory 5003, so as to realize the content shown in any one of the foregoing method embodiments.

其中,电子设备包括但不限于:移动电话、笔记本电脑、数字广播接收器、PDA(个人数字助理)、PAD(平板电脑)、PMP(便携式多媒体播放器)、车载终端(例如车载导航终端)等等的移动终端以及诸如数字TV、台式计算机等等的固定终端。Among them, electronic devices include but are not limited to: mobile phones, notebook computers, digital broadcast receivers, PDAs (personal digital assistants), PADs (tablet computers), PMPs (portable multimedia players), vehicle-mounted terminals (such as vehicle-mounted navigation terminals), etc. Mobile terminals such as digital TVs, desktop computers, etc. and fixed terminals.

本发明第四方面提供了一种计算机可读存储介质,计算机可读存储介质上存储有计算机程序,该程序被处理器执行时实现本发明第一方面所示的一种基于改进的LMS算法的自干扰抵消方法。The fourth aspect of the present invention provides a computer-readable storage medium, and a computer program is stored on the computer-readable storage medium. When the program is executed by a processor, an improved LMS algorithm based on the improved LMS algorithm shown in the first aspect of the present invention is implemented. Self-interference cancellation method.

本发明的又一实施例提供了一种计算机可读存储介质,该计算机可读存储介质上存储有计算机程序,当其在计算机上运行时,使得计算机可以执行前述方法实施例中相应内容。Another embodiment of the present invention provides a computer-readable storage medium, on which a computer program is stored, and when the computer program is run on a computer, the computer can execute the corresponding content in the foregoing method embodiments.

此外,本发明的实施例还提出一种计算机可读存储介质,其上存储有计算机程序,该程序被处理器执行时实现上述方法的步骤。In addition, an embodiment of the present invention also proposes a computer-readable storage medium on which a computer program is stored, and when the program is executed by a processor, the steps of the above method are implemented.

综上,为解决固定步长值固有缺陷无法合理协调收敛速度和稳态误差及现有改进的LMS算法中存在的问题,本发明提供一种基于改进的LMS算法的自干扰抵消方法、系统、设备和介质,通过改进的箕舌线函数与误差信号建立新的非线性关系来构造步长调整函数,进而更新权重系数的步长因子,可在算法收敛过程中动态地改变步长因子的大小,使得改进的LMS算法在保持较快收敛性的同时还能保持较低的计算复杂度,且在低信噪比条件下,仍能保持良好的性能,并在全双工通感一体化数字域自干扰抵消中,在数字域的ICR高于其他比较的算法且具有较快的收敛速度。In summary, in order to solve the inherent defects of the fixed step value that cannot reasonably coordinate the convergence speed and steady-state error and the problems existing in the existing improved LMS algorithm, the present invention provides a self-interference cancellation method, system, For equipment and media, the step size adjustment function is constructed by establishing a new nonlinear relationship between the improved skip line function and the error signal, and then the step size factor of the weight coefficient is updated, and the size of the step size factor can be dynamically changed during the algorithm convergence process , so that the improved LMS algorithm can maintain a low computational complexity while maintaining fast convergence, and can still maintain good performance under the condition of low signal-to-noise ratio. In domain self-interference cancellation, the ICR in the digital domain is higher than other compared algorithms and has a faster convergence speed.

本说明书中的各个实施例均采用递进的方式描述,各个实施例直接相同或相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于系统实施例而言,由于其基本相似于方法实施例,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。需要说明的是,上述实施例的各技术特征可以进行任意的组合,为使描述简洁,未对上述实施例中的各个技术特征所有可能的组合都进行描述,然而,只要这些技术特征的组合不存在矛盾,都应当认为是本说明书记载的范围。Each embodiment in this specification is described in a progressive manner, and the same or similar parts of each embodiment can be directly referred to each other, and each embodiment focuses on the differences from other embodiments. In particular, for the system embodiment, since it is basically similar to the method embodiment, the description is relatively simple, and for relevant parts, refer to part of the description of the method embodiment. It should be noted that the various technical features of the above-mentioned embodiments can be combined arbitrarily. For the sake of concise description, all possible combinations of the various technical features in the above-mentioned embodiments are not described. Where there is a contradiction, all should be deemed to be within the scope of this specification.

以上所述实施例仅表达了本发明的几种优选实施方式,其描述较为具体和详细,但并不能因此而理解为对发明专利范围的限制。应当指出的是,对于本技术领域的普通技术人员来说,在不脱离本发明技术原理的前提下,还可以做出若干改进和替换,这些改进和替换也应视为本发明的保护范围。因此,本发明专利的保护范围应以所述权利要求的保护范围为准。The above-mentioned examples only express several preferred implementation modes of the present invention, and the description thereof is relatively specific and detailed, but should not be construed as limiting the patent scope of the invention. It should be noted that those skilled in the art can make some improvements and substitutions without departing from the technical principle of the present invention, and these improvements and substitutions should also be regarded as the protection scope of the present invention. Therefore, the protection scope of the patent for the present invention should be based on the protection scope of the claims.

Claims (10)

1. A self-interference cancellation method based on an improved LMS algorithm, comprising the steps of:
calculating an output signal of a filter according to a transmission signal transmitted by a full duplex ISAC system transmitting end and a weight coefficient of the filter;
performing error calculation according to a baseband receiving signal and the output signal in the full duplex ISAC system to obtain an error signal;
constructing a step length adjusting function according to the error signal and the improved skip tongue function, and iteratively updating the weight coefficient according to the step length adjusting function;
and iteratively updating the error signal according to the updated weight coefficient, and outputting the updated error signal and the corresponding output signal when the power of the updated error signal reaches the minimum value.
2. The method for self-interference cancellation based on improved LMS algorithm as claimed in claim 1, wherein before said calculating the output signal of the filter according to the transmission signal transmitted from the transmitting end of the full duplex ISAC system and the weight coefficient of the filter, comprising:
initializing the filter to make the weight coefficient be 0.
3. A method of self-interference cancellation based on an improved LMS algorithm as claimed in claim 1, wherein said constructing a step size adjustment function based on said error signal and an improved skip tongue function comprises:
substituting the error function as an independent variable into an improved skip tongue function, adding an absolute value to the error signal according to a step length adjustment principle, and introducing an adjustment parameter to obtain the step length adjustment function.
4. A method of self-interference cancellation based on an improved LMS algorithm as claimed in claim 3, wherein said step size adjustment function is expressed as:
Figure FDA0004142398740000011
wherein μ (n) is a step size adjustment function, α is a first adjustment parameter, β is a second adjustment parameter, m is a third adjustment parameter, α, β, m are all greater than 0, and e (n) is an error function.
5. The method for self-interference cancellation based on improved LMS algorithm as claimed in claim 1, wherein said baseband received signal is a signal after antenna domain self-interference cancellation, analog domain self-interference cancellation and ADC quantization in said full duplex ISAC system, and the expression of said baseband received signal is:
r(n)=s I (n)+d(n)+ε(n)
where r (n) is the baseband received signal, s I (n) is the residual self-interference signal after the self-interference cancellation of the antenna domain and the self-interference cancellation of the analog domain, d (n) is the far-end expected signal, and epsilon (n) is the additive noise.
6. A method of self-interference cancellation based on an improved LMS algorithm as claimed in claim 1, further comprising: when the interference cancellation ratio reaches the maximum value, ending the iterative updating process, and outputting an updated error signal and a corresponding output signal thereof; wherein, the calculation formula of the interference cancellation ratio is as follows:
Figure FDA0004142398740000021
wherein ICR is interference cancellation ratio, P r For the power of the baseband received signal, P d Work for remote desired signalRate, sigma ε 2 E { E (++E) is the power of the additive noise signal 2 And n is the mean square error when n approaches infinity.
7. A self-interference cancellation system based on an improved LMS algorithm, comprising:
the output signal acquisition module is used for calculating the output signal of the filter according to the transmission signal transmitted by the transmitting end of the full-duplex ISAC system and the weight coefficient of the filter;
the error signal acquisition module is used for carrying out error calculation according to a baseband receiving signal and the output signal in the full-duplex ISAC system to obtain an error signal;
the weight coefficient updating module is used for constructing a step length adjusting function according to the error signal and the improved skip tongue function and iteratively updating the weight coefficient according to the step length adjusting function;
and the loop cut-off output module is used for carrying out iterative updating on the error signal according to the updated weight coefficient, and outputting the updated error signal and the corresponding output signal when the power of the updated error signal reaches the minimum value.
8. A self-interference cancellation system based on an improved LMS algorithm as claimed in claim 7, wherein said weight coefficient update module comprises;
the step length function construction module is used for substituting the error function as an independent variable into the improved skip tongue function, adding an absolute value to the error signal according to a step length adjustment principle, and introducing an adjustment parameter to obtain the step length adjustment function.
9. An electronic device comprising a processor, a memory, and a computer program stored in the memory and configured to be executed by the processor, the processor implementing the improved LMS algorithm-based self-interference cancellation method of any one of claims 1 to 6 when the computer program is executed.
10. A computer readable storage medium, characterized in that the computer readable storage medium comprises a stored computer program, wherein the computer program when run controls a device in which the computer readable storage medium is located to perform the self-interference cancellation method based on the modified LMS algorithm as claimed in any one of claims 1 to 6.
CN202310293398.1A 2023-03-23 2023-03-23 Self-interference cancellation method, system, device and medium based on improved LMS algorithm Pending CN116405129A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310293398.1A CN116405129A (en) 2023-03-23 2023-03-23 Self-interference cancellation method, system, device and medium based on improved LMS algorithm

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310293398.1A CN116405129A (en) 2023-03-23 2023-03-23 Self-interference cancellation method, system, device and medium based on improved LMS algorithm

Publications (1)

Publication Number Publication Date
CN116405129A true CN116405129A (en) 2023-07-07

Family

ID=87015202

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310293398.1A Pending CN116405129A (en) 2023-03-23 2023-03-23 Self-interference cancellation method, system, device and medium based on improved LMS algorithm

Country Status (1)

Country Link
CN (1) CN116405129A (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN118449672A (en) * 2024-06-21 2024-08-06 无锡瀚诺光电科技有限公司 A signal processing method for single-station full-duplex communication and perception integration
CN118473483A (en) * 2024-05-11 2024-08-09 深圳市中承科技有限公司 Self-adaptive beam forming method, system and device
CN119901512A (en) * 2025-04-02 2025-04-29 拾音汽车科技(上海)有限公司 A buffer block test bench compensation noise reduction method and system

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN118473483A (en) * 2024-05-11 2024-08-09 深圳市中承科技有限公司 Self-adaptive beam forming method, system and device
CN118449672A (en) * 2024-06-21 2024-08-06 无锡瀚诺光电科技有限公司 A signal processing method for single-station full-duplex communication and perception integration
CN118449672B (en) * 2024-06-21 2024-12-06 无锡学院 Single-station full duplex communication perception-oriented integrated signal processing method
CN119901512A (en) * 2025-04-02 2025-04-29 拾音汽车科技(上海)有限公司 A buffer block test bench compensation noise reduction method and system

Similar Documents

Publication Publication Date Title
CN116405129A (en) Self-interference cancellation method, system, device and medium based on improved LMS algorithm
CN110492868B (en) A Multi-parameter Variable Step Size LMS Adaptive Filtering Method
CN109741727A (en) Active noise reduction earphone, noise reduction method and storage medium based on active noise control algorithm
CN103716013B (en) Variable element ratio sef-adapting filter
CN107425929B (en) Unaided Data Equalization Method for Fading Channels under Alpha Stable Distributed Noise
Resende et al. A variable step-size NLMS algorithm with adaptive coefficient vector reusing
CN111884625B (en) Variable step-length sparse augmentation complex number self-adaptive filter
KR102190833B1 (en) Echo suppression
CN107342751A (en) A kind of variable step SA adaptive filter algorithms based on cross-correlation entropy
CN110535936A (en) A kind of energy efficient mist computation migration method based on deep learning
CN105027540A (en) Echo suppression
Li et al. Variable step-size LMS algorithm based on hyperbolic tangent function
CN115360999A (en) Adaptive variable step size LMS filter based on hyperbolic tangent function and its filtering method and computer equipment
CN112003588B (en) Adaptive signal filtering method based on polymorphic variable step length normalized mean square
WO2025020825A1 (en) Wind turbine-participating power system, and wind turbine torque oscillation suppression method and apparatus
CN118101391A (en) Interference cancellation method, device, terminal and storage medium
CN117077740B (en) Model quantification methods and equipment
Luo et al. Affine Projection Algorithms with Novel Schemes of Variable Projection Order
CN116165898A (en) Design and configuration method and system for robust power system stabilizer based on wolf algorithm
WO2023274285A1 (en) Filtering method and filtering control apparatus for load cell
CN113452350B (en) Variable step block sparse affine projection self-adaptive filter
CN113382081B (en) Time delay estimation adjusting method, device, equipment and storage medium
CN117240678A (en) Artificial intelligence model training method, system and related equipment
Shu et al. Set-membership LMS adaptive algorithms based on an error-estimation time-varying bound method
Apolinário Jr et al. Introduction to adaptive filters

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