CN110598348A - 一种基于Matlab的核信号产生及处理的仿真方法 - Google Patents
一种基于Matlab的核信号产生及处理的仿真方法 Download PDFInfo
- Publication number
- CN110598348A CN110598348A CN201910898345.6A CN201910898345A CN110598348A CN 110598348 A CN110598348 A CN 110598348A CN 201910898345 A CN201910898345 A CN 201910898345A CN 110598348 A CN110598348 A CN 110598348A
- Authority
- CN
- China
- Prior art keywords
- signal
- gamma ray
- matlab
- nuclear
- energy
- 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.)
- Granted
Links
Landscapes
- Measurement Of Radiation (AREA)
Abstract
本发明公开了一种基于Matlab的核信号产生及处理的仿真方法，用Matlab基于蒙特卡洛方法粒子输运模型，构造了只需输入γ射线光子的能量和数量便能得到探测器内沉积能量的简易仿真方法，通过高斯展宽模拟统计涨落和噪声影响，用双指数函数模拟探测器物理特性拟合核信号波形，并用Simulink构建以放大、滤波为主的信号处理系统对核信号进行电学处理，实现了对整个核信号产生与处理全过程的仿真。
Description
技术领域
本发明涉及核电子学领域，具体涉及一种基于Matlab的核信号产生及处理的仿真方法。
背景技术
随着核技术的进一步应用，社会对核技术开发及应用专业人才的需求量越来越大。目前对于核技术应用方面的本科教育而言，专业课尤其是《核电子学》难度大、实验环境复杂、仪器价格昂贵，使得学习成本大，学习内容抽象，实验复杂度高，因此相关教育工作阻力很大。
软件仿真是解决实验风险和设备成本的一条理想途径，但现有核信号仿真软件例如蒙特卡洛模拟软件，使用门槛高，操作较为复杂。教学的发展需要一种能够模拟整个核信号产生与处理的仿真工具或建模、仿真方法。
发明内容
针对现有技术中的上述不足，本发明提供的一种基于Matlab的核信号产生及处理的仿真方法解决了现有仿真工具或仿真方法使用门槛高、操作较为复杂，且不能同时具备模拟核信号生成过程和处理过程功能的问题。
为了达到上述发明目的，本发明采用的技术方案为：一种基于Matlab的核信号产生及处理的仿真方法，包括以下步骤：
S1、设定入射γ射线光子的数量N以及每个入射γ射线光子的能量Eγi，其中i为在区间(0,N]内的整数；
S2、采用Matlab基于蒙特卡洛方法粒子输运模型，根据每个入射γ射线光子的能量Eγi模拟每个入射γ射线光子在碘化钠探测器内的反应，得到每个入射γ射线光子在探测器内沉积的能量Ei；
S3、将每个入射γ射线光子在探测器内沉积的能量Ei进行高斯展宽，得到高斯展宽的γ射线光子能量Edi；
S4、判断入射γ射线光子的数量N是否大于20，若是，则跳转至步骤S5，否则跳转至步骤S6；
S5、记录从i＝1～N累加过程中Ei的数值以及该数值出现的次数，并以Ei的各数值为横坐标，各数值出现的次数为纵坐标建立二维坐标下的曲线，得到能谱，并用Matlab进行窗体显示，完成仿真过程；
S6、将每个高斯展宽的γ射线光子能量Edi设定为每个γ射线光子各自的脉冲信号曲线Vi(t)的峰值，采用双指数函数拟合方法生成完整的脉冲信号曲线Vi(t)；
S7、将每个γ射线光子的脉冲信号曲线Vi(t)在时间分布上按照泊松分布叠加，得到核信号V(t)；
S8、使用MATLAB中的Simulink工具构建信号处理系统，对核信号V(t)进行信号处理，得到终端核信号Vout(t)。
进一步地：步骤S2中入射γ射线光子在碘化钠探测器内的反应包括：光电效应、康普顿效应和电子对效应。
进一步地：步骤S6中双指数函数拟合方法遵循下式：
其中τ为时间常数，其值为模拟的碘化钠探测器等效电阻R与等效电容C的乘积，λ为荧光衰减常数的倒数，t为时间。
进一步地：步骤S8包括以下步骤：
S81、使用Simulink工具的放大器控件构建前置放大子系统，对核信号V(t)进行前置放大，得到放大信号VG(t)；
S82、使用Simulink工具的极零相消控件，对放大信号VG(t)进行极零相消，得到第一极零相消信号VX1(t)；
S83、使用Simulink工具的跟随器控件构建跟随器子系统，对第一极零相消信号VX1(t)进行阻抗匹配，得到跟随信号VZ(t)；
S84、使用Simulink工具的极零相消控件，对跟随信号VZ(t)进行第二次极零相消，得到第二极零相消信号VX2(t)；
S85、使用Simulink工具的滤波放大控件，构建滤波放大子系统，对第二极零相消信号VX2(t)进行滤波放大，得到终端核信号Vout(t)。
进一步地：步骤S85中滤波放大控件包括：积分反馈滤波放大控件和S-K低通滤波放大控件。
进一步地：积分反馈滤波放大控件对第二极零相消信号VX2(t)进行滤波放大，得到滤波放大后的终端核信号Vout(t)。
进一步地：S-K低通滤波放大控件对第二极零相消信号VX2(t)进行低通滤波放大，得到低通滤波放大后的终端核信号Vout(t)。
本发明的有益效果为：用Matlab基于蒙特卡洛方法粒子输运模型，构造了只需输入γ射线光子的能量和数量便能得到探测器内沉积能量的简易仿真方法，通过高斯展宽模拟统计涨落和噪声影响，用双指数函数模拟探测器物理特性拟合核信号波形，并用Simulink构建以放大、滤波为主的信号处理系统对核信号进行电学处理，实现了对整个核信号产生与处理全过程的仿真，解决了现有仿真工具或仿真方法使用门槛高、操作较为复杂，且不能同时具备模拟核信号生成过程和处理过程功能的问题。
附图说明
图1为基于Matlab的核信号产生及处理的仿真方法流程示意图
具体实施方式
下面对本发明的具体实施方式进行描述，以便于本技术领域的技术人员理解本发明，但应该清楚，本发明不限于具体实施方式的范围，对本技术领域的普通技术人员来讲，只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内，这些变化是显而易见的，一切利用本发明构思的发明创造均在保护之列。
如图1所示，一种基于Matlab的核信号产生及处理的仿真方法的步骤为：
S1、设定入射γ射线光子的数量N以及每个入射γ射线光子的能量Eγi，其中i为在区间(0,N]内的整数；
S2、采用Matlab基于蒙特卡洛方法粒子输运模型，根据每个入射γ射线光子的能量Eγi模拟每个入射γ射线光子在碘化钠探测器内的反应，得到每个入射γ射线光子在探测器内沉积的能量Ei；
步骤S2中入射γ射线光子在碘化钠探测器内的反应包括：光电效应、康普顿效应和电子对效应。
S3、将每个入射γ射线光子在探测器内沉积的能量Ei进行高斯展宽，得到高斯展宽的γ射线光子能量Edi；
S4、判断入射γ射线光子的数量N是否大于20，若是，则跳转至步骤S5，否则跳转至步骤S6，即入射γ射线光子的数量N大于20，则输出能谱；入射γ射线光子的数量N小于等于20，则输出终端核信号Vout(t)；
S5、记录从i＝1～N累加过程中Ei的数值以及该数值出现的次数，并以Ei的各数值为横坐标，各数值出现的次数为纵坐标建立二维坐标下的曲线，即为能谱，并用Matlab进行窗体显示，完成仿真过程；
S6、将每个高斯展宽的γ射线光子能量Edi设定为每个γ射线光子各自的脉冲信号曲线Vi(t)的峰值，采用双指数函数拟合方法生成完整的脉冲信号曲线Vii(t)；
其中步骤S6中双指数函数拟合方法遵循下式：
其中τ为时间常数，其值为模拟的碘化钠探测器等效电阻R与等效电容C的乘积，λ为荧光衰减常数的倒数，t为时间。
S7、将每个γ射线光子的脉冲信号曲线Vi(t)在时间分布上按照泊松分布叠加，得到核信号V(t)；
S8、使用MATLAB中的Simulink工具构建信号处理系统，对核信号V(t)进行信号处理，得到终端核信号Vout(t)。
其中步骤S8包括以下步骤：
S81、使用Simulink工具的放大器控件构建前置放大子系统，对核信号V(t)进行前置放大，得到放大信号VG(t)；
S82、使用Simulink工具的极零相消控件，对放大信号VG(t)进行极零相消，得到第一极零相消信号VX1(t)；
S83、使用Simulink工具的跟随器控件构建跟随器子系统，对第一极零相消信号VX1(t)进行阻抗匹配，得到跟随信号VZ(t)；
S84、使用Simulink工具的极零相消控件，对跟随信号VZ(t)进行第二次极零相消，得到第二极零相消信号VX2(t)；
S85、使用Simulink工具的滤波放大控件，构建滤波放大子系统，对第二极零相消信号VX2(t)进行滤波放大，得到终端核信号Vout(t)。
步骤S85中滤波放大控件包括：积分反馈滤波放大控件和S-K低通滤波放大控件。
积分反馈滤波放大控件对第二极零相消信号VX2(t)进行滤波放大，得到滤波放大后的终端核信号Vout(t)。
S-K低通滤波放大控件对第二极零相消信号VX2(t)进行低通滤波放大，得到低通滤波放大后的终端核信号Vout(t)。
各级电子学控件可选择输出，可选择电子学控件对核信号V(t)进行处理，得到用户需要的终端核信号Vout(t)。
本发明用Matlab基于蒙特卡洛方法粒子输运模型，构造了只需输入γ射线光子的能量和数量便能得到探测器内沉积能量的简易仿真方法，通过高斯展宽模拟统计涨落和噪声影响，用双指数函数模拟探测器物理特性拟合核信号波形，并用Simulink构建以放大、滤波为主的信号处理系统对核信号进行电学处理，实现了对整个核信号产生与处理全过程的仿真，解决了现有仿真工具或仿真方法使用门槛高、操作较为复杂，且不能同时具备模拟核信号生成过程和处理过程功能的问题。
Claims (7)
1.一种基于Matlab的核信号产生及处理的仿真方法，其特征在于，包括以下步骤：
S1、设定入射γ射线光子的数量N以及每个入射γ射线光子的能量Eγi，其中i为在区间(0，N]内的整数；
S2、采用Matlab基于蒙特卡洛方法粒子输运模型，根据每个入射γ射线光子的能量Eγi模拟每个入射γ射线光子在碘化钠探测器内的反应，得到每个入射γ射线光子在探测器内沉积的能量Ei；
S3、将每个入射γ射线光子在探测器内沉积的能量Ei进行高斯展宽，得到高斯展宽的γ射线光子能量Edi；
S4、判断入射γ射线光子的数量N是否大于20，若是，则跳转至步骤S5，否则跳转至步骤S6；
S5、记录从i＝1～N累加过程中Ei的数值以及该数值出现的次数，并以Ei的各数值为横坐标，各数值出现的次数为纵坐标建立二维坐标下的曲线，得到能谱，并用Matlab进行窗体显示，完成仿真过程；
S6、将每个高斯展宽的γ射线光子能量Edi设定为每个γ射线光子各自的脉冲信号曲线Vi(t)的峰值，采用双指数函数拟合方法生成完整的脉冲信号曲线Vi(t)；
S7、将每个γ射线光子的脉冲信号曲线Vi(t)在时间分布上按照泊松分布叠加，得到核信号V(t)；
S8、使用MATLAB中的Simulink工具构建信号处理系统，对核信号V(t)进行信号处理，得到终端核信号Vout(t)。
2.根据权利要求1所述的基于Matlab的核信号产生及处理的仿真方法，其特征在于，所述步骤S2中入射γ射线光子在碘化钠探测器内的反应包括：光电效应、康普顿效应和电子对效应。
3.根据权利要求1所述的基于Matlab的核信号产生及处理的仿真方法，其特征在于，所述步骤S6中双指数函数拟合方法遵循下式：
其中τ为时间常数，其值为模拟的碘化钠探测器等效电阻R与等效电容C的乘积，λ为荧光衰减常数的倒数，t为时间。
4.根据权利要求1所述的基于Matlab的核信号产生及处理的仿真方法，其特征在于，所述步骤S8包括以下步骤：
S81、使用Simulink工具的放大器控件构建前置放大子系统，对核信号V(t)进行前置放大，得到放大信号VG(t)；
S82、使用Simulink工具的极零相消控件，对放大信号VG(t)进行极零相消，得到第一极零相消信号VX1(t)；
S83、使用Simulink工具的跟随器控件构建跟随器子系统，对第一极零相消信号VX1(t)进行阻抗匹配，得到跟随信号VZ(t)；
S84、使用Simulink工具的极零相消控件，对跟随信号VZ(t)进行第二次极零相消，得到第二极零相消信号VX2(t)；
S85、使用Simulink工具的滤波放大控件，构建滤波放大子系统，对第二极零相消信号VX2(t)进行滤波放大，得到终端核信号Vout(t)。
5.根据权利要求4所述的基于Matlab的核信号产生及处理的仿真方法，其特征在于，所述步骤S85中滤波放大控件包括：积分反馈滤波放大控件和S-K低通滤波放大控件。
6.根据权利要求5所述的基于Matlab的核信号产生及处理的仿真方法，其特征在于，所述积分反馈滤波放大控件对第二极零相消信号VX2(t)进行滤波放大，得到滤波放大后的终端核信号Vout(t)。
7.根据权利要求5所述的基于Matlab的核信号产生及处理的仿真方法，其特征在于，所述S-K低通滤波放大控件对第二极零相消信号VX2(t)进行低通滤波放大，得到低通滤波放大后的终端核信号Vout(t)。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910898345.6A CN110598348B (zh) | 2019-09-23 | 2019-09-23 | 一种基于Matlab的核信号产生及处理的仿真方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910898345.6A CN110598348B (zh) | 2019-09-23 | 2019-09-23 | 一种基于Matlab的核信号产生及处理的仿真方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110598348A true CN110598348A (zh) | 2019-12-20 |
| CN110598348B CN110598348B (zh) | 2022-08-02 |
Family
ID=68862334
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910898345.6A Active CN110598348B (zh) | 2019-09-23 | 2019-09-23 | 一种基于Matlab的核信号产生及处理的仿真方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110598348B (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112201104A (zh) * | 2020-09-07 | 2021-01-08 | 上海培云教育科技有限公司 | 一种数值计算的放射线测量仿真方法及系统 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104316954A (zh) * | 2014-09-28 | 2015-01-28 | 中国石油大学(华东) | 核物理实验仿真系统及其进行能谱测量实验和强度测量实验的方法 |
| CN104537186A (zh) * | 2015-01-09 | 2015-04-22 | 南方电网科学研究院有限责任公司 | 一种模拟特高压直流输电线路雷击暂态过程的方法 |
| CN106021723A (zh) * | 2016-05-20 | 2016-10-12 | 成都理工大学 | 蒙特卡罗虚拟探测器 |
| CN107193036A (zh) * | 2017-06-26 | 2017-09-22 | 成都理工大学 | 一种改进型核信号梯形脉冲成形方法和装置 |
| US20170357012A1 (en) * | 2015-01-05 | 2017-12-14 | Raycan Technology Co., Ltd. (Suzhou) | Method for digitalizing scintillation pulse |
| CN109884891A (zh) * | 2015-08-28 | 2019-06-14 | 易良碧 | 一种高精度智能仿核信号发生系统及其工作方法 |
-
2019
- 2019-09-23 CN CN201910898345.6A patent/CN110598348B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104316954A (zh) * | 2014-09-28 | 2015-01-28 | 中国石油大学(华东) | 核物理实验仿真系统及其进行能谱测量实验和强度测量实验的方法 |
| US20170357012A1 (en) * | 2015-01-05 | 2017-12-14 | Raycan Technology Co., Ltd. (Suzhou) | Method for digitalizing scintillation pulse |
| CN104537186A (zh) * | 2015-01-09 | 2015-04-22 | 南方电网科学研究院有限责任公司 | 一种模拟特高压直流输电线路雷击暂态过程的方法 |
| CN109884891A (zh) * | 2015-08-28 | 2019-06-14 | 易良碧 | 一种高精度智能仿核信号发生系统及其工作方法 |
| CN106021723A (zh) * | 2016-05-20 | 2016-10-12 | 成都理工大学 | 蒙特卡罗虚拟探测器 |
| CN107193036A (zh) * | 2017-06-26 | 2017-09-22 | 成都理工大学 | 一种改进型核信号梯形脉冲成形方法和装置 |
Non-Patent Citations (5)
| Title |
|---|
| MD FOIEZ AHMED 等: "A Monte Carlo Model of a Benchtop X-Ray Fluorescence Computed Tomography System and Its Application to Validate a Deconvolution-Based X-Ray Fluorescence Signal Extraction Method", 《IEEE TRANSACTIONS ON MEDICAL IMAGING》 * |
| THORSTEN STEZELBERGER 等: "Modeling the energy and timing digital signalprocessing for the Gamma Ray Energy Tracking In-Beam Array (GRETINA)", 《2017 IEEE NUCLEAR SCIENCE SYMPOSIUM AND MEDICAL IMAGING CONFERENCE (NSS/MIC)》 * |
| 周伟 等: "基于双参数模型核脉冲信号数字高斯成形技术", 《物探与化探》 * |
| 张软玉 等: "核信号数值仿真方法的研究及应用", 《核电子学与探测技术》 * |
| 覃骏 等: "基于MATLAB的核信号仿真与成形研究", 《电子质量》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112201104A (zh) * | 2020-09-07 | 2021-01-08 | 上海培云教育科技有限公司 | 一种数值计算的放射线测量仿真方法及系统 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110598348B (zh) | 2022-08-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP2715942B1 (en) | Spike domain circuit and modeling method | |
| JP2013226005A (ja) | 配電網電力フロー分析システム及び方法 | |
| Howlett et al. | Galaxy two-point covariance matrix estimation for next generation surveys | |
| CN107251090A (zh) | 尖峰域卷积电路 | |
| CN106021723B (zh) | 蒙特卡罗虚拟探测器 | |
| CN105824817A (zh) | 一种闪烁脉冲的数字化方法 | |
| CN110598348B (zh) | 一种基于Matlab的核信号产生及处理的仿真方法 | |
| Köthur et al. | Visual analytics for correlation‐based comparison of time series ensembles | |
| Pál et al. | Comments on the stochastic characteristics of fission chamber signals | |
| Saxena et al. | Digital pulse deconvolution with adaptive shaping for real-time high-resolution high-throughput gamma spectroscopy | |
| CN104242877B (zh) | 基于模拟cr‑rc电路的核脉冲信号数字高斯成形方法 | |
| CN106580359A (zh) | 探测器信号处理方法及装置 | |
| Regadio et al. | Unfolding using deep learning and its application on pulse height analysis and pile-up management | |
| CN104382611B (zh) | 光电倍增管输出能量的校正方法和装置 | |
| CN107894860A (zh) | 精准模拟触摸控制方法及装置 | |
| CN111553111B (zh) | 一种基于mcnp的数字仿核信号发生器 | |
| CN112462676A (zh) | 一种可模拟重叠核脉冲信号发生装置及其控制方法 | |
| CN105720950B (zh) | 一种核脉冲信号的产生方法及装置 | |
| CN116432703A (zh) | 基于复合神经网络模型的脉冲高度估计方法、系统及终端 | |
| Kafaee et al. | Pile-up correction by genetic algorithm and artificial neural network | |
| CN201622351U (zh) | 核信号随机特性模拟器 | |
| JPH06274662A (ja) | 非線形演算ユニットおよびそれによる情報処理装置 | |
| Bobin et al. | Sparse BSS From Poisson Measurements | |
| Lee et al. | Monte Carlo based time-domain Hspice noise simulation for CSA-CRRC circuit | |
| Jenderny et al. | Wave digital model of calcium‐imaging‐based neuronal activity of mice |
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 |