電磁場-點電荷-電場線-電勢-MATLAB--仿真-中南大學(xué)(2024063506)

1、電磁場理論實驗一利用Matlab模擬點電荷電場的分布.實驗?zāi)康模? .熟悉單個點電荷及一對點電荷的電場分布情況；2 .學(xué)會使用Matlab進行數(shù)值計算，并繪出相應(yīng)的圖形；二.實驗原理:根據(jù)庫倫定律：在真空中，兩個靜止點電荷之間的作用力與這兩個電荷的電 量乘積成正比，與它們之間距離的平方成反比，作用力的方向在兩個電荷的連線上，兩電荷同號為斥力，異號為吸力，它們之間的力 F滿足：F k警RR2式 1由電場強度E的定義可知:式 2對于點電荷，根據(jù)場論根底中的定義，有勢場E的勢函數(shù)為式3式4kQRUU,電場強度E后，可以用Matlab在Matlab中，由以上公式算出各點的電勢 自帶的庫函數(shù)繪出相應(yīng)電荷

2、的電場分布情況三實驗內(nèi)容：1.單個點電荷點電荷的平面電力線和等勢線真空中點電荷的場強大小是 E=kq /rA2，其中k為靜電力恒量,q為電量,r為點電荷到場點Px,y的距離.電場呈球?qū)ΨQ分布，取電量q> 0,電力線是以電荷為 起點的射線簇.以無窮遠處為零勢點，點電荷的電勢為U=kq /r,當U取常數(shù)時，此 式就是等勢面方程 .等勢面是以電荷為中心以 r 為半徑的球面 .平面電力線的畫法在平面上,電力線是等角分布的射線簇,用MATLAB畫射線簇很簡單.取射 線的半徑為 ( 都取國際制單位 ) r0, 不同的角度用向量表示 ( 單位為弧度 ) th=linspace(0,2*pi,13).射

3、線簇的終點的直角坐標為：x,y=pol2cart(th,r0).插入 x 的 起始坐標x=x; *x.同樣插入y的起始坐標，y=y; *y, x和y都是二維數(shù)組，每 一列是一條射線的起始和終止坐標.用二維畫線命令plot(x,y)就畫出所有電力線.平面等勢線的畫法在過電荷的截面上 , 等勢線就是以電荷為中心的圓簇 , 用 MATLAB 畫等勢 線更加簡單 .靜電力常量為 k=9e9, 電量可取為 q=1e- 9; 最大的等勢線的半徑應(yīng) 該比射線的半徑小一點？ r0.其電勢為u0=k8q /r0.如果從外到里取7條等勢線，最 里面的等勢線的電勢是最外面的 3 倍, 那么各條線的電勢用向量表示為

4、： u=linspace(1,3,7)*u0.從-r0到r0取偶數(shù)個點，例如100個點，使最中心點的坐 標繞過 0, 各點的坐標可用向量表示 ： x=linspace(- r0,r0,100), 在直角坐標系中可 形成網(wǎng)格坐標：X,Y=meshgrid(x).各點到原點的距離為：r=sqrt(X.A2+Y.A2),在乘 方時, 乘方號前面要加點 , 表示對變量中的元素進行乘方計算 .各點的電勢為 U=k8q. /r, 在進行除法運算時 , 除號前面也要加點 , 同樣表示對變量中的元素進 行除法運算 .用等高線命令即可畫出等勢線 contour(X,Y,U,u), 在畫等勢線后一般 會把電力線擦

5、除 , 在畫等勢線之前插入如下命令 hold on 就行了.平面電力線和 等勢線如圖 1, 其中插入了標題等等 .越靠近點電荷的中心 , 電勢越高, 電場強度 越大, 電力線和等勢線也越密 .單個點電荷的電場線與等勢線0.15圖1點電荷的立體電力線和等勢面立體電力線的畫法先形成三維單位球面坐標，繞z軸一周有8條電力線X,Y,Z=sphere(8),每維 都是9X 9的網(wǎng)格矩陣,將X化為行向量,就形成各條電力線的終點x坐標 x=r0=X(:)',其他兩個坐標也可同樣形成終點坐標y=r0+Y(:)' , z=r0+Z(:)'.對x坐標插入原點x=x(zeros(size(x

6、),其他兩個坐標如下形成y=y(zeros(size(y), z=z(zeros(size(z),用三維畫線命令plot3(x,y,z),就畫出所有電力線.立體等勢面的畫法畫5條等勢面時,各面的電勢為u=linspace(1,3,5)+u0,各等勢面的半徑為 r=k6q. /u,其中第一個球面的半徑為rr=r(1).三維單位球面的坐標可由 X,Y,Z=sphere命令形成，每維都是21X 21的網(wǎng)格矩陣，由于外球會包圍內(nèi)球， 因此把球面的四分之一設(shè)為非數(shù),表示割去該局部Z(X<0&丫<0)=nan.用曲面命 令可畫出第一個曲面surf(rr6X,rr6Y,rr6Z),只要取

7、不同的半徑就能畫出不同的等 勢面.為了使等勢面好看，可設(shè)置一個顏色濃淡連續(xù)變化的命令 shadi ng in terp點 電荷的立體電力線和等勢面如圖2,旋轉(zhuǎn)圖片可從不同的角度觀察.正電荷電場線等勢面的三維圖形0.20.150.1 -0.05 .-0.05-0.1-0.15-0.20.20.10.20.1-0.1-0.1y圖2-0.2 -0.22 一對點電荷平面等勢線的畫法,電量不妨分仍然用MATLAB的等高線命令畫等勢線.對于正負兩個點電荷 別取q1=2e- 9,q2=- 1e- 9,正電荷在x軸正方,負電荷在x軸負方,它們到原點 的距離定為a=0.02;假設(shè)平面范圍為xx0=0.05,yy

8、0=0.04,兩個坐標向量分別 x=linspace(- xx0,xx0,20)和 y=linspace(- yy0,yy0,50).設(shè)置平面網(wǎng)格坐標為 X,Y=meshgrid(x),各點到兩電荷的距離分別為 r仁sqrt(X- a).A2+Y.A2)和 r2=sqrt(X+a).A2+Y.A2).各點的電勢為 U=k6q1. /r1+k6q2. /r2,取最高電勢為 u0=50, 最低電勢取其負值.在兩者之間取11個電勢向量u=linspace (u0,- u0,11),等高線 命令contour(X,Y,U,u,'k-')用黑實線,畫出等勢線如圖4所示,其中,左邊從里到

9、外的第6條包圍負電荷的等勢線為零勢線.平面電力線的畫法利用MATLAB的箭頭命令,可用各點的電場強度方向代替電力線.根據(jù)梯度可求各點的場強的兩個分量Ex,Ey=gradient(- U),合場強為E=sqrt(Ex.A2+Ey.A2).為了使箭頭等長,將場強Ex=Ex. /E,Ey=Ey. /E歸一化,用箭頭命令 quiver(X,Y,Ex,Ey)可標出各網(wǎng)點的電場強度的方向，異號點電荷對的場點方向如 圖3所示.為了畫出連續(xù)的電力線，先確定電力線的起點電荷的半徑可取為ro=O.OO2, 如圖4所示，假設(shè)第一條電力線的起始角為30度,其弧度為q=30+pi /180,起始 點到第一個點電荷的坐標

10、為 x仁r0+cos(q),y=r0+sin(q),到第二個點電荷的坐標只 有橫坐標x2=2+a+x1不同.用前面的方法可求出該點到兩個電荷之間的距離r1和r2,從而計算場強的兩個分量以及總場強Ex=q1+x1 /r1A3 +q2+x2 /r2A3,Ey=q1+y/r1A3+q2+y/r2A3, E=sqrt(Ex6Ex+Ey6Ey).下面只要用到場強分量與總場 強的比值，在計算場強分量時沒有乘以靜電力常量k.由于電力線的方向與場強的切線方向相同,取線段為s=0.0001 由此可求出終點的坐標為 x1=x1+s#Ex/E,y=y+s+Ey/E,從而計算x2.以終點為新的起點就能計算其他終點.

11、當終點出界時或者到達另一點電荷時，這個終點可作為最后終點.這種計算電力 線的方法稱為切線法.y 0對點電荷的電場分布圖0.050.040.030.020.01-0.01-0.02-0.03-0.04一對不相等的電荷的等勢線圖和電場線圖0.050.040.030.020.01-0.01-0.02-0.03-0.04-0.05-0.05,/2*? <-F- h*i雀' 時.Td !#、峯't十r*- -1KH.4 m廠 JV1 -Vxr«-i>-乂./%/*' v崔 L_ h -*"-w +- ： - - z、T1 q z、. x - *i&

12、quot;【-.-/9' - / 11 fl=一 A*-:.'-" Vu -.-!：i L1二 Ttjj 3(v 一11 1 1 1'、q .宰C' I-11 J.1 rr r -y * .'""Th' i.;耳M.1/ >"- . -K E=r r",-r*>1屮H >1/n/f1-廠薩* /*-r-_梓*/ A: X.-'“=-./ $ATl*- 甲;%, <直- 1%K ' *J.'-r 4.- -.Jt-a/、I、 >>1F Jf

13、r'l、 >(!氣 -0.04-0.03-0.02-0.0100.010.020.030.040.05-5點電荷電場分布的3-D圖100 一0.5-0.5-1.5-2-2-102 1.5局部 M-file;1. 點電荷的平面電力線和等勢線 %點電荷的平面電力線和等勢線 %平面電力線的畫法 q=1e-9;r0=0.12; th=linspace(0,2*pi,13); x,y=pol2cart(th,r0);x=x;0.1*x;y=y;0.1*y;plot(x,y);grid onhold on plot(0,0,'o','MarkerSize',1

14、2) xlabel('x','fontsize',16) ylabel('y','fontsize',16)title('單個點電荷的電場線與等勢線','fo ntsize',20)%平面等勢線的畫法k=9e9;r0=0.1;u0=k*q/r0;u=linspace(1,3,7)*u0;x=linspace(-r0,r0,100);X,Y=meshgrid(x);r=sqrt(X.A2+Y.A2);U=k*q./r;hold on;contour(X,Y,U,u)2);2. 一對電荷平面等勢線和電場

15、線圖%一對電荷平面等勢線和電場線clear all;clf;%平面等勢線的畫法q1=2e-9;q2=-1e-9;a=0.02;%到原點的距離xx0=0.05;yy0=0.04;k=9e9;x=linspace(-xx0,xx0,20);y=linspace(-yy0,yy0,50);X,Y=meshgrid(x);r11=sqrt(xx0/1.7-a)A2+(yy0/1.7)A r22=sqrt(xx0/1.7+a)A2+(yy0/1.7)A r1=sqrt(X-a).A2+Y.A2);%各點到點電荷的距離r2=sqrt(X+a).A2+Y.A2);U=k*q1./r1+k*q2./r2;%各

16、點的電勢u0=k*q1/r11+k*q2/r22; u=linspace(u0,-u0,11); %取21 個等勢向量contour(X,Y,U,u,'k-');hold ongrid on plot(a,0,'o','MarkerSize',12); plot(-a,0,'o','MarkerSize',12); xlabel('x','fontsize',16); ylabel('y','fontsize',16);%平面電力線的畫法Ex,Ey=g

17、radient(-U);E=sqrt(Ex.A2+Ey.A2);Ex=Ex./E;Ey=Ey./E;hold on;quiver(X,Y,Ex,Ey);title('對不相等的電荷的等勢線 圖和電場線圖 ','fontsize',20)clear;4.3. 立體電力線的畫法%立體電力線的畫法q=1e-9;X,Y ,Z=sphere(8);r0=0.18;r1=0.2;k=9e9;u0=k*q/r0;x=r1*X(:)'y=r1*Y(:)'z=r1*Z(:)'x=x;zeros(size(x);y=y;zeros(size(y);z=z;z

18、eros(size(z);plot3(x,y,z)hold on;%立體等勢線之畫法u=linspace(1,3,5)*u0;%畫 5 條等勢面時 , 各面的電勢 為 u=linspace(1,3,5)+u0, r=k*q./u;%各等勢面的半徑為 r=k6q. /u X,Y ,Z=sphere; Z(X<0&Y<0)=nan; surf(r(1)*X,r(1)*Y ,r(1)*Z);%第一到第五個球面 surf(r(2)*X,r(2)*Y ,r(2)*Z); surf(r(3)*X,r(3)*Y ,r(3)*Z); surf(r(4)*X,r(4)*Y ,r(4)*Z);

19、 surf(r(5)*X,r(5)*Y ,r(5)*Z); shading interp % 個顏色濃淡連續(xù)變化的命令 shading interp.xlabel('x','fontsize',16);ylabel('y','fontsize',16);zlabel('z','fontsize',16); title('正電荷電場線等勢面的三維 圖形','fo ntsize',20);clear all;clf;q1=1;q2=1;a=0.02;xx0=0.05;yy

20、0=0.04;k=9e9;x=linspace(-xx0,xx0,20);y=linspace(-yy0,yy0,50);X,Y=meshgrid(x);r1仁 sqrt(xx0/1.7-a)A2+(yy0/1.7)A2);r22=sqrt(xx0/1.7+a)A2+(yy0/1.7)Ar1 =sqrt(X-a).A2+Y.A2); r2=sqrt(X+a).A2+Y.A2); U=k*q1./r1+k*q2./r2; u0=k*q1/r11+k*q2/r22; u=linspace(u0,-u0,11); contour(X,Y,U,u,'k-');hold onEx,Ey=

21、gradient(-U);E=sqrt(Ex.A2+Ey.A2);Ex=Ex./E;Ey=Ey./E; dth1=20;th1=(dth1:dth1:180-dth1)*pi/180; r0=a/5;clear;x1=r0*cos(th1)+a; y1=r0*sin(th1); streamline(X,Y,Ex,Ey,x1,y1); streamline(-X,-Y,-Ex,-Ey,x1,-y1); q=abs(q1/q2);dth2=dth1/q; th2=(180-dth2:-dth2:dth2)*pi/180; x2=r0*cos(th2)-a;y2=r0*sin(th2); stre

22、amline(X,Y,Ex,Ey,x2,y2); streamline(X,-Y,Ex,-Ey,x2,-y2); grid on plot(a,0,'o','MarkerSize',12); plot(-a,0,'o','MarkerSize',12); xlabel('x','fontsize',16); ylabel('y','fontsize',16);title(' 一對點電荷的電場分布圖 '); clear;clear all;clf;q1=

23、1;q2=1;a=0.02;xx0=0.05;yy0=0.04;k=9e9;x=linspace(-xx0,xx0,20);y=linspace(-yy0,yy0,50);X,Y=meshgrid(x);r1仁 sqrt(xxO/1.7-aF2+(yyO/1.7F2);r22=sqrt(xxO/1.7+a)A2+(yyO/1.7)A2);r1=sqrt(X-a).A2+Y.A2); r2=sqrt(X+a).A2+Y.A2); U=k*q1./r1+k*q2./r2; uO=k*q1/r11+k*q2/r22; u=linspace(uO,-uO,11);contour(X,Y,U,u,'k-');