電磁場的數(shù)值計算方法

電磁場的數(shù)值計算方法摘要：數(shù)值計算方法是一種研究并解決數(shù)學問題數(shù)值近似解的方法，廣泛運用于電氣、軍事、經(jīng)濟、生態(tài)、醫(yī)療、天文、地質(zhì)等眾多領域。本文綜述了電磁場數(shù)值計算方法的發(fā)展歷史、分類，詳細介紹了三種典型的數(shù)值計算方法—有限差分法、有限元法、矩量法,對每種方法的解題思路、原理、步驟、特點、應用進行了詳細闡述,并就不同方法的區(qū)別進行了深入分析,最后對電磁場數(shù)值計算方法的應用前景作了初步探討。關(guān)鍵詞：電磁場；數(shù)值計算；有限差分法；有限元法；矩量法引言自從1864年Maxwell建立了統(tǒng)一的電磁場理論，并得出著名的Maxwell方程以來，經(jīng)典的數(shù)學分析方法是一百多年來電磁學學科發(fā)展中一個極為重要的手段,圍繞電磁分布邊值問題的求解國內(nèi)外專家學者做了大量的工作。在數(shù)值計算方法之前,電磁分布的邊值問題的研究方法主要是解析法，但其推導過程相當繁瑣和困難,缺乏通用性,可求解的問題非常有限。上個世紀六十年代以來,伴隨著電子計算機技術(shù)的飛速發(fā)展,多種電磁場數(shù)值計算方法不斷涌現(xiàn),并得到廣泛地應用,相對于解析法而言,數(shù)值計算方法受邊界形狀的約束大為減少,可以解決各種類型的復雜問題。但各種數(shù)值計算方法都有一定的局限性,一個復雜的問題往往難以依靠一種單一方法解決，因此如何充分發(fā)揮各種方法的優(yōu)勢,取長補短,將多種方法結(jié)合起來解決實際問題,即混合法的研究和應用已日益受到人們的關(guān)注。本文綜述電磁場的數(shù)值計算方法，對三種常用的電磁場數(shù)值計算方法進行分類和比較。1電磁場數(shù)值計算方法的發(fā)展歷史在上世紀四十年代，就有人試探用數(shù)值計算的方法來求解具有簡單邊界的電磁場問題，如采用Ritz法[1]，以多項式在整個求解場域范圍內(nèi)整體逼近二階偏微分方程在求解域中的解。五十年代，采用差分方程近似二階偏微分方程，誕生了有限差分數(shù)值計算方法，開始是人工計算，后來采用機械式的手搖計算機計算，使簡單、直觀的有限差分法得到應用和發(fā)展，該方法曾在歐、美風行一時。1964年美國加州大學學者Winslow以矢量位為求解變量，用有限差分法在計算機上成功地解算了二維非線性磁場[1]，此后有限差分法在工程電磁場計算領域大為發(fā)展。1965年，Winslow首先將有限元法從力學界引入電氣工程中，1969年加拿大MeGill大學P.P.Silvester運用有限元法成功地進行了波導的計算[2]；七十年代初，P.P.Silvester和M.V.K.Chari合作將有限元法應用于二維非線性磁場的計算，成功地計算了直流電機、同步電機的恒定磁場。此后有關(guān)有限元法探討的論文越來越多，有限元法運用的范圍由靜態(tài)場到渦流場到輻射場，由線性場到非線性場，由各向同性媒質(zhì)到各向異性、要考慮磁滯損耗，由工程電磁場到生物電磁場等等。有人認為有限元法是求解工程電磁場的最有效最成功的方法。有限元法和有限差分法都是求解邊值問題的方法，屬于微分方程法。對于開區(qū)域或要求求解連續(xù)分布場量的區(qū)域，這類方法就會受到自身的限制。1972年英國盧瑟福實驗室的C.W.Trowbridge等人提出了積分方程法的思想，給出了二維、三維場問題的離散形式[2]，由于此種方法只需離散源區(qū)，不需考慮邊界條件，所以它較好地解決了無界開域場和要求連續(xù)計算場量的問題。該方法計算精度高，但計算量很大。該實驗室Sinkin等人又在積分方程法基礎上提出了邊界積分方程法(又稱邊界元法)，用此解決線性場的計算，計算量大為減小。此后該室的學者們將積分方程與微分方程法結(jié)合起來，提出了求解三維靜磁場的雙標量位法等。在解決天線輻射場、散射場問題中，矩量法是一個很重要的數(shù)值計算方法。1968年R.F.Harrington發(fā)表了專著“FieldcomputationbyMomentMethod”，對散射場、天線輻射場、波導場等方面的問題起了很好的推進作用。除以上所介紹的方法外，隨著電磁場數(shù)值分析的不斷發(fā)展，各種新方法不斷涌現(xiàn)，如計算電場的模擬電荷法，最小二乘配點法，求解磁場的模擬電流法，以及計算場的圖論模型法，快速Fourier變換法、有限體元法、無網(wǎng)格計算法等等。各種方法互相配合，出現(xiàn)了一些混合方法，如：矩量法—模擬電荷法、模擬電荷法—有限元法、有限元法—邊界元法等，有效地解決了一些實際問題。近年來人工神經(jīng)網(wǎng)絡，小波理論[3]等也引入了電磁場的數(shù)值計算中，瞬態(tài)電磁場計算如時域有限差分法的應用有了長足的發(fā)展。總之隨著現(xiàn)有的電磁場數(shù)值計算方法的不斷深入發(fā)展、提高和完善，新的方法不斷產(chǎn)生。在電磁場的數(shù)值解法不斷發(fā)展的同時，人們并沒有忘記長期以來所運用的解析方法。解析法計算結(jié)果精確，且可以用解析式表達計算結(jié)果，受這些特點吸引，解析法與數(shù)值計算方法相結(jié)合形成的半解析法應運而生，也成為了一種主流解算方法，并還在不斷發(fā)展。電磁場數(shù)值計算方法發(fā)展走向成熟的一個重要標志是：成熟的方法越來越多地應用于工程實際問題中，商業(yè)化通用軟件包不斷出現(xiàn)[4]。一個商業(yè)化軟件包通常由下面幾部分組成：數(shù)據(jù)定義：幾何尺寸、材料性能參數(shù)、邊界條件前處理模擬化： 空調(diào)剖分、網(wǎng)格自動產(chǎn)生、節(jié)點形成網(wǎng)格圖形顯示離散方程組系數(shù)矩陣形成數(shù)據(jù)處理求解代數(shù)方程組 非線性疊代按要求輸出計算結(jié)果場圖顯示(含線性媒質(zhì)和非線性媒質(zhì)區(qū))后處理 受力和損耗計算與圖形顯示局部場域分布的精細計算與顯示以上三部分中前、后處理占用了軟件包語句的90%以上，編程的主要工作量在此，而數(shù)據(jù)處理，也就是我們目前正在學習的數(shù)值計算方法僅占軟件語句的10%以內(nèi)，但它卻是占用計算機內(nèi)存量和消耗CPU時間的主要部分。電磁場數(shù)值計算方法的分類求解電磁問題的最終要求就是獲得滿足實際條件的Maxwell方程的解，借助于計算數(shù)學中的數(shù)值算法能夠得到大多數(shù)電磁問題的近似解。數(shù)值算法的基本思想[5]就是把連續(xù)變量函數(shù)離散化，把微分方程化為差分方程；把積分方程化為有限和的形式，從而建立起收斂的代數(shù)方程組，然后利用計算機技術(shù)進行求解。數(shù)值計算方法從求解方程的形式看，主要分為積分方程法和微分方程法兩大類。積分方程法主要有矩量法和邊界元法，微分方程法主要有有限差分法和有限元法。對兩種方程法的比較，如表一所示。表一積分方程法和微分方程法的比較積分方程法微分方程法共性對場問題的處理是一致的，即需離散化場域，結(jié)果是數(shù)值解不同點離散域僅在場源區(qū)，無需對整個場域離散整個場域計算對象場量先求位函數(shù)，再求場量求解域可在場域內(nèi)某一局部區(qū)域求解，也可在全場域內(nèi)求解全場域內(nèi)求解計算程度較高較低應用不適用邊界區(qū)域復雜的場域邊界形狀復雜的場域較易處理聯(lián)系兩種方法的結(jié)合形式，可處理較復雜的電磁場問題幾種重要的數(shù)值計算方法有限差分法在電磁場數(shù)值計算方法中，有限差分法是應用最早的一種方法。有限差分法以其概念清晰，方法簡單、直觀，有大致固定的處理和計算模式，具有一定的通用性等特點，在電磁場數(shù)值分析領域內(nèi)得到了廣泛的應用。有限差分法的基本原理有限差分法的基本思想是把連續(xù)的定解區(qū)域用有限個離散點構(gòu)成的網(wǎng)格來代替，這些離散點稱作網(wǎng)格的節(jié)點；把連續(xù)定解區(qū)域上的連續(xù)變量的函數(shù)用在網(wǎng)格上定義的離散變量函數(shù)來近似；把原方程和定解條件中的微商用差商來近似，積分用積分和來近似，于是原微分方程和定解條件就近似地代之以代數(shù)方程組，即有限差分方程組，解此方程組就可以得到原問題在離散點上的近似解。然后再利用插值方法便可以從離散解得到定解問題在整個區(qū)域上的近似解。差分與差商設函數(shù)f(x)的自變量x有一小增量xh，則f(x)的增量為f(x)f(xh)f(x)（3.1）f(x)為函數(shù)f(x)的一階差分。當增量h足夠小，差分f與微分df之間的差才足夠小。一階差分f是自變量x的函數(shù)。按式（1）,計算f(x)的差分2f(x)稱二階差分，且2f(x)f(xh)f(x)（3.2）函數(shù)f(x)的一階導數(shù)f'(x)為fxdflimfx dx xx0應用差分，f'(x)可表示為f(x)f(xh)f(x)f'(x) （3.3） x h故f'(x)可表示為差分f(x)除以有限小差分x的商，稱為差商。同理，函數(shù)f(x)的二階導數(shù)f''(x)可表示為 d2f1df df ( )dx2xdxxxdxx1f(xh)f(x)f(x)f(xh) h h （3.4）hf(xh)2f(x)f(xh)h23.1.3差分方程的構(gòu)造現(xiàn)以二維靜態(tài)電、磁場泊松方程的第一類邊值問題為例，來具體闡明有限差分法的應用。設具有平行平面場特征的電磁場場域D,如圖1所示,為一由閉合邊界L所界定的平面域，其定解條件可表述為2ux,y2u2uFx,yx,yD（3.5）x2y2ux,yfx,y（3.6）L對于所給定的偏微分方程定解問題，應用有限差分法，首先需從網(wǎng)格剖分著手決定離散點的分布方式。原則上，可以采用任意的網(wǎng)格剖分方式，但這將直接影響所得差分方程的具體內(nèi)容，進而影響解題的經(jīng)濟性與計算精度。為簡化問題，通常采用完全有規(guī)律的分布方式，這樣在每個離散點上就能得出相同形式的差分方程，有效地提高解題速度，因而經(jīng)常采用正方形的網(wǎng)格的剖分方式?，F(xiàn)即以這種正方形網(wǎng)格剖分場域D，也就是說，用分別與x,y兩坐標軸平行的兩簇等距網(wǎng)格線來生成正方形網(wǎng)格，即L03ML03MDyhh124xx(i0,1,2)yyjh j (j0,1,2) yj1yh為步長，網(wǎng)格線的交點Ox,y稱為節(jié)yj i jj1點，這樣D域就離散化為由網(wǎng)格節(jié)點標成的離散點得集合。對場域D中節(jié)點Ox,y是一典型 i j節(jié)點，它與周圍的1，2，3和4點構(gòu)成一x個對稱星型。設這些離散點上待求函數(shù)的 i1 i1圖1正方形網(wǎng)格劃分近似值記為u0u(i,j)，u1u(i1,j)，u2u(i,j1)，u3u(i1,j)，u4u(i,j1)則式（6）可近似離散化為1u(i1,j)2u(i,j)u(i1,j)1u(i,j1)2u(i,j)u(i,j1)F（3.7）h2 h2即u(i1,j)u(i,j1)u(i1,j)u(i,j1)4u(i,j)h2F（3.8）若式（6）F=0,則節(jié)點O上函數(shù)u的值等于其四周相鄰點函數(shù)值的平均。因為差分方程（7），（8）只出現(xiàn)待求函數(shù)u在點Ox,y及其四個臨近點的值，故i j稱之為五點差分格式[6]，根據(jù)差分方程組解出各離散點處的待求函數(shù)值。3.2有限元法傳統(tǒng)的變分法在20世紀二三十年代為其新型時期，理論上發(fā)展很快，各種變分問題的最后求解都可歸結(jié)為解尤拉方程的邊值問題，然而只有在一些特殊情況下尤拉方程才能求出精確解，在大多數(shù)情況下，尤拉方程的精確解無法求出。四五十年代，隨著計算機的出現(xiàn)，使其在實際應用中逐漸為比較靈活、通用的有限差分法所替代。但是，有限差分法在理論上沒有以變分原理為基礎，因而其收斂性和數(shù)值穩(wěn)定性往往得不到保證。隨后發(fā)展形成的有限元法正是變分法與有限差分法相結(jié)合的成果，它取長補短地在理論上以變分原理為基礎，在具體方法構(gòu)造上又利用了有限差分法網(wǎng)格離散化處理的思想[7]。3.2.1有限元法的基本原理有限元法是以變分原理為基礎，將要求解的微分方程型數(shù)學模型—邊值問題，首先轉(zhuǎn)化為相應的變分問題，即泛函求極 A(x,y) xdxdsdxdsB(x2,y2)3.2.2泛函、變分問題簡介在微積分學形成初期，以數(shù)學物理問題為背景，與多元函數(shù)的極值問題相對應，已在幾y何、力學上提出了若干個求解泛函極值的問題。如圖2中的質(zhì)點最速降線問題所述，質(zhì)點A圖2最速降線問題從定點圖2最速降線問題 1 1 2 2使滑行時間最短的質(zhì)點下滑軌道yy(x)。圖示滑行弧段ds所需時間為dssecdx 1y2dxdt v 2gy 2gy滑行總時間為J[y(x)]T[y(x)]Tdtx21y2dx（3.9） 0 x 2gy1（3.9）式JJ[y(x)]不僅取決于積分端點x和x，而且取決于yy(x)的選取。 1 2J取決于y(x)，所以J是函數(shù)y(x)的函數(shù)，稱之為y(x)的泛函，記作J[y(x)]。于是所述之最速降線問題，在數(shù)學上就歸結(jié)為研究泛函J[y(x)]的極值問題，即x21y2dxminJ[y(x)] x 2gy （3.10） 1 y(x)0 y(x)y 1 2 2泛函的極值問題就稱為變分問題。對一般問題而言，可導出下列對應于一個自變量x、單個函數(shù)y(x)及其導數(shù)y(x)的已知函數(shù)J[y]x2F(x,y,y)dx（3.11）x1式中F為x、y和y的已．．．．知函數(shù)。泛函J[y]的自變量不是一般的自變量，而是一個或幾個函數(shù)所屬的函數(shù)族y(x)。在端點x和x上分別等于給定值的無．．數(shù) 1 2個．函數(shù)y(x)中，僅有一個y(x)能使定積分J[y]達到極小值，此函y(x)數(shù)稱為極值函數(shù)。因此，變分問題就在于尋求使泛函達到極值的該極值函數(shù)y(x)，即分析研究泛函的極值問題。3.2.3泛函的變分與尤拉方程泛函變分問題的經(jīng)典解法有兩種，一種稱之為直接解法，另一類是間接解法。直接解法是直接把泛函的極值問題近似地轉(zhuǎn)化為一般多元函數(shù)的極值問題，用有限維子空間中的函數(shù)去逼近無窮維空間中的極值函數(shù)，從而近似求得泛函的極值。間接解法是將變分問題轉(zhuǎn)化為尤拉方程（微分方程）的定解問題，即邊值問題來求解。以式（3.11）這種最簡形式來推導尤拉方程。設函數(shù)y(x)稍有變化，記作yy，y稱之為y(x)的變分，它反映了整個函數(shù)的變化量。這樣泛函J[y]的值也應隨之變動，相應于變分y的泛函增量為JJ[yy]J[y]x2[F(x,yy,yy)F(x,y,y)]dx（3.12）x1將（3.12）式由多元函數(shù)的泰勒公式展開J{[FyFy]1[2F(y)222Fyy2F(y)2]}dx y y 2y2 yyy2J2J3J式中作為泛函增量J的線性主部為（3.13）Jx2[FyFy]dxxy y（3.14）1J稱為泛函J[y]的一次變分（簡稱變分）。而2J、3J……分別是函數(shù)變分y及其導數(shù)y的二次、三次齊次式……等的積分，依次稱為二次變分，三次變分……令變分問題的解為yy(x)，且設極值解yy(x)稍有變動yy，令y(x)（3.15）式中為任意給定的微量實參數(shù)，值就確定了yy(x,)函數(shù)族中的某一曲線，進而確定泛函J[y(x,)]之值；而(x)是定義于區(qū)間[x,x]且滿足 1 2xx0齊次邊界條件的可微函數(shù)。于是泛函Jy=Jyx,=1 2就成為變量的函數(shù)，且當0時獲極值函數(shù)的解。()在0時取得極值的必要條件是()0(0)x12FyyFyydx0（3.16）x利用分部積分，并根據(jù)變分與微分順序可互換原理，(3.16）式可寫為x2[Fd(F)]ydx0（3.17）xydxy1由于（3.17）對任意y均成立，故有FdF()0（3.18）ydxy方程（3.18）就稱為泛函（3.11）的極值問題的尤拉方程。綜上所述，有限元法的基本特點是：離散化過程保持了明顯的物理意義。這是因為，變分原理描述了支配物理現(xiàn)象的物理學中的最小作用原理（如力學中的最小勢能原理、靜電學中的湯姆遜定理等）。優(yōu)異的解題能力。與其他數(shù)值計算方法相比較，有限元法在適應場域邊界幾何形狀，以及媒質(zhì)物理性質(zhì)變異情況復雜的問題求解上，有突出的優(yōu)點。從數(shù)學理論意義上講，有限元法作為應用數(shù)學的一個重要分支，很少有其他方法應用的這樣廣泛。它使微分方程的解法和理論面目一新，推動了泛函分析與計算方法的發(fā)展。3.3矩量法矩量法，是近年來在天線、微波技術(shù)和電磁波散射等方面廣泛應用的一種方法。從這些實際問題涉及開域、激勵場源分布形態(tài)較為復雜等特征出發(fā)，矩量法是將待求的積分方程問題轉(zhuǎn)化為一個矩陣方程的問題[7]，借助于計算機，求得其數(shù)值解，從而在所得激勵源分布的數(shù)值解基礎上，即可算出輻射場的分布及其波阻抗等特性參數(shù)。3.3.1矩量法的基本原理先選定基函數(shù)對未知函數(shù)進行近似展開，代入算子方程，再選取適當?shù)臋?quán)函數(shù)，使在加權(quán)平均的意義下方程的余量等于0，由此將連續(xù)的算子方程轉(zhuǎn)換為代數(shù)方程。原則上，矩量法可用于求解微分方程和積分方程，但用于微分方程時所得到的代數(shù)方程組的系數(shù)矩陣往往是病態(tài)的，故在電磁場問題中主要用于求解積分方程。3.3.2加權(quán)余量法設給定邊值問題的場方程統(tǒng)一表述為如下的算子方程，即L(f)g（3.19）已知邊界條件為uur（3.20）ssb1ur（3.21）qns2sb其中:L是線性算子,f是待求函數(shù),g是已知的源。若u為精確解，則方程(3.19)和邊界條件(3.20),(3.21)應該完全滿足。但大多數(shù)情況下，不能得到u的精確解，只能通過數(shù)值方法進行估計。構(gòu)造一個由有限個線性無關(guān)函數(shù)N(i=1,2,…,n)所組成的基函數(shù)集合iN，借以展開待求函數(shù)u的近似解為u~nNuNTu（3.22）ii1將u~代入式(3.19)中必然存在誤差，即Ru~Lu~g(3.23)取一個歸屬于試探函數(shù)的權(quán)函數(shù)集合W，令WLu~gdV0(j=1,2,…,n)（3.24）v式(3.24)由n個方程構(gòu)成的方程組，它等價于人為地強制近似解u~，使其因不能精確地滿足場方程而導致的誤差在平均的含義上等于零。按式(3.24)展開，所構(gòu)成的各種求解積分或微分方程近似解的方法可被統(tǒng)稱為加權(quán)余量法[8]。因為按給定權(quán)函數(shù)W展開式的式(3.24)，即意味著余量RLu~g對W取矩的一組j j平衡式，故式(3.24)的構(gòu)造亦就被稱為矩量法?；诩訖?quán)余量式(3.24)，進行移項處理，便得WLu~dVWgdV(j=1,2,…,n)（3.25） j j v v將式(3.22)代入式(3.25)的左端，有WLnNunuWLNdV（3.26） vji1 iii1ivj i為了書寫方便，令nuWLNdVW,LN和WgdVW,gijijijjvvi1代入(3.26)式，則可寫成nuWLNdVi j ivi1W,gj1,2,n（3.27）j這樣，即展開成含n個未知數(shù)u的n個方程。若用矩陣形式表示，則有ilug（3.28）綜上所述，矩量法的特點是：矩量法將連續(xù)方程離散化為代數(shù)方程組,既適用于求解微分方程,又適用于求解積分方程。它的求解過程簡單,求解步驟統(tǒng)一,應用起來比較方便，然而需要一定的數(shù)學技巧,如離散化的程度、基函數(shù)與權(quán)函數(shù)的選取,矩陣求解過程等。另外必須指出的是,矩量法可以達到所需要的精確度、解析部分簡單,可計算量很大,即使用高速大容量計算機,計算任務也很繁瑣。4電磁場數(shù)值計算的應用前景電磁場數(shù)值計算方法近六十年來發(fā)展如此之快，除由于從事這方面的科研人員的努力之外，主要是其研究成果迅速被電機、電器、變壓器、加速器、微波器件、計算機磁頭等領域采用[9]，對改善產(chǎn)品性能、降低生產(chǎn)成本，起著越來越大的作用。眾所周知，任何電磁器件，包括國民經(jīng)濟中應用極其廣泛的電機與變壓器，其能量轉(zhuǎn)換都是通過電磁場來實現(xiàn)的。但傳統(tǒng)的電磁產(chǎn)品的設計方式由于客觀條件與手段的限制，把場的實際分布參數(shù)當作集中參數(shù)處理，不可避免地帶來相當大的誤差。在迫不得已的情況下，只能用模擬、實驗等方法處理，其耗費大、周期長，可借鑒的經(jīng)驗不多，而實驗的模擬也不都是有條件的，此外，工農(nóng)業(yè)和日常生活中所用電磁裝置越來越多，生產(chǎn)和銷售競爭劇烈，因此有效的設計方法顯然受到重視。在采用電磁場數(shù)值計算以前的任何方法，即使是十分精巧的代數(shù)解析方法，也只能適用于特別簡單的幾何結(jié)構(gòu)以及一些特殊的假定模型，有時模型甚至簡化到不能容忍的程度，但除此之外別無它法。因為實際的電磁場問題，特別是電機電磁場，其邊界情況十分復雜，加上鐵的飽和以及導體中的渦流效應，解析方法是無能為力的。而數(shù)值方法可以模擬復雜的形狀，可以適應非線性問題，可以進行渦流的分布計算。電磁場的數(shù)值計算與流體力學分析、溫度分析、機械應力分析、電磁力分析和生產(chǎn)計劃等等有機聯(lián)系起來，由計算機完成全部設計，構(gòu)成所謂CAE[10]系統(tǒng)(ComputerAidedEngineering)。CAE系統(tǒng)基本可包括設計與制造的全部過程，產(chǎn)品的設計、制造方案、準備零件、草圖、計算、成本核算、生產(chǎn)、機床數(shù)控和試驗、模擬和產(chǎn)品的自動測試都可能列入到CAE系統(tǒng)中。目前CAE系統(tǒng)在西方國家發(fā)展很快，估計每年增長40%，尤其在機械行業(yè)中，形成了所謂CADMAT即計算機輔助設計、加工與試驗。從長遠觀點看，電磁場數(shù)值計算要發(fā)揮更實際的作用，必須與其他多種學科相結(jié)合以形成CAE系統(tǒng)，CAE的發(fā)展必然給工業(yè)結(jié)構(gòu)上帶來巨大的變化。結(jié)束語數(shù)值計算是一門計算的藝術(shù)，電磁場的數(shù)值計算橫跨了多個學科，是數(shù)學理論、電磁理論和計算機應用能力的完美組合。通過本次論文設計，我對差分、變分、泛函、加權(quán)余量法等數(shù)學知識有了深刻的認識，掌握了有限差分法、有限元法、矩量法的基本原理，并能進行簡單的運算。但若用數(shù)值計算方法解決復雜的電磁場問題，則在很大程度上依賴于數(shù)學知識及計算機編程能力，這就需要進一步系統(tǒng)學習相關(guān)的理論知識。參考文獻文舸一．計算電磁學的進展與展望[J]．電子學報，1995，23(10)：62-69．劉圣民．電磁場的數(shù)值方法[M]．武漢：華中理工大學出版社，1991：23-45．鄧東皋，彭立中．小波分析[J]．數(shù)學進展，1991，20(3)：294-310．洪偉．計算電磁學研究進展［J］．東南大學學報（自然科學版），2002，32（3）：335-339．倪光正，楊仕友，錢秀英等．工程電磁場數(shù)值計算[M]．北京：機械工業(yè)出版社，2003：123-161．盛劍霓．工程電磁場數(shù)值分析[M]．西安：西安交通大學出版社，1991：69-75．方靜，汪文秉．有限元法和矩量法結(jié)合分析背腔天線的輻射特性［J］．微波學報，2000，16(2)：139-143．連漢雄．電磁場理論的數(shù)學方法[M]．北京：北京理工大學出版社，1990：10-15．張平文，劉法啟，張宇．小波函數(shù)值的計算[J]．計算數(shù)學，1995(3)：173-185．樓仁海，符果行，袁敬閎．電磁理論［M］．成都：電子科技大學出版社，1996：73-102．ElectromagneticnumericalmethodAbstract：Numericalcalculation,whichiswidelyusedinmanyfields,suchaselectric,militaryaffairs,economy,ecology,medicaltreatment,astronomy,geologyandsoonis