A new representation of power spectral density and correlation function by means of fractional spectral moments

https://doi.org/10.1016/j.probengmech.2010.04.003Get rights and content

Abstract

In this paper, a new perspective for the representation of both the power spectral density and the correlation function by a unique class of function is introduced. We define the moments of order γ (γC) of the one sided power spectral density and we call them Fractional Spectral Moments (FSM). These complex quantities remain finite also in the case in which the ordinary spectral moments diverge, and are able to represent the whole Power Spectral Density and the corresponding correlation function.

Introduction

The Spectral Moments (SMs) introduced by Vanmarcke [1] are the moments of the one-sided Power Spectral Density (PSD) function. These quantities have been introduced because of their relevance in many cases of engineering interest such as the prediction of the first excursion failure, fatigue failure, and so on. In some previous works [2], [3], [4] it has been shown that in time domains the SMs are covariances of a complex process constructed in such a way that its real part is the given process, while the imaginary part is its Hilbert transform. Such a complex process is usually referred to in literature as analytic process. In particular the zero-th and second order SMs are the common variances of the process and of its derivative, respectively, while the first SM is the cross covariance of the analytic process and its first derivative. Higher order SMs have been also introduced for statistical distribution of peaks and ranges [5]. Other relevant information on the SM may be found in [6], [7], [8], [9], [10], [11], [12], [13].

In any cases, we can say that spectral moments up to the fourth order give some information on the process but higher order spectral moments are not in general useful. Moreover, higher SMs might be divergent quantities [3] and, in any case the knowledge of the complete set of SMs cannot be used to restore the whole PSD.

In this paper, we introduce a different perspective on the spectral moments problem, that allows us to formulate a representation formula for both the entire PSD and the correlation function, in all the ranges 0<ω< and 0<τ<, respectively. To this aim we define the moments of order γ(γC) of the one sided PSD and we call them as Fractional Spectral Moments (FSM). These quantities remain finite quantities also in the case in which the ordinary SMs diverge, by simply selecting a proper value of the real part of γ.

Then, a new generalized Taylor expansion in integral form in terms of Riesz fractional derivatives and integrals, involving FSMs, is introduced. By using such a Taylor integral form both the PSD and the correlation function are expressed in terms of the FSM function. Discretization of such a Taylor integral form give quite satisfactory results also in some pathological cases like the exponential type correlation function and Pierson–Moskowitz spectrum.

Section snippets

Prebelliminary concepts and definitions

In this section, some well known concepts on analytic processes and SMs are briefly outlined, for clarity’s sake as well as to introduce appropriate symbols.

Fractional Spectral Moments (FSM)

Even though the first few SMs give information on the distribution of the one-sided PSD or on the bandwidth, on the envelopes etc., they are not able to reconstruct the overall one-sided PSD. In this section we will define a new form of SM that are dubbed Fractional Spectral Moments (FSMs) defined as λXγ=0SX(ω)ωγdω;γC.

Since λXγ is a function of γ it will be denoted hereinafter as ΛX(γ). Obviously, for γ=0,1,2, the FSM defined in Eq. (7) coincide with the SM already defined in Eq. (5).

We now

Representation of PSD and correlation function by FSM

As previously described in Section 2, the classical spectral moments do not give back the PSD nor the correlation function, in general, and in many cases they might also diverge. We now will demonstrate that the fractional spectral moments previously introduced are very useful quantities to restore both the PSD and the correlation function in the whole domain. In order to achieve this result, we first recall that the Riesz fractional integral is related to the Riemann–Liouville fractional

Numerical applications

In order to show the consistency of the approximations previously discussed, two pathological examples are now shown in the next two paragraphs. In the first one, we consider the Pierson–Moskowitz spectrum for ocean wave elevation while in the second one the case RY(τ)=σX2exp(α|τ|) is examined. The former is selected because in zero the PSD is very flat, causing particular difficulties. Readers are referred to the paper of Spanos [22] for a deeper insight on the representation of such a

Conclusions

The spectral moments introduced by Vanmarcke are the moments of the one sided PSD. Even though these quantities are useful, for example in evaluating the bandwidth and the peak response, they do not give the PSD nor the correlation function. In this paper, a new concept is introduced giving a new prospective on the probabilistic characterization of stochastic processes. In particular, it has been shown that the Fractional Spectral Moment function is able to represent both the PSD and the

Acknowledgements

The authors wish to thank Pol D. Spanos for stimulating observations on the paper. The unknown referees are also acknowledged for their interesting scientific comments.

References (25)

  • M. Di Paola et al.

    Spectral moments and pre-envelope covariances of nonseparable processes

    Journal of Applied Mechanics, Transactions ASME

    (1990)
  • G. Petrucci et al.

    On the characterization of dynamic properties of random processes by spectral parameters

    Journal of Applied Mechanics

    (2000)
  • Cited by (0)

    View full text