Carlo模拟

摘要： 为了度量估计误差,该文在E-Bayes估计(expected Bayesian estimation)的基础上引入了E-MSE(expected mean square error)的定义,并推导不同损失函数(包括平方损失函数和LINEX损失函数)下失效概率的E-Bayes估计及其E-MSE的表达式.通过MonteCarlo模拟比较了所提出估计方法的性能(比较基于E-MSE).最后,分别采用E-Bayes方法和MCMC方法,结合发动机可靠性问题进行了计算和分析.在考虑评价不同损失函数下参数的E-Bayes估计时,该文提出用E-MSE作为评价标准.

摘要： 在海量数据的假设检验中,传统的假设检验方法和错误控制率会导致重要的信号被遗漏掉。文章提出利用数据结构进行分层检验,并结合较为宽松的错误发现率(FDR)进行控制;重点讨论了双层结构检验中检验阈值的选择,对比单层和双层结构检验方法,分析它们的优点并利用Monte Carlo数值模拟它们在错误控制率和检验灵敏度上的表现。

摘要： 在地表水地下水耦合模拟过程中,源汇项的不确定性会影响模型输出结果的可靠性,故对其进行不确定性分析十分必要。针对一假想算例建立了地表水地下水耦合模拟数学模型,并运用Hydro Geo Sphere软件对其进行同步并行求解。采用Monte Carlo方法分析源汇项的不确定性对于地表水地下水耦合模拟输出结果的影响。在不确定性分析过程中,为减少多次调用模拟模型所产生的计算负荷,采用拉丁超立方抽样和Kriging方法建立模拟模型的替代模型。结果表明:地表水地下水耦合模拟模型能够描述研究区地表水流与地下水流的水力联系;采用Kriging方法建立模拟模型的替代模型,能减少不确定性分析的计算负荷;结合风险评估,地下水水位低于地下水生态水位阈值的风险值为0.2;基于替代模型对地表水地下水耦合模拟模型进行不确定性分析,大幅度地减少了计算成本,为决策者的实际工作提供参考。

摘要： 目的:探讨新药在国际多中心临床试验(MRCT)中,以正态资料为主要疗效终点、目标区采用调整后的检验水准α'作为决策依据方法的可行性和α'的合理取值,确定目标区样本量的比例,为MRCT在目标区通过审批提供参考。方法:利用Monte Carlo法模拟优效性设计和目标区不同规模MRCT的数据模型,在检验水准α=0.05水平上显示试验组优于对照组的前提下,目标区采用调整后检验水准α'作为决策依据的条件Ⅰ型错误率(CFPR)、条件检验效能(CP)随目标区样本量占MRCT样本量占比K的变化情况。结果:目标区CFPR和CP随着α'的增加而增加,K值越大,相应的目标区CFPR也越高。当K0.8,则目标区CP可超过0.70;当K=25%,如果f>0.8,目标区CP可超过0.75;当K=30%,只要f=0.7,目标区CP就能维持在0.80以上;当K=50%时,f即使只有0.5,目标区CP依然可以达到0.80以上。结论:本研究提出的方法容易理解,便于操作,尤其对于f≥0.7的情形有较好表现。

摘要： Objective The performance of propensity score(PS)methods were compared through constructing different confounding structures and generating covariates with different correlations when the outcome model was linear.Methods Monte Carlo method was used to simulate the datasets by constructing four confounding structures from simple to complex.Then four PS-based methods including PS matching,covariate adjustment,inverse probability of weighting(IPW)and stratification were applied to estimate the treatment effect.The results were compared from different aspects including the point estimate,standard error,relative bias and mean square error.Results When the outcome model was linear,covariate adjustment showed the least biased and stable estimates among the four methods.PS matching with caliper 0.02 performed better than the other matching methods when the caliper is 0.2 of the standard deviation of the PS value.When there were nonlinear relationship in the treatment model,IPW showed biased results and largest standard error.PS stratification resulted in biased estimates in all settings.Conclusion Covariate adjustment by PS is robust to complex confounding structure and achieved the least biased estimates.We propose that when the relationships between confounding factors and treatment or outcome variable cannot be confirmed,using PS covariate adjustment seems a better choice.%目的 通过构建不同混杂结构的处理因素模型和结局模型、不同相关性的协变量,比较多种倾向性评分方法在结局模型为线性回归模型的情况下估计处理效应的优劣.方法 采用Monte Carlo模拟方法,通过构建四种由简单到复杂的不同结构的混杂模型,生成相应的数据集,再分别应用倾向性评分匹配、回归调整、加权以及分层的方法估计处理效应并进行比较.评价指标包括点估计、标准误、相对偏倚、均方误差.结果 在结局模型为线性回归模型情况下,倾向性评分回归调整法估计的相对偏倚最小,稳定性也最好.匹配法卡钳值取0.02较卡钳值取倾向性评分标准差的0.2倍估计的相对偏倚更小.当处理因素模型中含有非线性效应时,用逆概率加权法估计的偏倚较大,并且加权法估计的标准误也最大.倾向性评分分层法在各种情况下估计的相对偏倚都较大.结论 倾向性评分回归调整法能够较好地估计处理效应,并且在各种情况下估计都较为稳健.建议当协变量与处理因素和结局变量的关系无法确定时,这四种方法中可以考虑优先使用回归调整法.

摘要： 在金融风险评估、事故预测、保险索赔等领域的研究中,极值理论已发展成为一种重要的统计学方法。Gumbel分布是一种常用的极值分布函数,并逐渐成为了对于随机变量极端变异性建模的重要工具。文章将二项分布与Gumbel分布函数复合,提出了一种新的复合极值分布函数即二项-Gumbel分布。重点介绍了极值理论以及二项分布与Gumbel分布复合函数,运用极大似然估计(MLE)对二项-Gumbel复合分布的各种参数进行估计,并通过计算机模拟得KS检验统计量的临界值。

摘要： 利用三参数威布尔分布模型对采集终端的寿命分布进行建模,并用Monte Carlo模拟方法分别比较了由极大似然法,相关系数优化法和灰色估计法拟合得到的威布尔分布的三参数的拟合特性,得出了这3种方法的拟合精度和适用范围。相关结论对正确选用威布尔参数估计方法评估采集终端的可靠性具有一定的指导作用。

摘要： 测量了10—25 keV电子碰撞厚W,Au靶产生的韧致辐射谱,并与Monte Carlo程序PENELOPE模拟的X射线谱进行了比较,除在3 keV前实验谱略低于理论谱外,整体上两者符合得很好.在模拟电子与靶材料相互作用产生韧致辐射时,PENELOPE程序中只包含有普通韧制辐射的截面数据.实验结果表明,在电子与固体靶相互作用时,没有明显的极化韧致辐射产生,PENELOPE程序能够可靠地描述电子与固体厚靶相互作用产生的韧致辐射.

摘要： 研究模拟均匀、各向同性介质内的逆康普顿散射过程的Monte Carlo方法,分析电子速度分布乘抽样法和光子散射方向分布乘抽样法的抽样效率,完善电子速度分布的直接抽样法、光子散射能量分布的乘加抽样法.通过对这四种抽样方法抽样费用的分析和比较,得出各种方法的最适条件,该方法可以模拟更高能的辐射.

摘要： 空间误差分量模型(Spatial Error Components,SEC)传统的空间相关性LM检验存在严重的水平扭曲和较低的检验功效,导致检验统计量失效。文章将Bootstrap方法应用于SEC模型的空间相关性LM检验,提高检验统计量的有效性。Monte Carlo模拟实验表明,Bootstrap LM检验的水平受误差项分布、空间权重矩阵和样本量影响较小,并且远优于渐近LM检验,具有理想的检验水平;渐近LM检验和Bootstrap LM检验的功效均随着空间相关性的增强,及样本量的增大而增大,但Bootstrap LM检验在各种情形下均具有更高的检验功效,尤其是样本量较小时。简言之,Bootstrap LM检验是SEC模型更为优越的空间相关性检验方法。

摘要： 单粒子翻转效应的研究目的在于研究计算机CPU和存储器受高能中子、质子等粒子辐照后软错误发生的概率，为计算机系统抗辐射加固提供数值模拟方法和参考数据。本文从中子与硅原子相互作用的物理机理出发，引入随机数，利用MonteCarlo方法编制了中子引起的单粒子翻转截面的计算模拟程序。本文介绍的单粒子翻转的Montc Carlo计算模拟程序可为14 MeV中子环境下的16K位静态存储器硅片翻转提供截面和描述内部物理过程的参考数据。

摘要： 单粒子翻转效应的研究目的在于研究计算机CPU和存储器受高能中子、质子等粒子辐照后软错误发生的概率，为计算机系统抗辐射加固提供数值模拟方法和参考数据。本文从中子与硅原子相互作用的物理机理出发，引入随机数，利用MonteCarlo方法编制了中子引起的单粒子翻转截面的计算模拟程序。本文介绍的单粒子翻转的Montc Carlo计算模拟程序可为14 MeV中子环境下的16K位静态存储器硅片翻转提供截面和描述内部物理过程的参考数据。

摘要： 单粒子翻转效应的研究目的在于研究计算机CPU和存储器受高能中子、质子等粒子辐照后软错误发生的概率，为计算机系统抗辐射加固提供数值模拟方法和参考数据。本文从中子与硅原子相互作用的物理机理出发，引入随机数，利用MonteCarlo方法编制了中子引起的单粒子翻转截面的计算模拟程序。本文介绍的单粒子翻转的Montc Carlo计算模拟程序可为14 MeV中子环境下的16K位静态存储器硅片翻转提供截面和描述内部物理过程的参考数据。

摘要： 单粒子翻转效应的研究目的在于研究计算机CPU和存储器受高能中子、质子等粒子辐照后软错误发生的概率，为计算机系统抗辐射加固提供数值模拟方法和参考数据。本文从中子与硅原子相互作用的物理机理出发，引入随机数，利用MonteCarlo方法编制了中子引起的单粒子翻转截面的计算模拟程序。本文介绍的单粒子翻转的Montc Carlo计算模拟程序可为14 MeV中子环境下的16K位静态存储器硅片翻转提供截面和描述内部物理过程的参考数据。

摘要： 本发明公开了一种基于Copula和Monte‑Carlo模拟的路径行程时间估计方法，步骤为：S1、获取一条路径中各路段路径行程时间波动率；S2、通过最大似然估计对各路段路径行程时间波动率进行边缘分布拟合；S3、基于各路段路径行程时间波动率分布函数使用Copula理论进行路径行程时间波动率分布函数的拟合；S4、在Copula函数的基础上进行Monte‑Carlo模拟以获得路径行程时间波动率序列；S5、确定某一置信水平下的置信区间并计算路径行程时间波动率的期望；S6、把路径行程时间波动率期望值与上一时段的路径行程时间结合以估计当前时段的路径行程时间。本发明考虑了各路段行程时间在时序上的波动，还考虑了路段行程时间之间的空间关系，深度挖掘行程时间的时空特征，具有较高的精度和可靠性。

摘要： 本发明公开了一种适于研究ZnMgO/ZnO异质结中合金群散射的MonteCarlo模拟方法，包括：1)计算得ZnO和Zn1‑xMgxO的能带结构，采用五能谷解析带近似；2)拟合得ZnO和Zn1‑xMgxO最低五个能谷的有效质量；3)得到ZnMgO/ZnO异质结的电子波函数、量子化能级和电子面密度；4)建立模拟ZnMgO/ZnO异质结输运特性的MC模型；5)计算各种散射机制的散射率；6)初始化各粒子的波矢；7)设置粒子编号n和电场强度F；8)n加1；9)判断粒子编号与模拟粒子总数的关系；10)判断模拟时间是否为总仿真时间；11)计算电子的稳态漂移速度和电子迁移率；12)绘制合金群散射对电子输运特性影响的关系图。得到的电子输运特性更加精确，为在ZnMgO/ZnO异质结中降低合金群散射的影响和提高电子迁移率特性提供了参考。

摘要： 本发明涉及一种模拟导线边，属于教学工具领域，用于解决教学中对于导线模拟困难的技术问题，技术要点是：主要由滑动部和具有磁性的滑道组成，在所述滑道中成型有与所述滑动部相互配合的滑槽，该滑槽使滑动部能沿滑道的长度方向移动，在滑动部的一个端部上置有与滑动部一体成型的夹持部，与该夹持部相对的滑道的端部上延伸出用于安装连接装置的连接部。有益效果是得到了一种导线模拟装置。

摘要： 本发明提供一种通过Ewald sum的Monte Carlo分子模拟计算苯与水液相界面张力的方法，包括如下几个步骤：步骤A:初始化构象：准备两个模拟盒子，大小均为(30×30×45)–(35×35×55)。将1200‑2000个水放在第一个盒子的格点上，230‑380个苯分子放在第二个盒子的格点上；步骤B:挑选模拟所用力场：苯分子采用TraPPE力场，水分子使用TIP4P/2005力场；步骤C:在正则系综下对两个盒子分别进行Monte Carlo模拟；步骤D:将水盒子与苯盒子沿着z方向拼接在一起，形成一个（30×30×90Å）‑（35×5×110Å）的大盒子，并Monte Carlo cycle 50000‑80000次，以使体系重新达到平衡；步骤E:体系达到平衡后，进入数据收集阶段。本发明除了提供分子模拟常用的流程外，改进了Ewald sum部分的计算方式，使用本发明的Ewald sum计算方法后，总体模拟效率可以提升约20%。

摘要： 本发明公开了一种基于Copula和Monte‑Carlo模拟的路径行程时间估计方法，步骤为：S1、获取一条路径中各路段路径行程时间波动率；S2、通过最大似然估计对各路段路径行程时间波动率进行边缘分布拟合；S3、基于各路段路径行程时间波动率分布函数使用Copula理论进行路径行程时间波动率分布函数的拟合；S4、在Copula函数的基础上进行Monte‑Carlo模拟以获得路径行程时间波动率序列；S5、确定某一置信水平下的置信区间并计算路径行程时间波动率的期望；S6、把路径行程时间波动率期望值与上一时段的路径行程时间结合以估计当前时段的路径行程时间。本发明考虑了各路段行程时间在时序上的波动，还考虑了路段行程时间之间的空间关系，深度挖掘行程时间的时空特征，具有较高的精度和可靠性。

摘要： 本发明公开了一种适于研究ZnMgO/ZnO异质结中合金群散射的MonteCarlo模拟方法，包括：1)计算得ZnO和Zn1‑xMgxO的能带结构，采用五能谷解析带近似；2)拟合得ZnO和Zn1‑xMgxO最低五个能谷的有效质量；3)得到ZnMgO/ZnO异质结的电子波函数、量子化能级和电子面密度；4)建立模拟ZnMgO/ZnO异质结输运特性的MC模型；5)计算各种散射机制的散射率；6)初始化各粒子的波矢；7)设置粒子编号n和电场强度F；8)n加1；9)判断粒子编号与模拟粒子总数的关系；10)判断模拟时间是否为总仿真时间；11)计算电子的稳态漂移速度和电子迁移率；12)绘制合金群散射对电子输运特性影响的关系图。得到的电子输运特性更加精确，为在ZnMgO/ZnO异质结中降低合金群散射的影响和提高电子迁移率特性提供了参考。

摘要： 本发明公开了一种n型Si材料中电子输运问题的Monte Carlo模拟方法，其步骤如下：一、载流子散射机制的确定以及对应输入条件下各种散射率的计算；二、载流子漂移模型的建立以及载流子漂移后能量与波矢量的计算；三、载流子散射模型的建立以及散射类型的选择；四、Monte Carlo方法模拟n型半导体Si材料中电子的输运问题计算程序的实现。本发明使得计算n型Si材料的平均速率以及迁移率变得简单快捷，避免了之前使用实验测试方法受到半导体器件尺寸以及实验条件影响造成的难测试以及误差较大等一系列问题，而且该方法具有较强的可推广性，其他半导体材料载流子输运或微观粒子的碰撞问题也可以通过改变对应的输入参数来进行计算。

摘要： 利用Monte Carlo方法模拟双能X射线成像进行物质识别的方法，涉及一种双能X射线成像进行物质识别的方法。它是为了解决目前在安检过程中，物质识别准确率低的问题。其方法：步骤一、利用Monte Carlo方法模拟能量为9MeV和能量为6MeV的电子束轰击钨靶过程，产生轫致辐射，获得X射线谱；步骤二、对步骤一获得的X射线谱进行准直，准直后扫描待测样品，获得扫描图像；步骤三、分析步骤二获得的扫描图像，获得待测样品的识别曲线，根据该识别曲线与物质原子之间的关系，从而实现物质识别。本发明适用于安防检查。

摘要： 本发明涉及一种模拟导线边，属于教学工具领域，用于解决教学中对于导线模拟困难的技术问题，技术要点是：主要由滑动部和具有磁性的滑道组成，在所述滑道中成型有与所述滑动部相互配合的滑槽，该滑槽使滑动部能沿滑道的长度方向移动，在滑动部的一个端部上置有与滑动部一体成型的夹持部，与该夹持部相对的滑道的端部上延伸出用于安装连接装置的连接部。有益效果是得到了一种导线模拟装置。

摘要： 本发明涉及照明灯具光学设计技术领域，尤其涉及一种基于尺寸公差Monte‑Carlo分析的灯具配光设计方法，包括：步骤一、通过光学设计模拟建立关键间距与光束角的关联关系；步骤二、确定与关键间距相关的多个零件参数及多个零件参数的尺寸；步骤三、对多个零件参数的尺寸公差工艺控制能力进行评估，确定每个零件参数的尺寸公差分布；步骤四、在多个零件参数的尺寸公差分布范围内，通过Monte‑Carlo统计方法计算获得关键间距的概率分布；步骤五、根据关键间距与光束角的关联关系和关键间距的概率分布，推算光束角是否满足配光性能设计要求。可预测零件生产工艺及工艺控制能力能否达到配光性能设计要求，为零件尺寸调整提供理论分析指导。