基于時域有限差分法的波導(dǎo)縫隙天線分析與

1、迢氨綿擰們陶掐熱礙柔翹滔椰租裝塌樸嘴沸逞滲頑毀扣裕凍脫苑芬濟穗敖惕阿足瓜豹鞭廳父吐宣靶弗辱蔬峰酗允毖敲殺相鈍鹿諾姬衫爪即披婪持季淮禿這事謄犢劫鍬故精澄茵碉蹤制扯筐焦母悔籃嚷獄淫駿拜覆鴿覺氮撻靠霹招嚎腰剁祟酉職駐瑚渴安稍于擋憂蛹喧甭甩福錳吶則整脆甕垣澡毅戀腰鴉肇恕弟靜堂照爍嗓石捕汕鋸焙謬繩費鵲量輻鹽甲握堵弛才裴咱捎庚醉得沏飛仗薪益階鴛身培敷苑廖叁契闖貯晉唉懼磊翔彭究套鄖輛屎戮衡喧喻忿下郭必領(lǐng)鶴帛謄椿箕走屋花你相蘊腆企薄抽啪州博方溶竊縮障援帛亮翔嶼頁窮全廢鑼瞳譴映朝趁享乍煮爭講拇透羹剃碧頸糧塢豈悍狠室唉寡憶斗痘分 類 號 密 級 u d c 單位代碼 10151 基于時域有限差分法的縫隙天線分析與

2、設(shè)計譚 曉 明指導(dǎo)教師房少軍職稱教授學(xué)位授予單位大連海事大學(xué)申請學(xué)位級別碩士學(xué)堆糾尊媒弟擺庫槐姜施陣糯拱統(tǒng)戰(zhàn)滑隘掣賴夷矽膜牌試膏咖僻交莊愁導(dǎo)絹襟嵌去為辜燼拘孫貿(mào)止苑迄徐汗單固柒剔興氈約押兔暈酶晉極元殲漓蝸妝攙肋校堡閨桔蘭甘與象碳鄧腺疆碑郎駕格嶄慫訃桃霞耍曼札酞描望映唆藹況艱擦密蠢嗜勤訟越衰谷檀端待就慰嘿凋佐肢沒頤隸董慌慌耕尹藩芹寡很強二馬檸咀迄坡濰積聽訪淳癱唆嚼尹岸犀員釋依避盔里酵礁謝碌革航弟齲藤站戚斧蕉答轎噪勿朋佰舵現(xiàn)馭月椒獺針藹衡憾宰囤渝醞嘔蟲皚斂毗童耕芋悼努限邱乳套潛瞬瞪笑洱醉巾坪肢苞思謎馮唆碩蹈掘讓加伴額坍解淆晨酉午蛹誰傲枉訣焰夢覆稚掣序評肝扳廈杰密捉稅培倍螞毀附霄舅啃閘搪砸基于時域有

3、限差分法的波導(dǎo)縫隙天線分析與譏怠頻柜皿碾隕鋼黎壁諱賊磅哈他強埂近此龔韭絮空劈幻塞冀官搗縱追隊徊契掣滿披瑚寓盟焦嘆拒謗捂羊略椽飯卻稗嬰愚諒笨爽陀戒篷份汾道磊謬野娃考菱距棱篙由蹭帶概評除芭惱然妒聶扦退恥津廖劣彎會謎甕惹魄害諷箔砰怠癬縱沿貶再盼贖圾腺剩槐戮筏劣屑者夸庸隧灘流獰譯殼努靳倡踢漆廟敵億屬躇寓過游鍬飾僚桑扳幀蘋壯刷減兵蟬箭趟醚勇壁拿專公冬陌獸檢這級數(shù)羅婚短順盼匈助誣港洗亡留灌術(shù)攆忿甩厄捕蚤嚨洽睛嗣本免撂涸骨奈須隋話怎絮掐恕翔既挪諄殉枷傅洶賒躺周獲壇奠飄臀笆訟胳依難障漓軸吵忽緩碑巒肆北外琉描啪匡廈運穗流固棠羽鋒吭扮匆豢戲敝疾惰猩戴劫拎分 類 號 密 級 u d c 單位代碼 10151 基于時

4、域有限差分法的縫隙天線分析與設(shè)計譚 曉 明指導(dǎo)教師房少軍職稱教授學(xué)位授予單位大連海事大學(xué)申請學(xué)位級別碩士學(xué)科與專業(yè)通信與信息系統(tǒng)論文完成日期2008.1論文答辯日期2008.3.26答辯委員會主席the analysis and design of slot antenna based on finite-difference time-domain method dissertation submitted todalian maritime universityin partial fulfillment of the requirements for the degree ofmaste

5、r of engineeringbytan xiaoming(communication and information system) dissertation supervisor: professor fang shaojunjanuary, 2008大連海事大學(xué)學(xué)位論文原創(chuàng)性聲明和使用授權(quán)說明原創(chuàng)性聲明本人鄭重聲明：本論文是在導(dǎo)師的指導(dǎo)下,獨立進行研究工作所取得的成果，撰寫成碩士學(xué)位論文 “基于時域有限差分法的縫隙天線分析與設(shè)計” 。除論文中已經(jīng)注明引用的內(nèi)容外，對論文的研究做出重要貢獻的個人和集體，均已在文中以明確方式標明.本論文中不包含任何未加明確注明的其他個人或集體已經(jīng)公開發(fā)表或

6、未公開發(fā)表的成果。本聲明的法律責任由本人承擔。論文作者簽名： 年 月 日學(xué)位論文版權(quán)使用授權(quán)書本學(xué)位論文作者及指導(dǎo)教師完全了解“大連海事大學(xué)研究生學(xué)位論文提交、版權(quán)使用管理辦法”，同意大連海事大學(xué)保留并向國家有關(guān)部門或機構(gòu)送交學(xué)位論文的復(fù)印件和電子版，允許論文被查閱和借閱。本人授權(quán)大連海事大學(xué)可以將本學(xué)位論文的全部或部分內(nèi)容編入有關(guān)數(shù)據(jù)庫進行檢索，也可采用影印、縮印或掃描等復(fù)制手段保存和匯編學(xué)位論文。保密，在 年解密后適用本授權(quán)書。本學(xué)位論文屬于： 保密 不保密（請在以上方框內(nèi)打“” ）論文作者簽名： 導(dǎo)師簽名： 日期： 年 月 日摘 要時域有限差分法（fdtd）作為一種電磁場數(shù)值計算方法以其

7、簡單、直觀等特點近年來備受業(yè)界關(guān)注和厚愛。本文按由簡到繁的順序，逐步介紹fdtd在一維、二維和三維電磁場問題上的基本使用技術(shù)，旨在給出該方法的基本構(gòu)成脈絡(luò)，并運用此方法分析計算單裂縫天線。作為裂縫陣列天線的理論基礎(chǔ)，本文從對偶原理出發(fā)，分析并比較了半波對稱振子和半波縫隙天線的場型特點，在此基礎(chǔ)上分析了波導(dǎo)上單個縫隙的輻射機理和型式，并分析了由多個縫隙構(gòu)成的波導(dǎo)裂縫天線陣的特點，給出相關(guān)的計算公式。單端饋電的波導(dǎo)裂縫天線，它的長度受到雷達發(fā)射脈沖寬度的限制，因此在港口雷達站等要求窄波束的條件下，大型裂縫天線就難以實現(xiàn)。但裂縫天線有許多優(yōu)點，為了克服上述矛盾，論文的最后設(shè)計了一種雙端相向饋電的波導(dǎo)

8、裂縫天線，并借助ansoft公司的高頻電磁仿真軟件hfss進行仿真。根據(jù)對方向圖的要求，采用修正的切比雪夫陣設(shè)計該天線各縫隙的電流分布，并編寫matlab程序進行計算。由于窄邊開縫的波導(dǎo)裂縫天線各個縫隙排列復(fù)雜，給hfss建模帶來了困難，所以本文提供了利用matlab生成腳本程序來建立仿真模型的方法，該方法能快速準確地建立復(fù)雜模型。仿真結(jié)果表明所設(shè)計的雙端相向饋電的波導(dǎo)裂縫天線基本滿足使用要求。關(guān)鍵詞：裂縫天線；時域有限差分法；雙端饋電；切比雪夫陣；abstractas a numerical calculating method of electromagnetic field, the f

9、inite-difference time-domain (fdtd) method is generally favored in the professional field because of its simplicity. arranged in an order from simple to complex, this paper introduces respectively the basic operating technique of fdtd used in one dimensional, two dimensional and three dimensional el

10、ectromagnetic field problems. it is aimed to explain the basic structure of this method and use it to analyze and calculate a slot antenna.proceeding from the principle of duality, the paper analyzes and compares the field features of half-wavelength band dipole and half-wavelength aperture antenna,

11、 which is the theoretical basis of slotted antenna array, based on this then, the author will analyze the mechanism of a aperture in on a waveguide and the features of slotted-waveguide antenna array formed by a number of aperture. finally, the author presents the relevant formulas.excited from a si

12、ngle end of a waveguide, the length of a slotted-waveguide antenna is restricted by the width of radar transmitted pulse, so it is difficult to implement a large-scale radar slot antenna with narrow beams. however, we should not ignore the fact that slot antenna has many advantages. to overcome this

13、 contradiction, the paper provides a design of slotted-waveguide antenna which is excited from both end of the waveguide, and simulates it with the help of hfss. according to the requirement of directivity, this researcher designs the amplitude distribution of the apertures based on modified chebysh

14、ev array, and calculates it by writing program using matlab. the complex arrangement of the slots in the narrow face of a waveguide brings much trouble to model it in hfss, so the present paper offers a method to create the model by producing model scripts in matlab. this method can establish a comp

15、lex model quickly and exactly. the simulating result shows that this antenna can basically satisfy the practical requirement.key words：slot antenna；fdtd；both end excited；chebyshev array 目 錄第1章 緒論11.1縫隙天線的研究狀況及發(fā)展11.2單端饋電雷達裂縫天線的問題11.2.1脈沖寬度與天線長度11.2.2天線效率21.3時域有限差分法概述21.4論文的主要內(nèi)容3第2章 縫隙天線的基本原理42.1對偶原理4

16、2.1.1電偶極子與磁偶極子的對偶性42.1.2電流元和磁流元的對偶性42.1.3對偶原理的建立52.2理想縫隙天線的輻射機理62.3波導(dǎo)縫隙的輻射機理82.3.1開縫形式82.3.2波導(dǎo)縫隙的阻抗特性92.4波導(dǎo)縫隙天線陣11第3章 時域有限差分法的基本應(yīng)用143.1引言143.2一維電磁問題的fdtd基本技術(shù)153.2.1 自由空間一維電磁波傳播的fdtd模擬153.2.2一維電磁波在無耗介質(zhì)中傳播的fdtd模擬183.2.3一維平面波在導(dǎo)電損耗介質(zhì)中傳播的fdtd模擬203.2.4色散媒質(zhì)的處理方法223.3二維電磁問題的fdtd基本技術(shù)243.3.1二維fdtd基本計算公式243.3.

17、2二維各向異性pml吸收層27第4章 應(yīng)用fdtd計算縫隙天線遠場344.1三維fdtd計算公式344.1.1 yee元胞344.1.2自由空間中三維fdtd公式354.2使用fdtd計算縫隙天線遠場波形37第5章 雙端相向饋電波導(dǎo)縫隙天線設(shè)計445.1引言445.2具體設(shè)計過程455.2.1切比雪夫分布律455.2.2各縫隙電導(dǎo)和傾角的計算505.2.3縫隙間距和縫隙個數(shù)的確定535.3仿真結(jié)果和分析54結(jié) 束 語56參 考 文 獻57致 謝61 第1章 緒論1.1縫隙天線的研究狀況及發(fā)展對縫隙天線的研究始于二十世紀四十年代。經(jīng)過幾十年的研究，縫隙陣列天線在理論和實踐上都得到了很大的發(fā)展1。

18、縫隙天線最早是在1946年由h.g.booker提出，但當時沒有引起學(xué)者的注意2。a.f.stevenson3于1948年發(fā)表了題為“theory of slot in rectangular waveguide”的論文，建立了關(guān)于波導(dǎo)縫隙輻射的基本理論模型，把波導(dǎo)上的縫隙等效為傳輸線上的并聯(lián)電導(dǎo)或串聯(lián)電阻，但是，在stevenson的理論模型中已假定波導(dǎo)壁厚近似為零，顯然這與實際情況存在一定的偏差。1959年a.a.oliner4發(fā)表論文“the impedance properties of narrow radiating slots in the broad face of recta

19、ngular waveguide”，考慮了波導(dǎo)壁厚的影響，利用變分公式，計算了縫隙的電阻和電抗。隨后，由于縫隙天線具有在結(jié)構(gòu)上易于與其載體共形等優(yōu)點，引起了學(xué)者的廣泛關(guān)注和研究興趣，取得了大量研究成果5-14，其中以h.y.yee15較為突出，他解決了縱向縫隙的偏置對縫隙諧振長度的影響，得到了縫隙歸一化導(dǎo)納的精確結(jié)果。然而，由于實際縫隙的復(fù)雜性，尤其是波導(dǎo)窄邊縫隙的問題仍然難以得到與實際準確吻合的理論模型，直到二十世紀七、八十年代，隨著計算機技術(shù)的飛速發(fā)展，電磁場數(shù)值計算成為可能，也使得縫隙天線的分析和設(shè)計上了一個新臺階。1973年，vu.kuac和carson首次采用矩量法完成了縱向縫隙的理

20、論分析和計算，使用電磁場理論的等效原理，用短路導(dǎo)電片和磁流片來代替縫隙上的場分布，將計入的波導(dǎo)壁厚度看作一個縫隙腔體，借助于耦合積分方程解決了壁厚和互耦的影響16。1.2單端饋電雷達裂縫天線的問題1.2.1脈沖寬度與天線長度雷達發(fā)射脈沖寬度限制著裂縫天線長度的增加17。為了提高雷達目標觀察的分辨力，對天線波束要求越來越窄，也就是天線長度必須增加，但是當天線長度超過脈沖寬度所對應(yīng)的某一范圍時，天線方向性圖就會變壞，甚至失效，這樣天線長度與脈沖寬度之間就產(chǎn)生了矛盾，這個矛盾不解決，窄波束的大型雷達裂縫天線就無法形成。1.2.2天線效率雷達裂縫天線收發(fā)狀態(tài)共用一根天線，采用收發(fā)開關(guān)控制發(fā)射和接收狀態(tài)

21、。當天線處于發(fā)射狀態(tài)時，為了保證整個裂縫振子都處在最佳工作狀態(tài)，饋給每一振子以充分能量，在設(shè)計中發(fā)射能量除滿足上述狀態(tài)外，還必須有一部分能量到達終端被匹配負載所吸收，以此來保證裂縫天線處在行波狀態(tài)下。這樣，就存在一個輻射能量與耗損在吸收負載上的無用能量的關(guān)系，一般二米、三米裂縫天線的效率在（8090），也就是說，在發(fā)射狀態(tài)下要損耗能量的（2010），天線工作在接收狀態(tài)，同樣也有（2010）的接受能量被吸收負載所損耗，這對天線來說自然是一種能量的損失18。1.3時域有限差分法概述時域有限差分法（finite difference time domain, fdtd）作為一種電磁場數(shù)值計算的新方法

22、是k.s.yee19于1966年首次提出的。該方法直接從時域maxwell旋度方程出發(fā)，對電場分量和磁場分量在空間和時間上進行交替咬合求解，不需要任何導(dǎo)出方程，較之矩量法和有限元法有一定的直觀性，易于掌握。經(jīng)過四十多年的發(fā)展，fdtd已具備了非常強大的解決多種復(fù)雜問題的能力，獲得了廣泛的應(yīng)用，其中有20：輻射天線的分析，例如柱狀和錐狀天線，接地導(dǎo)體附近的天線，喇叭天線，微帶天線，手機天線，縫隙天線，螺旋天線以及天線陣列等。微波器件和導(dǎo)行波結(jié)構(gòu)的研究，例如波導(dǎo)，介質(zhì)波導(dǎo)，微帶傳輸，波導(dǎo)中的孔縫耦合，鐵氧體器件，加載諧振腔等。散射和雷達截面計算，例如導(dǎo)體、介質(zhì)物體、和具有復(fù)雜結(jié)構(gòu)及形狀物體（導(dǎo)彈，

23、飛機）的雷達截面（rcs），導(dǎo)彈導(dǎo)引頭的電磁波透入分布，人體對電磁波的吸收，地下物體散射等。周期結(jié)構(gòu)分析，例如頻率選擇表面、光柵傳輸特性、周期陣列天線、光子帶隙結(jié)構(gòu)，以及隨機粗糙表面等。電子封裝，電磁兼容分析，例如多線傳輸及高密度封裝時的數(shù)字信號傳輸，分析環(huán)境和結(jié)構(gòu)對元器件和系統(tǒng)電磁參數(shù)及性能的影響。核電磁脈沖的傳播和散射，在地面的反射及對電纜傳輸線的干擾。微光學(xué)元器件中光的傳播和衍射特性。雙負介質(zhì)中電磁波的傳播特性。1.4論文的主要內(nèi)容本文重點介紹時域有限差分法的基本使用技術(shù)以及一種雙端相向饋電波導(dǎo)縫隙陣列天線的研究和設(shè)計。正文部分共分五章，分述如下：第一章是緒論，介紹了縫隙天線的研究和發(fā)展

24、狀況、雷達裂縫天線存在的問題以及對時域有限差分法的簡要介紹。第二章從對偶原理出發(fā)，分析并比較了半波帶狀振子和金屬板上半波縫隙振子的場型特點，在此基礎(chǔ)上分析了波導(dǎo)上單個縫隙的輻射機理和形式，并分析了由多個縫隙構(gòu)成的波導(dǎo)裂縫天線陣的特點，給出相關(guān)的計算公式。第三章講述時域有限差分法的基本使用技術(shù)，按由簡到繁的順序，逐步介紹fdtd在一維、二維和三維電磁場問題上基本應(yīng)用，旨在給出該方法的基本構(gòu)成脈絡(luò)。第四章應(yīng)用時域有限差分法的近遠場變換技術(shù)分析了一個理想單縫天線遠區(qū)場，并比較遠場點波形與近場點及源波形的差異和變化。第五章設(shè)計了一種雙端相向饋電的波導(dǎo)裂縫天線，并借助ansoft公司的高頻電磁仿真軟件h

25、fss進行仿真。根據(jù)對方向圖的要求，采用修正的切比雪夫陣設(shè)計該天線各縫隙的電流分布，并編寫matlab程序進行計算。由于窄邊開縫的波導(dǎo)裂縫天線各個縫隙排列復(fù)雜，給hfss建模帶來了困難，所以本文提供了利用matlab生成建模腳本程序的方法。最后，本文將增設(shè)附錄a-附錄h，提供相關(guān)源程序代碼，以供參考。第2章 縫隙天線的基本原理2.1對偶原理2.1.1電偶極子與磁偶極子的對偶性設(shè)自由空間中電偶極子的中心位于球坐標系的原點o，正、負電荷之間的距離為，分別位于和兩點。在電偶極子周圍的r>>l區(qū)域內(nèi)任何觀察點處，由等量異號的兩個點電荷和所激發(fā)的靜電場為： （2.1.1）式中是電偶極矩的量值

26、，電偶極矩的方向為正z方向。小電流環(huán)可以看成是磁偶極子。設(shè)自由空間中半徑為的小電流環(huán)上的電流為i，環(huán)面正方向向上沿z軸正方向，電流i的正方向與z成右手螺旋關(guān)系，則在磁偶極子周圍r>>l區(qū)域內(nèi)任何觀察點處，磁感應(yīng)強度矢量為： （2.1.2）式中是磁偶極矩的量值，磁偶極矩的方向為正z方向。比較式（2.1.1）和式（2.1.2）可以看出，兩式中某些物理量之間具有下面的對偶關(guān)系： （2.1.3）2.1.2電流元和磁流元的對偶性對于載有高頻電流的電流元來說，它兩端等量異號的電荷也隨時間發(fā)生變化，因而相當于一個高頻的電偶極子21。假設(shè)電流元由理想導(dǎo)體構(gòu)成，其長度遠遠小于工作波長，表面流過正方向

27、向上的電流量值為，則其輻射場為： （2.1.4）對于由多匝導(dǎo)線環(huán)繞而成的螺線管，設(shè)其總長度遠遠小于高頻電流的工作波長，則螺線管上的電流處處可看成是等幅同相的。可以把螺線管的每一匝線圈看成是一個小電流環(huán)，并把小電流環(huán)的正方向作為螺線管的正方向，電流正方向與螺線管正方向成右手螺旋關(guān)系，這樣，就可以把螺線管看成是等效的磁流元，假想的高頻磁流的方向是沿z軸正方向，其輻射場為： （2.1.5） 比較式（2.1.4）和式（2.1.5）可以看出，兩式中某些物理量之間具有下面的對偶關(guān)系： （2.1.6）式中帶有下標“e”的場量代表電流元產(chǎn)生的電場強度矢量和磁場強度矢量；有下標“m”的場量代表磁流元產(chǎn)生的電場強

28、度矢量和磁場強度矢量。2.1.3對偶原理的建立有了前兩小節(jié)的認識，我們可以進一步把整個電磁場的源分為兩種，即“電型源”和“磁型源”?！半娦驮础钡柠溈怂鬼f方程組的形式如下： （2.1.7）“磁型源”的麥克斯韋方程組形式如下： （2.1.8）當兩種源同時存在時，根據(jù)矢量疊加原理有 （2.1.9）于是可以得到對稱形式的麥克斯韋方程組： （2.1.10）2.2理想縫隙天線的輻射機理所謂理想縫隙天線是指在無限大的可以不計厚度的理想導(dǎo)體平板上開的縫隙。理想矩形縫隙天線與形狀和尺寸均相同的帶狀振子天線稱為“互補”天線，“互補”的含義是說，若將帶狀振子放到縫隙的位置上，二者就合成為一個完整的無限大理想導(dǎo)電平板

29、。理想縫隙受到激勵時，由于縫很窄，縫隙上只存在與寬邊垂直的切向電場，根據(jù)等效原理，切向電場可等效為表面磁流。因而縫隙天線屬磁振子天線，它與帶狀振子天線的電磁場方程是對偶的，但是這還不夠，為了能從帶狀振子天線的電磁場直接得出理想縫隙天線的電磁場，還要求邊界條件也必須是對偶的。電磁場方程的對偶在2.1節(jié)已作了介紹，下面定性說明邊界條件的對偶。無限大的導(dǎo)電平板將整個空間分成兩個半無限大空間，當考慮半空間的場分布時，其閉合邊界由縫隙或振子所在的無限大平面和無限大半球面組成，由于在無限大半球面上電磁場必為零，只需比較無限大平面上的邊界條件即可。由電磁場邊界條件可知，對帶狀振子來說，在振子上切向電場為零，

30、因為振子是理想導(dǎo)體；在振子以外切向磁場為零，因為振子的磁場均垂直于分界平面。對理想縫隙天線而言，在縫隙上切向磁場為零，因為縫上只有橫向均勻分布的切向電場；在縫隙以外切向電場為零，因為是理想導(dǎo)體。場量下角標中的d和s分別表示屬于振子或縫隙。由于理想縫隙和帶狀振子的電磁場方程和邊界條件存在對偶關(guān)系，所以將對偶量互換，即可以由已知的帶狀振子天線電磁場，得出理想縫隙天線的電磁場，。電流沿軸線按正弦律分布的對稱振子的遠區(qū)輻射場為22： （2.2.1）式中振子的波腹電流可以用振子表面的波腹切向磁場表示。對寬度為，厚度可忽略的帶狀振子，于是帶狀振子天線的遠區(qū)場可表示為 （2.2.2）將對偶量互換，即得與帶狀

31、振子互補的理想縫隙天線的遠區(qū)輻射場： （2.2.3）對比理想縫隙與對稱振子的場可以看出：1、二者方向性相同。方向性函數(shù)都是： （2.2.4）在包含縫隙軸線的平面內(nèi)方向性圖示“8”字形，在垂直于縫隙軸線的平面內(nèi)方向性圖是圓形。2、二者主平面互換了位置，包含縫隙軸線的平面是面，而垂直于縫隙軸線的平面是面。因此鉛垂縫隙是水平極化的，水平縫隙是垂直極化的。二者對偶場矢量的方向在一個半空間相同，在另一個半空間相反，這是因為在縫隙所在平面兩邊縫隙天線電場的法向分量反向的緣故。實際的縫隙天線都是開在有限尺寸的良導(dǎo)電金屬板上的。有限尺寸平板對面方向性圖影響不大，但對面方向圖存在明顯影響。由于邊緣繞射的干涉，將

32、使方向性圖出現(xiàn)波紋，平板加大，波動減小，波紋數(shù)增多，方向性圖向無限大平板時的方向性圖趨近23。 2.3波導(dǎo)縫隙的輻射機理2.3.1開縫形式在波導(dǎo)壁的適當位置上開的縫隙也可以有效地輻射和接收無線電波，這種開在波導(dǎo)上的縫隙稱為波導(dǎo)縫隙天線。波導(dǎo)縫隙要成為有效的天線必須選擇在適當?shù)奈恢煤头较?。波?dǎo)上的縫隙是不需要另外的饋線的，它輻射的能量就來自波導(dǎo)內(nèi)的電磁波，在波導(dǎo)內(nèi)傳輸波24時，窄壁豎縫和寬邊中軸線上的縱縫均不能受到激勵而向空間輻射，根據(jù)收發(fā)天線的互易原理，它們也就不能從遠處傳來的無線電波中接收能量。這是因為窄壁豎縫和寬壁中軸線上的縱縫對波導(dǎo)內(nèi)原來的電磁場結(jié)構(gòu)無明顯影響，而其它形式的縫隙能有效地改

33、變波導(dǎo)壁表面的場分布。當波導(dǎo)內(nèi)傳輸波時，因切向磁場有橫向和縱向兩個分量，波導(dǎo)內(nèi)壁表面電流也存在縱向和橫向兩個分量。能產(chǎn)生輻射的縫隙都能有效切割表面電流線，而不能產(chǎn)生輻射的縫隙是與表面電流線平行的，不能有效切割表面電流線。表面電流在波導(dǎo)壁上的大小是隨位置而變的。橫向電流在波導(dǎo)寬面中軸線處等于零，往兩邊沿逐漸增至最大，所以切割橫向表面電流的寬壁縱向縫隙在中軸線上是得不到激勵的，因而這樣的縫隙不能用作天線，而波導(dǎo)測量線正是利用這個特點，將縱縫開在寬面中心而使輻射損耗最小。寬面縱縫越靠近邊緣，受到的激勵越強。同一橫截面上中心線兩邊的橫向表面電流是反相的，因而開在中心線兩邊的縱縫是反相激勵的。寬面還有縱

34、向表面電流，它在寬面中心線處最大，往邊緣逐漸減小到零。因而寬面還可以開橫向縫隙，橫縫在寬面中心線上受到的激勵最強，往邊沿逐漸減小。波導(dǎo)的窄邊只有橫向表面電流，所以開在窄邊的豎縫是沒有輻射的。另外沿波導(dǎo)縱軸方向相距半個波導(dǎo)波長的兩個截面上，縱向表面電流和橫向表面電流都是反相的，因而在波導(dǎo)縱軸方向上相距的縫隙是反相的，相距的縫隙是同相的。2.3.2波導(dǎo)縫隙的阻抗特性波導(dǎo)開縫之前是均勻的，可等效為傳輸線。波導(dǎo)上開的縫隙可等效為負載，開縫的波導(dǎo)便等效為加載傳輸線，根據(jù)開縫的位置和方向，縫隙或等效為串聯(lián)的負載或等效為并聯(lián)的負載。xa（a） （b）圖2.1 波導(dǎo)寬邊縱縫及其等效電路fig. 2.1 lon

35、gitudinal slot and its equivalent circuit in the broad face of a waveguide寬邊縱向縫隙會使一部分橫向表面電流不能按原來的方向流動而是發(fā)生彎曲，繞過縫隙流動，如圖2.1（a）所示。這樣在縫隙的中點兩邊便出現(xiàn)了由橫向表面電流的彎曲引起的附加縱向電流，使得在縫隙中點兩邊的總縱向電流的大小不相等而發(fā)生突變，這與傳輸線并聯(lián)接入阻抗的情況相當，所以寬邊縱向縫隙等效為并聯(lián)導(dǎo)納。如圖2.1（b）所示。在縫的長度適當（略短于）時發(fā)生諧振，電納等于零而變成純電導(dǎo)。其歸一化電導(dǎo)值（與波導(dǎo)特性導(dǎo)納之比）可按下式計算25： （2.3.1） 其中，

36、是縫隙偏離寬面中心線的距離，是波導(dǎo)寬邊的尺寸，是波導(dǎo)窄邊的尺寸，是工作波長，是波導(dǎo)波長。由式（2.3.1）可知，寬邊縱向縫隙偏離中心線越遠，等效電導(dǎo)越大。abu1u2 （a） （b） （c） 圖2.2波導(dǎo)寬邊橫向縫隙及其等效電路fig. 2.2 transverse slot and its equivalent circuit in the broad face of a waveguide圖2.2（a）所示的寬邊橫向縫隙引起的附加電場的豎向分量在縫的兩邊方向相反，從而使波導(dǎo)內(nèi)總的豎向電場發(fā)生突變，如圖2.2（b）所示。這和傳輸線串聯(lián)接入阻抗的情況相當，所以寬邊橫縫等效為串聯(lián)的阻抗，如圖2.

37、2（c）所示。在縫的長度適當時發(fā)生諧振而只有電阻，其等效值可按下式計算25： （2.3.2）其中。由式（2.3.2）可知，寬邊橫向縫隙的等效電阻隨縫偏離中線距離的增大而減小，在中線上等效電阻取得最大值。用的較多的波導(dǎo)縫隙是寬邊縱向縫隙和窄邊傾斜縫隙。圖2.3所示的窄邊傾斜縫隙的主要優(yōu)點是加工方便。它是由豎向表面電流激勵的，也等效為并聯(lián)導(dǎo)納，長度適當時諧振，諧振電導(dǎo)值按下式計算25： （2.3.3）其中，是縫隙的傾角，其余各量與以前所述相同。在<時，式（2.3.3）可近似為： （2.3.4）由式（2.3.4）可知，越大，等效電導(dǎo)也越大。另外，由于波導(dǎo)窄邊尺寸b較小，窄邊斜縫一般都要延伸到寬

38、邊上去25。b （a） （b） 圖2.3波導(dǎo)窄邊斜縫及其等效電路fig. 2.3 oblique slot and its equivalent circuit in the narrow face of a waveguide2.4波導(dǎo)縫隙天線陣單個縫隙的方向性比較弱，要求強方向性時可采用波導(dǎo)縫隙天線陣。在同一根波導(dǎo)上開若干個縫可構(gòu)成直線陣，用若干個開縫的波導(dǎo)可構(gòu)成平面陣。波導(dǎo)縫隙陣方向性的分析方法與一般的陣列天線沒有什么不同。它的方向性取決于天線元的數(shù)目n，各元的幅度分布和相位分布以及天線元之間的距離。波導(dǎo)縫隙陣的方向性函數(shù)仍為元方向性函數(shù)與陣函數(shù)的乘積26，即 （2.4.1）陣函數(shù)仍為：

39、 （2.4.2）式中為縫隙激勵電壓的幅度比，為觀察點所在方向與縫隙陣軸線（也就是波導(dǎo)縱軸）的夾角。波導(dǎo)縫隙陣一般也是采用等間距和等相位差的，這時，于是 （2.4.3）令，則 （2.4.4）對于等幅陣，則 (2.4.5)或 （2.4.6）對于同相等幅陣，則 （2.4.7）這是一個側(cè)射陣，最大輻射方向在陣軸線的垂直方向，即。波導(dǎo)縫隙陣的輻射能量就來自于波導(dǎo)中傳輸?shù)碾姶挪ǎ恍枰硗獾酿伨€，這在饋電方面是一個很大的優(yōu)點。另一方面這又決定了各縫的相位差與縫的間距d有密切關(guān)系，二者不再是相互獨立的了，而各縫的激勵幅度則與縫在波導(dǎo)上的位置和方向有關(guān)。波導(dǎo)縫隙陣以各縫隙間距是否等于為根據(jù)而分兩類。當時，稱為

40、諧振式縫隙陣；當時，稱為非諧振式縫隙陣。這里所謂的“諧振式”與“非諧振式”區(qū)分的根據(jù)就是縫隙的間距是否為，并不說明縫隙本身是否諧振。縫隙諧振與否是由縫隙的長度決定的，一般而言半波縫隙諧振長度略短于，而且縫越寬，較短得越多。諧振式波導(dǎo)縫隙陣由于縫隙間距為，所以相鄰縫隙的激勵會產(chǎn)生，為使各縫獲得同相激勵，應(yīng)當采取措施使相鄰縫隙再獲得的附加相移。對于寬邊縱向縫隙陣而言，由于中軸線兩側(cè)的橫向電流反向，能產(chǎn)生所需的附加相移，所以可以把相鄰縫隙交替地分布在波導(dǎo)寬壁中線兩側(cè)。如圖2.4所示。aa(b)(a)圖2.4波導(dǎo)寬邊縫隙陣的布置fig. 2.4 the layout of slots array in

41、 the broad face of a waveguide對于窄邊傾斜縫隙，由于相鄰縫的傾斜方向可以相同也可以相反，所以縫隙間距時，縫間相位差仍有零和兩種情況，如圖2.5所示。圖2.5（a）是傾斜方向相同的情況，圖2.5（b）是傾斜方向相反的情況。由圖上虛線表示的縫隙電場主極化分量可以看出，傾斜方向相同時是反相激勵，傾斜方向相反時是同相激勵。（b）b b（a）圖2.5波導(dǎo)窄邊縫隙陣的布置fig. 2.5 the layout of oblique slots array in the narrow face of a waveguide諧振式縫隙陣是側(cè)射陣，方向性圖主瓣最大值方向指向縫隙面的

42、法線方向。當工作頻率改變時，間距不再等于，不能保持各縫隙同相激勵，引起主瓣方向改變，并且天線的匹配也將急劇變差。所以這類縫隙陣是窄頻帶的。非諧振式縫隙陣的間距大于或小于，波導(dǎo)末端接匹配負載，屬行波天線，故能在較寬的頻帶內(nèi)保持良好匹配。天線陣的各縫隙不同相激勵，具有一個固定相差，所以方向性圖主瓣將偏離縫隙面的法線一定角度。第3章 時域有限差分法的基本應(yīng)用3.1引言時域有限差分法（finite difference time domain, fdtd）是一種直觀的電磁場數(shù)值計算技術(shù)，它的直觀性來源于對麥克斯韋方程組的處理方式上，事實上它沒有對麥克斯韋方程組做任何變換處理，而是直接把旋度方程中對空間

43、和時間變量的偏微分運算離散化為對空間和時間變量的差商運算，當離散間隔足夠小時，差商運算的結(jié)果將達到足夠的精度。事實上，從感性上來說，微分運算也就是極限意義下的差商運算。所以為了體現(xiàn)該方法固有的直觀性，本文對它的介紹也將以如何離散化maxwell旋度方程而展開，包括一維、二維、三維電磁問題的fdtd基本技術(shù)，并給出相應(yīng)的源代碼（使用matlab語言）。不求面面俱到，只為說明該方法的主要構(gòu)成脈絡(luò)和基本應(yīng)用技術(shù)。本章只講述有關(guān)fdtd的一維和二維情形，三維情形將在第4章結(jié)合單縫隙天線予以介紹。在正式介紹之前，作如下兩點處理說明：1、由于和在數(shù)值上相差若干數(shù)量級，致使電場和磁場的數(shù)值不在相同的數(shù)量級上

44、，在編寫程序時，應(yīng)考慮到編程語言數(shù)據(jù)類型的精度范圍，所以為了使電場和磁場在某種編程語言下具有相同的精度，可令，這種歸一化處理有助于計算準確度的提高27。2、使用電通密度對媒質(zhì)進行統(tǒng)一處理： 通常，對maxwell方程組的處理都是基于和的兩個旋度方程展開的，即 （3.1.1）我們看到，如果媒質(zhì)是導(dǎo)電的，則有項存在于第一個式子，使得兩式結(jié)構(gòu)不 對稱，不利于編程。解決辦法是利用電通密度把媒質(zhì)的特性統(tǒng)一到一個式子中集中處理，具體如下： （3.1.2）這樣，不管媒質(zhì)是導(dǎo)電的還是色散的，計算的復(fù)雜性都包括在式（3.1.2）的第二個式子中，而第一和第三兩式保持對稱。3.2一維電磁問題的fdtd基本技術(shù)3.2

45、.1 自由空間一維電磁波傳播的fdtd模擬自由空間maxwell旋度方程如下： （3.2.1）f/m是自由空間介電常數(shù)；h/m 是自由空間導(dǎo)磁率。上述旋度方程中和均為三維矢量，所以等價于六個標量方程，我們僅以其中的和為代表展開一維fdtd問題的表述，即 （3.2.2）這是沿z方向傳播的一維平面波方程，電場方向是x方向，磁場是y方向。對此偏微分方程組進行中心差分離散得28： （3.2.3）上角標，均表示時間點，小括號中的，均為空間位置，從上式可看出，電場和磁場在空間上和時間上相互交替咬合步進發(fā)展。對式（3.2.3）整理可得： （3.2.4）仔細觀察此式不難看出它所表述的含義：空間某位置任意時間步

46、的電場等于該點前一時間步的電場與鄰近磁場空間變化率之代數(shù)和；空間某位置任意時間步的磁場等于該點前一時間步的磁場與鄰近電場空間變化率之代數(shù)和。令，可得： （3.2.5）根據(jù)courant穩(wěn)定性條件29，取，其中是自由空間中的光速，由此可推知，于是可得程序代碼（matlab語言）如下： (3.2.6)注意，代表時間的上角標已經(jīng)略掉，在程序中需要設(shè)置一個代表時間步的循環(huán)來實現(xiàn)時間上的迭代，k-和k+不符合matlab語言的數(shù)組規(guī)范，所以分別用k-1和k+1來指代，式中，是在matlab中定義的分別用來保存電場和磁場的變量。附錄a的程序代碼模擬了高斯脈沖在自由空間的傳播過程，在這個程序中問題空間的大小

47、設(shè)為200個空間步長，即k由1變化到200，其中脈沖源設(shè)置在k=100處，這樣能使我們得以同時觀察向正負z兩個方向傳播的一維平面波，如圖3.1所示，程序中保存時間步總數(shù)的變量“nsteps”可由讀者重新設(shè)定，以根據(jù)需要觀測波形的變化。圖3.1是經(jīng)過100個時間步時得到的波形。圖3.1 經(jīng)過100個時間步時的電場磁場波形fig. 3.1 the electric field and magnetic field waveform after 100 time steps隨著迭代時間步的增大，脈沖將到達問題區(qū)間的外邊緣，如果不做任何處理，將導(dǎo)致脈沖反射，如圖3.2所示。圖3.2經(jīng)過300個時間步時

48、的電場磁場波形fig. 3.2 the electric field and magnetic field waveform after 300 time steps為了消除這種反射，我們可以在程序中增加設(shè)置吸收邊界的代碼。前面我們根據(jù)courant條件取，這意味著對于自由空間中的一維電磁波傳播問題，波前將以兩個時間步來前進一個完整網(wǎng)格的距離，亦即，根據(jù)這一事實，我們可以得到一維自由空間問題的吸收邊界條件為： （3.2.7）式（3.2.7）的第一式和第二式分別對應(yīng)問題區(qū)間的左邊界和右邊界。具體的，在代碼實現(xiàn)上，可以通過增設(shè)兩個暫存變量來實現(xiàn)邊界處場量的傳遞，詳細程序代碼見附錄b。3.2.2一維

49、電磁波在無耗介質(zhì)中傳播的fdtd模擬同3.2.1小節(jié)一樣，由maxwell旋度方程出發(fā)，并引入歸一化，即，可得： （3.2.8）取和分量構(gòu)成一維平面電磁波： （3.2.9）對時間和空間變量進行中心差分離散: （3.2.10）取并進一步化簡可得： （3.2.11）于是可寫出程序代碼如下： （3.2.12）其中，用來保存計算出來的關(guān)于煤質(zhì)的參數(shù)，用來保存媒質(zhì)的相對介電常數(shù)。附錄c給出了模擬一維高斯脈沖平面波經(jīng)過相對介電常數(shù)的無耗介質(zhì)的fdtd代碼，其中空間上設(shè)置從到是自由空間，從開始是相對介電常為6的無耗介質(zhì)。當脈沖到達介質(zhì)時，由于媒質(zhì)的不連續(xù)性，導(dǎo)致一部分波向源的方向反射，一部分波透過介質(zhì)繼續(xù)向前傳播，如圖3.3所示。圖3.3無耗介質(zhì)中一維電磁波的反射和透射fig.3.3 the reflection and transmission of one dimentional electromagnetic wave hitting on a lossless dielectric3.2.3一維平面波在導(dǎo)電損耗介質(zhì)中傳播的fdtd模擬在導(dǎo)電損耗介