CN109039725A - It is a kind of with the complex network optimal estimating method that couples at random - Google Patents
It is a kind of with the complex network optimal estimating method that couples at random Download PDFInfo
- Publication number
- CN109039725A CN109039725A CN201810814462.5A CN201810814462A CN109039725A CN 109039725 A CN109039725 A CN 109039725A CN 201810814462 A CN201810814462 A CN 201810814462A CN 109039725 A CN109039725 A CN 109039725A
- Authority
- CN
- China
- Prior art keywords
- node
- complex network
- time
- matrix
- ith node
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 55
- 238000010168 coupling process Methods 0.000 claims abstract description 43
- 238000005859 coupling reaction Methods 0.000 claims abstract description 43
- 238000005259 measurement Methods 0.000 claims abstract description 42
- 230000008878 coupling Effects 0.000 claims abstract description 41
- 239000011159 matrix material Substances 0.000 claims description 80
- 238000005457 optimization Methods 0.000 claims description 7
- 238000012512 characterization method Methods 0.000 claims description 3
- 230000017105 transposition Effects 0.000 claims 1
- 230000005540 biological transmission Effects 0.000 abstract description 9
- 230000001960 triggered effect Effects 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 5
- 238000013461 design Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L41/00—Arrangements for maintenance, administration or management of data switching networks, e.g. of packet switching networks
- H04L41/14—Network analysis or design
- H04L41/145—Network analysis or design involving simulating, designing, planning or modelling of a network
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L41/00—Arrangements for maintenance, administration or management of data switching networks, e.g. of packet switching networks
- H04L41/14—Network analysis or design
- H04L41/142—Network analysis or design using statistical or mathematical methods
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Physics & Mathematics (AREA)
- Algebra (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Probability & Statistics with Applications (AREA)
- Pure & Applied Mathematics (AREA)
- Data Exchanges In Wide-Area Networks (AREA)
Abstract
It is a kind of with the complex network optimal estimating method that couples at random, the present invention relates to Random Coupling complex network optimal estimating method.The present invention solves standing state estimation method cannot handle the time lag complex network coupled at random with measurement Loss and inaccurate probability of happening simultaneously, cause to estimate that performance accuracy rate is low, and there are loss of transmitted data, transmission failure, switching nodes can not receive in the simultaneous situation of other nodal informations, leads to the problem for estimating that performance accuracy rate is low.Process are as follows: one, the random generation coupling time lag complex network dynamic model of foundation;Two, state estimation is carried out to dynamic model under event-triggered protocols;Three, ∑ is calculatedi,k+1|k;Four, K is calculatedi,k+1;Five, it obtainsJudge whether k+1 reaches M, if k+1 < M, executes six, otherwise terminate;Six, ∑ is calculatedi,k+1|k+1;Another k=k+1 executes two, until meeting k+1=M.The present invention is used for complex network optimal estimating field.
Description
Technical Field
The invention relates to a random coupling complex network optimization estimation method.
Background
State estimation of a complex network is an important research problem in a control system, and is widely applied to signal estimation tasks in the fields of engineering, power grids, social networks and the like.
When the network is congested, the measurement loss phenomenon often occurs, and in practical application, the coupling between nodes occurs randomly due to network switching. Therefore, it is necessary to design a state estimation method suitable for these network-induced phenomena, especially when the probability of random coupling is uncertain;
the existing state estimation method can not process the randomly coupled time-lag complex network with measurement loss and inaccurate occurrence probability at the same time, so that the estimation performance accuracy is low;
the existing state estimation method has the problems of low estimation performance accuracy rate caused by the loss of transmission data, transmission failure and incapability of receiving information of other nodes by a coupling node.
Disclosure of Invention
The invention solves the problems that the existing state estimation method can not simultaneously process the time-lag complex network with random coupling and measurement loss and inaccurate occurrence probability, which causes low estimation performance accuracy, and the estimation performance accuracy is low under the condition that transmission data is lost, transmission fails, and coupling nodes can not receive other node information simultaneously, and provides the complex network optimization estimation method with random coupling.
A complex network optimization estimation method with random coupling comprises the following specific processes:
the complex network can be a network formed by satellites, a network formed by robots, a network formed by spacecrafts or a network formed by radars;
establishing a random coupling time-lag complex network dynamic model with measurement loss and inaccurate occurrence probability;
secondly, performing state estimation on the random coupling time-lag complex network dynamic model with the measurement loss phenomenon and the inaccurate occurrence probability established in the first step under an event trigger protocol;
step three, calculating the upper bound sigma of the one-step prediction error covariance matrix of the state estimation of each node of the complex networki,k+1|k;
Step four, according to the one-step prediction error covariance matrix upper bound sigma obtained in step threei,k+1|kCalculating an estimated gain moment for each node of the complex networkMatrix Ki,k+1;
Step five, obtaining the estimated gain matrix K of each node in the step fouri,k+1Substituting the state estimation formula 8 in the step two to obtain the state estimation of the ith node at the (k + 1) th momentThe state estimation of the randomly generated coupling time-lag complex network with the measurement loss phenomenon and the inaccurate generation probability is realized.
Judging whether k +1 reaches the total time length M of the complex network, if k +1 is less than M, executing a sixth step, and if k +1 is equal to M, ending the operation;
step six, according to the estimated gain matrix K of each node of the complex network calculated in the step fouri,k+1Calculating the upper bound sigma of the covariance matrix of the estimation error of each node of the complex networki,k+1|k+1;
And k is k +1, and the step two is executed until k +1 is M.
The invention has the following effects:
the invention provides a time-lag complex network state estimation method with random coupling of measurement loss phenomenon and inaccurate occurrence probability, which considers the influence of the random coupling of the measurement loss phenomenon and the inaccurate occurrence probability on state estimation performance, and utilizes an extended Kalman filtering method to comprehensively consider effective information of an estimation error covariance matrix; the method solves the problems that the existing state estimation method can not simultaneously process the randomly coupled time-lag complex network with the measurement loss phenomenon and the inaccurate occurrence probability, so that the estimation performance accuracy is low, and the estimation performance accuracy is improved.
The method utilizes an extended Kalman filtering method, obtains an estimated error covariance matrix by considering effective information of the estimated error covariance matrix, and then ensures that the trace of the estimated error covariance matrix can obtain the minimum value in each step by designing a gain matrix. A minimization of the estimation error is guaranteed. The performance estimation is not affected under the condition that the transmission data is lost, the transmission fails and the condition that the existing coupling node cannot receive other node information (coupling occurs randomly) occurs simultaneously, and the estimation accuracy is improved.
Solves the problem of low accuracy of estimation performance caused by the existing method under the condition that the transmission data is lost, the transmission fails and the coupling node can not receive the information of other nodes at the same time,
the relative error at all times of all nodes in case one is 24.82%, and the relative error at all times of all nodes in case two is 38.48% with reference to the figure. It can be seen that as the trigger threshold increases, the state estimation performance of the network gradually decreases.
Drawings
FIG. 1 is a flow chart of the method of the present invention;
FIG. 2a is a diagram of the actual state trace of the first node of the complex networkState estimation trajectory at two different trigger thresholdsIn the comparison of the figures of (a),a first component of a state variable at a kth time for a 1 st node of the complex network; whereinIs the trajectory of the state of the system,is the state estimation trajectory for the case that,is the state estimation trajectory in case two;
FIG. 2b is a diagram of the actual state trace of the first node of the complex networkState estimation trajectory at two different trigger thresholdsIn the comparison of the figures of (a),a second component of the state variable at the kth time for the 1 st node of the complex network;
FIG. 3a is a diagram of the actual state trace of a second node of a complex networkState estimation trajectory at two different trigger thresholdsIn the comparison of the figures of (a),a first component of a state variable of a 2 nd node of the complex network at a kth time;
FIG. 3b is a diagram of the actual state trace of a second node of a complex networkState estimation trajectory at two different trigger thresholdsIn the comparison of the figures of (a),a second component of the state variable at the kth time of the 2 nd node of the complex network;
FIG. 4a is a diagram of the actual state trace of the third node of the complex networkState estimation trajectory at two different trigger thresholdsIn the comparison of the figures of (a),a first component of a state variable of a 3 rd node of the complex network at a kth time;
FIG. 4b is the track of the actual state of the third node of the complex networkState estimation trajectory at two different trigger thresholdsIn the comparison of the figures of (a),a second component of the state variable at the kth time of the 3 rd node of the complex network;
FIG. 5a is a trace plot of the upper bound trace of the state estimation error covariance matrix of the first node at different measurement loss probabilities1,k|kUpper bound of covariance matrix of estimation error of 1 st node at k-th time2,k|kUpper bound of covariance matrix of estimation error of 2 nd node at k-th time, sigma3,k|kUpper bound, tra, of the covariance matrix of the estimation error of the 3 rd node at the k-th time instantce(Σ1,k|k) Is sigma1,k|kTrace of (d, trace (sigma)2,k|k) Is sigma2,k|kTrace of (d, trace (sigma)3,k|k) Is sigma3,k|kStep is iteration times, and k is time; whereinIs thatA trace map of the trace that is under,is thatA trace map of the trace that is under,is thatA trace map of the trace that is under,is thatA trace map of the traces;
FIG. 5b is a trace plot of the upper bound trace of the state estimation error covariance matrix of the second node at different measurement loss probabilities;
FIG. 5c is a trace plot of the upper bound of the state estimation error covariance matrix of the third node at different measurement loss probabilities.
Detailed Description
The first embodiment is as follows: the embodiment is described with reference to fig. 1, and a specific process of the complex network optimization estimation method with random occurrence coupling of the embodiment is as follows:
the complex network can be a network formed by satellites, a network formed by robots, a network formed by spacecrafts or a network formed by radars;
establishing a random coupling time-lag complex network dynamic model with measurement loss and inaccurate occurrence probability;
secondly, performing state estimation on the random coupling time-lag complex network dynamic model with the measurement loss phenomenon and the inaccurate occurrence probability established in the first step under an event trigger protocol;
step three, calculating the upper bound sigma of the one-step prediction error covariance matrix of the state estimation of each node of the complex networki,k+1|k;
Step four, according to the one-step prediction error covariance matrix upper bound sigma obtained in step threei,k+1|kCalculating an estimated gain matrix K for each node of the complex networki,k+1;
Step five, obtaining the estimated gain matrix K of each node in the step fouri,k+1Substituting the state estimation formula (8) in the step two to obtain the state estimation of the ith node at the (k + 1) th momentThe state estimation of the randomly generated coupling time-lag complex network with the measurement loss phenomenon and the inaccurate generation probability is realized.
Judging whether k +1 reaches the total time length M of the complex network, if k +1 is less than M, executing a sixth step, and if k +1 is equal to M, ending the operation;
step six, according to the estimated gain matrix K of each node of the complex network calculated in the step fouri,k+1Calculating the upper bound sigma of the covariance matrix of the estimation error of each node of the complex networki,k+1|k+1;
And k is k +1, and the step two is executed until k +1 is M.
The second embodiment is as follows: the first difference between the present embodiment and the specific embodiment is: establishing the measurement loss phenomenon (lambda) in the step onei,k) And probability of occurrence of inaccuracy (gamma)i,k) The random coupling time lag complex network dynamic model;
the state space form of the random coupling time-lag complex network dynamic model with the measurement loss phenomenon and the inaccurate occurrence probability is as follows:
yi,k=λi,kCi,kxi,k+νi,k,i=1,2,…,N, (2)
in the formula,
respectively representing the state variables of the ith node of the complex network at the k, k +1 and k-d moments; d is a fixed network time lag;is a real number domain of the state of the complex network dynamic model, and n is a dimension; x is the number ofj,kThe state variable of the jth node of the complex network at the kth moment;the measurement output of the ith node at the kth moment;a real number domain is output by a complex network dynamic model, and p is a dimension;is the initial value of the ith node at time k, k-d, -d +1, …, 0; Γ is a known connection matrix; w ═ Wij]N×NIs a known coupling matrix, wijCoupling weights for the ith node and the jth node;is a mean of zero and a variance of Qi,kThe process noise of (a) is generated,a real number domain of the process noise of the complex network dynamic model is formed, and q is a dimension;is a mean of zero and a variance of Ri,kThe measurement noise of (2); n is the number of nodes of the complex network; a. thei,kIs a system matrix of known time k, Bi,kIs a known noise distribution matrix at time k, Ci,kFor a known measurement matrix at time k,a system matrix of known dimensions and associated with a time lag k-d; gamma rayi,kAnd λi,kRandom variables obeying Bernoulli distribution respectively depict random coupling phenomenon and measurement loss phenomenon, and satisfy the following conditions:
wherein,
γi,k+1、λi,k+1random variables obeying Bernoulli distribution respectively depict random coupling phenomenon and measurement loss phenomenon,is a known constant, which is the expected probability of random coupling of the ith node at the time k + 1;the constant is a known constant and is the measurement loss probability of the ith node at the moment k + 1;andhas a value range of (0,1), Δ γi,k=1Inaccuracy of characterization probability; prob { } is a probability,as desired;
the above-mentionedThe inaccuracy of the characterization probability is such that,is an upper bound of probability inaccuracy, is a known constant, and satisfies
Other steps and parameters are the same as those in the first embodiment.
The third concrete implementation mode: the present embodiment differs from the first or second embodiment in that: performing state estimation on the random coupling time-lag complex network dynamic model with the measurement loss phenomenon and the inaccurate occurrence probability established in the step one under the event triggering protocol in the step two; the specific process is as follows:
firstly, aiming at the ith node, selecting the following event trigger formula:
whereinFor the measurement output of the ith node at the last trigger time,for the corresponding last value of the triggering time, δi0 is a known adjustment threshold, T is transposed; the next event trigger sequence for the ith nodeIteratively generated by:
whereinIs a positive integer set, inf { } is a lower limit taking function;
after the event triggering mechanism, the true measurement value delivered to the filter is
When the next event trigger sequence does not arrive, the real measurement value is always the value of the previous trigger time;
and aiming at the ith node of the complex network, constructing a state estimator:
in the formula,
is xi,kOne-step prediction at the time k (the specific process of one-step prediction is formula 7),for the state estimate of the ith node at time k +1,for the state estimation of the ith node at the K-d time, Ki,k+1For the estimated gain matrix at time k +1 for the ith node,for the state estimate of the ith node at the kth time,for the state estimate of the jth node at the kth instant,for the measurement output of the ith node at time k +1, Ci,k+1A measurement matrix at a known time k + 1;the expected probability of random coupling at time k for the ith node,is a known constant.
Other steps and parameters are the same as those in the first or second embodiment.
The fourth concrete implementation mode: the difference between this embodiment mode and one of the first to third embodiment modes is: computing a one-step prediction error covariance matrix upper bound sigma of a state estimate for each node of the complex network in step threei,k+1|k(ii) a The specific process is as follows:
calculating the upper bound sigma of the one-step prediction error covariance matrix of the state estimation of each node of the complex network according to the following formulai,k+1|k:
In the formula, τ1、τ2、τ4、τ5Is the intermediate variable(s) of the variable,
in the formula, τ1=1+ε1+ε2+ε3,
εsFor weighting factors greater than zero, s is 1,2, …,7 (for ∑ isi,k+1|kTrace of (d) is minimal);is epsilons1,2, …, 7;
wiis the intermediate variable(s) of the variable,Σi,k+1|kone-step prediction error covariance matrix upper bound, Σ, for the ith node at time ki,k|kThe covariance matrix of the estimation error at the kth time for the ith node is bounded,are respectively Ai,k,Γ,Bi,k,Transposing; sigmai,k-d|k-dUpper bound of covariance matrix of estimation error of ith node at k-dj,k|kAnd (4) the estimation error covariance matrix of the jth node at the kth moment is bounded.
Other steps and parameters are the same as those in one of the first to third embodiments.
The fifth concrete implementation mode: the difference between this embodiment and one of the first to fourth embodiments is: the fourth step is to obtain the upper bound sigma of the one-step prediction error covariance matrix according to the third stepi,k+1|kCalculating an estimated gain matrix K for each node of the complex networki,k+1(ii) a The specific process is as follows:
calculating an estimated gain matrix K of each node of the complex network according to the formula (9) as followsi,k+1:
In the formula,
is the intermediate variable(s) of the variable,ηlis a constant greater than zero (in order to make Σ in step six)i,k+1|k+1Trace of (d) is minimal), l ═ 1,2,3, 4;is etalThe inverse of (a) is,is λi,k+1The variance of (a) is determined,Ki,k+1for the estimated gain matrix of the ith node at time k +1,is Ci,k+1I is an identity matrix, Ri,k+1V isi,k+1The variance of (c) is.
Other steps and parameters are the same as in one of the first to fourth embodiments.
The sixth specific implementation mode: the difference between this embodiment and one of the first to fifth embodiments is: in the sixth step, the estimated gain matrix K of each node of the complex network calculated in the fourth stepi,k+1Calculating the upper bound sigma of the covariance matrix of the estimation error of each node of the complex networki,k+1|k+1(ii) a The specific process is as follows:
the formula is as follows:
wherein,
Σi,k+1|k+1estimate for the ith node at time k +1The upper bound of the covariance matrix of the error is counted,andare respectively asAnd Ki,k+1The transposing of (1).
Furthermore, it can be shown that with probability, other conditions are not changedIncrease of (1), estimation error covariance matrix upper bound ∑i,k+1|k+1The trace of (a) is not increased. Sigmai,k+1|k+1Carry over to step three (replace Σ in equation 9)i,k|k)。
The theory in the third step, the fourth step and the fifth step is as follows:
the minimum upper bound of the estimation error covariance matrix of each node is calculated. Find sigmai,k+1|k+1So that P isi,k+1|k+1≤∑i,k+1|k+1WhereinThe covariance matrix of the estimated error at the time k +1 for the ith node,is the estimation error at the time instant k +1,to the expectation of the element { · },is composed ofThe transposing of (1).
Because the estimation error covariance matrix has uncertain items, the true value of the estimation error covariance matrix cannot be obtained. Optimization estimation error covariance matrix upper bound sigmai,k+1|k+1The estimated gain matrix K of the ith node at the moment K +1 can be obtainedi,k+1。
The following examples were used to demonstrate the beneficial effects of the present invention:
the first embodiment is as follows:
the method of the invention is adopted for simulation:
system parameters:
B1,k=[0.8 0.7]T,B2,k=[7.5+0.1sin(2k) 0.55]T,B3,k=[0.65 0.85]T,
C1,k=[1.6 1.8],C2,k=[1.6 1.3],C3,k=[1.65 1.8],
Γ=diag{0.8,0.8},
d=3。
other simulation initial values are selected as follows:
x1,0=[0.5 1]T,x2,0=[-1 0.25]T,x3,0=[-0.5 -0.75]T,xi,j=[0.2 0.2]T,(i=1,2,3,j=-1,-2),Q1,k=0.2,Q2,k=0.15,Q3,k=0.1,R1,k=0.2,R2,k=0.1,R3,k=0.1,ε1=0.3,εi=0.1(i=2,3,4),εi=1(i=5,6,7),η1=0.2,η2=5,η3=1,η4=0.5,Σ1,0|0=2I2,Σ2,0|0=2.5I2,Σ3,0|0=3I2,Σi,j|j=10I2(i=1,2,3,j=-1,-2,-3)。
case one (Case I):
δ1=0.1,δ2=0.5,δ3=0.5;
case two (Case II):
δ1=5,δ2=5,δ3=5。
the state estimator effect:
as can be seen from fig. 2a, 2b, 3a, 3b, 4a, and 4b, the state estimator design method of the present invention can effectively estimate the target state for a randomly-coupled time-lapse complex network with measurement loss and inaccurate occurrence probability under an event trigger mechanism.
As can be seen from fig. 5a, 5b, 5c, for each node, the probability is followedIncrease of (1), estimation error covariance matrix upper bound ∑i,k+1|k+1Is decremented.
The present invention is capable of other embodiments and its several details are capable of modifications in various obvious respects, all without departing from the spirit and scope of the present invention.
Claims (6)
1. A complex network optimization estimation method with random coupling is characterized in that: the method comprises the following specific processes:
step one, establishing a complex network dynamic model with randomly generated coupling time lag;
secondly, performing state estimation on the randomly-generated coupling time-lag complex network dynamic model established in the first step under an event trigger protocol;
step three, calculating the upper bound sigma of the one-step prediction error covariance matrix of the state estimation of each node i of the complex networki,k+1|k;
Step four, according to the one-step prediction error covariance matrix upper bound sigma obtained in step threei,k+1|kCalculating an estimated gain matrix K for each node of the complex networki,k+1;
Step five, obtaining the estimated gain matrix K of each node in the step fouri,k+1Substituting the step two to obtain the state estimation of the ith node at the k +1 th moment
Judging whether k +1 reaches the total time length M of the complex network, if k +1 is less than M, executing a sixth step, and if k +1 is equal to M, ending the operation;
step six, according to the estimated gain matrix K of each node of the complex network calculated in the step fouri,k+1Calculating an estimation error covariance matrix upper bound sigma of each node of the complex networki,k+1|k+1;
And k is k +1, and the step two is executed until k +1 is M.
2. The method of claim 1, wherein the method comprises: establishing a random coupling time lag complex network dynamic model in the first step;
the state space form of the complex network dynamic model with the randomly generated coupling time lag is as follows:
yi,k=λi,kCi,kxi,k+νi,k,i=1,2,…,N, (2)
in the formula,
respectively representing the state variables of the ith node of the complex network at the k, k +1 and k-d moments; d is the network skew;is a real number domain of the state of the complex network dynamic model, and n is a dimension; x is the number ofj,kThe state variable of the jth node of the complex network at the kth moment;the measurement output of the ith node at the kth moment;a real number domain is output by a complex network dynamic model, and p is a dimension;is the initial value of the ith node at time k, k-d, -d +1, …, 0; Γ is the connection matrix; w is aijCoupling weights for the ith node and the jth node;is a mean of zero and a variance of Qi,kThe process noise of (a) is generated,a real number domain of the process noise of the complex network dynamic model is formed, and q is a dimension;is a mean of zero and a variance of Ri,kThe measurement noise of (2); n is the number of nodes of the complex network; a. thei,kIs a system matrix, Bi,kAs a noise distribution matrix, Ci,kIs the measurement matrix for the time instant k,of known dimensions and with a time lag k-d a related system matrix; gamma rayi,kAnd λi,kRandom variables obeying Bernoulli distribution respectively depict random coupling phenomenon and measurement loss phenomenon, and satisfy the following conditions:
wherein,
γi,k+1、λi,k+1random variables obeying Bernoulli distribution respectively depict random coupling phenomenon and measurement loss phenomenon,the expected probability of random coupling of the ith node at the moment k + 1;the measurement loss probability of the ith node at the moment k +1 is obtained;andhas a value range of (0,1), Δ γi,k=1Inaccuracy of characterization probability; prob { } is a probability,as desired.
3. The method of claim 2, wherein the method comprises: performing state estimation on the randomly generated coupling time-lag complex network dynamic model established in the step one under the event triggering protocol in the step two; the specific process is as follows:
firstly, aiming at the ith node, selecting the following event trigger formula:
whereinFor the measurement output of the ith node at the last trigger time,for the corresponding last value of the triggering time, δiIf > 0, the adjusting threshold is set, and T is transposition; the next event trigger sequence for the ith nodeIteratively generated by:
whereinIs a positive integer set, inf { } is a lower limit taking function;
after the event triggering mechanism, the true measurement value delivered to the filter is
And aiming at the ith node of the complex network, constructing a state estimator:
in the formula,
is xi,kA one-step prediction at the time instant k,for the state estimate of the ith node at time k +1,for the state estimation of the ith node at the K-d time, Ki,k+1For the estimated gain matrix at time k +1 for the ith node,for the state estimate of the ith node at the kth time,for the state estimate of the jth node at the kth instant,for the measurement output of the ith node at time k +1, Ci,k+1The measurement matrix is the measurement matrix at the moment k + 1;the expected probability of random coupling at time k for the ith node.
4. The method of claim 3, wherein the method comprises: calculating the one-step prediction error covariance matrix upper bound sigma of the state estimation of each node i of the complex network in the step threei,k+1|k(ii) a The specific process is as follows:
calculating each node of the complex network according to the following formulaUpper bound sigma of the covariance matrix of the one-step prediction error of the state estimationi,k+1|k:
In the formula,
τ1、τ2、τ4、τ5is an intermediate variable, wiIs the intermediate variable(s) of the variable,Σi,k+1|kone-step prediction error covariance matrix upper bound, Σ, for the ith node at time ki,k|kThe covariance matrix of the estimation error at the kth time for the ith node is bounded,ΓT,are respectively Ai,k,Γ,Bi,k,Transposing; sigmai,k-d|k-dUpper bound of covariance matrix of estimation error of ith node at k-dj,k|kAnd (4) the estimation error covariance matrix of the jth node at the kth moment is bounded.
5. The method of claim 4, wherein the method comprises: the fourth step is to obtain the upper bound sigma of the one-step prediction error covariance matrix according to the third stepi,k+1|kCalculating an estimated gain matrix K for each node of the complex networki,k+1(ii) a The specific process is as follows:
according to (A)9) Formula an estimated gain matrix K for each node of a complex network is calculated as followsi,k+1:
In the formula,
is the intermediate variable(s) of the variable,ηlis a constant greater than zero, 1,2,3, 4;is etalThe inverse of (a) is,is λi,k+1The variance of (a) is determined,Ki,k+1for the estimated gain matrix of the ith node at time k +1,is Ci,k+1I is an identity matrix, Ri,k+1V isi,k+1The variance of (c) is.
6. The method of claim 5, wherein the method comprises: in the sixth step, the estimated gain matrix K of each node of the complex network calculated in the fourth stepi,k+1Calculating the upper bound sigma of the covariance matrix of the estimation error of each node of the complex networki,k+1|k+1(ii) a The specific process is as follows:
the formula is as follows:
wherein, sigmai,k+1|k+1The covariance matrix of the estimation error at the (k + 1) th node is bounded,andare respectively asAnd Ki,k+1The transposing of (1).
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810814462.5A CN109039725B (en) | 2018-07-23 | 2018-07-23 | Complex network optimization estimation method with random occurrence coupling |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810814462.5A CN109039725B (en) | 2018-07-23 | 2018-07-23 | Complex network optimization estimation method with random occurrence coupling |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109039725A true CN109039725A (en) | 2018-12-18 |
CN109039725B CN109039725B (en) | 2021-07-27 |
Family
ID=64644328
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810814462.5A Active CN109039725B (en) | 2018-07-23 | 2018-07-23 | Complex network optimization estimation method with random occurrence coupling |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109039725B (en) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109088749A (en) * | 2018-07-23 | 2018-12-25 | 哈尔滨理工大学 | The method for estimating state of complex network under a kind of random communication agreement |
CN109728795A (en) * | 2018-12-24 | 2019-05-07 | 哈尔滨理工大学 | Time-varying event under unknown probability situation with loss of data triggers filtering method |
CN113242063A (en) * | 2021-04-29 | 2021-08-10 | 江南大学 | Large-scale MIMO channel model modeling method based on random coupling |
CN115859030A (en) * | 2022-11-29 | 2023-03-28 | 哈尔滨理工大学 | Two-step state estimation method under complex coupling |
CN115883408A (en) * | 2022-11-28 | 2023-03-31 | 哈尔滨理工大学 | Multi-rate complex network state estimation method based on compensation |
CN116088303A (en) * | 2022-12-12 | 2023-05-09 | 哈尔滨理工大学 | Uncertain complex dynamic network state time-varying recursion estimation method |
CN116542541A (en) * | 2023-02-17 | 2023-08-04 | 威海市天罡仪表股份有限公司 | Thermal dynamic estimation method and system for multi-user radiator of next ladder under coding and decoding strategy |
CN117544956A (en) * | 2023-11-06 | 2024-02-09 | 哈尔滨理工大学 | Multi-mobile robot security positioning method based on network communication |
CN117910622A (en) * | 2023-12-28 | 2024-04-19 | 哈尔滨理工大学 | Insect population dynamic forecast estimation method |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106054084A (en) * | 2016-07-06 | 2016-10-26 | 南京航空航天大学 | Power battery SOC estimation method |
CN109088749A (en) * | 2018-07-23 | 2018-12-25 | 哈尔滨理工大学 | The method for estimating state of complex network under a kind of random communication agreement |
-
2018
- 2018-07-23 CN CN201810814462.5A patent/CN109039725B/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106054084A (en) * | 2016-07-06 | 2016-10-26 | 南京航空航天大学 | Power battery SOC estimation method |
CN109088749A (en) * | 2018-07-23 | 2018-12-25 | 哈尔滨理工大学 | The method for estimating state of complex network under a kind of random communication agreement |
Non-Patent Citations (6)
Title |
---|
HONGXU ZHANG等: "Recursive state estimation for time-varying complex networks subject to missing measurements and stochastic inner coupling under random access protocol", 《NEUROCOMPUTING》 * |
JUN HU等: "One state estimation for nonlinear dynamical networks with random sensor delays and coupling strength under event-based communication mechanism", 《INFORMATION SCIENCES》 * |
WENLING LI等: "Recursive state estimation for complex networks with random coupling strength", 《NEUROCOMPUTING》 * |
WENLING LI等: "State estimation for stochastic complex networks with switching topolpgy", 《IEEE TRANSACTIONS ON AUTOMATIC CONTROL》 * |
王立成: "基于信息不完整模型的复杂网络系统控制与状态估计", 《中国优秀硕士学位论文全文数据库(电子期刊)》 * |
陈东彦等: "具有相关噪声和不确定观测系统的全局最优Kalman滤波", 《哈尔滨理工大学学报》 * |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109088749A (en) * | 2018-07-23 | 2018-12-25 | 哈尔滨理工大学 | The method for estimating state of complex network under a kind of random communication agreement |
CN109088749B (en) * | 2018-07-23 | 2021-06-29 | 哈尔滨理工大学 | State estimation method of complex network under random communication protocol |
CN109728795A (en) * | 2018-12-24 | 2019-05-07 | 哈尔滨理工大学 | Time-varying event under unknown probability situation with loss of data triggers filtering method |
CN109728795B (en) * | 2018-12-24 | 2022-11-18 | 哈尔滨理工大学 | Time-varying event trigger filtering method with data loss under unknown probability condition |
CN113242063A (en) * | 2021-04-29 | 2021-08-10 | 江南大学 | Large-scale MIMO channel model modeling method based on random coupling |
CN115883408A (en) * | 2022-11-28 | 2023-03-31 | 哈尔滨理工大学 | Multi-rate complex network state estimation method based on compensation |
CN115883408B (en) * | 2022-11-28 | 2023-05-30 | 哈尔滨理工大学 | Multi-rate complex network state estimation method based on compensation |
CN115859030A (en) * | 2022-11-29 | 2023-03-28 | 哈尔滨理工大学 | Two-step state estimation method under complex coupling |
CN116088303A (en) * | 2022-12-12 | 2023-05-09 | 哈尔滨理工大学 | Uncertain complex dynamic network state time-varying recursion estimation method |
CN116542541A (en) * | 2023-02-17 | 2023-08-04 | 威海市天罡仪表股份有限公司 | Thermal dynamic estimation method and system for multi-user radiator of next ladder under coding and decoding strategy |
CN116542541B (en) * | 2023-02-17 | 2024-03-26 | 威海市天罡仪表股份有限公司 | Thermal dynamic estimation method and system for multi-user radiator of next ladder under coding and decoding strategy |
CN117544956A (en) * | 2023-11-06 | 2024-02-09 | 哈尔滨理工大学 | Multi-mobile robot security positioning method based on network communication |
CN117544956B (en) * | 2023-11-06 | 2024-04-30 | 哈尔滨理工大学 | Multi-mobile robot safety positioning method based on network communication |
CN117910622A (en) * | 2023-12-28 | 2024-04-19 | 哈尔滨理工大学 | Insect population dynamic forecast estimation method |
Also Published As
Publication number | Publication date |
---|---|
CN109039725B (en) | 2021-07-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109039725B (en) | Complex network optimization estimation method with random occurrence coupling | |
CN109088749B (en) | State estimation method of complex network under random communication protocol | |
Liu et al. | Extended dissipative sliding mode control for nonlinear networked control systems via event-triggered mechanism with random uncertain measurement | |
CN109696669B (en) | Event-triggered multi-sensor fusion estimation method in correlated noise environment | |
CN109728796B (en) | Filtering method based on event trigger mechanism | |
CN107994885B (en) | Distributed fusion filtering method for simultaneously estimating unknown input and state | |
WO2021036414A1 (en) | Co-channel interference prediction method for satellite-to-ground downlink under low earth orbit satellite constellation | |
CN108153259B (en) | Multi-controller optimal state estimation control strategy design method based on Kalman filtering | |
CN108847828B (en) | Nonlinear event trigger filtering method with random modeling error | |
CN116088303B (en) | Uncertain complex dynamic network state time-varying recursion estimation method | |
CN109728795B (en) | Time-varying event trigger filtering method with data loss under unknown probability condition | |
CN107292410A (en) | tunnel deformation prediction method and device | |
CN110377955A (en) | Distributed set-membership filtering device design method based on event trigger mechanism | |
Liu et al. | Robust measurement fusion steady-state estimator design for multisensor networked systems with random two-step transmission delays and missing measurements | |
Yao et al. | Recursive distributed fusion estimation for multi-sensor systems with missing measurements, multiple random transmission delays and packet losses | |
Li et al. | Outlier-resistant interval observer design for multirate time-delayed systems under the adaptive event-triggered protocols | |
CN112269318B (en) | Finite time remote safety state estimation method for time delay uncertain system | |
CN115865702B (en) | Distributed fusion estimation method with data attenuation under network scheduling strategy | |
CN111462479A (en) | Traffic flow prediction method based on Fourier-recurrent neural network | |
Caballero-Águila et al. | Covariance‐Based Estimation from Multisensor Delayed Measurements with Random Parameter Matrices and Correlated Noises | |
Caballero-Águila et al. | Recursive least-squares quadratic smoothing from measurements with packet dropouts | |
CN115883408A (en) | Multi-rate complex network state estimation method based on compensation | |
CN115859030B (en) | Two-step state estimation method under complex coupling | |
CN118157630B (en) | Distributed filtering method of positioning system based on sensor network | |
CN117544956B (en) | Multi-mobile robot safety positioning method based on network communication |
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 |