物理学毕业论文-波动方程的数值解法要点分析.doc

学 院 本科毕业论文（设计） 题 目 波动方程的数值解法研究 院 系 物理与电子工程学院 专 业 物理学 姓 名 学 号 学习年限 2007年9月至2011年7月 指导教师 申请学位 理学 学士学位 二零一一年五月二十三日 波动方程的数值解法研究 摘 要: 本文利用有限差分法及其显示差分格式表示了波动方程的差分格式以及初始、边界条件的差分格式, 进而讨论了不同初始、边界条件下的波动方程的解的情况, 并利用MATLAB数学软件编程绘制出了各种解的图像. 为了方便对振动过程的进一步研究, 还给出了振动过程的动画演示. 对于方程初始、边界条件比较复杂因而无法求得解析解的波动方程, 有限差分方法的数值解法给出了一种较为普适而精确的计算方法. 关键词: 波动方程; 有限差分法; 显式差分格式; 数值解  Numerical Solution of Wave Equation ABSTRACT: In this paper, using finite difference method and show difference scheme showed that the difference scheme for wave equation and the initial and boundary conditions, then discussed the different initial and boundary conditions of the wave equation under the circumstances, and using mathematical software program gives MATLAB images of each solution. In order to facilitate the further research of vibration process, also gives a vibration process animation. For more complex boundary conditions and therefore should not be obtained by analytic solution wave equation, provides a relatively simple method. KEYWORDS：Wave equation; finite difference method; show difference scheme; numerical solution  目 录 引 言 1 1 有限差分法 2 1.1 有限差分法 2 2 显式差分格式 5 2.1 波动方程的差分格式 5 2.2 初始、边界条件的差分格式 5 3 计算示例 8 3.1 边界条件为零, 初始位移、初始速度不为零 8 3.2 初始条件为零, 边界条件不为零 12 3.3 自由项不为零 17 3.4 自由项、定解条件都不为零 21 总 结 26 致 谢 26 参考文献 26  引 言 在实际工作中,经常涉及求解常微分方程边值问题、偏微分方程边值问题. 通常,方程的精确解难以求解,而只能得到求解区域内某些点处解的近似值, 即数值解. 在数学上, 二阶线性偏微分方程分为双曲型、抛物型、椭圆型三种, 其中双曲型方程的典型代表是波动方程. 当满足合适的边界条件和初始条件时, 波动方程是可以求解的. 常用的求解方法有分离变量法[1,2,3,5,6,7]、行波法[1,3,5,6]、格林函数法[1,3,4]、积分变换法[1,2,3]等等. 其中最为常用的是分离变量法. 而当所讨论的问题的边界条件较为复杂时, 分离变量法将不再便于计算. 因此, 对于不易求解解析解的波动方程的数值解的研究和发展具有重要意义. 有限差分法是一种较为成熟的数值解法, 它的主要思想是将连续的问题离散化. 是将每一处导数由有限差分近似公式代替, 从而把求解偏微分方程的问题转换成求解代数方程的问题的方法. 本文在介绍了有限差分法的基础知识后给出了波动方程及其初始、边界条件的差分格式,列举了一些不同定解条件下的波动方程的数值求解方法及图像.  1 有限差分法 有限差分法 有限差分法是在采用数值