matlab结题报告(电偶极子的辐射场).doc

电偶极子的辐射场 背景与意义： 对于一个带电体来说，如果正负电荷呈电偶分布，正、负电荷的重心不重合，那么讨论这种带电体的电场时，可以把它模拟成两个相距很近的等量异号的点电荷+q 和?q，这样的带电系统称为电偶极子。实际生活中电偶极子的例子随处可见，例如，在研究电解质极化时，采用重心模型描述后电解质分子可等效为电偶极子；在电磁波的发射和吸收中电子做周期性运动形成振荡电偶极子；生物体所有的功能和活动都以生物电的形式涉及到电偶极子的电场等，当天线长度l远小于波长时，它的辐射就是电偶极辐射。因此，研究电偶极子在空间激发的电场问题具有重要意义。我们主要讨论宏观电荷系统在其线度远小于波长情形下的辐射问题。 基本内容介绍： 计算辐射场的一般公式 (1) （2） 其中 A(x 若电流J是一定频率的交变电流，有 Jx, 代入（3）式得 Ax 式中k=ω/c为波数。令 A 有 （6） 失势的展开 在失势公式（6）中，存在三个线度：电荷分布区域的线度l，它决定积分区x,的大小；波长λ=2π/k以及电荷到场点的距离r。我们研究分布于一个小区域的电流所产生的辐射。所谓小区域是指它的线度远小于波长λ 这种情况下，可以讲失势做展开得 Ax 电偶极辐射 我们研究展开式的第一项 （8） 先看电流密度体积分的意义。电流是有运动的带电粒子组成的。设单位体积内有ni个带电荷为qi,速度为vi 其中求和符号表示对各类带电粒子求和。上式也等于对单位体积内的所有带电粒子的qv求和。因此 式中求和符号表示对区域内所有带电粒子求和。但 式中P是电荷系统的电偶极矩。因此 ?l 如右图所示，当两个相距为?l的导体球组成，两个 导体之间由导线连接。当导线上有交变电流I时，两导体上的电荷±Q就交替变化，形成一个振荡电偶极子。这系统的电偶极矩为 当导线上有电流I时，Q的变化率为 因而体系的电偶极矩变化率为 （9） 由此可得，（8）式代表振荡电偶极矩产生的辐射 （10） 在计算电磁场时，需要对A作用算符?。我们只保留1/R 低次项，因而算符?不需作用到分母的R上，而仅需作用到因子eikr 由此得辐射场 (11) （12） 写成分量形式得 （13） (14) (15) 编程实现： 要实现电场的可视化操作，首先要得出电场线的方程 由电场个分量之间关系可得出 sin 2θ*cos 由式中K为积分常数，K取不同的值则得到不同的电力线。因此由(16)式可绘制出电偶极子的电力线族。在绘图时，需要将球坐标还原成直角坐标： 由于电场分布与?角无关，故电场分布关于z轴对称，因此可以只考虑某个过z轴的平面(如xoz平面)上电力线图，对于xoz平面，y=0，因此(4)式中球坐标 r=x2+z 且x、z的取值范围均为。 显然，(16)式可以写成 的形式，这其实是标量函数u(x,z)的等值线方程，因此电偶极子的电力线方程就是函数u(x,z)的等值线方程。MATLAB提供了一个专门的函数用于绘制标量函数u的等值线(或称等高线)图： [c, h] = contour (X, Z, U, V) 其中，X,Z,U为同维的矩阵，X，Z指定平面上点的x、z坐标，可由meshgrid命令取得，在本例中： x=–r:0.1:r; z=-r:0.1:r; [X,Z]=meshgrid(x, z); k是函数u(x,z)在坐标X,Z上的值，V 是向量，指定各条等高线的高度值（例如(16)式中的K值）。h 是返回的句柄值。 1. 以影片动画的方式仿真电偶极子辐射过程 要模拟电偶极子辐射场的动态过程，首先要绘制各个时刻的电力线图，即使用contour函数在?t取不同值的情况下绘制电力线方程式。绘制电力线图时应注意下面几个环节： 适当选取每个画面上电力线的根数，太多连成一片，太少没有真实感。有2个参数控制电力线的根数，一为K值，K每取一个值代表一条电力线（环形线，