电磁学
一、
1、点电荷的电场
研究真空中,两个带正电的点电荷,在电量相同和电量不同情况下的电场分布。 V=V1+V2=q1 4??0r1+q24??0r2,E=-▽V
2、程序实现
主程序文件名为point.m
clear all
ep0=8.85*le-12; %真空中的电容率
c0=1/(4*pi*ep0);
e=1.6e-10;
h=0.018;
x=-0.5:h:0.5;
y=-0.5:h:0.5;
str{1}=’两同号等量点电荷’;
str{2}=’两同号不等量点电荷’;
[X,Y]=meshgrid(x,y);
q=[e;1.9*e];
for i=1:2
V=c0*e./sqrt((X+0.2).^2+Y.^2)+c0.*q(i)./sqrt((X-0.2).^2+Y.^2); %求电势
[Ex,Ey]=gradient(-V,h); %求电场
figure(i)
counter(X(:,:,1),Y(:,:,1),V,… %等势面
[20,-20,19,-19,18,-18,17,-17,16,-16,15,-15,14,-14,13,-13,12,-12,11,-11,10,-10],’r’); Axis([-0.38,0.38,-0.28,0.28])
hold on
phi=0:pi/17:2*pi; %以下画电场线
sx1=0.2+0.01*cos(phi);
sy1=0.01*sin(phi);
streamline(X(:,:,1),Y(:,:,1),Ex,Ey,sx1,sy1);
hold on
sx2=-0.2+0.01*cos(phi);
sy2=0.01*sin(phi);
streamline(X(:,:,1),Y(:,:,1),Ex,Ey,sx2,sy2);
title(str(i))
text(-0.215,0,’+’,’fontsize’,20); %标示点电荷
text(0.185,0,’+’,’fontsize’,20);
end
电磁学
一、
1、点电荷的电场
研究真空中,两个带正电的点电荷,在电量相同和电量不同情况下的电场分布。 V=V1+V2=q1 4??0r1+q24??0r2,E=-▽V
2、程序实现
主程序文件名为point.m
clear all
ep0=8.85*le-12; %真空中的电容率
c0=1/(4*pi*ep0);
e=1.6e-10;
h=0.018;
x=-0.5:h:0.5;
y=-0.5:h:0.5;
str{1}=’两同号等量点电荷’;
str{2}=’两同号不等量点电荷’;
[X,Y]=meshgrid(x,y);
q=[e;1.9*e];
for i=1:2
V=c0*e./sqrt((X+0.2).^2+Y.^2)+c0.*q(i)./sqrt((X-0.2).^2+Y.^2); %求电势
[Ex,Ey]=gradient(-V,h); %求电场
figure(i)
counter(X(:,:,1),Y(:,:,1),V,… %等势面
[20,-20,19,-19,18,-18,17,-17,16,-16,15,-15,14,-14,13,-13,12,-12,11,-11,10,-10],’r’); Axis([-0.38,0.38,-0.28,0.28])
hold on
phi=0:pi/17:2*pi; %以下画电场线
sx1=0.2+0.01*cos(phi);
sy1=0.01*sin(phi);
streamline(X(:,:,1),Y(:,:,1),Ex,Ey,sx1,sy1);
hold on
sx2=-0.2+0.01*cos(phi);
sy2=0.01*sin(phi);
streamline(X(:,:,1),Y(:,:,1),Ex,Ey,sx2,sy2);
title(str(i))
text(-0.215,0,’+’,’fontsize’,20); %标示点电荷
text(0.185,0,’+’,’fontsize’,20);
end