運用數(shù)值微分方法研究核脈沖信號放大電路

周 偉, 周建斌, 肖永江, 周 靖, 趙 祥

(1.成都理工大學(xué) 核技術(shù)學(xué)院，四川 成都 610059; 2.萍鄉(xiāng)學(xué)院 機電系，江西 萍鄉(xiāng) 337055)

?

運用數(shù)值微分方法研究核脈沖信號放大電路

周 偉1, 周建斌1, 肖永江2, 周 靖1, 趙 祥1

(1.成都理工大學(xué) 核技術(shù)學(xué)院，四川 成都 610059; 2.萍鄉(xiāng)學(xué)院 機電系，江西 萍鄉(xiāng) 337055)

提出了一種嶄新的基于數(shù)值微分法的時域內(nèi)核脈沖信號放大電路的分析方法，彌補了頻域分析方法在脈沖信號的寬度、形狀和幅度特性分析上的不足。從實際信號放大電路出發(fā)，運用基爾霍夫電流定律，建立輸入信號和輸出信號之間的一階微分方程關(guān)系等式。方程解算過程中，利用數(shù)值微分算法求解，獲得輸入輸出之間的數(shù)學(xué)關(guān)系模型。通過反相放大電路和同相放大電路等實際電路的測試結(jié)果，表明時域內(nèi)運用數(shù)值微分算法研究核脈沖信號放大電路的方法是可行的。

信號放大電路； 數(shù)值微分法； 時域分析； VBA仿真

0 引 言

信號放大電路作為一種常用信號處理電路，在儀器儀表設(shè)計領(lǐng)域中得到了廣泛的應(yīng)用。放大電路常以運算放大器為基礎(chǔ)，通過將運算放大器、電阻和電容等進行組合，并引入反饋，組成如反相放大電路、同相放大電路等常用電路形式，以實現(xiàn)對核信號的放大、濾波以及成形等功能。核輻射測量中，放大器有兩種形式，一是前置放大器，二是主放大器[1-2]。前置放大器需要盡可能地靠近前級探測器，以減少探測器輸出端到放大器輸入端之間分布電容的影響，提高信噪比以及減少外界干擾[3]。主放大器用于將前置放大器的輸出信號作進一步放大和成形，同時保持探測器輸出的有用信息，如射線的能量信息和時間信息等，盡可能減少它們的失真[4-5]。

在信號放大電路的設(shè)計與分析上，傳統(tǒng)的方法有兩種。一種是運用純數(shù)學(xué)變換(拉氏變換、傅里葉變換等)的方式，將對信號的分析由時域轉(zhuǎn)換到頻域進行[5～7]，另一種是基于電路仿真軟件(Multisim，Pspice等)直接進行頻率響應(yīng)和基本性能分析[8～13]。上述兩種方法都忽略了信號放大電路在時域內(nèi)的特性，例如核脈沖信號的隨機性、寬度、幅度和形狀等。論文從核脈沖信號放大電路分析入手，基于基爾霍夫電流定律(Kirchhoff Current Law，KCL)構(gòu)建輸入與輸出之間一階微分方程關(guān)系等式，再運用數(shù)值微分法求解該等式[14-15]，最后得到時域內(nèi)信號放大電路的數(shù)學(xué)模型。同時，論文也對建立的數(shù)學(xué)模型進行了仿真驗證。

1 反相放大電路數(shù)學(xué)模型

圖1是一個典型的反相放大電路，常用于核脈沖信號的前置放大器和主放大器中。

圖1 反相放大電路

上述電路中，運用KCL以所標(biāo)位置為“節(jié)點”，可以建立輸入信號Uin和輸出信號Uout之間的關(guān)系等式，

(1)

式中，Un和Up分別是運算放大器的負端和正端電壓。根據(jù)運算放大器“虛短”特性可知，Un=Up=0，則式⑴可以變換為,

(2)

從形式上可知，式(2)是一個一階微分方程。由于核脈沖信號具有統(tǒng)計漲落的特性，對于式(2)，采用求解連續(xù)函數(shù)微分方程解的方法是不合適的。為了能夠在時域里面對式(2)進行求解，將輸入、輸出的模擬信號轉(zhuǎn)換成具有一定時間間隔(由ADC采樣頻率決定)的離散序列，即令Uin=x[n]，Uout=y[n]，dt=Δt。其中，Δt為ADC采樣間隔，則式⑵可以表示為：

(3)

(4)

由此可見，式⑷即為圖1反相放大電路的數(shù)學(xué)模型。輸出信號的脈沖寬度、形狀和幅度隨著參數(shù)k1和k2的變化而變化。實際應(yīng)用過程中，ADC選定后，Δt也隨之被確定，k1和k2只與電路中的電阻R1、Rf和電容C有關(guān)。此時，調(diào)整電阻和電容的值，可以得到不同的輸出信號，從而也實現(xiàn)時域內(nèi)對該放大電路的分析。1.2 數(shù)學(xué)模型的仿真

基于VBA平臺完成了反相放大電路數(shù)學(xué)模型的仿真測試。其中，VBA技術(shù)(Visual Basic for Application)是微軟公司為工程應(yīng)用而專門設(shè)計的一種可視化程序開發(fā)環(huán)境。將VBA技術(shù)和Excel應(yīng)用相結(jié)合，不僅能夠使用Excel提供的庫函數(shù)，訪問其特有的命令和結(jié)構(gòu)，使用統(tǒng)計分析和圖表繪制等功能；還可以創(chuàng)建自定義的用戶界面[16-17]。實際仿真過程中，選擇疊加噪聲的標(biāo)準(zhǔn)負指數(shù)信號作為輸入信號，可表示為：

瀝青性能優(yōu)劣直接影響到排水瀝青混凝土路用性能及耐久性。本研究選用A型和B型國產(chǎn)瀝青改性劑TAFPACK—Seper（簡稱TPS），改性基質(zhì)瀝青來制備高黏度改性瀝青，并參照《公路工程瀝青及瀝青混合料試驗規(guī)程》進行相應(yīng)試驗，主要技術(shù)性能檢測結(jié)果見表1。

(5)

其中：A表示標(biāo)準(zhǔn)負指數(shù)信號的幅度；τ表示時間常數(shù)；B(0.5-Rnd(1))表示疊加的噪聲。仿真時，為了觀察數(shù)學(xué)模型的信號放大和濾波等綜合處理性能，選擇的4組輸入信號在信號幅度和疊加噪聲上是不同的。其中，x1[n]為單一噪聲信號。

(6)

圖2是不同k1，k2參數(shù)反相放大電路數(shù)學(xué)模型的仿真波形。其中，取Δt=50 ns，對于輸出信號Output1，k1=k2=0.5，可以選擇R1=Rf=100 Ω，C=1 nF的電阻和電容組合。

圖2 不同參數(shù)時反相放大電路數(shù)學(xué)模型的仿真波形

從圖2中可以看出，在輸出信號的幅度上，根據(jù)電阻R1和Rf的組合進行控制。另外，在輸出信號的信噪比上，電容C起到低通濾波器的作用。綜合選擇濾波器的參數(shù)(Rf和C)，可以控制輸出信號的脈沖寬度和形狀。

2 同相放大電路數(shù)學(xué)模型

2.1 數(shù)學(xué)模型的建立

圖3是一種常見的同相放大電路。

圖3 同相放大電路

采用相同的分析方法，輸入信號Uin和輸出信號Uout之間的關(guān)系等式可以描述為：

(7)

其中，Un=Up=Uin，則式(7)可以轉(zhuǎn)換為：

(8)

令Uin=x[n]，Uout=y[n]，dt=Δt，式(8)經(jīng)過變換后，得到：

(8)

(9)

2.2 數(shù)學(xué)模型的仿真

圖4是不同k1，k2參數(shù)同相放大電路數(shù)學(xué)模型的仿真波形。其中，取Δt=50 ns，對于輸出信號Output1，k1=0.4，k2=0.2，可以選擇R2=125 Ω，Rf=250 Ω，C=1 nF的電阻和電容組合。

圖4 不同參數(shù)時反相放大電路數(shù)學(xué)模型的仿真波形

3 結(jié) 語

在時域內(nèi)引入數(shù)值微分法研究核脈沖信號放大電路的數(shù)學(xué)模型是一次成功的嘗試。時域內(nèi)實現(xiàn)信號的寬度、形狀和幅度等特性分析，這是對頻域分析方法是極大的補充。論文中，先后研究了反相放大電路和同相放大電路實際電路的數(shù)學(xué)模型，并完成了標(biāo)準(zhǔn)負指數(shù)信號(疊加噪聲信號)的仿真測試。測試結(jié)果顯示所研究數(shù)字成形模型的正確性?？梢灶A(yù)見，隨著實時高速采樣技術(shù)的快速發(fā)展，采樣的最小時間單元可以微分到Psec量級(部分可以達到fsec量級)，從而更加縮小離散系統(tǒng)與連續(xù)系統(tǒng)的時間界限。鑒于在時間域內(nèi)的模型簡單、概念清楚、易于實現(xiàn)、對比直觀等優(yōu)點，數(shù)值微分法對于完善離散系統(tǒng)的時域分析方法必將會有極大的促進。此外，研究的順利開展也得益于VBA軟件仿真平臺便捷的完成了數(shù)學(xué)模型的仿真驗證。最后，在前期工作基礎(chǔ)之上，下一步的研究包括2個方面：① 完成實際核脈沖信號測試，并進一步完善優(yōu)化數(shù)學(xué)模型。進而探討數(shù)學(xué)模型基于硬件載體(FPGA或DSP器件)的實現(xiàn)技術(shù)，以提升實際工程應(yīng)用中核儀器系統(tǒng)的綜合性能;② 繼續(xù)運用數(shù)值微分法研究濾波電路、成形電路等實際信號處理電路的時域內(nèi)數(shù)學(xué)模型，從而為研究電子電路仿真軟件的開發(fā)工作提出一些新的思路。

[1] 王經(jīng)瑾, 范天民, 錢永庚. 核電子學(xué)[M]. 北京：原子能出版社, 1983.

[2] 王芝英, 樓濱喬, 朱俊杰, 等. 核電子技術(shù)原理[M]. 北京：原子能出版社, 1989.

[3] 薛 楠. Protel DXP 2004 原理圖與PCB設(shè)計實用教程[M]. 北京：機械工業(yè)出版社, 2012.

[4] 湯 彬, 葛良全, 方 方, 等. 核輻射測量原理[M]. 哈爾濱：哈爾濱工程大學(xué)出版社, 2011.

[5] 薛 山. MATLAB基礎(chǔ)教程[M]. 北京：清華大學(xué)出版社, 2011.

[6] 劉 維. 精通Matlab與C/C++混合程序設(shè)計(第3版)[M]. 北京：北京航空航天大學(xué)出版社, 2012.

[7] 黃 喜, 張新亮, 董建績, 等. 半導(dǎo)體光放大器超快折射率變化動態(tài)特性的研究[J]. 物理學(xué)報，2009,58(5):3185-3192.

[8] 朱高中. 基于Multisim的高頻諧振功率放大器仿真實驗[J]. 實驗室研究與探索，2013,32(2):92-94.

[9] 鄭 鑫, 張曉潔. 虛擬實驗室立體化教學(xué)資源平臺建設(shè)簡[J]. 實驗室研究與探索，2014,33(8):109-111.

[10] 蒲永紅, 余 粟, 王維榮. Multisim輔助電工電子實驗教學(xué)的探討[J]. 實驗室研究與探索，2013,32(9): 174-177.

[11] 劉麗娜, 杜欽君. OrCAD PSpice在自動化專業(yè)綜合設(shè)計型實驗課程中的應(yīng)用探索[J]. 實驗技術(shù)與管理，2013,30(11):133-136.

[12] 王連英. 基于Multisim10的電子仿真試驗與設(shè)計[M]. 北京：北京郵電大學(xué)出版社, 2013.

[13] 周潤景. 基于PROTEUS的電路及單片機設(shè)計與仿真[M].第2版. 北京：北京航空航天學(xué)出版社, 2010.

[14] Zhou J B, Zhou W, Lei J R,etal. Study of time-domain digital pulse shaping algorithms for nuclear signals[J]. Nuclear Science & Techniques, 2012, 23(3):150-155.

[15] Zhou J B, Hong X, Wang R B,etal. Study of recursive model for pole-zero cancellation circuit[J]. Nuclear Science & Techniques, 2014(010403).

[16] 譚浩強. Visual Basic程序設(shè)計[M].第3版 北京：清華大學(xué)出版社, 2012.

[17] 譚浩強. QBasic語言教程[M].第2版. 北京：電子工業(yè)出版社, 2002.

Study of Nuclear Signal Amplifier Circuit Based on Numerical Differential Method

ZHOUWei1,ZHOUJian-bin1,XIAOYong-jiang2,ZHOUJing1,ZHAOXiang1

(1. Institute of Nuclear, Chengdu University of Technology, Chengdu 610059, China;2. School of Mechanical and Electronic Engineering, PingXiang Institute, Pingxiang 337055, China)

In the time-domain, a novel method of nuclear pulse signal amplifier circuit is presented based on numerical differential in this paper. The method can cover the many shortages of the analyzing methods in the frequency-domain, such as pulse width, shape and amplitude. According to the Kirchhoff's current law, an equation between input and output is established based on actual signal amplifier circuits. And after equation is solved, the mathematical model can also be acquired by using numerical differential. The simulating results of the inverting amplifying circuit and non-inverting amplifier circuit for nuclear pulse signal show that it is feasible for studying the amplifier circuits by using numerical differential method in the time-domain.

signal amplifier circuit; numerical differential method; time-domain analysis; simulation based on VBA

2014-08-11

周 偉(1979-)，男，江蘇南京人，博士，副教授，系主任，從事核信號獲取與數(shù)字處理的研究工作。

Tel.:13551805424,02884078773;E-mail:zhouwei@cdut.edu.cn。

TL 821

A

1006-7167(2015)03-0134-03