简谐振动电偶极子辐射场研究.pptVIP

Page ? * 做简谐振动的电偶极子辐射电场 组员： 1.物理模型 1．1振动性偶极子的电场 1．2振动性偶极子电场的电场线 2.制作振动性偶极子电场线图像、动画并进行研究 2．1对方程的处理 2．2K值在振动性偶极子电场线方程中的含义 2. 3振动性偶极子电场线的演化进程研究 2.3.1单电场线演化进程研究 2.3.2全电场线演化进程研究 2．4对振动性偶极子电场的分析 1.物理模型： 1．1振动性偶极子的电场： 设振动性偶极子的电矩为 采用球坐标可得到在任意时刻t，空间任意处r的辐射电场： (1) 在krl的远区，库仑电场比感应电场弱得多，故远区的电场以感应电场为主导。而在 krl的近区与kr≈l的过渡区，库仑电场和感应电场不仅大小有差别，而且二者相位不尽相同，使此区域的电场呈现比较复杂的情况，这是需要进行认真分析的。 在近区 (kr1)，此振动性偶极子的库仑电场为： (2) 直接从式（2）出发进行分析很难看清楚，用图示的方法则可以非常清晰明了。而要作图，首先要导出辐射的电场线所应满足的方程。 1．2振动性偶极子电场的电场线方程 引入 代入(1)式即有： (3) (4) 将式（4）代入式（5），整理得 上式表明的全微分为零，即 =恒量。将C的表达式代入上式，并将 并入恒量，设为-K。有 在 （定值）的平面内： (5) (6) 2.制作振动性偶极子电场线图像、动画并进行研究： 2.1对方程的处理： 再作代换： 得 2.2K值在振动性偶极子电场线方程中的含义： K=0时作出的图为闭合圆环，这些圆不是电场线。 图中取相同K值的电场线用同一种线条表示。 由图中可见： 1)方程中 取值越大代表离K=0的圆越远的电场线 2)两个K=0的圆之间的电场线K值只取正或只取负，正负交替。 另外，K值较大的曲线未在图中表示出，它们处于电场的核心区域，满足1）的情况。 2．3振动性偶极子电场线的演化进程研究： 2.3.1单电场线的演化进程研究： K取定值（K取值如图），分别在n取不同值时对电场线方程作图。程序如下： syms x y; n=linspace(0.00,1.00,m);%%n在0~1中取m个值 for i=1:m; z=sqrt((1/(x^2+y^2)/2/pi/2/pi)+1)*(x^2/(x^2+y^2))*(cos(2*pi*n(i)-2*pi*sqrt(x^2+y^2)+atan(2*pi*sqrt(x^2+y^2))))+K; %%在K的位置输入具体数值 ezplot(z,[-1,1,-1,1]); %%取相应坐标作图 hold on; end title(振动性偶极子的电场线辐射的进程); k取不同值的图像.docx 图为运行结果，为了便于观察，前三图取m=5，后三图取m=10。K值互为相反数的电场线作出的演化进程图像是一致的，因而合为同一张图来表示。 K=±0.2 K=±0.5 K=±0.8 K=±1.0 K=±1.1 K=±1.3 动画程序： syms x y; m=moviein(300)%%共300帧 t=linspace(0.00,30.00,300);%%在0~30之间取300个t值 for i=1:300; z=sqrt((1/(x^2+y^2)/2/pi/2/pi)+1)*(x^2/(x^2+y^2))*(cos(2*pi*t(i)-2*pi*sqrt(x^2+y^2)+atan(2*pi*sqrt(x^2+y^2))))+K;%%按所想研究的电场线取K值 ezplot(z,[-3,3,-3,3]); %%作图 m(:,i)=getframe; end movie(m,3,20)%%每秒播放20帧，播放3遍 oujizifushedonghua.m 3.3.2全电场线演化进程研究: 取K=0,±0.2,±0.4,±0.6,±0.8,±1,±1.1,±1.2, ±1.4,±1.6分别对不同t值做图。程序如下： syms x y T; K=linspace(-1.6,1.6,17); t=T/64 n= (t/T ); for i=1:17; z=sqrt((1/(x^2+y^2)/2/pi/2/pi)+1)*(x^2/(x^2+y^2))*(cos(2*pi*n-2*pi*sqrt(x^2+y^2)+atan(2*pi*sqrt(x^2+y^2))))+K(i); ezplot(z,[-1,1,