Error analysis in the approximation of Sturm-Liouville spectral density functions

DOI10.1006/JMAA.1996.0394MaRDI QIDQ1815460zbMATH OpenOpenAlexFDO

Authors Charles T. Fulton, Steven A. Pruess

Publication date 9 December 1996

Published in Journal of Mathematical Analysis and Applications (Search for Journal in Brave)

Full work available at URL https://doi.org/10.1006/jmaa.1996.0394

zbMATH Keywords

asymptotic behaviorsingular Sturm-Liouville problemspectral function

Mathematics Subject Classification ID

Eigenvalues, estimation of eigenvalues, upper and lower bounds of ordinary differential operators (34L15) Numerical solution of eigenvalue problems involving ordinary differential equations (65L15)

Review text

An analysis is made of the order of error in a method of numerical approximation to the spectral function for the singular Sturm-Liouville (SL) problem on a half infinite interval. Let \(\rho(t)\) be the spectral function of the singular problem on \([a,\infty)\) and let \(\rho_b(t)\) be the spectral function for the regular truncated problem defined on \([a,b]\). \(\widehat \rho_b(t)\) will be a computed approximation to \(\rho_b(t)\) and \(\widehat I\widehat \rho_b(t)\) an approximation to \(\rho(t)\) obtained by interpolation of \(\widehat \rho_b(t)\). We may write (i) \(|\rho(t)-\widehat I\widehat \rho_b(t)|=E_1(t)+E_2(t)\) where \(E_1=|I\rho_b(t)-\widehat I\widehat \rho_b(t)|\), \(E_2=|\rho(t)-I\rho_b(t)|\) and \(I\rho_b(t)\) is the interpolation scheme applied to \(\rho_b(t)\). Let \(h_i^\ell\), \(i=0,1,\dots,2^\ell N\), \(\ell=0,1,\dots,\) be mesh lengths associated with a sequence of meshes on \([a,b]\) which are generated recursively by bisection, \(H_\ell=\max_i h_i^\ell\). The authors prove that if the coefficients of the SL operators are \(C^4[a,b]\) then (i) \(|\rho_b(t)-\widehat\rho_b(t)|=O(H_\ell^2)\), \(\ell\to \infty\), (ii) \(E_1(t)=|I\rho_b^\ell(t)-\widehat I\widehat\rho_b^\ell(t)|= O(H_\ell^2)\), \(\ell\to\infty\), where \(I\) is one of three possible interpolation procedures. Asymptotic behavior of \(E_2(t)\) is studied for three SL systems. For the first no asymptotic formula is obtained. In the second \(E_2(t)=O(1/b)\), \(b\to \infty\), and in the third \(E_2(t)=O(1/\text{ln}^2b)\), \(b\to \infty\). The analysis applies to the author's package SLEDGE used for computation of spectral functions [ACM Trans. Math. Software, 19, 360-376 (1993)].

Recommendations

• Computing the spectral function for singular Sturm--Liouville problems
• The computation of spectral density functions for singular Sturm-Liouville problems involving simple continuous spectra
• Efficient calculation of spectral density functions for specific classes of singular Sturm-Liouville problems

Cited in

(8)

• Computing the spectral function for singular Sturm--Liouville problems
• Multilayer heat equations and their solutions via oscillating integral transforms
• Efficient calculation of spectral density functions for specific classes of singular Sturm-Liouville problems
• Error analysis of a spectral projection of the Benjamin-Ono equation
• scientific article; zbMATH DE number 645829 (Why is no real title available?)
• scientific article; zbMATH DE number 1498438 (Why is no real title available?)
• Computing Spectral Measures of Self-Adjoint Operators
• Estimating spectral density functions for Sturm-Liouville problems with two singular endpoints

Describes a project that uses

Uses Software

• SLEDGE

This page was built for publication: Error analysis in the approximation of Sturm-Liouville spectral density functions

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q1815460)