可视化计算离线作业(3)

3.（10分）在空间有一个带正电的点电荷，请画出此点电荷的空间电位分布与电场强度的空间分布图。

hh=figure('numbertitle','off','name','点电荷的空间电位分布与电场强度的空间分布图<2015年秋 >');

r0=0.12;

th=linspace(0,2*pi,13);

[x,y]=pol2cart(th,r0)

x=[x;0.1*x]

y=[y;0.1*y]

plot(x,y)

grid on

hold on

plot(0,0,'o','MarkerSize',12)

k=9e9;

q=1e-9;

r0=0.12;

u0=k*q/r0;

u=linspace(1,3,7)*u0

x=linspace(-r0,r0,100)

[X,Y]=meshgrid(x)

r=sqrt(X.^2+Y.^2);

U=k*q./r;

contour(X,Y,U,u)

xlabel('x','fontsize',16) 0 100 200 300 400 500 600

-200 -100

100

200

300

400 500

??