利用MATLAB仿真静电场.pdfVIP

电磁场基础第一次仿真报告 程序于MATLAB 7.11.0 测试通过 目录 题目一 1 问题描述 2 问题分析 2 问题的数值解决 3 画电力线分布图 4 数值法求电场强度 5 简化计算方法 6 实验结论 7 题目二 7 问题描述 7 理论分析 8 算法实现： 9 代码编写： 9 绘制方法的实现原理： 11 精度控制和编程体会 12 实验收获与体会 13 题目一 1 问题描述 截面为正方形的无限长线电荷如下图所示。设电荷面密度为20 ，边长a=2 。x 轴上有 A、B、C 三点，其坐标为1.5a, 5a, 10a。所有单位均取国际单位制。 图1 正方形无限长线电荷截面图 （1）画出电力线分布示意图； （2 ）采用数值方法计算A 、B 、C 三点处的电场强度。以A点为例，说明计算步骤； （3 ）对于远离正方形无限长电荷的观察点，是否有简化的计算方法，以C 点为例，予以说 明，给出简化方法的计算结果；讨论能够采用简化计算方法的条件。 问题分析 对于由无限长的均匀截面电荷，可以把它分解成无数无限长线电荷的集合体。从而，我们 可以认为它产生的电场是平行平面场。所以，只需要研究其在XOY 平面上的分布，就可以了解 其在整个空间的分布。那么，分别计算出每一条线电荷在其周围产生的电场，再进行叠加即可 得到整个区域内的电场分布。 ^  r 对于无限长线电荷，其场强公式E r 2 ，可以知道场强与从电荷到观察 2 r 2 r 0 0 点的矢径共线。又因为所给线电荷截面的对称性，可以得到电场在XOY 平面上的分布也是对称 的。通过定性分析，该电荷在XOY 平面上既关于x 轴和y 轴对称，同时也关于四个象限的角平 分线分别对称。由于截面为正方形，并不具有严格的旋转对称性，所以电场在靠近电荷的区域 2 分布并不均匀，具体表现为电场线不严格的放射状分布，而会具有一定的不均匀特性。但在远 离电荷的区域，可以将电荷看做一无限长线电荷，电场线沿放射状分布，可以近似认为是直 线。 欲求XOY 平面内任意点P 的电场矢量E，可以先求所有组成线电荷的“小截面”无限长线 电荷在P 点产生的电场矢量dE，再叠加即可。利用无限长线电荷的场强公式，可以知道坐标 为(x,y)的小截面对应线电荷dxdy 在P 点产生的场强 dxdy ^ r d E r 2 dxdy 2 r 2 r 0 0 设P(x ,y ),由于已知 2 ，那么 p p 0 r (xp x, yp y ) d E 2 dxdy 2 2 dxdy 2 r (x  x) ( y  y) 0 p p 将x