北理matlab实验报告,北理工电磁场实验报告.docx

北理工电磁场实验报告

、实验四实验报告　　实验名称：一条指令的执行过程　　学号姓名班级：实验时间：年月日　　实验报告表4-1一条指令执行过程记录表　　1　　2　　CENTRALSOUTHUNIVERSITY　　题目利用Matlab模拟点电荷电场的分布　　姓名刘畅学号班级电气试验1101班　　任课老师李志勇实验日期　　XX年11月10日　　一、实验目的：　　1、熟悉单个点电荷及一对点电荷的电场分布情况；2、学会使用Matlab进行数值计算，并绘出相应的图形；　　二、实验原理：　　根据库伦定律：在真空中，两个静止点电荷之间的作用力与这两个电荷的电量乘积成正比，与它们之间距离的平方成反比，作用力的方向在两个电荷的连线上，两电荷同号为斥力，异号为吸力，它们之间的力F满足：　　F?k　　Q1Q2?　　R2　　R(式1)　　由电场强度E的定义可知：　　E?　　kQ?　　R　　R2(式2)　　对于点电荷，根据场论基础中的定义，有势场E的势函数为　　U?　　kQ　　R(式3)　　而E???U(式4)在Matlab中，由以上公式算出各点的电势U，电场强度E后，可以用Matlab自带的库函数绘出相应电荷的电场分布情况。　　三、实验内容：　　1．根据库伦定律，利用Matlab强大的绘图功能画出单个点电荷的电场分布情况，包括电力线和等势面。　　实验代码：　　r0=;%取射线半径　　th=linspace(0,2*pi,13);%电力线的角度　　[x,y]=pol2cart(th,r0);%将极坐标转化为直角坐标x=[x;*x];%插入X的起始坐标y=[y;*y];%插入Y的起始坐标plot(x,y,'b')%用蓝色画出所有电力线gridon%加网格Holdon%保持图像　　plot(0,0,'o','MarkerSize',12)%画电荷　　xlabel('x','fontsize',16)%用16号字体标出X轴ylabel('y','fontsize',16)%用16号字体标出Y轴k=9e9;%设定K值　　q=1e-9;%设定电荷电量　　r1=;%设定最大等势线的半径u0=k*q/r1;%算出最小的电势　　u=linspace(1,3,7)*u0;%求出各条等势线的电势大小x=linspace(-r0,r0,100);%将X坐标分成100等份[X,Y]=meshgrid(x);%在直角坐标中形成网格坐标r=sqrt(X.^2+Y.^2);%各个网格点到电荷点的距离U=k*q./r;%各点的电势　　contour(X,Y,U,u);画出点电荷的点失眠　　title('单个正点电荷的电场线与等电势','fontsize',20);%显示标题　　截图：　　2．根据库伦定律,利用Matlab强大的绘图功能画出一对点电荷的电场分布情况,包括电力线的分布和等势面。　　实验代码：　　k=9e9;%设定K值　　q1=1e-9;%设定正电荷的电量q2=-1e-9;%设定负电荷的电量a=;%设置电荷到原点的距离xx0=;%设定X的范围yy0=;%设定Y的范围　　x=linspace(-xx0,xx0,20);%将X轴进行20等分y=linspace(-yy0,yy0,50);%将Y轴进行50等分[X,Y]=meshgrid(x);%形成网格坐标　　r1=sqrt((X-a).^2+Y.^2);%各点到正电荷的距离r2=sqrt((X+a).^2+Y.^2);%各点到负电荷的距离U=k*q1./r1+k*q2./r2;%各点的电势u0=500;%设定最大电势的大小　　u=linspace(u0,-u0,11);%计算各等势线的电势contour(X,Y,U,u,'k-');%画出所有等势线Gridon%画网格Holdon%保持图像　　plot(a,0,'o','MarkerSize',12)%画电荷plot(-a,0,'o','MarkerSize',12)%画电荷　　xlabel('x','fontsize',16)%用16号字体标出X轴ylabel('y','fontsize',16)%用16号字体标出Y轴[Ex,Ey]=gradient(-U);%各点的场强的两个分量E=sqrt(Ex.^2+Ey.^2);%各点的合场强　　Ex=Ex./E;%为使箭头等长，将场强归一化Ey=Ey./E;%%为使箭头等长，将场强归一化quiver(X,Y,Ex,Ey);%标出电场强度方向　　title('一对相异电荷的等势线图和电场线图','fontsize',20);%显示标题　　截图：　　3．实验内容1中，可以在正电荷和负电荷中任选一组画出其电场分布，实验内容