Projecting interval uncertainty through the discrete Fourier transform: An application to time signals with poor precision

• Published in: Mechanical systems and signal processing Vol. 172; p. 108920
• Main Authors: Behrendt, Marco, de Angelis, Marco, Comerford, Liam, Zhang, Yuanjin, Beer, Michael
• Format: Journal Article
• Language: English
• Published: Berlin Elsevier Ltd 01.06.2022; Elsevier BV
• Subjects: Amplitudes; Complex intervals; Dependency tracking; Discrete Fourier transform; Fourier transforms; Interval arithmetic; Monte Carlo simulation; Power spectral density; Power spectral density estimation; Time signals; Uncertainty; Uncertainty quantification; Power spectral density estimation; Uncertainty quantification; Dependency tracking; Complex intervals; Discrete Fourier transform; Interval arithmetic
• ISSN: 0888-3270; 1096-1216
• DOI: 10.1016/j.ymssp.2022.108920

The discrete Fourier transform (DFT) is often used to decompose a signal into a finite number of harmonic components. The efficient and rigorous propagation of the error present in a signal through the transform can be computationally challenging. Real data is always subject to imprecision because of measurement uncertainty. For example, such uncertainty may come from sensors whose precision is affected by degradation, or simply from digitisation. On many occasions, only error bounds on the signal may be known, thus it may be necessary to automatically propagate the error bounds without making additional artificial assumptions. This paper presents a method that can automatically propagate interval uncertainty through the DFT while yielding the exact bounds on the Fourier amplitude and on an estimation of the Power Spectral Density (PSD) function. The method allows technical analysts to project interval uncertainty – present in the time signals – to the Fourier amplitude and PSD function without making assumptions about the dependence and the distribution of the error over the time steps. Thus, it is possible to calculate and analyse system responses in the frequency domain without conducting extensive Monte Carlo simulations nor running expensive optimisations in the time domain. The applicability of this method in practice is demonstrated by a technical application. It is also shown that conventional Monte Carlo methods severely underestimate the uncertainty. •A novel algorithm for propagating interval signals through the DFT is presented.•Interval arithmetic is utilised to obtain exact bounds on the Fourier amplitude.•The problem of repeated variables arising in the amplitude calculation is addressed.•It is shown that MC underestimates the uncertainty in the Fourier amplitude bounds.•The efficiency is tested against an example involving dynamic structural analysis.