Cyber-Physical Systems are very vulnerable to sparse sensor attacks.But current protection mechanisms employ linear and deterministic models which cannot detect attacks precisely.Therefore,in this paper,we propose a n...Cyber-Physical Systems are very vulnerable to sparse sensor attacks.But current protection mechanisms employ linear and deterministic models which cannot detect attacks precisely.Therefore,in this paper,we propose a new non-linear generalized model to describe Cyber-Physical Systems.This model includes unknown multivariable discrete and continuous-time functions and different multiplicative noises to represent the evolution of physical processes and randomeffects in the physical and computationalworlds.Besides,the digitalization stage in hardware devices is represented too.Attackers and most critical sparse sensor attacks are described through a stochastic process.The reconstruction and protectionmechanisms are based on aweighted stochasticmodel.Error probability in data samples is estimated through different indicators commonly employed in non-linear dynamics(such as the Fourier transform,first-return maps,or the probability density function).A decision algorithm calculates the final reconstructed value considering the previous error probability.An experimental validation based on simulation tools and real deployments is also carried out.Both,the new technology performance and scalability are studied.Results prove that the proposed solution protects Cyber-Physical Systems against up to 92%of attacks and perturbations,with a computational delay below 2.5 s.The proposed model shows a linear complexity,as recursive or iterative structures are not employed,just algebraic and probabilistic functions.In conclusion,the new model and reconstructionmechanism can protect successfully Cyber-Physical Systems against sparse sensor attacks,even in dense or pervasive deployments and scenarios.展开更多
Stochastic modelling of hydrological time series with insufficient length and data gaps is a serious challenge since these problems significantly affect the reliability of statistical models predicting and forecasting...Stochastic modelling of hydrological time series with insufficient length and data gaps is a serious challenge since these problems significantly affect the reliability of statistical models predicting and forecasting skills.In this paper,we proposed a method for searching the seasonal autoregressive integrated moving average(SARIMA)model parameters to predict the behavior of groundwater time series affected by the issues mentioned.Based on the analysis of statistical indices,8 stations among 44 available within the Campania region(Italy)have been selected as the highest quality measurements.Different SARIMA models,with different autoregressive,moving average and differentiation orders had been used.By reviewing the criteria used to determine the consistency and goodness-of-fit of the model,it is revealed that the model with specific combination of parameters,SARIMA(0,1,3)(0,1,2)_(12),has a high R^(2) value,larger than 92%,for each of the 8 selected stations.The same model has also good performances for what concern the forecasting skills,with an average NSE of about 96%.Therefore,this study has the potential to provide a new horizon for the simulation and reconstruction of groundwater time series within the investigated area.展开更多
The dynamic and accurate forecasting of monthly streamflow processes of a river are important in the management of extreme events such as floods and drought, optimal design of water storage structures and drainage net...The dynamic and accurate forecasting of monthly streamflow processes of a river are important in the management of extreme events such as floods and drought, optimal design of water storage structures and drainage network. Many Rivers are selected in this study: White Nile, Blue Nile, Atbara River and main Nile. This paper aims to recommend the best linear stochastic model in forecasting monthly streamflow in rivers. Two commonly hydrologic models: the deseasonalized autoregressive moving average (DARMA) models and seasonal autoregressive integrated moving average (SARIMA) models are selected for modeling monthly streamflow in all Rivers in the study area. Two different types of monthly streamflow data (deseasonalized data and differenced data) were used to develop time series model using previous flow conditions as predictors. The one month ahead forecasting performances of all models for predicted period were compared. The comparison of model forecasting performance was conducted based upon graphical and numerical criteria. The result indicates that deasonalized autoregressive moving average (DARMA) models perform better than seasonal autoregressive integrated moving average (SARIMA) models for monthly streamflow in Rivers.展开更多
Short-term(up to 30 days)predictions of Earth Rotation Parameters(ERPs)such as Polar Motion(PM:PMX and PMY)play an essential role in real-time applications related to high-precision reference frame conversion.Currentl...Short-term(up to 30 days)predictions of Earth Rotation Parameters(ERPs)such as Polar Motion(PM:PMX and PMY)play an essential role in real-time applications related to high-precision reference frame conversion.Currently,least squares(LS)+auto-regressive(AR)hybrid method is one of the main techniques of PM prediction.Besides,the weighted LS+AR hybrid method performs well for PM short-term prediction.However,the corresponding covariance information of LS fitting residuals deserves further exploration in the AR model.In this study,we have derived a modified stochastic model for the LS+AR hybrid method,namely the weighted LS+weighted AR hybrid method.By using the PM data products of IERS EOP 14 C04,the numerical results indicate that for PM short-term forecasting,the proposed weighted LS+weighted AR hybrid method shows an advantage over both the LS+AR hybrid method and the weighted LS+AR hybrid method.Compared to the mean absolute errors(MAEs)of PMX/PMY sho rt-term prediction of the LS+AR hybrid method and the weighted LS+AR hybrid method,the weighted LS+weighted AR hybrid method shows average improvements of 6.61%/12.08%and 0.24%/11.65%,respectively.Besides,for the slopes of the linear regression lines fitted to the errors of each method,the growth of the prediction error of the proposed method is slower than that of the other two methods.展开更多
A patient co-infected with COVID-19 and viral hepatitis B can be atmore risk of severe complications than the one infected with a single infection.This study develops a comprehensive stochastic model to assess the epi...A patient co-infected with COVID-19 and viral hepatitis B can be atmore risk of severe complications than the one infected with a single infection.This study develops a comprehensive stochastic model to assess the epidemiological impact of vaccine booster doses on the co-dynamics of viral hepatitis B and COVID-19.The model is fitted to real COVID-19 data from Pakistan.The proposed model incorporates logistic growth and saturated incidence functions.Rigorous analyses using the tools of stochastic calculus,are performed to study appropriate conditions for the existence of unique global solutions,stationary distribution in the sense of ergodicity and disease extinction.The stochastic threshold estimated from the data fitting is given by:R_(0)^(S)=3.0651.Numerical assessments are implemented to illustrate the impact of double-dose vaccination and saturated incidence functions on the dynamics of both diseases.The effects of stochastic white noise intensities are also highlighted.展开更多
Several mathematical models have been developed to investigate the dynamics of tuberculosis(TB)and hepatitis B virus(HBV).Numerous current models for TB,HBV,and their co-dynamics fall short in capturing the important ...Several mathematical models have been developed to investigate the dynamics of tuberculosis(TB)and hepatitis B virus(HBV).Numerous current models for TB,HBV,and their co-dynamics fall short in capturing the important and practical aspect of unpredictability.It is crucial to take into account a stochastic co-infection HBV–TB epidemic model since different random elements have a substantial impact on the overall dynamics of these diseases.We provide a novel stochastic co-model for TB and HBV in this study,and we establish criteria on the uniqueness and existence of a nonnegative global solution.We also looked at the persistence of the infections as long its dynamics are governable by the proposed model.To verify the theoretical conclusions,numerical simulations are presented keeping in view the associated analytical results.The infections are found to finally die out and go extinct with certainty when L´evy intensities surpass the specified thresholds and the related stochastic thresholds fall below unity.The findings also demonstrate the impact of noise on the decline in the co-circulation of HBV and TB in a given population.Our results provide insights into effective intervention strategies,ultimately aiming to improve the management and control of TB and HBV co-infections.展开更多
The Mean First-Passage Time (MFPT) and Stochastic Resonance (SR) of a stochastic tumor-immune model withnoise perturbation are discussed in this paper. Firstly, considering environmental perturbation, Gaussian whiteno...The Mean First-Passage Time (MFPT) and Stochastic Resonance (SR) of a stochastic tumor-immune model withnoise perturbation are discussed in this paper. Firstly, considering environmental perturbation, Gaussian whitenoise and Gaussian colored noise are introduced into a tumor growth model under immune surveillance. Asfollows, the long-time evolution of the tumor characterized by the Stationary Probability Density (SPD) and MFPTis obtained in theory on the basis of the Approximated Fokker-Planck Equation (AFPE). Herein the recurrenceof the tumor from the extinction state to the tumor-present state is more concerned in this paper. A moreefficient algorithmof Back-Propagation Neural Network (BPNN) is utilized in order to testify the correction of thetheoretical SPDandMFPT.With the existence of aweak signal, the functional relationship between Signal-to-NoiseRatio (SNR), noise intensities and correlation time is also studied. Numerical results show that both multiplicativeGaussian colored noise and additive Gaussian white noise can promote the extinction of the tumors, and themultiplicative Gaussian colored noise can lead to the resonance-like peak on MFPT curves, while the increasingintensity of the additiveGaussian white noise results in theminimum of MFPT. In addition, the correlation timesare negatively correlated with MFPT. As for the SNR, we find the intensities of both the Gaussian white noise andthe Gaussian colored noise, as well as their correlation intensity can induce SR. Especially, SNR is monotonouslyincreased in the case ofGaussian white noisewith the change of the correlation time.At last, the optimal parametersin BPNN structure are analyzed for MFPT from three aspects: the penalty factors, the number of neural networklayers and the number of nodes in each layer.展开更多
Based on theWorld Health Organization(WHO),Meningitis is a severe infection of the meninges,the membranes covering the brain and spinal cord.It is a devastating disease and remains a significant public health challeng...Based on theWorld Health Organization(WHO),Meningitis is a severe infection of the meninges,the membranes covering the brain and spinal cord.It is a devastating disease and remains a significant public health challenge.This study investigates a bacterial meningitis model through deterministic and stochastic versions.Four-compartment population dynamics explain the concept,particularly the susceptible population,carrier,infected,and recovered.The model predicts the nonnegative equilibrium points and reproduction number,i.e.,the Meningitis-Free Equilibrium(MFE),and Meningitis-Existing Equilibrium(MEE).For the stochastic version of the existing deterministicmodel,the twomethodologies studied are transition probabilities and non-parametric perturbations.Also,positivity,boundedness,extinction,and disease persistence are studiedrigorouslywiththe helpofwell-known theorems.Standard and nonstandard techniques such as EulerMaruyama,stochastic Euler,stochastic Runge Kutta,and stochastic nonstandard finite difference in the sense of delay have been presented for computational analysis of the stochastic model.Unfortunately,standard methods fail to restore the biological properties of the model,so the stochastic nonstandard finite difference approximation is offered as an efficient,low-cost,and independent of time step size.In addition,the convergence,local,and global stability around the equilibria of the nonstandard computational method is studied by assuming the perturbation effect is zero.The simulations and comparison of the methods are presented to support the theoretical results and for the best visualization of results.展开更多
In the Saloum region of central-western Senegal, water needs are essentially met by tapping an underground aquifer associated with the sandy-clay formations of the Continental Terminal, in contact with both the ocean ...In the Saloum region of central-western Senegal, water needs are essentially met by tapping an underground aquifer associated with the sandy-clay formations of the Continental Terminal, in contact with both the ocean to the west and the highly saline waters of the Saloum River to the north. In this estuarine and deltaic zone with its very low relief, the hydraulic loads in the water tables are generally close to zero or even negative, creating a reversal of the natural flow and encouraging saline intrusion into this system, which makes it very vulnerable. This study concerns the implementation of a numerical model of saline intrusion to provide a better understanding of the vulnerability of the water table by analyzing the variability of the freshwater/saltwater interface. The Modflow-2005 code is used to simulate saline intrusion using the SWI2 module, coupled with the GRASS (Geographic Resources Analysis Support System) software under the Linux operating system with the steep interface approach. The probable expansion of the wedge is studied in three scenarios, taking into account its position relative to the bedrock at 1 m, 5 m and 10 m. Simulations carried out under imposed potential and river conditions, based on variations in groundwater reserves using two effective porosity values, 10−1 and 10−2, show that the water table is highly vulnerable in the northwest sector. The probable expansion of the wedge increases as the storage coefficient decreases and is more marked with river conditions in the areas surrounding the Saloum River, reaching 6 km with a probability of 1. The probability of the wedge reaching a certain degree of expansion decreases from 1 to 0.5, and then cancels out as it moves inland. The probable position of the wedge is limited to 500 m or even 1 km depending on the corner around the coast to the southwest and in the southern zone. This modelling, carried out under natural conditions, will be developed further, taking into account climatic parameters and pumping from wells and boreholes.展开更多
The research on spatial epidemic models is a topic of considerable recent interest. In another hand, the advances in computer technology have stimulated the development of stochastic models. Metapopulation models are ...The research on spatial epidemic models is a topic of considerable recent interest. In another hand, the advances in computer technology have stimulated the development of stochastic models. Metapopulation models are spatial designs that involve movements of individuals between distinct subpopulations. The purpose of the present work has been to develop stochastic models in order to study the transmission dynamics and control of infectious diseases in metapopulations. The authors studied Susceptible-Infected-Susceptible （SIS） and Susceptible-lnfected-Recovered （SIR） epidemic schemes, using the Gillespie algorithm, Computational numerical simulations were carried in order to explore the models. The results obtained show how the dynamics of transmission and the application of control measures within each subpopulation may affect all subpopulations of the system. They also show how the distribution of control measures among subpopulations affects the efficacy of these strategies. The dynamics of the stochastic models developed in the current study follow the trends observed in the classic deterministic designs. Also, the present models exhibit fluctuating behavior. This work highlights the importance of the spatial distribution of the population in spread and control of infectious diseases. In addition, it shows how chance could play an important role in these scenarios.展开更多
An optimal quota-share and excess-of-loss reinsurance and investment problem is studied for an insurer who is allowed to invest in a risk-free asset and a risky asset.Especially the price process of the risky asset is...An optimal quota-share and excess-of-loss reinsurance and investment problem is studied for an insurer who is allowed to invest in a risk-free asset and a risky asset.Especially the price process of the risky asset is governed by Heston's stochastic volatility(SV)model.With the objective of maximizing the expected index utility of the terminal wealth of the insurance company,by using the classical tools of stochastic optimal control,the explicit expressions for optimal strategies and optimal value functions are derived.An interesting conclusion is found that it is better to buy one reinsurance than two under the assumption of this paper.Moreover,some numerical simulations and sensitivity analysis are provided.展开更多
In this paper, data streams are classified into four types conforming to a standardized infrastructure of communication networks for a substation automation system （SAS） based on IEC61850 system. The data exchanged ...In this paper, data streams are classified into four types conforming to a standardized infrastructure of communication networks for a substation automation system （SAS） based on IEC61850 system. The data exchanged on the net are demonstrated to be stochastic according to investigation on the Ethemet communication principles. Four generalized stochastic Petri nets （GSPN） based models for performance analysis of communication networks of IEC61850 system are developed based on the three-level structure of SAS, different time requirements of the four data streams and different networks topology for different voltage level. The GSPN-based model associated with immediate and exponential transitions is proven to be theoretically isomorphic with Markov chain; hence we apply the mathematic methods of performance evaluation contained in Markov chain to the GSPN models proposed. The computer simulation of the model including only sample value data streams shows that it can meet performance evaluation needs of communication networks of IEC61850 system. Further researches should be focused on the pe^ormance of the other three models to explain clear how those different data streams are interrelated to and interact on each other.展开更多
Nonlinearity and randomness are both the essential attributes for the real world,and the case is the same for the models of infectious diseases,for which the deterministic models can not give a complete picture of the...Nonlinearity and randomness are both the essential attributes for the real world,and the case is the same for the models of infectious diseases,for which the deterministic models can not give a complete picture of the evolution.However,although there has been a lot of work on stochastic epidemic models,most of them focus mainly on qualitative properties,which makes us somewhat ignore the original meaning of the parameter value.In this paper we extend the classic susceptible-infectious-removed(SIR)epidemic model by adding a white noise excitation and then we utilize the large deviation theory to quantitatively study the long-term coexistence exit problem with epidemic.Finally,in order to extend the meaning of parameters in the corresponding deterministic system,we tentatively introduce two new thresholds which then prove rational.展开更多
<span style="font-family:Verdana;">In this paper we build and analyze two stochastic epidemic models with death. The model assume</span><span style="font-family:Verdana;"><span...<span style="font-family:Verdana;">In this paper we build and analyze two stochastic epidemic models with death. The model assume</span><span style="font-family:Verdana;"><span style="font-family:Verdana;"><span style="font-family:Verdana;">s</span></span></span><span style="font-family:Verdana;"><span style="font-family:Verdana;"><span style="font-family:Verdana;"> that only susceptible individuals (S) can get infected (I) and may die from this disease or a recovered individual becomes susceptible again (SIS model) or completely immune (SIR Model) for the remainder of the study period. Moreover, it is assumed there are no births, deaths, immigration or emigration during the study period;the community is said to be closed. In these infection disease models, there are two central questions: first it is the disease extinction or not and the second studies the time elapsed for such extinction, this paper will deal with this second question because the first answer corresponds to the basic reproduction number defined in the bibliography. More concretely, we study the mean-extinction of the diseases and the technique used here first builds the backward Kolmogorov differential equation and then solves it numerically using finite element method with FreeFem++. Our contribution and novelty </span></span></span><span style="font-family:Verdana;"><span style="font-family:Verdana;"><span style="font-family:Verdana;">are</span></span></span><span style="font-family:Verdana;"><span style="font-family:Verdana;"><span style="font-family:Verdana;"> the following: however the reproduction number effectively concludes the extinction or not of the disease, it does not help to know its extinction times because example with the same reproduction numbers has very different time. Moreover, the SIS model is slower, a result that is not surprising, but this difference seems to increase in the stochastic models with respect to the deterministic ones, it is reasonable to assume some uncertainly.</span></span></span>展开更多
The asymptotic stability of two species stochastic Lotka-Volterra model is explored in this paper. Firstly, the Lotka-Volterra model with random parameter is built and reduced into the equivalent deterministic system ...The asymptotic stability of two species stochastic Lotka-Volterra model is explored in this paper. Firstly, the Lotka-Volterra model with random parameter is built and reduced into the equivalent deterministic system by orthogonal polynomial approximation. Then, the linear stability theory and Routh-Hurwitz criterion for nonlinear deterministic systems are applied to the equivalent one. At last, at the aid of Lyapunov second method, we obtain that as the random intensity or statistical parameter of random variable is changed, the stability about stochastic Lotka-Volterra model is different from the deterministic system.展开更多
Process variations can reduce the accuracy in estimation of interconnect performance. This work presents a process variation based stochastic model and proposes an effective analytical method to estimate interconnect ...Process variations can reduce the accuracy in estimation of interconnect performance. This work presents a process variation based stochastic model and proposes an effective analytical method to estimate interconnect delay. The technique decouples the stochastic interconnect segments by an improved decoupling method. Combined with a polynomial chaos expression （PCE）, this paper applies the stochastic Galerkin method （SGM） to analyze the system response. A finite representation of interconnect delay is then obtained with the complex approximation method and the bisection method. Results from the analysis match well with those from SPICE. Moreover, the method shows good computational efficiency, as the running time is much less than the SPICE simulation＇s.展开更多
Weighted total least squares(WTLS)have been regarded as the standard tool for the errors-in-variables(EIV)model in which all the elements in the observation vector and the coefficient matrix are contaminated with rand...Weighted total least squares(WTLS)have been regarded as the standard tool for the errors-in-variables(EIV)model in which all the elements in the observation vector and the coefficient matrix are contaminated with random errors.However,in many geodetic applications,some elements are error-free and some random observations appear repeatedly in different positions in the augmented coefficient matrix.It is called the linear structured EIV(LSEIV)model.Two kinds of methods are proposed for the LSEIV model from functional and stochastic modifications.On the one hand,the functional part of the LSEIV model is modified into the errors-in-observations(EIO)model.On the other hand,the stochastic model is modified by applying the Moore-Penrose inverse of the cofactor matrix.The algorithms are derived through the Lagrange multipliers method and linear approximation.The estimation principles and iterative formula of the parameters are proven to be consistent.The first-order approximate variance-covariance matrix(VCM)of the parameters is also derived.A numerical example is given to compare the performances of our proposed three algorithms with the STLS approach.Afterwards,the least squares(LS),total least squares(TLS)and linear structured weighted total least squares(LSWTLS)solutions are compared and the accuracy evaluation formula is proven to be feasible and effective.Finally,the LSWTLS is applied to the field of deformation analysis,which yields a better result than the traditional LS and TLS estimations.展开更多
The issue of document management has been raised for a long time, especially with the appearance of office automation in the 1980s, which led to dematerialization and Electronic Document Management (EDM). In the same ...The issue of document management has been raised for a long time, especially with the appearance of office automation in the 1980s, which led to dematerialization and Electronic Document Management (EDM). In the same period, workflow management has experienced significant development, but has become more focused on the industry. However, it seems to us that document workflows have not had the same interest for the scientific community. But nowadays, the emergence and supremacy of the Internet in electronic exchanges are leading to a massive dematerialization of documents;which requires a conceptual reconsideration of the organizational framework for the processing of said documents in both public and private administrations. This problem seems open to us and deserves the interest of the scientific community. Indeed, EDM has mainly focused on the storage (referencing) and circulation of documents (traceability). It paid little attention to the overall behavior of the system in processing documents. The purpose of our researches is to model document processing systems. In the previous works, we proposed a general model and its specialization in the case of small documents (any document processed by a single person at a time during its processing life cycle), which represent 70% of documents processed by administrations, according to our study. In this contribution, we extend the model for processing small documents to the case where they are managed in a system comprising document classes organized in subclasses;which is the case for most administrations. We have thus observed that this model is a Markovian <i>M<sup>L×K</sup>/M<sup>L×K</sup>/</i>1 queues network. We have analyzed the constraints of this model and deduced certain characteristics and metrics. <span style="white-space:normal;"><i></i></span><i>In fine<span style="white-space:normal;"></span></i>, the ultimate objective of our work is to design a document workflow management system, integrating a component of global behavior prediction.展开更多
We consider state and parameter estimation for compartmental models having both timevarying and time-invariant parameters.In this manuscript,we first detail a general Bayesian computational framework as a continuation...We consider state and parameter estimation for compartmental models having both timevarying and time-invariant parameters.In this manuscript,we first detail a general Bayesian computational framework as a continuation of our previous work.Subsequently,this framework is specifically tailored to the susceptible-infectious-removed(SIR)model which describes a basic mechanism for the spread of infectious diseases through a system of coupled nonlinear differential equations.The SIR model consists of three states,namely,the susceptible,infectious,and removed compartments.The coupling among these states is controlled by two parameters,the infection rate and the recovery rate.The simplicity of the SIR model and similar compartmental models make them applicable to many classes of infectious diseases.However,the combined assumption of a deterministic model and time-invariance among the model parameters are two significant impediments which critically limit their use for long-term predictions.The tendency of certain model parameters to vary in time due to seasonal trends,non-pharmaceutical interventions,and other random effects necessitates a model that structurally permits the incorporation of such time-varying effects.Complementary to this,is the need for a robust mechanism for the estimation of the parameters of the resulting model from data.To this end,we consider an augmented state vector,which appends the time-varying parameters to the original system states whereby the time evolution of the time-varying parameters are driven by an artificial noise process in a standard manner.Distinguishing between time-varying and time-invariant parameters in this fashion limits the introduction of artificial dynamics into the system,and provides a robust,fully Bayesian approach for estimating the timeinvariant system parameters as well as the elements of the process noise covariance matrix.This computational framework is implemented by leveraging the robustness of the Markov chain Monte Carlo algorithm permits the estimation of time-invariant parameters while nested nonlinear filters concurrently perform the joint estimation of the system states and time-varying parameters.We demonstrate performance of the framework by first considering a series of examples using synthetic data,followed by an exposition on public health data collected in the province of Ontario.展开更多
A Markov chain-based stochastic model （MCM） is developed to simulate the movement of particles in a 2D bubbling fluidized bed （BFB）. The state spaces are determined by the discretized physical cells of the bed, an...A Markov chain-based stochastic model （MCM） is developed to simulate the movement of particles in a 2D bubbling fluidized bed （BFB）. The state spaces are determined by the discretized physical cells of the bed, and the transition probability matrix is directly calculated by the results of a discrete element method （DEM） simulation. The Markov property of the BFB is discussed by the comparison results calculated from both static and dynamic transition probability matrices. The static matrix is calculated based on the Markov chain while the dynamic matrix is calculated based on the memory property of the particle movement. Results show that the difference in the trends of particle movement between the static and dynamic matrix calculation is very small. Besides, the particle mixing curves of the MCM and DEM have the same trend and similar numerical values, and the details show the time averaged characteristic of the MCM and also expose its shortcoming in describing the instantaneous particle dynamics in the BFB.展开更多
基金supported by Comunidad de Madrid within the framework of the Multiannual Agreement with Universidad Politécnica de Madrid to encourage research by young doctors(PRINCE).
文摘Cyber-Physical Systems are very vulnerable to sparse sensor attacks.But current protection mechanisms employ linear and deterministic models which cannot detect attacks precisely.Therefore,in this paper,we propose a new non-linear generalized model to describe Cyber-Physical Systems.This model includes unknown multivariable discrete and continuous-time functions and different multiplicative noises to represent the evolution of physical processes and randomeffects in the physical and computationalworlds.Besides,the digitalization stage in hardware devices is represented too.Attackers and most critical sparse sensor attacks are described through a stochastic process.The reconstruction and protectionmechanisms are based on aweighted stochasticmodel.Error probability in data samples is estimated through different indicators commonly employed in non-linear dynamics(such as the Fourier transform,first-return maps,or the probability density function).A decision algorithm calculates the final reconstructed value considering the previous error probability.An experimental validation based on simulation tools and real deployments is also carried out.Both,the new technology performance and scalability are studied.Results prove that the proposed solution protects Cyber-Physical Systems against up to 92%of attacks and perturbations,with a computational delay below 2.5 s.The proposed model shows a linear complexity,as recursive or iterative structures are not employed,just algebraic and probabilistic functions.In conclusion,the new model and reconstructionmechanism can protect successfully Cyber-Physical Systems against sparse sensor attacks,even in dense or pervasive deployments and scenarios.
文摘Stochastic modelling of hydrological time series with insufficient length and data gaps is a serious challenge since these problems significantly affect the reliability of statistical models predicting and forecasting skills.In this paper,we proposed a method for searching the seasonal autoregressive integrated moving average(SARIMA)model parameters to predict the behavior of groundwater time series affected by the issues mentioned.Based on the analysis of statistical indices,8 stations among 44 available within the Campania region(Italy)have been selected as the highest quality measurements.Different SARIMA models,with different autoregressive,moving average and differentiation orders had been used.By reviewing the criteria used to determine the consistency and goodness-of-fit of the model,it is revealed that the model with specific combination of parameters,SARIMA(0,1,3)(0,1,2)_(12),has a high R^(2) value,larger than 92%,for each of the 8 selected stations.The same model has also good performances for what concern the forecasting skills,with an average NSE of about 96%.Therefore,this study has the potential to provide a new horizon for the simulation and reconstruction of groundwater time series within the investigated area.
文摘The dynamic and accurate forecasting of monthly streamflow processes of a river are important in the management of extreme events such as floods and drought, optimal design of water storage structures and drainage network. Many Rivers are selected in this study: White Nile, Blue Nile, Atbara River and main Nile. This paper aims to recommend the best linear stochastic model in forecasting monthly streamflow in rivers. Two commonly hydrologic models: the deseasonalized autoregressive moving average (DARMA) models and seasonal autoregressive integrated moving average (SARIMA) models are selected for modeling monthly streamflow in all Rivers in the study area. Two different types of monthly streamflow data (deseasonalized data and differenced data) were used to develop time series model using previous flow conditions as predictors. The one month ahead forecasting performances of all models for predicted period were compared. The comparison of model forecasting performance was conducted based upon graphical and numerical criteria. The result indicates that deasonalized autoregressive moving average (DARMA) models perform better than seasonal autoregressive integrated moving average (SARIMA) models for monthly streamflow in Rivers.
基金supported by National Natural Science Foundation of China,China(No.42004016)HuBei Natural Science Fund,China(No.2020CFB329)+1 种基金HuNan Natural Science Fund,China(No.2023JJ60559,2023JJ60560)the State Key Laboratory of Geodesy and Earth’s Dynamics self-deployment project,China(No.S21L6101)。
文摘Short-term(up to 30 days)predictions of Earth Rotation Parameters(ERPs)such as Polar Motion(PM:PMX and PMY)play an essential role in real-time applications related to high-precision reference frame conversion.Currently,least squares(LS)+auto-regressive(AR)hybrid method is one of the main techniques of PM prediction.Besides,the weighted LS+AR hybrid method performs well for PM short-term prediction.However,the corresponding covariance information of LS fitting residuals deserves further exploration in the AR model.In this study,we have derived a modified stochastic model for the LS+AR hybrid method,namely the weighted LS+weighted AR hybrid method.By using the PM data products of IERS EOP 14 C04,the numerical results indicate that for PM short-term forecasting,the proposed weighted LS+weighted AR hybrid method shows an advantage over both the LS+AR hybrid method and the weighted LS+AR hybrid method.Compared to the mean absolute errors(MAEs)of PMX/PMY sho rt-term prediction of the LS+AR hybrid method and the weighted LS+AR hybrid method,the weighted LS+weighted AR hybrid method shows average improvements of 6.61%/12.08%and 0.24%/11.65%,respectively.Besides,for the slopes of the linear regression lines fitted to the errors of each method,the growth of the prediction error of the proposed method is slower than that of the other two methods.
文摘A patient co-infected with COVID-19 and viral hepatitis B can be atmore risk of severe complications than the one infected with a single infection.This study develops a comprehensive stochastic model to assess the epidemiological impact of vaccine booster doses on the co-dynamics of viral hepatitis B and COVID-19.The model is fitted to real COVID-19 data from Pakistan.The proposed model incorporates logistic growth and saturated incidence functions.Rigorous analyses using the tools of stochastic calculus,are performed to study appropriate conditions for the existence of unique global solutions,stationary distribution in the sense of ergodicity and disease extinction.The stochastic threshold estimated from the data fitting is given by:R_(0)^(S)=3.0651.Numerical assessments are implemented to illustrate the impact of double-dose vaccination and saturated incidence functions on the dynamics of both diseases.The effects of stochastic white noise intensities are also highlighted.
文摘Several mathematical models have been developed to investigate the dynamics of tuberculosis(TB)and hepatitis B virus(HBV).Numerous current models for TB,HBV,and their co-dynamics fall short in capturing the important and practical aspect of unpredictability.It is crucial to take into account a stochastic co-infection HBV–TB epidemic model since different random elements have a substantial impact on the overall dynamics of these diseases.We provide a novel stochastic co-model for TB and HBV in this study,and we establish criteria on the uniqueness and existence of a nonnegative global solution.We also looked at the persistence of the infections as long its dynamics are governable by the proposed model.To verify the theoretical conclusions,numerical simulations are presented keeping in view the associated analytical results.The infections are found to finally die out and go extinct with certainty when L´evy intensities surpass the specified thresholds and the related stochastic thresholds fall below unity.The findings also demonstrate the impact of noise on the decline in the co-circulation of HBV and TB in a given population.Our results provide insights into effective intervention strategies,ultimately aiming to improve the management and control of TB and HBV co-infections.
基金National Natural Science Foundation of China(Nos.12272283,12172266).
文摘The Mean First-Passage Time (MFPT) and Stochastic Resonance (SR) of a stochastic tumor-immune model withnoise perturbation are discussed in this paper. Firstly, considering environmental perturbation, Gaussian whitenoise and Gaussian colored noise are introduced into a tumor growth model under immune surveillance. Asfollows, the long-time evolution of the tumor characterized by the Stationary Probability Density (SPD) and MFPTis obtained in theory on the basis of the Approximated Fokker-Planck Equation (AFPE). Herein the recurrenceof the tumor from the extinction state to the tumor-present state is more concerned in this paper. A moreefficient algorithmof Back-Propagation Neural Network (BPNN) is utilized in order to testify the correction of thetheoretical SPDandMFPT.With the existence of aweak signal, the functional relationship between Signal-to-NoiseRatio (SNR), noise intensities and correlation time is also studied. Numerical results show that both multiplicativeGaussian colored noise and additive Gaussian white noise can promote the extinction of the tumors, and themultiplicative Gaussian colored noise can lead to the resonance-like peak on MFPT curves, while the increasingintensity of the additiveGaussian white noise results in theminimum of MFPT. In addition, the correlation timesare negatively correlated with MFPT. As for the SNR, we find the intensities of both the Gaussian white noise andthe Gaussian colored noise, as well as their correlation intensity can induce SR. Especially, SNR is monotonouslyincreased in the case ofGaussian white noisewith the change of the correlation time.At last, the optimal parametersin BPNN structure are analyzed for MFPT from three aspects: the penalty factors, the number of neural networklayers and the number of nodes in each layer.
基金Deanship of Research and Graduate Studies at King Khalid University for funding this work through large Research Project under Grant Number RGP2/302/45supported by the Deanship of Scientific Research,Vice Presidency forGraduate Studies and Scientific Research,King Faisal University,Saudi Arabia(Grant Number A426).
文摘Based on theWorld Health Organization(WHO),Meningitis is a severe infection of the meninges,the membranes covering the brain and spinal cord.It is a devastating disease and remains a significant public health challenge.This study investigates a bacterial meningitis model through deterministic and stochastic versions.Four-compartment population dynamics explain the concept,particularly the susceptible population,carrier,infected,and recovered.The model predicts the nonnegative equilibrium points and reproduction number,i.e.,the Meningitis-Free Equilibrium(MFE),and Meningitis-Existing Equilibrium(MEE).For the stochastic version of the existing deterministicmodel,the twomethodologies studied are transition probabilities and non-parametric perturbations.Also,positivity,boundedness,extinction,and disease persistence are studiedrigorouslywiththe helpofwell-known theorems.Standard and nonstandard techniques such as EulerMaruyama,stochastic Euler,stochastic Runge Kutta,and stochastic nonstandard finite difference in the sense of delay have been presented for computational analysis of the stochastic model.Unfortunately,standard methods fail to restore the biological properties of the model,so the stochastic nonstandard finite difference approximation is offered as an efficient,low-cost,and independent of time step size.In addition,the convergence,local,and global stability around the equilibria of the nonstandard computational method is studied by assuming the perturbation effect is zero.The simulations and comparison of the methods are presented to support the theoretical results and for the best visualization of results.
文摘In the Saloum region of central-western Senegal, water needs are essentially met by tapping an underground aquifer associated with the sandy-clay formations of the Continental Terminal, in contact with both the ocean to the west and the highly saline waters of the Saloum River to the north. In this estuarine and deltaic zone with its very low relief, the hydraulic loads in the water tables are generally close to zero or even negative, creating a reversal of the natural flow and encouraging saline intrusion into this system, which makes it very vulnerable. This study concerns the implementation of a numerical model of saline intrusion to provide a better understanding of the vulnerability of the water table by analyzing the variability of the freshwater/saltwater interface. The Modflow-2005 code is used to simulate saline intrusion using the SWI2 module, coupled with the GRASS (Geographic Resources Analysis Support System) software under the Linux operating system with the steep interface approach. The probable expansion of the wedge is studied in three scenarios, taking into account its position relative to the bedrock at 1 m, 5 m and 10 m. Simulations carried out under imposed potential and river conditions, based on variations in groundwater reserves using two effective porosity values, 10−1 and 10−2, show that the water table is highly vulnerable in the northwest sector. The probable expansion of the wedge increases as the storage coefficient decreases and is more marked with river conditions in the areas surrounding the Saloum River, reaching 6 km with a probability of 1. The probability of the wedge reaching a certain degree of expansion decreases from 1 to 0.5, and then cancels out as it moves inland. The probable position of the wedge is limited to 500 m or even 1 km depending on the corner around the coast to the southwest and in the southern zone. This modelling, carried out under natural conditions, will be developed further, taking into account climatic parameters and pumping from wells and boreholes.
文摘The research on spatial epidemic models is a topic of considerable recent interest. In another hand, the advances in computer technology have stimulated the development of stochastic models. Metapopulation models are spatial designs that involve movements of individuals between distinct subpopulations. The purpose of the present work has been to develop stochastic models in order to study the transmission dynamics and control of infectious diseases in metapopulations. The authors studied Susceptible-Infected-Susceptible （SIS） and Susceptible-lnfected-Recovered （SIR） epidemic schemes, using the Gillespie algorithm, Computational numerical simulations were carried in order to explore the models. The results obtained show how the dynamics of transmission and the application of control measures within each subpopulation may affect all subpopulations of the system. They also show how the distribution of control measures among subpopulations affects the efficacy of these strategies. The dynamics of the stochastic models developed in the current study follow the trends observed in the classic deterministic designs. Also, the present models exhibit fluctuating behavior. This work highlights the importance of the spatial distribution of the population in spread and control of infectious diseases. In addition, it shows how chance could play an important role in these scenarios.
基金National Natural Science Foundation of China(No.62073071)Fundamental Research Funds for the Central Universities and Graduate Student Innovation Fund of Donghua University,China(No.CUSF-DH-D-2021045)。
文摘An optimal quota-share and excess-of-loss reinsurance and investment problem is studied for an insurer who is allowed to invest in a risk-free asset and a risky asset.Especially the price process of the risky asset is governed by Heston's stochastic volatility(SV)model.With the objective of maximizing the expected index utility of the terminal wealth of the insurance company,by using the classical tools of stochastic optimal control,the explicit expressions for optimal strategies and optimal value functions are derived.An interesting conclusion is found that it is better to buy one reinsurance than two under the assumption of this paper.Moreover,some numerical simulations and sensitivity analysis are provided.
文摘In this paper, data streams are classified into four types conforming to a standardized infrastructure of communication networks for a substation automation system （SAS） based on IEC61850 system. The data exchanged on the net are demonstrated to be stochastic according to investigation on the Ethemet communication principles. Four generalized stochastic Petri nets （GSPN） based models for performance analysis of communication networks of IEC61850 system are developed based on the three-level structure of SAS, different time requirements of the four data streams and different networks topology for different voltage level. The GSPN-based model associated with immediate and exponential transitions is proven to be theoretically isomorphic with Markov chain; hence we apply the mathematic methods of performance evaluation contained in Markov chain to the GSPN models proposed. The computer simulation of the model including only sample value data streams shows that it can meet performance evaluation needs of communication networks of IEC61850 system. Further researches should be focused on the pe^ormance of the other three models to explain clear how those different data streams are interrelated to and interact on each other.
基金supported by the National Natural Science Foundation of China(No.12172167)。
文摘Nonlinearity and randomness are both the essential attributes for the real world,and the case is the same for the models of infectious diseases,for which the deterministic models can not give a complete picture of the evolution.However,although there has been a lot of work on stochastic epidemic models,most of them focus mainly on qualitative properties,which makes us somewhat ignore the original meaning of the parameter value.In this paper we extend the classic susceptible-infectious-removed(SIR)epidemic model by adding a white noise excitation and then we utilize the large deviation theory to quantitatively study the long-term coexistence exit problem with epidemic.Finally,in order to extend the meaning of parameters in the corresponding deterministic system,we tentatively introduce two new thresholds which then prove rational.
文摘<span style="font-family:Verdana;">In this paper we build and analyze two stochastic epidemic models with death. The model assume</span><span style="font-family:Verdana;"><span style="font-family:Verdana;"><span style="font-family:Verdana;">s</span></span></span><span style="font-family:Verdana;"><span style="font-family:Verdana;"><span style="font-family:Verdana;"> that only susceptible individuals (S) can get infected (I) and may die from this disease or a recovered individual becomes susceptible again (SIS model) or completely immune (SIR Model) for the remainder of the study period. Moreover, it is assumed there are no births, deaths, immigration or emigration during the study period;the community is said to be closed. In these infection disease models, there are two central questions: first it is the disease extinction or not and the second studies the time elapsed for such extinction, this paper will deal with this second question because the first answer corresponds to the basic reproduction number defined in the bibliography. More concretely, we study the mean-extinction of the diseases and the technique used here first builds the backward Kolmogorov differential equation and then solves it numerically using finite element method with FreeFem++. Our contribution and novelty </span></span></span><span style="font-family:Verdana;"><span style="font-family:Verdana;"><span style="font-family:Verdana;">are</span></span></span><span style="font-family:Verdana;"><span style="font-family:Verdana;"><span style="font-family:Verdana;"> the following: however the reproduction number effectively concludes the extinction or not of the disease, it does not help to know its extinction times because example with the same reproduction numbers has very different time. Moreover, the SIS model is slower, a result that is not surprising, but this difference seems to increase in the stochastic models with respect to the deterministic ones, it is reasonable to assume some uncertainly.</span></span></span>
文摘The asymptotic stability of two species stochastic Lotka-Volterra model is explored in this paper. Firstly, the Lotka-Volterra model with random parameter is built and reduced into the equivalent deterministic system by orthogonal polynomial approximation. Then, the linear stability theory and Routh-Hurwitz criterion for nonlinear deterministic systems are applied to the equivalent one. At last, at the aid of Lyapunov second method, we obtain that as the random intensity or statistical parameter of random variable is changed, the stability about stochastic Lotka-Volterra model is different from the deterministic system.
文摘Process variations can reduce the accuracy in estimation of interconnect performance. This work presents a process variation based stochastic model and proposes an effective analytical method to estimate interconnect delay. The technique decouples the stochastic interconnect segments by an improved decoupling method. Combined with a polynomial chaos expression （PCE）, this paper applies the stochastic Galerkin method （SGM） to analyze the system response. A finite representation of interconnect delay is then obtained with the complex approximation method and the bisection method. Results from the analysis match well with those from SPICE. Moreover, the method shows good computational efficiency, as the running time is much less than the SPICE simulation＇s.
基金the financial support of the National Natural Science Foundation of China(Grant No.42074016,42104025,42274057and 41704007)Hunan Provincial Natural Science Foundation of China(Grant No.2021JJ30244)Scientific Research Fund of Hunan Provincial Education Department(Grant No.22B0496)。
文摘Weighted total least squares(WTLS)have been regarded as the standard tool for the errors-in-variables(EIV)model in which all the elements in the observation vector and the coefficient matrix are contaminated with random errors.However,in many geodetic applications,some elements are error-free and some random observations appear repeatedly in different positions in the augmented coefficient matrix.It is called the linear structured EIV(LSEIV)model.Two kinds of methods are proposed for the LSEIV model from functional and stochastic modifications.On the one hand,the functional part of the LSEIV model is modified into the errors-in-observations(EIO)model.On the other hand,the stochastic model is modified by applying the Moore-Penrose inverse of the cofactor matrix.The algorithms are derived through the Lagrange multipliers method and linear approximation.The estimation principles and iterative formula of the parameters are proven to be consistent.The first-order approximate variance-covariance matrix(VCM)of the parameters is also derived.A numerical example is given to compare the performances of our proposed three algorithms with the STLS approach.Afterwards,the least squares(LS),total least squares(TLS)and linear structured weighted total least squares(LSWTLS)solutions are compared and the accuracy evaluation formula is proven to be feasible and effective.Finally,the LSWTLS is applied to the field of deformation analysis,which yields a better result than the traditional LS and TLS estimations.
文摘The issue of document management has been raised for a long time, especially with the appearance of office automation in the 1980s, which led to dematerialization and Electronic Document Management (EDM). In the same period, workflow management has experienced significant development, but has become more focused on the industry. However, it seems to us that document workflows have not had the same interest for the scientific community. But nowadays, the emergence and supremacy of the Internet in electronic exchanges are leading to a massive dematerialization of documents;which requires a conceptual reconsideration of the organizational framework for the processing of said documents in both public and private administrations. This problem seems open to us and deserves the interest of the scientific community. Indeed, EDM has mainly focused on the storage (referencing) and circulation of documents (traceability). It paid little attention to the overall behavior of the system in processing documents. The purpose of our researches is to model document processing systems. In the previous works, we proposed a general model and its specialization in the case of small documents (any document processed by a single person at a time during its processing life cycle), which represent 70% of documents processed by administrations, according to our study. In this contribution, we extend the model for processing small documents to the case where they are managed in a system comprising document classes organized in subclasses;which is the case for most administrations. We have thus observed that this model is a Markovian <i>M<sup>L×K</sup>/M<sup>L×K</sup>/</i>1 queues network. We have analyzed the constraints of this model and deduced certain characteristics and metrics. <span style="white-space:normal;"><i></i></span><i>In fine<span style="white-space:normal;"></span></i>, the ultimate objective of our work is to design a document workflow management system, integrating a component of global behavior prediction.
基金the funding from the New Frontiers in Research Fund(NFRF)2022 Special Call e Research for Postpandemic Recovery(Grant no:NFRFR-2022-00395).
文摘We consider state and parameter estimation for compartmental models having both timevarying and time-invariant parameters.In this manuscript,we first detail a general Bayesian computational framework as a continuation of our previous work.Subsequently,this framework is specifically tailored to the susceptible-infectious-removed(SIR)model which describes a basic mechanism for the spread of infectious diseases through a system of coupled nonlinear differential equations.The SIR model consists of three states,namely,the susceptible,infectious,and removed compartments.The coupling among these states is controlled by two parameters,the infection rate and the recovery rate.The simplicity of the SIR model and similar compartmental models make them applicable to many classes of infectious diseases.However,the combined assumption of a deterministic model and time-invariance among the model parameters are two significant impediments which critically limit their use for long-term predictions.The tendency of certain model parameters to vary in time due to seasonal trends,non-pharmaceutical interventions,and other random effects necessitates a model that structurally permits the incorporation of such time-varying effects.Complementary to this,is the need for a robust mechanism for the estimation of the parameters of the resulting model from data.To this end,we consider an augmented state vector,which appends the time-varying parameters to the original system states whereby the time evolution of the time-varying parameters are driven by an artificial noise process in a standard manner.Distinguishing between time-varying and time-invariant parameters in this fashion limits the introduction of artificial dynamics into the system,and provides a robust,fully Bayesian approach for estimating the timeinvariant system parameters as well as the elements of the process noise covariance matrix.This computational framework is implemented by leveraging the robustness of the Markov chain Monte Carlo algorithm permits the estimation of time-invariant parameters while nested nonlinear filters concurrently perform the joint estimation of the system states and time-varying parameters.We demonstrate performance of the framework by first considering a series of examples using synthetic data,followed by an exposition on public health data collected in the province of Ontario.
基金The National Science Foundation of China(No.51276036,51306035)the Fundamental Research Funds for the Central Universities(No.KYLX_0114)
文摘A Markov chain-based stochastic model （MCM） is developed to simulate the movement of particles in a 2D bubbling fluidized bed （BFB）. The state spaces are determined by the discretized physical cells of the bed, and the transition probability matrix is directly calculated by the results of a discrete element method （DEM） simulation. The Markov property of the BFB is discussed by the comparison results calculated from both static and dynamic transition probability matrices. The static matrix is calculated based on the Markov chain while the dynamic matrix is calculated based on the memory property of the particle movement. Results show that the difference in the trends of particle movement between the static and dynamic matrix calculation is very small. Besides, the particle mixing curves of the MCM and DEM have the same trend and similar numerical values, and the details show the time averaged characteristic of the MCM and also expose its shortcoming in describing the instantaneous particle dynamics in the BFB.
