二维区域本征值问题的研究毕业论文.doc

二维区域上本征值问题的研究 第 9 页 共 9 页二维区域上本征值问题的研究徐荣超 ( 安庆师范学院物理与电气工程学院 安徽 安庆 246011)指导教师：章礼华摘要：本文是对二维区域上的本征值问题进行研究，二维本征值问题涉及二维振动问题、输运问题以及稳定场等问题。首先写出相关的定解问题，然后导出二维本征值问题，如矩形、圆形、三角形、扇环形等区域，并运用数学物理方法教程中所涉及的方法对其进行求解，例如分离变量法，拉普拉斯变换法，得到本征值和本征函数。最后对各种二维区域的本征值问题进行总结，并进行分析比较。关键词：二维区域，本征值问题，分离变量法1 引言在运用数学物理方程来研究物理问题时，我们都必须要解数学物理方程,而求解数学物理方程的过程中必然会产生本征值问题，需要我们对本征值问题进行求解。梁昆淼的数学物理方法、姚端正的数学物理方法以及四川大学出版的高等数学第四册均详细的向我们讲解了利用分离变量法、积分变换等方法来对各种有界问题进行求解，其中包括矩形区域和圆形区域1 2 3；章礼华的等腰三角形薄膜振动的解析解在矩形膜的特征解的基础上，求解了占据等腰三角形区域的均匀膜的横振动问题，并讨论了等腰三角形的均匀膜的横振动的某些性质4；彭芳麟的数学物理方程的MATLAB求解与可视化以及计算物理基础向我们介绍了利用MATLAB指令来求解二维区域上本征值问题，并可将结果用图形甚至动画表现出来5 6。在解决二维振动方程问题、输运问题或者稳定场等问题时免不了要面临对二维本征值问题的求解，因此研究二维区域上本征值问题对于解决各种二维数学物理方程有着非常重要的作用和意义。2 本征值问题的基本概念考虑长为两端固定的弦的自由振动 (1)由经典力学的知识我们清楚地知道，两端固定的弦的自由振动会产生驻波，这就指导我们应该考虑利用驻波的一些性质来解决这样的问题。在经典力学中，驻波的表达式为 (2)对于定解问题(1)式可设其分离变数的试探解为 (3)也就是说和分别只是变量和的函数，变量和可以分开互不影响。把试探解（3）式带入定解问题方程，分离变量得到，即 （4）由（4）式可以看出，等式左边只是的函数，右边只是的函数，它们是两个相互独立的变量，因此只有当两边都是常数时，上述等式才可以成立。我们令这个常数为，即（4）式可以写成 （5）由（5）式改写得到两个常微分方程 ， ， （6）其中被称为分离常数。现将分离变数的试探解代入边界条件，即可以得到 ， （7）因为我们考虑的是非零解，所即应不恒为零，所以我们可以得到齐次边界条件。这样一来，求偏微分方程定解问题的解，即是求常微分方程满足齐次边界条件的非零解。只有当取某些特定数值时，偏微分方程定解问题才有非零解，由偏微分方程和边界条件构成的问题，即 （8）称为本征值问题，称为本征值问题的本征值，相应的非零解为本征值问题的本征函数，求本征值和本征函数的问题即为本征值问题。3 二维区域上本征值问题的研究3.1矩形区域上本征值问题0考虑边界条件为四边固定，求解区域为0，0的矩形膜的横向振动问题，如图1所示。图1.矩形薄膜的横向振动与之相对应的定解问题为 （9）将定解问题分离变量得三个方程为 (10) (11) (12)方程（10）和其对应的边界条件构成一个本征值问题 (13) 方程（11）和其对应的边界条件构成又一个本征值问题 (14)分别求解得到本征值问题（13）和（14）的本征值和本征函数分别为 (15) (16)即得到矩形膜的振动的本征模为 3.2 圆形区域上本征值问题考虑一个半径为,上下两面绝热的薄圆盘上的稳定温度分布问题，如图2所示。由于边界是圆周，因此采用极坐标求解更为简便，又因为圆盘的上下两面绝热，所以没有热量流动，又因为圆盘很薄，因此可以看成是二维圆形区域的稳定场问题。在极坐标下该问题可以表示为 （17）其中为已知的函数，表示圆盘边缘的温度函数。将极坐标中的的表达式带入定解问题（17）式中可以得到O图2.圆形的薄圆上稳定的温度分布 （18） 以分离变数形式的试探解 （19）带入（18）式，运用分离变量法得到两个方程 （20） 由于定解问题的边界条件是并不是齐次的，因此我们无法运用分离变量法对其进行求解，但是我们可以经过分类讨论来求解本征值以及与其相对应的本征函数。由于在一般的物理问题中，函数为单值函数，因此我们可以得到周期性的边界条件由（19）式即可得也就是需要满足条件这就给关于的方程提出了周期性的边界条件，因此我们得到了一组本征值问题 （21）若则本征值问题（21）的解为其中为常数，由周期性边界条件可知所以可得，即若，记，则本征值问题（21）的通解为由于实指数函数不是周期函数，显然它不满足周期条件，因此此种情况并不存在。若，则本征值问题的通解为将边界条件带入即可得到，但是由于当取正整数和负整数时本征值方程是线性相关的，故只要取综上所述，本征值问题（21）的本征值和相应的本征函数为 （22） （23）3.3等腰直角三角形区域上本征值问题 图3.等腰直角三角形均匀膜的横向振动考虑如图3所示的直角边固定斜边固定（或自由）的等腰直角三角形均匀膜作横振动的问题时，可以写出与之对应的定解问题 （24）以分离变数的试探解，分离变量可得所满足的本征值问题为 （25）式子中为常数，由于斜边上边界条件显然不能直接采取分离变量法来进一步求解，但是根据已在数学物理方法教材1 2 3中讨论过的问题，边长为的正方形区域上相对应的本征值，相应的本征函数为 （26）显然满足本征值问题（25）中的方程以及两个直角边上的边界条件，但还是满足不了斜边上的边界条件。为了找到能满足斜边上边界条件的解，我们可以假想式（26）中的相应于同一本征值的两个不同函数的组合，即，这样假想给出的显然可以满足式（25）中的方程以及两个直角边上的边界条件，并且同属于本征值。将这样给出的带入斜边上的边界条件，得到 （27）因此只需要取就可以了，因此我们可以得到本征值问题（25）的本征值和本征函数分别为， (28)， (29)3.4扇环形区域上本征值问题O图4.扇环形薄膜的横向振动考虑如图4所示的夹角为周边固定的扇环形薄膜的横向振动，为了简便，我们采用极坐标系来表示，表示内外环的半径。可以写出其对应的定解问题为 （30）采用分离变数的试探解，令代回上述方程得 （31）即可得到两个方程如下 （32） （33）由方程（33）和式（30）中的边界条件构成了本征值问题 （34）当时，（33）式的通解为 （35）由边界条件可得即解得本征值和本征函数分别为 （36） （37）3 结语综合全文来看，在矩形和等腰直角三角形区域中，涉及的本征值问题的类型为 相应的本征值，本征函数为。在圆形和扇环形区域中，定解问题采用极坐标求解，其中在圆形区域中涉及的本征值问题为相应的本征值为，本征函数为，其边界条件为周期性的边界条件,而在扇环形区域中周期性的边界条件便不再成立。参考文献1梁昆淼，刘法，缪国庆，数学物理方法(第四版)，高等教育出版社,2010。2姚端正，梁家宝，数学物理方法（第三版），科学出版社，2010。3四川大学数学系高等数学教研室，高等数学第四册（第二版），四川大学出版社，1985。4章礼华，朱德权，王其申，等腰三角形薄膜横振动的解析解，大学物理，2012（7），31-33。5彭芳麟,数学物理方程的MATLAB求解与可视化，清华大学出版社,2004。6彭芳麟,计算物理基础，高等教育出版社,2010。Research on eigenvalue problem of the two-dimensional areaXU Rong-chao (School of Physics and Electrical Engineering of Anqing Normal College, Anqing 246011)Abstract : In this paper, we make a research on eigenvalue problem of the two-dimensional area. The two-dimensional eigenvalue problem relates to vibration problem, transport problem and stable field problem. First, we have to write out relevant definite solutions problem, then find out the two-dimensional eigenvalue problem. The two-dimensional area includes the range of rectangular, circular, triangular and sector annual area. Besides, we can use the methods what had be mentioned in the Methods of Mathematical physics, such as variables separation, method of Laplace transformation. Through solving the eigenvalue equation, we obtained the eigenvalue and