简谐振动电偶极子辐射场分析.docxVIP

2Pk3E - 2Pk3 E -4矶 P.k3 S签昭［爲—存°曲如+击叽心「+彳） (2-1) 研究简谐振动的电偶极子电场 【摘要】本文首先对振动性偶极子电场的物理模型进行简要的分析并推导出其电场线方程，然 后利用数学软件Matlab对隐函数直接作图的功能作出其电场线的演化进程图像，并用Matlab动画模 拟其电场线辐射过程，最后结合图像和动画对了振动性偶极子电场进行具体的分析，得出结论。特别 是，文中清楚地模拟了部分不闭合电场线“分裂”出闭合电场线的过程，这在一般论文和教材中较为 少见。 【关键字】振动性偶极子（振荡电偶极子 偶极振子）；Matlab;作图；动画；感应电场；库仑电场 引言 振动性偶极子是电磁波辐射理论的基础，对其电场辐射情况的研究具有重要的意义。但由于振动性 偶极子电场的概念抽象，理论计算过程又十分复杂，推导和掌握需要较深的数学基础，而图形绘制也要 考虑诸多因素，极其繁琐，致使这方面的研究较为困难。使用Matlab则可以轻松地应对这些问题，它 能够针对振动性偶极子电场的各个参量变化时的特点快速地绘制出其电场线图像。在图形的帮助下，就 很容易对其电场进行简明而清楚的分析。 物理模型 2?1振动性偶极子的电场 设振动性偶极子的电矩为 P = exPo cos M 采用球坐标可得到在任意时刻/,空间任意处r的辐射电场叫 cos（血一 Zr） + — cos（血一比厂 +兰） （kr）2 2 上式中k=~. c 在kr?l的远区，库仑电场比感应电场弱得多，故远区的电场以感应电场为主导。而在kr?l的近 区与kr~l的过渡区，库仑电场和感应电场不仅大小有差别，而且二者相位不尽相同，使此区域的电场 呈现比较复杂的情况，这是需要进行认真分析的，也是本文的重点。 由（1）式可直接推导出振动性偶极子在在kr?l的近区的电场为⑷： Er = -^^COS0COSM 4亦0厂Ee = — sin 0cos cot 4 Er = -^^COS0COSM 4亦0厂 Ee = — sin 0cos cot 4亦0厂 E(p = 0 (2-2) 直接从式（2）出发进行分析很难看清楚，用图示的方法则可以非常清晰明了。而要作图, 导出辐射的电场线所应满足的方程。 首先耍 2. 2振动性偶极子电场的电场线方程 由式（1）导出振动性偶极子电场的电场线方程⑵，根据此方程即画出电场线图。引入 C =初 PN [— 4-l]l/2 sin 0cos[cot -kr-\- arctan伙r)]121 4殆)厂(kr) (2-3) 由式（1） 可以验证E = VxCr1,即有: (2-4)Er = ——- (CsinO) r sin Odd (2-4) = (rC)r dr E(p = O 在(p-(pe（定值）的平面内，E线满足的微分方程:dr Errd 0 Eg 在(p-(pe （定值）的平面内，E线满足的微分方程: dr Er rd 0 Eg (2-5) 将式（4） 代入式(5),整理得 一丄?(厂 C)d 厂=亠 2(Csin d)dO r dr sin 0 d0 —(Crsin + ——(Crsin ff)d = 0 ⑵ dr (2-6) 上式表明Csin的全微分为零，即Csin二恒量。将C的表达式代入上式中， 并将 空并入 4 兀￡（） 恒量,设为有 I 1 忑)2 +1F sin2 cos[d-Q + arctan伙厂)]+ K = 0121 1 补? 2 t (2-7)(上)2 +1 卜 sirT 0cos[27T--kr + arctan(Z:r)| + K = 0 (2-7) 这就是振动性偶极子辐射的电场线所应满足的方程，当K取不同的值时得到不同的辐射电场线。 用MATLAB制作振动性偶极子电场线的图像、动画并进行研究 3. 1对方程的处理 首先将（7）式写成直角坐标式X 好写成 —子[（訐+（亍” 而 sin ^ = — = / 厂 a/x2 + z2 将以上两式代入式（7）得到4,4所满足的方程卩 A A +11(2兀)2吟2+(自2]X再作代换：y= x ,A.ZI1 1[————+1]2-（2龙凡午+于） +11 (2兀)2吟2+(自2] X 再作代换：y= x , A. Z I 1 1 [————+1]2- （2龙凡午+于） 现在就可以根据上式编写Matlab程序进行作图以进行研究了。 ?= y, — = n，得 T x~ x2 +z2 cos] 27vn 一 2^(x2 + y2 y + arctan [2^(x2 +)，)打]+ K = 0 ⑶(2-8) 3.2 K值在振动性偶极子电场线方程中的含义 振动性偶极子电场的电场线图象3-1 n=（）时不同 振动性偶极子电场的电场线图象 3-1 n=（）时