浅谈使用Matlab软件作为计算物理教学平台的优势

2023-02-05 06:05:36 发布在教学论文

　　浅谈使用Matlab软件作为计算物理教学平台的优势

浅谈使用Matlab软件作为计算物理教学平台的优势

　　文林，邱旭

　　（重庆师范大学物理与电子工程学院，重庆401331）

　　摘要：Matlab为数值计算提供了大量功能强大的应用工具箱，使数值计算变得更加的高效。作为数值计算的软件平台，Matlab在计算物理教学中的应用能降低学生学习计算物理的难度，提高学生对计算物理的学习兴趣。

　　关键词：Matlab；计算物理；线性方程；BCS理论

　　一、引言

　　计算物理是伴随着计算机的出现和发展而逐步形成的物理学科，是运用计算机对所要研究的复杂物理问题进行数值计算或模拟实验，并从中探索和发现新的物理规律。因此，计算物理是物理、数学和计算机三者相结合的产物[1]。在当今的物理学各个分支领域，大量物理问题的求解都离不开计算机的辅助。可以说，计算物理的发展对物理学起着重要的推动作用。

　　在高等教育中，许多高校都为物理学专业本科生开设了计算物理课程。通过对该门课程的学习，学生可掌握许多科学计算中所应用的方法与技术，培养他们对科研的兴趣，并为将来开展科学研究和应用打下良好的基础。然而，对于大学本科生来说，在计算物理的学习中存在着许多的困难。首先，计算物理课程通常涉及较为复杂的数学知识，如线性代数、偏微分方程等。在传统的计算物理教学中，大部分教师把数值计算方法和原理作为教学重点和难点，使得高等数学知识基础较差的学生学习起来较为困难[2]。另一方面，计算物理课程的教学离不开计算软件。在早期的教学过程中，大部分教师都以Fortran作为软件平台进行程序的编写，要求学生具有较强的编程能力。但学生通常很难掌握和应用复杂的.Fortran语言，使得计算物理教学过程中存在教与学的难度，降低了学生学习的效率[3]。随着科学的进步，许多新型计算软件的诞生为计算物理教学和学生的学习提供了便利，最为典型的就是Matlab软件在计算物理教学中的应用。

　　Matlab是美国MathWorks公司所开发的商业数学软件，其全称为Matrix Laboratory，意为矩阵实验室。

　　Matlab的主要功能包括数值分析、数字图像处理、数字信号处理、工程与科学绘图和通讯系统设计与仿真等[4]。与其他数值计算软件相比，Matlab的优势主要体现在如下几方面：（1 ）Matlab为用户提供了高效的数值计算和符号计算功能，特别是在处理矩阵方面具有明显的优势；（2）Matlab在编程语言上较为简单，语法结构更加接近人们常用的数学表达习惯，提高了编程效率。最重要的是，Maltab内嵌了许多功能强大的应用工具箱，为用户提供了大量方便且实用的处理工具，使得一些烦琐的数值计算过程能用单个简单的命令代替；（3）Matlab具有完美的图像处理功能，让复杂的计算结果变得可视化。因此，如果将Matlab作为计算物理教学的软件平台，将有效地改变学生学习计算物理困难的现状，提高学生学习积极性的同时，也能降低高校教师的教学难度。本文将通过一些计算物理实例，进一步分析Matlab软件作为计算物理教学平台的优势。

　　二、Matlab在求解线性方程组中的应用

　　在物理学中，许多实际物理问题可用线性方程组来描述，如电学网络、分子振动、应力分析等问题，而这些线性方程组又可用转化成如下的矩阵形式AX=B其中X=（x1，x2，…，xN）为未知数所构成的N×1矩阵，A为N×N的系数矩阵，而B为N×1矩阵。如果矩阵A是非奇异矩阵，则方程组的解为X=A-1B。所以，非线性方程组的求解转变为求逆A-1。在线性代数中，矩阵的逆A-1＝A*/|A|，这里A*为伴随矩阵。|A|为矩阵的行列式。如果利用标准的线性代数方法求N×N矩阵A的行列式值[5]，则必须要计算N!个乘积。对于较小的矩阵来说，通常可利用简单的数学公式求解行列式。当N较大时，尽管理论上仍然可以运用线性代数公式进行计算，但计算量非常大。以N=20的矩阵为例，如果我们的计算机每秒可做1010次运算，那么计算该矩阵的行列式要花费大约为7.7年的时间。这说明，在实际的数值计算中，不能用标准的线性代数公式做计算。此外，伴随矩阵A*的求解在数值上也不太容易实现。

　　Matlab为我们提供了一个计算矩阵行列式的命令———det。该命令所使用的算法为高斯消去法，其基本思想是将矩阵化为上三角或下三角矩阵，然后将对角元素相乘即可得到矩阵的行列式。我们以N=20的随机矩阵为例，发现计算该矩阵的行列式所需要的时间约为5×10-5秒。由此可见，Matlab所提供的命令det在计算行列式时，既高效又简单。特别的，基于高斯－约当消去法，Matlab也为我们提供了计算矩阵逆的命令———inv。因此，学生在利用这些命令求解线性方程组时，不仅可以节省计算量，而且还可以避免编写高斯消去法和高斯－约当消去法时所涉及的烦琐程序。

　　三、Matlab求解非线性方程的应用举例

　　在日常生活中，许多现象在理论上可用一些复杂的非线性方程或方程组来描述。通过求解这些非线性方程，可以深入研究这些现象背后所隐藏的物理规律。非线性方程的求解没有简单的求根公式，必须借助数值计算方法进行近似求解。对于单个非线性方程的求解，常用的数值方法有区间对分法、迭代法、搜索法及牛顿迭代法[1]。尽管这些方法能够有效地求解非线性方程，但在实际的应用中，学生很难直接利用这些算法求解较为复杂的非线性方程或方程组。Matlab为我们提供了一个求解非线性方程或方程组的命令———fsolve，该命令主要采用最小二乘法来求解非线性方程[1]，其用法比较简单，学生容易学习和掌握。

　　为了说明利用Matlab求解非线性方程的优势，我们以BCS超导理论中的能隙方程和粒子数方程为例，其基本形式如下：
浅谈使用Matlab软件作为计算物理教学平台的优势

　　这里为as原子相互作用的s波散射长度，kF为费米动量，配对序参量Δ和化学势μ需要自洽求解。该方程组为非线性方程组，可以通过牛顿迭代法进行求解。牛顿迭代法的基本思想是，将两个方程在初始值附近用二元泰勒级数做展开，并取线性部分，可得到求解方程时所需要的迭代格式[1]。然而，牛顿迭代法在每一步迭代的过程中都需要计算一次雅可比行列式。对于较为复杂的非线性方程组而言，雅可比行列式的求解较为困难，程序编写过程比较复杂，以至于学生在编程过程中容易出现错误。如果直接利用Matlab的fsolve命令，则可有效避地免该问题。以下是利用fsolve求解配对序参量和化学势的Matlab程序：

浅谈使用Matlab软件作为计算物理教学平台的优势