Locating discontinuities of a bounded function by the partial sums of its Fourier series (Q1580937)

In problems dealing with the reconstruction of a function with jump discontinuities, the accurate approximation of the jumps and their locations are of vital importance. The present authors develop an algorithm which enables one to approximate the locations of discontinuities and the magnitudes of jumps of a bounded function given by its truncated Fourier series. The first limit relation to determine the jumps of a function \(g\) of bounded variation by means of the partial sums \(S_n(g;\theta)\) of its Fourier series is the following: \[ \lim_{n\to\infty} {S_n'(g; \theta)\over n}= {1\over\pi} [g(\theta+)- g(\theta-)], \] where \(\theta\) is arbitrary, provided that \(g\) is periodic. The authors use more sophicated limit relations. The stability and complexity of the method suggested by the authors is investigated in details, and several numerical examples are also presented.