利用matlab计算与模拟点电荷对的电场

2013-08-12薛海国

河南科技 2013年12期

薛海国

(四川理工学院理学院，四川自贡 643000)

1 matlab 和点电荷系的电场

点电荷是电磁学中最简单的一种理想化模型。单个点电荷在空间某点所激发的电场强度大小和电势U =用电场线和等势面直观形象地来表示这两个物理量。电场线从正电荷(无穷远)出发终止于负电荷(无穷远)，电场线的切线方向代表该点的场强方向，电场线的疏密来表示电场的强弱。电场线和等势面处处垂直。

然而在真空中某点所激发的电场通常是由点电荷组或者连续带电体所激发。对于点电荷组在空间某点所激发的电场强度和电势，我们可以根据电场强度矢量叠加和电势标量叠加进行数值求解。比较有代表性的就是要求电偶极子在延长线上和中垂线上的E 和U。简单的点电荷系的模型还有两个等量同号点电荷系，两个不等量同号点电荷系等。但是如果空间带电体的数量多于两个，并且在空间有一定规律分布，我们虽然能求出E 和U 的数值，但由于电场线和等势面的抽象性和复杂性，此时就很难想象出空间的分布。

matlab 有强大的数学运算能力和绘图功能，能直接利用软件提供的函数，编写较简单的程序即可解决相应的物理问题。已有文献［1］［2］［3］对电偶极子、电四极子，多个点电荷系统的电势分布和矢量场图进行计算仿真。

2 等量同号点电荷电场中的E 和U

如图所示，等量同号点电荷在场点P(x，y)产生的电势为

其中

E=－▽U，

其中，劈形算符为

在xy 平面上，场强只有两个分量

3 利用matlab 计算和模拟电场强度和电势

3.1 计算距离两个等量同号点电荷远处任一点的电位和场强

clear

＞＞syms q C0 r1 r2 x y U1 U2 a;

＞＞r1 =sqrt((x－a).^2 ＋y.^2);

＞＞r2 =sqrt((x ＋a).^2 ＋y.^2);

＞＞U1 =C0.* q./r1;U2 =C0.* q./r2;

＞＞U=U1 ＋U2;

＞＞Ex=－diff(U，x)，Ey=－diff(U，y)，AE=sqrt(Ex.^2＋Ey.^2)

运行结果

Ex =

1/2* C0* q/(x^2－2* x* a ＋a^2 ＋y^2)^(3/2)* (2* x－2* a)＋1/2* C0* q/(x^2 ＋2* x* a ＋a^2 ＋y^2)^(3/2)* (2*x ＋2* a)

Ey =

C0* q/(x^2－2* x* a ＋a^2 ＋y^2)^(3/2)* y ＋C0* q/(x^2 ＋2* x* a ＋a^2 ＋y^2)^(3/2)* y

AE =

((1/2* C0* q/(x^2－2* x* a ＋a^2 ＋y^2)^(3/2)* (2*x－2* a)＋1/2* C0* q/(x^2 ＋2* x* a ＋a^2 ＋y^2)^(3/2)*(2* x ＋2* a))^2 ＋(C0* q/(x^2－2* x* a ＋a^2 ＋y^2)^(3/2)* y ＋C0* q/(x^2 ＋2* x* a ＋a^2 ＋y^2)^(3/2)* y)^2)^(1/2)

在中垂线上时，x=0.键入:

x=0，

Ex=subs(Ex)，Ey=subs(Ey)，

Ex =0

Ey =2* C0* q/(a^2 ＋y^2)^(3/2)* y

在坐标原点时，x=0，y=0，键入:

x=0，y=0，

Ex=subs(Ex)，Ey=subs(Ey)

Ex =0

Ey =0

3.2 绘制距离两等量同号点电荷远处任一点的电位和场强分布图

clear