matlab偶极矩电场强度分布图_MATLAB点电荷电偶极子程序

作业九：

点电荷

%点电荷电场的可视化

%by 小辣椒

%2015.4.15

%--参数设定及自变量取样--

q=1.6e-19;%点电荷电量

C0=1/(4*pi*8.85e-12);%合并系数，

%C0=1/(4*π*ε0)

%ε0：真空介电常数

x=-5:0.3:5;%自变量

y=x;%自变量

%——求解电势——

[X,Y]=meshgrid(x,y);%用于三维曲面的分格线坐标，自变量坐标矩阵

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

U=q.*C0.*(1./r);

%——绘制电势分布图——

subplot(2,2,1);

mesh(X,Y,U);%三维网线图

xlabel(x);

ylabel(y);

zlabel(U);

title(点电荷电势分布图);

%——对数坐标绘图——

subplot(2,2,3);

mesh(X,Y,log10(U));%三维网线图

xlabel(x);

ylabel(y);

zlabel(log10(U));

title(点电荷电势分布图);

%——求解电场强度——

[Ex,Ey]=gradient(-U,0.3,0.3);%gradient：梯度

AE=sqrt(Ex.^2+Ey.^2);

Ex=Ex./AE;

Ey=Ey./AE;

%——绘制场强分布图——<