На основании закона Кулона электрическое поле, создаваемое одним дискретным зарядом q на расстоянии r, определяется как:
E = q / (4 * pi * e0 * r. ^ 2);
Если у вас есть несколько зарядов, вы можете использовать принцип суперпозиции и добавить вклад каждого заряда.
Осталось только создать сетку для вычисления электрического поля.Для этого вы можете использовать функцию Matlab meshgrid
.
Простой пример в Matlab:
k=1/4/pi/8.854e-12;
d=2;
q=[-1 1];
x=[-d/2 d/2];
y=[0 0];
dx=0.01;
X=(-d:dx:d);
Y=(-d:dx:d);
[X Y]=meshgrid(X,Y);
E=zeros(size(X));
for i=1:numel(q)
r=sqrt((X-x(i)).^2+(Y-y(i)).^2);
E=E+k*q(i)./r.^2;
end
E(isinf(E))=NaN;
figure;
contourf(X,Y,E);
axis image;
Надеюсь, он вам поможет.Вы можете прочитать документацию contourf
, чтобы настроить сюжет под свои нужды.