基于matlab的電磁場仿真與分析

?

基于matlab的電磁場仿真與分析

宋應龍

（西安電子科技大學,陜西西安,710126）

摘要：電磁場問題往往較為抽象，并且計算也相對比較復雜。為此，借助計算機軟件對復雜問題進行求解的重要性也得到凸顯。本文通過兩個例子，著重介紹了matlab在電磁場仿真與數值分析中的應用。結果表明這一方法具有計算效率高，能使抽象的場的問題具體化的優點。

關鍵詞：Matlab；電磁場；仿真；數值分析

0　引言

電磁場問題一直是人類科學史上不斷探索和深究的問題，雖然Maxwell方程用簡潔的數學語言描述了電磁場所遵循的普遍規律，但其具體問題的求解往往是復雜的，并且很多時候往往不能直接精確求解，需要進行大量的近似計算來獲取結果。除此之外，對于數學表達式形式的求解結果，人們又往往難以形象地理解電磁場具體的空間分布。考慮到Matlab是功能強大的科學與工程計算軟件，擁有強大的數值分析和繪圖能力，被廣泛應用在信號處理、系統仿真、自動控制、小波分析以及圖形圖像分析等領域。我個人認為，Matlab的廣泛應用不僅得益于它良好的人機交互界面和簡潔易學的編程語言，更主要的是其函數庫的強大以及良好的擴充能力。為了解決實際電磁場問題求解時所面臨的窘境，Matlab的使用顯得尤為重要。其中，Matlab PDE工具箱為我們提供了方便。PDE即為偏微分方程（partial differential equation）的英文縮寫。該工具箱主要運用有限元法對偏微分方程來求解，方便的解決了電磁場偏微分方程難于求解以及求解過程較為繁瑣的問題。有限元法是將由偏微分方程表征的連續函數所在的封閉場域劃分成有限個小區域，每一個小區域用一個選定的近似函數來代替，所以整個場域上的函數被離散化，由此獲得一組近似的代數方程，聯立求解后即為該場域中函數的近似解。本文嘗試利用Matlab PDE工具箱來解決一些切實遇到的電磁場問題，進一步使得求解電磁場問題的方式更加多元化。

1　電磁場的邊值問題

我們知道，通過微分方程及相關邊界條件描述的問題稱為邊值問題。根據不同的邊界條件，邊值問題又可以分為三類：Dirichlet問題，Neumann問題以及混合問題。對于一些結構特殊的靜態場問題，我們往往根據Maxwell方程組結合相應的邊界條件來具體求解場的分布。

2　電荷密度為零的圓柱面內部的電位分布

圖1　圓柱面內部的電位分布

圖2　銅線橫截面電流

考慮如下問題：一半徑為a（任意常數）的圓柱面上，電位分布為：當時，；當時，。則其內部的電位分布用matlab PDE來計算仿真，首先在Matlab命令窗口輸入pdetool,啟動PDE工具箱。然后可以在Options菜單中調整坐標紙的柵欄格的稀疏。然后在Draw菜單下選擇Draw Mode來建模，得到圓形區域。接著進入Boundary Mode，鼠標雙擊邊界來根據已知設置邊界條件，都選擇Dirichlet條件，其中，上半圓的r=15,h=1,下半圓的r=0,h=1。然后在PDE Spacification中設置方程參數，選擇Elliptic(橢圓形)，c=1;a=0;f=0;點擊進行網格剖分，最后點擊，得到結果如圖1所示。圖1中箭頭方向即為內部電場的方向，與箭頭正交的線即為等電位線。另外，可以在Plot菜單下的Parameter中設置所畫圖形的參數，如：有無等勢線，有無矢量方向箭頭等等。

由（*）式可知，內部的電位分布滿足拉普拉斯方程，如果我們采用圓柱坐標系的分離變量法求解的話，可以得到內部的電位是：

其中r為圓柱面內部一點到圓柱面圓心的距離。從圖中可以看出，等位線對于該圓柱面與側放的地球儀上的經線相似。將仿真結果和數學表達式對比后可以發現，通過數學表達式我們不大容易聯想到內部具體的電位分布狀況，此時matlab的建模仿真對于將問題的解形象化發揮了重要作用，并且將繁瑣的計算求解過程交給了計算機程序，使得計算效率提高。

3　交流電磁場的計算仿真

又有問題如下：設在一段銅線上加有50MHz的交流電，其表面的電流密度，由可知，表面電場的Dirichlet條件為。在matlab PDE工具箱中選擇AC Power Electromagnetics模型。又知銅的電導率為，磁導率，則在Specify Boundary conditions的對話框中設置Dirichlet條件為,在PDE Specification對話框中選擇Elliptic,設置相應的電磁參數及頻率:Omega=2*pi*5E7,mu=4*pi1E-7,sigma=5.7E7,epsilon=8.8E-12。然后剖分網格，如果要進一步提高精度，可以點擊來加密網格。通過計算仿真得到內部電流密度的分布結果如圖2～3所示。

Matlab繪制的立體圖可以用鼠標拖動向各個方向旋轉，進而可以從不同角度來觀察分布的情況。從圖2～3也可以看出，頻率較高的情況下，電流主要分布在銅線表面，即為集膚效應。

假如我們通過波動方程來求其解析解的話，除了繁瑣的過程外，其結果將含有Bessel函數，這對于理解和想象其電流分布又

圖3　銅線橫截面電流密度分布立體圖

增加了困難。

4　結束語

綜上所述，本文主要研究了Matlab PDE工具箱在電磁場問題求解中的應用，以及對求解結果的分析等。通過上面兩個例子可以看到，matlab在求解關于電磁場的位分布，電流密度分布等的問題中具有高效、形象直觀的優點，并且借助其繪制的圖形，便于更深刻地理解相應的電磁場理論和解決實際的電磁場問題。該方法主要要求是對具體問題的建模，然后借助計算機程序來求解，相對較為簡潔，便于理解和掌握，對于在校學生、工程技術人員甚至科研人員都不失為一種良好的電磁場分析手段。

參考文獻

[1] 王沫然.Matlab與科學計算[M].電子工業出版社.2005.

[2] 尚濤，石端偉等.工程計算可視化與matlab實現[M].武漢大學出版社.2002.

[3] 金建銘.電磁場有限元法[M].西安電子科技大學出版社.1998.

[4] 趙曉峰.基于matlab的電磁輻射數據分析與仿真[D].上海：復旦大學，2008.

[5] 路宏敏，趙永久，朱滿座.電磁場與電磁波基礎[M].科學出版社.2011.

[6] 黃作英，闕沛文.基于matlab的電磁場數值分析[J].計算機工程與應用.2004,(36):196-198.

宋應龍 （1994- ），男，西安電子科技大學電子工程學院，曾連續多次獲得獎學金，電磁場與無線技術專業。

Electromagnetic field simulation and analysis based on matlab

Song Yinglong

(Xidian University,Xi’an,710126,China)

Abstract：Electromagnetic field problems are usually abstract,and the computation to it is also relatively complex.Therefore,the ability of solving the problems with the help of computers is becoming more important.This article introduces the application of matlab on simulation and numerical analysis.And the results shows that this method is efficient,and can make the abstract problems easy to handle.

Keywords：Matlab;Electromagnetic field;Simulation;Numerical Analysis

作者簡介

中圖分類號：TP391.9

文獻標識碼:A