Matlab在電磁場可視化教學(xué)和數(shù)值計(jì)算方面的應(yīng)用研究

1、MATLAB 課程設(shè)計(jì)報(bào)告Matlab 在電磁場與電磁波數(shù)值計(jì)算和可視化教學(xué)方面應(yīng)題 目用初探學(xué) 院專 業(yè) 電子信息工程學(xué)院 電子信息工程姓名和學(xué)號 劉偉 P11314008指導(dǎo)教師程鴻Matlab 在電磁場與電磁波數(shù)值計(jì)算和可視化教學(xué)方面應(yīng)用初探劉偉電子信息工程 P11314008摘要：本文利用Matlab 實(shí)現(xiàn)了SOR 超松弛法，繪制了單個(gè)點(diǎn)電荷的平面電場線與等勢線，和等量同號點(diǎn)電荷的電場線和等勢線，并利用GUI 對平面簡諧波進(jìn)行了仿真，展現(xiàn)了Matlab 在電磁場與電磁波數(shù)值計(jì)算和可視化教學(xué)方面應(yīng)用前景。關(guān)鍵詞：Matlab ；電磁場與電磁波；數(shù)值計(jì)算；可視化；GUI ；仿真1 引言電磁

2、場與電磁波是電子信息類專業(yè)的一門專業(yè)核心課。Maxwell 用簡潔的方程描繪了描述了電磁場所遵循的普遍規(guī)律，但其具體問題的求解往往是復(fù)雜的，并且很多時(shí)候往往不能直接精確求解，需要進(jìn)行大量的近似計(jì)算來獲取結(jié)果。除此之外，對于數(shù)學(xué)表達(dá)式形式的求解結(jié)果，人們又往往難以形象地理解電磁場具體的空間分布。1而Matlab 由于其強(qiáng)大的功能、簡單易學(xué)的編程語言和可視化的仿真環(huán)境，則為電磁場與電磁波的計(jì)算和可視化提供了一定條件。2 Matlab簡介Matlab 是美國Mathworks 公司出品的商業(yè)數(shù)學(xué)軟件，它將數(shù)值分析、矩陣計(jì)算、科學(xué)數(shù)據(jù)可視化、非線性動態(tài)系統(tǒng)的建模和仿真等諸多強(qiáng)大功能集成在一個(gè)易于使用的

3、視窗環(huán)境中，為科學(xué)研究、工程設(shè)計(jì)以及必須進(jìn)行有效數(shù)值計(jì)算的眾多科學(xué)領(lǐng)域提供了一種全面的解決方案，并在很大程度上擺脫了傳統(tǒng)非交互式程序設(shè)計(jì)語言（如C 、Fortran ）的編輯模式，代表了當(dāng)今國際科學(xué)計(jì)算軟件的先進(jìn)水平。3 Matlab在電磁場數(shù)值計(jì)算方面的應(yīng)用實(shí)例Matlab 在電磁場數(shù)值計(jì)算方法上的應(yīng)用主要包含了有限差分法、矩量法、有限元法、時(shí)域有限差分法等，通過利用軟件計(jì)算相同的問題，可以驗(yàn)證結(jié)果并加深對算法的理解2 我們以有限差分法作為例子。實(shí)際問題中，一般借助于數(shù)值法求電磁場的數(shù)值解3。然而在求解時(shí)，由于節(jié)點(diǎn)個(gè)數(shù)很多，聯(lián)立差分方程的個(gè)數(shù)往往可達(dá)幾百甚至幾千個(gè)4，因此我們需要采用逐次近似

4、的迭代方法解決問題。最常用的迭代法包括Gauss-Seidel 迭代法和SOR 迭代法，由于Gauss-Sidel 迭代法是SOR 迭代法的一種特殊形式，我們以后者作為代表。i -1n (k +1w (k (k +1(k x =x +b -a x -a x i i ij j ij j i a ii j =1j =i （1-1） (k (k (k (k T (k =0, 1, 2.; i =1, 2,. n x =x 1, x 2,., x n (其中參數(shù)稱作松弛因子，顯然，當(dāng)時(shí)，解式（1-1）的SOR 方法就是Gauss-Seidel 迭代法。5實(shí)現(xiàn)SOR 迭代法的自定義函數(shù)sor 在附件中。

5、4 Matlab在電磁場和電磁波教學(xué)可視化實(shí)例4.1Matlab 與繪圖有關(guān)的函數(shù)(1 plot語句：二維圖形繪圖命令中最基本的指令就是plot 。如果輸入兩個(gè)矢量 x 、y 則plot(x ,y產(chǎn)生的是y 相對于x 的圖形。(2 fplot語句：當(dāng)已知函數(shù)關(guān)系時(shí)，采用fplot 語句可以更加快速、精確地繪制出指定區(qū)間上的函數(shù)圖像。(1plot3(x,y,z,s語句：X 、Y 、Z 是同維數(shù)的矢量或矩陣。每組X 、Y 、Z 構(gòu)成一個(gè)點(diǎn)的坐標(biāo)，各點(diǎn)依次相連，形成一條曲線。如果是矩陣，則它們相應(yīng)的列構(gòu)成一條三維曲線的數(shù)據(jù)點(diǎn)坐標(biāo)，所以用矩陣可以同時(shí)畫多條空間曲線。(2網(wǎng)線和表面圖：Matlab 中，

6、曲面是用xy 平面的各個(gè)格點(diǎn)上的z 坐標(biāo)來定義，相鄰點(diǎn)用直線連接。因此，為了描繪二元函數(shù)形成的空間曲面，首先要建立平面的數(shù)據(jù)網(wǎng)格。即根據(jù)x 、y 的定義域,，生成平面上的坐標(biāo)格點(diǎn), 每個(gè)格點(diǎn)對應(yīng)一組 x，y 。生成數(shù)據(jù)網(wǎng)格的命令是meshgrid ，格式為 =meshgrid(x，y, 使用 meshgrid 命令，將兩個(gè)獨(dú)立的矢量 x 、y 組成數(shù)據(jù)網(wǎng)格。mesh 是三維網(wǎng)格作圖命令，mesh(x,y,z畫出了每一個(gè)格點(diǎn)(x, y上對應(yīng)的z 值。surfl(X ,Y，Z 是畫帶光照模式的三維曲面圖 。 該命令顯示一個(gè)帶陰影的曲面, 融合了周圍的、散射的和鏡面反射的光照模式。要獲得較平滑的顏色

7、過度，應(yīng)該使用有線性強(qiáng)度變化的色圖，如 gray 、copper 、bone 、pink 等。64.2單個(gè)點(diǎn)電荷的平面電場線與等勢線點(diǎn)電荷是帶電體的一種理想模型。如果在研究的問題中，帶電體的形狀、大小以及電荷分布可以忽略不計(jì)，即可將它看作是一個(gè)幾何點(diǎn)，則這樣的帶電體就是點(diǎn)電荷。一個(gè)實(shí)際的帶電體能否看作點(diǎn)電荷，不僅和帶電體本身有關(guān)，還取決于問題的性質(zhì)和精度的要求。與質(zhì)點(diǎn)、剛體等概念一樣，點(diǎn)電荷是實(shí)際帶電體的抽象和近似，它是建立具有普遍意義的基本規(guī)律的不可或缺的理想模型，又是把復(fù)雜多樣的實(shí)際問題轉(zhuǎn)化或分解為基本問題時(shí)必不可少的分析手段。設(shè)有一個(gè)帶有+q電荷的點(diǎn)電荷，則該點(diǎn)電荷在P(x，y 處產(chǎn)生的

8、電勢為：（2-1）由電場強(qiáng)度可得, 在xOy 平面上, 電場強(qiáng)度的公式為:（2-2）（2-3）在Matlab 中, 應(yīng)用gradient 函數(shù)可輕松求出電場強(qiáng)度。為了簡單起見，我們令，演示結(jié)果如圖所示。單個(gè)點(diǎn)電荷的平面電場線與等勢線10.8 0.60.40.2E (U 0-0.2-0.4-0.6-0.8r4.3 等量同號點(diǎn)電荷的電場線和等勢線設(shè)有兩個(gè)帶有等量同號電荷+q的兩個(gè)點(diǎn)電荷，根據(jù)疊加原理，我們可以分別利用（1-1）、（1-2）、（1-3）求出xOy 平面中任意一點(diǎn)的場強(qiáng)。不妨設(shè)Q(x,y為平面中任意一點(diǎn)，距電荷為R 1，距電荷為R 2，則有（3-1）（3-2）（3-3）同樣地，我們令，

9、分別代入x1,y 1、x2,y 2，演示結(jié)果如圖所示。等量同號點(diǎn)電荷的電場線和等勢線2 1.510.5E (U 0-0.5-1-1.5r4.4 基于GUI 的簡諧波仿真用戶界面（或接口）是指：人與機(jī)器（或程序）之間交互作用的工具和方法。如鍵盤、鼠標(biāo)、跟蹤球、話筒都可成為與計(jì)算機(jī)交換信息的接口。圖形用戶界面(Graphical User Interfaces,GUI則是由窗口、光標(biāo)、按鍵、菜單、文字說明等對象(Objects構(gòu)成的一個(gè)用戶界面。用戶通過一定的方法（如鼠標(biāo)或鍵盤）選擇、激活這些圖形對象，使計(jì)算機(jī)產(chǎn)生某種動作或變化，比如實(shí)現(xiàn)計(jì)算、繪圖等。Matlab 為表現(xiàn)其基本功能而設(shè)計(jì)的演示程序

10、demo 是使用圖形界面的最好范例。使用Matlab 時(shí)，在指令窗中運(yùn)行demo 打開圖形界面后，只要用鼠標(biāo)進(jìn)行選擇和點(diǎn)擊，就可瀏覽那豐富多彩的內(nèi)容。一些基本的操作，簡單地制作了對簡諧波的仿真。當(dāng)然，我對 GUI 的操作和熟練還相距甚 遠(yuǎn)。 Matlab 是一個(gè)功能很強(qiáng)大的軟件，還有更多的功能等待發(fā)掘。 參考文獻(xiàn) 1 宋應(yīng)龍.基于 Matlab 的電磁場仿真與分析 J.電子測試，2015(9 2 Dikshitulu K.Kalluri 著，馬西奎等譯. 電磁場與波電磁材料及 MATLAB 計(jì)算 M.機(jī) 械工業(yè)出版社，2014. 3 孫玉發(fā)等.電磁場與電磁波M.合肥工業(yè)大學(xué)出版社，2014 4 馮慈璋 馬西奎