The objective of this paper is to evaluate the reliability of a system in its different states (absence of failures, partial failure and total failure) and to propose actions to improve this reliability by an approach...The objective of this paper is to evaluate the reliability of a system in its different states (absence of failures, partial failure and total failure) and to propose actions to improve this reliability by an approach based on Monte Carlo simulation. It consists of a probabilistic evaluation based on Markov Chains. In order to achieve this goal, the functionalities of Markov Chains and Monte Carlo simulation steps are deployed. The application is made on a production system. .展开更多
A single set of vertically aligned cracks embedded in a purely isotropic background may be con- sidered as a long-wavelength effective transversely iso- tropy (HTI) medium with a horizontal symmetry axis. The crack-...A single set of vertically aligned cracks embedded in a purely isotropic background may be con- sidered as a long-wavelength effective transversely iso- tropy (HTI) medium with a horizontal symmetry axis. The crack-induced HTI anisotropy can be characterized by the weakly anisotropic parameters introduced by Thomsen. The seismic scattering theory can be utilized for the inversion for the anisotropic parameters in weakly aniso- tropic and heterogeneous HTI media. Based on the seismic scattering theory, we first derived the linearized PP- and PS-wave reflection coefficients in terms of P- and S-wave impedances, density as well as three anisotropic parameters in HTI media. Then, we proposed a novel Bayesian Mar- kov chain Monte Carlo inversion method of PP- and PS- wave for six elastic and anisotropic parameters directly. Tests on synthetic azimuthal seismic data contaminated by random errors demonstrated that this method appears more accurate, anti-noise and stable owing to the usage of the constrained PS-wave compared with the standards inver- sion scheme taking only the PP-wave into account.展开更多
Covering a wide range of bulk densities, density profiles for hard-sphere chain fluids (HSCFs) with chain length of 3,4,8,20,32 and 64 confined between two surfaces were obtained by Monte Carlo simulations using exten...Covering a wide range of bulk densities, density profiles for hard-sphere chain fluids (HSCFs) with chain length of 3,4,8,20,32 and 64 confined between two surfaces were obtained by Monte Carlo simulations using extended continuum configurational-bias (ECCB) method. It is shown that the enrichment of beads near surfaces is happened at high densities due to the bulk packing effect, on the contrary, the depletion is revealed at low densities owing to the configurational entropic contribution. Comparisons with those calculated by density functional theory presented by Cai et al. indicate that the agreement between simulations and predictions is good. Compressibility factors of bulk HSCFs calculated using volume fractions at surfaces were also used to test the reliability of various equations of state of HSCFs by different authors.展开更多
This paper addresses the issues of channel estimation in a Multiple-Input/Multiple-Output (MIMO) system. Markov Chain Monte Carlo (MCMC) method is employed to jointly estimate the Channel State Information (CSI) and t...This paper addresses the issues of channel estimation in a Multiple-Input/Multiple-Output (MIMO) system. Markov Chain Monte Carlo (MCMC) method is employed to jointly estimate the Channel State Information (CSI) and the transmitted signals. The deduced algorithms can work well under circumstances of low Signal-to-Noise Ratio (SNR). Simulation results are presented to demonstrate their effectiveness.展开更多
目的 针对妇幼卫生纵向数据的任意缺失模式,采用多重填补方法进行填补,探求最佳填补结果,以便对数据作进一步分析与研究。方法 运用SAS9.0 ,采用多重填补方法Markov China Monte Carlo(MCMC)模型对缺失数据进行多次填补并综合分析。...目的 针对妇幼卫生纵向数据的任意缺失模式,采用多重填补方法进行填补,探求最佳填补结果,以便对数据作进一步分析与研究。方法 运用SAS9.0 ,采用多重填补方法Markov China Monte Carlo(MCMC)模型对缺失数据进行多次填补并综合分析。结果 填补5次所得结果最优。结论 多重填补方法可以处理有缺失数据资料中的许多普遍问题,可提高统计效率,尤其是MCMC模型在处理复杂的缺失数据上,优势明显。展开更多
With increasing complexity of today’s electromagnetic problems, the need and opportunity to reduce domain sizes, memory requirement, computational time and possibility of errors abound for symmetric domains. With sev...With increasing complexity of today’s electromagnetic problems, the need and opportunity to reduce domain sizes, memory requirement, computational time and possibility of errors abound for symmetric domains. With several competing computational methods in recent times, methods with little or no iterations are generally preferred as they tend to consume less computer memory resources and time. This paper presents the application of simple and efficient Markov Chain Monte Carlo (MCMC) method to the Laplace’s equation in axisymmetric homogeneous domains. Two cases of axisymmetric homogeneous problems are considered. Simulation results for analytical, finite difference and MCMC solutions are reported. The results obtained from the MCMC method agree with analytical and finite difference solutions. However, the MCMC method has the advantage that its implementation is simple and fast.展开更多
目的探讨基于Markov Chain Monte Carlo(MCMC)模型的多重估算法在处理医院调查资料缺失数据中的应用。方法运用SAS9.2编写程序,在分析数据的分布类型和缺失机制的基础上,采用MCMC法对缺失数据进行多次填补和联合统计推断,分析多重估算...目的探讨基于Markov Chain Monte Carlo(MCMC)模型的多重估算法在处理医院调查资料缺失数据中的应用。方法运用SAS9.2编写程序,在分析数据的分布类型和缺失机制的基础上,采用MCMC法对缺失数据进行多次填补和联合统计推断,分析多重估算法的优势。结果数据服从多元正态分布与随机缺失,采用MCMC法填补10次所得的结果最佳。结论多重估算既可反映缺失数据的不确定性,又可充分利用现有资料的信息、提高统计效率、对模型的估计结果更加可信,是处理缺失数据的有效方法。展开更多
The magnetic properties of a mixed spin-2 and spin-1/2 ferromagnetic diamond chain are studied by effective-field theory and Monte Carlo(MC) simulation based on the Ising model.The temperature dependences of magneti...The magnetic properties of a mixed spin-2 and spin-1/2 ferromagnetic diamond chain are studied by effective-field theory and Monte Carlo(MC) simulation based on the Ising model.The temperature dependences of magnetization,magnetic susceptibility,internal energy,and specific heat are studied,respectively.The exchange interaction dependences of magnetization and the critical temperature are calculated by MC simulation.The changes of magnetization depending on the field increasing and then the field decreasing under steady-static conditions are also given.展开更多
A program of Monte Carlo simulation of binary copolymerization for E-SBR (emulsion polymn. SB rubber) was made according to the terminal model. The simulation results obtained by this program were in good agreement wi...A program of Monte Carlo simulation of binary copolymerization for E-SBR (emulsion polymn. SB rubber) was made according to the terminal model. The simulation results obtained by this program were in good agreement with those experimental ones. A detail microstructure information of E-SBR molecular chain has been provided.展开更多
This paper proposes a new technique based on inverse Markov chain Monte Carlo algorithm for finding the smallest generalized eigenpair of the large scale matrices. Some numerical examples show that the proposed method...This paper proposes a new technique based on inverse Markov chain Monte Carlo algorithm for finding the smallest generalized eigenpair of the large scale matrices. Some numerical examples show that the proposed method is efficient.展开更多
A part of a long DNA chain was driven into a confined environment by an electric field, while the rest remains in the higher-entropy region. Upon removal of the field, the chain recoils to the higher-entropy region sp...A part of a long DNA chain was driven into a confined environment by an electric field, while the rest remains in the higher-entropy region. Upon removal of the field, the chain recoils to the higher-entropy region spontaneously. This dynamical process was investigated by Monte Carlo simulations. The simulation reproduces the experimentally-observed phenomenon that the recoil of the DNA chain is initially slow and gradually increases in speed due to the presence of the confinement-entropic force. The results show that with increasing the dimension or decreasing the spacing of the nanopillars the recoil velocity of the DNA chain will increase. Further analysis suggests that the characteristic entropy per monomer in the confinement is proportional to the area fraction of the free part in the confinement.展开更多
The configurational properties of tail-like polymer chains with one end attached to a flat surface are studied by using dynamic Monte Carlo technique. We find that the probability distribution of the free end in z dir...The configurational properties of tail-like polymer chains with one end attached to a flat surface are studied by using dynamic Monte Carlo technique. We find that the probability distribution of the free end in z direction P(Rz) and the density profile p(z) can be scaled approximately by a factor β to be a length independent function for both random walking (RW) and self-avoiding walking (SAW) tail-like chains, where the factor β is related to the mean square end-to-end distance 〈RE〉. The scaled P(Rz) of the SAW chain roughly overlaps that of the RW chain, but the scaled p(z) of the SAW chain locates at smaller βz than that of the RW chain.展开更多
We introduce the potential-decomposition strategy (PDS), which can be used in Markov chain Monte Carlo sampling algorithms. PDS can be designed to make particles move in a modified potential that favors diffusion in...We introduce the potential-decomposition strategy (PDS), which can be used in Markov chain Monte Carlo sampling algorithms. PDS can be designed to make particles move in a modified potential that favors diffusion in phase space, then, by rejecting some trial samples, the target distributions can be sampled in an unbiased manner. Furthermore, if the accepted trial samples are insumcient, they can be recycled as initial states to form more unbiased samples. This strategy can greatly improve efficiency when the original potential has multiple metastable states separated by large barriers. We apply PDS to the 2d Ising model and a double-well potential model with a large barrier, demonstrating in these two representative examples that convergence is accelerated by orders of magnitude.展开更多
Recently, a new soft-in soft-out detection algorithm based on the Markov Chain Monte Carlo (MCMC) simulation technique for Multiple-Input Multiple-Output (MIMO) systems is proposed, which is shown to perform significa...Recently, a new soft-in soft-out detection algorithm based on the Markov Chain Monte Carlo (MCMC) simulation technique for Multiple-Input Multiple-Output (MIMO) systems is proposed, which is shown to perform significantly better than their sphere decoding counterparts with relatively low complexity. However, the MCMC simulator is likely to get trapped in a fixed state when the channel SNR is high, thus lots of repetitive samples are observed and the accuracy of A Posteriori Probability (APP) estimation deteriorates. To solve this problem, an improved version of MCMC simulator, named forced-dispersed MCMC algorithm is proposed. Based on the a posteriori variance of each bit, the Gibbs sampler is monitored. Once the trapped state is detected, the sample is dispersed intentionally according to the a posteriori variance. Extensive simulation shows that, compared with the existing solution, the proposed algorithm enables the markov chain to travel more states, which ensures a near-optimal performance.展开更多
目的:为急性冠脉综合征(ACS)患者抗血小板治疗提供决策依据。方法:收集PLATO和TREAT研究中的相关数据建立Markov模型,计算ACS患者应用替格瑞洛或氯吡格雷的总出血风险、主要出血风险、次要出血风险及致命出血风险等,依据既往文献收集、...目的:为急性冠脉综合征(ACS)患者抗血小板治疗提供决策依据。方法:收集PLATO和TREAT研究中的相关数据建立Markov模型,计算ACS患者应用替格瑞洛或氯吡格雷的总出血风险、主要出血风险、次要出血风险及致命出血风险等,依据既往文献收集、计算各状态间的转移概率、各状态的成本及效用值,运用TreeAge Pro 2011软件计算不同方案的医疗费用、质量调整生命年(QALYs)以及增量-成本效果比(ICER),并对模型进行单因素敏感性分析及概率敏感性分析。结果:替格瑞洛组平均总成本是66449.38元,获得7.34 QALYs;氯吡格雷组平均总成本是53846.03元,获得6.68 QALYs。替格瑞洛组与氯吡格雷组相比的ICER为19095.98元/QALYs,即每多获得1个QALYs,替格瑞洛组患者需多花19095.98元,低于意愿支付阈值(64644元)。敏感性分析结果与上述结果一致。结论:相较于氯吡格雷,替格瑞洛对ACS患者尤其是发生出血不良事件的ACS患者具有更小的经济负担。展开更多
This paper first applies the sequential cluster method to set up the classification standard of infectious disease incidence state based on the fact that there are many uncertainty characteristics in the incidence cou...This paper first applies the sequential cluster method to set up the classification standard of infectious disease incidence state based on the fact that there are many uncertainty characteristics in the incidence course.Then the paper presents a weighted Markov chain,a method which is used to predict the future incidence state.This method assumes the standardized self-coefficients as weights based on the special characteristics of infectious disease incidence being a dependent stochastic variable.It also analyzes the characteristics of infectious diseases incidence via the Markov chain Monte Carlo method to make the long-term benefit of decision optimal.Our method is successfully validated using existing incidents data of infectious diseases in Jiangsu Province.In summation,this paper proposes ways to improve the accuracy of the weighted Markov chain,specifically in the field of infection epidemiology.展开更多
文摘The objective of this paper is to evaluate the reliability of a system in its different states (absence of failures, partial failure and total failure) and to propose actions to improve this reliability by an approach based on Monte Carlo simulation. It consists of a probabilistic evaluation based on Markov Chains. In order to achieve this goal, the functionalities of Markov Chains and Monte Carlo simulation steps are deployed. The application is made on a production system. .
基金sponsorship of the National Natural Science Foundation of China (No.41674130)the National Basic Research Program of China (973 Program,Nos.2013CB228604,2014CB239201)+1 种基金the National Oil and Gas Major Projects of China (Nos.2016ZX05027004-001,2016ZX05002005-009)the Fundamental Research Funds for the Central Universities (15CX08002A) for their funding in this research
文摘A single set of vertically aligned cracks embedded in a purely isotropic background may be con- sidered as a long-wavelength effective transversely iso- tropy (HTI) medium with a horizontal symmetry axis. The crack-induced HTI anisotropy can be characterized by the weakly anisotropic parameters introduced by Thomsen. The seismic scattering theory can be utilized for the inversion for the anisotropic parameters in weakly aniso- tropic and heterogeneous HTI media. Based on the seismic scattering theory, we first derived the linearized PP- and PS-wave reflection coefficients in terms of P- and S-wave impedances, density as well as three anisotropic parameters in HTI media. Then, we proposed a novel Bayesian Mar- kov chain Monte Carlo inversion method of PP- and PS- wave for six elastic and anisotropic parameters directly. Tests on synthetic azimuthal seismic data contaminated by random errors demonstrated that this method appears more accurate, anti-noise and stable owing to the usage of the constrained PS-wave compared with the standards inver- sion scheme taking only the PP-wave into account.
基金Supported by the National Science Foundation of China (No. 29736170, No. 20025618) and the Doctoral Research Foundation by Ministry of Education of China (No. 1999025103). Additional support provided by the Visiting Researcher Foundation of University La
文摘Covering a wide range of bulk densities, density profiles for hard-sphere chain fluids (HSCFs) with chain length of 3,4,8,20,32 and 64 confined between two surfaces were obtained by Monte Carlo simulations using extended continuum configurational-bias (ECCB) method. It is shown that the enrichment of beads near surfaces is happened at high densities due to the bulk packing effect, on the contrary, the depletion is revealed at low densities owing to the configurational entropic contribution. Comparisons with those calculated by density functional theory presented by Cai et al. indicate that the agreement between simulations and predictions is good. Compressibility factors of bulk HSCFs calculated using volume fractions at surfaces were also used to test the reliability of various equations of state of HSCFs by different authors.
文摘This paper addresses the issues of channel estimation in a Multiple-Input/Multiple-Output (MIMO) system. Markov Chain Monte Carlo (MCMC) method is employed to jointly estimate the Channel State Information (CSI) and the transmitted signals. The deduced algorithms can work well under circumstances of low Signal-to-Noise Ratio (SNR). Simulation results are presented to demonstrate their effectiveness.
文摘目的 针对妇幼卫生纵向数据的任意缺失模式,采用多重填补方法进行填补,探求最佳填补结果,以便对数据作进一步分析与研究。方法 运用SAS9.0 ,采用多重填补方法Markov China Monte Carlo(MCMC)模型对缺失数据进行多次填补并综合分析。结果 填补5次所得结果最优。结论 多重填补方法可以处理有缺失数据资料中的许多普遍问题,可提高统计效率,尤其是MCMC模型在处理复杂的缺失数据上,优势明显。
文摘With increasing complexity of today’s electromagnetic problems, the need and opportunity to reduce domain sizes, memory requirement, computational time and possibility of errors abound for symmetric domains. With several competing computational methods in recent times, methods with little or no iterations are generally preferred as they tend to consume less computer memory resources and time. This paper presents the application of simple and efficient Markov Chain Monte Carlo (MCMC) method to the Laplace’s equation in axisymmetric homogeneous domains. Two cases of axisymmetric homogeneous problems are considered. Simulation results for analytical, finite difference and MCMC solutions are reported. The results obtained from the MCMC method agree with analytical and finite difference solutions. However, the MCMC method has the advantage that its implementation is simple and fast.
文摘目的探讨基于Markov Chain Monte Carlo(MCMC)模型的多重估算法在处理医院调查资料缺失数据中的应用。方法运用SAS9.2编写程序,在分析数据的分布类型和缺失机制的基础上,采用MCMC法对缺失数据进行多次填补和联合统计推断,分析多重估算法的优势。结果数据服从多元正态分布与随机缺失,采用MCMC法填补10次所得的结果最佳。结论多重估算既可反映缺失数据的不确定性,又可充分利用现有资料的信息、提高统计效率、对模型的估计结果更加可信,是处理缺失数据的有效方法。
基金Project supported by the Shanghai Leading Academic Disciplines Project,China (Grant No. T0104)the National Natural Science Foundation of China(Grant No. 10674092)
文摘The magnetic properties of a mixed spin-2 and spin-1/2 ferromagnetic diamond chain are studied by effective-field theory and Monte Carlo(MC) simulation based on the Ising model.The temperature dependences of magnetization,magnetic susceptibility,internal energy,and specific heat are studied,respectively.The exchange interaction dependences of magnetization and the critical temperature are calculated by MC simulation.The changes of magnetization depending on the field increasing and then the field decreasing under steady-static conditions are also given.
文摘A program of Monte Carlo simulation of binary copolymerization for E-SBR (emulsion polymn. SB rubber) was made according to the terminal model. The simulation results obtained by this program were in good agreement with those experimental ones. A detail microstructure information of E-SBR molecular chain has been provided.
文摘This paper proposes a new technique based on inverse Markov chain Monte Carlo algorithm for finding the smallest generalized eigenpair of the large scale matrices. Some numerical examples show that the proposed method is efficient.
文摘A part of a long DNA chain was driven into a confined environment by an electric field, while the rest remains in the higher-entropy region. Upon removal of the field, the chain recoils to the higher-entropy region spontaneously. This dynamical process was investigated by Monte Carlo simulations. The simulation reproduces the experimentally-observed phenomenon that the recoil of the DNA chain is initially slow and gradually increases in speed due to the presence of the confinement-entropic force. The results show that with increasing the dimension or decreasing the spacing of the nanopillars the recoil velocity of the DNA chain will increase. Further analysis suggests that the characteristic entropy per monomer in the confinement is proportional to the area fraction of the free part in the confinement.
基金Project (No. 20204014) supported by the National Natural ScienceFoundation of China
文摘The configurational properties of tail-like polymer chains with one end attached to a flat surface are studied by using dynamic Monte Carlo technique. We find that the probability distribution of the free end in z direction P(Rz) and the density profile p(z) can be scaled approximately by a factor β to be a length independent function for both random walking (RW) and self-avoiding walking (SAW) tail-like chains, where the factor β is related to the mean square end-to-end distance 〈RE〉. The scaled P(Rz) of the SAW chain roughly overlaps that of the RW chain, but the scaled p(z) of the SAW chain locates at smaller βz than that of the RW chain.
基金Supported by the National Natural Science Foundation of China under Grant Nos.10674016,10875013the Specialized Research Foundation for the Doctoral Program of Higher Education under Grant No.20080027005
文摘We introduce the potential-decomposition strategy (PDS), which can be used in Markov chain Monte Carlo sampling algorithms. PDS can be designed to make particles move in a modified potential that favors diffusion in phase space, then, by rejecting some trial samples, the target distributions can be sampled in an unbiased manner. Furthermore, if the accepted trial samples are insumcient, they can be recycled as initial states to form more unbiased samples. This strategy can greatly improve efficiency when the original potential has multiple metastable states separated by large barriers. We apply PDS to the 2d Ising model and a double-well potential model with a large barrier, demonstrating in these two representative examples that convergence is accelerated by orders of magnitude.
文摘Recently, a new soft-in soft-out detection algorithm based on the Markov Chain Monte Carlo (MCMC) simulation technique for Multiple-Input Multiple-Output (MIMO) systems is proposed, which is shown to perform significantly better than their sphere decoding counterparts with relatively low complexity. However, the MCMC simulator is likely to get trapped in a fixed state when the channel SNR is high, thus lots of repetitive samples are observed and the accuracy of A Posteriori Probability (APP) estimation deteriorates. To solve this problem, an improved version of MCMC simulator, named forced-dispersed MCMC algorithm is proposed. Based on the a posteriori variance of each bit, the Gibbs sampler is monitored. Once the trapped state is detected, the sample is dispersed intentionally according to the a posteriori variance. Extensive simulation shows that, compared with the existing solution, the proposed algorithm enables the markov chain to travel more states, which ensures a near-optimal performance.
文摘目的:为急性冠脉综合征(ACS)患者抗血小板治疗提供决策依据。方法:收集PLATO和TREAT研究中的相关数据建立Markov模型,计算ACS患者应用替格瑞洛或氯吡格雷的总出血风险、主要出血风险、次要出血风险及致命出血风险等,依据既往文献收集、计算各状态间的转移概率、各状态的成本及效用值,运用TreeAge Pro 2011软件计算不同方案的医疗费用、质量调整生命年(QALYs)以及增量-成本效果比(ICER),并对模型进行单因素敏感性分析及概率敏感性分析。结果:替格瑞洛组平均总成本是66449.38元,获得7.34 QALYs;氯吡格雷组平均总成本是53846.03元,获得6.68 QALYs。替格瑞洛组与氯吡格雷组相比的ICER为19095.98元/QALYs,即每多获得1个QALYs,替格瑞洛组患者需多花19095.98元,低于意愿支付阈值(64644元)。敏感性分析结果与上述结果一致。结论:相较于氯吡格雷,替格瑞洛对ACS患者尤其是发生出血不良事件的ACS患者具有更小的经济负担。
基金supported in part by"National S&T Major Project Foundation of China"(2009ZX10004-904)Universities Natural Science Foundation of Jiangsu Province(09KJB330004),National Science Foundation Grant DMS-9971405National Institutes of Health Contract N01-HV-28183
文摘This paper first applies the sequential cluster method to set up the classification standard of infectious disease incidence state based on the fact that there are many uncertainty characteristics in the incidence course.Then the paper presents a weighted Markov chain,a method which is used to predict the future incidence state.This method assumes the standardized self-coefficients as weights based on the special characteristics of infectious disease incidence being a dependent stochastic variable.It also analyzes the characteristics of infectious diseases incidence via the Markov chain Monte Carlo method to make the long-term benefit of decision optimal.Our method is successfully validated using existing incidents data of infectious diseases in Jiangsu Province.In summation,this paper proposes ways to improve the accuracy of the weighted Markov chain,specifically in the field of infection epidemiology.