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

《电磁场基础》 第一次仿真作业 电81 井泽佳 2008010621 《电磁场基础》第一次仿真作业 电81 井泽佳 2008010621 一、截面为正方形的无限长线电荷如下图所示。设电荷面密度为2??0; 边长a = 2。x轴上有A、B、C三点，其坐标为1.5a, 5a, 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轴分别进行叠加，代入题干数据，ρ=2 ?q=ρ*(?a)2 由于对称性，y方向的场强应当为0。 因此，对于A点，d=1.5a=3，?a=an=2/n 。此即用数值法求解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; end end Ex,Ey %显示结果 MATLAB输出结果为： Ex = 0.23374 0.02010 0.00501 Ey = 1.0e-019 * 0.7284 -0.0126 0.0025 其中，C点E=EX值与理论计算符合的很好，并且所有EY的值都极其近似为零。说明在此数值划分（n=1000）下，结果精度已经相当高，结果可以采用。 MATLAB结果整理如下： A B C Ex(V/m) 0.23374 0.02010 0.00501 Ey(V/m)( 1.0e-019 *) 0.7284 -0.0126 0.0025 3、结果分析： 简化计算方法的精度： 调整输出结果的显示精度，我们可以看到，当取C点时简化计算法得到结果误差仅仅为0.00501-0.0050.005×100%=0.2%（简化算法0.005，而MATLAB为0.00501）。事实上，取点越远，该简化算法精度越高。对于精度要求不高的情形，距离 二、真空中无限长细线如下图所示，? = 2??0，a = 2，在，的范围内画出不少于10条起点和终点分别在＋τ和－τ附近的电场线。可以任意选择起点和终点。用Matlab自行编写画电场线的程序，必须说明起点和终点的定义。内容包括：