CN110034747B - Robust complex scale symbol adaptive filter - Google Patents
Robust complex scale symbol adaptive filter Download PDFInfo
- Publication number
- CN110034747B CN110034747B CN201910259804.6A CN201910259804A CN110034747B CN 110034747 B CN110034747 B CN 110034747B CN 201910259804 A CN201910259804 A CN 201910259804A CN 110034747 B CN110034747 B CN 110034747B
- Authority
- CN
- China
- Prior art keywords
- adaptive filter
- filter
- complex
- vector
- matrix
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F7/00—Methods or arrangements for processing data by operating upon the order or content of the data handled
- G06F7/38—Methods or arrangements for performing computations using exclusively denominational number representation, e.g. using binary, ternary, decimal representation
- G06F7/48—Methods or arrangements for performing computations using exclusively denominational number representation, e.g. using binary, ternary, decimal representation using non-contact-making devices, e.g. tube, solid state device; using unspecified devices
- G06F7/4806—Computations with complex numbers
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
- H03H17/00—Networks using digital techniques
- H03H17/02—Frequency selective networks
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Mathematical Optimization (AREA)
- Computing Systems (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Filters That Use Time-Delay Elements (AREA)
- Soundproofing, Sound Blocking, And Sound Damping (AREA)
Abstract
The invention discloses a robust complex scale symbol adaptive filter, and belongs to the field of digital filter design. The filter adopts a symbol operation and proportion self-adaptive method to enhance the anti-impulse noise interference capability and accelerate the convergence speed. The filter can be applied to instruments and communication equipment which are interfered by noise in an electronic and communication system and have sparse estimation systems.
Description
Technical Field
The invention discloses a robust complex scale symbol self-adaptive filter, and belongs to the field of digital filter design.
Background
Complex adaptive filters have been widely used in recent years because the amount of information that complex numbers can carry exceeds real numbers. Commonly used complex adaptive filters are complex least mean square filter (abbreviated as CLMS filter), complex normalized least mean square filter (abbreviated as CNLMS filter), etc. since these two filters are low in computational complexity and easy to implement. In some special circumstances, the output signal of an unknown system may be disturbed by impulse noise. Because the two filters adopt the minimum mean square error criterion, the pulse noise resistance is poor, namely the robustness is not strong.
System identification is an important branch of adaptive signal processing, and many problems of traditional adaptive channel equalization, adaptive noise cancellation, adaptive echo cancellation, active noise control and the like can be summarized as system identification problems. In some applications, most elements in the coefficient vector of an unknown system are zero or near zero, and such systems are generally referred to as sparse systems. Sparse system identification problems often involve theoretical and engineering practices, such as zero attraction theory and the application of proportional adaptive strategies in satellite transmission channels as well as echo cancellation channels. Duttweiler D L proposes a proportional normalized least mean square filter (abbreviated as PNLMS filter) [ Proportate normalized least mean square adaptation in echo cancellers, IEEE Transactions on spech and Audio Processing,2000,8 (5): 508-518]. However, the filter can only be used in the real number domain, and cannot be used for complex sparse system identification.
Therefore, to estimate a complex sparse system in an impulse interference environment, a robust sparse adaptive filter needs to be designed.
Disclosure of Invention
In order to solve the above problems, the present invention discloses a robust complex scale symbol adaptive filter (abbreviated as RCPSA filter). The RCPSA filter adopts an anti-impulse interference method and a proportional adaptive strategy to update the weight vector of the filter. The updating weight vector of the RCPSA filter mainly comprises the following steps:
1) By input signal x at n times n And a desired signal d n Calculating an error signal e n I.e. byWherein x is n =[x n ,x n-1 ,...x n-L+1 ] T Representing the first L samples x of the input signal n ,x n-1 ,...x n-L+1 Is formed of an input vector, w n =[w 0,n ,w 2,n ,...w L-1,n ] T The weight vector is formed by L weights of the adaptive filter, T represents transposition operation, and H represents conjugate transposition operation.
2) From an input vector x n Conjugation of error signalModulus e (n) of the error signal and a complex scaling matrix G n-1 Calculating the intermediate vector updated by the weight vectorWherein: delta is a small normal number to avoid zero denominator in the iteration, which in a typical environment ranges from 10 -4 ,10 -1 ](ii) a The complex scale matrix is denoted G n-1 =diag{g 0,n-1 ,g 1,n-1 ...,g L-1,n-1 Denotes the element g 0,n-1 ,g 1,n-1 ...,g L-1,n-1 Write into a diagonal matrix, taking the diagonal elements of the matrix as||w n-1 || 1 Representing the coefficient vector l to the adaptive filter 1 Norm, parameter alpha is in the range of-1, 1) for balancing between sparse and non-sparse systems, and parameter epsilon is a small normal number for avoidingThe denominator is zero in the iteration-free process, and the value range is [10 ] -6 ,10 -3 ]。
3) By iterative formula w n =w n-1 +μG n-1 f n Updating weight vector of the adaptive filter, wherein mu represents the step length of the filter and the value range is [0.0001,0.05 ]]。
In the above step 2), if the matrix G is n-1 Taking the identity matrix as the basis matrix, the RCPSA filter degenerates to a robust complex symbol adaptive filter (abbreviated as RCSA filter).
Advantageous effects
Compared with the existing symbol adaptive filter scheme, the adaptive filter disclosed by the invention has stronger robustness and can accelerate the convergence speed of the identification of the complex sparse system.
Drawings
The invention is further described with reference to the following figures and examples:
FIG. 1 is a block diagram of a robust complex scale symbol adaptive filter architecture;
fig. 2 is a comparison of mean square deviation curves of the adaptive filter under the action of a round gaussian white signal.
FIG. 3 is a comparison of the mean square deviation curves of the adaptive filter under the action of the round Gaussian colored signal.
Detailed Description
Examples
The principle of the technical scheme of the invention is as follows:
firstly, expanding a cost function of the established normalized symbol adaptive filter into a complex domain and correcting, then solving the gradient of the corrected cost function, taking the gradient as an intermediate vector for reducing the pulse interference sensitivity, and then introducing a complex proportion matrix into a fastest descent method to obtain a weight vector updating formula of the adaptive filter disclosed by the application.
The embodiment of the invention adopts a computer experiment method to verify the performance of the RCPSA filter. The RCPSA filter disclosed by the invention is used for identifying an unknown complex sparse system in an impulse noise interference environment, and the performance of the RCPSA filter is compared with the performance of a complex normalized least mean square adaptive filter (abbreviated as a CNLMS filter) and an RCSA filter.
In order to make the experimental result more general, the round Gaussian white signal and the round Gaussian colored signal are selected to be respectively subjected to the experiment, and the variance of the round Gaussian white signalVariance of round Gaussian colored signalWherein, the round Gaussian colored signal is generated by an AR (1) model, and the transfer function is F (z) = 1/(1-0.5 z) -1 ). In the experiment, mean Square Deviation (MSD) was used as a measure of algorithm performance, which is defined asIn dB, where wo represents the unknown system to be identified,express vector w n -w o L of 2 The square of the norm, log denotes the common logarithm. The simulated MSD curves in the figure are obtained by averaging 100 independent experiments.
Pulse signal v used in the experiment n Comprises a zero mean, a variance ofWhite gaussian noise gamma of n And an impulse noise z n I.e. v n =γ n +z n . Impulse noise z n Generated by the Bernoulli Gaussian process, i.e. z n =ξ n ψ n Wherein xi is n Is a Bernoulli process, and P [ delta ] n =1]=0.01,P[δ n =0]=0.99,ψ n White gaussian noise with zero mean.
The method for identifying the unknown complex sparse system by the RCPSA adaptive filter comprises the following steps:
2) From an input vector x n Conjugation of error signalsModulus e (n) of the error signal and a complex scaling matrix G n-1 Calculating the intermediate vector updated by the weight vectorIn this example, δ is taken to be 0.025. In other embodiments, the delta value lies in the interval [10 ] -4 ,10 -1 ]Internal; the complex scale matrix is denoted G n-1 =diag{g 0,n-1 ,g 1,n-1 ...,g L-1,n-1 Indicates to put the element g 0,n-1 ,g 1,n-1 ...,g L-1,n-1 Write into a diagonal matrix, taking the diagonal elements of the matrix asε is used to avoid zero denominator in the iteration, and ε is taken to be 0.0001 in this example. In other embodiments, the value of ε is located at [10 ] -6 ,10 -3 ]And (4) the following steps.
3) By iterative formula w n =w n-1 +0.01G n-1 f n And updating the formula.
As can be seen from fig. 2 and 3, the RCPSA filter disclosed by the present invention has good impulse noise resistance under both sparse conditions and non-sparse conditions, and has the fastest convergence rate when estimating a complex sparse system compared with an RCSA that does not adopt a scaling strategy.
The above embodiments are merely illustrative of the technical ideas and features of the present invention, and the purpose thereof is to enable those skilled in the art to understand the contents of the present invention and implement the present invention, and not to limit the protection scope of the present invention. All equivalent changes and modifications made according to the spirit of the present invention should be covered in the protection scope of the present invention.
Claims (1)
1. A robust complex scale sign adaptive filter characterized by: in the complex domain, the adaptive filter updates its weight vector based on a proportional-adaptive strategy;
the weight vector updating comprises the following steps:
1) By input signal x at n times n And a desired signal d n Calculating an error signal e n I.e. byWherein x is n =[x n ,x n-1 ,...x n-L+1 ] T Representing the first L samples x of the input signal n ,x n-1 ,...x n-L+1 Formed of input vectors, w n =[w 0,n ,w 2,n ,...w L-1,n ] T The weight vector is formed by L weights of the adaptive filter, T represents transposition operation, and H represents conjugate transposition operation;
2) From an input vector x n Conjugation of error signalModulo e (n) of the error signal and a complex scaling matrix G n-1 Calculating the intermediate vector updated by the weight vectorWherein: delta is a small normal number, which is used to avoid the occurrence of zero denominator in the iteration, and the value range is [10 ] -4 ,10 -1 ](ii) a The complex scale matrix is denoted G n-1 =diag{g 0,n-1 ,g 1,n-1 ...,g L-1,n-1 Indicates to put the element g 0,n-1 ,g 1,n-1 ...,g L-1,n-1 Write into a diagonal matrix, taking the diagonal elements of the matrix as||w n-1 || 1 Expressing the coefficient vector l to the adaptive filter 1 Norm with parameter alpha in the range of-1, 1) for balancing between sparse and non-sparse systems, and parameter epsilon, a very small normal number for avoiding denominator zero in the iterative process, with a value in the range of [10 ] -6 ,10 -3 ];
3) By iterative formula w n =w n-1 +μG n-1 f n Updating weight vector of the adaptive filter, wherein mu represents the step length of the filter and the value range is [0.0001,0.05 ]]。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910259804.6A CN110034747B (en) | 2019-04-02 | 2019-04-02 | Robust complex scale symbol adaptive filter |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910259804.6A CN110034747B (en) | 2019-04-02 | 2019-04-02 | Robust complex scale symbol adaptive filter |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110034747A CN110034747A (en) | 2019-07-19 |
CN110034747B true CN110034747B (en) | 2023-02-03 |
Family
ID=67237214
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910259804.6A Active CN110034747B (en) | 2019-04-02 | 2019-04-02 | Robust complex scale symbol adaptive filter |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110034747B (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113225045B (en) * | 2021-03-25 | 2023-06-23 | 苏州大学 | Sparse-facilitated affine projection adaptive filter with low computational complexity |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104158512A (en) * | 2014-08-08 | 2014-11-19 | 西南交通大学 | self-adaptive sparse system identification method based on impact-interference-resistance of independent activity factor |
CN106059531A (en) * | 2016-05-24 | 2016-10-26 | 苏州大学 | Nonnegative adaptive filter |
CN107592096A (en) * | 2017-09-29 | 2018-01-16 | 苏州大学 | A kind of robust deviation compensation sef-adapting filter and its filtering method |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8995594B2 (en) * | 2012-08-22 | 2015-03-31 | Intel Corporation | Baseband cancellation of platform radio interference |
-
2019
- 2019-04-02 CN CN201910259804.6A patent/CN110034747B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104158512A (en) * | 2014-08-08 | 2014-11-19 | 西南交通大学 | self-adaptive sparse system identification method based on impact-interference-resistance of independent activity factor |
CN106059531A (en) * | 2016-05-24 | 2016-10-26 | 苏州大学 | Nonnegative adaptive filter |
CN107592096A (en) * | 2017-09-29 | 2018-01-16 | 苏州大学 | A kind of robust deviation compensation sef-adapting filter and its filtering method |
Also Published As
Publication number | Publication date |
---|---|
CN110034747A (en) | 2019-07-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107592096B (en) | Filtering method of robust deviation compensation adaptive filter | |
CN111884625B (en) | Variable step-length sparse augmentation complex number self-adaptive filter | |
CN111277244A (en) | Variable-step zero-attraction normalization double-symbol adaptive filter | |
CN104158512B (en) | The adaptive sparse system identifying method that a kind of shock resistance based on the independent activities factor is disturbed | |
CN114420147A (en) | Echo eliminator based on proportion total adaptive method | |
CN114172490B (en) | Ekblom norm-based robust self-adaptive noise elimination method | |
CN110034747B (en) | Robust complex scale symbol adaptive filter | |
CN117896468B (en) | Deviation compensation echo cancellation method and system for telephone communication | |
CN109818888B (en) | Group sparse underwater acoustic channel estimation method in pulse interference environment | |
CN109089004B (en) | Collective member self-adaptive echo cancellation method based on correlation entropy induction | |
CN110190831B (en) | Mixed norm non-negative adaptive filter | |
CN109151237B (en) | Affine projection self-adaptive echo cancellation method based on zero attraction | |
CN112462352B (en) | Line spectrum enhancement method suitable for low signal-to-noise ratio condition | |
CN107800403B (en) | Robust spline self-adaptive filter | |
CN112886947A (en) | Variable-step robust affine projection adaptive filter | |
CN110190832B (en) | Regularization parameter multi-task adaptive filter network | |
CN113452350B (en) | Variable step block sparse affine projection self-adaptive filter | |
CN112003588A (en) | Adaptive signal filtering method based on polymorphic variable step size normalized mean square | |
CN111181531B (en) | Variable regularization deviation compensation symbol sub-band adaptive filter | |
CN113541650B (en) | Sparse linear constraint recursive maximum correlation entropy adaptive filter | |
CN112929006A (en) | Variable step size selection updating kernel minimum mean square adaptive filter | |
Albu | A proportionate affine projection algorithm using fast recursive filtering and dichotomous coordinate descent iterations | |
CN110492869A (en) | A kind of improved segmentation area block LMS adaptive filter algorithm | |
CN117040489B (en) | Spline self-adaptive filter with sparse constraint | |
Xu et al. | A simplified RLS algorithm and its application in acoustic echo cancellation |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |