电磁场实验指导书及实验报告

CENTRALSOUTHUNIVERSITY题目利用MATLAB模拟点电荷电场的分布姓名XXXX学号XXXXXXXXXX班级电气XXXX班任课老师XXXX实验日期201010电磁场理论实验一利用MATLAB模拟点电荷电场的分布1实验目的1熟悉单个点电荷及一对点电荷的电场分布情况；2学会使用MATLAB进行数值计算，并绘出相应的图形；2实验原理根据库伦定律在真空中，两个静止点电荷之间的作用力与这两个电荷的电量乘积成正比，与它们之间距离的平方成反比，作用力的方向在两个电荷的连线上，两电荷同号为斥力，异号为吸力，它们之间的力F满足RQKF21式1由电场强度E的定义可知RKE2式2对于点电荷，根据场论基础中的定义，有势场E的势函数为式3KQU而式4在MATLAB中，由以上公式算出各点的电势U，电场强度E后，可以用MATLAB自带的库函数绘出相应电荷的电场分布情况。3实验内容1单个点电荷点电荷的平面电力线和等势线真空中点电荷的场强大小是EKQ/R2,其中K为静电力恒量,Q为电量,R为点电荷到场点PX,Y的距离。电场呈球对称分布,取电量Q0,电力线是以电荷为起点的射线簇。以无穷远处为零势点,点电荷的电势为UKQ/R,当U取常数时,此式就是等势面方程等势面是以电荷为中心以R为半径的球面。平面电力线的画法在平面上,电力线是等角分布的射线簇,用MATLAB画射线簇很简单。取射线的半径为都取国际制单位R0012,不同的角度用向量表示单位为弧度THLINSPACE0,2PI,13。射线簇的终点的直角坐标为X,YPOL2CARTTH,R0。插入X的起始坐标XX01X同样插入Y的起始坐标,YY01Y,X和Y都是二维数组,每一列是一条射线的起始和终止坐标。用二维画线命令PLOTX,Y就画出所有电力线。平面等势线的画法在过电荷的截面上,等势线就是以电荷为中心的圆簇,用MATLAB画等势线更加简单。静电力常量为K9E9,电量可取为Q1E9最大的等势线的半径应该比射线的半径小一点R001。其电势为U0K8Q/R0。如果从外到里取7条等势线,最里面的等势线的电势是最外面的3倍,那么各条线的电势用向量表示为ULINSPACE1,3,7U0。从R0到R0取偶数个点,例如100个点,使最中心点的坐标绕过0,各点的坐标可用向量表示XLINSPACER0,R0,100,在直角坐标系中可形成网格坐标X,YMESHGRIDX。各点到原点的距离为RSQRTX2Y2,在乘方时,乘方号前面要加点,表示对变量中的元素进行乘方计算。各点的电势为UK8Q/R,在进行除法运算时,除号前面也要加点,同样表示对变量中的元素进行除法运算。用等高线命令即可画出等势线CONTOURX,Y,U,U,在画等势线后一般会把电力线擦除,在画等势线之前插入如下命令HOLDON就行了。平面电力线和等势线如图1,其中插入了标题等等。越靠近点电荷的中心,电势越高,电场强度越大,电力线和等势线也越密。02015010050005010150201501005000501015XY一一一一一一一一一一一一一图1源程序点电荷的平面电力线和等势线平面电力线的画法Q1E9R0012THLINSPACE0,2PI,13X,YPOL2CARTTH,R0XX01XYY01YPLOTX,YGRIDONHOLDONPLOT0,0,O,MARKERSIZE,12XLABELX,FONTSIZE,16YLABELY,FONTSIZE,16TITLE单个点电荷的电场线与等势线,FONTSIZE,20平面等势线的画法K9E9R001U0KQ/R0ULINSPACE1,3,7U0XLINSPACER0,R0,100X,YMESHGRIDXRSQRTX2Y2UKQ/RHOLDONCONTOURX,Y,U,UCLEAR点电荷的立体电力线和等势面立体电力线的画法先形成三维单位球面坐标,绕Z轴一周有8条电力线X,Y,ZSPHERE8,每维都是99的网格矩阵,将X化为行向量,就形成各条电力线的终点X坐标XR0X,其他两个坐标也可同样形成终点坐标YR0Y,ZR0Z。对X坐标插入原点XXZEROSSIZEX,其他两个坐标如下形成YYZEROSSIZEY,ZZZEROSSIZEZ,用三维画线命令PLOT3X,Y,Z,就画出所有电力线。立体等势面的画法画5条等势面时,各面的电势为ULINSPACE1,3,5U0,各等势面的半径为RK6Q/U,其中第一个球面的半径为RRR1。三维单位球面的坐标可由X,Y,ZSPHERE命令形成,每维都是2121的网格矩阵,由于外球会包围内球,因此把球面的四分之一设为非数,表示割去该部分ZX0X,Y,ZSPHERE8R0018R102K9E9U0KQ/R0XR1XYR1YZR1ZXXZEROSSIZEXYYZEROSSIZEYZZZEROSSIZEZPLOT3X,Y,ZHOLDON立体等势线之画法ULINSPACE1,3,5U0画5条等势面时,各面的电势为ULINSPACE1,3,5U0,RKQ/U各等势面的半径为RK6Q/UX,Y,ZSPHEREZX0SURFR1X,R1Y,R1Z第一到第五个球面SURFR2X,R2Y,R2ZSURFR3X,R3Y,R3ZSURFR4X,R4Y,R4ZSURFR5X,R5Y,R5ZSHADINGINTERP个颜色浓淡连续变化的命令SHADINGINTERP。XLABELX,FONTSIZE,16YLABELY,FONTSIZE,16ZLABELZ,FONTSIZE,16TITLE正电荷电场线等势面的三维图形,FONTSIZE,20CLEAR2一对点电荷平面等势线的画法仍然用MATLAB的等高线命令画等势线。对于正负两个点电荷,电量不妨分别取Q12E9,Q21E9,正电荷在X轴正方,负电荷在X轴负方,它们到原点的距离定为A002假设平面范围为XX0005,YY0004,两个坐标向量分别XLINSPACEXX0,XX0,20和YLINSPACEYY0,YY0,50。设置平面网格坐标为X,YMESHGRIDX,各点到两电荷的距离分别为R1SQRTXA2Y2和R2SQRTXA2Y2。各点的电势为UK6Q1/R1K6Q2/R2,取最高电势为U050,最低电势取其负值。在两者之间取11个电势向量ULINSPACEU0,U0,11,等高线命令CONTOURX,Y,U,U,K用黑实线,画出等势线如图2所示,其中,左边从里到外的第6条包围负电荷的等势线为零势线。平面电力线的画法利用MATLAB的箭头命令,可用各点的电场强度方向代替电力线。根据梯度可求各点的场强的两个分量EX,EYGRADIENTU,合场强为ESQRTEX2EY2。为了使箭头等长,将场强EXEX/E,EYEY/E归一化,用箭头命令QUIVERX,Y,EX,EY可标出各网点的电场强度的方向,异号点电荷对的场点方向如图3所示。为了画出连续的电力线,先确定电力线的起点。电荷的半径可取为R00002,如图4所示,假设第一条电力线的起始角为30度,其弧度为Q30PI/180,起始点到第一个点电荷的坐标为X1R0COSQ,YR0SINQ,到第二个点电荷的坐标只有横坐标X22AX1不同。用前面的方法可求出该点到两个电荷之间的距离R1和R2,从而计算场强的两个分量以及总场强EXQ1X1/R13Q2X2/R23,EYQ1Y/R13Q2Y/R23,ESQRTEX6EXEY6EY。下面只要用到场强分量与总场强的比值,在计算场强分量时没有乘以静电力常量K。由于电力线的方向与场强的切线方向相同,取线段为S00001,由此可求出终点的坐标为X1X1SEX/E,YYSEY/E,从而计算X2。以终点为新的起点就能计算其他终点。当终点出界时或者到达另一点电荷时,这个终点可作为最后终点这种计算电力线的方法称为切线法。源程序一对电荷平面等势线和电场线图CLEARALLCLF平面等势线的画法Q12E9Q21E9A002到原点的距离XX0005YY0004K9E9XLINSPACEXX0,XX0,20YLINSPACEYY0,YY0,50X,YMESHGRIDXR11SQRTXX0/17A2YY0/172R22SQRTXX0/17A2YY0/172R1SQRTXA2Y2各点到点电荷的距离R2SQRTXA2Y2UKQ1/R1KQ2/R2各点的电势U0KQ1/R11KQ2/R22ULINSPACEU0,U0,11取21个等势向量CONTOURX,Y,U,U,KHOLDONGRIDONPLOTA,0,O,MARKERSIZE,12PLOTA,0,O,MARKERSIZE,12XLABELX,FONTSIZE,16YLABELY,FONTSIZE,16平面电力线的画法EX,EYGRADIENTUESQRTEX2EY2EXEX/EEYEY/EHOLDONQUIVERX,Y,EX,EYTITLE一对不相等的电荷的等势线图和电场线图,FONTSIZE,20CLEAR图3源程序一对电荷平面等势线和电场线图CLEARALLCLF平面等势线的画法Q11Q21A002XX0005YY0004K9E9XLINSPACEXX0,XX0,20YLINSPACEYY0,YY0,50X,YMESHGRIDXR11SQRTXX0/17A2YY0/172R22SQRTXX0/17A2YY0/172R1SQRTXA2Y2R2SQRTXA2Y2UKQ1/R1KQ2/R2U0KQ1/R11KQ2/R22ULINSPACEU0,U0,11CONTOURX,Y,U,U,KHOLDON平面电力线的画法EX,EYGRADIENTUESQRTEX2EY2EXEX/EEYEY/EDTH120TH1DTH1DTH1180DTH1PI/180R0A/5X1R0COSTH1AY1R0SINTH1STREAMLINEX,Y,EX,EY,X1,Y1STREAMLINEX,Y,EX,EY,X1,Y1QABSQ1/Q2DTH2DTH1/QTH2180DTH2DTH2DTH2PI/180X2R0COSTH2AY2R0SINTH2STREAMLINEX,Y,EX,EY,X2,Y2STREAMLINEX,Y,EX,EY,X2,Y2GRIDONPLOTA,0,O,MARKERSIZE,12PLOTA,0,O,MARKERSIZE,12XLABELX,FONTSIZE,16YLABELY,FONTSIZE,16TITLE一对点电荷的电场分布图CLEARXY一一一一一一一一一一一00500400300200100010020030040050050040030020010001002003004005图420221510500511521050510一一一一一一一一3D一图54实验心得本次电磁场实验是利用MATLAB模拟点电荷电场的分布，刚收到实验指导书时，并不知道该怎么做，由于我们并没有正式学过MATLAB，只是在部分课程如信号，自控等课上对该软件有所接触。接到实验指导书后，我去图书馆借阅了有关MATLAB基础的书籍，了解了其中基础部分以及和实验有关的部分，并在网上搜约了相关资料，最后在同组同学的努力以及老师的指导下终于成功完成了实验。通过这次试验，我了解了MATLAB的基本操作以及其中部分函数的应用，学会用利用MATLAB模拟点电荷电场的分布。另外，我还体会到自学的重要性，大学有很多东西需要自学，只有通过自学更多的知识才能更好的应用所学的课程。电磁场理论实验二利用MATLAB模拟带电粒子在磁场中的运动1实验目的1理解数值模拟研究物理问题的思路，能独立地运用此方法研究物理问题，掌握数值模拟的编程。2运用MATLAB数值模拟的方法研究三维空间中带电粒子在复杂磁场环境下的运动行为。3实验原理带电粒子在磁场中运动时会受到洛伦兹力的作用，且随着初始运动方向和磁场分布的不同，其运动轨迹会发生不同的变化。由洛伦兹力的推导公式可知，它垂直于粒子的运动速度，不对运动粒子作功，只改变其运动方向，其大小为SINBQF；因此，综合牛顿运动定律就可以精确确定带电粒子在磁场中的运动轨迹。4实验内容1用MATLAB数值模拟的方法模拟带电粒子在恒定磁场中的螺旋运动，即带电粒子进入磁场的方向与磁场方向的角度（90）。2用MATLAB数值模拟的方法模拟磁聚焦现象，即在均匀磁场中某点引入一发散角不大的带电粒子束，并使束中粒子的速度V大致相同。3有兴趣的同学可以尝试模拟磁镜现象，即从带电粒子束进入方向，磁场逐渐增强。5实验步骤一1带电粒子在均匀稳定电磁场中受力分析BVQAM2带电粒子在均匀稳定电磁场中的运动微分方程为DTR2可将上式分解在直角坐标系展成标量式22XYZYZBDTTMQDTZYTDTTXD令1MQBX,2Y,3BZ则化简为DTYTXDTZYTDZTTXD1231223令DTXY21DTY43DTZY65则得出可以用MATLAB数值积分的一次微分方程组4126521634624321YDTYYTDYYTY3根据上述方程进行MATLAB编程建立微分方程函数实验微分方程电磁场中带电粒子FUNCTIONYDOTMF1T,Y,FLAG,Q,M,B1,B2,B3YDOTY2QB3Y4/MQB2Y6/MY4QB3Y2/MQB1Y6/MY6QB2Y2/MQB1Y4/M设置各参数的初值,并在COMMANDWINDOWS中输入相关命令,B10B20B32C0,5,0,9,0,8Q16E2M002FIGURESTRD1EXNEQ1,EYNEQ1,EZNEQ1,BXNEQO,BYNEQO,BZNEQ1T,YODE23MF1,0000120,C,Q,M,B1,B2,B3TITLESTRD1,FONTSIZE,12,FONTWEIGHT,DEMIXLABELXYLABELYZLABELZVIEW51,18COMET3Y,1,Y,3,Y,5PLOT3Y,1,Y,3,Y,5GRIDON程序中利用了ODE23求解数值微分,事实上,如果利用符号积分计算运动轨迹,由于计算机速度的限制,得不到结果利用COMET3绘制3D动态图,利用PLOT3绘制3D静止图实验结果如图505101510505050100150200二一个带电粒子进入磁场时速度的方向不与磁场垂直,则可将入射速度分解为沿着磁场方向的速度V1和垂直磁场的速度V2。在垂直磁场方向,由于粒子受到洛伦兹力的作用,做圆周运动,运动周期为QBM2T,粒子平行于磁场方向的分速度不受磁场的影响,因而将具有沿磁场方向的匀速运动。上述两种运动的合成是一个沿磁场方向的螺旋运动,这一螺旋轨迹的螺距为HV1T。如果在均匀磁场中某点引入一发散角不太大的带