约束路径Monte Carlo方法的倍维对称表示

为了研究我们提出的附相t-U Hubbard模型中一般相角的选取对模型电子配对关联性质的影响,在传统的约束路径Monte Carlo方法的基础上我们发展了用来研究体系一般厄米哈密顿量基态性质的这种模拟方法的倍维对称表示．     

附相t-t''''-U Hubbard 模型的约束路径Monte Carlo方法研究

利用约束路径Monte Carlo方法的双维表示，对附相t-t’-U Hubbard模型进行了研究．在研究中，我们发现，当在t-t’-U Hubbard模型中加入一体现内部非均匀性的相因子时，不同的相因子对扩展的s-波，d-波的配对关联函数和顶角贡献有一定的影响．而在相同的相角时，负的t‘值减小了关联函数顶角贡献的短程关联性．     

Tobit方差分量模型的极大似然估计

本文对非负的且含有大量零的混合类型数据提出了Tobit方差分量模型,许多很有用的Tobit模型是我们模型的特例.我们运用MCEM算法给出了模型的极大似然估计,其中E-步运用了Gibbs抽样的Monte Carlo模拟,并用Louis方法得到参数的标准误差估计.     

体素模型在Monte Carlo模拟计算中的描述

体素模型与Monte Carlo模拟计算方法相结合已开始应用于辐射防护的研究中.然而,通常的模拟计算需要很长的时间,这限制了体素模型的广泛应用.该文对模拟计算中体素模型的描述进行研究,实现了把三维体素合并算法用在MCNP(Monte Carlo N-particle)程序描述体素模型的方法.结果表明,使用该方法描述体素模型,粒子输运过程要比传统的方法快32%,但是对粒子在组织器官中沉积能量的纪录需要更多时间.     

体素模型在Monte Carlo模拟计算中的描述

体素模型与Monte Carlo模拟计算方法相结合已开始应用于辐射防护的研究中.然而,通常的模拟计算需要很长的时间,这限制了体素模型的广泛应用.该文对模拟计算中体素模型的描述进行研究,实现了把三维体素合并算法用在MCNP(Monte Carlo N-particle)程序描述体素模型的方法.结果表明,使用该方法描述体素模型,粒子输运过程要比传统的方法快32%,但是对粒子在组织器官中沉积能量的纪录需要更多时间.     

Poisson自回归模型的平稳性检验

利用经验似然方法研究Poisson自回归模型的平稳性检验问题, 通过构造相关检验的经验似然比检验统计量, 在零假设下获得了检验统计量的极限分布, 并通过Monte Carlo数值模拟考察该检验方法的有限样本特征.     

基于M-H抽样算法的贝叶斯Probit分位回归模型研究

针对Probit分位回归在参数随机化条件下的建模问题,提出基于Metropolis-Hastings算法的贝叶斯Probit分位回归模型.通过分析Probit分位回归模型结构,选择模型的先验分布,运用M-H算法进行参数估计.利用Monte Carlo仿真技术,得到不同分位点模型参数后验分布,同时用贝叶斯probit分位回归与分位回归方法和光滑分位回归方法对模型参数估计进行比较分析.研究结果表明:贝叶斯Probit分位回归模型可以更全面描述离散选择变量的影响,能够得到更加准确有效的参数估计.     

基于MCMC算法的贝叶斯面板单位根检验

针对面板单位根检验存在检验势不稳定和原假设设置主观选择的问题,提出基于面板数据分位自回归模型,选择非对称Laplace分布的似然函数对模型进行贝叶斯分位回归分析.结合参数的完全条件分布设计MCMC抽样算法,进行贝叶斯分位单位根检验,并利用Monte Carlo模拟实验研究了贝叶斯分位单位根检验的有效性与可行性.研究结果表明,基于面板数据分位自回归模型的贝叶斯单位根检验方法解决了检验势不稳定以及原假设主观设置的问题,能够给出更全面稳健的单位根检验判断.     

带有协变量动态可靠性模型参数的估计

考虑带有协变量的k-out-of-n系统, 假设系统中正常工作变量剩余寿命的条件风险函数为成比例的风险模型, 讨论该模型参数、 系统的可靠度函数和分位数函数的最大似然估计以及这些估计量的渐近正态性, 给出模型参数的置信区间, 并对点估计和区间估计进行Monte Carlo模拟研究. 结果表明, 估计量的均方误差随样本容量的增加而减少.     

一维亚铁磁无序混合自旋链的磁特性的Monte Carlo模拟

基于一维自旋链模型,采用Monte Carlo方法模拟研究了具有反铁磁相互作用的自旋S-1和自旋S-3/2混合的亚铁磁系统的磁特性,重点研究了自旋的无序排列对此混合系统磁特性的影响.研究发现,基态时,当自旋S-1和自旋S-3/2交替混合(有序排列)时,系统存在M-H磁化曲线的阶梯效应;而当自旋S-1和自旋S-3/2无序混合时,系统的阶梯效应出现涨落,台阶甚至消失.最后,通过对系统自旋组态和能量的研究,解释了该混合系统中所产生的阶梯效应和相图.     

用极限分析方法计算边坡可靠度上限值

为了减小边坡可靠度计算在模型建立、参数估计和算法选取上的误差,建立了基于极限分析上限定理的边坡可靠度分析模型,采用最大熵原理推断边坡结构参数的统计特征,并采用Monte Carlo模拟法计算边坡的结构可靠度上限值。最后将结果与其他方法进行了比较。研究结果表明：利用最大熵原理、Monte Carlo模拟法和基于上限定理的边坡可靠度分析模型所计算的可靠度值是边坡可靠度的一个严格上限值,该模型与极限平衡法模型、R.Jimenez-Rodriguez等模型计算的结果比较接近。     

基于三参数Logistic模型Gibbs抽样方法的敏感度分析

通过蒙特卡罗(Monte Carlo,MC)方法研究了三参数Logistic(three-parameter Logistic,3PL)模型的Gibbs抽样方法的估计效果.首先,介绍了 3PL模型的定义和参数的先验分布;其次,介绍了潜变量的引进及Gibbs抽样过程;最后,进行了模拟实验,考虑了不同的测试长度n=10,2...     

零温Hohenberg-Mermin-Wagner定理和Goldstone 定理的等效性

应用零温情形下的Bogoliubov不等式,研究了一维和二维Hubbard模型基态中电子配对和磁性(铁磁和反铁磁)长程序的可能性．发现如果在一维和二维Hubbard模型的自旋(电荷)激发谱中存在能隙,则该系统在基态中不可能呈现出磁性(或电子配对)的长程序,并由此从低维Hub—bard模型说明了零温Hohenberg—Mermin—Wagner定理与Goldstone定理是等效的．     

Monte Carlo方法在放射治疗计划剂量计算中的应用

为了在放射治疗计划系统中的临床应用,探讨Monte Carlo剂量计算方法。应用已经建立Monte Carlo虚拟源模型,针对一例食道肿瘤病例,完成该食道癌患者放射治疗计划各个射野的剂量计算,并进行了计划设计。研究发现,高剂量区与计划靶区(PTV)保持了很好的适形,关键器官的照射量很小,得到了很好的保护。结果表明,用Monte Carlo虚拟源剂量计算模型生成的剂量分布符合临床放射治疗的实施原则,治疗计划靶区得到了高剂量照射,危险器官受照剂量远低于耐受剂量。由于该模型的精确性,完全可以替代现有的解析型剂量计算模型。     

C11集装箱船参数横摇运动极值响应分析

基于窄带随机过程理论和Hermite变换法,本文提出了一种预测船舶参数横摇极值响应的新方法.以C11集装箱船为例,预测了其随机参数横摇极值响应的均值.将该结果与基于Monte Carlo法的数值结果对比,误差小于1%,证明了本文方法的正确性.同时,本文方法使用20条时间历程样本的预测精度与Monte Carlo法使用10~4条时间历程样本的预测精度相同,由此可以证明本文方法的高效性.本文还利用传统Gumbel模型对C11船参数横摇极值响应进行了预测.通过比较可以发现传统Gumbel模型的极值预测误差很大,证明传统Gumbel模型不适用于参数横摇这类非线性较强的运动极值预测.但是,即使利用新方法预测极值,也和模型试验结果之间存在一定的偏差.通过分析认为,这种预测误差是由于忽略了极大值之间的相关性导致的.     

左截断右删失数据下伽玛分布参数多变点的贝叶斯估计

首先通过添加数据得到了左截断右删失数据下伽玛分布的完全数据似然函数,然后研究了变点位置和其它参数的满条件分布,接着利用Gibbs抽样与Metropolis-Hastings算法相结合的MCMC(Markov Chain Monte Carlo)方法得到了参数的Gibbs样本,把Gibbs样本的均值作为各参数的贝叶斯估计,随机模拟试验的结果表明各参数贝叶斯估计的精度都较高.     

中子时间常数的Monte Carlo计算方法

为了研究核装置内中子密度随时间的瞬态变化规律,在Monte Carlo程序DSMC(dynamic system Monte Carlo)的基础上发展了一种新的计算核装置中子时间常数方法,该方法先求得本征分布源,然后求中子时间常数。通过与MCNP(Monte Carlo N-particle)程序计算α本征值的方法及两个模型的计算结果比对分析发现:在超临界和临界情况下两者结果一致;在较深次临界情况下,新方法克服了MCNP程序方法的不足,计算结果更为可靠,特别是在MCNP4C不能顺利完成计算的情况下,新方法依然有效。     

附相t-U Hubbard模型的约束路径Monte Carlo方法研究

作为研究强关联电子体系的一个重要数值模型,单带t UHubbard模型已经被人们广泛地用来研究材料的磁性,相变以及高温超导性等各种性质.借助于各种计算方法,人们在单带以及三带Hubbard模型中发现了关联函数的存在,然而人们同时也发现这些关联函数随着晶格尺度的增加而迅速衰减[1].考虑到目前一个被人们普遍接受的观点:即高温超导性与强关联效应和强磁性耦合相关,我们在传统的单带t UHubbard模型中引入一相因子θ,并将传统的Hubbard模型调整为附相t UHubbard模型.在对晶格施加了周期性边界条件之后,我们利用数值模拟技术—约束路径MonteCar…     

gl(2|1)超代数的齐次微分实现和玻色-费米实现

研究了gl(2/1)超代数在齐次多项式空间上的参数齐次微分实现及其相应的玻色－费米实现。该参数与关联电子Hubbard模型的Hubbard作用参数U相关联。     

基于Monte Carlo模拟设计MicroSPECT针孔准直器

为使小动物成像专用的MicroSPECT(micro single photon emission computed tomography)达到预期的分辨率和灵敏度指标,基于Monte Carlo方法对针孔准直器和成像系统参数进行设计并确定了两套设计方案。方案A保证灵敏度不低于9.19×10-6s-1.Bq-1,且视野中心位置的分辨率不差于1.0 mm;方案B在保证分辨率好于1.2 mm的前提下,可达到5.86×10-5s-1.Bq-1的灵敏度。进一步采用Monte Carlo方法估计了γ光子的准直器孔壁穿透效应和Compton散射效应对投影数据的影响。结果表明:穿透效应在方案A和B中对投影数据的贡献分别为28.6%和40.4%,而散射光子的贡献分别为12.4%和12.9%。