基于iir濾波器的數(shù)字核脈沖信號高斯成形方法
【專利摘要】本發(fā)明公開了一種基于IIR濾波器的數(shù)字核脈沖信號高斯成形方法，包括以下步驟：根據(jù)模擬高斯成形系統(tǒng)的頻率響應的幅度譜，確定IIR數(shù)字濾波器的設計指標，即IIR數(shù)字濾波器的通帶截止頻率、阻帶截止頻率、通帶最大衰減、阻帶最小衰減；確定所需濾波器的類型（巴特沃斯濾波器、切比雪夫I型濾波器、切比雪夫II型濾波器、橢圓濾波器）；根據(jù)選擇的濾波器的類型以及設計指標，在MATLAB中計算出對應的IIR數(shù)字濾波器的系統(tǒng)函數(shù)的系數(shù)，得到IIR數(shù)字濾波器的系統(tǒng)函數(shù)；用IIR數(shù)字濾波器對數(shù)字核脈沖信號進行處理，實現(xiàn)數(shù)字核脈沖信號的高斯成形。該方法在濾除噪聲的同時，用較低的階數(shù)將數(shù)字核脈沖信號成形為準高斯信號，成形后的波形具有較好的準高斯特性。
【專利說明】基于MR濾波器的數(shù)字核脈沖信號高斯成形方法
【技術領域】
[0001]本發(fā)明涉及放射性測量中數(shù)字核脈沖信號的高斯成形，尤其涉及一種基于無限沖激響應(IIR)濾波器的數(shù)字核脈沖信號高斯成形方法。
【背景技術】
[0002]核能譜測量技術作為一種對物質(zhì)成分進行分析的重要方法，由于準確、靈敏、無損等特點，在眾多領域中得到廣泛應用。在傳統(tǒng)的核能譜測量儀器中，為了提高信噪比以及滿足后級電路對信號波形的需要，常用模擬Sallen-Key濾波器將模擬核脈沖信號濾波成形為準高斯波形。由于數(shù)字濾波成形可避免模擬濾波器固有的溫度漂移、噪聲、電壓漂移等問題，數(shù)字濾波成形技術的研究引起了科研工作者的密切關注?；贗IR濾波器實現(xiàn)對數(shù)字核脈沖信號的濾波成形處理，在濾除噪聲的同時，可用較低的階數(shù)將數(shù)字核脈沖信號濾波成形為準高斯信號。

【發(fā)明內(nèi)容】

[0003]本發(fā)明的目的在于公開一種基于IIR濾波器的數(shù)字核脈沖信號高斯成形方法，該方法克服了核脈沖信號模擬高斯成形的不足，在濾除噪聲的同時，用較低的階數(shù)將數(shù)字核脈沖信號濾波成形為準高斯信號，解決了核脈沖的數(shù)字高斯成形需求。
[0004]本發(fā)明是通過以下技術方案實現(xiàn)的，具體包括以下步驟:
根據(jù)模擬高斯成形系統(tǒng)的電路，得到系統(tǒng)電路輸入信號與輸出信號的微分方程，將微分方程在時域中進行求解，獲得模擬高斯成形系統(tǒng)的單位沖激響應，對單位沖激響應進行傅里葉變換，得到模擬高斯成形系統(tǒng)的頻率響應；
根據(jù)模擬高斯成形系統(tǒng)的頻率響應的幅度譜，確定IIR數(shù)字濾波器的設計指標，即IIR數(shù)字濾波器的通帶截止頻率、阻帶截止頻率、通帶最大衰減、阻帶最小衰減；
確定所需濾波器的類型(巴特沃斯濾波器、切比雪夫I型濾波器、切比雪夫II型濾波器、橢圓濾波器)；
根據(jù)選擇的濾波器的類型以及設計指標，在MATLAB中計算出對應的IIR數(shù)字濾波器的系統(tǒng)函數(shù)的系數(shù)，得到IIR數(shù)字濾波器的系統(tǒng)函數(shù)；
用IIR數(shù)字濾波器對數(shù)字核脈沖信號進行處理，實現(xiàn)數(shù)字核脈沖信號高斯成形。
[0005]與現(xiàn)有技術相比，本發(fā)明的一個或多個實施例可以具有如下優(yōu)點:
有效克服模擬高斯成形系統(tǒng)的不足，用較低的階數(shù)實現(xiàn)對數(shù)字核脈沖信號的濾波成形，成形后的波形具有較好的準高斯特性。
[0006]本發(fā)明的其它特征和優(yōu)點將在隨后的說明書中闡述，并且，部分地從說明書中變得顯而易見，或者通過實施本發(fā)明而了解。本發(fā)明的目的和其他優(yōu)點可通過在說明書、權利要求書以及附圖中所特別指出的結構來實現(xiàn)和獲得。
【專利附圖】

【附圖說明】[0007]附圖用來提供對本發(fā)明的進一步理解，并且構成說明書的一部分，與本發(fā)明的實施例共同用于解釋本發(fā)明，并不構成對本發(fā)明的限制。在附圖中:
圖1是基于IIR濾波器的數(shù)字核脈沖信號高斯成形方法流程圖；
圖2是模擬高斯成形系統(tǒng)的電路原理圖；
圖3是模擬高斯成形系統(tǒng)頻率響應的幅度譜；
圖4是實測核脈沖信號濾波成形后的波形圖。
【具體實施方式】
[0008]容易理解，根據(jù)本發(fā)明的技術方案，在不變更本發(fā)明的實質(zhì)精神下，本領域的一般技術人員可以提出本發(fā)明的多個結構方式和制作方法。因此以下【具體實施方式】以及附圖僅是本發(fā)明的技術方案的具體說明，而不應當視為本發(fā)明的全部或者視為本發(fā)明技術方案的限定或限制。
[0009]下面結合實施例及附圖對本發(fā)明作進一步詳細的描述。
[0010]如圖1所示，為基于IIR濾波器的數(shù)字核脈沖信號高斯成形方法，該方法包括以下步驟:
步驟10根據(jù)模擬高斯成形系統(tǒng)的電路，得到系統(tǒng)電路輸入信號與輸出信號的微分方程，在時域中推導出模擬高斯成形系統(tǒng)的沖激響應，在對其進行傅里葉變換，得到模擬高斯成形系統(tǒng)的頻率響應，即模擬Sallen-Key濾波器的頻率響應，包括如下步驟Al-Cl:
Al根據(jù)模擬Sallen-Key濾波器的電路原理圖(如圖2所示)，列出圖中輸入信號f(t)與輸出信號y(t)之間的數(shù)學關系為:
【權利要求】
1.基于IIR濾波器的數(shù)字核脈沖信號高斯成形方法，其特征在于，所述方法包括: (1)根據(jù)模擬高斯成形系統(tǒng)的電路，得到系統(tǒng)電路輸入信號與輸出信號的微分方程，在時域中求解出模擬高斯成形系統(tǒng)的單位沖激響應，對單位沖激響應進行傅里葉變換，推導出模擬高斯成形系統(tǒng)的頻率響應； (2)根據(jù)模擬高斯成形系統(tǒng)的頻率響應的幅度譜，確定IIR數(shù)字濾波器的設計指標，SPIIR數(shù)字濾波器的通帶截止頻率、阻帶截止頻率、通帶最大衰減、阻帶最小衰減； (3)確定所需濾波器的類型(巴特沃斯濾波器、切比雪夫I型濾波器、切比雪夫II型濾波器、橢圓濾波器)； (4)根據(jù)選擇的濾波器的類型以及設計指標，在MATLAB中計算出對應的IIR數(shù)字濾波器的系統(tǒng)函數(shù)的系數(shù)，得到IIR數(shù)字濾波器的系統(tǒng)函數(shù)； (5)用IIR數(shù)字濾波器對數(shù)字核脈沖信號進行處理，實現(xiàn)數(shù)字核脈沖信號高斯成形。
2.根據(jù)權利要求1所述的基于IIR濾波器的數(shù)字核脈沖信號高斯成形方法，其特征在于，所述(3)確定所需濾波器的類型，可以選擇巴特沃斯濾波器、切比雪夫I型濾波器、切比雪夫II型濾波器或橢圓濾波器，不同類型的濾波器，其幅頻響應具有不同的特點，巴特沃斯濾波器的幅頻響應在通帶內(nèi)具有最大平坦的特性，并且在通帶和阻帶內(nèi)，隨著頻率的增加而單調(diào)地下降；切比雪夫I型濾波器的幅頻響應在通帶內(nèi)等波紋變化，阻帶內(nèi)單調(diào)下降；切比雪夫II型濾波器的幅頻響應在通帶內(nèi)單調(diào)下降，阻帶內(nèi)有等波紋變化；橢圓濾波器的幅頻響應在通帶和阻帶內(nèi)均為等波紋的，選擇不同類型的濾波器，推導出的IIR數(shù)字濾波器的系統(tǒng)函數(shù)不同，但都可實現(xiàn)對數(shù)字核脈沖信號的濾波成形。
【文檔編號】G01T1/36GK103777228SQ201410066240
【公開日】2014年5月7日 申請日期:2014年2月26日 優(yōu)先權日:2014年2月26日
【發(fā)明者】葛青, 葛良全, 羅耀耀 申請人:成都理工大學