SEMI-BLIND CHANNEL ESTIMATION OF MULTIPLE-INPUT/MULTIPLE-OUTPUT SYSTEMS BASED ON MARKOV CHAIN MONTE CARLO METHODS

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.

Seismic fluid identification using a nonlinear elastic impedance inversion method based on a fast Markov chain Monte Carlo method

Elastic impedance inversion with high efficiency and high stability has become one of the main directions of seismic pre-stack inversion. The nonlinear elastic impedance inversion method based on a fast Markov chain Monte Carlo （MCMC） method is proposed in this paper, combining conventional MCMC method based on global optimization with a preconditioned conjugate gradient （PCG） algorithm based on local optimization, so this method does not depend strongly on the initial model. It converges to the global optimum quickly and efficiently on the condition that effi- ciency and stability of inversion are both taken into consid- eration at the same time. The test data verify the feasibility and robustness of the method, and based on this method, we extract the effective pore-fluid bulk modulus, which is applied to reservoir fluid identification and detection, and consequently, a better result has been achieved.

多重填补法Markov Chain Monte Carlo模型在有缺失值的妇幼卫生纵向数据中的应用

目的 针对妇幼卫生纵向数据的任意缺失模式,采用多重填补方法进行填补,探求最佳填补结果,以便对数据作进一步分析与研究。方法 运用SAS9.0 ,采用多重填补方法Markov China Monte Carlo(MCMC)模型对缺失数据进行多次填补并综合分析。结果 填补5次所得结果最优。结论 多重填补方法可以处理有缺失数据资料中的许多普遍问题,可提高统计效率,尤其是MCMC模型在处理复杂的缺失数据上,优势明显。

AN IMPROVED MARKOV CHAIN MONTE CARLO METHOD FOR MIMO ITERATIVE DETECTION AND DECODING

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.

Evaluation of the Reliability of a System: Approach by Monte Carlo Simulation and Application

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. .

基于Markov Chain Monte Carlo模型对医院调查资料中缺失数据的多重估算

目的探讨基于Markov Chain Monte Carlo(MCMC)模型的多重估算法在处理医院调查资料缺失数据中的应用。方法运用SAS9.2编写程序,在分析数据的分布类型和缺失机制的基础上,采用MCMC法对缺失数据进行多次填补和联合统计推断,分析多重估算法的优势。结果数据服从多元正态分布与随机缺失,采用MCMC法填补10次所得的结果最佳。结论多重估算既可反映缺失数据的不确定性,又可充分利用现有资料的信息、提高统计效率、对模型的估计结果更加可信,是处理缺失数据的有效方法。

Bayesian Markov chain Monte Carlo inversion for anisotropy of PP-and PS-wave in weakly anisotropic and heterogeneous media

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.

On Finding the Smallest Generalized Eigenpair Using Markov Chain Monte Carlo Algorithm

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.

Potential-Decomposition Strategy in Markov Chain Monte Carlo Sampling Algorithms

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.

Markov Chain Monte Carlo Solution of Laplace’s Equation in Axisymmetric Homogeneous Domain

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.

基于Monte Carlo与Markov法参数不确定性条件下SIL评估

为了避免由参数不确定性因素导致较大的Markov法SIL评估偏差,且减少计算工作量和复杂性,采用Monte Carlo(MC)仿真方法处理含不确定性参数的Markov模型,并借助Matlab GUI编程开发MC仿真处理参数不确定性条件下Markov法SIL评估可视化仿真计算软件。在理论研究的基础上,为说明该研究方法与计算软件的可行性,以石油天然气工业高完整性压力保护系统(HIPPS)为算例进行SIL评估。结果表明:MC仿真方法可以有效处理Markov法SIL评估中参数不确定性问题;基于Matlab GUI编程设计出的仿真计算软件在一定程度上可以提高计算效率。

结合贝叶斯-MCMC更新的RC梁抗剪承载力概率模型的可靠度分析

为分析钢筋混凝土(RC)梁抗剪承载力,首先,基于国内外100组RC梁抗剪承载力试验数据作为先验信息,结合我国GB 50010-2010规范的抗剪承载力,综合考虑混凝土与箍筋的材料强度、截面尺寸、剪跨比、配箍率等因素的影响,通过贝叶斯-马尔科夫链蒙特卡洛(MCMC)方法对我国规范模型进行更新与修正;随后,结合Monte Carlo模拟,对所建立的RC梁概率抗剪承载力模型进行可靠度分析,验证了该模型具有较好的计算精度与可靠性。结果表明:概率模型均值与试验值比值K的均值与标准差分别为1.013与0.171,试验值落在概率模型的95%置信区间范围内,且可靠指标分布在4.0左右,说明建立的概率模型具有较好的预测效果与可靠性。

基于爬坡方向状态划分的MCMC风电功率序列建模方法

由于电网弃风或者灵活性资源不足往往发生在风电大量发电时,故提高风电时间序列模型对大出力状态的建模-抽样精度,有助于后续的电网灵活性资源相关研究。在传统马尔科夫链蒙特卡洛(MCMC)法和持续与波动蒙特卡罗(PV-MC)法基础上,提出一种考虑爬坡方向状态划分的改进方法,以更准确地描述风电出力连续爬坡至大出力状态的过程。该方法以累积分布概率而不是以功率大小均匀划分状态区间,使各个状态区间的样本分布更均匀,提高了风电时间序列模型对大出力状态的建模-抽样精度。通过算例比较所提方法、MCMC法及PV-MC法生成风电功率序列与历史数据的分布特性和统计特性指标,结果表明,所提方法的拟合度较好,且能够有效解决MCMC法和PV-MC法高出力、样本偏少的问题。

基于MCMC的电网安全稳定控制系统动态可靠性评估方法

现有安全稳定控制系统(简称稳控系统)的可靠性评估方法本质上属于静态建模,由于未能体现系统内各装置老化和检修等动态过程,在一定程度上影响了评估结果的准确性。为此,文中提出一种基于马尔可夫链蒙特卡洛(Markov chain Monte Carlo,MCMC)的稳控系统动态可靠性评估方法。首先针对失效过程,构建四状态非齐次马尔可夫模型来模拟装置老化过程,并给出各状态评判方法;其次针对修复过程,分析不同检修策略对装置状态转移的影响以体现状态检修的差异性;最后考虑稳控装置状态转移过程的时序或条件相关性,对稳控系统可靠性进行动态建模。以实际稳控系统为例,仿真对比不同检修策略下的可靠性,并对模型参数进行灵敏度分析。评估结果表明,该方法可以求解稳控系统的时变可用度,用于指导稳控装置现场合理检修。

Monte Carlo Simulations of Density Profiles for Hard-Sphere Chain Fluids Confined Between Surfaces

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.

Multi-Scale Variation Prediction of PM2.5 Concentration Based on a Monte Carlo Method

Haze concentration prediction,especially PM2.5,has always been a significant focus of air quality research,which is necessary to start a deep study.Aimed at predicting the monthly average concentration of PM2.5 in Beijing,a novel method based on Monte Carlo model is conducted.In order to fully exploit the value of PM2.5 data,we take logarithmic processing of the original PM2.5 data and propose two different scales of the daily concentration and the daily chain development speed of PM2.5 respectively.The results show that these data are both approximately normal distribution.On the basis of the results,a Monte Carlo method can be applied to establish a probability model of normal distribution based on two different variables and random sampling numbers can also be generated by computer.Through a large number of simulation experiments,the average monthly concentration of PM2.5 in Beijing and the general trend of PM2.5 can be obtained.By comparing the errors between the real data and the predicted data,the Monte Carlo method is reliable in predicting the PM2.5 monthly mean concentration in the area.This study also provides a feasible method that may be applied in other studies to predict other pollutants with large scale time series data.

利用Monte Carlo技术模拟研究不同缺失值处理方法对完全随机缺失数据的处理效果

目的 以医疗卫生机构年报资料为数据来源,采用成组删除法、极大似然估计法、多重填补法分别对模拟的完全随机缺失数据集缺失值进行处理,比较不同缺失率下三种方法的缺失处理效果。方法 运用SAS9.3,采用Monte Carlo技术模拟完整数据集及不同缺失比例数据集,利用成组删除法、EM算法、MCMC算法对缺失数据进行处理,得到不同处理方法后的参数估计结果,与完整数据集参数估计进行比较。结果 对于完全随机缺失数据,不同缺失率下,成组删除法的准确率均比较好;缺失率小于10%,三种方法处理效果差异不大;缺失率在10%-30%,成组删除法精确度逐渐降低,EM与MCMC准确度与精确度较好,缺失率大于30%,MCMC准确度与精确度相对较好。结论 对于不同缺失率的数据,综合考虑准确度和精确度,采用不同的方法进行处理。

基于Monte Carlo方法的一对二马尔可夫随机格斗模型

讨论了一对二马尔可夫随机格斗双方获胜概率计算问题。提出了一种新颖的一对二马尔可夫随机格斗任意对抗回合双方获胜概率的计算方法,该方法首先基于Monte Carlo仿真计算各个对抗回合中双方发射次序的概率分布,再利用全概率公式确定马尔可夫链的状态转移概率矩阵,从而克服了马尔可夫随机格斗模型往往只能提供无限对抗回合之后格斗双方获胜概率的缺点,为运用马尔可夫随机格斗研究火力运用和弹药分配提供了新途径,并用实例说明了该方法的有效性。

利用Monte Carlo方法求解线性抛物型问题(英文)

提出一种一维线性抛物型偏微分方程的温度分布函数的数值解法,数值算法是基于在空间和时间上采用紧有限差分法(CFD)得到离散化的控制方程进而利用Monte Carlo(MC)随机模拟方法求解所得的方程.通过比较由CFD方法和有限差分法(FD)得到的数值解与精确解的误差的计算结果说明了所提方法的效率和精度.

基于MCMC填补的SSA-SVM煤与瓦斯突出预测模型

为提升煤与瓦斯突出预测准确度,减小数据缺失对煤与瓦斯突出预测的不利影响,提出1种基于链式多重填补马尔科夫链蒙特卡罗(MCMC)的麻雀搜索算法(SSA)优化支持向量机(SVM)预测模型。根据突出影响因素选取模型参数,运用MCMC对突出事故缺失值进行数据填补,采用SSA优化SVM,建立MCMC-SSA-SVM模型对填补后数据集进行预测,验证MCMC填补有效性和SSA优化性能;分别构建SVM、SSA-SVM、PSO-SVM、GAM-SVM、CMC-SVM、MCMC-PSO-SVM和MCMC-GA-SVM这7种模型进行突出预测,对比预测准确度,分析MCMC-SSA-SVM、MCMC-PSO-SVM和MCMC-GA-SVM的适应度。研究结果表明:MCMC填补后准确度均提升7.89个百分点以上,SSA的优化性能强于PSO和GA,MCMC-SSA-SVM预测准确度最高,为97.37%,泛化能力优于对比模型。研究结果可为煤与瓦斯突出预测研究提供借鉴和参考。