云南大学 软件学院 大学物理实验4

课程：  大学物理实验    学期：   2012-2013学年 第一学期    任课教师：             

班级：          学号：               序号：        姓名：           成绩：       

实验4  静电场

问题1：

(1):

(2):

(3):

clear;

k=1;q=1;l=1;

x=-6:0.1:6;

x1p=l:0.1:6;

x1n=-6:0.1:-l;

x11=-l:0.1:l;

E1p=(4*k*q*l.*x1p)./(x1p.^2-l^2).^2;%X大于L

E1n=(4*k*q*l.*x1n)./(x1n.^2-l^2).^2;%X小于L

E11=-k*q./(x11+l).^2-k*q./(x11-l).^2;%-LE2=4*k*q*l./x.^3;%电偶极子

figure;

plot(x1p,E1p);

hold on

plot(x11,E11);

hold on

plot(x1n,E1n);

hold on

plot(x,E2);

grid on

plot(-l, 0, 'ro',  -l, 0, 'r-', 'MarkerSize', 8); % 标出负电荷

hold on

plot(l, 0, 'ro',  l, 0, 'r+', 'MarkerSize', 8); % 标出正电荷

ylabel('电场强度：N/C');

xlabel('距离：m');

title('电场强度分布图');

text(1,-50,'点电荷-Ltext(1,50,'点电荷X>L曲线');

text(-2.5,-50,'点电荷Xaxis([-6 6 -100 100]);

(4):

问题2：

(1):

(2):

(3):

clear;

k=1;q=1;l=1;

y=-6:0.1:6;

y2p=0:0.1:6;

y2n=-6:0.1:0;

E1=-2*k*q./(y.^2+l^2)*l./sqrt(y.^2+l^2);%点电荷方程

E2n=4*k*q*l./y2n.^3;%电偶极子Y>0方程

E2p=-4*k*q*l./y2p.^3;%电偶极子Y<0方程

figure;

plot(y,E1);

hold on

plot(y2p,E2p);

hold on

plot(y2n,E2n);

grid on

ylabel('电场强度：N/C');

xlabel('距离：m');

title('Y轴电场强度分布图');

text(-1,-1,'点电荷曲线');

text(-2,-2,'电偶极子曲线');

text(2,-2,'电偶极子曲线');

axis([-6 6 -10 10]);

(4):