CN117895920B - Distributed consistency Kalman filtering method for sensor network under communication link fault - Google Patents
Distributed consistency Kalman filtering method for sensor network under communication link fault Download PDFInfo
- Publication number
- CN117895920B CN117895920B CN202410282300.7A CN202410282300A CN117895920B CN 117895920 B CN117895920 B CN 117895920B CN 202410282300 A CN202410282300 A CN 202410282300A CN 117895920 B CN117895920 B CN 117895920B
- Authority
- CN
- China
- Prior art keywords
- sensor
- error covariance
- covariance matrix
- estimation
- 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
- 238000001914 filtration Methods 0.000 title claims abstract description 37
- 238000000034 method Methods 0.000 title claims abstract description 37
- 238000004891 communication Methods 0.000 title claims abstract description 29
- 239000011159 matrix material Substances 0.000 claims abstract description 77
- 230000004927 fusion Effects 0.000 claims abstract description 26
- 238000005259 measurement Methods 0.000 claims abstract description 25
- 230000008569 process Effects 0.000 claims description 18
- 230000005540 biological transmission Effects 0.000 claims description 4
- 230000003416 augmentation Effects 0.000 abstract description 2
- 238000005516 engineering process Methods 0.000 abstract description 2
- 238000013461 design Methods 0.000 description 4
- 239000013598 vector Substances 0.000 description 4
- 238000004088 simulation Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000002411 adverse Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 238000012938 design process Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000007499 fusion processing Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
- H03H21/00—Adaptive networks
- H03H21/0012—Digital adaptive filters
- H03H21/0025—Particular filtering methods
- H03H21/0029—Particular filtering methods based on statistics
- H03H21/003—KALMAN filters
-
- 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/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/25—Fusion techniques
- G06F18/251—Fusion techniques of input or preprocessed data
-
- 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)
- Data Mining & Analysis (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Computational Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Bioinformatics & Computational Biology (AREA)
- Software Systems (AREA)
- Evolutionary Biology (AREA)
- Probability & Statistics with Applications (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- Life Sciences & Earth Sciences (AREA)
- Artificial Intelligence (AREA)
- Evolutionary Computation (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Operations Research (AREA)
- Computing Systems (AREA)
- Arrangements For Transmission Of Measured Signals (AREA)
Abstract
The invention discloses a distributed consistency Kalman filtering method of a sensor network under communication link faults, which comprises the following steps: step 1, based on a standard Kalman filtering algorithm, carrying out local measurement update according to measurement information of a sensor network under a communication link fault to obtain local optimal estimation, wherein the local optimal estimation comprises an optimal estimation value and an estimation error covariance matrix of the sensor network; and 2, carrying out information consistency fusion on the local optimal estimation, and updating an optimal estimation value and an estimation error covariance matrix of the sensor network after the consistency fusion is completed. The invention can compensate the influence of the loss of the measured data on the filtering precision, calculates the error covariance matrix and the state estimation variance matrix of the augmentation system by means of matrix theory and random analysis technology, further obtains the updating equation of the distributed filter, can ensure the integral working performance of the whole distributed sensing network, and is suitable for the problem of distributed state estimation when the communication link fails.
Description
Technical Field
The invention belongs to the technical field of signal processing, and particularly relates to a distributed consistency Kalman filtering method of a sensor network under a communication link fault.
Background
The Kalman filter is used as a classical online optimal recursive filter, is widely applied in the fields of target tracking, industrial system state estimation and the like, and is widely focused at home and abroad when being popularized to a distributed sensor network.
It should be noted that, in the actual implementation process, some sensor nodes in the sensor network may fail due to energy exhaustion or environmental interference. Meanwhile, in order to compensate for the failure node and improve the monitoring precision, some new sensor nodes are supplemented into the network. In addition, routing algorithms are commonly employed to reduce power consumption, and these factors can cause changes in the communication links between network sensor nodes. Obviously, the communication link between the sensors is extremely easy to be broken down due to the influence of various complex factors, so that transmission data is lost, distributed resource sharing cannot be realized, the overall performance of the network is further destroyed, filtering precision is reduced, even estimation errors are dispersed, and therefore, the influence of the communication link fault on the distributed filtering performance of the sensor network is required to be studied deeply.
Disclosure of Invention
Aiming at the defects of the prior art, the invention provides a distributed consistency Kalman filtering method of a sensor network under the fault of a communication link, develops the research of a distributed filtering algorithm with robustness to the fault of the communication link, eliminates adverse effects caused by data loss as much as possible, deeply improves the existing filtering method to ensure that the existing filtering method is suitable for the problem of estimating the distributed state when the fault of the communication link exists, particularly compensates the influence of the loss of measured data on the filtering precision by introducing a new online data iterative updating strategy, solves the problem of difficult deduction of a time updating equation and a measuring updating equation, successfully realizes the design of relevant parameters which can meet the robustness requirement to the uncertainty factor of the link, and ensures that the distributed consistency Kalman filter has good performance under the influence of the fault of the communication link, thereby ensuring the whole working performance of the whole distributed sensing network.
In order to achieve the technical purpose, the invention adopts the following technical scheme:
The distributed consistency Kalman filtering method of the sensor network under the fault of the communication link comprises the following steps:
Step 1, based on a standard Kalman filtering algorithm, carrying out local measurement update according to measurement information of a sensor network under a communication link fault to obtain local optimal estimation, wherein the local optimal estimation comprises an optimal estimation value and an estimation error covariance matrix of the sensor network;
And 2, carrying out information consistency fusion on the local optimal estimation, and updating an optimal estimation value and an estimation error covariance matrix of the sensor network after the consistency fusion is completed.
In order to optimize the technical scheme, the specific measures adopted further comprise:
The local measurement update process in the step 1 is as follows:
;
;
wherein, Representing the optimal estimated value of the ith sensor at the k moment; /(I)Representing an estimated error covariance matrix of an ith sensor at a k moment; /(I)Representing the measured value of the ith sensor at time k; /(I)Representing a measurement matrix of an ith sensor; /(I)Representing a one-step predicted value of an ith sensor at a time k; /(I)Representing a one-step prediction error estimation error covariance matrix of an ith sensor at a k moment; /(I)The Kalman filtering gain matrix of the ith sensor at the k moment.
Above mentioned;
Wherein T represents the transpose of the matrix,Representing the estimated error covariance matrix of the noise signal in the ith sensor measurement.
Above mentioned;
Wherein,Representing an estimated value of an ith sensor at time k-1; /(I)Representing an estimated error covariance matrix of an ith sensor at a k-1 moment; /(I)Representing a process noise estimation error covariance matrix in a target system model; a represents a state matrix in the target system model.
The step 2 specifically includes:
setting initial parameters according to local optimal estimation And/>;
Based on initial parametersAnd (3) fusing the information consistency, and after finishing the information consistency fusion, recalculating and updating the optimal estimated value and the estimated error covariance matrix of the sensor network.
First, theThe information consistency fusion formula of the steps is as follows:
;
wherein, Representing the pass/>The inverse matrix of the local optimal estimation error covariance matrix of the ith sensor at the k moment after step fusion; /(I)Equal to/>Multiplying the local optimal estimation of the ith sensor at the k moment; representing the pass/> The inverse matrix of the local optimal estimation error covariance matrix of the ith sensor at the k moment after step fusion; /(I)Equal to/>Multiplying the local optimal estimation of the ith sensor at the k moment; /(I)Representing a set of neighbor nodes of an ith sensor in the sensor network; /(I)Is the weight.
Above mentioned;
Wherein,The weight of direct information transmission between an ith sensor and a jth sensor in the sensor network is represented; Indicating random packet loss of the link.
The optimal estimated value and the estimated error covariance matrix of the updated sensor network are as follows:
;
Wherein, Representing the pass/>The inverse matrix of the local optimal estimation error covariance matrix of the ith sensor at the k moment after step fusion; /(I)Equal to/>Multiplying by the locally optimal estimate of the ith sensor at time k.
The invention has the following beneficial effects:
The distributed consistency filtering algorithm robust to communication link faults provided by the invention firstly compensates the influence of measured data loss on filtering precision by introducing a new online data iterative updating strategy in the structural design process of a filter; secondly, the current state estimation is used to update the estimation error variance matrix in each iteration, which essentially linearizes the system equation and the measurement equation around the state estimation, and in this process, linearization errors are necessarily generated, and also the norm-bounded uncertainty term brought by the logarithmic quantizer is introduced. Therefore, compared with the existing distributed Kalman filtering scheme, the method not only considers the difficulty brought by the introduction of a new compensation strategy to the filter design and analysis, but also further considers the robustness of design parameters to link uncertainty factors. In the specific research process, calculating an optimal estimated value and an estimated error covariance matrix of the augmentation system by means of matrix theory and a random analysis technology, and further obtaining an updating equation of the distributed filter; finally, the theory and algorithm are further perfected by utilizing numerical simulation, and an evaluation result is given.
Drawings
FIG. 1 is a flow chart of the method of the present invention;
FIG. 2 is a diagram of a sensor network topology;
fig. 3 is a graph of the mean square estimation error of each sensor.
Detailed Description
The present invention will be described in further detail with reference to the drawings and examples, in order to make the objects, technical solutions and advantages of the present invention more apparent. It should be understood that the specific embodiments described herein are for purposes of illustration only and are not intended to limit the scope of the invention.
Although the steps of the present invention are arranged by reference numerals, the order of the steps is not limited, and the relative order of the steps may be adjusted unless the order of the steps is explicitly stated or the execution of a step requires other steps as a basis. It is to be understood that the term "and/or" as used herein relates to and encompasses any and all possible combinations of one or more of the associated listed items.
As shown in fig. 1, the distributed consistency kalman filtering method of the sensor network under the fault of the communication link of the invention comprises the following steps:
Step 1, based on a standard Kalman filtering algorithm, carrying out local measurement update according to measurement information of a sensor network under a communication link fault to obtain local optimal estimation, wherein the local optimal estimation comprises an optimal estimation value and an estimation error covariance matrix of the sensor network;
And 2, carrying out information consistency fusion on the local optimal estimation, and updating an optimal estimation value and an estimation error covariance matrix of the sensor network after the consistency fusion is completed.
When the algorithm of the invention is used for filtering, the filtering result is as follows: the estimated values of all the sensors can accurately track the state of the target system after distributed fusion. And obtaining the state estimation of the sensor network.
The covariance matrix in the algorithm refers to covariance of the estimation error of the sensor network, and because the estimation error is a random quantity, statistical information, covariance, is needed, and because the estimation error is a vector, the covariance matrix is needed. The covariance matrix is used for completing the updating of the subsequent estimated value and is also an evaluation index of the estimation accuracy.
In an embodiment, the specific implementation process of the invention comprises the following steps: the sensor obtains local optimal estimation based on a standard Kalman filtering algorithm according to self measurement information; setting iteration initial parameters; carrying out local information fusion in an information consistency mode; after the consistency is completed, the updated estimation error covariance matrix and the estimation equation are recalculated and recorded as the estimation error covariance matrix and the estimation value of the sensor at the moment. The finite nature of the estimation error at finite coherence steps and infinite coherence steps is demonstrated. The concrete introduction is as follows:
Firstly, constructing and analyzing a sensor network architecture model:
Consider the following type of linear system described by a state space model: ;
Wherein the method comprises the steps of Is a state vector of the system. First/>The measurement model of the individual sensor nodes is described as follows:
;
Wherein the method comprises the steps of Is/>Individual sensor at/>Measurement of time of day, disturbance term/>And/>The zero-mean and mutually independent Gaussian white noise has the estimated error covariance matrixes of/>, respectivelyAnd/>. Initial state/>For covariance matrix/>, with zero mean and estimation errorAnd/or > Gaussian noiseAnd/>Independent of each other; /(I)And/>Is a known matrix with the appropriate dimensions.
In fact, directed graphsA communication topology of a network may be described in whichFor node set,/>For edge set,/>Is a weighted adjacency matrix. Edge(s)Represents the/>The individual node may receive the data from the/>Information of individual nodes. Weighted adjacency matrix/>Is a row random matrix (i.e., the sum of each row of elements is 1), where each element/>Non-negative; if/>Then/>Otherwise. In addition, since the node can obtain its own information without communication, it/>. At this time, considering that the network node and the link mode may change randomly in a complex environment, a markov process is used to describe the communication topology change of the sensor network. Graph/>Evolution of (1) is defined by a homogeneous Markov process/>Decision, the Markov process is in a finite setAnd (3) taking the value. Accordingly, we can get the node set/>Edge set/>Weighted adjacency matrix/>And contiguous node set/>The transition probability matrix of the markov process is/>。
Then on the basis, a distributed consistency Kalman filtering algorithm meeting the Markov switching topological structure is designed, and the accurate and reliable estimation of the target state is completed under the condition of dynamic topology. The design of the distributed consistency Kalman filtering algorithm is specifically as follows:
each sensor node adopts a Kalman distributed filtering algorithm to estimate the state of the target system.
Firstly, the sensor obtains a local optimal estimation based on a standard Kalman filtering algorithm according to self measurement information. The one-step state prediction process of the Kalman filtering algorithm is as follows:
;
wherein, Representing an estimated value of an ith sensor at time k-1; /(I)Representing an estimated error covariance matrix of an ith sensor at a k-1 moment; /(I)Representing a process noise estimation error covariance matrix in a target system model; a represents a state matrix in the target system model.
The local measurement update procedure is as follows:
;
;
wherein, Representing the optimal estimated value of the ith sensor at the k moment; /(I)Representing an estimated error covariance matrix of an ith sensor at a k moment; /(I)Representing the measured value of the ith sensor at time k; /(I)Representing a measurement matrix of an ith sensor; /(I)Representing a one-step predicted value of an ith sensor at a time k; /(I)Representing a one-step prediction error estimation error covariance matrix of an ith sensor at a k moment; /(I)The Kalman filtering gain matrix of the ith sensor at the k moment.
The Kalman filter gain matrix is given by:
;
Wherein T represents the transpose of the matrix, Representing the estimated error covariance matrix of the noise signal in the ith sensor measurement.
The collaborative estimation is given using the distributed consensus kalman filter algorithm as shown in table 1:
;
In Table 1, step 3, Representing the pass/>The inverse matrix of the local optimal estimation error covariance matrix of the ith sensor at the k moment after step fusion; /(I)Equal to/>Multiplying the local optimal estimation of the ith sensor at the k moment; /(I)Representing the pass/>The inverse matrix of the local optimal estimation error covariance matrix of the ith sensor at the k moment after step fusion; /(I)Equal to/>Multiplying the local optimal estimation of the ith sensor at the k moment; /(I)Representing a set of neighbor nodes of an ith sensor in the sensor network; /(I)Is the weight.
In Table 1, step 4,Representing the pass/>The inverse matrix of the local optimal estimation error covariance matrix of the ith sensor at the k moment after step fusion; /(I)Equal to/>Multiplying by the locally optimal estimate of the ith sensor at time k.
In table 1, step 3, the adjacent sensing nodes need to implement information interaction through a network and acquire information vectors and information matrixes, that is, the sensors need to mutually transmit the information vectors and the information matrixes with the adjacent sensing nodes through the network, so that local information fusion estimation is completed. The process is carried out in an information consistency mode, L steps of consistency coordination steps are shared in consideration of time k, and information is required to be transmitted outwards in each step of sensor network. It is not difficult to find that the consistency process information transfer at this time depends on the real-time topology of the sensor network.
In Table 1, step 3 thus completes the local information fusion process. That is, the consistency process is affected by random packet loss during communication, and thus the weight is defined as follows:
;
wherein, The weight of direct information transmission between an ith sensor and a jth sensor in the sensor network is represented; Indicating random packet loss of the link, which obeys an independent co-distributed Bernoulli distribution. Specifically, at time k/> In the step of consistent filtering, nodes/>And/>Communication betweenAnd if not, packet loss occurs. For each communication link, it is assumed that the packet loss rate is less than 1, namely:
;
Wherein the method comprises the steps of For all/>And/>All are true. Because the node can obtain its own information regardless of the network conditions, there is/>,/>. Let us assume random variable/>Independent of process noise/>Measurement noise/>And initial state/>. Thus far, a distributed consistency Kalman filter algorithm under the condition of link failure is given.
Finally, deep performance analysis is carried out on the algorithm, and the convergence condition of the algorithm is discussed.
The convergence analysis of the distributed consistency kalman filter algorithm is as follows:
To analyze convergence, the algorithm proposed by the present invention will be described as having consistency. "compatibility" here means that the resulting covariance of the algorithm must be the upper bound of the system's true covariance matrix. First, the true one-step prediction error covariance is defined as follows:
;
The true estimation error covariance is as follows:
;
the true estimated error covariance after multi-step consistency synergy/fusion is as follows:
;
It is easy to verify that the consistent kalman filtering algorithm under the markov switching topology is coherent, i.e.:
;
The sensor network is assumed to be globally observable, i.e Considerable, where/>;
In applications where sensor nodes are limited in functionality, it is difficult for a single node to be able to perform the task of target tracking, and the assumption of local observability is often not reasonable. As a whole, the sensor network comprises different types of sensors, and the defects of the performance of a single sensor can be generally overcome through cooperation among the sensors, so that global observability is realized, and the assumption is reasonable.
Under globally observable conditions/concepts, the following results can be demonstrated:
(1) Estimation error bounded analysis under finite consistency step size: if the sensing network is globally considerable and the link fault probability is not constant 1, then the estimated error covariance is randomly bounded, namely:
;
(2) Estimation error bounded analysis under infinite consistency step size: if the sensing network is globally considerable, the link fault probability is not constant 1, and the number of the consistency steps is large enough, a positive scalar exists So that the following consistent bounded condition is established,
;
For a pair ofAlmost certainly.
Therefore, the system convergence is influenced by link faults, and has deep requirements on the network topology structure in addition to the requirements on the observation performance of the sensors in the sensor network. If the topology structure cannot meet the condition of communication in the mean value sense, the network part nodes are likely to be isolated nodes, and the information of the network part nodes cannot be shared globally, so that the task fails.
In order that the practice of the invention may be better understood by those skilled in the art, the invention is simulated using Matlab software.
The linear system model considered is as follows:
;
Wherein the method comprises the steps of For sampling period,/>Zero mean and variance of/>Is a gaussian distribution of (c). Let us assume initial state/>Satisfying a mean of 0 variance as/>Is a gaussian distribution of (c). For convenience, for all/>Weights are all set to/>. The probability of no packet loss of the communication link is given to be 10%. The topology of the hypothetical network is already given.
The individual sensor measurement matrices are as follows:
;
Wherein the method comprises the steps of Obeying the mean to be 0 variance/>Is a gaussian distribution of (c).
The network is not difficult to verify to meet global observability, and the mean square estimation error (MSE) obtained by the proposed distributed consistency filtering algorithm can be verified to have consistency through the simulation result graph, and the overall requirement of error convergence is met. In fact, for a single sensor, the estimation error curve will diverge because of not satisfying observability itself, but the estimation error remains converged as a whole for the sensor network, thus it can be seen that the performance of the present invention is significantly better than that of a single sensor algorithm.
Fig. 2 is a network topology structure diagram arbitrarily given, and fig. 3 is a simulation result of mean square estimation error (MSE) of each sensor, from which MSE can be found to have consistency and successfully meet the overall requirement of error convergence, which verifies the effectiveness of the present invention.
It will be evident to those skilled in the art that the invention is not limited to the details of the foregoing illustrative embodiments, and that the present invention may be embodied in other specific forms without departing from the spirit or essential characteristics thereof. The present embodiments are, therefore, to be considered in all respects as illustrative and not restrictive, the scope of the invention being indicated by the appended claims rather than by the foregoing description, and all changes which come within the meaning and range of equivalency of the claims are therefore intended to be embraced therein. Any reference sign in a claim should not be construed as limiting the claim concerned.
Furthermore, it should be understood that although the present disclosure describes embodiments, not every embodiment is provided with a separate embodiment, and that this description is provided for clarity only, and that the disclosure is not limited to the embodiments described in detail below, and that the embodiments described in the examples may be combined as appropriate to form other embodiments that will be apparent to those skilled in the art.
Claims (2)
1. The distributed consistency Kalman filtering method of the sensor network under the communication link fault is characterized by comprising the following steps:
Step 1, based on a standard Kalman filtering algorithm, carrying out local measurement update according to measurement information of a sensor network under a communication link fault to obtain local optimal estimation, wherein the local optimal estimation comprises an optimal estimation value and an estimation error covariance matrix of the sensor network;
The local measurement updating process comprises the following steps:
wherein, Representing the optimal estimated value of the ith sensor at the k moment; /(I)Representing an estimated error covariance matrix of an ith sensor at a k moment; y k,i represents the measurement of the ith sensor at time k; c i represents the measurement matrix of the ith sensor; /(I)Representing a one-step predicted value of an ith sensor at a time k; /(I)Representing a one-step prediction error estimation error covariance matrix of an ith sensor at a k moment; /(I)A Kalman filtering gain matrix of the ith sensor at the k moment;
Wherein T represents a matrix transpose, and R i represents an estimated error covariance matrix of the noise signal in the ith sensor measurement;
wherein, Representing an estimated value of an ith sensor at time k-1; /(I)Representing an estimated error covariance matrix of an ith sensor at a k-1 moment; q represents a process noise estimation error covariance matrix in a target system model; a represents a state matrix in a target system model;
Step 2, information consistency fusion is carried out on the local optimal estimation, and after consistency fusion is completed, an optimal estimation value and an estimation error covariance matrix of the sensor network are updated;
The step 2 specifically includes:
setting initial parameters according to local optimal estimation And/>
Carrying out L-step information consistency fusion based on initial parameters, and after finishing the information consistency fusion, recalculating and updating an optimal estimated value and an estimated error covariance matrix of the sensor network;
the information consistency fusion formula of the step L is as follows:
wherein, Representing the inverse matrix of the local optimal estimation error covariance matrix of the ith sensor at the k moment after l+1 steps of fusion; /(I)Equal to/>Multiplying the local optimal estimation of the ith sensor at the k moment; /(I)Representing the inverse matrix of the local optimal estimation error covariance matrix of the ith sensor at the k moment after l-step fusion; /(I)Equal toMultiplying the local optimal estimation of the ith sensor at the k moment; n i represents a set of neighbor nodes of an ith sensor in the sensor network; /(I)The formula is as follows:
wherein pi i,j represents the weight of direct information transmission between the ith sensor and the jth sensor in the sensor network; Indicating random packet loss of the link.
2. The distributed consistency kalman filtering method of sensor network under communication link failure according to claim 1, wherein the updated optimal estimated value and estimated error covariance matrix of the sensor network are:
wherein, Representing the inverse matrix of the local optimal estimation error covariance matrix of the ith sensor at the k moment after the L-step fusion; /(I)Equal to/>Multiplying by the locally optimal estimate of the ith sensor at time k.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410282300.7A CN117895920B (en) | 2024-03-13 | 2024-03-13 | Distributed consistency Kalman filtering method for sensor network under communication link fault |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410282300.7A CN117895920B (en) | 2024-03-13 | 2024-03-13 | Distributed consistency Kalman filtering method for sensor network under communication link fault |
Publications (2)
Publication Number | Publication Date |
---|---|
CN117895920A CN117895920A (en) | 2024-04-16 |
CN117895920B true CN117895920B (en) | 2024-05-17 |
Family
ID=90642594
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202410282300.7A Active CN117895920B (en) | 2024-03-13 | 2024-03-13 | Distributed consistency Kalman filtering method for sensor network under communication link fault |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117895920B (en) |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106878375A (en) * | 2016-12-22 | 2017-06-20 | 中国民航大学 | A kind of cockpit pollutant monitoring method based on distribution combination sensor network |
CN108427656A (en) * | 2018-01-03 | 2018-08-21 | 东南大学 | A kind of multiple target encirclement control method based on matrix weights consistency-Kalman filter |
CN108645415A (en) * | 2018-08-03 | 2018-10-12 | 上海海事大学 | A kind of ship track prediction technique |
CN111601269A (en) * | 2020-05-15 | 2020-08-28 | 中国民航大学 | Event trigger Kalman consistency filtering method based on information freshness judgment |
CN114024524A (en) * | 2021-11-01 | 2022-02-08 | 北京科技大学 | Distributed asynchronous consistent filtering method, electronic device and computer storage medium |
CN115603845A (en) * | 2022-09-19 | 2023-01-13 | 安徽理工大学(Cn) | Optimal distributed filter cooperative target tracking method based on consistency strategy |
CN116186643A (en) * | 2023-04-27 | 2023-05-30 | 北京航空航天大学 | Multi-sensor collaborative target tracking method, system, equipment and medium |
CN116699598A (en) * | 2023-07-10 | 2023-09-05 | 哈尔滨工业大学 | Multi-sensor target tracking method, system and equipment based on consistency filtering |
-
2024
- 2024-03-13 CN CN202410282300.7A patent/CN117895920B/en active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106878375A (en) * | 2016-12-22 | 2017-06-20 | 中国民航大学 | A kind of cockpit pollutant monitoring method based on distribution combination sensor network |
CN108427656A (en) * | 2018-01-03 | 2018-08-21 | 东南大学 | A kind of multiple target encirclement control method based on matrix weights consistency-Kalman filter |
CN108645415A (en) * | 2018-08-03 | 2018-10-12 | 上海海事大学 | A kind of ship track prediction technique |
CN111601269A (en) * | 2020-05-15 | 2020-08-28 | 中国民航大学 | Event trigger Kalman consistency filtering method based on information freshness judgment |
CN114024524A (en) * | 2021-11-01 | 2022-02-08 | 北京科技大学 | Distributed asynchronous consistent filtering method, electronic device and computer storage medium |
CN115603845A (en) * | 2022-09-19 | 2023-01-13 | 安徽理工大学(Cn) | Optimal distributed filter cooperative target tracking method based on consistency strategy |
CN116186643A (en) * | 2023-04-27 | 2023-05-30 | 北京航空航天大学 | Multi-sensor collaborative target tracking method, system, equipment and medium |
CN116699598A (en) * | 2023-07-10 | 2023-09-05 | 哈尔滨工业大学 | Multi-sensor target tracking method, system and equipment based on consistency filtering |
Also Published As
Publication number | Publication date |
---|---|
CN117895920A (en) | 2024-04-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111178385B (en) | Target tracking method for robust online multi-sensor fusion | |
CN108896047B (en) | Distributed sensor network collaborative fusion and sensor position correction method | |
CN105740203A (en) | Multi-sensor passive synergic direction finding and positioning method | |
CN110958639A (en) | Target state estimation method and system | |
CN110011879B (en) | Sensor network safety real-time online monitoring system based on parallel filtering | |
CN116186643B (en) | Multi-sensor collaborative target tracking method, system, equipment and medium | |
CN116088303B (en) | Uncertain complex dynamic network state time-varying recursion estimation method | |
CN108347694B (en) | Node positioning method and system based on boundary conditions | |
Jin et al. | Distributed Kalman filters with random sensor activation and noisy channels | |
CN103313386B (en) | Based on the radio sensing network method for tracking target of consistency on messaging right-value optimization | |
CN104535963B (en) | Cooperative target positioning achievement method of multiple mobile nodes based on distance measurement | |
De Silva et al. | Efficient distributed multi-robot localization: A target tracking inspired design | |
CN117895920B (en) | Distributed consistency Kalman filtering method for sensor network under communication link fault | |
CN103313384B (en) | A kind of radio sensing network method for tracking target based on consistency on messaging | |
CN109362049B (en) | Square root volume filtering method based on mixed information weighting consistency | |
CN111883265A (en) | Target state estimation method applied to fire control system | |
CN107966697B (en) | Moving target tracking method based on progressive unscented Kalman | |
CN108445445B (en) | Distributed passive positioning method and device | |
CN113110321B (en) | Distributed estimation method of networked industrial control system based on event trigger | |
CN109474892B (en) | Strong robust sensor network target tracking method based on information form | |
CN115859030B (en) | Two-step state estimation method under complex coupling | |
CN110807478A (en) | Cooperative target tracking method under condition of observing intermittent loss | |
Marinakis et al. | Hybrid inference for sensor network localization using a mobile robot | |
Wang et al. | Multi-rate data fusion for wireless sensor networks with time-delay based on improved cubature kalman filter | |
CN114430525B (en) | Sensor network-based distributed positioning method for closed space |
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 |