J4 ›› 2012, Vol. 25 ›› Issue (3): 80-87.doi: 10.3976/j.issn.1002-4026.2012.03.017

• 目录 • 上一篇    下一篇

波动问题的高精度重心有理插值配点法

马燕,王兆清,唐炳涛   

  1. 马燕,王兆清,唐炳涛
  • 收稿日期:2012-02-01 出版日期:2012-06-20 发布日期:2012-06-20
  • 作者简介:马燕(1985-),女,硕士研究生,研究方向计算力学及其应用软件。

High-precison barycentric rational interpolation collocation method of wave problems

 MA Yan, WANG Zhao-Qing, TANG Bing-Chao   

  1.  MA Yan, WANG Zhao-Qing, TANG Bing-Chao
  • Received:2012-02-01 Online:2012-06-20 Published:2012-06-20

摘要:

       提出一种数值求解波动问题的高精度重心有理插值配点法。对于给定的时间和空间上的计算节点,采用重心有理插值近似未知函数,建立未知函数关于时间和空间变量导数的微分矩阵。将未知函数的重心有理插值近似函数代入波动问题的控制方程,得到波动问题方程和定解条件的离散代数方程组。利用微分矩阵的记号,将离散后的代数方程组写成简洁的矩阵形式。通过置换法施加边界条件和初始条件,求解代数方程组,得到波动问题在计算节点处的位移值。数值算例表明,重心有理插值配点法具有计算公式简单、计算节点适应性好、程序实施方便和计算精度高的优点。

关键词: 波动问题, 重心有理插值, 微分矩阵, 配点法

Abstract:

       We construct the differentiation matrices of an unknown function regarding temporal and spatial variables for the given computational nodes in temporal and spatial fields with the approaximation of barycentric rational interpolation to the function. We initially acquire discrete algebraic equations of a wave equation and its definite conditions by inserting barycentric rational interpolation of an unknown function into the governing equation of the wave equation. We then denote the discrete algebraic equations as a concise matrix with the notation of differentiation matrices. We eventually obtain the displacements of the wave equation on the nodes by replacement method and applying boundary and initial conditions. Numerical examples demonstrate that the approach has such advantages as simple computation, easy programming and high precision.

Key words: wave motion problems, barycentric rational interpolation, differentiation matrix, collocation method

中图分类号: 

  • TB115.1