王曉明，閆偉，徐澤林，陳榮，崔少康

中山市中醫(yī)院 設(shè)備科，廣東 中山528401

X射線管燈絲發(fā)射特性曲線的擬合算法

王曉明，閆偉，徐澤林，陳榮，崔少康

中山市中醫(yī)院 設(shè)備科，廣東 中山528401

X射線管燈絲發(fā)射特性曲線是X射線管的重要參數(shù)，在發(fā)生器里儲(chǔ)存當(dāng)前使用的X射線管燈絲發(fā)射特性曲線，X射線曝光準(zhǔn)備時(shí)就可以將燈絲電流加載到預(yù)定的值，以實(shí)現(xiàn)快速曝光。X射線管燈絲發(fā)射特性曲線的儲(chǔ)存方式有數(shù)據(jù)表格式和曲線方程式，前者獲取速度快但儲(chǔ)存容量大，而后者是相反的特點(diǎn)，但是以目前計(jì)算機(jī)的速度，使用曲線方程式來(lái)進(jìn)行存儲(chǔ)已經(jīng)可以滿足實(shí)際需求。本文采用VB軟件界面，對(duì)X射線管燈絲發(fā)射特性曲線以描點(diǎn)的方式采樣，再利用最小二乘法進(jìn)行曲線擬合，以獲得燈絲發(fā)射特性曲線，為發(fā)生器提供X射線管燈絲加熱電流參數(shù)值。

X射線管；燈絲發(fā)射特性曲線；最小二乘法；曲線擬合

0 引言

X射線管作為X射線的產(chǎn)生裝置，在X射線機(jī)和CT等設(shè)備的工作中起著關(guān)鍵的作用。X射線管產(chǎn)生的X射線質(zhì)量是由高壓發(fā)生器控制，曝光時(shí)為了能獲得用戶設(shè)定的曝光參數(shù)（管電壓和管電流），發(fā)生器在曝光之前就需要對(duì)燈絲按一定要求加熱，即在X射線管的燈絲發(fā)射特性曲線上獲得相應(yīng)的燈絲電流值并發(fā)送給燈絲控制回路進(jìn)行燈絲加熱。不同X射線管有不同的燈絲發(fā)射特性曲線，故發(fā)生器中必須存儲(chǔ)當(dāng)前使用的X射線管的燈絲發(fā)射特性曲線。存儲(chǔ)曲線的方法有數(shù)據(jù)表格式和曲線方程式，前者需要較大存儲(chǔ)區(qū)域，以不同的管電壓和管電流建立數(shù)據(jù)表格，計(jì)算機(jī)通過(guò)查詢的方式獲得數(shù)據(jù)，由于表格數(shù)據(jù)的不連續(xù)性，所以往往還需要通過(guò)插值算法獲得其它數(shù)據(jù)，但這種方法的優(yōu)點(diǎn)是速度快。而曲線方程式存儲(chǔ)量較小，只需要存貯方程的系數(shù)即可，最終通過(guò)函數(shù)值的計(jì)算獲得燈絲電流值，但這種方法的缺點(diǎn)是速度慢。不過(guò)在計(jì)算機(jī)高速發(fā)展下，這種函數(shù)求值的方法已經(jīng)能夠保證發(fā)生器實(shí)際控制的需要[1-2]。

本文設(shè)計(jì)了1種X射線管燈絲發(fā)射特性曲線的擬合算法軟件，能夠方便地對(duì)燈絲發(fā)射特性曲線的圖片通過(guò)描點(diǎn)采樣的方式獲得采樣點(diǎn)，利用最小二乘法原理進(jìn)行曲線擬合，從而得出燈絲發(fā)射特性曲線的近似方程，提供給發(fā)生器使用。

1 選取樣點(diǎn)

本系統(tǒng)是在VB6.0平臺(tái)上進(jìn)行開發(fā)的，燈絲發(fā)射特性曲線的樣點(diǎn)選取界面，見圖1。左側(cè)則顯示X射線管的燈絲發(fā)射特性曲線，右側(cè)為獲得的一系列采樣數(shù)據(jù)點(diǎn)坐標(biāo)，該曲線可由X射線管參數(shù)手冊(cè)掃描后導(dǎo)入。進(jìn)入采集界面后，根據(jù)界面頂部的提示進(jìn)行操作。首先設(shè)置其坐標(biāo)軸原點(diǎn)位置，即將鼠標(biāo)“十”字定位到坐標(biāo)的原點(diǎn)，然后按鼠標(biāo)左鍵確認(rèn)，界面右邊就會(huì)顯示坐標(biāo)點(diǎn)的位置。然后用同樣的方法確定縱坐標(biāo)（管電流坐標(biāo)）最大坐標(biāo)點(diǎn)位置和橫坐標(biāo)（燈絲電流）最大坐標(biāo)點(diǎn)位置。最后在圖中選定某一根燈絲發(fā)射特性曲線，用鼠標(biāo)“十”字定位到曲線上，一般選取10個(gè)采樣點(diǎn)左右，且具有均勻性和一定的代表性，每選一個(gè)點(diǎn)，按鼠標(biāo)左鍵進(jìn)行一次“樣點(diǎn)記錄”，并顯示在界面的右邊，最后將所選點(diǎn)的坐標(biāo)存入數(shù)據(jù)庫(kù)，供曲線擬合時(shí)調(diào)用[3]。

在曲線擬合之前，還需要輸入坐標(biāo)原點(diǎn)對(duì)應(yīng)的最小管電流值和最小燈絲電流值，以及縱坐標(biāo)和橫坐標(biāo)對(duì)應(yīng)的最大管電流值和燈絲加熱電流值。最后按曲線擬合按鈕就可以獲得所需要的曲線擬合方程的系數(shù)值和擬合效果。

圖1 燈絲發(fā)射特性曲線采樣點(diǎn)選取界面

2 曲線擬合

曲線擬合就是從一組已知的數(shù)據(jù)（xi,yi）（i = 1,2,3,……,m），求出自變量x和因變量y之間的近似函數(shù)關(guān)系：y=p(x)。從圖形上看，就是由幾個(gè)已知的點(diǎn)去求出曲線函數(shù)，在實(shí)際工程中應(yīng)用廣泛。不要求曲線p(x)經(jīng)過(guò)每1個(gè)數(shù)據(jù)點(diǎn)，只要能盡量反應(yīng)數(shù)據(jù)點(diǎn)的變化趨勢(shì)即可[4]。

2.1 最小二乘原理

曲線擬合的算法常用的是最小二乘法，最小二乘法的數(shù)據(jù)擬合中，最簡(jiǎn)單也是最常用的數(shù)學(xué)模型是多項(xiàng)式擬合。

最小二乘法由勒讓德和高斯分別在18世紀(jì)初期分別創(chuàng)立，其基本原理是根據(jù)給定的 [a,b]區(qū)間上的m個(gè)數(shù)據(jù)點(diǎn)（xi,yi），找到一個(gè)多項(xiàng)式 (1)，其中m＞n，使得其誤差平方和 [式(2)]最小，這時(shí)的Pn(x)稱為最小二乘法擬合多項(xiàng)式，對(duì)應(yīng)的曲線為擬合曲線[5]。

2.2 曲線擬合的具體實(shí)現(xiàn)

進(jìn)行曲線擬合時(shí)，首先確定多項(xiàng)式次數(shù)，根據(jù)已知數(shù)據(jù)點(diǎn)構(gòu)造多項(xiàng)式方程，從而解出多項(xiàng)式系數(shù)。系統(tǒng)中建立一個(gè)類模塊，用于擬合實(shí)現(xiàn)。主要步驟如下：

（1）首先確定擬合多項(xiàng)式的次數(shù)n，n必須小于數(shù)據(jù)點(diǎn)數(shù)m。

（2）計(jì)算并建立兩個(gè)數(shù)組SumX(0 to 2n) 和SumYX (0 to n)；其中即SumX( )=m，∑x2，∑x3，……，∑x2n；即SumYX( )=∑y，∑xy，

（3）構(gòu)造矩陣（二維數(shù)組）Matrix (0 to n+1)，其值為：

（4）根據(jù)矩陣求解方程組，獲得一組最終系數(shù)，部分代碼參考文獻(xiàn)[6]。

3 結(jié)果與討論

以球管燈絲發(fā)射特性曲線為例，管電流最大上限為1.2 A，燈絲加熱電流范圍是4.0～5.5 A，選定40 kV曲線進(jìn)行采樣，為確定曲線次數(shù)以及數(shù)據(jù)樣點(diǎn)選取的影響，本文選取了兩組數(shù)據(jù)，每組10個(gè)樣點(diǎn)，并將采樣點(diǎn)坐標(biāo)值轉(zhuǎn)換成實(shí)際的數(shù)據(jù)值，見表1。設(shè)計(jì)2～9次的多項(xiàng)式進(jìn)行擬合，擬合曲線見圖2。

圖2的擬合曲線表明：對(duì)于選取的兩組不同的樣點(diǎn)，在2～7次多項(xiàng)式擬合時(shí)，曲線都基本一致；在8次和9次多項(xiàng)式擬合時(shí)，顯示出較明顯的差異，但這個(gè)差異僅體現(xiàn)在邊界。對(duì)于多項(xiàng)式的次數(shù)，2次的拋物線擬合誤差較大，3次曲線也存在一定量的誤差，從4次開始，擬合誤差很小，曲線基本都能經(jīng)過(guò)樣點(diǎn)，到了8次和9次，曲線邊界開始有微微振蕩，與實(shí)際中曲線單調(diào)遞增有所不同。所以最好選取4～6次多項(xiàng)式曲線，可以滿足應(yīng)用要求[7]。

另外，多項(xiàng)式次數(shù)的選擇還與多項(xiàng)式計(jì)算的執(zhí)行時(shí)間有關(guān)，即多項(xiàng)式次數(shù)越高，執(zhí)行時(shí)間越長(zhǎng)。若燈絲發(fā)射特性的多項(xiàng)式計(jì)算是在單片機(jī)系統(tǒng)中完成的，一般建議采用低次多項(xiàng)式，這是由于單片機(jī)的時(shí)鐘和浮點(diǎn)計(jì)算能力的限制；若由PC機(jī)完成多項(xiàng)式計(jì)算，則運(yùn)算速度較快，可選用較高次數(shù)的多項(xiàng)式。本文采用Microchip的PIC18F系列8位單片機(jī)，系統(tǒng)晶振32 MHz，機(jī)器周期頻率為8 MHz，在選取4次多項(xiàng)式擬合曲線時(shí)，其運(yùn)算時(shí)間在1～2 ms之間，這一速度完全可以滿足實(shí)際使用的要求。若系統(tǒng)的運(yùn)算速度較快，則可選取5或6次曲線，精度更高些。

除了曲線擬合算法外，還可采用樣點(diǎn)間線性插值的方法。線性插值算法最突出的優(yōu)點(diǎn)就是運(yùn)算速度高，適合單片機(jī)系統(tǒng)，但為了達(dá)到一定的精度，選取的樣點(diǎn)數(shù)也較多，造成手工操作的繁瑣，同時(shí)還需占用較大的ROM存儲(chǔ)空間[8]。

4 結(jié)論

本文設(shè)計(jì)了1種X射線管燈絲發(fā)射特性曲線的擬合算法軟件，能夠?qū)崿F(xiàn)對(duì)所選樣點(diǎn)進(jìn)行曲線擬合，從而得出燈絲發(fā)射特性曲線的近似方程，便于曝光前給出合適的燈絲加熱值。根據(jù)實(shí)驗(yàn)數(shù)據(jù)，可以得出如下結(jié)論：

圖2 2～9次方程的擬合曲線

（1）不同的樣點(diǎn)選取對(duì)擬合結(jié)果基本無(wú)影響。

（2）基于最小二乘原理的4～6次多項(xiàng)式擬合曲線基本滿足所需燈絲加熱值的精度要求。

（3）多項(xiàng)式次數(shù)可根據(jù)硬件系統(tǒng)需求事先進(jìn)行選擇。

（4）本系統(tǒng)操作簡(jiǎn)便，自動(dòng)化程度高，減輕了調(diào)試人員的繁瑣工作。

本系統(tǒng)還存在進(jìn)一步優(yōu)化之處，如曲線圖片導(dǎo)入后的旋轉(zhuǎn)、放大功能等，以提高系統(tǒng)的人性化操作。

[1]王保華.生物醫(yī)學(xué)測(cè)量與儀器[M].第二版.上海:復(fù)旦大學(xué)出版社,2009:356-357.

[2]徐躍,梁碧玲.醫(yī)學(xué)影像設(shè)備學(xué)[M].第二版.北京:人民衛(wèi)生出版社,2009:27-29.

[3]Varah JM.Least squares data ftting with implicit functions[J].BIT Numer Mathemat,1996,36(4):842-854.

[4]施妙根,顧麗珍.科學(xué)和工程計(jì)算基礎(chǔ)[M].北京:清華大學(xué)出版社,1999:146-147.

[5]李蓓蕾.多次自適應(yīng)最小二乘曲線擬合方法及其應(yīng)用[D].武漢:長(zhǎng)江大學(xué),2014.

[6]Solev VN.The accuracy of the method of least squares[J].J Mathemat Sci,1999,93(3):443-446.

[7]陳嵐峰,楊靜瑜,崔崧,等.基于MATLAB 的最小二乘曲線擬合仿真研究[J].沈陽(yáng)師范大學(xué)學(xué)報(bào)(自然科學(xué)版),2014,32(1):75-79.

[8]蔡山,張浩,陳洪輝,等.基于最小二乘法的分段三次曲線擬合方法研究[J].科學(xué)技術(shù)與工程,2007,7(3):352-355.

Fitting Algorithm of the Filament Emission Characteristic Curve for X-Ray Tube

As one of the important parameters of the X-ray tube,the filament emission characteristic curve used currently shall be stored in the generator,and then the flament current can be loaded to a predetermined value when the X-ray exposure is prepared. There are two methods to store the curve in the generator:data table format and curve equation method. The former has rapid access to the data and more storage space;and the latter is just opposite. Under the high calculating speed of the modern computers,the generator has adequate speed to acquire the data from the curve equation. This paper uses the VB software interface,and sample the flament emission curve of X-ray tube by point-by-point manner,then applies linear least-squares method to perform curve ftting. Finally,the flament emission characteristic curve is acquired so as to provide the value of X-ray tube flament heating current for the generator.

X-ray tube;flament mission characteristic curve;the least squares;curve ftting

WANG Xiao-ming,YAN Wei,XU Ze-lin,CHEN Rong,CUI Shao-kang
Department of Equipment,Zhongshan Hospital of Traditional Chinese Medicine,Zhongshan Guangdong 528401,China

TH774

A

10.3969/j.issn.1674-1633.2016.09.012

1674-1633(2016)09-0050-04

2016-05-27

2016-06-21

作者郵箱：king580315@163.com