[發明專利]一種中子脈沖序列的頻譜分析方法

無效

• 申請號: 200810237026.2
• 申請日: 2008-12-31
• 公開（公告）號: CN101533045A
• 公開（公告）日: 2009-09-16
• 發明（設計）人: 魏彪;唐躍林;任勇;馮鵬;米德伶;潘英俊
• 申請（專利權）人: 重慶大學
• 主分類號: G01R23/16
• 分類號: G01R23/16;G01T3/00
• 代理公司: 重慶弘旭專利代理有限責任公司
• 代理人: 周韶紅

1.一種中子脈沖序列的頻譜分析方法，包括以下步驟：

(1)對中子脈沖數據進行采集，得到中子源以及被中子源激發的被探測體產生的中子探測計數的時間分布，所述時間分布的形式是由“0”或“1”組成的中子脈沖序列，所述中子脈沖序列按脈沖出現的位置進行存儲，沒有脈沖的不存儲；

(2)對采集后的數據包按一定的長度進行分塊，得到各塊中的中子脈沖原序列；各數據塊之間連續無間隔；

(3)對各數據塊的數據進行相關計算，得到各數據塊或塊間的相關序列；

(4)對各塊或塊間的相關序列進行進一步計算得到功率譜；

其特征是：當所述步驟(3)中的相關計算為自相關計算時；所述自相關計算使用快速偏移的方法，包括以下步驟：

a、將上述步驟(2)中所得的中子脈沖原序列x(i)，i＝0，1，…，n-1復制一份，得到偏移序列y(i)，y(i)＝x(i)，i＝0，1，…，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，…，n-1與序列x(i)，i＝0，…，n-1-k中數據“1”間的間距，得到此時的最小間距d，如果d＞0，轉至步驟d；如果k≥n-1或d＝0，則計算完畢，轉至下步；

f、未計算的其他相關值為0，預先初始化為0；

g、進行計算結果的處理，得到自相關序列?j＝0，1，......，n-1；

當所述步驟(3)中的相關計算為互相關計算時；所述互相關計算使用快速偏移的方法，包括以下步驟：

a、用其中兩個通道的序列x₁(i)與x₂(i)計算偏移0時x₁(i)與?x₂?(i)的相關值，即?

b、計算x₁(i)與x₂(i)中數據“1”間的間距，得到最小間距d(d＞0)；設偏移計數為k，k＝0；

c、計算x₂(i)偏移最小間距d時與x₁(i)的相關值，即?同時，累計偏移計數k＝k+d；

d、如果k＜n-1，則計算x₂(i)，i＝k，…，n-1與x₁(i)，i＝0，…，n-1-k中數據“1”間的間距，得到此時的最小間距d；如果d＞0，轉至步驟c；如果k≥n-1或d＝0，則計算完畢，轉至下步；

e、未計算的其他相關值為0，預先初始化為0；

f、改變移位方向，重復步驟c、d、e，計算出另一半結果從而得到完整的互相關序列?j＝-n+1，...，-1，0，1，...，n-1。

2.如權利要求1所述的中子脈沖序列的頻譜分析方法，其特征是：所述步驟(4)按照如下步驟進行：

a.當各數據塊中的數累計達到預先設定的閾值時，對各數據塊中計算所得的相關序列進行求和得到?

b.對求和結果?進行快速傅立葉變換進而得到功率譜。

3.如權利要求1所述的中子脈沖序列的頻譜分析方法，其特征是：對中子脈沖數據進行采集時，中子源以及被中子源激發的被探測體產生的中子由中子探測器探測，利用置于計算機中的高速數據采集卡對源通道信號和被探測的二個通道共三路信號進行高速數據采集，從而得到所述時間分布。

4.如權利要求3所述的中子脈沖序列的頻譜分析方法，其特征是：所述高速數據采集卡的采樣率為1GHz。

5.如權利要求2所述中子脈沖序列的頻譜分析方法，其特征是：所述數據塊的長度為1024個字節。?