波动方程的差分反演模型.pdfVIP

应用数学和力学，第29卷第3期 Applied Mathematics and Mechanics 2008年3月15日出版 Vo1．29，No．3，Mar．15，2008 文章编号：1000-0887(2008)03-0325-06 ⑥ 应用数学和力学编委会，ISSN 1000-0887 波动方程的差分反演模型 王德明 (哈尔滨工业大学数学系，哈尔滨150001) (王彪推荐) 摘要： 为了反演波动方程的系数函数，利用差分离散及扰动假设，推导出一个适合迭代的数值模 型．解决了以往方法中正反演模型数值精度不一致问题，以及由此带来的一系列问题．经数值模 拟计算说明，该方法是可行的和有效的． 关 键 词：反问题；波动方程； 扰动迭代 中图分类号：O241．8 文献标识码：A 引 言 波动方程反问题是指那些在正问题中的一个或几个原来的已知量，现在变为未知．而原 问题的未知量可能仍是未知的．我们只知道与这些未知量有关的某些附加信息，用方程、定解 条件与附加的某些其他条件来确定这些未知量． 在以往的反问题求解方法中，脉冲谱技术(PST)_l J是一个典型的方法之一．它的基本思想 是，为了识别某一系统，可向该系统传送某些信息，如脉冲信息，然后从接收的信息中进行分析 与综合．信息的传送与接收在时域中进行，而分析与综合则是借助于一个迭代过程在频域中 进行． 脉冲谱技术的分析与综合过程是借助于Green函数的理论，将每个迭代过程中的关于未 知参数的微分方程转化为第一类Fredholm积分方程．理论上讲，这个转化过程是严格的．然 而，Green函数的解析表达式往往是不易得到的．它们的数值求解不仅计算量大，而且精度低， 这是脉冲谱技术在应用时遇到的一个重要问题．如果迭取特殊的附加条件及未知参数的特殊 边界条件，有时也可以避开Green函数的计算_2 J，这无疑又限制了测量条件，因为所迭取的特 殊附加条件，在工程中未必是能够实现的．更何况在迭代过程中，未知参数的特殊边界条件， 直接影响迭代的收敛性．而且，微分方程的离散和积分方程的离散，其精度是不一致的．这也 是以往在方法中的计算结果只能取一步迭代的重要原因之一． 根据前人的扰动假设思想，本文在频域中构造一个扰动迭代过程，利用中心差分格式进行 离散，将波动方程用离散值的矩阵和向量模型来代替．进而利用矩阵分析原理，将模型转化为 关于迭代中未知量的数值反演模型．这样，避开了Green函数的计算，而将Green函数的计算 收稿日期： 2007—09—24；修订日期： 2008一O1—24 作者简介： 王德明(196o__)，男，黑龙江人，副教授，博士(Tel：+86-451E-mail：wangdeming@ hit．edu．cn)． 325 326 王 德 明 隐含在矩阵的转换技巧之中．这就去掉了实际工程中对澳0量条件的限制．并且保证了正、反 演模型中的数值精度的一致性，去掉了由于正、反演模型中的数值精度的不一致，给迭代带来 的较大误差的问题．这是迭代不收敛的一个重大因素． 应该指出的是，反问题往往是不适定的．数值求解时，一般需要加人一些对解的限制．例 如光滑化假设 ]等．本文中采用Phillips和Twemey[ 4]的光滑化方法求解反问题，结果是令 人满意的． 1差分反演模型 考虑波动方程(杆的振动微分方程)的定解问题 ( ) 3u)= ) ， ∈ 6]，t u(。，￡)：O， 3 u(6 ￡)= (￡)， (1) ， ( ，￡)：0， ( ，0)：0， 其中，