matlab 电场线公式,基于Matlab的静电场电场线和等势线(面)的数值模拟

教育教学论坛 EDUCATIONTEACHINGFORUM 2017 年 12 月第 49 期 Dec.2017 NO.49 基于Matlab的静电场电场线和等势线(面)的数值模拟 收稿日期：2017-05-27 基金项目：南京工程学院校级科研基金项目(CKJC201602) 在大学物理的教学中，除理论分析、实验研究外，计算机数值模拟也是一种很重要的研究方法。Matlab是一套高性能的数值计算和可视化数学软件，集数值分析、矩阵运算、信号处理和图形显示等功能于一体[1]。其丰富的数学函数使得用Matlab进行数值模拟时程序变得简洁易懂，不像其他高级语言那样要从底层开始编写复杂的程序[2]。在研究如电磁学等的同类问题时，只需将Matlab代码稍加改动即可，类似的代码可重复使用，减轻了编程工作量。Matlab强大的功能为大学物理教学的数值化和可视化提供了强有力的手段。用这种语言编程可以方便地用图形、图像、动画等形象直观地描述物理过程和物理现象，有助于学生对物理问题的深刻理解[3]。本文用Matlab软件在大学物理电磁学的教学中进行了一些有益的尝试，用这种软件以电偶极子、均匀带电线段和等量异号平行直线电荷的静电场为例分析了数值模拟的过程并给出了模拟的结果。 一、电偶极子的电场如图1所示为两个相距为2a的等量异号点电荷+q 和-q组成的点电荷系，称为电偶极子。由点电荷的电势公式及电势的叠加原理，可求得场点P(x，y)处 产生 的电势为V= 1 4πε0 q r1 - 1 4πε0 q r2 ，式中，r 1和r2是场点 P到两点电荷的距离r1= (x+a) 2 +y 2 姨 ，r2= (x-a) 2 +y 2 姨 。 若取a为坐标单位，且令V0= q 4πε0a ，x′= x a ，y′= y a (V0为点电荷q在原点处产生的电势)， 则 电 势 可 表 示 为 V=V0( 1 (x′+1) 2 +y′ 2 姨 - 1 (x′-1) 2 +y′ 2 姨 )，当 场点 P 到电偶极子的距离比两点电荷的间距大得多时， 即x>>a，则有r1≈r+acosθ，r 2≈r-acosθ。于是可得电偶 极子在远处产生的电势为V=V0 2x′ (x′ 2 +y′ 2 ) 3/2 ，这就是 电偶极子的等势面方程，由 此 式可模拟描绘电偶极子的三维等势面。 王红兵 (南京工程学院数理部，江苏 南京 211167) 摘要：利用Matlab软件模拟描绘了静电场的电场线和等势线(面)，以电偶极子、均匀带电线段和等量异号平行直线电荷为例，导出了电场线和等势线方程并给出了数值模拟的相应结果，使教学内容更形象和直观，也更深刻地揭示了物理现象的内涵。 关键词：大学物理；Matlab；电场线；等势线(面)；数值模拟中图分类号：G642.0 文献标志码：A 文章编号：1674-9324(2017)49-0193-03 图1 图2 由以上数值分析，在 进 行 数 值 模 拟 时 只 需 用 几 个Matlab函数即可描绘电场线和等势线。Matlab的梯度函数gradient可以直接计算电场强度的数值分量，电偶 极子等势面可由网线曲面函数mesh描绘，等势线可用等高线函数contour描绘，电场线可用流线函数stream-line描绘。数值模拟的结果如图3、图4所示。 193- - 教育教学论坛 EDUCATIONTEACHINGFORUM 2017 年 12 月第 49 期 Dec.2017 NO.49 二、均匀带电线段的电场 如图4所示为长为2a的均匀带电线段，电荷线密度为λ。在线段l处取线元