资源描述：
《第十组-MATLAB仿真镜像法.docx》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、电磁场与电磁波大作业姓名：马杰学号学号：刘万康秦旺潘尚伟张一鸣黄璞一，实验要求：运用镜像电荷方法计算和模拟无穷大直角导体内部点电荷电位，电场计算及分布图。二，什么是镜像法镜像法是解静电边值问题的一种特殊方法。它主要用来求分布在导体附近的电荷（点电荷、线电荷）产生的场。如在实际工程中，要遇到水平架设的双线传输线的电位、电场计算问题。当传输线离地面距离较小的时候，要考虑地面的影响，地面可以看作一个无穷大的导体平面。由于传输线上所带的电荷靠近导体平面，导体表面会出现感应电荷。此时地面上方的电场由原电荷和感应电荷共同产生。镜像法是应用唯一性
2、定理的典型范例。在镜像法应用中应注意以下几点：1)镜像电荷位于待求场域边界之外。2)将有边界的不均匀空间处理为无限大均匀空间，该均匀空间中媒质特性与待求场域中一致。3)实际电荷(或电流)和镜像电荷(或电流)共同作用保持原边界上的边界条件不变。三，点电荷对相互正交的两个无限大接地导体平面的镜像设在自由空间有一点电荷q位于无限大接地导体平面上方，且与x导体平面距离为a,与z导体平面距离为d.点电荷无限大导体平面左上半空间的电位分布和电场强度计算可用镜像法解决。待求场域为z>0，x>0空间，边界为x=0,z=0的无限大导体平面，边界条件为
3、在边界上电位为零，即∅x,y,z=0(2.1)设想将无限大平面导体撤去，整个空间为自由空间。在原边界之外(-a,0,d)放置一镜像电荷q2，q2=-q，在(-a,0,-d)放置一镜像电荷q3，当q3=+q，在(a,0,-d)放置一镜像电荷q4，当q4=-q，如图2所示。点电荷q1和镜像电荷q2,q3.q4在边界上产生的电位满足式(2.1)所示的边界条件。图2镜像法图示根据镜像法原理，在x>0,z>0空间的电位为点电荷q1和镜像电荷q2,q3,q4所产生的电位叠加，即ϕ=14πε0{qr1-qr2+qr3-qr4}ϕ=q4πε0{1(
4、x-a2+y2+(z-d)2)12-1(x+a2+y2+(z-d)2)12+1(x+a2+y2+(z+d)2)12-1(x-a2+y2+(z+d)2)12}(2.2)上半空间任一点的电场强度为E=-∇∅电场强度E的三个分量分别为EX=q4πε0{x-a(x-a2+y2+(z-d)2)32-x+a(x+a2+y2+(z-d)2)32+x+a(x+a2+y2+(z+d)2)32-x-a(x-a2+y2+(z+d)2)32}(2.3.1)Ey=q4πε0{y(x-a2+y2+(z-d)2)32-y(x+a2+y2+(z-d)2)32+y(
5、x+a2+y2+(z+d)2)32-y(x-a2+y2+(z+d)2)32}(2.3.2)Ez=q4πε0{z-d(x-a2+y2+(z-d)2)32-z-d(x+a2+y2+(z-d)2)32+z+d(x+a2+y2+(z+d)2)32-z+d(x-a2+y2+(z+d)2)32}(2.3.3)可见，在导体表面x=0,处，Ey=Ez=0，只有Ex存在,在导体表面z=0处Ex=Ey=0，即导体表面上法向电场存在。导体表面感应电荷分布可由边界条件决定，即ρx=0=aq2π-1a2+y2+z-d232+1a2+y2+z+d232(2.4
6、.1)ρz=0=dq2π-1(x-a)2+y2+d232+1(x+a)2+y2+d232(2.4.2)由式(2.4)可以看出，导体表面上感应电荷分布是不均匀的，导体表面上感应电荷为qin=(ρx=0+ρz=0)dS=aq2π(-∞+∞-∞+∞(-1a2+y2+z-d232+1a2+y2+z+d232)+dq2π-∞+∞-∞+∞(-1x-a2+y2+d232+1x+a2+y2+d232)dxdy=-q四，MATLAB实现电场线和电位（镜像电荷理想等效）clearq1=1;q2=-1;q3=1;q4=-1;x1=2;x2=-2;y1=2
7、;y2=-2;xm=5;ym=5;x=linspace(-xm,xm);y=linspace(-ym,ym);[X,Y]=meshgrid(x,y);R1=sqrt((X-x1).^2+(Y-y1).^2)R2=sqrt((X-x2).^2+(Y-y1).^2);R3=sqrt((X-x2).^2+(Y-y2).^2);R4=sqrt((X-x1).^2+(Y-y2).^2);U=q1./R1+q2./R2+q3./R3+q4./R4;u=-4:0.5:4;figurecontour(X,Y,U,u)gridonlegend(num
8、2str(u'))holdonplot([-xm;xm],[0;0])plot([0;0],[-ym;ym])plot(x1,y1,'o','MarkerSize',12)plot(x2,y1,'o','MarkerSize',12)
