CN101533045A - 一种中子脉冲序列的频谱分析方法 - Google Patents
一种中子脉冲序列的频谱分析方法 Download PDFInfo
- Publication number
- CN101533045A CN101533045A CN200810237026A CN200810237026A CN101533045A CN 101533045 A CN101533045 A CN 101533045A CN 200810237026 A CN200810237026 A CN 200810237026A CN 200810237026 A CN200810237026 A CN 200810237026A CN 101533045 A CN101533045 A CN 101533045A
- Authority
- CN
- China
- Prior art keywords
- neutron
- sequence
- data
- pulse
- neutron pulse
- 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
Images
Landscapes
- Measurement Of Radiation (AREA)
Abstract
本发明涉及一种中子脉冲序列的频谱分析方法,包括以下步骤:(1)对中子脉冲数据进行采集,得到中子源以及被中子源激发的被探测体产生的中子探测计数的时间分布;(2)对采集后的数据包按一定的长度进行分块,得到各块中的中子脉冲原序列;各数据块之间连续无间隔;(3)对各数据块的数据进行相关计算,得到各数据块或块间的相关序列;(4)对各块或块间的相关序列进行进一步计算得到功率谱;其特征是:所述步骤(3)中对各块数据进行相关计算时使用快速偏移的方法。这种频谱分析方法利用可实时计算的快速偏移相关算法以及功率谱计算步骤的优化,使得该方法能实现中子脉冲序列的实时频谱分析。
Description
技术领域
本发明涉及一种中子脉冲序列频谱分析方法。
背景技术
中子源驱动噪声分析测量法,也称功率谱密度分析法(power spectral density analysismeasurements),它的基本原理是:首先创立目标参数的计算公式,其次依托实验测得原函数,并快速计算相关函数和功率谱,最后将这些功率谱代入所创立的目标参数公式,由此便可得到中子辐射环境中相应的目标参数。通过构建实验测量系统进行高达1GHz的采样率进行高速数据采集,可以获得原函数,即中子脉冲序列。
对中子脉冲序列的频谱进行经典分析的方法包括自相关法和周期图法。自相关法是建立在维纳-幸钦定理的基础上的,其基本步骤是先对随机信号的N点数据作自相关估计然后对其作傅立叶变换,便得到功率谱估计。其中,相关函数的计算有两种方法,一是利用公式直接计算,二是通过FFT来实现快速计算。周期图法是先求N点数据的离散时间傅立叶变换,得到X(ejω),然后取其幅频特性的平方乘以1/N作为功率谱估计,称为周期图。但是无论是自相关法还是周期图法的计算速度都很慢,不能满足实时分析的要求。例如:当采样率为1GHz,计数率为3×106S-1(即每秒有3×106个脉冲出现),块(Block)长度为1024,采集1×109个Block时,要求频谱分析在1小时内完成,也即单个Block的计算分析必须在3.6μs内完成。而在现有PC机平台(Pentium D3.4G,单线程)下,对于单个Block(长度为1024)而言,采用目前使用较为广泛且被公认为最快速的快速傅立叶变换(FFT)程序包之一的FFTW的计算时间为13.5μs,而自相关计算时间则需要130μs。显然,传统的分析方法完全无法满足系统在3.6μs内实时地进行分析的要求。
发明内容
本发明的目的在于提供一种能够对中子脉冲序列进行实时分析的方法。
本发明的目的是通过以下技术方案实现的:一种中子脉冲序列的频谱分析方法,包括以下步骤:(1)对中子脉冲数据进行采集,得到中子源以及被中子源激发的被探测体产生的中子探测计数的时间分布,所述时间分布的形式是由“0”或“1”组成的中子脉冲序列,所述中子脉冲序列按脉冲出现的位置进行存储,没有脉冲的不存储;(2)对采集后的数据包按一定的长度(设为n)进行分块(共m块),得到各块中的中子脉冲原序列;各数据块之间连续无间隔;(3)对各数据块的数据进行相关计算,得到各数据块或块间的相关序列;(4)对各块或块间的相关序列进行进一步计算得到功率谱;其特征是:所述步骤(3)中对各块数据进行相关计算时使用快速偏移的方法。
若需要得到自相关功率谱,则上述步骤(3)中的相关计算为自相关,所述快速偏移的方法包括以下步骤:a、将上述步骤(2)中所得的中子脉冲原序列x(i),i=0,1,...n-1复制一份,得到偏移序列y(i),y(i)=x(i),i=0,1,L,n-1;b、计算偏移0时中子脉冲原序列x(i)与偏移序列y(i)的相关值,即r(0)=p,p为中子脉冲原序列x(i)中脉冲为“1”的个数;c、计算中子脉冲原序列x(i)与偏移序列y(i)中数据“1”间的间距,得到最小间距d(d>0);设偏移计数为k,k=0;d、计算偏移序列y(i)偏移最小间距d时与中子脉冲原序列x(i)的相关值,即 偏移计数k=k+d;e、如果k<n-1,则计算序列y(i),i=k,L,n-1与序列x(i),i=0,L,n-1-k中数据“1”间的间距,得到此时的最小间距d,如果d>0,转至步骤d;如果k≥n-1或d=0,则计算完毕,转至下步;f、未计算的其他相关值为0,可预先初始化为0;g、进行计算结果的处理,得到相关序列
对中子脉冲数据进行采集时,中子源以及被中子源激发的被探测体产生的中子由中子探测器探测,利用置于计算机中的高速数据采集卡对源通道信号和被探测的二个通道共三路信号进行高速数据采集,从而得到所述时间分布。当然,被探测的通道也可以为三个或四个。
若需要得到互相关功率频谱,则上述步骤(3)中的相关计算为互相关计算,所述快速偏移的方法包括以下步骤:a、用其中两个通道的序列x1(i)与x2(i)计算偏移0时x1(i)与x2(i)的相关值,即 b、计算x1(i)与x2(i)中数据“1”间的间距,得到最小间距d(d>0);设偏移计数为k,k=0;c、计算x2(i)偏移最小间距d时与x1(i)的相关值,即 同时,累计偏移计数k=k+d;d、如果k<n-1,则计算x2(i),i=k,L,n-1与x1(i),i=0,L,n-1-k中数据“1”间的间距,得到此时的最小间距d;如果d>0,转至步骤d;如果k≥n-1或d=0,则计算完毕,转至下步;e、未计算的其他相关值为0,可预先初始化为0;f、改变移位方向,重复步骤c、d、e,可算出另一半结果。由此,可以计算得到完整的互相关序列
所述高速数据采集卡的采样率可以为1GHz,当然也可以为500MHz或者其他频率。
所述数据块的长度可以为1024个字节,当然,块的长度也可以是256、512、2048等。
本发明的有益效果是:
(1)由于在进行相关计算时,利用了中子脉冲数据计数率低的特点,使用快速偏移的方法,省去了没有必要的偏移和相关计算,因此大大减少了相关计算的时间,从而能够对中子脉冲序列进行实时分析;当采样率为1GHz,计数率为3×106S-1(即每秒有3×106个脉冲出现),数据块长度为1024,采集1×109个数据块时,利用现有的PC平台(Pentium D3.4G),频谱分析能在1小时内完成,也即单个Block的计算分析可以在3.6μs内完成。与相应的英特尔数学内核库(Intel Math Kernel Library,MKL)相比较,快速偏移的自相关计算方法的计算效率提高了249倍,具体测试数据见表1。
表1
(2)由于对中子脉冲序列进行数据分块,对各块数据先做相关计算,然后对各块数据相关计算的结果累计求和,再对累计求和结果进行傅立叶变换得到功率谱;而现有的算法是先对各块数据进行相关计算,然后对各相关计算的结果进行傅立叶变化,再对各傅立叶变换的结果求和,因此本发明所采用的方法大大节约了计算的时间,从而能够对中子脉冲序列进行实时分析。
附图说明
图1为本发明实施例1优化功率谱计算示意图;
图2为基于PC平台的频谱分析结构示意图;
图3为中子脉冲序列的频谱测量与分析流程图;
图4为本发明实施例1测试中自相关计算结果图;
图5为本发明实施例1测试中自功率频谱计算结果图;
图6为本发明实施例2优化功率谱计算示意图;
图7为本发明实施例2测试中互相关计算结果图;
图8为本发明实施例2测试中互功率频谱计算结果图。
具体实施方式
下面通过实施例的方式进一步说明本发明,并不因此将本发明限制在所述的实施例范围之中。
实施例1:如图1至图3所示,一种中子脉冲序列的频谱分析方法,包括以下步骤:(1)对中子脉冲数据进行采集,得到中子源以及被中子源激发的被探测体产生的中子探测计数的时间分布,所述时间分布的形式是由“0”或“1”组成的中子脉冲序列,所述中子脉冲序列按脉冲出现的位置进行存储,没有脉冲的不存储;(2)对采集后的数据包按1024个字节的长度(设为n)进行分块(共m块),得到各块中的中子脉冲原序列x(i),i=0,1,...n-1;各数据块之间连续无间隔;(3)对各数据块的数据进行相关计算,得到各数据块或块间的相关序列;(4)对各块或块间的相关序列进行进一步计算得到功率谱。
所述步骤(3)中对各块数据进行的相关计算为自相关计算,使用快速偏移的方法,包括以下步骤:a、将上述步骤(2)中所得的中子脉冲原序列x(i),i=0,1,...n-1复制一份,得到偏移序列y(i),y(i)=x(i),i=0,1,L,n-1;b、计算偏移0时中子脉冲原序列x(i)与偏移序列y(i)的相关值,即r(0)=p,p为中子脉冲原序列x(i)中脉冲为“1”的个数;c、计算中子脉冲原序列x(i)与偏移序列y(i)中数据“1”间的间距,得到最小间距d(d>0);设偏移计数为k,k=0;d、计算偏移序列y(i)偏移最小间距d时与中子脉冲原序列x(i)的相关值,即 偏移计数k=k+d;e、如果k<n-1,则计算序列y(i),i=k,L,n-1与序列x(i),i=0,L,n-1-k中数据“1”间的间距,得到此时的最小间距d,如果d>0,转至步骤d;如果k≥n-1或d=0,则计算完毕,转至下步;f、未计算的其他相关值为0,可预先初始化为0;g、进行计算结果的处理,得到相关序列
对中子脉冲数据进行采集时,中子源以及被中子源激发的被探测体产生的中子由中子探测器探测,利用置于计算机中的高速数据采集卡对源通道信号和被探测的二个通道共三路信号进行高速数据采集,从而得到所述时间分布。所述高速数据采集卡的采样率为1GHz,计数率为3×106S-1。
上述计算完毕后,显示终端的显示如图4和图5所示。
实施例2:如图6所示,本实施例与实施例1所不同的是,则上述步骤(3)中的相关计算为互相关计算,所述快速偏移的方法包括以下步骤:a、用其中两个通道的序列x1(i)与x2(i)计算偏移0时x1(i)与x2(i)的相关值,即 b、计算x1(i)与x2(i)中数据“1”间的间距,得到最小间距d(d>0);设偏移计数为k,k=0;c、计算x2(i)偏移最小间距d时与x1(i)的相关值,即 同时,累计偏移计数k=k+d;d、如果k<n-1,则计算x2(i),i=k,L,n-1与x1(i),i=0,L,n-1-k中数据“1”间的间距,得到此时的最小间距d;如果d>0,转至步骤d;如果k≥n-1或d=0,则计算完毕,转至下步;e、未计算的其他相关值为0,可预先初始化为0;f、改变移位方向,重复步骤c、d、e,可算出另一半结果。由此,可以计算得到完整的互相关序列 所得功率谱为互功率谱。本实施例的其他内容与实施例1相同。
上述计算完毕后,显示终端的显示如如图7和图8所示。
Claims (7)
1、一种中子脉冲序列的频谱分析方法,包括以下步骤:
(1)对中子脉冲数据进行采集,得到中子源以及被中子源激发的被探测体产生的中子探测计数的时间分布,所述时间分布的形式是由“0”或“1”组成的中子脉冲序列,所述中子脉冲序列按脉冲出现的位置进行存储,没有脉冲的不存储;
(2)对采集后的数据包按一定的长度进行分块,得到各块中的中子脉冲原序列;各数据块之间连续无间隔;
(3)对各数据块的数据进行相关计算,得到各数据块或块间的相关序列;
(4)对各块或块间的相关序列进行进一步计算得到功率谱;
其特征是:所述步骤(3)中对各块数据进行相关计算时使用快速偏移的方法。
2、如权利要求1所述的中子脉冲序列的频谱分析方法,其特征是:步骤(2)中的相关计算为自相关计算;所述快速偏移的方法包括以下步骤:
a、将上述步骤(2)中所得的中子脉冲原序列x(i),i=0,1,...n-1复制一份,得到偏移序列y(i),y(i)=x(i),i=0,1,L,n-1;
b、计算偏移0时中子脉冲原序列x(i)与偏移序列y(i)的相关值,即r(0)=p,p为中子脉冲原序列x(i)中脉冲为“1”的个数;
c、计算中子脉冲原序列x(i)与偏移序列y(i)中数据“1”间的间距,得到最小间距d(d>0);设偏移计数为k,k=0;
d、计算偏移序列y(i)偏移最小间距d时与中子脉冲原序列x(i)的相关值,即 偏移计数k=k+d;
e、如果k<n-1,则计算序列y(i),i=k,L,n-1与序列x(i),i=0,L,n-1-k中数据“1”间的间距,得到此时的最小间距d,如果d>0,转至步骤d;如果k≥n-1或d=0,则计算完毕,转至下步;
f、未计算的其他相关值为0,可预先初始化为0;
g、进行计算结果的处理,得到自相关序列
3、如权利要求1所述的中子脉冲序列的频谱分析方法,其特征是:对中子脉冲数据进行采集时,中子源以及被中子源激发的被探测体产生的中子由中子探测器探测,利用置于计算机中的高速数据采集卡对源通道信号和被探测的二个通道共三路信号进行高速数据采集,从而得到所述时间分布;步骤(2)中的相关计算为互相关计算;所述快速偏移的方法包括以下步骤:
a、用其中两个通道的序列x1(i)与x2(i)计算偏移0时x1(i)与x2(i)的相关值,即
b、计算x1(i)与x2(i)中数据“1”间的间距,得到最小间距d(d>0);设偏移计数为k,k=0;
c、计算x2(i)偏移最小间距d时与x1(i)的相关值,即 同时,累计偏移计数k=k+d;
d、如果k<n-1,则计算x2(i),i=k,L,n-1与x1(i),i=0,L,n-1-k中数据“1”间的间距,得到此时的最小间距d;如果d>0,转至步骤d;如果k≥n-1或d=0,则计算完毕,转至下步;
e、未计算的其他相关值为0,可预先初始化为0;
f、改变移位方向,重复步骤c、d、e,计算出另一半结果从而得到完整的互相关序列
5、如权利要求1或2所述的中子脉冲序列的频谱分析方法,其特征是:对中子脉冲数据进行采集时,中子源以及被中子源激发的被探测体产生的中子由中子探测器探测,利用置于计算机中的高速数据采集卡对源通道信号和被探测的二个通道共三路信号进行高速数据采集,从而得到所述时间分布。
6、如权利要求4所述的中子脉冲序列的频谱分析方法,其特征是:所述高速数据采集卡的采样率为1GHz。
7、如权利要求4所述中子脉冲序列的频谱分析方法,其特征是:所述数据块的长度为1024个字节。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008102370262A CN101533045B (zh) | 2008-12-31 | 2008-12-31 | 一种中子脉冲序列的频谱分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008102370262A CN101533045B (zh) | 2008-12-31 | 2008-12-31 | 一种中子脉冲序列的频谱分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101533045A true CN101533045A (zh) | 2009-09-16 |
CN101533045B CN101533045B (zh) | 2011-04-13 |
Family
ID=41103772
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2008102370262A Expired - Fee Related CN101533045B (zh) | 2008-12-31 | 2008-12-31 | 一种中子脉冲序列的频谱分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101533045B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101713830B (zh) * | 2009-11-17 | 2011-11-09 | 重庆大学 | 一种瞬发中子衰减常数的综合在线测定方法 |
CN102788990A (zh) * | 2012-08-09 | 2012-11-21 | 中国科学院合肥物质科学研究院 | 一种基于量子点的中子能谱测量方法 |
CN109545403A (zh) * | 2018-11-16 | 2019-03-29 | 深圳中广核工程设计有限公司 | 一种反应堆堆内构件的健康管理方法及系统 |
CN109584137A (zh) * | 2018-10-24 | 2019-04-05 | 北京大学 | 一种脉冲序列格式转换方法及系统 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4837442A (en) * | 1988-03-10 | 1989-06-06 | The United States Of America As Represented By The United States Department Of Energy | Neutron range spectrometer |
DE10325147A1 (de) * | 2003-05-28 | 2004-12-16 | Friedrich-Schiller-Universität Jena | Verfahren zur Zeit-Frequenz-Analyse von Signalfolgen |
US7430479B1 (en) * | 2004-08-17 | 2008-09-30 | Science Applications International Corporation | System and method for analyzing content data |
-
2008
- 2008-12-31 CN CN2008102370262A patent/CN101533045B/zh not_active Expired - Fee Related
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101713830B (zh) * | 2009-11-17 | 2011-11-09 | 重庆大学 | 一种瞬发中子衰减常数的综合在线测定方法 |
CN102788990A (zh) * | 2012-08-09 | 2012-11-21 | 中国科学院合肥物质科学研究院 | 一种基于量子点的中子能谱测量方法 |
CN102788990B (zh) * | 2012-08-09 | 2014-04-16 | 中国科学院合肥物质科学研究院 | 一种基于量子点的中子能谱测量方法 |
CN109584137A (zh) * | 2018-10-24 | 2019-04-05 | 北京大学 | 一种脉冲序列格式转换方法及系统 |
CN109584137B (zh) * | 2018-10-24 | 2021-02-02 | 北京大学 | 一种脉冲序列格式转换方法及系统 |
CN109545403A (zh) * | 2018-11-16 | 2019-03-29 | 深圳中广核工程设计有限公司 | 一种反应堆堆内构件的健康管理方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN101533045B (zh) | 2011-04-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109061589B (zh) | 随机跳频雷达的目标运动参数估计方法 | |
EP0132807B1 (en) | Signal representation generator | |
CN101533045B (zh) | 一种中子脉冲序列的频谱分析方法 | |
CN101900761B (zh) | 一种高准确度非整周期采样谐波分析测量方法 | |
CN104462695A (zh) | 双耦合Duffing振子与变尺度相结合的微弱信号检测方法 | |
CN103293515A (zh) | 一种舰船线谱噪声源纵向分布特性测量方法 | |
CN101874744B (zh) | 用于长骨分析的超声导波参数测量方法 | |
CN103543026A (zh) | 基于振动传递率函数和支持向量机的结构损伤识别方法 | |
Sharkova et al. | Wavelet transform-based cross-correlation in the time-delay estimation applications | |
CN105092022B (zh) | 基于抹香鲸叫声的水下目标测距单脉冲筛选方法及装置 | |
CN101713830B (zh) | 一种瞬发中子衰减常数的综合在线测定方法 | |
CN102073065B (zh) | 一种消除地震数据单频干扰的方法 | |
CN104849551B (zh) | 一种谐相角分析方法 | |
CN114019236A (zh) | 一种电网谐波单通道混叠目标信号检测方法和装置 | |
Bogomolov et al. | Entropy-based technique for denoising of acoustic emission signals | |
Koo et al. | A possible application of the PD detection technique using electro-optic Pockels cell with nonlinear characteristic analysis on the PD signals | |
CN104950168A (zh) | 一种基于二次平均的低信噪比正弦信号高精度频率测量方法 | |
RU2371736C2 (ru) | Способ формирования текущего энергетического спектра выходного сигнала приемника, устройство для его осуществления и способ измерения дальности | |
RU2306592C1 (ru) | Способ активной идентификации линейных объектов управления | |
EP0425095A1 (en) | Method and apparatus for frequency measurement | |
CN115824394B (zh) | 一种噪声测量中的微弱信号相位对齐方法 | |
CN113805228B (zh) | 基于面波频散的地面微地震定位方法 | |
RU2660320C2 (ru) | Способ дискретизации и восстановления непрерывного сигнала | |
RU2399060C1 (ru) | Способ анализа многочастотных сигналов, содержащих скрытые периодичности | |
CN101256552B (zh) | 矢量信号处理方法及采用该方法的矢量信号处理系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20110413 Termination date: 20111231 |