基于重心有理插值函数的微分求积法及其应用研究-计算数学专业论文.docxVIP

摘妥 论文题目z 基于重心有理插值函数的微分求积法及其应用研究 学科专业 z 计算数学 研究生z 李成澄 指导教师E 赵凤群教授  签名z 李戒渔 签名z 去;t 叫层 摘要 在微分方程的诸多的数值方法中，微分求积法 (DQM)具有原理简单、计算精度高、程 序实施方便等优点，这使得 DQ 法在许多领域得到了广泛应用。但是 DQ 法的应用具有局 限性，当节点取的较多时，系数矩阵会出现病态以及数值不稳定等缺陷。为了改善传统 DQ 法的缺陷，本文提出了一种基于重心有理插值函数的微分求积法。该方法采用高阶重 心高理插值函数作为基函数，通过数值算例表明该方法比传统微分求积法计算精度高，收 敛性好，同时具有实际的应用价值。下面将本文的主要工作介绍如下: 1、根据DQ 法的基本思想，选择重心有理插值函数作为 DQ 法的插值基函数，形成 了求解一维边值问题的重心有理插值微分求积法。该方法具有适用性强，计算精度高，数 值稳定性好的优点。以二阶变系数常微分方程两点边值问题为例，运用直接法处理边界条 件，验证了该方法的有效性。并将该方法用于求解 高阶 (囚阶)常微分方程，配合δ 法处 理边界条件，仍然取得了令人满意的效果。 2、将所形成的重心有理插值DQ 法推广到二维，在二维函数中，对 x 方向和 y 方向 分别用重心有理插值 DQ 法离散并与有限差分法结合，形成了求解二维抛物型问题的重心 有理插值 DQ 法。将该方法应用到求解变系数抛物型初边值问题和 Burgers 方程组的初边 值问题，通过数值实验与现有文献的数值结果比较，验证了该方法的优越性。 3、将本文方法应用到轴向运动梁的稳定性分析中。由D Alembert 原理，建立了轴 向运动蒙古弹性变截面梁的运动微分方程，采用重心有理插值 DQ 法求解。对于简支黠弹性 变截面梁，用该方法得到了特征方程，获得了变截面梁前两阶无量纲复频率与无 量纲轴向 运动速度的变化关系。分析了梯形截面梁和抛物形截面梁随轴向运动速度变化的失稳形 式，并与等截面梁进行了比较，同时分析了变截面梁的高度比和茹弹性系数对梁动力稳定 性的影响。 关键词E 微分方程: 微分求积法:重心有理插值基函数:变截面梁;稳定性 西安 西安理工大学硕士学位论文 ii ii Abstract Title: THE DQ METHOD BASED ON BARYCENTRIC RATIONAL INTERPOLATION FUNCTION DND ITS APPLlCATION RESEARCH Major: Computational Mathematics Name: Chengcheng Li  Signature:.QJ8flJeJ四li Supervisor: Prof. Fengqun Zhao Sig驯n…n :句句ω Z也 Abstract Among the many numerical methods of differential equation，differential quadrature method(DQM) has many advantages ，such as high calculation precision ，convenient program implementation，which makes the DQ method has been widely used in many fields. While 也e DQ method also has its disadvantages ，such as numerical instability and coefficient matrix will be abnormal when it takes more nodes. In order to solve these problems，a DQ method based on rational interpolation function is proposed by selecting high order barycen创c rational interpolation function to be as basis function in this paper，and numerical experiments show 由at the proposed method not only has higher calculation precision and good convergence than traditional method，but also has the actual application v