重心插值及其在求解高阶微分方程中的应用

摘要

将计算区间采用第二类Chebyshev点离散,利用数值稳定性好、计算精度高的重心插值近似未知函数,建立未知函数各阶导数在计算节点上的微分矩阵,提出数值求解微分方程初边值问题的重心插值法.采用重心插值法将微分方程及其初边值条件离散为线性代数方程.利用微分矩阵直接计算得到未知函数在节点上的各阶导数值.数值算例表明本文方法具有计算公式简单、程序实施方便和计算精度高等优点.