Projecting interval uncertainty through the discrete Fourier transform: An application to time signals with poor precision. (1st June 2022)

• Record Type: Journal Article
• Authors: Behrendt, Marco; de Angelis, Marco; Comerford, Liam; Zhang, Yuanjin; Beer, Michael
• Is Part Of: Mechanical systems and signal processing. Volume 172(2022)
• Publication Date: 2022-06-01
• Subjects: Discrete Fourier transform -- Complex intervals -- Dependency tracking -- Interval arithmetic -- Power spectral density estimation -- Uncertainty quantification; Structural dynamics -- Periodicals; Vibration -- Periodicals; Constructions -- Dynamique -- Périodiques; Vibration -- Périodiques; Structural dynamics; Vibration; Periodicals; 621
• DOI: 10.1016/j.ymssp.2022.108920
• Languages: English
• ISSNs: 0888-3270

Abstract

Abstract: 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. Highlights: 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.