等势线动图

clc;clear;%%等势线动图

for a=2:1:14;

q1=-1.602e-19;q2=1.602e-19;radius = 2;   %点电荷半径

x = 25; y = 25; %x和y轴显示边界值

k = 9e+9; %k = 1/4πε

delta = 1;%引入距离误差以消除点电荷处点位无限大的问题

phi = linspace(-pi,pi,20);  %电场线发射范围

[X,Y]=meshgrid(-x:0.5:x,-y:0.5:y);%获取左半部分矩阵

r1=sqrt((X+a).^2 + Y.^2 + delta);

r2=sqrt((X-a).^2 + Y.^2 + delta);

V1=q1*(k./r1) + q2*(k./r2);

subplot(1,2,1);meshc(V1);axis([0 100 0 100 -2e-9 2e-9]);

xlabel('x轴');ylabel('y轴');zlabel('z轴');

subplot(1,2,2);meshc(V1);axis([0 100 0 100 -2e-9 2e-9]);view(0,0);

xlabel('x轴');ylabel('y轴');

drawnow;%刷新图窗

end