CN108614804B - Regularization Kalman filtering method based on signal-to-noise ratio test - Google Patents
Regularization Kalman filtering method based on signal-to-noise ratio test Download PDFInfo
- Publication number
- CN108614804B CN108614804B CN201810379450.4A CN201810379450A CN108614804B CN 108614804 B CN108614804 B CN 108614804B CN 201810379450 A CN201810379450 A CN 201810379450A CN 108614804 B CN108614804 B CN 108614804B
- Authority
- CN
- China
- Prior art keywords
- state
- parameters
- parameter
- signal
- estimated
- 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.)
- Expired - Fee Related
Links
- 238000001914 filtration Methods 0.000 title claims abstract description 58
- 238000000034 method Methods 0.000 title claims abstract description 44
- 238000012360 testing method Methods 0.000 title claims abstract description 14
- 239000011159 matrix material Substances 0.000 claims description 59
- 238000012937 correction Methods 0.000 claims description 8
- 239000000126 substance Substances 0.000 claims description 7
- 230000008569 process Effects 0.000 claims description 3
- 230000007704 transition Effects 0.000 claims description 3
- 238000004422 calculation algorithm Methods 0.000 description 24
- 230000000694 effects Effects 0.000 description 8
- 230000006835 compression Effects 0.000 description 4
- 238000007906 compression Methods 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 230000002411 adverse Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000001010 compromised effect Effects 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000001575 pathological effect Effects 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 201000004569 Blindness Diseases 0.000 description 1
- 241000764238 Isis Species 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 238000010827 pathological analysis Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 239000013598 vector Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Life Sciences & Earth Sciences (AREA)
- Operations Research (AREA)
- Probability & Statistics with Applications (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Algebra (AREA)
- Evolutionary Biology (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Bioinformatics & Computational Biology (AREA)
- Filters That Use Time-Delay Elements (AREA)
- Complex Calculations (AREA)
Abstract
The invention provides a regularization Kalman filtering method based on signal-to-noise ratio test. The method comprises the following steps: step 1, obtaining t k+1 State parameter X to be estimated at time k+1 Initial state estimation ofAnd a typical parameter theta k+1 Least squares estimation ofStep 2, according to the state parameter X to be estimated k+1 The ith parameter ofState estimation ofDeterminingThe signal-to-noise ratio statistic of (1), 2 …, t-1, t, t is the number of state parameters to be estimated; step 3, according to the aboveThe signal-to-noise ratio statistic of (2), the state parameter X to be estimated k+1 Dividing into interference-related state parameters and non-interference-related state parameters; step 4, estimating according to least squareDetermining ridge parameters of disturbance-related state parametersRidge parameters with non-interference state parametersStep 5, according to the ridge parametersAnddetermining a state parameter X to be estimated k+1 To estimate said initial stateAnd (6) correcting. The invention effectively reduces the introduction of deviation while reducing the estimation variance of the state parameters, so that the estimation result of the state parameters is better in the mean square error sense.
Description
Technical Field
The invention relates to the technical field of dynamic data processing, in particular to a regularization Kalman filtering method based on signal-to-noise ratio test.
Background
Kalman filtering is one of the most common methods for dynamic data processing, and has been extensively studied and widely applied in geodetic surveying, satellite navigation, and satellite orbit determination. The ill-conditioned nature of the observation matrix in the discrete dynamic system can produce very big influence to Kalman filtering state estimation, in order to overcome the influence of the ill-conditioned nature of the observation matrix, improve the precision of parameter estimation, many scholars have given improved algorithm.
Currently, there are some methods proposed by scholars to solve the ill-posed problem in discrete dynamic systems from the perspective of biased estimation: (1) the Biased estimation is combined with Kalman filtering to provide a Biased Kalman Filter. (2) The ridge regression and Kalman filtering are combined, and the gain array is corrected, so that the adverse effect of observation matrix ill-conditioned on a filtering value is overcome. (3) Combining Stein compression estimation and ridge regression with Kalman filtering, providing corresponding compression type Kalman filtering and ridge type Kalman filtering and their algorithms, and providing a method for selecting compression coefficients and ridge parameters.
In fact, the estimated hazard size of each state parameter is different according to the ill-condition of the observation matrix, and the hazard size is related to the magnitude of the state parameter itself and the degree of the complex collinearity of the observation matrix data column corresponding to the parameter.
Disclosure of Invention
The invention provides a two-parameter ridge Kalman filtering method based on signal-to-noise ratio measurement, which measures the influence of the ill-conditioned property of parameter state estimation by using signal-to-noise ratio statistics, takes targeted measures on a Kalman filtering recursion process according to a measurement result, improves a ridge Kalman filtering algorithm, and further reduces the influence of the ill-conditioned property of an observation matrix on the state estimation.
The invention provides a regularization Kalman filtering method based on signal-to-noise ratio test, which comprises the following steps:
step 4, estimating according to least squareDetermining ridge parameters of disturbance-related state parametersRidge parameters with non-interference state parameters
Step 5, according to the ridge parametersAnddetermining a state parameter X to be estimated k+1 To estimate said initial stateAnd (6) correcting.
Further, the method further comprises:
obtaining t k+1 Normal matrix N of time instants k+1 ;
If the judgment result is positive matrix N k+1 Is greater than the preset condition number threshold, then for t k+1 And carrying out signal-to-noise ratio statistics on each state parameter at each moment.
Further, the step 2 specifically includes:
according to the formula
DeterminingSignal to noise ratio statistic ofWherein the content of the first and second substances,is composed ofThe variance of (c).
Further, the step 3 specifically includes:
Further, the step 4 specifically includes:
Determining; wherein the content of the first and second substances,to representThe maximum value of (a) is,s is the number of interference-related state parameters.
Further, the state parameter X to be estimated k+1 The correction matrix of (a) is:
the invention has the beneficial effects that:
the regularization Kalman filtering method based on signal-to-noise ratio test provided by the invention divides the state parameters to be estimated into interference parameters and non-interference parameters according to the magnitude of the signal-to-noise ratio statistic of each state parameter; and then estimate based on least squaresDetermining two ridge parametersAndso as to perform different intensity ridge-type corrections on the state estimates of the two part parameters. For interference parameters, the interference parameters are correctedRelatively large, for non-involved parameters, make it correct ridge parametersIs relatively small. The refined processing effectively reduces the introduction of deviation in ridge type Kalman filtering while reducing the state parameter estimation variance, so that the state parameter estimation result is better in the mean square error sense.
Drawings
FIG. 1 is a schematic flow chart of a regularization Kalman filtering method based on signal-to-noise ratio test according to an embodiment of the present invention;
FIG. 2 is a diagram illustrating a variation of condition numbers of a normal matrix in each filtering method according to an embodiment of the present invention;
fig. 3 is a schematic diagram illustrating changes of mean square errors in filtering methods according to an embodiment of the present invention;
fig. 4 is a schematic diagram illustrating changes in state estimation errors in each filtering method according to an embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly described below with reference to the accompanying drawings in the embodiments of the present invention, and it is obvious that the described embodiments are some, but not all, embodiments of the present invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
Before describing the regularized Kalman filtering method based on signal-to-noise ratio test provided by the invention, a discrete dynamic system, a Kalman filtering elementary equation, a ridge-type Kalman filtering state estimation and the influence of the ill-conditioned observation matrix in the discrete dynamic system on the Kalman filtering state estimation are introduced in the following.
The embodiment of the invention adopts the following state space model to describe a discrete dynamic system:
X k+1 =Φ k+1,k X k +W k (1)
Y k =H k X k +V k (2)
in the formula, X k ∈R t For discrete dynamic systems at time t k The state of (1); y is k ∈R n Is the corresponding observed signal; phi k+1,k Is a t × t nonsingular matrix, is t k Time to t k+1 A one-step state transition matrix of the moment; h k Is an nxt observation matrix; w k ∈R t Input noise that is normally distributed; v k ∈R n Observed noise that follows a normal distribution. Equation (1) is called a state equation, and equation (2) is called an observation equation.
At the same time, assume W k And V k Satisfy the requirement of
E(W k )=0,Cov(W k ,W j )=E(W k W j T )=Q k δ kj
E(V k )=0,Cov(V k ,V j )=E(V k V j T )=R k δ kj
Cov(W k ,V j )=E(W k V j T )=0 (3)
In the formula, Q k For a covariance matrix of the input noise, assume Q k Is a non-negative array, R k To observe the covariance matrix of the noise, let R be assumed k Is a positive definite matrix; delta kj Is a Kronecher-delta function, which is defined as
As can be seen from formula (3), W k And V k Mean value of zero and covariance matrix of Q k And R k Is white uncorrelated noise.
The following relevant introduction is made to the basic equation of the kalman filter algorithm:
and (3) state one-step prediction:
one-step prediction covariance matrix:
a filter gain matrix:
and (3) updating the state:
estimating an error variance matrix:
equations (4) to (8) are basic equations of recursive kalman filtering. Given an initial valueAnd P 0 According to t k+1 Observed quantity Y of time k+1 Then t can be calculated by recursion k+1 State estimation of time of day
From the basic Kalman filtering equation, t k+1 The state estimate for a time instant can in turn be expressed as:
is provided withN k+1 A normal matrix called kalman filtering. N is a radical of k+1 Is non-negative array, for N k+1 The decomposition of the characteristic value is carried out,
Then
wherein, the first and the second end of the pipe are connected with each other,
if the observation matrix H k+1 Presence of pathological conditions, then H k+1 Andis likely to be such that N k+1 The existence of the ill-conditioned character, which is shown by the practical work,the control effect on the ill-conditioned nature of the observation matrix is weak and cannot eliminate the adverse effect of the ill-conditioned nature of the observation matrix on the state estimation.
From the formula (12), N can be seen k+1 Presence of one or more small eigenvaluesAt this time, if Y k+1 Andthere is a small observed error or bias to which the inverse of the small eigenvalue in the above equation amplifies, thereby causing the estimate to deviate far from the true value.
In order to reduce the influence of small characteristic values on the observation result, a ridge type Kalman filtering method is provided on the basis of the algorithm, t k+1 The ridge-type kalman filter state estimation at time is:
ridge-type Kalman filtering utilizes a ridge parameter α k+1 For small eigenvalueAnd (4) suppressing is carried out, and the variance of state estimation is reduced, so that the amplification effect of small characteristic values on observation errors is weakened.
However, the ridge-type kalman filter has two drawbacks, one is that no ill-conditioned information is utilized, so that correction of ridge parameters is blindness, and all parameters are corrected completely and consistently; secondly, deviation is introduced by the ridge-type Kalman filtering, and the deviation introduced by the ridge-type Kalman filtering is possibly amplified in the continuous recursion process, so that the estimation precision is influenced.
From the above analysis, t k+1 State estimation of time of dayAssociated with the normal matrix (see equation (10)), if normal matrix N k+1 If a pathological condition exists, t is k+1 The state estimate at the time will become very unstable, which is where the ill-conditioned nature of the observation matrix exerts an influence on the state estimate. Normal matrix N k+1 The reason for the morbidity is that the complex collinearity relationship exists between the data columns, thereby causing t k+1 The state estimation at the moment is not good, but not all the state estimation is not good, the ill-conditioned property of the normal matrix only has larger influence on the state parameter estimation corresponding to the data columns participating in the complex collinearity, and has smaller influence on the state parameter estimation corresponding to the data columns not participating in the complex collinearity.
The invention provides a regularization Kalman filtering method (DPRTKF for short) based on signal-to-noise ratio test on the basis of a ridge Kalman filtering method, which corrects the state parameters of two parts with smaller and larger ill-conditioned observed matrixes with different strengths. Fig. 1 is a schematic flow chart of a regularization kalman filtering method based on signal-to-noise ratio test according to an embodiment of the present invention. As shown in fig. 1, the method comprises the steps of:
s101, acquiring t k+1 State parameter X to be estimated at time k+1 Initial state estimation ofAnd a typical parameter theta k+1 Least squares estimation of
Specifically, the method provided by the invention needs to acquire t in advance k+1 State parameter X to be estimated at time k+1 Initial state estimation ofAnd a typical parameter theta k+1 Least squares estimation ofSince the invention aims at estimating the state obtained by the existing Kalman filtering methodThe correction is performed, so that the initial state estimation obtained by using the existing Kalman filtering method needs to be obtained in advanceAnd a typical rule parameter theta k+1 Least squares estimation of
For example, given initial values of state parametersAnd mean square error thereofUsing the kalman filter primitive equation:
S102, according to the state parameter X to be estimated k+1 The ith parameter ofState estimation ofDeterminingThe signal-to-noise ratio statistic of (1), 2 …, t-1, t, t is the number of state parameters to be estimated;
s103, according to theThe signal-to-noise ratio statistic of the state parameter X to be estimated k+1 Dividing into interference-related state parameters and non-interference-related state parameters;
s104, estimating according to least squaresDetermining ridge parameters of disturbance-related state parametersRidge parameters with non-interference state parameters
S105, according to the ridge parametersAnddetermining a state parameter X to be estimated k+1 To estimate said initial stateAnd (6) correcting.
The regularization Kalman filtering method based on signal-to-noise ratio test provided by the invention divides the state parameters to be estimated into interference parameters and non-interference parameters according to the magnitude of the signal-to-noise ratio statistic of each state parameter; and then estimate based on least squaresDetermining two ridge parametersAndso as to perform different intensity ridge-type corrections on the state estimates of the two part parameters. For interference parameters, the parameters are correctedRelatively large, for non-involved parameters, make it correct ridge parametersIs relatively small. The refined processing effectively reduces the introduction of deviation in ridge Kalman filtering while reducing the state parameter estimation variance, so that the state parameter estimation result is better in mean square error sense.
On the basis of the above embodiment, the method further comprises the steps of:
obtaining t k+1 Normal matrix N of time instants k+1 ;
If judging the learning method matrix N k+1 Is greater than the preset condition number threshold, then for t k+1 Each state parameter at a momentThe numbers are counted for signal to noise ratio.
In particular, the Fa matrix N k+1 If the condition number of the observation matrix is greater than the preset condition number threshold, the observation matrix is considered to be a sick matrix, and the ill-condition of the observation matrix can adversely affect the state estimation result of the state parameter to be estimated, so that the t needs to be estimated k+1 And carrying out signal-to-noise ratio statistics on each state parameter at each moment. Therefore, the state parameters to be estimated can be grouped for the following signal-to-noise ratio statistics based on each state parameter, and the state parameters to be estimated are divided into non-interference state parameters with small influence and interference state parameters with large influence.
On the basis of the foregoing embodiments, the step S102 in the method specifically includes: according to the formula
DeterminingSignal to noise ratio statistic ofWherein the content of the first and second substances,is composed ofThe variance of (c).
On the basis of the foregoing embodiments, the step S103 in the method specifically includes:
In particular, whenConsidering the complex collinearity to estimate the corresponding stateThe damage of the (A) is serious, and the estimation effect is not good; when in useConsidering the complex collinearity to estimate the corresponding stateThe damage of (2) is small, and the estimation effect is good. By calculating the SNR statistic, t can be calculated k+1 The state parameter of the moment is divided into two partsThe s parameters for which the SNR statistic is smaller areThe complex collinearity is greatly damaged and is called as interference-related state parameter; the t-s parameters with larger signal-to-noise ratio statistics areThey are less compromised by complex collinearity and are referred to as non-interfering state parameters. For theBecause the estimation is poor, the parameters of the part are greatly corrected; for theWith minor modifications thereto. In practice, the threshold valueCan be determined flexibly according to the actual situation.
On the basis of the foregoing embodiments, step S104 in the method specifically includes:
Determining; wherein the content of the first and second substances,to representThe maximum value of (a) is,s is interferenceThe number of state parameters.
Specifically, all signal-to-noise ratio statistics are combinedRearranged in descending order and numbered asOrder to
Is the signal to noise ratio, then its inverseThe larger the noise-to-signal ratio, the greater the degree to which the estimate of the corresponding state parameter is compromised by complex collinearity. By usingNoise signal ratio meanAndnoise signal ratio meanQuantitatively reacting with the ratio ofAndthe proportional relation between the complex collinearity harmfulness degrees of the parameter estimation is adopted to determine c k+1 。
On the basis of the above embodiments, the state parameter X to be estimated in the method k+1 The correction matrix of (a) is:
in particular, the correction matrix Z k+1 For t rows and t columns matrix, the first s rows of the matrix are all main diagonalsFor correcting disturbance-related state parametersThe front and back t-s rows of the matrix are main diagonalsFor correcting non-interference state parameters
The invention is further illustrated by the following further specific examples.
and 3, updating the state:
step 4. judging method matrixIf the condition number is greater than 500, if so, go to step 5, otherwise return to step 2.
And 5, measuring the signal-to-noise ratio, and determining interference-related state parameters and non-interference-related state parameters.
And 7, correcting the state estimation by utilizing the two-parameter ridge Kalman filtering, and calculating a mean square error matrix
And estimating the state parameters.
From the above, it can be seen that:
and because
at the same time
Therefore, it is not only easy to useIs a state parameter X k A compressed biased estimation of (2).
In addition to this, the present invention is,
And because, when selecting the appropriate initial valueAnd mean square error thereofTime, Kalman filteringCan be an unbiased estimation, while the ordinary ridge-type Kalman filteringAnd regularized filtering provided by the present inventionAre all Kalman filteringIs estimated by compression ofFiltering KalmanHas a compressive strength greater thanFiltering KalmanThe compressive strength of (a), so from a deflection point of view,is greater thanTo be smaller than
Due to the fact that
So in the condition number senseAlgorithm condition number less than Kalman filter estimationLarger than ridge Kalman filteringThe algorithm condition number.
In conclusion, the regularization Kalman filtering method based on signal-to-noise ratio test provided by the invention not only has the anti-attitude property which is not possessed by the common Kalman filtering, but also overcomes the blind and excessive anti-attitude property of the common ridge Kalman filtering, thereby reducing the damage of excessive deviation caused by excessive anti-attitude property and achieving more accurate and proper effect.
The invention is further illustrated by the following further specific examples.
To verify the correctness of the pathological analysis herein and the validity of the new algorithm, the following algorithm was used to perform the simulation test. Setting a one-step State transition TorqueMatrix phi k+1,k Is composed of
Design observation matrix H k Of 5 column vectors of
a k1 =[15.57 44.02 20.42 18.74 49.20 44.92 55.48 59.28 94.39 128.02 96.00 131.42 127.21 252.90 409.20 463.70 510.22] T
a k2 =[2643 2048 3940 6505 5723 11520 5779 5969 8461 20106 13313 10771 45543 36194 34703 39204 86533] T
a k4 =[18.0 9.5 12.8 36.7 35.7 24.0 43.3 46.7 76.7 180.5 60.9 103.7 126.8 157.7 169.4 331.4 371.6] T
a k5 =[4.45 6.92 4.28 3.90 5.50 4.60 5.62 5.15 6.18 6.15 5.88 4.68 4.88 5.57 10.78 7.05 6.35] T
a k3 =2a k1 +0.5a k4 +e k ,e k ~N 17 (0,0.05 2 I)
H k =[a k1 a k2 a k3 a k4 a k5 ]
Q k =I 5 ,R k =0.5 2 ×I 17 Covariance matrices, initial values, of state noise sequence and observation noise sequence, respectivelyIs composed ofWherein X 0 =[200,15,35,16,-2.8,6] Τ In the true value, the value of,is the initial error matrix.
In the embodiment of the invention, the data columns of the observation matrix are set to have approximate linear relation, so that the normal matrix has serious ill-conditioned property. Three filtering algorithms are used for calculation and their condition numbers, mean square errors and deviations are compared, and the effect pair is shown in fig. 2 to 4.
(1) As can be seen from fig. 2, the normal Kalman filter algorithm has a condition number of 10 for the normal Kalman filter algorithm 4 The order of magnitude of (2) exceeds a threshold value of 500, which shows that complex collinearity of an observation matrix can really cause ill-conditioned law matrix, and the condition numbers of law matrices of an RTKF algorithm (namely, ordinary ridge-type Kalman filtering in the figure) and a DPRTKF algorithm (namely, two-parameter ridge-type Kalman filtering in the figure) are effectively controlled.
(2) It can be seen from fig. 3 that the RTKF algorithm and the DPRTKF algorithm effectively weaken the ill-posed nature in the Kalman filtering, and the DPRTKF algorithm is superior to the general Kalman filtering algorithm and the RTKF algorithm in the mean square error sense.
(3) As can be seen from fig. 4, compared with the Kalman filter algorithm, the RTKF algorithm properly sacrifices the unbiased property of the state estimation in exchange for the great reduction of the variance, so that the estimation result is relatively stable, while the DPRTKF algorithm reduces the introduction of the deviation while reducing the state estimation variance, so that the euclidean distance between the solution and the true value is smaller than the Kalman filter algorithm and the RTKF algorithm.
Finally, it should be noted that: the above examples are only intended to illustrate the technical solution of the present invention, but not to limit it; although the present invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some technical features may be equivalently replaced; and such modifications or substitutions do not depart from the spirit and scope of the corresponding technical solutions of the embodiments of the present invention.
Claims (4)
1. A regularization Kalman filtering method based on signal-to-noise ratio test is applied to a state parameter estimation process of a pathologically observed matrix in a discrete dynamic system, and a state space model is adopted to describe the discrete dynamic system:
X k+1 =Φ k+1,k X k +W k (1)
Y k =H k X k +V k (2)
in the formula, X k ∈R t For discrete dynamic systems at time t k The state of (1); y is k ∈R n Is the corresponding observed signal; phi k+1,k Is a t × t nonsingular matrix, is t k Time to t k+1 A one-step state transition matrix of the moment; h k Is an nxt observation matrix; w k ∈R t Input noise that is normally distributed; v k ∈R n Observed noise that is normally distributed; wherein the discrete dynamic system is any one of a geodetic surveying system, a satellite navigation system and a satellite orbit determination system, and the method comprises:
step 1, obtaining t k+1 State parameter X to be estimated at time k+1 Initial state estimation ofAnd a typical parameter theta k+1 Least squares estimation of
Step 2, according to the state parameter X to be estimated k+1 The ith parameter ofState estimation ofDeterminingThe signal-to-noise ratio statistic of (1), 2 …, t-1, t, t is the number of state parameters to be estimated; the step 2 specifically comprises the following steps:
according to the formula
DeterminingSignal to noise ratio statistic ofWherein the content of the first and second substances,is composed ofThe variance of (a);
step 3, according to the aboveThe signal-to-noise ratio statistic of (2), the state parameter X to be estimated k+1 Dividing into interference-related state parameters and non-interference-related state parameters; the method specifically comprises the following steps:
whereinOmega is the level of significance and is,is center x 2 Distribution ofUpper omega quantile of (a);
step 4, estimating according to least squareDetermining ridge parameters of disturbance-related state parametersRidge parameters with non-interference state parameters
2. The method of claim 1, further comprising, prior to step 1:
obtaining t k+1 Normal matrix N of time instants k+1 ;
If judging the learning method matrix N k+1 Is greater than the preset condition number threshold, then for t k+1 And carrying out signal-to-noise ratio statistics on each state parameter at each moment.
3. The method according to claim 1, wherein step 4 is specifically:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810379450.4A CN108614804B (en) | 2018-04-25 | 2018-04-25 | Regularization Kalman filtering method based on signal-to-noise ratio test |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810379450.4A CN108614804B (en) | 2018-04-25 | 2018-04-25 | Regularization Kalman filtering method based on signal-to-noise ratio test |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108614804A CN108614804A (en) | 2018-10-02 |
CN108614804B true CN108614804B (en) | 2022-09-02 |
Family
ID=63660835
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810379450.4A Expired - Fee Related CN108614804B (en) | 2018-04-25 | 2018-04-25 | Regularization Kalman filtering method based on signal-to-noise ratio test |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108614804B (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102020115919A1 (en) * | 2020-06-17 | 2021-12-23 | Fette Compacting Gmbh | Method for operating a mixing device in a plant |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102064798A (en) * | 2010-12-17 | 2011-05-18 | 北京大学 | Negative-feedback self-adaption on-line and real-time filtering method and system |
WO2012177335A1 (en) * | 2011-06-21 | 2012-12-27 | Exxonmobil Upstream Research Company | Improved dispersion estimation by nonlinear optimization of beam-formed fields |
CN103684349A (en) * | 2013-10-28 | 2014-03-26 | 北京理工大学 | Kalman filtering method based on recursion covariance matrix estimation |
CN103902834A (en) * | 2014-04-14 | 2014-07-02 | 重庆大学 | Structural damage identification method based on ridge estimation and L curve method |
CN104376383A (en) * | 2014-11-27 | 2015-02-25 | 东北大学 | Grid voltage monitoring and prediction system and method based on geographic information system |
CN107425548A (en) * | 2017-09-11 | 2017-12-01 | 河海大学 | A kind of interpolation H ∞ EKFs generator dynamic state estimator method |
CN107729845A (en) * | 2017-10-20 | 2018-02-23 | 开沃新能源汽车集团有限公司 | A kind of frequency respond noise-reduction method decomposed based on sub-space feature value |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
TW200803262A (en) * | 2006-06-23 | 2008-01-01 | Pixart Imaging Inc | Bumping algorithm used to estimate signal frequency and displacement-estimating method and device |
-
2018
- 2018-04-25 CN CN201810379450.4A patent/CN108614804B/en not_active Expired - Fee Related
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102064798A (en) * | 2010-12-17 | 2011-05-18 | 北京大学 | Negative-feedback self-adaption on-line and real-time filtering method and system |
WO2012177335A1 (en) * | 2011-06-21 | 2012-12-27 | Exxonmobil Upstream Research Company | Improved dispersion estimation by nonlinear optimization of beam-formed fields |
CN103684349A (en) * | 2013-10-28 | 2014-03-26 | 北京理工大学 | Kalman filtering method based on recursion covariance matrix estimation |
CN103902834A (en) * | 2014-04-14 | 2014-07-02 | 重庆大学 | Structural damage identification method based on ridge estimation and L curve method |
CN104376383A (en) * | 2014-11-27 | 2015-02-25 | 东北大学 | Grid voltage monitoring and prediction system and method based on geographic information system |
CN107425548A (en) * | 2017-09-11 | 2017-12-01 | 河海大学 | A kind of interpolation H ∞ EKFs generator dynamic state estimator method |
CN107729845A (en) * | 2017-10-20 | 2018-02-23 | 开沃新能源汽车集团有限公司 | A kind of frequency respond noise-reduction method decomposed based on sub-space feature value |
Non-Patent Citations (4)
Title |
---|
A novel GLM-based method for the Automatic IDentification of functional Events (AIDE) in fNIRS data recorded in naturalistic environments;Pinti Paola 等;《Neuroimage》;20170731;第155卷;291-304 * |
Nonparametric hemodynamic deconvolution of fMRI using homomorphic filtering;Sreenivasan K. R. 等;《IEEE Transactions on Medical Imaging》;20141212;第34卷(第5期);1155-1163 * |
一种改进的高斯近似滤波方法;黄玉龙 等;《自动化学报》;20151231;第42卷(第3期);385-401 * |
基于信噪比检验的双参数岭型Kalman滤波及其在BDS星地联合定轨中的应用;李豪 等;《大地测量与地球动力学》;20181115;第38卷(第11期);1143-1148 * |
Also Published As
Publication number | Publication date |
---|---|
CN108614804A (en) | 2018-10-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11245464B2 (en) | Direction-of-arrival estimation and mutual coupling calibration method and system with arbitrary sensor geometry and unknown mutual coupling | |
CN111985093B (en) | Adaptive unscented Kalman filtering state estimation method with noise estimator | |
Brüggemann et al. | Inference in VARs with conditional heteroskedasticity of unknown form | |
CN101251926B (en) | Remote sensing image registration method based on local configuration covariance matrix | |
US10754922B2 (en) | Method and apparatus for sensor fusion | |
CN110061716B (en) | Improved kalman filtering method based on least square and multiple fading factors | |
EP2784537A1 (en) | Inversion method and apparatus based on polarimetric interferometric synthetic aperture radar | |
CN107870001A (en) | A kind of magnetometer bearing calibration based on ellipsoid fitting | |
CN109239653B (en) | Multi-radiation source passive direct time difference positioning method based on subspace decomposition | |
CN110826021B (en) | Robust identification and output estimation method for nonlinear industrial process | |
CN113552565B (en) | Phase unwrapping method for SAR data high noise and large gradient change area | |
Metref et al. | A non-Gaussian analysis scheme using rank histograms for ensemble data assimilation | |
CN117236084B (en) | Intelligent management method and system for woodworking machining production | |
CN110032709B (en) | Positioning and estimation method for abnormal point in geographic coordinate conversion | |
CN104463876B (en) | Adaptive-filtering-based rapid multi-circle detection method for image under complex background | |
CN108614804B (en) | Regularization Kalman filtering method based on signal-to-noise ratio test | |
CN113238227B (en) | Improved least square phase unwrapping method and system combined with deep learning | |
CN111398523B (en) | Sensor data calibration method and system based on distribution | |
CN111142062B (en) | Grid-free target direction-of-arrival estimation method utilizing Toeplitz characteristic | |
CN110672127B (en) | Real-time calibration method for array type MEMS magnetic sensor | |
CN103065021A (en) | multistep effective method for finite element model updating | |
CN113569393B (en) | Robust absolute orientation algorithm | |
US20220308150A1 (en) | Method for direction-of-arrival estimation based on sparse reconstruction in the presence of gain-phase error | |
CN112881971B (en) | Direction finding method for coherent interference source under electromagnetic directional mutual coupling effect | |
Weng et al. | Why might relative fit indices differ between estimators? |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20220902 |
|
CF01 | Termination of patent right due to non-payment of annual fee |