[论文解读] Spectral method for substantial fractional differential equations

Can Huang, Qingshuo Song+1 more|arXiv (Cornell University)|Aug 26, 2014

Fractional Differential Equations Solutions参考文献 17被引用 5

一句话总结

本文提出了一种针对阶数为 $\nu$ 和 $1+\nu$（其中 $0<\nu<1$）的实质分数阶微分方程（FDEs）的非多项式谱Petrov-Galerkin方法及其相关配点技术。通过将广义拉盖尔多项式扩展至 $\alpha \leq -1$，并精心选择试函数空间与测试函数空间，该方法生成了对角且条件良好的线性系统，并通过基函数选择中的参数调节实现了高收敛速率。

ABSTRACT

In this paper, a non-polynomial spectral Petrov-Galerkin method and associated collocation method for substantial fractional differential equations (FDEs) are proposed, analyzed, and tested. We extend a class of generalized Laguerre polynomials to form our basis. By a proper scaling of trial basis and test basis, our Petrov-Galerkin method results in a diagonal and thus well-conditioned linear systems for both fractional advection equation and fractional diffusion equation. In the meantime, we construct substantial fractional differential collocation matrices and provide explicit forms for both type of equations. Moreover, the proposed method allows us to adjust a parameter in basis selection according to different given data to maximize the convergence rate. This fact has been proved in our error analysis and confirmed in our numerical experiments.

研究动机与目标

• 开发一种针对尚未被充分研究的实质分数阶微分方程的高阶数值方法。
• 通过基函数与函数空间的选择，构建条件良好且对角的线性系统，以解决谱方法在FDEs中出现的病态线性系统问题。
• 在保持多项式结构的前提下，将广义拉盖尔多项式扩展至 $\alpha \leq -1$，从而实现对非多项式解的精确逼近。
• 为对流与扩散方程提供实质分数阶微分配点矩阵的显式构造。
• 通过可调参数调节基函数，实现收敛速率的可调性，并通过误差分析与数值实验加以验证。

提出的方法

• 该方法采用基于基函数 $e^{-\sigma x} x^{\alpha} L_n^{\alpha}(x)$ 的非多项式谱Petrov-Galerkin框架，其中 $L_n^{\alpha}(x)$ 为扩展至 $\alpha \leq -1$ 的广义拉盖尔多项式。
• 通过不同缩放方式构建试函数空间与测试函数空间，确保所得线性系统为对角且条件良好，从而最小化条件数的增长。
• 该方法利用广义拉盖尔多项式及其分数阶导数的性质，实现在弱形式中对实质分数阶导数的精确计算。
• 为分数阶对流与扩散方程推导出实质分数阶微分配点矩阵的显式表达式。
• 通过调节基函数中的参数以最大化收敛速率，并通过误差分析提供理论依据。
• 通过严格的误差分析与数值实验验证该方法，结果表明其具有高收敛速率与良好的条件性。

实验结果

研究问题

• RQ1 能否为实质FDEs设计一种谱Petrov-Galerkin方法，使其生成对角且条件良好的线性系统？
• RQ2 如何在保持多项式结构的同时，将广义拉盖尔多项式扩展至 $\alpha \leq -1$，并实现对解的精确逼近？
• RQ3 基函数中可调参数对实质FDEs谱方法收敛速率的影响是什么？
• RQ4 能否为对流与扩散方程中的实质分数阶导数构造显式的配点矩阵？
• RQ5 在此新框架下，所得线性系统的条件数如何随多项式阶数变化？

主要发现

• 所提出的Petrov-Galerkin方法由于试函数与测试函数空间的精心选择，导致线性系统为对角形式，显著改善了条件性并提升了计算效率。
• 该方法实现了高收敛速率，且通过调节基函数中的参数可使收敛速率最大化，该结论在误差分析中得到证明，并通过数值实验得到验证。
• 将广义拉盖尔多项式扩展至 $\alpha \leq -1$ 后，对所有实数 $\alpha$ 均产生多项式结果，与以往扩展方式存在本质不同。
• 为分数阶对流与扩散方程推导出实质分数阶微分配点矩阵的显式形式，从而实现高效的基于配点的计算。
• 由于系统被对角化，其条件数增长缓慢，与标准谱配点方法中出现的 $\mathcal{O}(N^{2\nu})$ 增长形成鲜明对比。
• 数值实验验证了理论收敛速率，并展示了通过参数调节提升精度的有效性。