matlab 等位线,用MATLAB模拟双“点电荷”的电场和等位线

所谓双″点电荷″，就是空间中有一定间隔的两个单″点 电荷″。在电磁学中，经常要模拟点电荷的场和等位线，如果用人工绘图法来模拟，则模拟单个点电荷的场和等位线并不很难，但模拟双″点电荷″的场和等位线比较难。在本文中，笔者介绍了用 MATLAB 模拟双″点电荷″的场和等位线的方法，它包括数学建模和编写 MATLAB 程序代码。 1 数 学 建 模 设在(q1x，q1y)处有点电荷 q1，在(q2x，q2y)处有点电荷 q2。那么双″点电荷″所在平面上任何一点的电势和场强分别为： V(x,y)= q1 4πε0 3 1 r1 + q2 4πε0 6 1 r2 E=- !V 其中： r1= (x- q1x) 2 +(y- q1y) 2 " r2= (x- q2x) 2 +(y- q2y) 2 " 1 4πε0 =9×109 q1 和 q2 分别代表双″点电荷″中两个点电荷的电量。 (q1x，q1y)代表点电荷 q1 的坐标位置，(q2x，q2y)代表点电荷 q2 的坐标位置。 r1 为平面上任意一点到点电荷 q1 的距离，r2 为平面上任意一点到点电荷 q2 的距离。 2 编写 MATLAB 程序代码 %MATLAB程序 clear@ clfB %q1，q2 的坐标值和电量由键盘输入 q1x=input('q1x=') B q1y=input('q1y=')B q2x=input('q2x=') B q2y=input('q2y=') B q1=input('q1=') B q2=input('q2=') B %k 为比例系数 k=9e9B %建立平面坐标网点 x=- 6:0.6:6B y=xB [X,Y]=meshgrid(x,y)B %场强和等位线的计算代码 r1=sqrt((X- q1x).^2+(Y- q1y).^2) B %任意一点到点电荷 q1 的 距离 r2=sqrt((X- q2x).^2+(Y- q2y).^2) B %任意一点到点电荷 q2 的 距离 V=q13k3(1./r1)+q23k3(1./r2) B %计算电势 [Ex,Ey]=gradient(- V)B %计算场强 AE=sqrt(Ex.^2+Ey.^2)B Ex=Ex./AEB Ey=Ey./AEB cv=linspace(min(min(V)),max(max(V)),49)B %产生 49 个电位 contourf(X,Y,V,cv,'k- ') %画等位线axis('square') %图形标题 title('\fontname^黑体_\fontsize^16_场与等位线 ')；hold on quiver(X,Y,Ex,Ey,0.5) %场强箭头plot(q1x,q1y,'wo') %小圆圈代表点电荷 q1 2007 年 4 月 第 2 期 用MATLAB模拟双“点电荷”的电场和等位线 电 脑 学 习 莫 照* 胡奇光 摘 要 探讨了用 MATLAB 模拟双“点电荷”的场和等位线的方法, 详细介绍了其数学建模过程和 MATLAB 程序代码。关键词 MATLAB 场 等位线 中图分类号 TP391.7 文献标识码 A 文章编号:1002-2422( 2007) 02-0042-02 Simulating the Electric Field and Isoelectric Line of ″Double Point Charges″in MATLAB Mo Zhao Hu Qiguang Abstract The method of simulating the e