发布网友 发布时间:2022-05-15 19:07
共1个回答
热心网友 时间:2024-02-27 05:16
因为等高线全集中在一起了,改成contour(X,Y, V, linspace(min(V(:)), max(V(:)), 50));追问改了,还是不行啊。。。追答
更正:
等高线集中在一起的原因是靠近电荷的电场强度都接近无穷大了, 所以怎么分都不行. 所以*一个最大值.
clear all;
clc;
q=2e-6;k=9e9;x=0,y=0;
[X,Y]=meshgrid(-3:0.3:3);
r=sqrt((X-x).^2+(Y-y).^2);
V=k.*q.*(1./r);
V_max = 50*min(V(:));
V(V==Inf) = V_max;
V(V>V_max) = V_max;
[Ex,Ey]=gradient(-V);
E=sqrt(Ex.^2+Ey.^2);
Ex=Ex./E;
Ey=Ey./E;
hold on;
quiver(X,Y,Ex,Ey,0.4)
plot(x,y,'bo',x,y,'g+')
contour(X,Y,V,50);