麦克斯韦方程组的威名可谓如雷贯耳,是经典电磁学的最高成就。先来看看麦克斯韦方程组的样子:

方程组写成以上形式,需要一个前提:空间中的媒质是各向同性的(即,媒质中的每一点的物理性质不随方向改变)。

方程组的一式描述了电荷产生电场的高斯定律,二式描述了变化的磁场产生电场的法拉第电磁感应定律,三式描述了磁单极子不存在的高斯磁定律,四式描述了电流和变化的电场产生磁场的麦克斯韦-安培环路定律。

有物理意义上的场必定存在物理意义上的源,因此,我们可以发现,麦克斯韦方程组等式的左边是场,右边是源。由于,变化的电场产生磁场,变化的磁场产生电场;又由于,只存在电荷而不存在磁荷(即不存在磁单极子,或者理解为,电力线是发散的,而磁力线是闭合的,使得具有终始方向的磁感应强度B与闭合曲面S的积分总为0);所以,根据方程组,只要在已知电荷分布和电流分布(即运动电荷)的情况下,就可以得到电场和磁场的唯一分布。

在如今物理学的教材中,标准的麦克斯韦方程组都是以上述四个积分形式出现的,不过,方程组的最初形式却并不是这样。麦克斯韦本人当年(1865年)在论文《电磁场动力论》中写下的是20个方程,都是分量形式,不是矢量形式;上述的这四个矢量积分形式的方程组是由另一位物理学家赫兹在1890年写出的[1]。

关于麦克斯韦方程组,物理学教材都告诉学生:麦克斯韦从他的方程组中推导出了电磁波的存在,并同时做出两项预言,电磁波的传播速度可以达到光速、光是一种电磁波,后来被赫兹用实验所证实。虽说很多人都知道这一事实,但很少有人去思考:如何从麦克斯韦方程组推导出电磁波的存在。

要从麦克斯韦方程组推导出电磁波,最关键的一个步骤是:采用什么方法,把方程组的积分形式演变到微分形式。令人感到有些奇怪的是,对于这个关键的步骤,几乎所有的物理学教材中都没有展示,而这正是写下本文的意义所在。

为解决这个问题,需要引进数学上的矢量分析。

大家都知道,在物理学里,E和B代表着电场强度和磁感应强度(想一下,为什么不把B对称地叫做磁场强度?然后去翻查物理史资料,会发现里面的历史其实很有趣,也可以顺便了解E和B之间的联系),所以就理所当然地认为麦克斯韦方程组中的E和B就是矢量电场强度和矢量感应强度。其实这样的理解并不完全正确,方程组中的E和B表示的是矢量场,即:空间中任意一点(x,y,z),在时刻t的电场强度(或磁感应强度)的大小和方向,用数学的方式表示:

要研究一个矢量场的性质,可完全用散度和旋度来表明[2]。因此,为了解矢量场E、B,我们需要知道E、B的散度和旋度。

先来看散度。所谓散度,是指:在一个矢量场A中,空间上存在一个点p,围绕该点作一闭合曲面S(曲面S上每一处的法线方向都是向外),矢量场A与该曲面S的积分定义为矢量场A穿过闭合曲面S的通量,再设闭合曲面S包围的体积为ΔV,当ΔV趋向于0时(此时S的也趋向于0,记为dS),我们把矢量场A穿过闭合曲面dS的通量与ΔV的比值称为矢量场A在p点的散度[3],散度的符号记为div,用数学公式表示如下:

散度是一个标量,没有方向,在一个矢量场A中,空间每一点的散度构成一个散度标量场。在直角坐标系中,散度的运算表达式为:

根据上面两式作一番变换:

这个结果公式也称为高斯-奥斯特罗格拉茨基公式(也可称为矢量场的散度定理)。有了这个公式之后,将此式与麦克斯韦方程组中的一式三式比较,可以容易得到:

再来看旋度。比起散度,旋度的概念有些复杂。先给出矢量场的环量面密度概念:在一个矢量场A中,空间上存在一个点p,围绕该点作一闭合回路L,回路的绕行方向和回路所在面的法线方向符合右手螺旋关系,矢量场A与该回路L的积分定义为矢量场A经过闭合回路L的环量,再设闭合回路L所包围的面积为ΔS,当ΔS趋向于0时(此时L也趋向于0,记为dl),我们把矢量场A经过闭合回路dl的环量与ΔS的比值称为矢量场A在p点沿ΔS法线方向的环量面密度[4],用数学公式表示如下:

由此可知,环量面密度是一个与方向有关的概念。有了环量面密度的概念后,现在可以给出旋度的定义了:在一个矢量场A中,空间上有一点p,如果在该点处存在一个矢量P,使得矢量场A在点p处沿着矢量P方向的环量面密度为最大,而这个最大值正好等于矢量P的绝对值,我们就把矢量P称为矢量场A在点p处的旋度,旋度的符号记为rot,用数学公式表示[5]:

旋度是一个矢量,有方向也有数值,在一个矢量场A中,空间每一点的旋度构成了一个旋度矢量场。那么,在知道了一个矢量场A的情况下,如何求得矢量场A的旋度呢?在直角坐标系下,旋度的运算表达式为:

我们现在再来回顾一下旋度的定义,其中最关键的是:在点p处,如何找到该处环量面密度最大的方向。结合环量面密度的定义公式来看,如果dl的方向处处与矢量场A相同,则求得的环量面密度即为最大,因此,在点p处,根据矢量场A的方向,选择合适的闭合回路dl与绕行方向,就可求得该点的旋度,我们将环量面密度的定义式作一番变换:

上式中的上标M表示:选择了合适的闭合回路dl之后,计算出来的该处最大的环量面密度。根据上式,再进行推导:

这个结果公式也称为斯托克斯公式(也可称为矢量场的旋度定理)。有了这个公式之后,将此式与麦克斯韦方程组中的二式四式比较,可以容易得到:

是否还记得前面的两个微分形式?我们现在已经得到了四个微分形式:

这四个微分形式正是我们非常熟悉的麦克斯韦方程组的微分公式,文章开头列出来的是积分公式。

现在我们从这四个微分公式推导出电磁波方程。

在一个自由的各向同性的媒质空间中(即,没有电荷分布和电流分布,ρ=0,矢量场J=0),麦克斯韦方程组的微分形式就可以写成:

根据以上四式和矢量分析的运算法则,可以得到:

同样的方法,可以得到

这两个结果式再进行整理:

最后可得

以上两式就是电磁波方程,可以看出,电磁波方程的形式和经典物理学的机械波的波动方程完全一致:

波动方程中的u是波函数,表示媒质上各点(x,y,z)在t时刻的位移,v是波的传播速度。波动方程是物理学里最基本的方程之一,任何物理量对时间和空间坐标的关系只要满足上式,该物理量就按波的形式传播。如果想知道波动方程是怎么推导出来的,可以阅读梁昆淼的《数学物理方法》一书,里面对波动方程的来龙去脉有详尽的解释。

电磁波方程中的c是光速,带有下标r的希腊字母ε、μ分别表示:介质的相对介电常数、介质的相对磁导率,现代物理学认为在真空中,两者的数值都等于1,在其他介质中都大于1,换而言之,真空中的光速大于任何介质中的光速。变化的电场在邻近的空间产生磁场,这个变化的磁场在较远的空间又会产生新的电场,电场磁场相互交变产生,由近及远向周围的空间传播出去,形成电磁波。

通过计算,麦克斯韦预言,电磁波的传播速度等于光速,又断言光是一种电磁波,也就是说,光波的波动方程就是电磁波方程。1862年,麦克斯韦在论文《论物理力线之第三部分-分子涡旋理论应用于静电学》中写到:光是某种媒质的横向波动,这种媒质正是产生电磁现象的同一种媒质;1868年,麦克斯韦在论文《关于光的电磁理论》中明确提出:光是一种按照电磁规律在场内传播的电磁扰动[7]。

麦克斯韦这位物理学大师,因身体健康原因于1879年去世,年仅48岁,假如他能再活30年,也许他能解决以太之谜、解释迈克尔逊-莫雷实验,如此,后来也不会出现洛伦兹变换、洛伦兹协变,20世纪的物理学可能会走上一条完全不同于现在的道路。当然,历史不能假设。

matlab麦克斯韦电磁方程组,从麦克斯韦方程组到电磁波动方程相关推荐

  1. UA PHYS515 电磁理论I 麦克斯韦方程组基础5 电动力学的四类问题与对应的麦克斯韦方程

    UA PHYS515 电磁理论I 麦克斯韦方程组基础5 电动力学的四类问题与对应的麦克斯韦方程 Electrostatics Magnetostatics Electromagnetic wave 一 ...

  2. UA PHYS515 电磁理论I 麦克斯韦方程组基础4 介质中的麦克斯韦方程

    UA PHYS515 电磁理论I 麦克斯韦方程组基础4 介质中的麦克斯韦方程 推导介质中的麦克斯韦方程 电位移矢量与辅助磁场强度 推导介质中的麦克斯韦方程 前三讲我们介绍了真空中的麦克斯韦方程的建立, ...

  3. UA PHYS515 电磁理论I 麦克斯韦方程组基础3 麦克斯韦方程的势能形式

    UA PHYS515 电磁理论I 麦克斯韦方程组基础3 麦克斯韦方程的势能形式 定义电磁场的potential 改写Maxwell方程 上一讲我们基于实验定律导出了真空中电磁场的Maxwell方程: ...

  4. UA PHYS515 电磁理论I 麦克斯韦方程组基础1 库仑定律与毕奥-萨伐尔定律

    UA PHYS515 电磁理论I 麦克斯韦方程组基础1 库仑定律与毕奥-萨伐尔定律 Coulomb定律与电场 Biot-Savart定律与磁场 写在前面 电磁理论是物理学中最优美的理论之一,因为它完备 ...

  5. matlab麦克斯韦速率,MATLAB编程首先建立计算麦克斯韦分布律的子程序mxwm-Read.PPT...

    MATLAB编程首先建立计算麦克斯韦分布律的子程序mxwm-Read * * 第5章 综合实例 用 MATLAB 求解问题时,一般要经历建模和编程两个过程,只有在建模正确的前提下,方能得出正确的结果. ...

  6. matlab中 三种方法计算 Ax b,在MATLAB中,方程Ax=B的解可以用哪个命令求得? matlab 求助 解方程组...

    matlab中解方程组还是很方便的,例如,对于代数方程组Ax=b(A为系数矩阵,非奇异)的求解,MATLAB中有两种方法: (1)x=inv(A)*b - 采用求逆运算解方程组: (2)x=A\B - ...

  7. matlab多元方程整数解,matlab多元非线性方程组解法

    流放之路一天赚几百是真的吗 你想太多了.如果是搬砖收入的话一天能有2E就不错了.当然脸正的时候也可以一次就赚上10E.5173上1E=3+. 用MATLAB解多元非线性方程组,CSS布局HTML小编今 ...

  8. matlab二元方程区间求解,matlab求解二元方程组

    陈星似 魔法师 matlab求解二元方程组 悬赏分:0 提问时间:2010-11-30 23:29回答数:1浏览量:241问题指向:全国 t1=(q1+q2+q3+q4-q5-q6-q7)/g1/c1 ...

  9. MATLAB怎么解方程解,怎么用MATLAB求解微分方程组并画出解函数图?

    !using["XSLSF"];                //使用命名空间XSLSF //数组xArray存放x的值:ti为当前有效值的个数:tmax为ti对应的时间:tmi ...

  10. 【matlab】解方程组

    在matlab中解方程组是很方便的 例如,对于代数方程组Ax=b(A为系数矩阵,非奇异)的求解,MATLAB中有两种方法: (1)x=inv(A)*b - 采用求逆运算解方程组: (2)x=A\b - ...

最新文章

  1. 13.挂载点、文件系统、磁盘的命名、分区的命名
  2. 音视频技术开发周刊 93期
  3. Fedora 20 配置
  4. FPGA+NIOS2 Verilog+C 实现一个计算器
  5. 在浙学大学计算机基础,艺术设计-浙江大学本科生院.PDF
  6. #控制台大学课堂点名问题_你对大学生活的5大误解!看完我想静静......
  7. sql 存储过程 盲注入_一次非常规 SQL 注入(informixsql)的利用过程
  8. i-p2psearcher开源软件源码下载
  9. django解决使用DateTimeField添加、修改记录时不动态更新时间的问题
  10. java list树 广度_java树的广度优先遍历思路
  11. Flash Player9.0 跟Flash Player8.0区别
  12. 360笔试题-字符置换
  13. mac下的c语言程序开发,mac VS Code配置C语言开发环境(小白简单Clang版)
  14. printf二进制输出
  15. 手机软件测试规范(含具体用例)
  16. 前端基础之CSS盒子模型
  17. 【UOS】如何用Ventoy工具制作UOS启动盘操作
  18. Android开发者的寒冬:轮回之下,几人破茧,几人沉沦?
  19. linux是专业的gis系统,首款支持红旗Linux操作系统的GIS基础平台问世
  20. 易中天品汉代风云人物09:刘邦崛起之谜

热门文章

  1. 西门子PLC s7-1200学习之路
  2. (一)使用RED5搭建流媒体服务器
  3. idea的pom文件变灰色
  4. ERP库存管理 华夏
  5. linux ora -03113,ORA-03113解决方法
  6. 课设 stm32的8路电压表采集系统
  7. Android和ios的区别
  8. 小甲鱼 C语言 15课
  9. 小甲鱼 c语言 17课
  10. 【ANSYS命令流】模型的建立相关命令(一)