CN109271707B - 对核能谱线进行仿真的模拟能谱曲线仿真方法 - Google Patents
对核能谱线进行仿真的模拟能谱曲线仿真方法 Download PDFInfo
- Publication number
- CN109271707B CN109271707B CN201811085398.8A CN201811085398A CN109271707B CN 109271707 B CN109271707 B CN 109271707B CN 201811085398 A CN201811085398 A CN 201811085398A CN 109271707 B CN109271707 B CN 109271707B
- Authority
- CN
- China
- Prior art keywords
- spectrum curve
- energy spectrum
- nuclear
- actual
- curve
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000001228 spectrum Methods 0.000 title claims abstract description 180
- 238000000034 method Methods 0.000 title claims abstract description 113
- 238000004088 simulation Methods 0.000 title claims abstract description 37
- 230000008569 process Effects 0.000 claims abstract description 38
- 238000000342 Monte Carlo simulation Methods 0.000 claims abstract description 23
- 238000005070 sampling Methods 0.000 claims description 34
- 238000009826 distribution Methods 0.000 claims description 26
- 238000001914 filtration Methods 0.000 claims description 19
- 238000000605 extraction Methods 0.000 claims description 5
- 230000011218 segmentation Effects 0.000 claims description 5
- 230000009467 reduction Effects 0.000 claims description 3
- 230000008439 repair process Effects 0.000 claims description 3
- 238000007781 pre-processing Methods 0.000 claims 1
- 230000000694 effects Effects 0.000 description 12
- 238000010586 diagram Methods 0.000 description 8
- 230000003595 spectral effect Effects 0.000 description 5
- 238000005311 autocorrelation function Methods 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000005315 distribution function Methods 0.000 description 2
- 238000013139 quantization Methods 0.000 description 2
- 238000009827 uniform distribution Methods 0.000 description 2
- 241000764238 Isis Species 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 241000894007 species Species 0.000 description 1
- 230000005654 stationary process Effects 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Complex Calculations (AREA)
- Measurement Of Radiation (AREA)
Abstract
本发明涉及对核能谱线进行仿真的模拟能谱曲线仿真方法,包括如下步骤:步骤S1,获取实际核能谱曲线图;步骤S2,对实际核能谱曲线图进行处理以获得模拟能谱曲线;本发明的模拟能谱曲线仿真方法通过对实际核能谱曲线图进行曲线识别并将能谱曲线数值化,从而得到能谱曲线各点的数值(即核能谱的能级和各能级的计数率),再通过蒙特卡罗方法随机直接抽样这一组数值以得到关于各核能级的随机数,从而模拟核衰变过程的随机性,最终再对所述随机数进行统计处理得到模拟能谱曲线,通过反演对比模拟能谱曲线与实际能谱曲线以确定仿核信号发生器的可靠性与精确性。
Description
技术领域
本发明涉及一种核能领域,尤其涉及对核能谱线进行仿真的模拟能谱曲线仿真方法。
背景技术
核衰变过程在时间上是随机发生的,衰变过程释放射线(能量)也是随机的,但对其发生的时间间隔及能量值作统计分析,可得知核衰变过程有着以下特性:在发生核衰变的时间间隔上近似服从指数分布;核衰变过程对外释放的能量(即能谱)近似服从高斯分布。
基于核衰变过程存在以上特性,传统方式的仿核信号发生器是以服从不同分布的随机数来模拟核信号特征,即以服从指数分布随机数模拟核信号在时间间隔统计特性;以高斯分布随机数模拟核信号在幅度上的统计特性。但仅以服从高斯分布的随机数来模拟核信号在幅度上的统计特性是不确切的,模拟得到的高期分布曲线与实际能谱曲线存在着较大的误差,不能准确地反映核信号特性;与此同时,各核素的能谱统计特性不尽相同,因此对不同种类核素的幅度特性模拟则需要产生不同参数的高斯分布随机数与之匹配,这是不实际的,在际操作过程中难以实现。
鉴于传统意义上的仿核信号发生器存在的种种弊端,本文提出一种全新方法用以解决以上问题。
发明内容
本发明的目的是提供一种模拟能谱曲线仿真方法,以实现对核能谱线进行仿真。
为了解决上述技术问题,本发明提供了一种模拟能谱曲线仿真方法,其特征在于,包括如下步骤:步骤S1,获取实际核能谱曲线图;以及步骤S2,对实际核能谱曲线图进行处理以获得模拟能谱曲线。
进一步,所述模拟能谱曲线仿真方法还包括:步骤S3,通过反演比较模拟能谱曲线与实际能谱曲线,以获得模拟能谱曲线与实际能谱曲线之间的误差。
进一步,步骤S2中对实际核能谱曲线图进行处理以获得模拟能谱曲线的方法包括:步骤S21,对实际核能谱曲线图进行曲线识别并将能谱曲线数值化,以得到能谱曲线各点的数值;步骤S22,通过蒙特卡罗方法随机直接抽样这一组数值得到关于各核能级的随机数,以模拟核衰变过程的随机性;以及步骤S23,对所述随机数进行统计处理得到所述模拟能谱曲线。
进一步,所述步骤S21中对实际核能谱曲线图进行曲线识别并将能谱曲线数值化以得到能谱曲线各点的数值的方法包括:将实际核能各谱曲线图的图像经滤波去噪处理后再将实际能谱曲线图进行显示,并根据显示的实际能谱曲线图临摹该核能谱曲线的各关键点获得能谱曲线数据,以建立能谱曲线数据库;或将实际核能各谱曲线图的图像经滤波、降噪预处理、曲线识别、曲线特征提取,以及进行插值处理以完善及修复缺失能谱曲线的各点数据,以建立能谱曲线数据库。
进一步,步骤S22中通过蒙特卡罗方法随机直接抽样这一组数值得到关于各核能级的随机数以模拟核衰变过程的随机性的方法包括:核信号时间统计特性模拟;以及核信号幅度统计特性模拟。
进一步,所述核信号时间统计特性模拟的方法包括:通过服从指数分布的随机数实现核信号时间统计特性模拟,其中
指数分布的随机数由(0,1]均匀分布的随机数通过反函数法变换得到,且(0,1]均匀分布随机数适于通过线性同余法求得。
进一步,所述核信号幅度统计特性模拟的方法包括:
通过对实际核能谱曲线识别并数字化得到各能级幅值及计数率,再通过蒙特卡罗方法直接抽样并输出所述随机数;其中
对实际核能谱曲线的识别并数字化的过程包括:
步骤S221,对实际能谱曲线图进行滤波、降噪;
步骤S222,通过最大类间分割法求出阈值,并将能谱曲线图进行二值化处理,再由像素点扫描法提取出能谱曲线上各点的数值即坐标;
步骤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,对实际核能谱曲线图进行处理以获得模拟能谱曲线。
可选的,所述模拟能谱曲线仿真方法还包括:
步骤S3,通过反演比较模拟能谱曲线与实际能谱曲线,以获得模拟能谱曲线与实际能谱曲线之间的误差。
进一步,如图3所示,所述步骤S2中对实际核能谱曲线图进行处理以获得模拟能谱曲线的方法包括:
步骤S21,对实际核能谱曲线图进行曲线识别并将能谱曲线数值化,以得到能谱曲线各点的数值;步骤S22,通过蒙特卡罗方法随机直接抽样这一组数值得到关于各核能级的随机数,以模拟核衰变过程的随机性;以及步骤S23,对所述随机数进行统计处理得到所述模拟能谱曲线。
具体的,所述步骤S21中对实际核能谱曲线图进行曲线识别并将能谱曲线数值化以得到能谱曲线各点的数值的方法包括:
将实际核能各谱曲线图的图像经滤波去噪处理后再将实际能谱曲线图进行显示,并根据显示的实际能谱曲线图临摹该核能谱曲线的各关键点获得能谱曲线数据,以建立能谱曲线数据库;或将实际核能各谱曲线图的图像经滤波、降噪预处理、曲线识别、曲线特征提取,以及进行插值处理以完善及修复缺失能谱曲线的各点数据,以建立能谱曲线数据库。
其中,所述步骤S22中通过蒙特卡罗方法随机直接抽样这一组数值得到关于各核能级的随机数以模拟核衰变过程的随机性的方法包括:核信号时间统计特性模拟和核信号幅度统计特性模拟。
所述核信号时间统计特性模拟的方法包括:通过服从指数分布的随机数实现核信号时间统计特性模拟,其中指数分布的随机数由(0,1]均匀分布的随机数通过反函数法变换得到,且(0,1]均匀分布随机数适于通过线性同余法求得。
具体的,通过线性同余法求得(0,1]均匀分布随机数的方法如下:
线性同余法的递推公式如下:
xi+1≡λxi+c(mod M) (1)
其中λ、c为常数。选取的初始x1称为种子,对随机数的生成质量有一定影响,其值分别在1~216=65535之间选取。为了便于在计算机上使用,通常取
M=2S,其中S为计算机中二进制的最大可能有效位数。
图4是取10000个(0,1]随机数分布情况
指数分布随机数的产生方法,即指数分布随机数可由反函数法实现,具体过程如下:
设随机变量X的分布函数服从指数分布:
F(x)=1-e-ax,x≥0 (3)
其中,a是一个时间常数,e是自然底数。
由上式可以,F(x)∈[0,1),且在定义域内单调递减,因此函数F(x)在0~+∞必有反函数,求其反函数:
由于0<1-F(x)≤1,因此上式可以简化为
由式(5)可得知由符合(0,1]均匀分布的随机数抽样得到服从指数分布随机数x。
所述核信号幅度统计特性模拟的方法包括:通过对实际核能谱曲线识别并数字化得到各能级幅值及计数率,再通过蒙特卡罗方法直接抽样并输出所述随机数;其中对实际核能谱曲线的识别并数字化的过程包括:
步骤S221,对实际能谱曲线图进行滤波、降噪;步骤S222,通过最大类间分割法求出阈值,并将能谱曲线图进行二值化处理,再由像素点扫描法提取出能谱曲线上各点的数值即坐标;步骤S223,对能谱曲线进行修补及数值化。
具体的,所述通过蒙特卡罗方法直接抽样并输出所述随机数的方法,即通过蒙特卡罗方法直接抽样能谱曲线及曲线上的各点数值,以获得一系列随机的随机数,从而以模拟核衰变过程的随机性。
所述步骤S221中对实际能谱曲线图进行滤波的方法,即对实际能谱曲线图进行维纳滤波处理,以滤除能谱曲线图中的高斯噪声,从而尽可能的减小噪声带来的干扰。
所述核信号幅度统计特性模拟的方法的具体实施过程如下:
所述步骤S221中对实际能谱曲线图进行滤波、降噪的具体实施步骤如下:
通过维纳滤波对实际能谱曲线图进行滤波、降噪处理,即所述维纳滤波器为一种线性滤波器,并且还是一种基于最小均方误差准则、对平稳过程的最优估计器。
假设维纳滤波器输入信号为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中通过最大类间分割法求出阈值,并将能谱曲线图进行二值化处理,再由像素点扫描法提取出能谱曲线上各点的数值即坐标;
最大类间方差法的具体算法过程如下:
令灰度值大于阈值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)=ω1ω2(μ2-μ1)2
那么最佳阈值
T*=Arg max{σ2(k)},0≤k<m-1 (18)
求得分割阈值T*=0.6353。
所述步骤S223中对能谱曲线进行修补及数值化的具体步骤如下:
实际核能谱曲线图经滤波去噪、二值化后为提取出核能谱曲线上各点的数值即坐标,需提取能谱曲线特征,并将曲线数值化。具体过程如下:
首先,直线识别,即通过扫描核能谱曲线二值图的行和列,识别出核能谱图中的直线;
其次,定点,由识别的直线判断出能谱曲线所在坐标系的横、纵坐标,并定位原点,一般由上至下,从左往右扫描,识别出的第一条直线就为横、纵坐标;
第三,能谱曲线特征提取。为减少图像中边框及坐标对曲线的影响,需将边框滤除。滤除边框后再由像素点扫描方法逐行或逐列扫描像素点为0的点(二值图像中黑色为0,白色为1)。
最后,曲线数值化。提取出曲线后,通过计算扫描到的能谱曲线有效点到扫描原点的横行与纵行距离确定该像素点在图中的位置,最后通过乘以扩大坐标的比例因子得到该像素点的坐标值。
最终提取能谱曲线特征效果如图7所示。
以及能谱曲线初步模拟效果如图8所示。
进一步,从图7和图8可看到,得到的模拟的能谱曲线图较之原能谱曲线在某些点的数据发生缺失。为尽可能真实地反映实际能谱曲线特性,需要对缺失的数据进行填补修复。
具体的,所述通过三次样条插值法来填补在能谱曲线特征提取的过程中缺失的数据点,并通过坐标的比例扩伸以获得能谱曲线图上各点的数值,以实现对缺失的数据有效地填补与修复。
所述三次样条插值法来填补在能谱曲线特征提取的过程中缺失的数据点的具体算法如下:
定义区间[a,b]上的分段函数S(x),若满足:
①S(x)在每个子区间[xi,xi+1]上是一个三次多项式函数;
②S(x)在整个区间[a,b]上有连续的二阶导数。
则称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
②连续性条件,即
③一阶导数连续条件,即
④二阶导数连续条件,即
通过三次样条插值后的效果图,如图9所示,由其局部放大图可以看到,三次样条插值后的数据点较为平滑,较为逼近实际值。
实际核能谱曲线模拟效果,即对实际核能谱曲线图经上述图像处理后,最终得到该能谱曲线模拟效果图如图10所示。
具体的,通过蒙特卡罗方法直接抽样能谱曲线及曲线上的各点数值,以获得一系列随机的随机数,从而以模拟核衰变过程的随机性。
图11示出了模拟核信号随机发生过程的效果图;
图12示出了采用蒙特卡罗方法直接抽样最终效果图(本图由实际能谱曲线图数值化后得到能级和计数率这一数组后,随机抽样过程并统计得到的。本图适于证明由蒙特卡罗抽样的合理性和准确性)。
即通过以上数字图像处理过程得到了模拟能谱曲线及曲线上各点的数值(横坐标为道址Channel,纵坐标为计数率Count),再以蒙特卡罗方法直接抽样这一组数据就能得到一系列随机的能级(能级例如但不限于通过多道分析器量化得到的,所述道址是指核衰变过程释放的能量经多道分析器量化后得到的)随机数,从而以模拟核衰变过程的随机性。最后再对所述随机数进行统计,可以得到模拟能谱曲线图,这样一方面可以验证系统的可靠性与精确性,另一方面也可以反演于多道分析器,用以标定多道分析器的精确性。
采用蒙特卡罗方法,对某事件A发生的概率P(A)=p(未知),进行模拟计算,具体计算方法:
(1)进行N次重复独立抽样试验,计算事件A发生次数为nA。
引入随机变量Xi,表示第i次试验中事件A发生次数,令
(2)计算事件A在N次重复独立抽样试验中的发生频率fN,为
以及直接抽样,即对核信号在时间和幅度上的特性是以两组服从不同分布的随机数来模拟的,而随机数是离散的,不连续的。对于离散随机序列的抽样,直接抽样法非常理想的。
离散型分布直接抽样法具体抽样过程如下:
(1)产生(0,1)区间上均匀分布的随机数r;
(3)抽取离散随机变量X的抽样值为X=Xn。而当0<r≤P0时X=X0;
(4)重复步骤(1)、(2)、(3)直到抽取n个样本值。
又因为随机数r服从(0,1)上的均匀分布,其概率密度函数为
其分布函数如下:
故产生的随机数r抽中样本值为X=xn的概率为
由此可知,由直接抽样法抽取到(X=Xn)的概率等价于随机数Xn在随机数序列X1,X2,...Xn出现的频率。
对于直接抽样法可靠性可由以下证明:
由于事件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分布的随机变量,且
有
由(4.31)式可推得
而
因此可推得
化简得
即当抽取的次数n越大,抽样后事件A出现的次数与抽样总数的频率比越接近于事件A发生的概率。
由直接抽样抽取随机数其误差为:
以上述依据本发明的理想实施例为启示,通过上述的说明内容,相关工作人员完全可以在不偏离本项发明技术思想的范围内,进行多样的变更以及修改。本项发明的技术性范围并不局限于说明书上的内容,必须要根据权利要求范围来确定其技术性范围。
Claims (1)
1.一种模拟能谱曲线仿真方法,其特征在于,包括如下步骤:
步骤S1,获取实际核能谱曲线图;
步骤S2,对实际核能谱曲线图进行处理以获得模拟能谱曲线;
所述模拟能谱曲线仿真方法还包括:
步骤S3,通过反演比较模拟能谱曲线与实际能谱曲线,以获得模拟能谱曲线与实际能谱曲线之间的误差;
步骤S2中对实际核能谱曲线图进行处理以获得模拟能谱曲线的方法包括:
步骤S21,对实际核能谱曲线图进行曲线识别并将能谱曲线数值化,以得到能谱曲线各点的数值;
步骤S22,通过蒙特卡罗方法随机直接抽样这一组数值得到关于各核能级的随机数,以模拟核衰变过程的随机性;
步骤S23,对所述随机数进行统计处理得到所述模拟能谱曲线;
所述步骤S21中对实际核能谱曲线图进行曲线识别并将能谱曲线数值化以得到能谱曲线各点的数值的方法包括:
将实际核能各谱曲线图的图像经滤波、降噪预处理、曲线识别、曲线特征提取,以及进行插值处理以完善及修复缺失能谱曲线的各点数据,以建立能谱曲线数据库;
步骤S22中通过蒙特卡罗方法随机直接抽样这一组数值得到关于各核能级的随机数以模拟核衰变过程的随机性的方法包括:
核信号时间统计特性模拟;以及核信号幅度统计特性模拟;
所述核信号时间统计特性模拟的方法包括:通过服从指数分布的随机数实现核信号时间统计特性模拟,其中
指数分布的随机数由(0,1]均匀分布的随机数通过反函数法变换得到,且(0,1]均匀分布随机数适于通过线性同余法求得;
所述核信号幅度统计特性模拟的方法包括:通过对实际核能谱曲线识别并数字化得到各能级幅值及计数率,再通过蒙特卡罗方法直接抽样并输出所述随机数;其中
对实际核能谱曲线的识别并数字化的过程包括:
步骤S221,对实际能谱曲线图进行滤波、降噪;
步骤S222,通过最大类间分割法求出阈值,并将能谱曲线图进行二值化处理,再由像素点扫描法提取出能谱曲线上各点的数值即坐标;
步骤S223,对能谱曲线进行修补及数值化;
所述通过蒙特卡罗方法直接抽样并输出所述随机数的方法,即
通过蒙特卡罗方法直接抽样能谱曲线及曲线上的各点数值,以获得一系列随机数,从而以模拟核衰变过程的随机性;
所述步骤S221中对实际能谱曲线图进行滤波的方法,即对实际能谱曲线图进行维纳滤波处理,以滤除能谱曲线图中的高斯噪声;
所述步骤S223中对能谱曲线进行修补及数值化的方法包括:通过三次样条插值法来填补在能谱曲线特征提取的过程中缺失的数据点,并通过坐标的比例扩伸以获得能谱曲线图上各点的数值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811085398.8A CN109271707B (zh) | 2015-08-28 | 2015-08-28 | 对核能谱线进行仿真的模拟能谱曲线仿真方法 |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510539124.1A CN106483550B (zh) | 2015-08-28 | 2015-08-28 | 一种模拟能谱曲线仿真方法 |
CN201811085398.8A CN109271707B (zh) | 2015-08-28 | 2015-08-28 | 对核能谱线进行仿真的模拟能谱曲线仿真方法 |
Related Parent Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510539124.1A Division CN106483550B (zh) | 2015-08-28 | 2015-08-28 | 一种模拟能谱曲线仿真方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109271707A CN109271707A (zh) | 2019-01-25 |
CN109271707B true CN109271707B (zh) | 2022-04-19 |
Family
ID=58234736
Family Applications (4)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510539124.1A Active CN106483550B (zh) | 2015-08-28 | 2015-08-28 | 一种模拟能谱曲线仿真方法 |
CN201811085398.8A Active CN109271707B (zh) | 2015-08-28 | 2015-08-28 | 对核能谱线进行仿真的模拟能谱曲线仿真方法 |
CN201811086619.3A Active CN109239763B (zh) | 2015-08-28 | 2015-08-28 | 模拟核衰变过程的模拟能谱曲线仿真方法 |
CN201811085406.9A Active CN109086555B (zh) | 2015-08-28 | 2015-08-28 | 采用蒙特卡罗方法的模拟能谱曲线仿真方法 |
Family Applications Before (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510539124.1A Active CN106483550B (zh) | 2015-08-28 | 2015-08-28 | 一种模拟能谱曲线仿真方法 |
Family Applications After (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811086619.3A Active CN109239763B (zh) | 2015-08-28 | 2015-08-28 | 模拟核衰变过程的模拟能谱曲线仿真方法 |
CN201811085406.9A Active CN109086555B (zh) | 2015-08-28 | 2015-08-28 | 采用蒙特卡罗方法的模拟能谱曲线仿真方法 |
Country Status (1)
Country | Link |
---|---|
CN (4) | CN106483550B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108416355B (zh) * | 2018-03-09 | 2021-07-30 | 浙江大学 | 一种基于机器视觉的工业现场生产数据的采集方法 |
CN111276195B (zh) * | 2020-01-20 | 2023-09-19 | 深圳大榆树科技有限公司 | 利用最大类间方差法计算凝胶能谱图化合物数量的方法 |
CN111273337B (zh) * | 2020-02-27 | 2022-11-29 | 成都工贸职业技术学院 | 一种基于蒙特卡罗脉冲插值方法的核能谱处理方法 |
CN111553111B (zh) * | 2020-04-30 | 2023-03-28 | 成都航空职业技术学院 | 一种基于mcnp的数字仿核信号发生器 |
CN112462676B (zh) * | 2021-01-27 | 2021-05-07 | 泛华检测技术有限公司 | 一种可模拟重叠核脉冲信号发生装置及其控制方法 |
CN112462675B (zh) * | 2021-01-27 | 2021-05-07 | 泛华检测技术有限公司 | 一种级联型双指数核脉冲信号发生装置及其控制方法 |
CN114422041B (zh) * | 2021-12-23 | 2023-03-24 | 中国原子能科学研究院 | 一种核信号模拟方法、装置、终端及存储介质 |
CN114241846A (zh) * | 2021-12-24 | 2022-03-25 | 西安恒律模训科技发展有限公司 | 模拟γ放射性核素识别训练方法及其系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102298652A (zh) * | 2010-06-23 | 2011-12-28 | 成都理工大学 | 放射性测量中能谱漂移的模拟方法 |
CN102916683A (zh) * | 2012-10-18 | 2013-02-06 | 成都理工大学 | 一种可调参数核脉冲模拟方法 |
WO2014090297A1 (en) * | 2012-12-12 | 2014-06-19 | Qatar Foundation | System and method for the simulation of metabolic profiles |
CN104316954A (zh) * | 2014-09-28 | 2015-01-28 | 中国石油大学(华东) | 核物理实验仿真系统及其进行能谱测量实验和强度测量实验的方法 |
Family Cites Families (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4633088A (en) * | 1985-04-08 | 1986-12-30 | Packard Instrument Co., Inc. | Reverse sum quench measurement using a liquid scintillation counter |
AUPO427796A0 (en) * | 1996-12-20 | 1997-01-23 | University Of Queensland, The | Computer simulation of magnetic resonance spectra employing homotopy |
JP4309733B2 (ja) * | 2003-09-29 | 2009-08-05 | 株式会社東芝 | 燃焼計算方法および燃焼計算プログラム |
US7411188B2 (en) * | 2005-07-11 | 2008-08-12 | Revera Incorporated | Method and system for non-destructive distribution profiling of an element in a film |
CN201233446Y (zh) * | 2008-07-14 | 2009-05-06 | 成都理工大学 | 任意核能谱发生器 |
EP2347282A2 (en) * | 2008-10-10 | 2011-07-27 | Koninklijke Philips Electronics N.V. | Practical spect calibration method for quantification of nuclides with high-energy contributions |
CN201622351U (zh) * | 2009-11-18 | 2010-11-03 | 成都理工大学 | 核信号随机特性模拟器 |
CN102073060B (zh) * | 2009-11-24 | 2013-01-23 | 成都理工大学 | 核信号随机特性模拟方法 |
WO2011123837A2 (en) * | 2010-04-01 | 2011-10-06 | University Of Georgia Research Foundation, Inc. | Method and system using computer simulation for the quantitative analysis of glycan biosynthesis |
JP2012037305A (ja) * | 2010-08-05 | 2012-02-23 | Fujita Corp | 地盤の逐次非線形地震応答解析方法及び解析プログラムを記憶した記憶媒体 |
DE102011055075B4 (de) * | 2010-11-05 | 2013-04-18 | Stefan Brühl | Verfahren und Vorrichtung zur vorbereitenden Planung einer medizinischen Ionenstrahltherapie und Vorrichtung zur Durchführung einer medizinischen Ionenstrahltherapie |
CN102353972B (zh) * | 2011-07-01 | 2013-04-10 | 成都理工大学 | 多种模式的数字化多道谱仪 |
CN102928866B (zh) * | 2011-08-09 | 2015-05-20 | 中国辐射防护研究院 | 一种利用被动式探测器测量中子能谱和累积剂量的方法 |
CN103091701B (zh) * | 2011-10-28 | 2015-09-30 | 中国原子能科学研究院 | 多用途冷中子束流品质测量用飞行时间设备 |
CN102902958A (zh) * | 2012-09-19 | 2013-01-30 | 四川大学 | 基于图像识别的数字核信号处理和多参数核信息获取方法 |
CN103076622B (zh) * | 2012-10-31 | 2016-08-17 | 成都理工大学 | 一种稳谱用随机信号的产生方法 |
US10593436B2 (en) * | 2013-11-21 | 2020-03-17 | Terrapower, Llc | Method and system for generating a nuclear reactor core loading distribution |
CN103853929B (zh) * | 2014-03-17 | 2016-06-15 | 东华理工大学 | 一种基于蒙卡响应矩阵的低分辨率γ能谱反演解析系统及方法 |
-
2015
- 2015-08-28 CN CN201510539124.1A patent/CN106483550B/zh active Active
- 2015-08-28 CN CN201811085398.8A patent/CN109271707B/zh active Active
- 2015-08-28 CN CN201811086619.3A patent/CN109239763B/zh active Active
- 2015-08-28 CN CN201811085406.9A patent/CN109086555B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102298652A (zh) * | 2010-06-23 | 2011-12-28 | 成都理工大学 | 放射性测量中能谱漂移的模拟方法 |
CN102916683A (zh) * | 2012-10-18 | 2013-02-06 | 成都理工大学 | 一种可调参数核脉冲模拟方法 |
WO2014090297A1 (en) * | 2012-12-12 | 2014-06-19 | Qatar Foundation | System and method for the simulation of metabolic profiles |
CN104316954A (zh) * | 2014-09-28 | 2015-01-28 | 中国石油大学(华东) | 核物理实验仿真系统及其进行能谱测量实验和强度测量实验的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109086555B (zh) | 2022-04-19 |
CN106483550B (zh) | 2018-09-18 |
CN109271707A (zh) | 2019-01-25 |
CN109086555A (zh) | 2018-12-25 |
CN109239763A (zh) | 2019-01-18 |
CN106483550A (zh) | 2017-03-08 |
CN109239763B (zh) | 2022-04-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109271707B (zh) | 对核能谱线进行仿真的模拟能谱曲线仿真方法 | |
CN109884891B (zh) | 一种高精度智能仿核信号发生系统及其工作方法 | |
CN106483551B (zh) | 一种仿核信号发生器及其工作方法 | |
Venema et al. | A stochastic iterative amplitude adjusted Fourier transform algorithm with improved accuracy | |
CN109858408B (zh) | 一种基于自编码器的超声信号处理方法 | |
CN107480619B (zh) | 基于eemd和排列熵的探地雷达b扫描图像的降噪方法及系统 | |
CN106127688A (zh) | 一种超分辨率图像重建方法及其系统 | |
CN112415078A (zh) | 一种质谱数据谱图信号校准方法及装置 | |
US10396911B1 (en) | Noise analysis to reveal jitter and crosstalk's effect on signal integrity | |
CN107886113B (zh) | 一种基于卡方检验的电磁频谱噪声提取和滤波方法 | |
CN106483867B (zh) | 仿核信号发生系统用处理器的工作方法 | |
CN104111109B (zh) | 一种基于不同阶次统计量及支持向量机的机械振动状态识别方法 | |
CN116432703A (zh) | 基于复合神经网络模型的脉冲高度估计方法、系统及终端 | |
CN114118177B (zh) | 一种基于奇异谱分析的位平面降噪方法、系统及存储介质 | |
CN111091194A (zh) | 一种基于cavwnb_kl算法的操作系统识别方法 | |
Trigano et al. | Intensity estimation of spectroscopic signals with an improved sparse reconstruction algorithm | |
CN112508862B (zh) | 一种通过改进gan增强裂纹磁光图像的方法 | |
CN107622035B (zh) | 一种基于模拟退火的多项式相位信号自适应时频变换方法 | |
CN112733775A (zh) | 基于深度学习的高光谱影像分类方法 | |
CN111008356A (zh) | 一种基于WTSVD算法扣除背景的γ能谱集分析方法 | |
US9851201B2 (en) | Method and apparatus for using cepstrum and wavelet based algorithms for wall thickness measurement | |
Ramey et al. | Real-time RFI mitigation in radio astronomy | |
Grimaldi et al. | Multivariate linear parametric models applied to daily rainfall time series | |
CN116977227B (zh) | 基于局部结构变差的图像平滑方法及装置 | |
CN111160185B (zh) | 多尺度时间序列遥感影像趋势和断点检测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20221118 Address after: 201100 Room 201-202, Building E, No. 1618, Yishan Road, Minhang District, Shanghai Patentee after: ARRAY MICROELECTRONICS LTD. Address before: 213000 Hohai University, Xinbei District, Changzhou, Jiangsu Patentee before: Yi Liangbi |