电磁场第一次仿真作业探索.docx

《电磁场基础》第一次仿真作业电81井泽佳2008010621《电磁场基础》第一次仿真作业电81 井泽佳 2008010621一、截面为正方形的无限长线电荷如下图所示。设电荷面密度为2??0; 边长a = 2。x轴上有A、B、C三点，其坐标为1.5a, 5a, 10a。所有单位均取国际单位制。a）、画出电力线分布示意图；b）、采用数值方法计算A、B、C三点处的电场强度。以A点为例，说明计算步骤；c）、对于远离正方形无限长电荷的观察点，是否有简化的计算方法，以C点为例予以说明，给出简化方法的计算结果；讨论能够采用简化计算方法的条件和简化方法的精度。解：本题求解截面为正方形的无限长线电荷问题。对于无限长线电荷，其电场围绕该线电荷对称。并且对于线电荷的任意垂直于该线的截平面，电场线及等势线分布都是相同的。所以我们可以通过研究某一截平面的电场性质，来研究该无限长线电荷电场性质。故以下题解中，研究的场源电荷为正方形。(a)应用PDETool工具箱，画出电力线分布示意图：应用MATLAB工具箱中的PDE工具，可以仿真出电力线及其等势线如图所示。电力线呈围绕线电荷发散状，由图中可见，电力线与等势线关于原点对称，并且电力线垂直于等势线，越远处，等势线越稀疏，场强越小，这点符合理论预期。(b)本题要求采用数值法计算电场强度（若不要求用数值法，可以套用高斯通量定理求解更为容易）以A点为例，说明数值法计算步骤。若要应用数值法，必然要对正方形进行细分割。比如将正方形面电荷均匀地分割为n2个边长为Δa的小正方形，令每个小正方形中心坐标为xk、yk，带电量为Δq。设观察点坐标A坐标为（d,0），则每个单位小正方形，在距其距离为r 的x轴上A点x、y方向的电场强度分别为：其中：θ为两点连线与x轴负半轴夹角。由场强的叠加性质，将各面积微元形成的场强x、y轴分别进行叠加，代入题干数据，,a=2;容易得到：由于对称性，y方向的场强应当为0。因此，对于A点，d=1.5a=3，，带入（1）式得，。此即用数值法求解A点电场强度公式。（c）对于较远处的点，如C点。从Ex表达式可以看出当x、y相对d为小量时,表达式变为，x,y近似为0：其物理意义即为将正方形的电荷忽略面积，视作点电荷来处理。代入C点坐标数据，a=2,d=10a=20，简化计算的结果为E=0.005V/m。2、MATLAB模型及其求解：在将其分割为10002个小正方形由上式写出计算场强的Matlab程序：clear all;close all;clc;n=1000;a=2;da=a/n;d=[1.5,5,10]*a;x=[-((n-1)/2):((n-1)/2)]*da%每一个微元正方形中心的x,y坐标y=x;Ex=[0,0,0];Ey=[0,0,0]; %初始化for k=1:n%数值法，叠加for l=1:n dEx=0.5*(da.^2)*(d-x(k))./(y(l)^2+(d-x(k)).^2).^1.5; //EX公式 dEy=0.5*(da.^2)*y(l)./(y(l)^2+(d-x(k)).^2).^1.5; //EY公式 Ex=Ex+dEx; Ey=Ey+dEy;endendEx,Ey%显示结果MATLAB输出结果为：Ex = 0.23374 0.02010 0.00501Ey = 1.0e-019 * 0.7284 -0.0126 0.0025其中，C点E=EX值与理论计算符合的很好，并且所有EY的值都极其近似为零。说明在此数值划分（n=1000）下，结果精度已经相当高，结果可以采用。MATLAB结果整理如下：ABCEx(V/m)0.233740.020100.00501Ey(V/m)( 1.0e-019 *)0.7284-0.01260.00253、结果分析：简化计算方法的精度：调整输出结果的显示精度，我们可以看到，当取C点时简化计算法得到结果误差仅仅为=0.2%（简化算法0.005，而MATLAB为0.00501）。事实上，取点越远，该简化算法精度越高。对于精度要求不高的情形，距离10a时，可以采用简化算法求解电场。二、真空中无限长细线如下图所示，? = 2??0，a = 2，在，的范围内画出不少于10条起点和终点分别在＋τ和－τ附近的电场线。可以任意选择起点和终点。用Matlab自行编写画电场线的程序，必须说明起点和终点的定义。内容包括：（1）程序实现原理；（2）绘制出的电场线图（画出导线的示意位置，指出起点和终点）；（3）需要说明的内容（如程序精度控制，编程体会等）。解：应用老师课上给的画电场线的小软件，可以大概看出如题干所示极性的3个点电荷电场线分布情况。应该注意到小软件中存在一些问