MATLAB靜電場電場電勢

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

2、b在圖像方面問題處理的應(yīng)用,問題一：等量同號點(diǎn)電荷的電場線的繪制,問題分析：根據(jù)庫侖定律：在真空中，兩個靜止的點(diǎn)電荷之間的作用力與這兩個電荷的電量乘積成正比，與它們之間距離的平方成反比，作用力的方向在兩個電荷的連線上，兩電荷同號為斥力，異號為吸力，他們之間的力F滿足： 由電場強(qiáng)度E的定義可知,對于點(diǎn)電荷，根據(jù)場論基礎(chǔ)中的定義，有勢場E的勢函數(shù)為： 在matlab中，由以上公式算出各點(diǎn)的電勢U，電場強(qiáng)度E后，可以用matlab自帶的庫函數(shù)繪出一對點(diǎn)電荷的電場線和等勢線,源程序,q=0.5; xm=2.5; ym=2; x=linspace(-xm,xm);%橫坐標(biāo)向量 y=linspace(-y

3、m,ym); %縱坐標(biāo)向量 X,Y=meshgrid(x,y); %產(chǎn)生自變量網(wǎng)絡(luò)坐標(biāo) R1=sqrt(X+1).2+Y.2);%第一個正電荷到場點(diǎn)的距離 R2=sqrt(X-1).2+Y.2);%第二個正電荷到場點(diǎn)的距離 U=1./R1+q./R2;%計算電勢 u=1:0.5:4;%等勢線的電場向量 figure,contour(X,Y,U,u)%畫等勢線 gridon%有網(wǎng)格 holdon%設(shè)置圖形保持狀態(tài) plot(-xm;xm,0;0)%畫水平線 plot(0;0,-ym;ym)%畫豎直線 plot(-1,0,o,MarkerSize,12) plot(1,0,o,MarkerSize

4、,12) Ex,Ey=gradient(-U,x(2)-x(1),y(2)-y(1);%用等勢梯度求場強(qiáng)的兩個分量 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

5、; 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)%畫右下電場線 axisequaltight%縱橫坐標(biāo)軸采用等長刻度 title%題目 xlabel(x,fontsize,16)%X軸說明 ylabel(y,fontsize,16)%Y軸說明 text(-xm,-ym+0.5,txt,fontsize,16,結(jié)論,可以清楚的看出同號電荷相互排斥的現(xiàn)象，各電荷的等勢線呈圈狀圍在電荷周圍，且越靠近電荷電場

6、線分布越密集,問題二：等量異種電荷的電勢，電場的分布,物理理論：平面上在x=2,y=0處有一正電荷，x=-2,y=0處有一負(fù)電荷 根據(jù),源程序,x,y=meshgrid(-5:0.2:5,-4:0.2:4);%建立數(shù)據(jù)網(wǎng)格 z=1./sqrt(x-2).2+y.2+0.01)-1./sqrt(x+2).2+y.2+0.01);%電勢的表達(dá)式 mesh(x,y,z)%三維曲面繪圖,可以看出，在不同位置的兩正負(fù)電荷由于所帶電性的不同電勢分布也大不相同，也可以說是在某些角度看來是完全相反的,電場分布,x,y =meshgrid(-2:0.1:2,-2:0.1:2); %以0.1為步長建立平面數(shù)據(jù)網(wǎng)格

7、.x,y的取值范圍都是-2,2 z=1./sqrt(x-1).2+y.2+0.01)-1./sqrt(x+1).2+y.2+0.01); %電勢公式 px,py=gradient(z); %求電勢在x,y方向的梯度即電場強(qiáng)度 contour(x,y,z,-12,-8,-5,-3,-1,-0.5,-0.1,0.1,0.5,1,3,5,8,12) %畫出等勢線 hold on %保留上一次作圖痕跡 quiver(x,y,px,py, k) %畫出各點(diǎn)上電場的大小和方向,k為圖形參數(shù)控制%quiver是繪制點(diǎn)x,y處的矢量px,py,即畫出各點(diǎn)電場的大小和方向,問題三 ：帶電粒子在電磁場中的運(yùn)動,如

8、果一個帶電粒子在既有電場又有磁場的區(qū)域里運(yùn)動，則其受到的電磁力為 這里，我們進(jìn)一步討論既有電場又有磁場情況下帶電粒子的運(yùn)動情況,題目,設(shè)質(zhì)量為m，帶電量為q的粒子在磁感應(yīng)強(qiáng)度為B，電場強(qiáng)度為E的電磁場中運(yùn)動，建立和求解該問題的微分方程，并分E0，B0；E=0，B0和E0，B=0三種情況畫出帶電粒子在電磁場中的運(yùn)動軌跡,源程序,syms w x y z t B E m q; %定義變量 E=input(E=);B=input(B=); %輸入E和B值 x,y,z=dsolve(D2x=q*B/m*Dy,D2y=q*E/m-q*B/m*Dx,D2z=0,x(0)=0,y(0)=0,z(0)=0,D

9、x(0)=0.01,Dy(0)=6,Dz(0)=0.01) ; %初始條件取x(0)=y(0)=z(0)=0,Dx(0)=0.01,Dy(0)=6,Dz(0)=0.01 q=1.6e-2; m=0.02; X=subs(x y z); x=X(1),y=X(2),z=X(3), %賦值函數(shù) ezplot3(X(1),X(2),X(3) 運(yùn)行上述程序，例如，取E=4, B=8可得下列特解并給出圖（a）。 x =-15/16*cos(32/5*t)-49/640*sin(32/5*t)+1/2*t+15/16 y =15/16*sin(32/5*t)-49/640*cos(32/5*t)+49/640 z =1/100*t,感想與總結(jié),經(jīng)過一學(xué)期一學(xué)期MATLAB課程學(xué)習(xí)，我們對一些數(shù)學(xué)題目有了新的認(rèn)識。就像這次我們做的有關(guān)電場分布以及帶電粒子在電磁場中運(yùn)動的問題，通過編寫程序我們可以直觀地看到電荷周圍的電場分布，以及帶電粒子在電磁場中的運(yùn)動軌跡。讓問題變得簡單化，可視化，更加易于理解，也加