Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userIdentifying source term in the subdiffusion equation with L^2-TV regularization
117
0
0.0
(
0
)
Added by
Bin Fan
Publication date
2021
fields
Informatics Engineering
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
In this paper, we consider the inverse source problem for the time-fractional diffusion equation, which has been known to be an ill-posed problem. To deal with the ill-posedness of the problem, we propose to transform the problem into a regularized problem with L^2 and total variational (TV) regularization terms. Differing from the classical Tikhonov regularization with L^2 penalty terms, the TV regularization is beneficial for reconstructing discontinuous or piecewise constant solutions. The regularized problem is then approximated by a fully discrete scheme. Our theoretical results include: estimate of the error order between the discrete problem and the continuous direct problem; the convergence rate of the discrete regularized solution to the target source term; and the convergence of the regularized solution with respect to the noise level. Then we propose an accelerated primal-dual iterative algorithm based on an equivalent saddle-point reformulation of the discrete regularized model. Finally, a series of numerical tests are carried out to demonstrate the efficiency and accuracy of the algorithm.
rate research
Read More
In this article, we are concerned with the analysis on the numerical reconstruction of the spatial component in the source term of a time-fractional diffusion equation. This ill-posed problem is solved through a stabilized nonlinear minimization system by an appropriately selected Tikhonov regularization. The existence and the stability of the optimization system are demonstrated. The nonlinear optimization problem is approximated by a fully discrete scheme, whose convergence is established under a novel result verified in this study that the $H^1$-norm of the solution to the discrete forward system is uniformly bounded. The iterative thresholding algorithm is proposed to solve the discrete minimization, and several numerical experiments are presented to show the efficiency and the accuracy of the algorithm.
We consider $L^1$-TV regularization of univariate signals with values on the real line or on the unit circle. While the real data space leads to a convex optimization problem, the problem is non-convex for circle-valued data. In this paper, we derive exact algorithms for both data spaces. A key ingredient is the reduction of the infinite search spaces to a finite set of configurations, which can be scanned by the Viterbi algorithm. To reduce the computational complexity of the involved tabulations, we extend the technique of distance transforms to non-uniform grids and to the circular data space. In total, the proposed algorithms have complexity $mathscr{O}(KN)$ where $N$ is the length of the signal and $K$ is the number of different values in the data set. In particular, the complexity is $mathscr{O}(N)$ for quantized data. It is the first exact algorithm for TV regularization with circle-valued data, and it is competitive with the state-of-the-art methods for scalar data, assuming that the latter are quantized.
We consider the Vlasov-Fokker-Planck equation with random electric field where the random field is parametrized by countably many infinite random variables due to uncertainty. At the theoretical level, with suitable assumption on the anisotropy of the randomness, adopting the technique employed in elliptic PDEs [Cohen, DeVore, 2015], we prove the best N approximation in the random space breaks the dimension curse and the convergence rate is faster than the Monte Carlo method. For the numerical method, based on the adaptive sparse polynomial interpolation (ASPI) method introduced in [Chkifa, Cohen, Schwab, 2014], we develop a residual-based adaptive sparse polynomial interpolation (RASPI) method which is more efficient for multi-scale linear kinetic equation, when using numerical schemes that are time-dependent and implicit. Numerical experiments show that the numerical error of the RASPI decays faster than the Monte-Carlo method and is also dimension independent.
We propose a set of iterative regularization algorithms for the TV-Stokes model to restore images from noisy images with Gaussian noise. These are some extensions of the iterative regularization algorithm proposed for the classical Rudin-Osher-Fatemi (ROF) model for image reconstruction, a single step model involving a scalar field smoothing, to the TV-Stokes model for image reconstruction, a two steps model involving a vector field smoothing in the first and a scalar field smoothing in the second. The iterative regularization algorithms proposed here are Richardsons iteration like. We have experimental results that show improvement over the original method in the quality of the restored image. Convergence analysis and numerical experiments are presented.
We show stability of the $L^2$-projection onto Lagrange finite element spaces with respect to (weighted) $L^p$ and $W^{1,p}$-norms for any polynomial degree and for any space dimension under suitable conditions on the mesh grading. This includes $W^{1,2}$-stability in two space dimensions for any polynomial degree and meshes generated by newest vertex bisection. Under realistic but conjectured assumptions on the mesh grading in three dimensions we show $W^{1,2}$-stability for all polynomial degrees. We also propose a modified bisection strategy that leads to better $W^{1,p}$-stability. Moreover, we investigate the stability of the $L^2$-projection onto Crouzeix-Raviart elements.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
