題目 姓名 學號 班級任課老師實驗日期電磁場理論實驗二-利用Matlab模擬點電荷電場的分布?實驗目的：熟悉單個點電荷及一對點電荷的電場分布情況；學會使用Matlab進行數值計算，并繪出相應的圖形;二?實驗原理:根據庫倫定律：在真空中，兩個靜止點電荷之間的作用力與這兩個電荷的電量乘積成正比，與它們之間距離的平方成反比，作用力的方向在兩個電荷的連線上，兩電荷同號為斥力，異號為吸力，它們之間的力F滿足:（式1）（式2）（式3）（式4）F二（式1）（式2）（式3）（式4）R2由電場強度E的定義可知:E二怨R

R2對于點電荷，根據場論基礎中的定義，有勢場E的勢函數為u二kQR在Matlab中，由以上公式算出各點的電勢U,電場強度E后，可以用Matlab自帶的庫函數繪出相應電荷的電場分布情況。三、實驗內容：根據庫倫定律，利用Matlab強大的繪圖功能畫出單個點電荷的電場分布情況，包括電力線和等勢面。根據庫倫定律，利用Matlab強大的繪圖功能畫出一對點電荷的電場分布r0=0.12;th=linspace(0,2*pi,13);[x,y]=pol2cart(th,r0);x=[x;0.1*x];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,'b')gridonHoldonplot(0,0,'o','MarkerSize',12)xlabel('x','fontsize',16)ylabel('y','fontsize',16)k=9e9;q=1e-9;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;contour(X,YU,u)情況，包括電力線的分布和等勢面。實驗代碼1:%射線的半徑%電力線的角度%將極坐標轉化為直角坐標%插入x的起始坐標%插入y的起始坐標%用藍色畫出所有電力線%加網格%保持圖像%畫電荷%用16號字體標出X軸%用16號字體標出Y軸%設定k值%設定電荷電量%設定最大等勢線的半徑%算出最小的電勢%求出各條等勢線的電勢大小%將X坐標分成100等份%在直角坐標中形成網格坐標%各個網格點到電荷點的距離%各點的電勢%畫出點電荷的電勢面title（'正電荷的電場線和等勢線','fontsize',20） %顯示標題單個正點電荷的電場線等勢線0.15實驗代碼.2：0.10.05-0.1-0.150>1-0.05-0.15 -0.1 -0.05 0 0.05 0.1 0.15>>k=9e9q1=1e-9;q2=-1e-9;a=0.02;q1=1e-9;q2=-1e-9;a=0.02;xx0=0.05;yy0=0.04;x=linspace(-xx0,xx0,20);y=linspace(-yy0,yy0,50);[X,Y]=meshgrid(x);r1=sqrt((X-a).A2+Y.A2);r2=sqrt((X+a).A2+Y.A2);U=k*q1./r1+k*q2./r2;u0=500;u=linspace(u0,-u0,11);contour(X,YU,u,'k-');GridonHoldonplot(a,0,'o','MarkerSize',12)plot(-a,0,'o','MarkerSize',12)xlabel('x','fontsize',16)ylabel('y','fontsize',16)[Ex,Ey]=gradient(-U);E=sqrt(Ex.A2+Ey.A2);Ex=Ex./E;Ey=Ey./E;quiver(X,YEx,Ey);% 設置正電荷電量%設置負電荷電量% 設置電荷到原點的距離%設置X軸的范圍%設置Y軸的范圍%將X軸進行20等分%將Y軸進行50等分%形成網格坐標%各點到正電荷的距離%各點到負電荷的距離%各點的電勢% 設定最大電勢的大小% 計算各等勢線的電勢% 畫出所有的等勢線% 形成網格% 保持圖形%畫電荷%畫電荷%用16號字體標出X軸%用16號字體標出Y軸%各點的場強的兩個分量%各點的合場強%為使箭頭等長，將場強歸一化%標出各網點的電場強度方向title('一對相異電荷的等勢線圖和電場線圖',title('一對相異電荷的等勢線圖和電場線圖',—對相異電荷的等勢線圖和電場線圖0.050.040.030.020.010-0.01-0.02-0.03-0.04-0.05-0.05-0.04-0.03-0.02-0.01 0 0.01 0.02 0.03 0.04 0.05四、實驗總結通過這次的matlab仿真點電荷電場的實驗,我們對于點電荷周圍電