CN112583381A - Self-adaptive filtering method based on deviation compensation auxiliary variable - Google Patents

Self-adaptive filtering method based on deviation compensation auxiliary variable Download PDF

Info

Publication number
CN112583381A
CN112583381A CN202011481813.9A CN202011481813A CN112583381A CN 112583381 A CN112583381 A CN 112583381A CN 202011481813 A CN202011481813 A CN 202011481813A CN 112583381 A CN112583381 A CN 112583381A
Authority
CN
China
Prior art keywords
time
noise
estimated
estimation
input
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202011481813.9A
Other languages
Chinese (zh)
Other versions
CN112583381B (en
Inventor
贾丽娟
李妍
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Institute of Technology BIT
Original Assignee
Beijing Institute of Technology BIT
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beijing Institute of Technology BIT filed Critical Beijing Institute of Technology BIT
Priority to CN202011481813.9A priority Critical patent/CN112583381B/en
Publication of CN112583381A publication Critical patent/CN112583381A/en
Application granted granted Critical
Publication of CN112583381B publication Critical patent/CN112583381B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03HIMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
    • H03H21/00Adaptive networks
    • H03H21/0012Digital adaptive filters
    • H03H21/0043Adaptive algorithms

Landscapes

  • Image Processing (AREA)
  • Filters That Use Time-Delay Elements (AREA)

Abstract

The invention relates to a self-adaptive filtering method based on deviation compensation auxiliary variables, and belongs to the technical field of signal estimation and digital filters. The method comprises the following steps: presetting iteration times M and constructing a variable error-containing model; constructing a class auxiliary variable strongly related to an input signal vector; calculating an estimated value and an estimated deviation of a similar auxiliary variable method under the condition of colored output noise; calculating an estimated value of unknown input noise variance under the colored output noise condition; and compensating the deviation caused by the noise by using a deviation compensation principle to obtain unbiased estimation of unknown system parameters. The method can stably work when the input signal is a white Gaussian process or a colored Gaussian process and the output noise signal is colored noise; the influence of the output noise variance is eliminated by introducing the class auxiliary variable, only the input noise variance is estimated, and the complexity of the method is reduced; the method for estimating the input noise variance in real time is provided, and the unknown input noise variance can be accurately estimated.

Description

Self-adaptive filtering method based on deviation compensation auxiliary variable
Technical Field
The invention relates to a self-adaptive filtering method based on deviation compensation auxiliary variables, and belongs to the technical field of signal estimation and digital filters.
Background
The traditional RLS adaptive filter has wide application in the field of adaptive filtering, such as adaptive equalization, echo cancellation, antenna array beam forming, parameter estimation, noise cancellation, spectrum estimation, etc. in the field of communication.
Convergence speed, steady state imbalance and robustness are three important performance indicators of the adaptive filter. The convergence rate determines the time taken for the adaptive filter to approach an unknown system, the steady state maladjustment height determines the estimation precision of the proposed method for the unknown system, and the robustness determines the application range and the effectiveness of the proposed method, and the three indexes simultaneously influence the quality of signal processing.
When the output end of the adaptive filter is interfered by colored noise, the traditional RLS method cannot realize a good filtering effect, and the problem that the least square estimation value is biased when the output noise is interfered by the colored noise under some conditions is solved by introducing auxiliary variables. However, in a common situation, when the input signal is a white gaussian process, the auxiliary variable method is unstable, and cannot accurately estimate unknown parameters of the system, and the application range needs to be expanded.
In some applications, the adaptive filter is disturbed by input noise in addition to output noise. If the parameters of an unknown system are estimated by adopting the traditional RLS adaptive filtering method, the estimation deviation introduced by noise exists, and a more accurate estimation result cannot be obtained.
Wherein RLS is Recursive Least square;
disclosure of Invention
The invention aims to overcome the defect that the conventional deviation compensation RLS method cannot normally work under the condition of colored output noise, provides a self-adaptive filtering method based on deviation compensation auxiliary variables, and realizes unbiased estimation of an unknown system under the condition that additive noise interference of unknown variance exists at the input end and the output end.
In order to achieve the technical purpose, the invention adopts the following technical scheme:
the adaptive filtering method based on the deviation compensation auxiliary variable comprises the following steps:
step A, presetting iteration times M, and constructing a variable error-containing model;
wherein, the variable error-containing model, namely the EIV-IIR filter model, can be expressed as:
Figure BDA0002837805710000021
wherein, EIV is Errors in Variables, meaning that the Variables contain Errors; IIR is Infine Impulse Response which is Infinite Impulse Response; y (i) is the noisy output signal at time i, piAn input signal vector at the ith moment of the EIV-IIR filter model is represented as h, an L-order unknown system parameter vector to be estimated is represented as v (i), composite noise at the ith moment is represented as v (i), superscript T represents transposition operation, and L is the order of the EIV-IIR filter model;
vector piAnd v (i) in particular:
pi=[y(i-1)y(i-2)…y(i-L)x(i-1)x(i-2)…x(i-L)]T (2)
Figure BDA0002837805710000022
ni=[e(i-1)e(i-2)…e(i-L)n(i-1)n(i-2)…n(i-L)]T (4)
wherein y (i-1) y (i-2) … y (i-L) is the delay of the noisy output signal at time i-1 and at an interval of 1 up to time i-L, respectively; x (i-1) x (i-2) … x (i-L) are the delay of the noisy input signal at time i-1 and at an interval of 1 up to time i-L, respectively; e (i) e (i-1) e (i-2) … e (i-L) are the delay of the output noise at time i and at intervals of 1 up to time i-L, respectively; n (i-1) n (i-2) … n (i-L) are the delay of the input noise at time i-1 and at an interval of 1 up to time i-L, respectively;
step B, constructing a class auxiliary variable, specifically: aiming at the condition that the input signal is a white Gaussian process or a colored Gaussian process and the output noise signal is colored noise, adjusting the parameters of the auxiliary variables to form a structural auxiliary variable based on the auxiliary variables, so that the class auxiliary variable xi at the j momentjInput signal vector p at the j time of EIV-IIR filter modeljStrong correlation;
wherein, the class auxiliary Variable, i.e. the instrument Variable-like, abbreviated as IV-like, must satisfy the following conditions to ensure the consistency of parameter estimation:
Figure BDA0002837805710000023
is a non-singular matrix
E[ξjv(j)]=0
Wherein v (j) is the composite noise at the j time;
aiming at the condition that the input signal is a white Gaussian process or a colored Gaussian process and the output noise signal is colored noise, constructing a j-th time auxiliary variable xijAs shown in (5):
ξj=[x(j-L-1)x(j-L-2)…x(j-2L)
x(j-1)x(j-2)…x(j-L)]T (5)
wherein x (j-L-1) x (j-L-2) … x (j-2L) is the delay of the noisy input signal at the j-L-1 th time and the interval of 1 to the j-2L th time respectively;
and C, calculating the parameter estimation value of IV-like under the EIV-IIR filter model through the step (6):
Figure BDA0002837805710000031
wherein,
Figure BDA0002837805710000032
representing an estimate of IV-like at time i,
Figure BDA0002837805710000033
representing an estimated value, y (j) representing a noisy output signal at the j-th time, y (j) calculating by (1), and replacing i in formula (1) by j;
step D, calculating the estimation deviation of parameter estimation by using IV-like under the condition of colored output noise, wherein the calculation process comprises the following steps: taking the limit on (6), we get (7):
Figure BDA0002837805710000034
wherein,
Figure BDA0002837805710000035
for the estimation of the IV-like method at the ith time
Figure BDA0002837805710000036
The deviation from the unknown system parameter vector h is recorded as: Δ h;
Figure BDA0002837805710000037
for unknown input noise variance, ILAn identity matrix of dimension L × L; from (7), it can be seen that, due to the presence of noise in the EIV-IIR filter model,
Figure BDA0002837805710000038
therefore,. DELTA.h.not equal to 0; since h is unknown, the deviation Δ h can be estimated by using the unbiased estimate at time i-1
Figure BDA0002837805710000039
Instead of obtaining the unknown system parameter vector h in the Δ h expression, the estimated value of the deviation Δ h can be recorded as
Figure BDA00028378057100000310
Step E, calculating the variance of the unknown input noise
Figure BDA00028378057100000311
Comprises the following substeps:
step E1, calculating the i-th backward input estimation vector aiIs estimated value of
Figure BDA00028378057100000312
Figure BDA00028378057100000313
Wherein,
Figure BDA00028378057100000314
the estimated values of the estimated vectors are input backwards at the ith and the i-1 time respectively,
Figure BDA00028378057100000315
is a vector of 0 s as an initial value,
Figure BDA00028378057100000316
as an i-th time class auxiliary variable xiiTransposing;
step E2, calculating the error estimated by the IV-like method
Figure BDA0002837805710000041
And backward input estimation error
Figure BDA0002837805710000042
The estimated value of the cross-correlation function g (i) at the ith moment obtained by multiplication
Figure BDA0002837805710000043
Figure BDA0002837805710000044
Wherein,
Figure BDA0002837805710000045
the estimated values of the cross-correlation function of the IV-like method estimated error and the backward input estimated error at the ith and the ith-1 time respectively,
Figure BDA0002837805710000046
the initial value is 0;
step E3, calculating the variance of the unknown input noise
Figure BDA0002837805710000047
Real-time estimation ofThe value:
Figure BDA0002837805710000048
wherein,
Figure BDA0002837805710000049
for unknown input noise variance
Figure BDA00028378057100000410
Real-time estimate of (a);
and F, compensating the deviation caused by the noise in the IV-like estimated value by using a deviation compensation principle, and calculating the unbiased estimation of the unknown system parameter vector h, wherein the unbiased estimation specifically comprises the formula (11):
Figure BDA00028378057100000411
wherein,
Figure BDA00028378057100000412
unbiased estimation of an unknown system parameter vector h at the ith moment;
and G, circulating the step B to the step F, executing iterative updating until a preset iteration number M is reached, and ending the method.
Advantageous effects
Compared with the prior art, the self-adaptive filtering method based on the deviation compensation auxiliary variable has the following beneficial effects that:
1. in the step B of the method, a class auxiliary variable is introduced, so that the problem that the method cannot work normally due to instability of the auxiliary variable when an input signal is a white Gaussian process is solved, the proposed class auxiliary variable is strongly correlated with an input signal vector when the input signal is the white Gaussian process or a colored Gaussian process, the proposed class auxiliary variable method is stable, and the application range of the auxiliary variable in practice is expanded;
2. in the step D of the method, the estimation deviation of parameter estimation is carried out by using an IV-like method in an EIV-IIR filter model through calculation, and the influence of the output noise variance on the estimation value of the unknown parameter is eliminated by the method shown in the step (7). compared with the deviation compensation RLS method, the method can work under the condition of colored output noise, only the input noise variance needs to be estimated, the unbiased estimation of the unknown parameter can be obtained, the complexity of the method is reduced, and the application range of the method is expanded;
3. the method can estimate the unknown input noise variance in real time, can accurately estimate the input noise variance under the condition of no need of input noise prior knowledge, and further realizes unbiased estimation of unknown parameters through deviation compensation.
Drawings
FIG. 1 is an EIV-IIR filter model used in step A of an adaptive filtering method based on bias compensation type auxiliary variables according to the present invention;
fig. 2 is a diagram of an adaptive filtering method based on offset compensation-like auxiliary variables, which is implemented specifically by comparing a recursive offset compensation-like auxiliary variable method with an RLS method, a recursive auxiliary variable method, and a recursive classification-like auxiliary variable method without offset compensation under the condition that an input signal is a white gaussian process and output noise is colored noise;
fig. 3 is a diagram of an adaptive filtering method based on bias compensation type auxiliary variables, which is implemented specifically, when an input signal is a colored gaussian process and an output noise is colored noise under the conditions, a recursive bias compensation type auxiliary variable method is performed, and a comparison is performed between an RLS method without bias compensation, a recursive auxiliary variable method, and a recursive classification type auxiliary variable method.
Detailed Description
The technical solution of the present invention will be described in detail below with reference to the accompanying drawings and examples. The present embodiment is implemented on the premise of the technical solution of the present invention, and a detailed implementation manner and a specific operation process are given, but the protection scope of the present invention is not limited to the following embodiments.
Examples
The present embodiment uses the mean square error criterion as a performance indicator, and is used for system identification for noisy inputs in a colored output noise environment.
Fig. 1 shows an EIV-IIR filter model adopted by the present invention, which is an adaptive filter model under a system identification framework. The following describes a specific implementation of the adaptive filtering method based on the offset compensation auxiliary variable proposed in the present invention with reference to fig. 1, which is summarized as follows:
step A, presetting iteration times M, and constructing a variable error-containing model;
wherein, the variable error-containing model, namely the EIV-IIR filter model, can be expressed as:
Figure BDA0002837805710000051
wherein, EIV is Errors in Variables, meaning that the Variables contain Errors; IIR is Infine Impulse Response which is Infinite Impulse Response; y (i) is the noisy output signal at time i, piAn input signal vector at the ith moment of the EIV-IIR filter model is represented as h, an L-order unknown system parameter vector to be estimated is represented as v (i), composite noise at the ith moment is represented as v (i), superscript T represents transposition operation, and L is the order of the EIV-IIR filter model;
vector piAnd v (i) in particular:
pi=[y(i-1)y(i-2)…y(i-L)x(i-1)x(i-2)…x(i-L)]T (2)
Figure BDA0002837805710000061
ni=[e(i-1)e(i-2)…e(i-L)n(i-1)n(i-2)…n(i-L)]T (4)
wherein y (i-1) y (i-2) … y (i-L) is the delay of the noisy output signal at time i-1 and at an interval of 1 up to time i-L, respectively; x (i-1) x (i-2) … x (i-L) are the delay of the noisy input signal at time i-1 and at an interval of 1 up to time i-L, respectively; e (i) e (i-1) e (i-2) … e (i-L) are the delay of the output noise at time i and at intervals of 1 up to time i-L, respectively; n (i-1) n (i-2) … n (i-L) are the delay of the input noise at time i-1 and at an interval of 1 up to time i-L, respectively;
step B, constructing a class auxiliary variable, specifically: aiming at the condition that the input signal is a white Gaussian process or a colored Gaussian process and the output noise signal is colored noise, adjusting the parameters of the auxiliary variables to form a structural auxiliary variable based on the auxiliary variables, so that the class auxiliary variable xi at the j momentjInput signal vector p at the j time of EIV-IIR filter modeljStrong correlation;
wherein, the class auxiliary Variable, i.e. the instrument Variable-like, abbreviated as IV-like, must satisfy the following conditions to ensure the consistency of parameter estimation:
Figure BDA0002837805710000062
is a non-singular matrix
E[ξjv(j)]=0
Wherein v (j) is the composite noise at the j time;
aiming at the condition that the input signal is a white Gaussian process or a colored Gaussian process and the output noise signal is colored noise, constructing a j-th time auxiliary variable xijAs shown in (5):
ξj=[x(j-L-1)x(j-L-2)…x(j-2L)
x(j-1)x(j-2)…x(j-L)]T (5)
wherein x (j-L-1) x (j-L-2) … x (j-2L) is the delay of the noisy input signal at the j-L-1 th time and the interval of 1 to the j-2L th time respectively;
and C, calculating the parameter estimation value of IV-like under the EIV-IIR filter model through the step (6):
Figure BDA0002837805710000071
wherein,
Figure BDA0002837805710000072
is shown asAn estimate of the IV-like at time i,
Figure BDA0002837805710000073
representing an estimated value, y (j) representing a noisy output signal at the j-th time, y (j) calculating by (1), and replacing i in formula (1) by j;
step D, calculating the estimation deviation of parameter estimation by using IV-like under the condition of colored output noise, wherein the calculation process comprises the following steps: taking the limit on (6), we get (7):
Figure BDA0002837805710000074
wherein,
Figure BDA0002837805710000075
for the estimation of the IV-like method at the ith time
Figure BDA0002837805710000076
The deviation from the unknown system parameter vector h is recorded as: Δ h;
Figure BDA0002837805710000077
for unknown input noise variance, ILAn identity matrix of dimension L × L; from (7), it can be seen that, due to the presence of noise in the EIV-IIR filter model,
Figure BDA0002837805710000078
therefore,. DELTA.h.not equal to 0; since h is unknown, the deviation Δ h can be estimated by using the unbiased estimate at time i-1
Figure BDA0002837805710000079
Instead of obtaining the unknown system parameter vector h in the Δ h expression, the estimated value of the deviation Δ h can be recorded as
Figure BDA00028378057100000710
Step E, calculating the variance of the unknown input noise
Figure BDA00028378057100000711
Comprises the following substeps:
step E1, calculating the i-th backward input estimation vector aiIs estimated value of
Figure BDA00028378057100000712
Figure BDA00028378057100000713
Wherein,
Figure BDA00028378057100000714
the estimated values of the estimated vectors are input backwards at the ith and the i-1 time respectively,
Figure BDA00028378057100000715
is a vector of 0 s as an initial value,
Figure BDA00028378057100000716
as an i-th time class auxiliary variable xiiTransposing;
step E2, calculating the error estimated by the IV-like method
Figure BDA00028378057100000717
And backward input estimation error
Figure BDA00028378057100000718
The estimated value of the cross-correlation function g (i) at the ith moment obtained by multiplication
Figure BDA00028378057100000719
Figure BDA00028378057100000720
Wherein,
Figure BDA00028378057100000721
the estimated values of the cross-correlation function of the IV-like method estimated error and the backward input estimated error at the ith and the ith-1 time respectively,
Figure BDA0002837805710000081
the initial value is 0;
step E3, calculating the variance of the unknown input noise
Figure BDA0002837805710000082
Real-time estimation of (d):
Figure BDA0002837805710000083
wherein,
Figure BDA0002837805710000084
for unknown input noise variance
Figure BDA0002837805710000085
Real-time estimate of (a);
and F, compensating the deviation caused by the noise in the IV-like estimated value by using a deviation compensation principle, and calculating the unbiased estimation of the unknown system parameter vector h, wherein the unbiased estimation specifically comprises the formula (11):
Figure BDA0002837805710000086
wherein,
Figure BDA0002837805710000087
unbiased estimation of an unknown system parameter vector h at the ith moment;
and G, circulating the step B to the step F, executing iterative updating until a preset iteration number M is reached, and ending the method.
Simulation experiment
The effect of the invention is verified by the following simulation experiment:
the output noise is colored noise, e (i) ═ u (i) -0.3u (i-1), where u (i) is white noise with a mean of 0 and a variance of 1. The input noise n (i) is zero mean,
Figure BDA0002837805710000088
white gaussian noise. The input signal is divided into two cases: (1) input signal
Figure BDA0002837805710000089
A white gaussian process with zero mean with a variance of 1; (2) input signal
Figure BDA00028378057100000810
In order to be a colored gaussian process,
Figure BDA00028378057100000811
Figure BDA00028378057100000812
the adaptive IIR filter can be described by the following model:
Figure BDA00028378057100000813
the iteration number M of the simulation experiment is 20000, the independent experiment number is 100, and the mean square error criterion is used as the performance index.
Fig. 2 shows the comparison of the unknown system parameter estimation accuracy based on the recursive deviation compensation type auxiliary variable method and the RLS method, the recursive auxiliary variable method and the recursive classification type auxiliary variable method which do not perform deviation compensation in the case that the input signal is a white gaussian process and the output noise is colored noise in the embodiment of the present invention. It is seen from the figure that after the deviation caused by the noise is compensated, the estimation accuracy of the deviation compensation type auxiliary variable adaptive filtering method provided by the invention is obviously improved, and the unbiased estimation of unknown parameters can be realized.
Fig. 3 shows a comparison of the estimation accuracy of unknown system parameters based on the recursive bias compensation type auxiliary variable method, the RLS method without bias compensation, the recursive auxiliary variable method, and the recursive classification type auxiliary variable method in the embodiment of the present invention, when the input signal is a colored gaussian process and the output noise is colored noise. As seen from the figure, the deviation compensation type auxiliary variable adaptive filtering method provided by the invention has the highest estimation precision, and can realize unbiased estimation on unknown system parameters under the condition of colored output noise.
The above-mentioned embodiments are merely preferred embodiments of the present invention, and the present invention is not limited thereto, and the technical means disclosed in the present invention is not limited to the technical means disclosed in the above-mentioned embodiments, but also includes technical means formed by any combination of the above technical features. It should be understood that any modification, equivalent replacement, and improvement made by those skilled in the art without departing from the principle of the present invention shall fall within the protection scope of the present invention.

Claims (6)

1. An adaptive filtering method based on deviation compensation auxiliary variables is characterized in that: the method comprises the following steps:
step A, presetting iteration times M, and constructing a variable error-containing model;
step B, constructing a class auxiliary variable, specifically: aiming at the condition that the input signal is a white Gaussian process or a colored Gaussian process and the output noise signal is colored noise, adjusting the parameters of the auxiliary variables to form a structural auxiliary variable based on the auxiliary variables, so that the class auxiliary variable xi at the j momentjInput signal vector p at the j time of EIV-IIR filter modeljStrong correlation;
wherein, the auxiliary Variable of class, namely the instrument Variable-like, is abbreviated as IV-like;
aiming at the condition that the input signal is a white Gaussian process or a colored Gaussian process and the output noise signal is colored noise, constructing a j-th time auxiliary variable xijAs shown in (5):
ξj=[x(j-L-1)x(j-L-2)…x(j-2L)
x(j-1)x(j-2)…x(j-L)]T (5)
wherein x (j-L-1) x (j-L-2) … x (j-2L) is the delay of the noisy input signal at the j-L-1 th time and the interval of 1 to the j-2L th time respectively;
step C, calculating the parameter estimation value of IV-like under the EIV-IIR filter model
Figure FDA0002837805700000011
Step D, calculating the estimation deviation of parameter estimation by using IV-like under the condition of colored output noise, wherein the calculation process comprises the following steps: taking the limit on (6), we get (7):
wherein,
Figure FDA0002837805700000012
for the estimation of the IV-like method at the ith time
Figure FDA0002837805700000013
The deviation from the unknown system parameter vector h is recorded as: Δ h;
Figure FDA0002837805700000014
for unknown input noise variance, ILAn identity matrix of dimension L × L; from (7), it can be seen that, due to the presence of noise in the EIV-IIR filter model,
Figure FDA0002837805700000015
therefore,. DELTA.h.not equal to 0; since h is unknown, the deviation Δ h can be estimated by using the unbiased estimate at time i-1
Figure FDA0002837805700000016
Instead of obtaining the unknown system parameter vector h in the Δ h expression, the estimated value of the deviation Δ h can be recorded as
Figure FDA0002837805700000017
Step E, calculating the variance of the unknown input noise
Figure FDA0002837805700000018
Comprises the following substeps:
step E1, calculating the i-th backward input estimation vector aiIs estimated value of
Figure FDA0002837805700000019
Figure FDA00028378057000000110
Wherein,
Figure FDA0002837805700000021
the estimated values of the estimated vectors are input backwards at the ith and the i-1 time respectively,
Figure FDA0002837805700000022
is a vector of 0 s as an initial value,
Figure FDA0002837805700000023
as an i-th time class auxiliary variable xiiTransposing;
step E2, calculating the error estimated by the IV-like method
Figure FDA0002837805700000024
And backward input estimation error
Figure FDA0002837805700000025
The estimated value of the cross-correlation function g (i) at the ith moment obtained by multiplication
Figure FDA0002837805700000026
Figure FDA0002837805700000027
Wherein,
Figure FDA0002837805700000028
the estimated values of the cross-correlation function of the IV-like method estimated error and the backward input estimated error at the ith and the ith-1 time respectively,
Figure FDA0002837805700000029
the initial value is 0;
step E3, calculating the variance of the unknown input noise
Figure FDA00028378057000000210
Real-time estimation of (d):
Figure FDA00028378057000000211
wherein,
Figure FDA00028378057000000212
for unknown input noise variance
Figure FDA00028378057000000213
Real-time estimate of (a);
f, compensating the deviation caused by the noise in the IV-like estimated value by using a deviation compensation principle, namely calculating the unbiased estimation of the parameter vector h of the unknown system;
and G, circulating the step B to the step F, executing iterative updating until a preset iteration number M is reached, and ending the method.
2. The adaptive filtering method based on the bias compensation auxiliary-like variable according to claim 1, characterized in that: in step a, the variable error-containing model, i.e., the EIV-IIR filter model, may be expressed as:
Figure FDA00028378057000000214
wherein, EIV is Errors in Variables, meaning that the Variables contain Errors; IIR is Infine Impulse Response which is Infinite Impulse Response; y (i) is the noisy output signal at time i, piThe method comprises the steps of obtaining an input signal vector of an EIV-IIR filter model at the moment i, h an unknown system parameter vector of an L order to be estimated, v (i) composite noise of the moment i, superscript T representing transposition operation, and L an order of the EIV-IIR filter model.
3. The adaptive filtering method based on the bias compensation auxiliary-like variable according to claim 1, characterized in that: in step A, vector piAnd v (i) in particular:
pi=[y(i-1)y(i-2)…y(i-L)x(i-1)x(i-2)…x(i-L)]T(2)
Figure FDA00028378057000000215
wherein,
ni=[e(i-1)e(i-2)…e(i-L)n(i-1)n(i-2)…n(i-L)]T (4)
y (i-1) y (i-2) … y (i-L) are the delay of the noisy output signal at time i-1 and at an interval of 1 up to time i-L, respectively; x (i-1) x (i-2) … x (i-L) are the delay of the noisy input signal at time i-1 and at an interval of 1 up to time i-L, respectively; e (i) e (i-1) e (i-2) … e (i-L) are the delay of the output noise at time i and at intervals of 1 up to time i-L, respectively; n (i-1) n (i-2) … n (i-L) are the delay of the input noise at time i-1 and at intervals of 1 up to time i-L, respectively.
4. The adaptive filtering method based on the bias compensation auxiliary-like variable according to claim 3, characterized in that: in step B, the class auxiliary variables must satisfy the following conditions to ensure consistency of parameter estimation:
Figure FDA0002837805700000031
is a non-singular matrix
E[ξjv(j)]=0
Where v (j) is the composite noise at time j.
5. The adaptive filtering method based on the bias compensation auxiliary-like variable according to claim 4, wherein: and C, calculating the parameter estimation value of IV-like under the EIV-IIR filter model through the step (6)
Figure FDA0002837805700000032
Figure FDA0002837805700000033
Wherein,
Figure FDA0002837805700000034
representing an estimate of IV-like at time i,
Figure FDA0002837805700000037
representing an estimated value, y (j) representing a noisy output signal at the j-th time, y (j) calculating by (1), and replacing i in formula (1) by j;
6. the adaptive filtering method based on the bias compensation auxiliary-like variable according to claim 5, characterized in that: in step F, calculating an unbiased estimate of the unknown system parameter vector h by (11):
Figure FDA0002837805700000035
wherein,
Figure FDA0002837805700000036
And carrying out unbiased estimation on the unknown system parameter vector h at the ith moment.
CN202011481813.9A 2020-12-15 2020-12-15 Self-adaptive filtering method based on deviation compensation auxiliary variable Active CN112583381B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011481813.9A CN112583381B (en) 2020-12-15 2020-12-15 Self-adaptive filtering method based on deviation compensation auxiliary variable

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011481813.9A CN112583381B (en) 2020-12-15 2020-12-15 Self-adaptive filtering method based on deviation compensation auxiliary variable

Publications (2)

Publication Number Publication Date
CN112583381A true CN112583381A (en) 2021-03-30
CN112583381B CN112583381B (en) 2022-07-15

Family

ID=75135396

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011481813.9A Active CN112583381B (en) 2020-12-15 2020-12-15 Self-adaptive filtering method based on deviation compensation auxiliary variable

Country Status (1)

Country Link
CN (1) CN112583381B (en)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20150323232A1 (en) * 2012-12-14 2015-11-12 Daikin Industries, Ltd. Active filter control device
CN106656881A (en) * 2016-11-22 2017-05-10 北京理工大学 Adaptive blind equalization method based on deviation compensation
CN107592096A (en) * 2017-09-29 2018-01-16 苏州大学 A kind of robust deviation compensation sef-adapting filter and its filtering method

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20150323232A1 (en) * 2012-12-14 2015-11-12 Daikin Industries, Ltd. Active filter control device
CN106656881A (en) * 2016-11-22 2017-05-10 北京理工大学 Adaptive blind equalization method based on deviation compensation
CN107592096A (en) * 2017-09-29 2018-01-16 苏州大学 A kind of robust deviation compensation sef-adapting filter and its filtering method

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
贾丽娟 等: "有噪输入输出动态系统辨识方法", 《北京理工大学学报》 *

Also Published As

Publication number Publication date
CN112583381B (en) 2022-07-15

Similar Documents

Publication Publication Date Title
CN109163720A (en) Kalman filter tracking method based on fading memory exponent
CN107592096B (en) Filtering method of robust deviation compensation adaptive filter
CN109871508A (en) There are the centralized two stages kalman estimate methods of correlated measurement noise
CN108512528B (en) Ratio control and normalization LMP filtering method under a kind of CIM function
US20100030529A1 (en) System identification method and program, storage medium, and system identification device
CN110398942B (en) Parameter identification method for industrial production process control
CN113467236B (en) Method for time lag compensation of error signal
CN117896468B (en) Deviation compensation echo cancellation method and system for telephone communication
CN113670315B (en) Variation iterative Kalman filtering-based method for estimating attitude of Liqun heavy tail interference noise dynamic aircraft
CN111444474A (en) Multiplicative noise correlation self-adaptive CKF-based target tracking method
CN112583381B (en) Self-adaptive filtering method based on deviation compensation auxiliary variable
CN114221638A (en) Maximum cross-correlation entropy Kalman filtering method based on random weighting criterion
CN114357370A (en) Solving time-varying matrix square root for noise-suppressed nulling neural networks with non-convex activation functions
CN117789691A (en) Echo cancellation method and system for telephone communication
CN112866149B (en) Adaptive estimation method based on bias compensation block sparse proportional NLMS
CN116680500B (en) Position estimation method and system of underwater vehicle under non-Gaussian noise interference
CN113503879A (en) Dynamic adaptive Kalman filter method based on ensemble empirical mode decomposition
CN102724152A (en) Multinomial adaptive active noise cancellation method based on Laguerre structure
CN110989353B (en) Design method of periodic disturbance observer
CN109582914B (en) Parallel fusion estimation method of noise-related deviation system
CN115800957A (en) Deviation compensation adaptive filtering method based on matrix eigenvalue solution
CN109840069B (en) Improved self-adaptive fast iterative convergence solution method and system
CN110034746B (en) Kalman filtering method based on maximum collaborative entropy
Vergés et al. Optimal linear mean square filter for Markov jump linear systems with hidden Markov chain and partial observations of the state
CN108512529A (en) Ratio control under a kind of μ rule function and normalization LMP filtering methods

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