(論文)matlab在電磁場中的應(yīng)用

Matlab在電場中的應(yīng)用

[摘要]Matlab是一種用于算法開發(fā)，數(shù)據(jù)可視化，數(shù)值分析及數(shù)值圖形生成的高級工具語言，它主要被應(yīng)用于信號和圖像處理，通訊，控制系統(tǒng)設(shè)計(jì)，測試和測量等廣泛領(lǐng)域�；�于MATLAB強(qiáng)大的繪圖和仿真功能,對物理學(xué)中電磁學(xué)中的等量同號點(diǎn)電荷的電場線的繪制和帶電粒子在均勻電磁場中的運(yùn)動等問題進(jìn)行仿真, 來簡便、直觀、高效分析物理問題。在本文中，用Matlab的功能使靜電場里的某些模型(電場強(qiáng)度電勢、電場線、等勢線、等勢面)可視化，方便了我們對有關(guān)靜電場的知識的學(xué)習(xí)，提高了我們對知識的理解和運(yùn)用能力，本文主要是從圖像處理功能方面介紹了Matlab語言在靜電場一些問題中的應(yīng)用。

【關(guān)鍵字】：

一、引言

靜電場中的電場線，等勢線，等勢面等圖形是一種抽象的模型，在現(xiàn)實(shí)世界不具可視化的空間場的物體。所以，形象的模擬出以上問題的圖形，對于更進(jìn)一步學(xué)習(xí)與研究電場知識有很大的意義。靜電場的問題學(xué)習(xí)與理解起來具有一定的特殊性：它既有理論數(shù)值的計(jì)算，又有圖形圖像的輔助處理與理解。例如：形象的模擬出電場線，等勢線，等勢面，這能在教學(xué)中解決教師的授課難題，又能解決學(xué)生的理解上的困難。近年來，一直有人在不斷的探索這方面的問題，并且取得一定的成績。但還存在一定的缺陷，而Matlab恰好解決了這些問題！這使得這些抽象問題能有一門精確的工具軟件來處理完成。這正是Matlab在圖像方面問題處理的應(yīng)用。

二、Matlab在電場中的應(yīng)用

（1）等量點(diǎn)電荷的電場線的繪制

根據(jù)庫侖定律：在真空中，兩個(gè)靜止的點(diǎn)電荷之間的作用力與這兩個(gè)電荷的電量乘積成正比，與它們之間距離的平方成反比，作用力的方向在兩個(gè)電荷的連線上，兩電荷同號維斥力，異號為吸力，他們之間的力F滿足：

QQ F?k1

22 （1） R

由電場強(qiáng)度E的定義可知：

E?kQ （2） R2

建立電場線的微分方程（二維情況）。 因?yàn)殡妶鲋腥我稽c(diǎn)的電場方向都沿該點(diǎn)電場線的切線方向，所以滿足：

dyEy? (3) dxEx

引入?yún)⒆兞縯得到： dydx??dt (4) ExEy

設(shè)二點(diǎn)電荷位于(-2,0)和(2,0)，二點(diǎn)電荷“電量”為q1和q2（均等于10）, 由庫倫定律和電場的疊加原理，得出下列微分方程： x?dx?Ex?q1?x?2?

3?q2?x?2?3 (5) dt??x?2?2?y222??????x?2?2?y2??

y?dyq1ydt?Ey?3?q2y3 (6)

??x2??2?2?y2????x?2?2?y22???

解此方程就可以繪制出電場線。

clc

clear

q=0.5;

xm=2.5;

ym=2;

x=linspace(-xm,xm); %橫坐標(biāo)向量

y=linspace(-ym,ym); %縱坐標(biāo)向量

[X,Y]=meshgrid(x,y);%產(chǎn)生自變量網(wǎng)絡(luò)坐標(biāo)

R1=sqrt((X+1).^2+Y.^2); %第一個(gè)正電荷到場點(diǎn)的距離

R2=sqrt((X-1).^2+Y.^2); %第二個(gè)正電荷到場點(diǎn)的距離

U=1./R1+q./R2; %計(jì)算電勢

u=1:0.5:4; %等勢線的電場向量

figure

contour(X,Y,U,u) %畫等勢線

grid on %有網(wǎng)格

hold o

大學(xué)網(wǎng) n %設(shè)置圖形保持狀態(tài)

plot([-xm;xm],[0;0]) %畫水平線

plot([0;0],[-ym;ym]) %畫豎直線

plot(-1,0,'o','MarkerSize',12)

plot(1,0,'o','MarkerSize',12)

[Ex,Ey]=gradient(-U,x(2)-x(1),y(2)-y(1)); %用等勢梯度求場強(qiáng)的兩個(gè)分量

dth1=20;

th1=(dth1:dth1:180-dth1)*pi/180; %電場線的起始角度

r0=0.1;

x1=r0*cos(th1)-1; %電場線的起點(diǎn)橫坐標(biāo)

y1=r0*sin(th1); %電場線的起點(diǎn)縱坐標(biāo)

streamline(X,Y,Ex,Ey,x1,y1) %畫左上電場線

streamline(X,-Y,Ex,-Ey,x1,-y1)%畫左下電場線

dth2=dth1/q;%右邊電場線角度間隔

th2=(180-dth2:-dth2:dth2)*pi/180;

x2=r0*cos(th2)+1; %電場線的起點(diǎn)橫坐標(biāo)

y2=r0*sin(th2); %電場線的起點(diǎn)縱坐標(biāo)

streamline(X,Y,Ex,Ey,x2,y2) %畫右下電場線

streamline(X,-Y,Ex,-Ey,x2,-y2) %畫右下電場線

axis equal tight %縱橫坐標(biāo)軸采用等長刻度

title %題目

xlabel('x','fontsize',16) %X軸說明

ylabel('y','fontsize',16) %Y軸說明

txt=['電荷比:\itQ\rm_2/\itQ\rm_1='num2str(q)];

text(-xm,-ym+0.5,txt,'fontsize',16)

（2）、線電荷產(chǎn)生的電位：

設(shè)電荷均勻分布在從z=-L到z=L,通過原點(diǎn)的線段上，其密度為q(單位C/m),求在xy平面上的電位分布。

點(diǎn)電荷產(chǎn)生的電位可表示為 V?Q 是一個(gè)標(biāo)量。其中r為電荷到測量點(diǎn)的距離。線電4?R?0

荷所產(chǎn)生的電位可用積分或疊加的方法來求。為此把線電荷分為N段，每段長為dL。每段上電荷為q*dL,看作集中在中點(diǎn)的點(diǎn)電荷，它產(chǎn)生的電位為 dV?qdL然后對全部電荷求和即可。 4?R?0

把xy平面分成網(wǎng)格，因?yàn)閤y平面上的電位僅取決于離原點(diǎn)的垂直距離R，所以可以省略一維，只取R為自變量。把R從0到10米分成Nr+1點(diǎn)，對每一點(diǎn)計(jì)算其電位。

clc

clear all;

L=input(‘線電荷長度L＝ ’);

N=input(‘分段數(shù)N＝ ’);

Nr=input(‘分段數(shù)Nr＝ ’);

q=input(‘電荷密度q= ’);

E0=8.85e-12;

C0=1/4/pi/E0;

L0=linspace(-L,L,N+1);

L1=L0(1:N);L2=L0(2:N+1);

Lm=(L1+L2)/2;dL=2*L/N;

R=linspace(0,10,Nr+1);

for k=1:Nr+1

Rk=sqrt(Lm.^2+R(k)^2);

Vk=C0*dL*q./Rk;

V(k)=sum(Vk);

end

[max(V),min(V)]

plot(R,V),grid

(3)、計(jì)算平面上N個(gè)電荷之間的庫倫引力 由庫侖定律：F?q1q2 24??0r

Fx?q1q2?x2?x1? 4??0r3

q1q2?y2?y1? 34??0rFy?

r??x2?x1?2?y2?y12

先輸入電荷的數(shù)目，各電荷的坐標(biāo)及電荷量，再選一個(gè)電荷，求其它電荷對它的作用力，疊加求合力。再選下一個(gè)電荷，依次類推。

clc

clear all;

disp('計(jì)算平面上n個(gè)點(diǎn)電荷之間的庫侖力的有關(guān)程序')

n=input('輸入電荷數(shù)目n=');

for ic = 1:n %輸入給定條件

fprintf('對電荷q%g\n',ic);%自動命名電荷q1,q2,?qn,分行顯示

rc = input('輸入電荷位置[x,y]（米）:');

x(ic) = rc(1); %電荷ic的x坐標(biāo)

y(ic) = rc(2); %電荷ic的y坐標(biāo)

q(ic) = input('輸入電荷量（庫侖）：');

end

E0 = 8.85e-12; %真空中的介電常數(shù)

C0 = 1/(4*pi*E0); %合并常數(shù)

for ic = 1:n %循環(huán)計(jì)每個(gè)電荷所受的力

Fx = 0.0;Fy = 0.0;

for jc = 1:n

if(ic ~= jc)

xij = x(ic)-x(jc);yij = y(ic)-y(jc);

Rij = sqrt(xij^2+yij^2);

Fx = Fx+C0*q(ic)*q(jc)*xij/Rij^3;

Fy = Fy+C0*q(ic)*q(jc)*yij/Rij^3;

F=sqrt(Fx^2+Fy^2);

end

end

fprintf('其它電荷作用在電荷q%g上的合力與分力分別為：\n',ic);

fprintf('F=%gN\n',F);

fprintf('Fx=%gN\n',Fx);

fprintf('Fy=%gN\n',Fy);

end

三、結(jié)論

從全文可以了解Matlab是一種有著強(qiáng)大的圖像繪制功能的可視化軟件。是一個(gè)為科學(xué)，教育，和工程數(shù)值計(jì)算等多方面設(shè)計(jì)的高級語言。它的特點(diǎn)是程序設(shè)計(jì)過程簡單，這使我們將我們將主要的，最難的工作還是放在工程問題上，這正是體現(xiàn)了它的優(yōu)越性。圖像處理僅僅是Matlab功能一個(gè)方面的體現(xiàn)，他還有更強(qiáng)大的功能值得我們?nèi)�學(xué)習(xí)與應(yīng)用。

通過Matlab在靜電場中的實(shí)際運(yùn)用，可以看出利用Matlab解決靜電場問題有以下優(yōu)點(diǎn)：

(1) 介入Matlab，使得靜電場中抽象的問題清晰與明朗化。這在教學(xué)與研究當(dāng)中是問題變得簡單。相比以往，學(xué)生更能形象的理解靜電場這方面的知識，將這一方法與工具應(yīng)用到中學(xué)教學(xué)中可以更好的幫助學(xué)生學(xué)好電學(xué)知識。

(2) 通過這一問題的處理。由當(dāng)初的學(xué)習(xí)電場到現(xiàn)在更深一步的探討，使得我對電場方面的知識有了更進(jìn)一步的學(xué)習(xí)與理解。Matlab還有很多方面的應(yīng)用，我們要繼續(xù)深入的學(xué)習(xí)Matlab來處理更多，更為復(fù)雜的問題。

參考文獻(xiàn)：

[1] 朱漢敏．MATLAB在靜電場教學(xué)中的優(yōu)越性[M]．中國期刊網(wǎng),2005年第3期．

[2] 劉衛(wèi)國．MATLAB程序設(shè)計(jì)與應(yīng)用[M]．北京：高等教育出版社，2006．

[3] 馬文蔚．物理學(xué)中冊（第四版）[M]．北京：高等教育出版社，1999．

[4] 陳鍾賢．計(jì)算物理學(xué)[M]．哈爾濱：哈爾濱工業(yè)大學(xué)出版社，2001．

[5] 張智星．MATLAB程序設(shè)計(jì)與應(yīng)用[M]．北京：清華大學(xué)出版社, 2002．

【(論文)matlab在電磁場中的應(yīng)用】相關(guān)文章：

Matlab在遙感圖像雜散光處理中的應(yīng)用04-30

Matlab在北京飛機(jī)增雨航跡分析中的應(yīng)用04-29

基于MATLAB的DCT變換在JPEG圖像壓縮中的應(yīng)用04-27

Matlab和Word的鏈接在水文計(jì)算中的應(yīng)用05-02

MATLAB回歸分析在測繪數(shù)據(jù)處理中應(yīng)用04-28

電磁場數(shù)值方法在光波導(dǎo)課程教學(xué)中的應(yīng)用探討04-29

VC++和Matlab混合編程在圖像處理中的應(yīng)用04-30

基于Matlab的數(shù)據(jù)處理方法在GPS高程擬合中的應(yīng)用05-01

基于MATLAB的BP神經(jīng)網(wǎng)絡(luò)在清潔生產(chǎn)審核中的應(yīng)用04-25

基于Matlab的數(shù)據(jù)處理方法在GPS高程擬合中的應(yīng)用04-30