電磁場(chǎng)與電磁波研究性學(xué)習(xí)

電磁場(chǎng)與電磁波研究性學(xué)習(xí)報(bào)告組長(zhǎng)：組員：北京交通大學(xué)電子信息工程學(xué)院2013年5月靜電場(chǎng)特性研究第二題：針對(duì)以下給定的電荷分布，用matlab仿真畫(huà)出對(duì)應(yīng)的電位和電場(chǎng)分布。并對(duì)結(jié)果進(jìn)行分析。（1）電荷為Q、相距d的電偶極子放置在真空中。（2）兩個(gè)接地的半無(wú)限大導(dǎo)體板分別放置在x軸和y軸上，形成900夾角，正電荷

放置在點(diǎn)（a，a）處。（3）一個(gè)兩維的電位分布近似用二次方表示如下：

為電荷分布。證明上述V函數(shù)滿足泊松方程。畫(huà)出電荷圖形和電位分布。（1）電荷為Q、相距d的電偶極子放置在真空中。【題目分析】由真空中靜電場(chǎng)點(diǎn)電荷公式：其中：（1）電荷為Q、相距d的電偶極子放置在真空中。【仿真程序及仿真結(jié)果】：q=1.6e-19;k=9e9;a=4;b=0;x=-12:0.6:12[X,Y]=meshgrid(x,x);r1=sqrt((X-a).^2+(Y-b).^2);r2=sqrt((X+a).^2+(Y+b).^2);V=q*k*(1./r1-1./r2);[Ex,Ey]=gradient(V);E=sqrt(Ex.^2+Ey.^2);Ex=Ex./E;Ey=Ey./E;cv=linspace(min(min(V)),max(max(V)),40);contour(X,Y,V,cv,'r','LineWidth',2);holdonquiver(X,Y,Ex,Ey,'b','LineWidth',0.01);plot(a,b,'w-')plot(-a,-b,'w+')xlabel('x');ylabel('y');title('電偶極子電位電場(chǎng)二維分布圖')holdoff（1）電荷為Q、相距d的電偶極子放置在真空中。（1）電荷為Q、相距d的電偶極子放置在真空中?！痉抡娉绦蚣胺抡娼Y(jié)果】：q=1.6e-19;k=9e9;x=-10:0.5:10;a=4;b=0;[X,Y]=meshgrid(x,x);r1=sqrt((X-a).^2+(Y-b).^2);r2=sqrt((X+a).^2+(Y+b).^2);V=q*k*(1./r1-1./r2);meshc(X,Y,V);title('電偶極子電勢(shì)三維分布圖');靜電場(chǎng)特性研究（2）兩個(gè)接地的半無(wú)限大導(dǎo)體板分別放置在x軸和y軸上，形成900夾角，正電荷

放置在點(diǎn)（a，a）處?！绢}目分析】利用真空中靜電場(chǎng)鏡像法可以將四個(gè)點(diǎn)電荷等效為空間四個(gè)點(diǎn)電荷在第一象限產(chǎn)生的電勢(shì)疊加和電場(chǎng)矢量和。如下圖所示：

Q2(-a,a)Q1(a,a)

Q3(-a,-a)Q4(a,-a)因此，由于電荷對(duì)稱(chēng)，電場(chǎng)應(yīng)垂直與x=0和y=0軸。其中：（2）兩個(gè)接地的半無(wú)限大導(dǎo)體板分別放置在x軸和y軸上，形成900夾角，正電荷

放置在點(diǎn)（a，a）處?！痉抡娉绦蚣胺抡娼Y(jié)果】：a=4;b=4;x=-12:1.2:12;[X,Y]=meshgrid(x,x);r1=sqrt((X-a).^2+(Y-b).^2);r2=sqrt((X+a).^2+(Y-b).^2);r3=sqrt((X+a).^2+(Y+b).^2);r4=sqrt((X-a).^2+(Y+b).^2);V=(1./r1-1./r2+1./r3-1./r4);[Gx,Gy]=gradient(V);G=sqrt(Gx.^2+Gy.^2);Gx=Gx./G;Gy=Gy./G;cv=linspace(min(min(V)),max(max(V)),40);contour(X,Y,V,cv,'r','LineWidth',1.45);holdonquiver(X,Y,Gx,Gy,'b','LineWidth',0.01)plot(a,b,'w-')plot(-a,-b,'w+')xlabel('x');ylabel('y');title('電場(chǎng)電勢(shì)二維分布圖')holdoff（2）兩個(gè)接地的半無(wú)限大導(dǎo)體板分別放置在x軸和y軸上，形成900夾角，正電荷

放置在點(diǎn)（a，a）處。（2）兩個(gè)接地的半無(wú)限大導(dǎo)體板分別放置在x軸和y軸上，形成900夾角，正電荷

放置在點(diǎn)（a，a）處?！痉抡娉绦蚣胺抡娼Y(jié)果】：x=-10:0.5:10;a=4;b=4;[X,Y]=meshgrid(x,x);r1=sqrt((X-a).^2+(Y-b).^2);r2=sqrt((X+a).^2+(Y-b).^2);r3=sqrt((X+a).^2+(Y+b).^2);r4=sqrt((X-a).^2+(Y+b).^2);E=(1./r1-1./r2+1./r3-1./r4);meshc(X,Y,E);title('電勢(shì)三維分布圖');（2）兩個(gè)接地的半無(wú)限大導(dǎo)體板分別放置在x軸和y軸上，形成900夾角，正電荷

放置在點(diǎn)（a，a）處。（3）一個(gè)兩維的電位分布近似用二次方表示如下：

證明上述V函數(shù)滿足泊松方程。畫(huà)出電荷圖形和電位分布?！绢}目分析】證明：∴滿足泊松方程（3）一個(gè)兩維的電位分布近似用二次方表示如下：

證明上述V函數(shù)滿足泊松方程。畫(huà)出電荷圖形和電位分布?！痉抡娉绦蚣胺抡娼Y(jié)果】：pv=1.6e-16;t=8.85e-12;x=-5:0.5:5;[X,Y]=meshgrid(x,x);V=-pv.*(X.^2+Y.^2)./(4*t);figure(1)mesh(X,Y,V);title('電位分布圖');figure(2)surf(X,Y,V);view(0,0);xlabel('y(x)')text(-2,-5,'V');title('截面圖');（3）一個(gè)兩維的電位分布近似用二次方表示如下：

證明上述V函數(shù)滿足泊松方程。畫(huà)出電荷圖形和電位分布。（3）一個(gè)兩維的電位分布近似用二次方表示如下：

證明上述V函數(shù)滿足泊松方程。畫(huà)出電荷圖形和電位分布。靜電場(chǎng)特性研究心得體會(huì)：我們小組選擇這個(gè)題目主要是想學(xué)習(xí)一下通過(guò)Matlab來(lái)制作電場(chǎng)線和電位線。大家都知道，電偶極子之間的電場(chǎng)、電位線美妙而且形象，能夠很好地描述點(diǎn)電荷之間的電場(chǎng)和電位分布。但是它也十分復(fù)雜，通過(guò)手繪則會(huì)帶來(lái)很大的工作量，不太現(xiàn)實(shí)。而Matlab則可以很好地解決這個(gè)問(wèn)題，用它繪出的二維電場(chǎng)分布和三維電勢(shì)分布圖形精確、形象，能很好地反映點(diǎn)偶極子之間的電場(chǎng)關(guān)系。靜電場(chǎng)特性研究這次的研究性學(xué)習(xí)，在理論上并沒(méi)有多復(fù)雜，因?yàn)轭}目所涉及內(nèi)容在教材上都能找到，我們的時(shí)間主要花費(fèi)在查詢(xún)函數(shù)的使用和編寫(xiě)程序上。由于我們接觸的MATLAB還是很淺的層次，所以需要大量的查找資料和函數(shù)設(shè)計(jì)，雖然費(fèi)時(shí)費(fèi)力，但我們感覺(jué)受益匪淺。通過(guò)這次研討，我們對(duì)點(diǎn)偶極子電場(chǎng)和電勢(shì)的分布有了更加清晰、深刻的認(rèn)識(shí)。我們學(xué)會(huì)了微分計(jì)算、梯度計(jì)算以及多種作圖函數(shù)的使用，也為以后的專(zhuān)業(yè)學(xué)習(xí)和設(shè)計(jì)豐富了知識(shí)，開(kāi)拓了思維。相信付出的汗水絕對(duì)不會(huì)白費(fèi)。靜電場(chǎng)特性研究與此同時(shí)，在研討查找資料的過(guò)程中，我們發(fā)現(xiàn)Malab遠(yuǎn)比想象中的強(qiáng)大，它還能做出很多復(fù)雜形狀，在我們今后的學(xué)習(xí)中如果能用好Matlab這個(gè)工具，我們的學(xué)習(xí)和對(duì)知識(shí)的理解將變得更加輕松。這需要我們繼續(xù)努力學(xué)好Matlab。靜電場(chǎng)特性研究遇到的問(wèn)題及建議：研討過(guò)程中我們也遇到了一些問(wèn)題，比如：Matlab的相關(guān)函數(shù)很復(fù)雜，種類(lèi)也很多，不同的函數(shù)會(huì)有不同的效果，對(duì)這些函數(shù)的不熟悉讓我們?cè)谔暨x和使用函數(shù)的時(shí)候進(jìn)度相當(dāng)慢，所以希望老師可以提供給我們相關(guān)的函數(shù)，這樣我們就可以不用花費(fèi)大量的時(shí)間去找合適的函數(shù)了。還有就是在描繪二維圖時(shí)，我們想將電荷的極性標(biāo)在圖上，這樣就更加一目了然，然而加上去的效果并不是太好，如下圖，直到最后還是沒(méi)能實(shí)現(xiàn)，希望之后老師能教我們解決方案。靜電場(chǎng)特性研究箭頭所示，只有一個(gè)電荷極性畫(huà)出來(lái)了，但顯示效果也并不