李新獻+李偉勤+施岱松

摘要：矩量法是將連續(xù)方程離散為代數(shù)方程組的方法，此法對于求解微分方程和積分方程均適用，本文以半波振子天線為例，詳細推導(dǎo)了半波振子天線的Pocklington積分方程，利用分域三角基和伽遼金方法的矩量法求解Pocklington積分方程，得到了半波振子天線上的電流分布，分析了對應(yīng)的電場和磁場的分布情況，并簡單介紹了半波陣子天線在電子藥丸上實現(xiàn)無線通信及能量傳輸?shù)膽?yīng)用。

關(guān)鍵詞：矩量法；半波陣子天線；分域三角基；伽遼金方法

中圖分類號：TM144

文獻標(biāo)識碼：A

DOI：10.3969/j.issn.1003-6970.2015.06.026

本文著錄格式：李新獻，李偉勤，施岱松，矩量法求解半波偶極子天線上的電流分布陰，軟件，2015，36（6）：139-145

CalculationoftheSurfaceCurrentDistributionoftheDipoleAntennaBasedonMoM

LIXin-xianl，LIWei-qin2，SHIDai-song2

[Abstract]：Suitableforsolvingthedifferentialequationsandintegralequations，MethodofMomentsisthediscretecontinuityequationintoalgebraicequations.Takingthehalf-wavedipoleantennaasanexample，detailedhalf-wavedipoleantennaforthePocklingtonintegralequationisdeduced，andthemomentmethodtosolvePocklingtonintegralequationsbyusingthetrianglebasisandgammaLiao.jinmethodandthecurrentdistributionhalf-wavedipoleantennawereobtained.Thedistributionofthecorrespondingelectricfieldandmagneticfieldareanalyzedandtheapplicationofhalfwavedipolearrayantennaforwirelesscommunicationandenergytransmissioninelectronicpillareintroduced.

[Keywords]：Methodofmoment；Half-wavedipoleantenna；Piecewisetrianglefunction；Galerkinmethod

0引言

天線是用來進行電磁波輻射和接收裝置，主要用于通信、探測和能量輸送，結(jié)構(gòu)上主要有由導(dǎo)線或金屬棒構(gòu)成的用于長波、短波、超短波線天線和由金屬面或介質(zhì)面構(gòu)成的面天線，主要是用于微波波段[8]；對天線的求解方法有三種：一種為精確求解的解析方法，是通過求解滿足邊界條件下的麥克斯韋方程組實現(xiàn)，第二種是以逐步逼近法、微擾法、變分法和迭代法為代表的近似法（準(zhǔn)解析法），其求解結(jié)果一般表示為級數(shù)解，多用于第一種方法不能解決的問題；第三種方法是數(shù)值方法，其核心思想是將微分方程化為差分方程，或?qū)⒎e分方程化為有限求和，從而建立代數(shù)方程組，也可以將微分方程或積分方程用矩量求解。

半波偶極子天線是一種典型的低增益的輻射組件，其感應(yīng)電流分布在零階近似的情況下數(shù)學(xué)處理非常方便，其傳輸函數(shù)、頻域和時域的輻射場均可以用解析函數(shù)表達出來[7]。因此，通過分析其天線輻射機理和特性，對弄清楚寬帶和超寬帶天線天線如何工作具有重要意義。而基于數(shù)值方法的矩量法不僅可以用來求解積分方程，而且可以用來求解其他算子方程，是一種誤差最小的方法，具有廣泛的應(yīng)用。矩量法是一種采用基函數(shù)和檢驗函數(shù)離散化積分方程的數(shù)值方法，它的基本思想是將一個泛函方程化為矩陣方程，然后求解該矩陣方程。對于天線問題可以建立描述天線表面感應(yīng)電流的積分方程，求出該電流即可進一步得到天線的輻射特性。本文應(yīng)用矩量法求解偶極子天線電流的Pocklington積分方程，進而由電流分布求出該天線的其他特性。

1矩量法基本原理

根據(jù)線性空間理論，N個線性方程的聯(lián)立方程組、微分方程、差分方程及積分方程均屬于希爾伯特空間中的算子方程，它們可化作矩陣方程予以求解，在求解過程中需計算廣義矩量，故此法稱為矩量法。

如果非齊次方程為：

L（f）=g

（1）

式中L是線性算子，g為已知激勵函數(shù)，f為未知響應(yīng)函數(shù)，算子L的定義域為算子作用于其上的函數(shù)f的集合，算子L的值域為算子在其定義域上運算而得到的函數(shù)g。L取不同形式，人們便可描繪不同的電磁場工程問題，如：

1.1矩量法（MoM）的基本步驟

1）展開未知函數(shù)廠為有限個線性無關(guān)的已知簡單函數(shù)fn之和。式中an是系數(shù)，fn被稱為N線性無關(guān)的展開函數(shù)或基函數(shù)。N個未知數(shù)al.a2，a3aN，此時算子L作用域已知函數(shù)fn上，待求系數(shù)a移到算子L之外，使求解方便。

2）作內(nèi)積

在L的值域內(nèi)選一組線性無關(guān)的函數(shù)Wm（權(quán)函數(shù)），分別與Lf和g作內(nèi)積式（4）即為三∽_g的近似算式

3）變換為矩陣方程：lmm=gm=則Llmn][an]=[gn]

4）矩陣求逆如果矩陣Llmn]是非奇異，則

注意：fn必須是線性無關(guān)，選擇恰當(dāng)可使∑anfn很快逼近fo權(quán)函數(shù)wm的選擇適當(dāng)也會是計算方便，

二1

當(dāng)fn=Wm時稱仂Ⅱ略金法（GarlerKin'smethord）

1.2基函數(shù)和權(quán)函數(shù)的選擇

基函數(shù)：MoM法的一個重要問題是基函數(shù)fn的選取，理論上有許多基函數(shù)可供選擇，只要它們是線性獨立的即可。但實際上，人們往往只能有少量的某些函數(shù)可較好地逼近待求地f通常選取地基函數(shù)應(yīng)使矩陣有較少的階次，求逆陣方便，收斂快等要求。

1）全域基函數(shù)：在待求函數(shù)廠的全部定義域中存在且不為“0”，如：

（1）冪級數(shù)

fn=x-xn+1

（3）麥克勞林級數(shù)n=InX2（n-l）

2）分域基函數(shù)：在待求函數(shù)廠的全部定義域中存在，但部分為“0”。如

（1）脈沖函數(shù)（分段均勻）

如圖1所示，其表達式為

（2）分段線性函數(shù)

如圖2所示，其表達式為

（4）二次插值函數(shù)

（5）正旋插值函數(shù)

權(quán)函數(shù)：作為對展開基函數(shù)的檢驗函數(shù)，常用的有以下兩類。

1）全部基函數(shù)均可作為權(quán)函數(shù)（或階段函數(shù)）

Wm的作用是對Lf-g系統(tǒng)加權(quán)（激勵）。通常Lf=g通過選取基函數(shù)已經(jīng)簡化為anlfn=g在N個待求的ai的方程，此時選N個激勵去列N個方程從而解出N個ai。一般這種激勵選δ函數(shù)是方便的。

當(dāng)選擇基函數(shù)等于權(quán)函數(shù)時，（厶=%時為伽略金法）稱為伽略金法。在特定的問題中，主要任務(wù)是選擇展開基函數(shù)和權(quán)函數(shù)，它們必須是線性無關(guān)的。展開基函數(shù)可以分為全域基和分域基，分域三角型基函數(shù)如圖1所示；權(quán)函數(shù)一般有點匹配法，伽略金法，最小二乘法等，伽遼金方法是選用和展開基函數(shù)相同的函數(shù)作為檢驗函數(shù)的方法。應(yīng)用矩量法需注意誤差分析、方程收斂性和積分奇異點處理等問題。δ2矩量法在線天線分析中的應(yīng)用

由于線天線在實際中應(yīng)用廣泛，而且是分析其它天線的基礎(chǔ)，所以用矩量法分析線天線有重要意義，多數(shù)文獻采用脈沖基函數(shù)和點匹配的方法分析線天線的Hallen積分方程[7]，本文詳細推導(dǎo)了線天線的Pocklington積分方程，應(yīng)用分域三角基和伽遼金的矩量法對半波對稱振子天線上電流分布特性進行了分析分析研究。δ2.1Pocklington積分方程的建立

考察如圖3所示的l=λ/2半波偶極天線，并作以下近似：電流沿導(dǎo)線軸流動；體電流密度，可以用線電流，近似；體電荷密度p用線荷密度σ近似；忽略天線端面的周向電流IΦ和和徑向電流Ip，認(rèn)為線上電流僅為長度變量z的函數(shù)，即I=I（z）與φ，p無關(guān)[6]。

當(dāng)半波振子天線上在加上饋電信號V時，或在接收電場作用時，近似有Eo=V/△z，陣子上將產(chǎn)生感應(yīng)電流，我們將電流分布按分域三角基函數(shù)展開。下面推導(dǎo)對應(yīng)的波克林頓積分方程。

令天線線長為，，線徑為a，且滿足>>a，天線端面處之Je和天線側(cè)面處之Jω略之不計，由于a<<λ，所以Je貢獻很小，可以忽略不計，又由于圓對稱，所以J=0，Jsz沿導(dǎo)體周界均勻分布，所以Jsz面電流分布可視為集中于軸線上的線電流，其表達式如下：根據(jù)以上假設(shè)可知動態(tài)矢量位A只有Az分量，即邊界條件EsIt+E1|l，=0（電場切向分量=0）導(dǎo)體表面處E|t=2Ez；Z向分量即為切向分量

利用洛淪滋條件

代人（12）式

在導(dǎo)體表面處

導(dǎo)體表面上

將Az代人（14）式，并交換微積分次序，可得

由導(dǎo)體邊界條件

式（16）即為Pocklington'sEquation，寫成算子形式為：

上述方程即是Pocklington，Pocklington積分方程激勵為E。對于一段d的電流元，其矢量位為

z位觀測點位置，z為源點位置，如圖1所示，對長度為2，的載流導(dǎo)體得Pocklington積分方程為

E，是由源場Ez，在陣子上產(chǎn)生感應(yīng)電流輻射產(chǎn)生的場，源場則來自施加在天線饋電端的電壓，由邊界

條件知在導(dǎo)體上的電場之和趨于零，則有

對Ez的表達式進行微分重組可得到如下形式

2.2矩量法求解Pocklington方程

根據(jù)第3節(jié)的模型，將電流函數(shù)，（z7）按圖所示的分域三角型基函數(shù)厶（z7）展開，并取有限近似，即選伽遼金方法，即取和展開基函數(shù)相同的檢驗函數(shù)記為厶（z），得

得到矩陣方程

根據(jù)以上推導(dǎo)，進行MATLAB仿真，得到偶極子天線的各項特性如圖4—圖7所示，從圖可以看出，該天線在E面的輻射方向圖呈8字形，在H面是一個全向輻射。

3偶極子天線在生物醫(yī)學(xué)上應(yīng)用

2013年6月24日出版的美國《紐約時報》印刷版刊文稱，正當(dāng)人們探討可穿戴計算設(shè)備時，形似藥丸的內(nèi)服式計算設(shè)備也逐漸開始興起。通過在藥丸中整合各種傳感器和發(fā)射器，便可以實現(xiàn)監(jiān)控人體內(nèi)部狀況的目的。如今也有一種新的可吞服電子產(chǎn)品，可以當(dāng)藥品或醫(yī)療器械。

“我給你開了兩枚電子藥丸。請每隔12小時服用一次。你吃下它之后，我會收到它發(fā)來的電子郵件，獲得你最新的內(nèi)臟器官運行數(shù)據(jù)。分析完這些數(shù)據(jù)，我會通知你詳細結(jié)果?！?/p>

如圖9所示的電子藥丸，包裹在膠囊內(nèi)部的主要部件是超微傳感器和發(fā)射器。和普通藥丸一樣，它可以用水或牛奶沖服。服用它就是監(jiān)控你的身體狀況，并將數(shù)據(jù)無線發(fā)送到外部設(shè)備中。由于偶極子天線的小型化和高集成性能，使得無線藥丸或膠囊中無線收發(fā)多采用偶極子天線來實現(xiàn)。4結(jié)論

本文分析了矩量法在半波偶極子線天線計算中應(yīng)用，通過采用分域三角基函數(shù)和伽遼金方法求解了Pocklington方程，用Matlab仿真實現(xiàn)偶極子天線的各項特性，計算結(jié)果與解析解比較吻合，表明矩量法在線天線計算中能得到較好的效果，當(dāng)然矩量法還可以用于其他天線的分析，如八木.宇田天線，陣列天線等，文中也簡單介紹了偶極子線在電子藥丸上的應(yīng)用。

參考文獻

[1]王秉中.計算電磁學(xué)[M].北京：科學(xué)出版社，2002.

[2]（美）哈林頓（R.F.Harrington）.計算電磁場的矩量法[M]王爾杰譯.北京：國防工業(yè)出版社，1981.

[3]姜光興，曹偉，朱洪波.基于不同積分方程的線天線矩量法分析[J].南京郵電大學(xué)學(xué)報，2005，25（1）：46-50.

[4]呂英華.計算電磁學(xué)的數(shù)值方法[M].北京：清華大學(xué)ILIJ版社，2006.

[5]（美）CleveB.Moler著.MATLAB數(shù)值計算[M]喻文健譯.北京：機械工業(yè)出版社，2006.

[6]黃惠芬，張少芳.使用新型頻率選擇表面反射板來提高超寬帶天線的性能[J].新型工業(yè)化，2013，3（12）：9-14.

[7]劉曉志，吳永剛.基于雙曲余弦函數(shù)的智能天線白適應(yīng)波束形成算法[J].新型工業(yè)化，2014，4（3）：74-78.

[8]劉寧，呂英華.一種用于人體中心網(wǎng)絡(luò)的超寬帶可穿戴織物天線[J].新型工業(yè)化，2011，1（6）：49-5.