中南大学电磁场理论实验2.docVIP

题 目 姓名 学号 班级 任课老师 实验日期 电磁场理论实验二 -利用Matlab模拟点电荷电场的分布 实验目的： 熟悉单个点电荷及一对点电荷的电场分布情况； 学会使用Matlab进行数值计算，并绘出相应的圈形； 实验原理： 根据库伦定律：在真空中，两个静止点电荷之间的作用力与这两个电荷的电 量乘积成正比，与它们之间距离的平方成反比，作用力的方向在两个电荷的连线 上，两电荷同号为斥力，异号为吸力，它们之间的力F满足： TOC \o 1-5 \h \z F = k^^R （式 1） R- 由电场强度E的定义可知： E二學灸 （式2 对于点电荷，根据场论基础中的定义，有势场E的势函数为 - （式 3） R 而 E = -VU （式 4） 在Matlab中，由以上公式算出各点的电势U,电场强度E后，可以用Matlab 自带的库函数绘出相应电荷的电场分布情况。 三、实验内容： 根据库伦定律，利用Matlab强大的绘图功能画出单个点电荷的电场分布 情况，包括电力线和等势面。 根据库伦定律，利用Matlab强大的绘图功能画出一对点电荷的电场分布 r0=0.12;th=linspace(0,2*pi/13);[x,y]=pol2cart(th,rO);x=[x;0.1*x];y=[y;0.1*y];plot(x,y/b)grid onHold onplot(0/0/o/MarkerSize r0=0.12; th=linspace(0,2*pi/13); [x,y]=pol2cart(th,rO); x=[x;0.1*x]; y=[y;0.1*y]; plot(x,y/b) grid on Hold on plot(0/0/o/MarkerSize/12) xlabel(x?fontsize’,16) ylabel(lyl/fontsize,16) k=9e9; q=le-9; r0=0.1; uO=k*q/rO; u=linspace(l,3,7)*u0; x=linspace(-r0,r0,100); [X,Y]=meshgrid(x); r=sqrt(X.A2+Y.A2); U=k*q./r; contour(X,Y,U,u) 单个正点电荷的电场线等势线 0.15 0.1 0.05 -0.05 -0.15 ?0.2 -0.2 -0.15 -0.1 ?0.05 0.05 0.15 情况，包括电力线的分布和等势面。 实验代码1: %射线的半径 %电力线的角度 %将极坐标转化力直角叱标 %插入x的起始坐标 %插入y的起始坐标 %用蓝色画山所有电力线 %加网格 %保持图像 % 画电荷 %用16号字体标山X轴 %用16号字体标出Y轴 %设定k值 %设定电荷电量 %没定最大等势线的半径 %算出最小的电势 %求出各条等势线的电势大小 %将X坐标分成100等份 %在直角坐标中形成网格來标 %各个网格点到电荷点的距离 % 各点的电势 % rt出点电荷的电势面 titlef正电荷的电场线和等势线Yfontsjze,20 %显示标题 实验代码2: ? k=9e9 E=sqrt(Ex. E=sqrt(Ex.A2+Ey.A2); Ex=Ex./E; ql=le-9; % 设罝正电荷电量 q2=-le-9; % 设置负电荷电量 a=0.02; % 设置电荷到原点的距离 xx0=0.05; % 没置X轴的范围 yy0=0.04; % 设置Y轴的范围 x=linspace(-xx0,xx0,20); % 将X轴进行20等分 y=linspace(-yy0,yy0,50); % 将Y轴进行50等分 [X,Y]=meshgrid(x); % 形成网格坐标 rl=sqrt((X-a).A2+Y.A2); % 各点到正电荷的距离 r2=sqrt((X+a).A2+Y.八 2); % 各点到负电荷的距离 U=k*ql./rl+k*q2./r2; % 各点的电势 u0=500; % 设定S大电势的大小 u=linspace(uO,-uO/ll); % 计算各等势线的电势 contourfXXU^/k-); % 画出所有的等势线 Grid on % 形成网格 Hold on % 保持图形 plotfa^/o/MarkerSize.l?) % 画电荷 plotf-aAo/MarkerSize1,!^) °/c )画电荷 xlabelfx\fontsizei,16) % 用16号字体标出X轴 ylabel(Y/fontsize\16) % 用16号字体称出Y轴 [Ex,Ey]=gradient(-U); % 各点的场强的两个分量 %各点的合场强 %为使箭头等长，将场强归一化 Ey=Ey./E; quiver(X,Y,Ex,Ey); %标出各网点的电场强度方向 一对相异电荷的等势钱图