平码五不中公式规律
  • / 21
  • 下载费用:30 金币  

一种模拟能谱曲线仿真方法.pdf

关 键 ?#21097;?/dt>
一种 模拟 谱曲 仿真 方法
  专利查询网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
摘要
申请专利号:

CN201510539124.1

申请日:

2015.08.28

公开号:

CN106483550A

公开日:

2017.03.08

当前法律状态:

授权

?#34892;?#24615;:

有权

法?#19978;?#24773;: 授权|||专利申请权的转移IPC(主分类):G01T 1/36登记生效日:20180821变更事项:申请人变更前权利人:易良碧变更后权利人:成都理工大学变更事项:地址变更前权利人:213000 江苏省常州?#34892;?#21271;区河海大学变更后权利人:610000 四川省成?#38469;?#20108;仙桥东三路1号|||实质审查的生效IPC(主分类):G01T 1/36申请日:20150828|||公开
IPC分类号: G01T1/36 主分类号: G01T1/36
申请人: 易良碧
发明人: 易良碧; 余国刚; 王礼; 涂小芳
地址: 213000 江苏省常州?#34892;?#21271;区河海大学
优?#28909;ǎ?/td>
专利代理机构: 南京同泽专利事务所(特殊普通合伙) 32245 代理人: 蒋全强
PDF完整版下载: PDF下载
法律状态
申请(专利)号:

CN201510539124.1

授权公告号:

|||||||||

法律状态公告日:

2018.09.18|||2018.09.07|||2017.04.05|||2017.03.08

法律状态类型:

授权|||专利申请权、专利权的转移|||实质审查的生效|||公开

摘要

本发明涉及一种模拟能谱曲线仿真方法,包括如下步骤:步骤S1,获取实际核能谱曲线图;步骤S2,对实际核能谱曲线图进行处理以获得模拟能谱曲线;本发明的模拟能谱曲线仿真方法通过对实际核能谱曲线图进行曲线识别并将能谱曲线数值化,从而得到能谱曲线各点的数值(即核能谱的能级和各能级的计数率),再通过蒙特卡罗方法随机直接抽样这一组数值以得到关于各核能级的随机数,从而模拟核衰变过程的随机性,最终再对所述随机数进行统计处理得到模拟能谱曲线,通过反演对比模拟能谱曲线与实际能谱曲线以确定仿核信号发生器的可靠性与精确性。

权利要求书

1.一种模拟能谱曲线仿真方法,其特征在于,包括如下步骤:
步骤S1,获取实际核能谱曲线图;
步骤S2,对实际核能谱曲线图进行处理以获得模拟能谱曲线。
2.根据权利要求1所述的模拟能谱曲线仿真方法,其特征在于,所述模拟
能谱曲线仿真方法还包括:
步骤S3,通过反演比较模拟能谱曲线与实际能谱曲线,以获得模拟能谱曲
线与实际能谱曲线之间的误差。
3.根据权利要求1或2所述的模拟能谱曲线仿真方法,其特征在于,
步骤S2中对实际核能谱曲线图进行处理以获得模拟能谱曲线的方法包括:
步骤S21,对实际核能谱曲线图进行曲线识别并将能谱曲线数值化,以得到
能谱曲线各点的数值;
步骤S22,通过蒙特卡罗方法随机直接抽样这一组数值得到关于各核能级的
随机数,以模拟核衰变过程的随机性;
步骤S23,对所述随机数进行统计处理得到所述模拟能谱曲线。
4.根据权利要求3所述的模拟能谱曲线仿真方法,其特征在于,
所述步骤S21中对实际核能谱曲线图进行曲线识别并将能谱曲线数值化以
得到能谱曲线各点的数值的方法包括:
将实际核能各谱曲线图的图像经滤波去噪处理后再将实际能谱曲线图进行
显示,并根据显示的实际能谱曲线图临摹该核能谱曲线的各关键点获得能谱曲
线数据,以建立能谱曲线数据库;或
将实际核能各谱曲线图的图像经滤波、降噪预处理、曲线识别、曲线特征
提取,以及进行插值处理以完善及修复缺失能谱曲线的各点数据,以建立能谱
曲线数据库。
5.根据权利要求4所述的模拟能谱曲线仿真方法,其特征在于,
步骤S22中通过蒙特卡罗方法随机直接抽样这一组数值得到关于各核能级
的随机数以模拟核衰变过程的随机性的方法包括:
核信号时间统计特性模拟;以及核信号幅度统计特性模拟。
6.根据权利要求5所述的模拟能谱曲线仿真方法,其特征在于,
所述核信号时间统计特性模拟的方法包括:通过服从指数分布的随机数实
现核信号时间统计特性模拟,其中
指数分布的随机数由(0,1]均匀分布的随机数通过反函数法变换得到,且
(0,1]均匀分布随机数适于通过线性同余法求得。
7.根据权利要求6所述的模拟能谱曲线仿真方法,其特征在于,所述核信
号幅度统计特性模拟的方法包括:
通过对实际核能谱曲线识别并数字化得到各能级幅值及计数率,再通过蒙
特卡罗方法直接抽样并输出所述随机数;其中
对实际核能谱曲线的识别并数字化的过程包括:
步骤S221,对实际能谱曲线图进行滤波、降噪;
步骤S222,通过最大类间分割法求出阈值,并将能谱曲线图进行二值化处
理,再?#19978;?#32032;点扫描法提取出能谱曲线上各点的数值即坐标;
步骤S223,对能谱曲线进行修补及数值化。
8.根据权利要求7所述的模拟能谱曲线仿真方法,其特征在于,所述通过
蒙特卡罗方法直接抽样并输出所述随机数的方法,即
通过蒙特卡罗方法直接抽样能谱曲线及曲线上的各点数值,以获得一系列
随机数,从而以模拟核衰变过程的随机性。
9.根据权利要求8所述的模拟能谱曲线仿真方法,其特征在于,所述步骤
S221中对实际能谱曲线图进行滤波的方法,即对实际能谱曲线图进行维纳滤波
处理,以滤除能谱曲线图中的高斯噪声。
10.根据权利要求9所述的模拟能谱曲线仿真方法,其特征在于,
所述步骤S223中对能谱曲线进行修补及数值化包括:通过三次样条插值法
来填补在能谱曲线特征提取的过程中缺失的数据点,并通过坐标的比例扩伸以
获得能谱曲线图上各点的数值。

说明书

一种模拟能谱曲线仿真方法

?#38469;?#39046;域

本发明涉及一种核能领域,尤其涉及一种模拟能谱曲线仿真方法。

背景?#38469;?br />

核衰变过程在时间上是随机发生?#27169;?#34928;变过程释放射线(能量)也是随机
?#27169;?#20294;对其发生的时间间隔及能量值作统计分析,可得知核衰变过程有着以下
特性:在发生核衰变的时间间隔上近?#21697;?#20174;指数分布;核衰变过程对外释放的
能量(即能谱)近?#21697;?#20174;高斯分布。

基于核衰变过程存在以上特性,传统方式的仿核信号发生器是以服从不同
分布的随机数来模拟核信号特征,即以服从指数分布随机数模拟核信号在时间
间隔统计特性;以高斯分布随机数模拟核信号在幅度上的统计特性。但仅以服
从高斯分布的随机数来模拟核信号在幅度上的统计特性是不确?#26800;模?#27169;拟得到
的高期分布曲线与实际能谱曲线存在着较大的误差,不能准确地?#20174;?#26680;信号特
性;与此同时,各核素的能谱统计特性不尽相同,因此对不同种类核素的幅度
特性模拟则需要产生不同参数的高斯分布随机数与之匹配,这是不实?#23454;模?#22312;
际操作过程中难以实现。

鉴于传统意义上的仿核信号发生器存在的种种弊端,本文提出一种全新方
法用以解决以上问题。

发明内容

本发明的目的是提供一种模拟能谱曲线仿真方法,以实现对核能谱线进行
仿真。

为?#31169;?#20915;上述?#38469;?#38382;题,本发明提供了一种模拟能谱曲线仿真方法,其特
征在于,包括如下步骤:步骤S1,获取实际核能谱曲线图;以及步骤S2,对实
际核能谱曲线图进行处理以获得模拟能谱曲线。

进一步,所述模拟能谱曲线仿真方法还包括:步骤S3,通过反演比较模拟
能谱曲线与实际能谱曲线,以获得模拟能谱曲线与实际能谱曲线之间的误差。

进一步,步骤S2中对实际核能谱曲线图进行处理以获得模拟能谱曲线的方
法包括:步骤S21,对实际核能谱曲线图进行曲线识别并将能谱曲线数值化,以
得到能谱曲线各点的数值;步骤S22,通过蒙特卡罗方法随机直接抽样这一组数
值得到关于各核能级的随机数,以模拟核衰变过程的随机性;以及步骤S23,对
所述随机数进行统计处理得到所述模拟能谱曲线。

进一步,所述步骤S21中对实际核能谱曲线图进行曲线识别并将能谱曲线
数值化以得到能谱曲线各点的数值的方法包括:将实际核能各谱曲线图的图像
经滤波去噪处理后再将实际能谱曲线图进行显示,并根据显示的实际能谱曲线
图临摹该核能谱曲线的各关键点获得能谱曲线数据,以建立能谱曲线数据库;
或将实际核能各谱曲线图的图像经滤波、降噪预处理、曲线识别、曲线特征提
取,以及进行插值处理以完善及修复缺失能谱曲线的各点数据,以建立能谱曲
线数据库。

进一步,步骤S22中通过蒙特卡罗方法随机直接抽样这一组数值得到关于
各核能级的随机数以模拟核衰变过程的随机性的方法包括:核信号时间统计特
性模拟;以及核信号幅度统计特性模拟。

进一步,所述核信号时间统计特性模拟的方法包括:通过服从指数分布的
随机数实现核信号时间统计特性模拟,其中

指数分布的随机数由(0,1]均匀分布的随机数通过反函数法变换得到,且
(0,1]均匀分布随机数适于通过线性同余法求得。

进一步,所述核信号幅度统计特性模拟的方法包括:

通过对实际核能谱曲线识别并数字化得到各能级幅值及计数率,再通过蒙
特卡罗方法直接抽样并输出所述随机数;其中

对实际核能谱曲线的识别并数字化的过程包括:

步骤S221,对实际能谱曲线图进行滤波、降噪;

步骤S222,通过最大类间分割法求出阈值,并将能谱曲线图进行二值化处
理,再?#19978;?#32032;点扫描法提取出能谱曲线上各点的数值即坐标;

步骤S223,对能谱曲线进行修补及数值化。

进一步,所述通过蒙特卡罗方法直接抽样并输出所述随机数的方法,即

通过蒙特卡罗方法直接抽样能谱曲线及曲线上的各点数值,以获得一系列
随机数,从而以模拟核衰变过程的随机性。

进一步,所述步骤S221中对实际能谱曲线图进行滤波的方法,即对实际能
谱曲线图进行维纳滤波处理,以滤除能谱曲线图中的高斯噪声。

进一步,所述步骤S223中对能谱曲线进行修补及数值化的方法包括:通过
三次样条插值法来填补在能谱曲线特征提取的过程中缺失的数据点,并通过坐
标的比例扩伸以获得能谱曲线图上各点的数值。

本发明的有益效果是,本发明的模拟能谱曲线仿真方法通过对实际核能谱
曲线图进行曲线识别并将能谱曲线数值化,从而得到能谱曲线各点的数值(即
核能谱的能级和各能级的计数率),再通过蒙特卡罗方法随机直接抽样这一组数
值以得到关于各核能级的随机数,从而模拟核衰变过程的随机性,最终再对所
述随机数进行统计处理得到模拟能谱曲线,通过反演对比模拟能谱曲线与实际
能谱曲线以确定仿核信号发生器的可靠性与精确性。

附图说明

下面结合附图和实施例对本发明进一步说明。

图1是本发明的模拟能谱曲线仿真方法的原理框图;

图2是本发明的模拟能谱曲线仿真方法流程图;

图3是本发明的所述步骤S2中对实际核能谱曲线图进行处理以获得模拟能
谱曲线的方法流程图;

图4是本发明所涉及的生成n=10000个(0,1)均匀分布随机数分布图;

图5是本发明的指数分布随机数分布图;

图6是本发明的对以上指数分布随机数取值范围均匀划分1000个组矩并进
行统计的统计图;

图7是本发明的提取到的能谱曲线特征图;

图8是本发明的能谱曲线初步模拟效果图;

图9是本发明的三次样条插值后的效果图;

图10是本发明的最终得到该能谱曲线模拟效果图;

图11示出了模拟核信号随机发生过程的效果图;

图12示出了采用蒙特卡罗方法直接抽样最终效果图。

具体实施方式

现在结合附图对本发明作进一步详细的说明。这些附图均为简化的示意图,
仅以示意方式说明本发明的基本结构,因此其仅显示与本发明有关的构成。

如图1所示,本发明的模拟能谱曲线仿真方法通过对实际核能谱曲线图进
行曲线识别并将能谱曲线数值化,从而得到能谱曲线各点的数值(即核能谱的
能级和各能级的计数率),再通过蒙特卡罗方法随机直接抽样这一组数值以得到
关于各核能级的随机数,从而模拟核衰变过程的随机性,最终再对所述随机数
进行统计处理得到模拟能谱曲线,通过反演对比模拟能谱曲线与实际能谱曲线
以确定仿核信号发生器的可靠性与精确性。

本发明的具体实施方式如以下实施例所示。

如图2所示,本发明的一种模拟能谱曲线仿真方法,包括如下步骤:

步骤S1,获取实际核能谱曲线图;

步骤S2,对实际核能谱曲线图进行处理以获得模拟能谱曲线。

可选?#27169;?#25152;述模拟能谱曲线仿真方法还包括:

步骤S3,通过反演比较模拟能谱曲线与实际能谱曲线,以获得模拟能谱曲
线与实际能谱曲线之间的误差。

进一步,如图3所示,所述步骤S2中对实际核能谱曲线图进行处理以获得
模拟能谱曲线的方法包括:

步骤S21,对实际核能谱曲线图进行曲线识别并将能谱曲线数值化,以得到
能谱曲线各点的数值;步骤S22,通过蒙特卡罗方法随机直接抽样这一组数值得
到关于各核能级的随机数,以模拟核衰变过程的随机性;以及步骤S23,对所述
随机数进行统计处理得到所述模拟能谱曲线。

具体?#27169;?#25152;述步骤S21中对实际核能谱曲线图进行曲线识别并将能谱曲线
数值化以得到能谱曲线各点的数值的方法包括:

将实际核能各谱曲线图的图像经滤波去噪处理后再将实际能谱曲线图进行
显示,并根据显示的实际能谱曲线图临摹该核能谱曲线的各关键点获得能谱曲
线数据,以建立能谱曲线数据库;或将实际核能各谱曲线图的图像经滤波、降
噪预处理、曲线识别、曲线特征提取,以及进行插值处理以完善及修复缺失能
谱曲线的各点数据,以建立能谱曲线数据库。

其中,所述步骤S22中通过蒙特卡罗方法随机直接抽样这一组数值得到关
于各核能级的随机数以模拟核衰变过程的随机性的方法包括:核信号时间统计
特性模拟和核信号幅度统计特性模拟。

所述核信号时间统计特性模拟的方法包括:通过服从指数分布的随机数实
现核信号时间统计特性模拟,其中指数分布的随机数由(0,1]均匀分布的随机
数通过反函数法变换得到,且(0,1]均匀分布随机数适于通过线性同余法求得。

具体?#27169;?#36890;过线性同余法求得(0,1]均匀分布随机数的方法如下:

线性同余法的递推公式如下:

xi+1≡λxi+c(mod M) (1)


其中λ、c为常数。选取的初始x1称为种子,对随机数的生成质量有一定影
响,其值分别在1~216=65535之间选取。为了便于在计算机上使用,通常取

M=2S,其中S为计算机中二进制的最大可能?#34892;?#20301;数。

图4是取10000个(0,1]随机数分布情况

指数分布随机数的产生方法,即指数分布随机数可由反函数法实现,具体
过程如下:

设随机变量X的分布函数服从指数分布:

F(x)=1-e-ax,x≥0 (3)

其中,a是一个时间常数,e是自然底数。

由上式可以,F(x)∈[0,1),?#20197;?#23450;义域内单调递减,因此函数F(x)在0~+∞
必有反函数,求其反函数:


由于0<1-F(x)≤1,因此上式可以简化为


由式(5)可得知由符合(0,1]均匀分布的随机数抽样得到服从指数分布随
机数x。

取由以上单位均分布随机数通过反函数法产生的指数分布随机
数分布图如图5。对以上指数分布随机数取值范围均匀划分1000个组矩并进行
统计,最终统计图如图6所示。

所述核信号幅度统计特性模拟的方法包括:通过对实际核能谱曲线识别并
数字化得到各能级幅值及计数率,再通过蒙特卡罗方法直接抽样并输出所述随
机数;其中对实际核能谱曲线的识别并数字化的过程包括:

步骤S221,对实际能谱曲线图进行滤波、降噪;步骤S222,通过最大类间
分割法求出阈值,并将能谱曲线图进行二值化处理,再?#19978;?#32032;点扫描法提取出
能谱曲线上各点的数值即坐标;步骤S223,对能谱曲线进行修补及数值化。

具体?#27169;?#25152;述通过蒙特卡罗方法直接抽样并输出所述随机数的方法,即通
过蒙特卡罗方法直接抽样能谱曲线及曲线上的各点数值,以获得一系列随机的
随机数,从而以模拟核衰变过程的随机性。

所述步骤S221中对实际能谱曲线图进行滤波的方法,即对实际能谱曲线图
进行维纳滤波处理,以滤除能谱曲线图中的高斯噪声,从而尽可能的减小噪声
带来的干扰。

所述核信号幅度统计特性模拟的方法的具体实施过程如下:

所述步骤S221中对实际能谱曲线图进行滤波、降噪的具体实施步骤如下:

通过维纳滤波对实际能谱曲线图进行滤波、降噪处理,即所述维纳滤波器
为一种线性滤波器,并?#19968;?#26159;一种基于最小均方误差准则、对平稳过程的最优
估计器。

假设维纳滤波器输入信号为s(t),叠加噪声n(t)。输出信号x(t)通过滤波器g(t)
由下面卷积运算得到:

x(t)=g(t)*(s(t)+n(t)) (6)

对于估计出的信号x(t),期望等同于s(t)。

其误差为:e(t)=s(t+d)-x(t) (7)

方差为:e2(t)=s2(t+d)-2s(t+d)x(t)+x2(t) (8)

其中s(t+d)是所期望的滤波器输出。

将x(t)写成卷积积分,即


即可计算平方误差为:


其中Rs是s(t)的自相关函数,Rx是x(t)的自相关函数,Rxs是x(t)和s(t)的自相
关函数。维纳滤波的最终目的就是求最优的g(t),使得E(e2)最小。

所述步骤S222中通过最大类间分割法求出阈值,并将能谱曲线图进行二值
化处理,再?#19978;?#32032;点扫描法提取出能谱曲线上各点的数值即坐标;

最大类间方差法的具体算法过程如下:

设一幅图像的灰度值为1~m,其中灰度值为i的像素点个数为ni,N表示图
像像素点总数,那么灰度值为i出现的概率为:

令灰度值大于阈值k为c1组,即c1={1~k},灰度值大于阈值k的则为c2组,
c2={k+1~m},那么C1和C2出现的概率分别为:



计算得到C1和C2的灰度均值为:



其中,那么可得:

μr=ω1·μ1+ω2·μ2 (16)

由此可计算两组间的方差σ2为:σ2(k)=ω1(μ1-μr)2+ω2(μ2-μr)2 (17)

将式(16)代入式(17)可得:σ2(k)=ω2ω2(μ2-μ1)2

那么最佳阈值T*=Arg max{σ2(k)},0≤k<m-1 (18)

求得分割阈值T*=0.6353。

所述步骤S223中对能谱曲线进行修补及数值化的具体步骤如下:

实际核能谱曲线图经滤波去噪、二值化后为提取出核能谱曲线上各点的数
值即坐标,需提取能谱曲线特征,并将曲线数值化。具体过程如下:

首先,直线识别,即通过扫描核能谱曲线二值图的行和列,识别出核能谱
图中的直线;

其次,定点,由识别的直线判断出能谱曲线所在坐标系的横、纵坐标,并
定位原点,一般由上至下,从左往右扫描,识别出的第一条直线就为横、纵坐
标;

第三,能谱曲线特征提取。为减少图像中边框及坐标对曲线的影响,需将
边框滤除。滤除边框后再?#19978;?#32032;点扫描方法逐行或逐列扫描像素点为0的点(二
值图像中黑色为0,白色为1)。

最后,曲线数值化。提取出曲线后,通过计算扫描到的能谱曲线?#34892;?#28857;到
扫描原点的横行与纵行距离确定该像素点在图中的位置,最后通过乘以扩大坐
标的比例因子得到该像素点的坐标值。

最终提取能谱曲线特征效果如图7所示。

以及能谱曲线初步模拟效果如图8所示。

进一步,从图7和图8可看到,得到的模拟的能谱曲线图较之原能谱曲线
在某些点的数据发生缺失。为尽可能真实地?#20174;?#23454;际能谱曲线特性,需要对缺
失的数据进行填补修复。

具体?#27169;?#25152;述通过三次样条插值法来填补在能谱曲线特征提取的过程中缺
失的数据点,并通过坐标的比例扩伸以获得能谱曲线图上各点的数值,以实现
对缺失的数据?#34892;?#22320;填补与修复。

所述三次样条插值法来填补在能谱曲线特征提取的过程中缺失的数据点的
具体算法如下:

定义区间[a,b]上的分段函数S(x),若满足:

①S(x)在每个子区间[xi,xi+1]上是一个三次多项式函数;

②S(x)在整个区间[a,b]上有连续的二阶?#38469;?br />

则称S(x)为区间[a,b]上关于a=x0<x1<…<xn=b的一个三次样条函数。
从而三次样条插值问题为:给定函数g(x)的n+1个节点x0,x1,…,xn得函数
y0,y1,…,yn,求一个三次样条函数S(x),使其满足:

S(xj)=yj,j=0,1,...,n (19)

其中,函数S(x)称为g(x)的三次样条插值函数。

如果S(x)是f(x)的三次样样条插值函数,那么必须满足以下条件:

①插值条件,即

S(xj)=yj,j=0,1,...,n-1

②连续性条件,即


③一阶?#38469;?#36830;续条件,即


④二阶?#38469;?#36830;续条件,即


通过三次样条插值后的效果图,如图9所示,由其局部放大图可以看到,
三次样条插值后的数据点较为?#20132;?#36739;为逼近实际值。

实际核能谱曲线模拟效果,即对实际核能谱曲线图经上述图像处理后,最
终得到该能谱曲线模拟效果图如图10所示。

具体?#27169;?#36890;过蒙特卡罗方法直接抽样能谱曲线及曲线上的各点数值,以获
得一系列随机的随机数,从而以模拟核衰变过程的随机性。

图11示出了模拟核信号随机发生过程的效果图;

图12示出了采用蒙特卡罗方法直接抽样最终效果图(本图由实际能谱曲线
?#38469;?#20540;化后得到能级?#22270;?#25968;率这一数组后,随机抽样过程并统计得到的。本图
适于证明由蒙特卡罗抽样的合理性和准确性)。

即通过以上数字图像处理过程得到了模拟能谱曲线及曲线上各点的数值
(横坐标为道址Channel,纵坐标为计数率Count),再以蒙特卡罗方法直接抽
样这一组数据就能得到一系列随机的能级(能级例如但不限于通过多道分析器
量化得到?#27169;?#25152;述道址是指核衰变过程释放的能量经多道分析器量化后得到的)
随机数,从而以模拟核衰变过程的随机性。最后再对所述随机数进行统计,可
以得到模拟能谱曲线图,这样一方面可以验证系统的可靠性与精确性,另一方
面也可以反演于多道分析器,用以标定多道分析器的精确性。

采用蒙特卡罗方法,对某事件A发生的概率P(A)=p(未知),进行模拟计
算,具体计算方法:

(1)进行N次重复独立抽样试验,计算事件A发生次数为nA。

引入随机变量Xi,表?#38236;趇次试验中事件A发生次数,令


则有

(2)计算事件A在N次重复独立抽样试验中的发生频率fN,为


(3)当N充分大时,以概率fN作为概率P(A)=p的估计值为


(4)要求估计值为概率P(A)=p的无偏估计,即


以及直接抽样,即对核信号在时间和幅度上的特性是以两组服从不同分布
的随机数来模拟?#27169;?#32780;随机数是离散?#27169;?#19981;连续的。对于离散随机序列的抽样,
直接抽样法非常理想的。

离散型分布直接抽样法具体抽样过程如下:

设离散随机变量X的取值范围为Xi(i=0,1,2,3…),其概率分布为

P(X=Xi)=Pi(i=0,1,2,3……)。其中Pi≥0,

(1)产生(0,1)区间上均匀分布的随机数r;

(2)求得正整数n=0,1,2...,使得r满足

(3)抽取离散随机变量X的抽样值为X=Xn。而当0<r≤P0时X=X0;

(4)重复步骤(1)、(2)、(3)直到抽取n个样本值。

由于产生(0,1)均匀分布的随机数r若在区间的概率为


即事件出现的概?#23454;?#21516;于事件X=Xn发生的概率。

又因为随机数r服从(0,1)上的均匀分布,其概率密度函数为


其分布函数如下:


故产生的随机数r抽中样本值为X=xn的概率为


由此可知,由直接抽样法抽取到(X=Xn)的概?#23454;?#20215;于随机数Xn在随机数序
列X1,X2,...Xn出现的频率。

对于直接抽样法可靠性可由以下证明:

设X为离散型随机变量,其概率分布为Pi=P{X=Xi},其中i=1,2,…。X分别以
Pi取得Xi,则事件|X-E(X)|≥ε表示随机变量X取得所有满足
不等式|Xi-E(X)|≥ε的可能值Xi,则






由于事件X=Xi(i=0,1,2,…N)发生的概率为pi(0<pi<1),则X≠Xi的概
率则为1-pi,而每次X=Xi发生的概率是不变,且每次抽样结果同其它各次抽取
结果无关。因此X=Xi单个事件是一次贝努利试验,那么抽样n次,则为n重伯
努利试验。若令事件A(X=Xi)发生的次数为nA,即nA~B(n,p)。由于X1,X2,…,
Xn是n个相互独立且服从参数为p的0-1分布的随机变量,且

有D(Xi)=p(1-p),i=1,2,...n。给定?#25105;猞牛?则有


由(4.31)式可推得




因此可推得


化简得


即当抽取的次数n越大,抽样后事件A出现的次数与抽样总数的频率比越
接近于事件A发生的概率。

由直接抽样抽取随机数其误差为:

令因此


即是p的无偏估计,


即抽样n的次数越大,估计值越接近理论值p。

以上述依据本发明的理想实施例为启示,通过上述的说明内容,相关工作
人员完全可以在不偏离本项发明?#38469;?#24605;想的范围内,进行多样的变更以及修改。
本项发明的?#38469;?#24615;范围并不局限于说明书上的内容,必须要根据权利要求范围
来确定其?#38469;?#24615;范围。

关于本文
本文标题:一种模拟能谱曲线仿真方法.pdf
链接地址:http://www.pqiex.tw/p-5994910.html
关于我们 - 网站声明 - 网?#38236;?#22270; - 资源地图 - 友情链接 - 网站客服 - 联系我们

[email protected] 2017-2018 zhuanlichaxun.net网站版权所有
经营许可证编号:粤ICP备17046363号-1 
 


收起
展开
平码五不中公式规律 重庆欢乐生肖五星走势图 时时彩后二对应码 新闻联播 后三组胆是什么玩法 正版双喜大厅属于哪个公司 彩盈娱乐平台app 欢乐捕鱼人安卓版下载 聚星娱乐入口 老奶毛栗加盟赚钱嘛 北京pk10定位计划软件