
Citation: Luca Bertelli, Frédéric Gibou. Fast two dimensional to three dimensional registration of fluoroscopy and CT-scans using Octrees on segmentation maps[J]. Mathematical Biosciences and Engineering, 2012, 9(3): 527-537. doi: 10.3934/mbe.2012.9.527
[1] | Chun-Chao Yeh, Ke-Jia Jhang, Chin-Chun Chang . An intelligent indoor positioning system based on pedestrian directional signage object detection: a case study of Taipei Main Station. Mathematical Biosciences and Engineering, 2020, 17(1): 266-285. doi: 10.3934/mbe.2020015 |
[2] | Xiaohao Chen, Maosheng Fu, Zhengyu Liu, Chaochuan Jia, Yu Liu . Harris hawks optimization algorithm and BP neural network for ultra-wideband indoor positioning. Mathematical Biosciences and Engineering, 2022, 19(9): 9098-9124. doi: 10.3934/mbe.2022423 |
[3] | Mingju Chen, Hongyang Li, Hongming Peng, Xingzhong Xiong, Ning Long . HPCDNet: Hybrid position coding and dual-frquency domain transform network for low-light image enhancement. Mathematical Biosciences and Engineering, 2024, 21(2): 1917-1937. doi: 10.3934/mbe.2024085 |
[4] | Diguo Zhai, Xinqi Bao, Xi Long, Taotao Ru, Guofu Zhou . Precise detection and localization of R-peaks from ECG signals. Mathematical Biosciences and Engineering, 2023, 20(11): 19191-19208. doi: 10.3934/mbe.2023848 |
[5] | Guohui Zhang, Gulbahar Tohti, Ping Chen, Mamtimin Geni, Yixuan Fan . SFEMM: A cotton binocular matching method based on YOLOv7x. Mathematical Biosciences and Engineering, 2024, 21(3): 3618-3630. doi: 10.3934/mbe.2024159 |
[6] | Kezhou Chen, Xu Lu, Rongjun Chen, Jun Liu . Wireless wearable biosensor smart physiological monitoring system for risk avoidance and rescue. Mathematical Biosciences and Engineering, 2022, 19(2): 1496-1514. doi: 10.3934/mbe.2022069 |
[7] | Bin Zhang, Linkun Sun, Yingjie Song, Weiping Shao, Yan Guo, Fang Yuan . DeepFireNet: A real-time video fire detection method based on multi-feature fusion. Mathematical Biosciences and Engineering, 2020, 17(6): 7804-7818. doi: 10.3934/mbe.2020397 |
[8] | Yongmei Ren, Xiaohu Wang, Jie Yang . Maritime ship recognition based on convolutional neural network and linear weighted decision fusion for multimodal images. Mathematical Biosciences and Engineering, 2023, 20(10): 18545-18565. doi: 10.3934/mbe.2023823 |
[9] | Wen Li, Xuekun Yang, Guowu Yuan, Dan Xu . ABCNet: A comprehensive highway visibility prediction model based on attention, Bi-LSTM and CNN. Mathematical Biosciences and Engineering, 2024, 21(3): 4397-4420. doi: 10.3934/mbe.2024194 |
[10] | Yujun Zhang, Yufei Wang, Shuijia Li, Fengjuan Yao, Liuwei Tao, Yuxin Yan, Juan Zhao, Zhengming Gao . An enhanced adaptive comprehensive learning hybrid algorithm of Rao-1 and JAYA algorithm for parameter extraction of photovoltaic models. Mathematical Biosciences and Engineering, 2022, 19(6): 5610-5637. doi: 10.3934/mbe.2022263 |
Along with the rapid development of biology and biotechnology, the regulation and interaction of mRNAs and the produced proteins in genetic regulatory networks (GRNs) have attracted great research attention, and various related applications have been available to address the burgeoning areas [1,2,3]. Furthermore, much effort has been devoted to the biochemical dynamics of GRNs, which can better describe the expression mechanisms of genes from a dynamical system perspective [4,5]. The dynamical behaviors of GRNs can considerably provide deeper insights into the processing conditions. Generally speaking, differential equation and Boolean network models are employed and investigated to characterize complex properties of mRNA and protein concentrations in GRNs [6,7,8]. In particular, biological evidences have proved the existences of jumping mechanisms of GRNs. This means that the dynamical parameters of GRNs are varying according to different modes rather than kept constant owing to random abrupt changes during gene expressions. This means that the parameters of GRNs should be modified according to certain modes. With this context, mathematical models of Markovian jumping systems with Markov chains have been widely employed to mimic the intrinsic variate of GRNs [9,10,11]. As a result, many remarkable results on Markovian jumping GRNs have been reported in the literature, which includes stochastic state estimation [12], stability [13] and synchronization issues [14].
In addition, the accessibility of true jumping mode information of GRNs has also drawn much interest, especially for some actual applications. It is noteworthy that the exact mode information of GRNs in practical applications is always difficult or expensive to acquire, which is mostly because of detection delays or unpredictable disturbances during the transcription and translation processes [15,16]. Especially, the mode information characterized by Markov chains for state estimation or control issues is also challenging to main synchronous design. On the other hand, nonsynchronous analysis and synthesis methods are more meaningful and practical for Markovian jumping GRNs with mode-dependent designs. Lately, hidden Markovian jumping systems have been initially investigated by external mode observations and encouraging results of the so-called nonsynchronous designs have been reported accordingly. More precisely, on the basis of conditional probability between true modes and observed modes, the observed models for mode information are effectively utilized without synchronous mode operations [17,18,19].
On another active research line, distributed strategies for control systems have been extensively developed due to their increasing requirements for efficiency and robustness in engineering areas. Compared with traditional control techniques, the distributed architecture can process and fuse the information via multiple sensors with each domain data. With respect to its advances, many remarkable analysis and synthesis of distributed control systems can be found in recent works [20,21,22]. To name a few, in [23], a novel distributed-observer-based distributed control design is developed for affine nonlinear systems and is effectively applied to interconnected cruise control of intelligent vehicles. In [24], the distributed observer and controller design issue is discussed for spatially interconnected systems. In [25], the distributed state observers and disturbance observers are proposed for cooperative output regulation of LTI system to ensure the reference tracking properly. In [26], a distributed Nash equilibrium seeking scheme is designed with nonlinear dynamical system. In [27], the leader-following output consensus problem is concerned with distributed nonlinear observers. These notable works have proven that distributed frameworks can effective improve the system performance for complex dynamical systems. However, the problems of network communication burden must be faced among the nodes. It should be emphasized that traditional time-triggered periodic communication schedules have certain drawbacks in computation and network resources. Fortunately, the emerging event-triggered schemes are developed aiming at necessary transmission only when speechified triggering conditions are satisfied [28,29,30]. Hence, the event-triggered mechanism could lead to less network consumptions and could be more reliable in line with practical network implementation. Furthermore, in view of distributed control systems, distinguishing advantages can also been made by adopting event-triggered communications [31,32,33]. Nevertheless, there are challenging problems for distributed event-triggered control in corresponding algorithm complexity. For the consideration of GRNs with HMJPs, it is reasonable and significant to apply the distributed event-triggered design. Unfortunately, until now, distributed event-triggered state estimation problem of GRNs has not been entirely studied, not to mention the case of HMJPs with NSE. This motivates the authors for current study to shorten such a gap.
To deal with the aforementioned issue, the key idea in this paper is to investigate the distributed event-triggered state estimation problem for GRNs with HMJPs based on the NSE scheme. Compared with most existing literature, the primary novelties of our work are highlighted as follows:
1) The attempt to deal with the distributed NSE design of GRNs with HMJPs is made for the first time, which can better model GRNs dynamics and mode jumping phenomena during the state estimation procedures.
2) The distributed nonsynchronous state estimation framework is also introduced for GRNs with nonsynchronous event-triggered communications, such that more estimation efficiency for concentrations of mRNAs and proteins can be achieved with effective mode information utilization.
3) The sufficient mode-dependent analysis criteria are established by novel Lyapunov-Krasovski functions, and then the desired nonsynchronous mode-dependent estimator gains can be derived with convex optimization.
The remainder of this paper is stated by the following lines. Section 2 introduces some background information regarding GRNs with HMJPs and formulates the state estimation problem. Subsequent session of Section 3 presents the main theoretical developments with details. Section 4 approves the correctness of our developed results through a numerical example. Section 5 concludes the paper and addresses future challenges.
All referred notations are standardly utilized. $ \mathscr{R}^{n} $ stands for $ n $-Euclidean space. Matrix $ \mathbb{P} > 0 $ means that $ \mathbb{P} $ is symmetric positive definite, and vise versa. $ (\mathbb{O}, \mathbb{F}, \mathbb{P}) $ corresponds a probability space, where $ \mathbb{O} $ denotes sample space, $ \mathbb{F} $ represents $ \delta $-algebra of subsets within $ \mathbb{O} $ and $ \mathbb{P} $ stands for probability measurement. $ \mathbb{E} $ means the mathematical expectation and $ * $ implies a induced term in symmetry matrices, respectively.
Consider the equilibrium points of GRNs model shifted to the origin and the following GRNs model with HMJPs can be described by:
$ (ΣM):{˙xm(t)=−A(r(t))xm(t)+B(r(t))g(xp(t−dp(t)))+Em(r(t))ω1(t),˙xp(t)=−C(r(t))xp(t)+D(r(t))xm(t−dm(t))+Ep(r(t))ω2(t),zm(t)=Um(r(t))xm(t),zp(t)=Up(r(t))xp(t),xm(t)=φm(t),xp(t)=φp(t),∀t∈[−max(ˉdp,ˉdm),0), $ | (2.1) |
where $ x_{m}\left(t\right) = [x_{m1}\left(t\right), x_{m2}\left(t\right), \ldots, x_{mn}\left(t\right) ]^{T}\in \mathscr{R}^{n} $, $ x_{p}\left(t\right) = [x_{p1}\left(t\right), x_{p2}\left(t\right), \ldots, x_{pn}\left(t\right) ]^{T}\in \mathscr{R}^{n} $ with $ x_{mi}\left(t\right) $ and $ x_{pi}\left(t\right), $ $ i = 1, 2, \ldots, n $ representing the concentrations of mRNAs and proteins, respectively; $ z_{m}\left(t\right) = [z_{m1}\left(t\right), z_{m2}\left(t\right), \ldots, z_{ml}\left(t\right) ]^{T}\in \mathscr{R} ^{l} $, $ z_{p}\left(t\right) = [z_{p1}\left(t\right), z_{p2}\left(t\right), \ldots, z_{pl}\left(t\right) ]^{T}\in \mathscr{R}^{l} $ with $ z_{mk}\left(t\right) $ and $ z_{pk}\left(t\right), $ $ k = 1, 2, \ldots, l $ being the concentrations of certain mRNAs and proteins to be estimated through expression level; $ g\left(x_{p}\left(t-d_{p}(t)\right) \right) = [g_{1}\left(x_{p1}\left(t-d_{p}(t)\right) \right), g_{2}\left(x_{p2}\left(t-d_{p}(t)\right) \right), \ldots, g_{n}\left(x_{pn} \left(t-d_{p}(t)\right) \right) ]^{T}\in \mathscr{R}^{n} $ stands for the nonlinear feedback regulation function of the protein on transcription, which satisfies the following assumption
$ 0≤gi(xi)xi≤ςi,xi≠0;gi(0)=0,i=1,2,…,n,ς=diag{ς1,ς2,…,ςn}>0. $ |
i.e.,
$ gT(x)(g(x)−ςx)≤0; $ |
$ 0\leq d_{p}(t)\leq \bar{d}_{p} $ and $ 0\leq d_{m}(t)\leq \bar{d}_{m} $ are the translation delay and the feedback regulation delay, respectively$; $ $ \omega _{1}(t) $ and $ \omega _{2}(t) $ refers to the external disturbances that belongs to $ \mathcal{L}_{2}[0, \infty) $.
The degradation or dilution rates of mRNAs and proteins can be described by known matrices $ A\left(r\left(t\right) \right) = diag\{a_{1}\left(r\left(t\right) \right), a_{2}\left(r\left(t\right) \right), \ldots, a_{n}\left(r\left(t\right) \right) \} $, $ C\left(r\left(t\right) \right) = diag\{c_{1}\left(r\left(t\right) \right), c_{2}\left(r\left(t\right) \right), \ldots, c_{n}\left(r\left(t\right) \right) \} $, $ D\left(r\left(t\right) \right) = diag\{d_{1}\left(r\left(t\right) \right), d_{2}\left(r\left(t\right) \right), \ldots, d_{n}\left(r\left(t\right) \right) \} $, and the coupling matrix $ B\left(r\left(t\right) \right) = (b_{ij}\left(r\left(t\right) \right))\in \mathscr{R}^{n} $ can be defined by the following
$ bij(r(t))={ˉbij(r(t))0−ˉbij(r(t))iftranscriptionfactorjisanactivatorofgenei,ifthereisnolinkfromnodejtoi,iftranscriptionfactorjisanrepressorofgenei. $ |
Furthermore, $ U_{m}\left(r\left(t\right) \right) $, $ U_{p}\left(r\left(t\right) \right), $ $ E_{m}\left(r\left(t\right) \right) $ and $ E_{p}\left(r\left(t\right) \right) $ are also known matrices for a certain mode $ r\left(t\right). $
Moreover, the continuous-time discrete-state Markov process $ r\left(t\right) \in \mathcal{S} = \{1, \ldots, \mathcal{N}\} $ on $ (\mathbb{O}, \mathbb{F}, \mathbb{P}) $ is given to model the HMJPs, which can be defined by $ \Pi = (\pi _{ij})_{\mathcal{N}\times \mathcal{N}} $, with
$ Pr(r(t+Δ)=j|r(t)=i)={πijΔ+o(Δ), i≠j,1+πiiΔ+o(Δ), i=j,, $ | (2.2) |
and $ \pi _{ii} = -\underset{j = 1, \;\;i\neq j}{\overset{\mathcal{N}}{\sum \limits}} \pi _{ij} $.
In this paper, a direct graph $ \mathbb{G = }\left\{ \mathbb{V}, \mathbb{E}, \mathbb{A} \right\} $ is given to describe the communication topology of $ N $ sensor nodes, where $ \mathbb{V = }\left\{ v_{1}, v_{2}, \cdots, v_{N}\right\} $ and $ \mathbb{E}\subseteq \mathbb{V}\times \mathbb{V} $ stand for the set of nodes and edges, respectively, $ \mathbb{A = }\left[ a_{ij}\right] \in \mathbb{R}^{N\times N} $ represent the adjacency matrix with $ a_{ii} = 0 $ for any $ i $. $ \mathbb{A} $ is associated with the edges of $ \mathbb{G} $ are positive, i.e., $ a_{ij} > 0\Longleftrightarrow \varepsilon _{ij}\in \mathbb{E} $. In addition, the corresponding Laplacian matrix of $ \mathbb{G} $ is defined as $ \mathbb{L} $.
Moreover, the distributed measurement is achieved by $ N $ mode-dependent sensor nodes as follows:
$ (ΣS):{ylm(t)=Hlm(r(t))xm(t)+Glm(r(t))ω3(t),ylp(t)=Hlp(r(t))xp(t)+Glp(r(t))ω4(t), $ | (2.3) |
where $ y_{m}^{l}(t) $ and $ y_{p}^{l}(t) $, $ l = 1, 2, \ldots, N, $ represent the measurements of mRNAs and proteins on sensor $ l $, respectively; $ \omega _{3}(t) $ and $ \omega _{4}(t) $ are measurement disturbances; $ H_{m}^{l}(r(t)), $ $ H_{p}^{l}(r(t)), $ $ G_{m}^{l}(r(t)) $, $ G_{p}^{l}(r(t)) $ are known constant matrices.
Then, define the sensor estimation errors $ \tilde{y}_{m}^{l}(t) $ and $ \tilde{y}_{p}^{l}(t) $ as follows
$ {˜ylm(t)=ylm(t)−Hlm(r(t))ˆxlm(t),˜ylp(t)=ylp(t)−Hlp(r(t))ˆxlp(t), $ | (2.4) |
where $ \hat{x}^{l}_{m}(t) $ and $ \hat{x}^{l}_{p}(t) $ denote the filter state for sensor $ l $.
By considering the event-triggered strategy for data transmission with ZOHs, it is assumed that the sampling period $ h_{k} $ with sampling instants $ t_{k} $ is assumed to be $ h_{k} = t_{k+1}-t_{k}\leq \bar{h}, $ $ \bar{h} > 0. $ Furthermore, the transmission instant $ t_{\delta }^{l} $ for sensor $ l $ is determined by the event generation functions as follows:
$ tlδ+1=mintk>tlδ{tk‖˜yl(tk)−˜yl(tlδ)‖≥κ‖˜yl(tlδ)‖}, $ | (2.5) |
where $ \tilde{y}^{l}(t_{k}) = [(\tilde{y}_{m}^{l}(t))^{T}, (\tilde{y} _{p}^{l}(t))^{T}]^{T} $ and $ 0 < \kappa < 1 $ is a scalar parameter.
Consequently, the following distributed mode-dependent state estimation filter $ (\Sigma _{F}) $ can be designed for system $ (\Sigma _{M}) $:
$ (ΣF):{˙ˆxlm(t)=−A(r(t))ˆxm(t)+B(r(t))g(ˆxp(t−dp(t)))+Flm(σ(t))˜ylm(tlδ)+Klm(σ(t))N∑s=1als[˜ysm(tsδ∗)−˜ylm(tlδ)],˙ˆxlp(t)=−C(r(t))ˆxp(t)+D(r(t))ˆxm(t−dm(t))+Flp(σ(t))˜ylp(tlδ)+Klp(σ(t))N∑s=1als[˜ysp(tsδ∗)−˜ylp(tlδ)],ˆzlm(t)=Um(r(t))ˆxlm(t),ˆzlp(t)=Up(r(t))ˆxlp(t), $ | (2.6) |
where $ t\in \left[ t_{\delta }^{l}, t_{\delta +1}^{l}\right), $ $ \delta ^{\ast } = \arg \min_{\mu }\left\{ t-t_{\mu }^{s}|t\geq t_{\mu }^{s}\right\}, $ $ \hat{z}_{m}^{l}\left(t\right) $ and $ \hat{z}_{p}^{l}\left(t\right) \ $ stand for estimation of $ z_{m}\left(t\right) $ and $ z_{p}\left(t\right), $ respectively; $ F_{m}^{l}(\sigma (t)), F_{p}^{l}(\sigma (t)) $ and $ K_{m}^{l}(\sigma (t)), K_{p}^{l}(\sigma (t)) $ represent the nonsynchronous mode-dependent filter gains to be determined later. More precisely, $ \sigma (t)\in \mathcal{F} = \{1, \ldots, F\} $ represents another Markov process corresponds to $ r\left(t\right) $ with
$ Pr{σ(t)=ρ|r(t)=i}=λiρ, $ |
and $ \sum_{\rho = 1}^{F}\lambda _{i\rho } = 1. $
Remark 1. As event-triggered control strategy for networked control systems becomes more important, it is reasonable to investigate the analysis and synthesis issues for GRNs over sensor networks with event-triggered schemes. In comparison with common time-triggered strategies with extensive data transmissions, the event-triggered framework can considerably reduce the communication burden by generated events. Moreover, our developed event-triggered mechanism is related with both measurements of mRNAs and proteins, which is more efficient for potential applications.
Then, it can be deduced by dividing $ \left[ t_{\delta }^{l}, t_{\delta +1}^{l}\right) $ into $ \cup _{\iota_{k} = t_{\delta }^{l}}^{t_{\delta +1}^{l}-h_{k+1}}\left[ \iota_{k}, \iota_{k+1}\right) $ that
$ {˙ˆxlm(t)=−A(r(t))ˆxlm(t)+B(r(t))g(ˆxlp(t−dp(t)))+Flm(σ(t))˜ylm(ιk)−Flm(σ(t))εlm(ιk)+Klm(σ(t))N∑s=1als[˜ysm(ιk)−˜ylm(ιk)]−Klm(σ(t))N∑s=1als[εsm(ιk)−εlm(ιk)],˙ˆxlp(t)=−C(r(t))ˆxlp(t)+D(r(t))ˆxlm(t−dm(t))+Flp(σ(t))˜ylp(ιk)−Flp(σ(t))εlp(ιk)+Klp(σ(t))N∑s=1als[˜ysp(ιk)−˜ylp(ιk)]−Klp(σ(t))N∑s=1als[εsp(ιk)−εlp(ιk)],ˆzlm(t)=Um(r(t))ˆxlm(t),ˆzlp(t)=Up(r(t))ˆxlp(t), $ | (2.7) |
where $ \varepsilon _{m}^{l}(\iota _{k}) = \tilde{y}_{m}^{l}\left(\iota _{k}\right) - \tilde{y}_{m}^{l}(t_{\delta }^{l}) $ and $ \varepsilon _{p}^{l}(\iota _{k}) = \tilde{y}_{p}^{l}\left(\iota _{k}\right) -\tilde{y}_{p}^{l}(t_{\delta }^{l}), $ $ t_{\delta }^{l}\leq \iota _{k} < t_{\delta +1}^{l}, $ represent the event-triggered measurement errors, respectively.
As a result, by denoting $ e_{m}^{l}(t) = x(t)-\hat{x}_{m}^{l}(t), $ $ e_{p}^{l}(t) = x(t)-\hat{x}_{p}^{l}(t), $ $ \tilde{z}_{m}^{l}(t) = z(t)-\hat{z} _{m}^{l}(t) $ and $ \tilde{z}_{p}^{l}(t) = z(t)-\hat{z}_{p}^{l}(t) $, one can deduce that
$ {˙elm(t)=−A(r(t))elm(t)+B(r(t))g(xp(t−dp(t)))−B(r(t))g(ˆxlp(t−dp(t)))−Flm(σ(t))Hlm(r(t))elm(ιk)+Flm(σ(t))εlm(ιk)−Klm(σ(t))N∑s=1als[Hsm(r(t))esm(ιk)−Hlm(r(t))elm(ιk)]+Km(σ(t))N∑s=1als[εsm(ιk)−εlm(ιk)]+Em(r(t))ω1(t)−Flm(σ(t))Glm(r(t))ω3(ιk)−Klm(σ(t))N∑s=1als[Gsm(r(t))ω3(ιk)−Glm(r(t))ω3(ιk)],˙elp(t)=−C(r(t))elp(t)+D(r(t))elm(t−dm(t))−Flp(σ(t))Hlp(r(t))elp(ιk)+Flp(σ(t))εlp(ιk)−Klp(σ(t))N∑s=1als[Hsp(r(t))esp(ιk)−Hlp(r(t))elp(ιk)]+Klp(σ(t))N∑s=1als[εsp(ιk)−εlp(ιk)]+Ep(r(t))ω2(t)−Flp(σ(t))Glp(r(t))ω4(ιk)−Klp(σ(t))N∑s=1als[Gsp(r(t))ω4(ιk)−Glp(r(t))ω4(ιk)],˜zlm(t)=Um(r(t))elm(t),˜zlp(t)=Up(r(t))elp(t), $ | (2.8) |
For simplicity, the resulting closed-loop filter system can be rewritten as follows:
$ {˙em(t)=−ˉA(i)em(t)+ˉB(i)ˉg(ep(t−dp(t)))−(ˉFm(ρ)−ˉKm(ρ)ˉL)ˉHm(i)em(ιk)+(ˉFm(ρ)−ˉKm(ρ)ˉL)εm(ιk)+ˉEm(i)ω1(t)−(ˉFm(ρ)−ˉKm(ρ)ˉL)ˉGm(i)ω3(ιk),˙ep(t)=−ˉC(i)ep(t)+ˉD(i)em(t−dm(t))−(ˉFp(ρ)−ˉKp(ρ)ˉL)ˉHp(i)ep(ιk)+(ˉFp(ρ)−ˉKp(ρ)ˉL)εp(ιk)+ˉEp(i)ω2(t)−(ˉFp(ρ)−ˉKp(ρ)ˉL)ˉGp(i)ω4(ιk),˜zm(t)=ˉUm(i)em(t),˜zp(t)=ˉUp(i)ep(t), t∈[ιk,ιk+1), $ | (2.9) |
where
$ em(t)=[eTm1(t),eTm2(t),…,eTmN(t)]T,ep(t)=[eTp1(t),eTp2(t),…,eTpN(t)]T,˜zm(t)=[˜zTm1(t),˜zTm2(t),…,˜zTmN(t)]T,˜zp(t)=[˜zTp1(t),˜zTp2(t),…,˜zTpN(t)]T,εm(t)=[εTm1(t),εTm2(t),…,εTmN(t)]T,εp(t)=[εTp1(t),εTp2(t),…,εTpN(t)]T,ˉA(i)=IN⊗A(i),ˉB(i)=IN⊗B(i),g(ep(t))=[(g(ep(t))−g(ˆep1(t)))T,(g(ep(t))−g(ˆep2(t)))T,…,(g(ep(t))−g(ˆep2(t)))T]T,ˉFm(ρ)=diag{F1m(ρ),F2m(ρ),…,FNm(ρ)},ˉFp(ρ)=diag{F1p(ρ),F2p(ρ),…,FNp(ρ)},ˉKm(ρ)=diag{K1m(ρ),K2m(ρ),…,KNm(ρ)},ˉKp(ρ)=diag{K1p(ρ),K2p(ρ),…,KNp(ρ)},ˉHm(i)=diag{H1m(i),H2m(i),…,HNm(i)},ˉHp(i)=diag{H1p(i),H2p(i),…,HNp(i)},ˉEm(i)=1N⊗Em(i),ˉEp(i)=1N⊗Ep(i),ˉUm(i)=1N⊗Um(i),ˉUp(i)=1N⊗Up(i),ˉGm(i)=[(G1m(i))T,(G2m(i))T,…,(GNm(i))T]T,ˉGp(i)=[(G1p(i))T,(G2p(i))T,…,(GNp(i))T]T,ˉL=L⊗In. $ |
Remark 2. It is noticed that the adoption of conditional probability has been proved to be effective to cope with the Markov or semi-Markov jump systems with limited mode information. For the GRNs with hidden Markovian jumping parameters, lack of true operation modes of Markov chains would lead to control performance degradation or even unstable with mode-dependent designs. In this paper, the nonsynchronous mode information corresponds to the distributed sensor mode information, which implies that mode information mismatch between true GRNs model and our developed distributed state estimators is considered during the design. As a consequence, we have developed a novel distributed nonsynchronous filter framework, where the observed Markovian jumping modes are utilized for more practical distributed event-triggered state estimators.
Consequently, by applying the input-delay approach, it can be obtained that
$ {˙em(t)=−ˉA(i)em(t)+ˉB(i)ˉg(ep(t−dp(t)))−(ˉFm(ρ)−ˉKm(ρ)ˉL)ˉHm(i)em(t−τ(t))+(ˉFm(ρ)−ˉKm(ρ)ˉL)εm(t−τ(t))+ˉEm(i)ω1(t)−(ˉFm(ρ)−ˉKm(ρ)ˉL)ˉGm(i)ω3(t−τ(t)),˙ep(t)=−ˉC(i)ep(t)+ˉD(i)em(t−dm(t))−(ˉFp(ρ)−ˉKp(ρ)ˉL)ˉHp(i)ep(t−τ(t))+(ˉFp(ρ)−ˉKp(ρ)ˉL)εp(t−τ(t))+ˉEp(i)ω2(t)−(ˉFp(ρ)−ˉKp(ρ)ˉL)ˉGp(i)ω4(t−τ(t)),˜zm(t)=ˉUm(i)em(t),˜zp(t)=ˉUp(i)ep(t), t∈[ιk,ιk+1), $ | (2.10) |
where $ \tau (t): = t-\iota _{k} $ with $ 0\leq \tau (t)\leq \bar{\tau}. $
The aim of our paper is to design the desired nonsynchronous state estimator gains, such that the modified $ H_{\infty } $ performance $ \gamma $ can be satisfied in the mean-square sense for augmented state estimation error system with zero initial conditions, if it can hold that
$ E{∫∞0(˜zTm(t)˜zm(t)+˜zTp(t)˜zp(t))dt}≤γ∫∞0(ωT1(t)ω1(t)+ωT2(t)ω2(t)+ωT3(t−τ(t))ω3(t−τ(t))+ωT4(t−τ(t))ω4(t−τ(t))dt. $ | (2.11) |
In the sequel, the following lemmas are employed for later theoretical derivations.
Lemma 1. [34] For any matrix $ \mathbb{M} > 0 $, scalars $ \bar{\xi} > 0 $, $ \xi (t) $ satisfying $ 0\leq \xi (t)\leq \bar{\xi} $, vector function $ \dot{x} (t):[-\bar{\xi}, 0]\rightarrow \mathbb{R}^{n} $ such that the concerned integrations are well defined in (2.12), then
$ −τ∫tt−ˉξ˙xT(s)M˙x(s)ds≤ζT(t)Uζ(t), $ | (2.12) |
where
$ ζ(t)=[xT(t),xT(t−ξ(t)),xT(t−ˉξ)]T,U=[−MM0∗−2MM∗∗−M]. $ |
Lemma 2. [35] If there exists a parameter $ \varepsilon $ and real matrices $ \mathbb{X} $, $ \mathbb{Y}, $ $ \mathbb{Z} $, $ \mathbb{W} $, and it holds that
$ [XZ+εWT∗−εYT−εY]<0, $ | (2.13) |
then one has
$ X+ZTε−TWT+Wε−1Z<0. $ | (2.14) |
In this section, sufficient modified $ H_{\infty } $ performance conditions are first established for the augmented state estimation error system and the distributed state estimator gains are further designed in terms of convex optimization.
Theorem 1. For given $ \bar{d}_{p} $, $ \bar{d}_{m} $, $ \bar{\tau} $ and mode-dependent filter gains $ F_{m}^{l}(\rho), F_{p}^{l}(\rho) $ and $ K_{m}^{l}(\rho), K_{p}^{l}(\rho) $, the modified $ H_{\infty} $ filtering can be satisfied over distributed sensor networks according to Definition 1, if there exist mode-dependent matrices $ P_{m}(i) > 0 $, $ P_{p}(i) > 0 $, and matrices $ Q_{m} > 0 $, $ Q_{p} > 0 $, $ Q_{\tau} > 0 $, $ R_{m} > 0 $, $ R_{p} > 0 $, $ R_{\tau} > 0 $, such that $ \Xi _{i} < 0 $, for all $ i\in \mathcal{N} $, $ k = 1, 2, \ldots, N $ and $ \rho\in \mathcal{F} $, $ k = 1, 2, \ldots, F $, where
$ Ξ(i)=[Ξ1(i)Ξ2(i)∗Ξ3(i)],Ξ1(i)=[Ξ11(i)Ξ12(i)∗Ξ13(i)],Ξ11(i)=[Ξ111(i)Ξ112(i)00Rm00∗Ξ113(i)Rτ0000∗∗−Qτ−Rτ0000∗∗∗−Qm−RmRTm00∗∗∗∗−2Rm00∗∗∗∗∗−Qp−RpRp∗∗∗∗∗∗−2Rp],Ξ111(i)=−2Pm(i)ˉA(i)+Qm+Qτ−Rm−Rτ+ˉUTm(i)ˉUm(i)+∑j∈NπijPm(j),Ξ112(i)=−F∑ρ=1λiρPm(i)ˉFm(ρ)ˉHm(i)+F∑ρ=1λiρPm(i)ˉKm(ρ)ˉLˉHm(i)+Rτ,Ξ113(i)=−2Rτ+κ2ˉHTm(i)ˉHm(i),Ξ12(i)=[000Pm(i)ˉB(i)Ξ111(i)0000−κ2ˉHTm(i)0000000000ˉDT(i)Pp(i)000000000Rp00ς0],Ξ121(i)=F∑ρ=1λiρPm(i)ˉFm(ρ)−F∑ρ=1λiρPm(i)ˉKm(ρ)ˉL,Ξ13(i)=[Ξ131(i)Ξ132(i)000∗Ξ133(i)Rτ00∗∗−Rτ00∗∗∗−2I0∗∗∗∗κ2I−I],Ξ131(i)=−2Pp(i)ˉC(i)+Qp+Qτ−Rp−Rτ+ˉUTp(i)ˉUp(i)+∑j∈NπijPp(i),Ξ132(i)=−F∑ρ=1λiρPp(i)ˉFp(ρ)ˉHp(i)+F∑ρ=1λiρPp(i)ˉKp(ρ)ˉLˉHp(i)+Rτ,Ξ133(i)=−Qτ−2Rτ+κ2ˉHTp(i)ˉHp(i),Ξ2(i)=[Ξ21(i),Ξ22(i)],Ξ21(i)=[0Pm(i)ˉEm(i)0Ξ211(i)0000κ2ˉHTm(i)ˉGm(i)00000000000000000000000000],Ξ22(i)=[−ˉdmˉAT(i)0−ˉτˉAT(i)0Ξ212(i)0Ξ213(i)0000000000ˉdpˉDT(i)0ˉτˉDT(i)00000000],Ξ211(i)=−F∑ρ=1λiρPm(i)ˉFm(ρ)ˉGm(i)+F∑ρ=1λiρPm(i)ˉKm(ρ)ˉLˉGm(i),Ξ212(i)=−ˉdmF∑ρ=1λiρˉHTm(i)ˉFTm(ρ)+ˉdmF∑ρ=1λiρˉHTm(i)ˉLTˉKTm(ρ),Ξ213(i)=−ˉτF∑ρ=1λiρˉHTm(i)ˉFTm(ρ)+ˉτF∑ρ=1λiρˉHTm(i)ˉLTˉKTm(ρ), $ |
and
$ Ξ3(i)=[Ξ31(i)Ξ32(i)∗Ξ33(i)],Ξ31(i)=[Ξ311(i)0Pp(i)ˉEp(i)0Ξ312(i)−κ2ˉHTp(i)000κ2ˉHTp(i)ˉGp(i)0000000000000−ˉGm(i)0κ2I−I000−ˉGp(i)∗−γI000∗∗−γI00∗∗∗κ2ˉGTm(i)ˉGm(i)−γI0∗∗∗∗κ2ˉGTp(i)ˉGp(i)−γI],Ξ311(i)=F∑ρ=1λiρPp(i)ˉFp(ρ)−F∑ρ=1λiρPp(i)ˉKp(ρ)ˉL,Ξ312(i)=−F∑ρ=1λiρPp(i)ˉFp(ρ)ˉGp(i)+F∑ρ=1λiρPp(i)ˉKp(ρ)ˉLˉGp(i),Ξ32(i)=[0−ˉdpˉCT(i)0−ˉτˉCT(i)0Ξ321(i)0Ξ322(i)0000ˉdmˉBT(i)0ˉτˉBT(i)0Ξ323(i)0Ξ324(i)00Ξ325(i)0Ξ326(i)ˉdmˉEm(i)0ˉτˉEm(i)00ˉdpˉETp(i)0ˉτˉETp(i)],Ξ321(i)=−ˉdpF∑ρ=1λiρˉHTp(i)ˉFTp(ρ)+ˉdpF∑ρ=1λiρˉHTp(i)ˉLTˉKTp(ρ),Ξ322(i)=−ˉτF∑ρ=1λiρˉHTp(i)ˉFTp(ρ)+ˉτF∑ρ=1λiρˉHTp(i)ˉLTˉKTp(ρ),Ξ323(i)=ˉdmF∑ρ=1λiρˉFTm(ρ)−ˉdmF∑ρ=1λiρˉLTˉKTm(ρ),Ξ324(i)=ˉτF∑ρ=1λiρˉFTm(ρ)−ˉτF∑ρ=1λiρˉLTˉKTm(ρ),Ξ325(i)=ˉdpF∑ρ=1λiρˉFTp(ρ)−ˉdpF∑ρ=1λiρˉLTˉKTp(ρ),Ξ326(i)=ˉτF∑ρ=1λiρˉFTp(ρ)−ˉτF∑ρ=1λiρˉLTˉKTp(ρ),Ξ33(i)=[Ξ331(i)0Ξ332(i)00Ξ333(i)0Ξ334(i)−R−1m000∗−R−1p00∗∗−R−1τ0∗∗∗−R−1τ],Ξ331(i)=−ˉdmF∑ρ=1λiρˉGTm(i)ˉFTm(ρ)+ˉdmF∑ρ=1λiρˉGTm(i)ˉLTˉKTm(ρ),Ξ332(i)=−ˉτF∑ρ=1λiρˉGTm(i)ˉFTm(ρ)+ˉτF∑ρ=1λiρˉGTm(i)ˉLTˉKTm(ρ),Ξ333(i)=−ˉdpF∑ρ=1λiρˉGTp(i)ˉFTp(ρ)+ˉdpF∑ρ=1λiρˉGTp(i)ˉLTˉKTp(ρ),Ξ334(i)=−ˉτF∑ρ=1λiρˉGTp(i)ˉFTp(ρ)+ˉτF∑ρ=1λiρˉGTp(i)ˉLTˉKTp(ρ). $ |
Proof. Construct the following mode-dependent Lyapunov-Krasovski functions:
$ V(i,t)=V1(i,t)+V2(i,t)+V3(i,t), $ | (3.1) |
where
$ V1(i,t)=eTm(t)Pm(i)em(t)+eTp(t)Pp(i)ep(t),V2(i,t)=∫tt−ˉdmeTm(s)Qmem(s)ds+∫tt−ˉdpeTp(s)Qpep(s)ds+∫tt−ˉτ(eTm(s)Qτem(s)+eTp(s)Qτep(s))ds,V3(i,t)=ˉdm∫0−ˉdm∫tt+φeTm(φ)Rmem(φ)dφds+ˉdp∫0−ˉdp∫tt+φeTp(φ)Rpep(φ)dφds+ˉτ∫0−ˉτ∫tt+φ(eTm(φ)Rτem(φ)+eTp(φ)Rτep(φ))dφds. $ |
Then, for all $ r(t) = i $ and $ \sigma (t) = \rho $, by defining the weak infinitesimal operator $ \mathcal{L} $ of $ V(i, t) $ as follows
$ LV(i,t)=limΔ→0+1Δ{E{V(r(t+Δ),t+Δ)|r(t)=i,t}−V(i,t)}, $ | (3.2) |
it can be deduced that
$ LV1(i,t)=E{2eTm(t)Pm(i)˙em(t)+2eTp(t)Pp(i)˙ep(t)+∑j∈NπijeTm(t)Pm(j)em(t)+∑j∈NπijeTp(t)Pp(i)ep(t)}=2eTm(t)[−Pm(i)ˉA(i)]em(t)+2eTm(t)[Pm(i)ˉB(i)]ˉg(ep(t−dp(t)))+2F∑ρ=1λiρeTm(t)[−Pm(i)(ˉFm(ρ)+ˉKm(ρ)ˉL)ˉHm(i)]em(t−τ(t))+2F∑ρ=1λiρeTm(t)[Pm(i)(ˉFm(ρ)−ˉKm(ρ)ˉL)]εm(t−τ(t))+2eTm(t)[Pm(i)ˉEm(i)]ω1(t)+2F∑ρ=1λiρeTm(t)[−Pm(i)(ˉFm(ρ)−ˉKm(ρ)ˉL)ˉGm(i)]ω3(t−τ(t))+2eTp(t)[−Pp(i)ˉC(i)]ep(t)+2eTp(t)[Pp(i)ˉD(i)]em(t−dm(t))+2eTp(t)[Pp(i)ˉEp(i)]ω2(t)+2F∑ρ=1λiρeTp(t)[−Pp(i)(ˉFp(ρ)−ˉKp(ρ)ˉL)ˉHp(i)]ep(t−τ(t))+2F∑ρ=1λiρeTp(t)[Pp(i)(ˉFp(ρ)−ˉKp(ρ)ˉL)]εp(t−τ(t))+2F∑ρ=1λiρeTp(t)[−Pp(i)(ˉFp(ρ)−ˉKp(ρ)ˉL)ˉGp(i)]ω4(t−τ(t))+∑j∈NπijeTm(t)Pm(j)em(t)+∑j∈NπijeTp(t)Pp(i)ep(t). $ | (3.3) |
Similarly, it can be deduced that
$ LV2(i,t)=eTm(t)Qmem(t)−eTm(t−ˉdm)Qmem(t−ˉdm)+eTp(t)Qpep(t)−eTp(t−ˉdp)Qpep(t−ˉdp)+eTm(t)Qτem(t)+eTp(t)Qτep(t)−eTm(t−ˉτ)Qτem(t−ˉτ)−eTp(t−ˉτ)Qτep(t−ˉτ)=eTm(t)(Qm+Qτ)em(t)+eTp(t)(Qp+Qτ)ep(t)−eTm(t−ˉdm)Qmem(t−ˉdm)−eTp(t−ˉdp)Qmep(t−ˉdp)−eTm(t−ˉτ)Qτem(t−ˉτ)−eTp(t−ˉτ)Qτep(t−ˉτ), $ | (3.4) |
and
$ LV3(i,t)=E{ˉd2m˙eTm(t)Rm˙em(t)+ˉd2p˙eTp(t)Rp˙ep(t)+ˉτ2˙eTm(t)Rτ˙em(t)+ˉτ2˙eTp(t)Rτ˙ep(t)−ˉdm∫tt−ˉdm˙eTm(φ)Rm˙em(φ)dφ−ˉdp∫tt−ˉdp˙eTp(φ)Rp˙ep(φ)dφ−ˉτ∫tt−ˉτ˙eTm(φ)Rτ˙em(φ)dφ−ˉτ∫tt−ˉτ˙eTp(φ)Rτ˙ep(φ)dφ}=ψT(t)(ϝmRmϝTm+ϝpRpϝTp+ϝτRτϝTτ)ψ(t)−ˉdm∫tt−ˉdm˙eTm(φ)Rm˙em(φ)dφ−ˉdp∫tt−ˉdp˙eTp(φ)Rp˙ep(φ)dφ−ˉτ∫tt−ˉτ˙eTm(φ)Rτ˙em(φ)dφ−ˉτ∫tt−ˉτ˙eTp(φ)Rτ˙ep(φ)dφ, $ | (3.5) |
where $ \psi ^{T}(t) = [e_{m}^{T}(t), e_{m}^{T}(t-\tau (t)), e_{m}^{T}(t-\bar{\tau }), e_{m}^{T}\left(t-\bar{d}_{m}\right), e_{m}^{T}\left(t-d_{m}(t)\right), e_{p}^{T}\left(t-\bar{d}_{p}\right), e_{p}^{T}\left(t-d_{p}(t)\right), e_{p}^{T}(t), e_{p}^{T}(t-\tau (t)), e_{p}^{T}(t-\bar{\tau}), \bar{g} ^{T}\left(e_{p}\left(t-d_{p}(t)\right) \right), \varepsilon _{m}^{T}(t-\tau (t)), \varepsilon _{p}^{T}(t-\tau (t)), \omega _{1}^{T}(t), \omega _{2}^{T}(t), \omega _{3}^{T}(t-\tau (t)), \omega _{4}^{T}(t-\tau (t))]^{T} $ and
$ ϝm=[−ˉdmˉAT(i)−ˉdm∑Fρ=1λiρˉHTm(i)(ˉFTm(ρ)−ˉLTˉKTm(ρ))00000000ˉdmˉBT(i)ˉdm∑Fρ=1λiρ(ˉFTm(ρ)−ˉLTˉKTm(ρ))0ˉdmˉEm(i)0−ˉdm∑Fρ=1λiρˉGTm(i)(ˉFTm(ρ)−ˉLTˉKTm(ρ))0],ϝp=[0000ˉdpˉDT(i)00−ˉdpˉCT(i)−ˉdp∑Fρ=1λiρˉHTp(i)(ˉFTp(ρ)−ˉLTˉKTp(ρ))000ˉdp∑Fρ=1λiρ(ˉFTp(ρ)−ˉLTˉKTp(ρ))0ˉdpˉETp(i)0−ˉdp∑Fρ=1λiρˉGTp(i)(ˉFTp(ρ)−ˉLTˉKTp(ρ))], $ |
$ ϝτ=[−ˉτˉAT(i)−ˉτ∑Fρ=1λiρˉHTm(i)(ˉFTm(ρ)−ˉLTˉKTm(ρ))00000000ˉτˉBT(i)ˉτ∑Fρ=1λiρ(ˉFTm(ρ)−ˉLTˉKTm(ρ))0ˉτˉEm(i)0−ˉτ∑Fρ=1λiρˉGTm(i)(ˉFTm(ρ)−ˉLTˉKTm(ρ))0]+[0000ˉτˉDT(i)00−ˉτˉCT(i)−ˉτ∑Fρ=1λiρˉHTp(i)(ˉFTp(ρ)−ˉLTˉKTp(ρ))000ˉτ∑Fρ=1λiρ(ˉFTp(ρ)−ˉLTˉKTp(ρ))0ˉτˉETp(i)0−ˉτ∑Fρ=1λiρˉGTp(i)(ˉFTp(ρ)−ˉLTˉKTp(ρ))]. $ |
In the sequel, by virtue of Lemma 1, one can obtain that
$ −ˉdm∫tt−ˉdm˙eTm(φ)Rm˙em(φ)dφ≤[eTm(t)eTm(t−dm(t))eTm(t−ˉdm)]T[−RmRm0∗−2RmRm∗∗−Rm][em(t)em(t−dm(t))em(t−ˉdm)], $ | (3.6) |
$ −ˉdp∫tt−ˉdp˙eTp(φ)Rp˙ep(φ)dφ≤[eTp(t)eTp(t−dp(t))eTp(t−ˉdp)]T[−RpRp0∗−2RpRp∗∗−Rp][ep(t)ep(t−dp(t))ep(t−ˉdp)] $ | (3.7) |
$ −ˉτ∫tt−ˉτ˙eTm(φ)Rτ˙em(φ)dφ≤[eTm(t)eTm(t−τ(t))eTm(t−ˉτ)]T[−RτRτ0∗−2RτRτ∗∗−Rτ][em(t)em(t−τ(t))em(t−ˉτ)], $ | (3.8) |
$ −ˉτ∫tt−ˉτ˙eTp(φ)Rτ˙ep(φ)dφ≤[eTp(t)eTp(t−τ(t))eTp(t−ˉτ)]T[−RτRτ0∗−2RτRτ∗∗−Rτ][ep(t)ep(t−τ(t))ep(t−ˉτ)]. $ | (3.9) |
Then, it can be deduced that
$ LV(i,t)≤ηT(t)Λη(t) $ |
where $ \eta ^{T}\left(t\right) = [\eta ^{T}_{1}\left(t\right), \eta ^{T} _{2}\left(t\right), \eta ^{T} _{3}\left(t\right), \eta ^{T} _{4}\left(t\right) ] $ with
$ ηT1(t)=[eTm(t)eTm(t−τ(t))eTm(t−ˉτ)eTm(t−ˉdm)eTm(t−dm(t))],ηT2(t)=[eTp(t−ˉdp)eTp(t−dp(t))eTp(t)eTp(t−τ(t))eTp(t−ˉτ)],ηT3(t)=[ˉgT(ep(t−dp(t)))εTm(t−τ(t))εTp(t−τ(t))],ηT4(t)=[ωT1(t)ωT2(t)ωT3(t−τ(t))ωT4(t−τ(t))], $ |
and $ \Lambda = \Lambda _{1}+\Lambda _{2}R_{m}\Lambda _{2}^{T}+\Lambda _{3}R_{p}\Lambda _{3}^{T}+\Lambda _{4}R_{\tau }\Lambda _{4}^{T}+\Lambda _{5}R_{\tau }\Lambda _{5}^{T} $ with
$ Λ1=[Λ11Λ12∗Λ22],Λ11=[Λ111Λ11200Rm00∗−2RτRτ0000∗∗−Qτ−Rτ0000∗∗∗−Qm−RmRTm00∗∗∗∗−2Rm00∗∗∗∗∗−Qm−RpRp∗∗∗∗∗∗−2Rp],Λ111=−2Pm(i)ˉA(i)+Qm+Qτ−Rm−Rτ+∑j∈NπijPm(j),Λ112=−F∑ρ=1λiρPm(i)ˉFm(ρ)ˉHm(i)+F∑ρ=1λiρPm(i)ˉKm(ρ)ˉLˉHm(i)+RτΛ12=[000Pm(i)ˉB(i)Λ1210Pm(i)ˉEm(i)0Λ1220000000000000000000000000000000ˉDT(i)Pp(i)0000000000000000000Rp000000000],Λ121=F∑ρ=1λiρPm(i)ˉFm(ρ)−F∑ρ=1λiρPm(i)ˉKm(ρ)ˉLΛ122=−F∑ρ=1λiρPm(i)ˉFm(ρ)ˉGm(i)+F∑ρ=1λiρPm(i)ˉKm(ρ)ˉLˉGm(i)Λ22=[Λ221Λ222000Λ2230Pp(i)ˉEp(i)0Λ224∗−Qτ−2RτRτ0000000∗∗−Rτ0000000∗∗∗0000000∗∗∗∗000000∗∗∗∗∗00000∗∗∗∗∗∗0000∗∗∗∗∗∗∗000∗∗∗∗∗∗∗∗00∗∗∗∗∗∗∗∗∗0],Λ221=−2Pp(i)ˉC(i)+Qp+Qτ−Rp−Rτ+∑j∈NπijPp(i),Λ222=−F∑ρ=1λiρPp(i)ˉFp(ρ)ˉHp(i)+F∑ρ=1λiρPp(i)ˉKp(ρ)ˉLˉHp(i)+RτΛ223=F∑ρ=1λiρPp(i)ˉFp(ρ)−F∑ρ=1λiρPp(i)ˉKp(ρ)ˉLΛ224=−F∑ρ=1λiρPp(i)ˉFp(ρ)ˉGp(i)+F∑ρ=1λiρPp(i)ˉKp(ρ)ˉLˉGp(i) $ |
$ Λ2=[−ˉdmˉAT(i)−ˉdm∑Fρ=1λiρˉHTm(i)(ˉFTm(ρ)−ˉLTˉKTm(ρ))00000000ˉdmˉBT(i)ˉdm∑Fρ=1λiρ(ˉFTm(ρ)−ˉLTˉKTm(ρ))0ˉdmˉEm(i)0−ˉdm∑Fρ=1λiρˉGTm(i)(ˉFTm(ρ)−ˉLTˉKTm(ρ))0],Λ3=[0000ˉdpˉDT(i)00−ˉdpˉCT(i)−ˉdp∑Fρ=1λiρˉHTp(i)(ˉFTp(ρ)−ˉLTˉKTp(ρ))000ˉdp∑Fρ=1λiρ(ˉFTp(ρ)−ˉLTˉKTp(ρ))0ˉdpˉETp(i)0−ˉdp∑Fρ=1λiρˉGTp(i)(ˉFTp(ρ)−ˉLTˉKTp(ρ))],Λ4=[−ˉτˉAT(i)−ˉτ∑Fρ=1λiρˉHTm(i)(ˉFTm(ρ)−ˉLTˉKTm(ρ))0000000ˉτˉBT(i)ˉτ∑Fρ=1λiρ(ˉFTm(ρ)−ˉLTˉKTm(ρ))0ˉτˉEm(i)0−ˉτ∑Fρ=1λiρˉGTm(i)(ˉFTm(ρ)−ˉLTˉKTm(ρ))0],Λ5=[0000ˉτˉDT(i)0−ˉτˉCT(i)−ˉτ∑Fρ=1λiρˉHTp(i)(ˉFTp(ρ)−ˉLTˉKTp(ρ))000ˉτ∑Fρ=1λiρ(ˉFTp(ρ)−ˉLTˉKTp(ρ))0ˉτˉETp(i)0−ˉτ∑Fρ=1λiρˉGTp(i)(ˉFTp(ρ)−ˉLTˉKTp(ρ))]. $ |
Furthermore, considering the nonlinear properties for the GRNs, it yields that
$ [ˉg(ep(t−dp(t)))ep(t−dp(t))]T[−2Iςς0][ˉg(ep(t−dp(t)))ep(t−dp(t))]≥0. $ | (3.10) |
In addition, by recalling the event-triggering conditions, one has
$ N∑l=1κ2[˜yl(ιk)−εl(ιk)]T[˜yl(ιk)−εl(ιk)]−N∑l=1[εl(ιk)]Tεl(ιk)≥0, $ | (3.11) |
where $ \tilde{y}^{l}(\iota _{k}) = [(\tilde{y}_{m}^{l}(t))^{T}, (\tilde{y} _{p}^{l}(t))^{T}]^{T} $, $ \varepsilon ^{l}(\iota _{k}) = [(\varepsilon _{m}^{l}(\iota _{k}))^{T}, (\varepsilon _{p}^{l}(\iota _{k}))^{T}]^{T} $ and it can be further verified that
$ N∑l=1κ2[Hlm(i)elm(ιk)+Glm(i)ω3(ιk)−εlm(ιk)]T[Hlm(i)elm(ιk)+Glm(i)ω3(ιk)−εlm(ιk)]+N∑l=1κ2[Hlp(i)elp(ιk)+Glp(i)ω4(ιk)−εlp(ιk)]T[Hlp(i)elp(ιk)+Glp(i)ω4(ιk)−εlp(ιk)]−(εm(t−τ(t)))Tεm(t−τ(t))−(εp(t−τ(t)))Tεp(t−τ(t))≥0. $ | (3.12) |
As a consequence, it can be derived by Schur complement that the modified $ H_{\infty} $ performance can be achieved under the mean-square sense when $ \Xi(i) < 0 $ is satisfied according to Theorem 1, which implies that
$ LV(i,t)+eTm(t)ˉUTm(i)ˉUm(i)em(t)+eTp(t)ˉUTp(i)ˉUp(i)ep(t)−γωT1(t)ω1(t)−γωT2(t)ω2(t)−γωT3(t−τ(t))ω3(t−τ(t))−γωT4(t−τ(t))ω4(t−τ(t))<0 $ | (3.13) |
and therefore completes the proof.
Remark 3. The above established sufficient convex optimization conditions are derived with aid of linear matrix inequality technique, which can be easily solved by Matlab LMI toolbox or other mathematical softwares. It should be pointed out that the matrix dimension dependents on both system mode information and the sensor networks structure, which means that the more system modes and sensors are involved, the more computational consumption is needed. As a result, the distributed state estimators should be well designed for practical system implements with considerable computational efficiency.
Theorem 2. For given $ \bar{d}_{p} $, $ \bar{d}_{m} $, $ \bar{\tau} $ and the modified $ H_{\infty } $ filtering can be satisfied over distributed sensor networks according to Definition 1, if there exist mode-dependent matrices $ P_{m}(i) > 0 $, $ P_{p}(i) > 0 $, $ X_{m}(i) > 0 $, $ X_{p}(i) > 0 $, $ Y_{F, m}\left(\rho \right) $, $ Y_{F, p}\left(\rho \right) $, $ Y_{K, m}\left(\rho \right) $, $ Y_{K, p}\left(\rho \right), $ and matrices $ Q_{m} > 0 $, $ Q_{p} > 0 $, $ Q_{\tau } > 0 $, $ R_{m} > 0 $, $ R_{p} > 0 $, $ R_{\tau } > 0 $, such that $ \tilde{\Xi} _{i} < 0 $, for all $ i\in \mathcal{N} $, $ k = 1, 2, \ldots, N $ and $ \rho \in \mathcal{F} $, $ k = 1, 2, \ldots, F $, where
$ ˜Ξ(i)=[˜Ξ1(i)˜Ξ2(i)∗˜Ξ3(i)],˜Ξ1(i)=[˜Ξ11(i)˜Ξ12(i)∗˜Ξ13(i)],˜Ξ11(i)=[˜Ξ111(i)˜Ξ112(i)∗˜Ξ113(i)],˜Ξ111(i)=[˜Ξ1111(i)˜Ξ1112(i)00Rm00∗˜Ξ1113(i)Rτ0000∗∗−Qτ−Rτ0000∗∗∗−Qm−RmRTm00∗∗∗∗−2Rm00∗∗∗∗∗−Qp−RpRp∗∗∗∗∗∗−2Rp],˜Ξ112(i)=[000Pm(i)ˉB(i)˜Ξ1121(i)00000−κ2ˉHTm(i)0000000000000ˉDT(i)Pp(i)00000000000Rp00ς00],˜Ξ113(i)=[˜Ξ1131(i)˜Ξ1132(i)000˜Ξ1133(i)∗˜Ξ1134(i)Rτ00−κ2ˉHTp(i)∗∗−Rτ000∗∗∗−2I00∗∗∗∗κ2I−I0∗∗∗∗∗κ2I−I],˜Ξ1111(i)=−2Pm(i)ˉA(i)+Qm+Qτ−Rm−Rτ+ˉUTm(i)ˉUm(i)+∑j∈NπijPm(j),˜Ξ1112(i)=−F∑ρ=1λiρYF,m(ρ)ˉHm(i)+F∑ρ=1λiρYK,m(ρ)ˉLˉHm(i)+Rτ,˜Ξ1113(i)=−2Rτ+κ2ˉHTm(i)ˉHm(i),˜Ξ1121(i)=F∑ρ=1λiρYF,m(ρ)−F∑ρ=1λiρYK,m(ρ)ˉL,˜Ξ1131(i)=−2Pp(i)ˉC(i)+Qp+Qτ−Rp−Rτ+ˉUTp(i)ˉUp(i)+∑j∈NπijPp(i),˜Ξ1132(i)=−F∑ρ=1λiρYF,p(ρ)ˉHp(i)+F∑ρ=1λiρYK,p(ρ)ˉLˉHp(i)+Rτ,˜Ξ1133(i)=F∑ρ=1λiρYF,p(ρ)−F∑ρ=1λiρYK,p(ρ)ˉL,˜Ξ1134(i)=−Qτ−2Rτ+κ2ˉHTp(i)ˉHp(i), $ |
$ ˜Ξ12(i)=[˜Ξ121(i)˜Ξ122(i)˜Ξ123(i)˜Ξ124(i)],˜Ξ121(i)=[Pm(i)ˉEm(i)0˜Ξ1211(i)000κ2ˉHTm(i)ˉGm(i)0000000000000000000000Pp(i)ˉEp(i)0˜Ξ1212(i)],˜Ξ1211(i)=−F∑ρ=1λiρYF,m(ρ)ˉGm(i)+F∑ρ=1λiρYK,m(ρ)ˉLˉGm(i),˜Ξ1212(i)=−F∑ρ=1λiρYF,p(ρ)ˉGp(i)+F∑ρ=1λiρYK,p(ρ)ˉLˉGp(i),˜Ξ122(i)=[−ˉdmˉAT(i)Pm(i)0−ˉτˉAT(i)Pm(i)0˜Ξ1221(i)0˜Ξ1222(i)0000000000ˉdpˉDT(i)Pp(i)0ˉτˉDT(i)Pp(i)000000000−ˉdpˉCT(i)Pp(i)0−ˉτˉCT(i)Pp(i)],˜Ξ1221(i)=−ˉdmF∑ρ=1λiρˉHTm(i)YTF,m(ρ)+ˉdmF∑ρ=1λiρˉHTm(i)ˉLTYTK,m(ρ),˜Ξ1222(i)=−ˉτF∑ρ=1λiρˉHTm(i)YTF,m(ρ)+ˉτF∑ρ=1λiρˉHTm(i)ˉLTYTK,m(ρ),˜Ξ123(i)=[000κ2ˉHTp(i)ˉGp(i)0000000000ˉdmˉBT(i)Pm(i)00−ˉGm(i)0˜Ξ1231(i)000−ˉGp(i)0],˜Ξ1231(i)=ˉdmF∑ρ=1λiρYTF,m(ρ)−ˉdmF∑ρ=1λiρˉLTYTK,m(ρ),˜Ξ124(i)=[˜Ξ1241(i)0˜Ξ1242(i)0000ˉτˉBT(i)Pm(i)00˜Ξ1243(i)0˜Ξ1244(i)0˜Ξ1245(i)],˜Ξ1241(i)=−ˉdpF∑ρ=1λiρˉHTp(i)YTF,p(ρ)+ˉdpF∑ρ=1λiρˉHTp(i)ˉLTYTK,p(ρ),˜Ξ1242(i)=−ˉτF∑ρ=1λiρˉHTp(i)YTF,p(ρ)+ˉτF∑ρ=1λiρˉHTp(i)ˉLTYTK,p(ρ),˜Ξ1243(i)=ˉτF∑ρ=1λiρYTF,m(ρ)−ˉτF∑ρ=1λiρˉLTYTK,m(ρ),˜Ξ1244(i)=ˉdpF∑ρ=1λiρYTF,p(ρ)−ˉdpF∑ρ=1λiρˉLTYTK,p(ρ),˜Ξ1245(i)=ˉτF∑ρ=1λiρYTF,p(ρ)−ˉτF∑ρ=1λiρˉLTYTK,p(ρ), $ |
$ ˜Ξ13(i)=[˜Ξ131(i)˜Ξ132(i)∗˜Ξ133(i)],˜Ξ131(i)=[−γI000∗−γI00∗∗κ2ˉGTm(i)ˉGm(i)−γI0∗∗∗κ2ˉGTp(i)ˉGp(i)−γI],˜Ξ132(i)=[ˉdmˉEm(i)Pm(i)0ˉτˉEm(i)Pm(i)00ˉdpˉETp(i)Pp(i)0ˉτˉETp(i)Pp(i)˜Ξ1311(i)0˜Ξ1312(i)00˜Ξ1313(i)0˜Ξ1314(i)],˜Ξ133(i)=[Rm−2Pm(i)000∗Rm−2Pp(i)00∗∗Rτ−2Pm(i)0∗∗∗Rτ−2Pp(i)],˜Ξ1311(i)=−ˉdmF∑ρ=1λiρˉGTm(i)YTF,m(ρ)+ˉdmF∑ρ=1λiρˉGTm(i)ˉLTYTK,m(ρ),˜Ξ1312(i)=−ˉτF∑ρ=1λiρˉGTm(i)YTF,m(ρ)+ˉτF∑ρ=1λiρˉGTm(i)ˉLTYTK,m(ρ),˜Ξ1313(i)=−ˉdpF∑ρ=1λiρˉGTp(i)YTF,p(ρ)+ˉdpF∑ρ=1λiρˉGTp(i)ˉLTYTK,p(ρ),˜Ξ1314(i)=−ˉτF∑ρ=1λiρˉGTp(i)YTF,p(ρ)+ˉτF∑ρ=1λiρˉGTp(i)ˉLTYTK,p(ρ),˜Ξ2(i)=[˜Ξ21(i)˜Ξ22(i)˜Ξ23(i)˜Ξ24(i)],˜Ξ21(i)=[˜Ξ211(i),˜Ξ212(i)],˜Ξ211(i)=[λi1(Xm(1)−Pm(i))ϵmˉHTm(i)(YTF,m(1)−ˉLTYTK,m(1))0],˜Ξ212(i)=[λi2(Xm(2)−Pm(i))⋯λiF(Xm(F)−Pm(i))ϵmˉHTm(i)(YTF,m(2)−ˉLTYTK,m(2))⋯ϵmˉHTm(i)(YTF,m(F)−ˉLTYTK,m(F))0⋯0],˜Ξ22(i)=[˜Ξ221(i),˜Ξ222(i)],˜Ξ221(i)=[0λi1(Xp(1)−Pp(i))ϵpˉHTp(i)(YTF,p(1)−ˉLTYTK,p(1))0],˜Ξ222(i)=[0⋯0λi2(Xp(2)−Pp(i))⋯λiF(Xp(F)−Pp(i))ϵpˉHTp(i)(YTF,p(2)−ˉLTYTK,p(2))⋯ϵpˉHTp(i)(YTF,p(F)−ˉLTYTK,p(F))0⋯0],˜Ξ23(i)=[˜Ξ231(i),˜Ξ232(i)],˜Ξ231(i)=[ϵm(ˉLTYTK,m(1)−YTF,m(1))000ϵmˉGTm(i)(YTF,m(1)−ˉLTYTK,m(1))0ˉdmλi1(Xm(1)−Pm(i))0ˉdmλi1(Xm(1)−Pm(i))0],˜Ξ232(i)=[ϵm(ˉLTYTK,m(2)−YTF,m(2))⋯ϵm(ˉLTYTK,m(F)−YTF,m(F))0⋯00⋯00⋯0ϵmˉGTm(i)(YTF,m(2)−ˉLTYTK,m(2))⋯ϵTmˉGTm(i)(YTF,m(F)−ˉLTYTK,m(F))0⋯0ˉdmλi2(Xm(2)−Pm(i))⋯ˉdmλiF(Xm(F)−Pm(i))0⋯0ˉdmλi2(Xm(2)−Pm(i))⋯ˉdmλiF(Xm(F)−Pm(i))0⋯0],˜Ξ24(i)=[˜Ξ241(i),˜Ξ242(i)],˜Ξ241(i)=[0ϵp(ˉLTYTK,p(1)−YTF,p(1))000ϵpˉGTp(i)(YTF,p(1)−ˉLTYTK,p(1))0ˉdpλi1(Xp(1)−Pp(i))0ˉdpλi1(Xp(1)−Pp(i))],˜Ξ242(i)=[0⋯0ϵp(ˉLTYTK,p(2)−YTF,p(2))⋯ϵp(ˉLTYTK,p(F)−YTF,p(F))0⋯00⋯00⋯0ϵpˉGTp(i)(YTF,p(2)−ˉLTYTK,p(2))⋯ϵTpˉGTp(i)(YTF,p(F)−ˉLTYTK,p(F))0⋯0ˉdpλi2(Xp(2)−Pp(i))⋯ˉdpλiF(Xp(F)−Pp(i))0⋯0ˉdpλi2(Xp(2)−Pp(i))⋯ˉdpλiF(Xp(F)−Pp(i))], $ |
and
$ ˜Ξ3(i)=[˜Ξ31(i)˜Ξ32(i)∗˜Ξ33(i)], $ | (3.14) |
$ ˜Ξ31(i)=[−2ϵmXm(1)0⋯0∗−2ϵmXm(2)⋯0⋮⋮⋱⋮∗∗⋯−2ϵmXm(F)] $ | (3.15) |
$ \begin{align} \tilde{\Xi}_{32}(i) = &0, \end{align} $ | (3.16) |
$ \begin{align} \tilde{\Xi}_{33}(i) = &\left[ \begin{array}{cccc} -2\epsilon _{p}X_{p}(1) & 0 & \cdots & 0 \\ \ast & -2\epsilon _{p}X_{p}(2) & \cdots & 0 \\ \vdots & \vdots & \ddots & \vdots \\ \ast & \ast & \ast & -2\epsilon _{p}X_{p}(F) \end{array} \right] . \end{align} $ | (3.17) |
With the aforementioned feasible solutions, the desired mode-dependent filter gains $ F_{m}^{l}(\sigma (t)), F_{p}^{l}(\sigma (t)) $ and $ K_{m}^{l}(\sigma (t)), K_{p}^{l}(\sigma (t)) $ can be obtained by the following calculation:
$ \begin{align*} Y_{F,m}\left( \rho \right) X_{m}^{-1}\left( \rho \right) = &\bar{F} _{m}(\rho ), \\ Y_{F,p}\left( \rho \right) X_{p}^{-1}\left( \rho \right) = &\bar{F} _{p}(\rho ), \\ Y_{K,m}\left( \rho \right) X_{m}^{-1}\left( \rho \right) = &\bar{K} _{m}(\rho ), \\ Y_{K,p}\left( \rho \right) X_{p}^{-1}\left( \rho \right) = &\bar{K} _{p}(\rho ). \end{align*} $ |
Proof. By pre- and post-multiplying $ diag\{I, I, I, I, I, I, I, I, I, I, I, P_{m}\left(i\right), P_{p}\left(i\right), P_{m}\left(i\right), P_{p}\left(i\right) \} $ to $ \Xi(i) < 0 $, Lemma 2 is employed accordingly.
Eventually, define $ Y_{F, m}\left(\rho \right) = X_{m}\left(\rho \right) \bar{F}_{m}(\rho) $, $ Y_{F, p}\left(\rho \right) = X_{p}\left(\rho \right) \bar{F}_{p}(\rho) $, $ Y_{K, m}\left(\rho \right) = X_{m}\left(\rho \right) \bar{K}_{m}(\rho) $, $ Y_{K, p}\left(\rho \right) = X_{p}\left(\rho \right) \bar{K}_{p}(\rho) $, then the remainder of proof will follow directly from Theorem 1.
Remark 4. The proposed modified $ H_{\infty} $ performance is more applicable for the distributed sensor networks design. In addition, the performance index $ \gamma $ can be further optimized by solving the following convex optimization problem:
$ \begin{align} &{\rm{minimize }}\; \gamma , \\ & {\rm{subject\; to } }\;\tilde{\Xi}_{i} < 0, i\in \mathcal{N}, \rho \in \mathcal{F}. \end{align} $ | (3.18) |
This section is dedicated to validate the applicability of our proposed distributed nonsynchronous designs, while illustrative studies have been performed with simulation results.
Consider the following HMJP GRNs model with a set of biological simulated parameters as
$ \begin{align*} A\left( 1\right) = &\left[ \begin{array}{cc} 1 & 0 \\ 0 & 1 \end{array} \right] ,A\left( 2\right) = \left[ \begin{array}{cc} 1.3 & 0 \\ 0 & 1.1 \end{array} \right] , \\ B\left( 1\right) = &\left[ \begin{array}{cc} 0.2 & 0.1 \\ 0.8 & 0.3 \end{array} \right] ,B\left( 2\right) = \left[ \begin{array}{cc} 0.3 & 0.2 \\ 0.1 & 0.2 \end{array} \right] , \\ C\left( 1\right) = &\left[ \begin{array}{cc} 1 & 0 \\ 0 & 1.2 \end{array} \right] ,C\left( 2\right) = \left[ \begin{array}{cc} 1 & 0 \\ 0 & 1 \end{array} \right] , \\ D\left( 1\right) = &\left[ \begin{array}{cc} 0.8 & 0 \\ 0 & 0.8 \end{array} \right] ,D\left( 2\right) = \left[ \begin{array}{cc} 1.1 & 0 \\ 0 & 1.3 \end{array} \right] , \\ E_{m}\left( 1\right) = &\left[ \begin{array}{l} 1.1 \\ 1 \end{array} \right] ,E_{m}\left( 2\right) = \left[ \begin{array}{l} 1.2 \\ 1.1 \end{array} \right] , \\ E_{p}\left( 1\right) = &\left[ \begin{array}{l} 1 \\ 1.1 \end{array} \right] ,E_{p}\left( 2\right) = \left[ \begin{array}{l} 1.2 \\ 1 \end{array} \right] , \\ U_{m}\left( 1\right) = &\left[ \begin{array}{cc} 1 & 0 \\ 0 & 1 \end{array} \right] ,U_{m}\left( 2\right) = \left[ \begin{array}{cc} 1.02 & 0 \\ 0 & 1.02 \end{array} \right] , \\ U_{p}\left( 1\right) = &\left[ \begin{array}{cc} 1 & 0 \\ 0 & 1 \end{array} \right] ,U_{p}\left( 2\right) = \left[ \begin{array}{cc} 1.01 & 0 \\ 0 & 1.01 \end{array} \right] , \end{align*} $ |
where the transition rate matrix with two modes of hidden Marko chain are assumed by
$ \begin{equation*} \Pi = \left[ \begin{array}{cc} -0.2 & 0.2 \\ 0.8 & -0.8 \end{array} \right] . \end{equation*} $ |
For the sensor networks, it is assumed that two networked sensors are deployed for mRNAs as well as proteins, such that the communication topology matrix for each network is given by
$ \begin{equation*} \mathbb{L} = \left[ \begin{array}{cc} 1 & -1 \\ -1 & 1 \end{array} \right] . \end{equation*} $ |
Moreover, the mode-dependent sensor parameters are set as follows:
$ \begin{align*} H_{m}^{1}\left( 1\right) & = \left[ \begin{array}{cc} 1 & 0 \\ 0 & 1 \end{array} \right] ,H_{m}^{1}\left( 2\right) = \left[ \begin{array}{cc} 1.1 & 0 \\ 0 & 1.1 \end{array} \right] , \\ H_{m}^{2}\left( 1\right) & = \left[ \begin{array}{cc} 1 & 0 \\ 0 & 1 \end{array} \right] ,H_{m}^{2}\left( 2\right) = \left[ \begin{array}{cc} 1.2 & 0 \\ 0 & 1.1 \end{array} \right] , \\ H_{p}^{1}\left( 1\right) & = \left[ \begin{array}{cc} 1 & 0 \\ 0 & 1 \end{array} \right] ,H_{p}^{1}\left( 2\right) = \left[ \begin{array}{cc} 0.9 & 0 \\ 0 & 1.1 \end{array} \right] , \\ H_{p}^{2}\left( 1\right) & = \left[ \begin{array}{cc} 1 & 0 \\ 0 & 1 \end{array} \right] ,H_{p}^{2}\left( 2\right) = \left[ \begin{array}{cc} 1 & 0 \\ 0 & 1.1 \end{array} \right] , \\ G_{m}^{1}\left( 1\right) & = \left[ \begin{array}{l} 1 \\ 1 \end{array} \right] ,G_{m}^{1}\left( 2\right) = \left[ \begin{array}{l} 1 \\ 1.1 \end{array} \right] , \\ G_{m}^{2}\left( 1\right) & = \left[ \begin{array}{l} 1 \\ 1 \end{array} \right] ,G_{m}^{2}\left( 2\right) = \left[ \begin{array}{l} 1.1 \\ 1 \end{array} \right] , \\ G_{p}^{1}\left( 1\right) & = \left[ \begin{array}{l} 1 \\ 1 \end{array} \right] ,G_{p}^{1}\left( 2\right) = \left[ \begin{array}{l} 0.9 \\ 1 \end{array} \right] , \\ G_{p}^{2}\left( 1\right) & = \left[ \begin{array}{l} 1 \\ 1 \end{array} \right] ,G_{p}^{2}\left( 2\right) = \left[ \begin{array}{l} 1 \\ 1.2 \end{array} \right] . \end{align*} $ |
In the simulation example, the sampling interval of networked sensors is set by $ h = 0.1 $s. Furthermore, the probability matrix of observed modes is supposed to be
$ \begin{equation*} \Lambda = \left[ \begin{array}{cc} 0.7 & 0.3 \\ 0.2 & 0.8 \end{array} \right] , \end{equation*} $ |
which indicates that the nonsynchronous mode information is utilized instead of true system modes.
Consequently, with these simulation parameters, the optimized nonsynchronous mode-dependent filter gains are solved by Theorem 2 as follows:
$ \begin{align*} F_{m}^{1}\left( 1\right) & = \left[ \begin{array}{ll} 0.0071 & 0.0153 \\ -0.0353 & 0.0983 \end{array} \right] ,F_{m}^{1}\left( 2\right) = \left[ \begin{array}{ll} 0.0063 & 0.0147\\ -0.0387 & 0.1017 \end{array} \right] , \\ F_{m}^{2}\left( 1\right) & = \left[ \begin{array}{ll} 0.0130 & 0.0075\\ -0.0256 & 0.0718 \end{array} \right] ,F_{m}^{2}\left( 2\right) = \left[ \begin{array}{ll} 0.0141 & 0.0061\\ -0.0295 & 0.0782 \end{array} \right] , \\ F_{p}^{1}\left( 1\right) & = \left[ \begin{array}{ll} 0.0116 & 0.0051\\ -0.0122 & 0.0406 \end{array} \right] ,F_{p}^{1}\left( 2\right) = \left[ \begin{array}{ll} 0.0107 & 0.0088\\ -0.0133 & 0.0464 \end{array} \right] , \\ F_{p}^{2}\left( 1\right) & = \left[ \begin{array}{ll} 0.0106 & 0.0050\\ -0.0151 & 0.0461 \end{array} \right] ,F_{p}^{2}\left( 2\right) = \left[ \begin{array}{ll} 0.0089 & 0.0044\\ -0.0208 & 0.0444 \end{array} \right] , \\ K_{m}^{1}\left( 1\right) & = \left[ \begin{array}{ll} -0.0134 & -0.0113 \\ -0.0387 & -0.0323 \end{array} \right] ,K_{m}^{1}\left( 2\right) = \left[ \begin{array}{ll} -0.0135 & -0.0127\\ -0.0378 & -0.0351 \end{array} \right] , \\ K_{m}^{2}\left( 1\right) & = \left[ \begin{array}{ll} -0.0118 & -0.0082\\ -0.0263 & -0.0177 \end{array} \right] ,K_{m}^{2}\left( 2\right) = \left[ \begin{array}{ll} -0.0135 & -0.0115\\ -0.0283 & -0.0232 \end{array} \right] , \\ K_{p}^{1}\left( 1\right) & = \left[ \begin{array}{ll} -0.0080 & -0.0079\\ -0.0134 & -0.0149 \end{array} \right] ,K_{p}^{1}\left( 2\right) = \left[ \begin{array}{ll} -0.0079 & -0.0075\\ -0.0145 & -0.0153 \end{array} \right] , \\ K_{p}^{2}\left( 1\right) & = \left[ \begin{array}{ll} -0.0079 & -0.0066\\ -0.0145 & -0.0097 \end{array} \right] ,K_{p}^{2}\left( 2\right) = \left[ \begin{array}{ll} -0.0086 & -0.0061\\ -0.0170 & -0.0102 \end{array} \right] . \end{align*} $ |
As a result, the concentrations of mRNAs and proteins estimation errors over sensor networks are depicted in Figures 1–4, where the jumping modes of GRNs and sensors are illustrated in Figure 5, respectively. It can be seen that all the sensors can display effective estimations without true GRNs system mode information and thus provides an improved nonsynchronous distributed state estimation approach. Furthermore, the triggering events of the sensors are shown in Figures 6 and 7, respectively. Compared with common time-triggered strategies, one can clearly found that the event-triggered scheme can achieve a considerable improvement of communication efficiency. Therefore, these simulation results well demonstrate our distributed state estimation approach and supports our theoretical findings.
This paper is centered on the state estimation for GRNs with HMJPs over distributed sensor networks. More precisely, a novel event-triggered sensing scheme is developed for more effective estimation. In addition, since the modes of HMJPs might be nonsynchronous with the sensors, the nonsynchronous distributed state estimation is developed with aid of observed information, which is more practical for the implementations. On the basis of model transformation, sufficient mode-dependent conditions for ensuring the modified $ H_{\infty} $ performance are established and the desired nonsynchronous state estimation gains are designed by matrix manipulation techniques. In the end, numerical simulations are performed to show the usefulness and advantages of our theoretical results. It is noted that one interesting research issue of our future work is extending the current method to the cases with more complicated mode jumping features.
This work is supported by the National Natural Science Foundation of China under Grant 62173326 and the National Key Research and Development Plan of China under Grant 2020AAA0105900, and supported by Beijing Natural Science Foundation under Grant L211023.
The authors declare there is no conflict of interest.
1. | Ahmad Tavana, Sam Ziamanesh, Hadi S. Salimi, Amir Aminzadeh Ghavifekr, Paolo Visconti, 2023, A GOA-Optimized Visible Light Communication System For Indoor High-Precision 3-D Positioning Service, 979-8-3503-1140-2, 2585, 10.1109/CoDIT58514.2023.10284069 | |
2. | 汝贵 Ru Gui, 秦岭 Qin Ling, 王凤英 Wang Fengying, 胡晓莉 Hu Xiaoli, 赵德胜 Zhao Desheng, 基于光电二极管阵列传感器的可见光定位系统, 2024, 44, 0253-2239, 2106005, 10.3788/AOS241001 | |
3. | Jiahui Wang, Xiangqing Wang, Chaochuan Jia, Cui Yang, Indoor Visible Light 3D Localization System Based on Black Wing Kite Algorithm, 2025, 13, 2169-3536, 20960, 10.1109/ACCESS.2025.3535528 | |
4. | Martin Dratnal, Lukas Danys, Radek Martinek, Bio-inspired optimization methods for visible light communication: a comprehensive review, 2025, 58, 1573-7462, 10.1007/s10462-025-11251-5 | |
5. | Xiaotong Mao, Lei Jing, Zhengrong Tong, Weihua Zhang, Peng Li, Xue Wang, Hao Wang, Tianyi Cao, Zhen Sun, Indoor visible light positioning system driven by deep neural network based on Kalman filtering and clustering optimization., 2025, 00304018, 132090, 10.1016/j.optcom.2025.132090 |