A QSC method for fractional subdiffusion equations with fractional boundary conditions and its application in parameters identification
Jun Liu,Hongfei Fu,Jiansong Zhang
DOI: https://doi.org/10.1016/j.matcom.2020.02.019
IF: 3.601
2020-08-01
Mathematics and Computers in Simulation
Abstract:<p>A quadratic spline collocation (QSC) method combined with <span class="math"><math>L1</math></span> time discretization, named QSC-<span class="math"><math>L1</math></span>, is proposed to solve fractional subdiffusion equations with artificial boundary conditions. A novel norm-based stability and convergence analysis is carefully discussed, which shows that the QSC-<span class="math"><math>L1</math></span> method is unconditionally stable in a discrete space–time norm, and has a convergence order <span class="math"><math>O(τ2−α+h2)</math></span>, where <span class="math"><math>τ</math></span> and <span class="math"><math>h</math></span> are the temporal and spatial step sizes, respectively. Then, based on fast evaluation of the Caputo fractional derivative (see, Jiang et al., 2017), a fast version of QSC-<span class="math"><math>L1</math></span> which is called QSC-F<span class="math"><math>L1</math></span> is proposed to improve the computational efficiency. Two numerical examples are provided to support the theoretical results. Furthermore, an inverse problem is considered, in which some parameters of the fractional subdiffusion equations need to be identified. A Levenberg–Marquardt (L-M) method equipped with the QSC-F<span class="math"><math>L1</math></span> method is developed for solving the inverse problem. Numerical tests show the effectiveness of the method even for the case that the observation data is contaminated by some levels of random noise.</p>
mathematics, applied,computer science, software engineering, interdisciplinary applications