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

云南大学软件学院 实验报告 学号： 序号： 姓名： 成绩： 实验4 静电场 问题1： : : : 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;%-LXl E2=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,点电荷-LXL曲线); text(1,50,点电荷XL曲线); text(-2.5,-50,点电荷XL曲线); axis([-6 6 -100 100]); : 问题2： : : : 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;%电偶极子Y0方程 E2p=-4*k*q*l./y2p.^3;%电偶极子Y0方程 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]); : 2