一类弦振动问题的计算机辅助分析

【摘要】 利用MATLAB软件及有限差分方法求解一类有界弦振动问题，并将数值解可视化，使对该类问题的解的物理意义有了直观的认识，将数值解与解析解比较，说明所选择数值方法的可靠性和计算机辅助分析的重要性.

【关键词】 弦振动问题；MATLAB；数值分析

0引言

弦振动问题是数学物理方程中的经典问题，弦振动理论在生产、生活中有着较为广泛的应用.如对于斜拉桥、高空观光缆车运行安全和游客人身安全而言，检测并控制钢索工作过程中的张力是至关重要的，而此类问题最终归结为弦振动问题的研究.另外，工厂车间进行络纱时，可以观察到纱会以相当快的速度沿着中心轴进行退绕，达到一定条件时便会形成多节气圈，这种现象也可以通过弦振动理论加以说明.弦振动是波动的一种特殊形式，因此弦振动理论在声学中也占有举足轻重的地位，弦振动理论与现代先进的软件和硬件相结合，在许多领域诸如音乐物理学、 材料学和系统分析中也都得到了广泛的应用[1].弦振动问题是由弦振动方程附加初边值条件构成的，弦振动方程是一类重要的偏微分方程，它关于时间和空间变量的导数均是二阶的，在解决实际问题的过程中，遇到的弦振动问题通常既有初始条件又有边界条件，是一类混合问题.对于一些简单的混合弦振动问题可以通过分离变量法、叠加原理、齐次化原理及傅里叶变换等方法求其解析解或者是广义解，但是对于一些复杂的弦振动问题，传统的方法求解起来就比较困难，或者求出的解的形式相对复杂，与实际相结合较为困难.这时，就可以通过数值方法进行求解，借助计算机进行辅助分析，从而使求解过程变的轻松，并且能够将所得结果以图形等形式直观地显示出来[3].

文中考虑了以下两端固定的有界弦振动问题

1数值分析

该节首先给出求解问题（1）的有效的数值方法，然后利用这些方法及Matlab软件求出数值解，并且比较几种数值方法的优劣.

1.1数值方法

求解偏微分方程的数值方法有很多，比较常见的方法就是有限元法和差分法，文献[4]中利用有限元法对（1）进行了求解.文中将利用差分法对（1）进行求解，主要采用显式、隐式两种差分格式.

1.1.1区域剖分

首先，对定解问题（1）的求解区域W={（x，t）|0 0}进行剖分.在区间[0，1]内取n-1个等分点xi=ih，i=1，…，n-1将区间n等分，其中x0=0，xn=l，h为空间步长.同样，在时间区间内插入若干等分点tj=jτ，j=1，2，…，其中，t0=0，τ为时间步长.用两族平行直线x=xi，1

1.1.2方程离散

用差分方法求解偏微分方程的关键是对方程的离散，下面采用显示、隐式差分格式对问题（1）中方程进行离散，得到如下形式的差分方程.

（i）显示差分方程

1.2数值实验

下面将利用计算机实验对问题（1）的数值解进行辅助分析，并将其可视化，实验主要采用的软件是MATLAB软件.

首先，根据显示差分格式编制MATLAB程序模拟问题（1）的数值解，获得动态结果如图1所示. 图1在图1实验中取初始位移为

φ（x）=sin7πx（37

0其他

初始速度ψ（x）取为零，时间步长τ=0.0005，空间步长h=0.0024.通过对上述结果进行观察，发现初始状态的波平均分成两个波向相反方向传播，在它们到达弦的端点之前，和波在无限长弦上的传播是一致的，到达弦的端点以后，两个波同时发生反射，出现了半波损失，这是有别于无限长弦的振动的，接下来，被反射回来的两个波进行相向运动并在弦的中间相遇，相遇后两个波叠加继续传播，整个波的传播过程清晰明了，实现了问题（1）的数值解的可视化，有利于进一步研究更为复杂的弦振动问题.

对于问题（1）还可以根据隐式差分格式编制MATLAB程序模拟它的数值解，相应程序的编码要比显示差分格式复杂的多，涉及到三对角方程组的求解等技术性问题，具体程序如下

2结束语

以上利用MATLAB软件及有限差分方法求解两端固定的有界弦振动问题，首先利用区域转化的思想对方程进行离散，然后编制MATLAB程序进行求解，并将数值解可视化，使对问题（1）的解的物理意义有了更深刻的认识，另外，通过与解析解比较，也充分说明了所选择数值方法的可靠性和计算机辅助分析的重要性.

参考文献

[1]李韵，郭怡文，吕郁文.基于 Matlab 环境的弦振动方程的图像与音效模拟.科协论坛[J].2009（7）：74-75.

[2]谷超豪，李大潜，等.数学物理方程：第三版[M].北京：高等教育出版社，2012.

[3]孙丽男，张馥菊.基于MATLAB的一类输运问题的数值分析.哈尔滨师范大学自然科学学报，2012，28（1）：14-17.

[4]高峰，陈君若，等.有界弦振动方程的有限元方法求解[J].机械设计，2008，25（11）：15-18.Computer Auxiliary Analysis on a Class of String Vibration ProblemsSun Linan1，2，Zhang Fuju3

（1. Heihe University； 2.Dalian University of Technology；3.Harbin Normal University）Abstract：A class of bounded string vibration problem is soved by MATLAB software and finite difference methods and it makes the numerical solution visualization， makes the physical meaning of the solution of this kind of problem with intuitive understanding. Compared numerical solution with analytical solution， the reliability of the selected numerical method and the importance of computer auxiliary analysis are illustrated.

Keywords： String vibration problem； MATLAB： Numerical analysis.

（责任编辑：黄永辉）