Mathematical Problems in Engineering

Volume 2015 (2015), Article ID 870420, 10 pages

http://dx.doi.org/10.1155/2015/870420

## Wavelet-Based Methodology for Evolutionary Spectra Estimation of Nonstationary Typhoon Processes

^{1}College of Civil and Transportation Engineering, Hohai University, Nanjing 210098, China^{2}China Key Laboratory of Concrete and Prestressed Concrete Structures of Ministry of Education, Southeast University, Nanjing 210096, China

Received 14 August 2014; Accepted 9 October 2014

Academic Editor: Ting-Hua Yi

Copyright © 2015 Guang-Dong Zhou et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

Closed-form expressions are proposed to estimate the evolutionary power spectral density (EPSD) of nonstationary typhoon processes by employing the wavelet transform. Relying on the definition of the EPSD and the concept of the wavelet transform, wavelet coefficients of a nonstationary typhoon process at a certain time instant are interpreted as the Fourier transform of a new nonstationary oscillatory process, whose modulating function is equal to the modulating function of the nonstationary typhoon process multiplied by the wavelet function in time domain. Then, the EPSD of nonstationary typhoon processes is deduced in a closed form and is formulated as a weighted sum of the squared moduli of time-dependent wavelet functions. The weighted coefficients are frequency-dependent functions defined by the wavelet coefficients of the nonstationary typhoon process and the overlapping area of two shifted wavelets. Compared with the EPSD, defined by a sum of the squared moduli of the wavelets in frequency domain in literature, this paper provides an EPSD estimation method in time domain. The theoretical results are verified by uniformly modulated nonstationary typhoon processes and non-uniformly modulated nonstationary typhoon processes.

#### 1. Introduction

The wind-induced dynamic responses are important factors that must be taken into account during the design and assessment of long-span bridges, tall buildings, and other large-scale structures above ground or sea level with low stiffness and damping. In order to exactly predict structural wind-induced vibration, comprehensively understanding wind properties is the prior work. Traditional concepts treat the wind as a realization of a stationary stochastic process, whose statistical moments and frequency content keep constants during the entire duration. However, due to the energy dissipation when it attacks a structure, sudden and sporadic behaviors exhibit in many field monitoring typhoon processes. Therefore, it is more realistic to model typhoon winds as nonstationary processes, whose amplitudes, mean values, and frequency content change in time [1, 2]. Clearly, nonstationary features cannot be captured by conventional spectra estimation methods such as the Fourier transform, the Welch method, and the modern spectra estimation technology since these methods only calculate the mean composition of the power spectral density during the whole process.

In the past decades, intensive research efforts have been made to describe the time-dependent spectral characteristics of different nonstationary processes and many methods have been developed, such as the Wigner-Ville method (WVM), the short-time Fourier transform (STFT), and the S-transform (ST). The WVM defines a time-dependent spectrum at time as the Fourier transform of an instantaneous correlation function , which corresponds to the standard correlation function with time lag centered at time [3–5]. The STFT applies a moveable narrow window to establish a local spectrum for any time instant based on the fast Fourier transform [2]. The ST, which is a generalization of the STFT, employs the scalable Gaussian window to view the local frequency content of nonstationary signals [6, 7]. Unfortunately, these methods enclose some intrinsic drawbacks. The WVM cannot reflect the actual local behavior of the process at time since the variable must be integrated over an infinite range to compute the Fourier transform of the instantaneous correlation function . The STFT cannot achieve high resolutions in both time and frequency domains simultaneously due to the inherent limitation of the Fourier transform. And the ST may miss high frequency energy and it requires highly complex computation.

In practical application, such as investigating nonstationary features of the typhoon process, high resolution in both time and frequency domains is desired. To overcome the intrinsic shortcomings of existing nonstationary signal processing technologies, the wavelet transform has been developed. The wavelet transform has the capability of focusing on multiscale details of a signal by a series of moveable and scalable time-frequency windows, which are named as wavelets, and provides a promising tool to explore signal features in both time and frequency domains [8, 9]. Applications of the wavelet transform to engineering problems can be found in a large number of studies that relate to signal processing [10–12], dynamic analysis of structures [13, 14], system identification [15–17], and damage detection [18–20]. Recently, wavelet analysis is getting utilized as a method of transforming time-sequential data to data on a time-frequency plane in civil engineering. Basu and Gupta have introduced closed-form expressions of wavelet-based spectra by modified Littlewood-Paley wavelets for pursuing accurate time-dependent spectra of seismic records [21, 22]. Iyama and Kuwamura have also proposed a time-dependent spectrum based on nonorthogonal wavelets under restrictive assumptions [23]. However, in these studies, the computed time-dependent spectrum has been introduced in a somewhat heuristic manner as the local contribution to the mean-square value of the Fourier transform of the process and no rigorous treatment has been presented. In 2004, an innovated wavelet-based method was developed by Spanos and Failla to estimate the evolutionary power spectral density (EPSD) of a specific class of Priestley’s oscillatory processes, such as seismic motions and fluctuating winds [24]. In this work, the rigorous derivation has been conducted and the EPSD is represented as a sum of time-independent shape functions, which is equal to the squared moduli of the Fourier transforms of the wavelets associated with the selected scales, modulated by time-dependent coefficients. Subsequently, the method has been extended to arbitrary nonstationary processes by Failla et al. [5]

The method developed by Spanos and Failla provides a desirable accuracy to estimate the EPSD of Priestley’s oscillatory processes and is effective in a large number of wavelet functions, including those which are not orthogonal in frequency domain. However, this method has an intrinsic drawback that it is inapplicable to those wavelets whose analytical expressions in frequency domain are not existent since the EPSD is formulated by the closed-form expressions of wavelet in frequency domain. In this paper, the EPSD of nonstationary typhoon processes is deduced and constructed by closed-form representations of wavelets in time domain, which can break the restriction in the method proposed by Spanos and Failla that the analytical frequency-dependent wavelet functions are required. And then, the reliability of the proposed method is demonstrated by uniformly modulated nonstationary typhoon processes and non-uniformly modulated nonstationary typhoon processes.

#### 2. Wavelet-Based Evolutionary Spectra Estimation

##### 2.1. Brief Description of the Wavelet Transform

The wavelet transform is a linear transform, which decomposes a signal via basis functions that are simply generated by scaling and shifting the “mother” wavelet through the convolution of the signal [25]. Scaling allows the time duration of the nonzero portion of the wavelet to be adjusted according to the local frequency content of the signal. And shifting enables the wavelet to observe the local frequency content of the signal at different time instant. The basis functions originate from the “mother” wavelet by where and represent the scale parameter and shift parameter, respectively.

The continuous wavelet transform (CWT) of a finite energy stochastic process is defined as where is the wavelet coefficient, means the convolution, and the symbol represents the complex conjugate.

In numerical analysis, the discrete version of wavelet transform is generally employed. The scale parameter and the shift parameter are discretized as where and are constants, represents the total number of discrete scales, and denotes the total number of time intervals retained in the discrete process .

##### 2.2. Brief Description of the Evolutionary Power Spectral Density

Because of the complexity of nonstationary stochastic processes, defining a universal time-dependent spectrum for arbitrary nonstationary processes is impracticable. The EPSD, which is introduced by Priestley [26], is a time-dependent analogue of the spectrum of a stationary stochastic process and valid for a class of nonstationary oscillatory processes. The EPSD describes the local power-frequency distribution at each time instant, so it preserves the physical interpretation of the spectrum as a power-frequency distribution of a stationary stochastic process and is applied in many fields [27]. After extracting the time-varying mean wind speed, the remaining fluctuating wind speed of the typhoon process can be represented by a zero-mean nonstationary oscillatory process [1]. Therefore, the EPSD is introduced to describe the time-frequency features of the fluctuating component of the typhoon process.

Regarding an arbitrary stationary stochastic process , its spectral decomposition equation can be expressed as where represents the imaginary unit and denotes a complex random process with orthogonal increments.

Since the spectral decomposition equation as (4) of nonstationary processes is not existent, the nonstationary stochastic process can be represented by the following Riemann-Stieltjes (R-S) integration: where denotes the slowly varying time- and frequency-dependent modulating function, namely, non-uniformly modulating function.

Then, the EPSD of the nonstationary stochastic process is defined as where is the EPSD of the process and is the time-independent power spectral density of . It can be seen from (6) that, by using the modulating function, the time-independent power spectral density defined only in frequency domain is extended to three-dimensional space with time, frequency, and amplitude.

In engineering application, the non-uniformly modulating function is often replaced by a slowly varying uniformly modulating function , and (5) and (6) are simplified to

##### 2.3. Wavelet-Based EPSD Estimation

According to the definition of the Fourier transform, the Fourier transform of is where represents the Fourier transform of . Substituting (5) into (8) results in Employing the definition of CWT in frequency domain, the CWT of is Since the uniformly modulating functions of practical signals in civil engineering vary slowly and the wavelets in frequency domain decay quickly at the vicinity of the central frequency, it can be assumed that the scales have one-to-one corresponding relation with the frequency in the wavelet transform. Thus, (5) and (6) can be approximately rewritten in the form Then, the following equation can be obtained as follows:

It is worth noting that (12) becomes more accurate if the change rate of the modulation function is smaller than that of the wavelet functions in frequency domain. At a given time instant , the wavelet transform of the process is It can be found that the item in (13) is the inverse Fourier transform of the wavelet function at time instant . That is, Substituting (14) into (13), the following equation can be derived: Define a new nonstationary stochastic process , whose R-S integration is where represents the non-uniformly modulating function of , which is Then, (15) can be simplified as Therefore, the wavelet coefficients, at a given instant , , can be interpreted as the Fourier transform of , which is a new nonstationary stochastic process defined by the non-uniformly modulating function . The is determined by (17), which indicates that is equal to the wavelet function in time domain multiplied by the non-uniformly modulating function of the original nonstationary stochastic process . According to (6), the EPSD of the nonstationary stochastic process is Using the marginal property of the EPSD, that the integral of EPSD with time is equal to the square of the time-independent power spectral density, the following equation can be obtained: Considering several samples of the process , the EPSD of different samples may be different. The mathematical expectation, which can roughly represent the average value of a large number of samples with the same probability distribution, is employed to describe different samples belonging to a typhoon process. Setting mathematical expectation in the left-hand side of (20) derives The general solution to (21) takes the following form: where is the frequency-dependent coefficient. Substituting (22) into (21), the following equation can be made:

To sum up, the EPSD of the process can be calculated by (22). This equation can be explained as that the EPSD is represented as a sum of time-dependent shape functions modulated by frequency-dependent coefficients. The shape functions are the squared moduli of wavelet functions in time domain during the whole duration, and the coefficients relate to the squared values of wavelet transform results and two wavelet functions with different time shifts. Comparing (22) with in [24], the difference is obvious although the format of the two expressions is similar. In this paper, the frequency-dependent coefficients are determined at each frequency point and the EPSD is constructed by wavelet functions in time domain. But, in [24], the time-dependent coefficients are determined at each time instant and the EPSD is constructed by wavelet functions in frequency domain. Therefore, combining the method proposed in this paper with the equation in [24], the EPSD can be estimated by various wavelets whether their analytical expressions exist in time or frequency domain. The procedure of the EPSD estimation is as follows: implementing wavelet transform to obtain wavelet coefficients with selected scales and different time instants; establishing the frequency-dependent coefficients by solving (23) for all scales; assembling the frequency-dependent coefficients and time-dependent wavelet functions to calculate the EPSD of the process . It should be noted that all theoretical formulations are originated from the definition of the EPSD, which is introduced to describe a class of nonstationary oscillatory processes. Therefore, the results in this paper are only applicable in oscillatory processes. When employing the wavelet transform, only the definition of the wavelet transform is used and the function of wavelet is not specified. So, in theory, the method proposed in this paper is applicable to any wavelet whose function in time domain exists.

#### 3. Numerical Example

##### 3.1. Wavelet Function

From (22) and (23), it can be concluded that the accuracy of the wavelet-based EPSD relies heavily on the wavelet function. So, an appropriate wavelet function should be selected before carrying out EPSD estimation. Since large-scale structures are susceptible to the energy distribution of excitation in both time and frequency domains, the high time resolution and frequency resolution of wavelet functions are required. However, according to the Heisenberg uncertainty principle, the time resolution and the frequency resolution are inconsistent. It is impossible that a wavelet function has highest resolution in both time and frequency domains simultaneously, so only the compromise solution is available. In this paper, the Morlet wavelet, which has acceptable resolution in both time and frequency domains, is adopted to estimate the EPSD of the typhoon process. Its expressions in time domain and frequency domain are where represents a parameter that involves the center frequency of the time-frequency window.

##### 3.2. Uniformly Modulated Nonstationary Typhoon Process

Early in 1982, Lin and Yang have suggested that the nonstationary turbulence wind can be described by a stationary stochastic process multiplied with a determined modulating function [28]. In this paper, the concept is employed to simulate the nonstationary wind. The Simiu spectra are introduced as the target power spectra of the stationary stochastic wind . The process can be synthesized by the classical spectral representation method [29, 30]. Since there is no widely accepted uniformly modulating function describing the typhoon process, at present, the modulating function in earthquake engineering is referred to and the modified formula is supposed to be as follows: where and are problem specific parameters. In this example, the selected and are 200 and 18000, respectively; then the sketch of is shown in Figure 1.