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

Authors: Guang-Dong Zhou; You-Liang Ding; Ai-Qun Li
Journal: Mathematical Problems in Engineering
 (2015)

---

## 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.


---

## Body

## 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 timet as the Fourier transform of an instantaneous correlation function R(τ,t), which corresponds to the standard correlation function R(τ) with time lag τ centered at time t [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 t since the variable τ must be integrated over an infinite range to compute the Fourier transform of the instantaneous correlation function R(τ,t). 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 ψa,b(t) originate from the “mother” wavelet ψ(t) by
(1)ψa,bt=1aψt-baa,b∈R,a>0,
where a and b represent the scale parameter and shift parameter, respectively.The continuous wavelet transform (CWT) of a finite energy stochastic processf(t) is defined as
(2)WTfa,b=ft,ψa,bt=1a∫-∞+∞ftψ*t-badt,
where WTf(a,b) 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 parametera and the shift parameter b are discretized as
(3)ar=a0rr=1,2,…mbk=k-1b0k=1,2,…n,
where a0 and b0 are constants, m represents the total number of discrete scales, and n denotes the total number of time intervals retained in the discrete process f(t).


### 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 processg(t), its spectral decomposition equation can be expressed as
(4)gt=∫-∞∞eiωtdZω,
where i represents the imaginary unit and Z(ω) 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 f(t) can be represented by the following Riemann-Stieltjes (R-S) integration:
(5)ft=∫-∞+∞Hω,teiωtdZω,
where H(ω,t) denotes the slowly varying time- and frequency-dependent modulating function, namely, non-uniformly modulating function.Then, the EPSD of the nonstationary stochastic processf(t) is defined as
(6)Sffω,t=Hω,t2Sggω,
where Sff(ω,t) is the EPSD of the process f(t) and Sgg(ω) is the time-independent power spectral density of g(t). It can be seen from (6) that, by using the modulating function, the time-independent power spectral density Sgg(ω) 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 functionh(t), and (5) and (6) are simplified to
(7)ft=∫-∞+∞Hω,teiωtdZω=htgt,Sffω,t=Hω,t2Sggω=ht2Sggω.


### 2.3. Wavelet-Based EPSD Estimation

According to the definition of the Fourier transform, the Fourier transform off(t) is
(8)Fω=∫-∞∞fte-iωtdt,
where F(ω) represents the Fourier transform of f(t). Substituting (5) into (8) results in
(9)Fω=∫-∞∞∫-∞∞Hω,teiωtdZωe-iωtdt.
Employing the definition of CWT in frequency domain, the CWT of f(t) is
(10)WTf(a,b)=a2π∫-∞∞∫-∞∞∫-∞∞Hω,teiωtdZωe-iωtdtllllllllllllllllllllllll×Ψ*(aω)eiωbdω.
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 a have one-to-one corresponding relation with the frequency ω in the wavelet transform. Thus, (5) and (6) can be approximately rewritten in the form
(11)ft≈∫-∞∞Ha,teiatdZa,Sffa,t≈Ha,t2Sgga.
Then, the following equation can be obtained as follows:(12)WTf(a,b)≈a2π∫-∞+∞∫-∞+∞∫-∞+∞H(a,t)eiatdZ(a)e-iωtdtlllllllllllllllllllllll×Ψ*(aω)eiωbdω.It is worth noting that (12) becomes more accurate if the change rate of the modulation function H(a,t) is smaller than that of the wavelet functions in frequency domain. At a given time instant b=bk, the wavelet transform of the process f(t) is
(13)WTf(a,bk)≈a2π∫-∞+∞∫-∞+∞∫-∞+∞H(a,t)eiatdZ(a)e-iωtdtlllllllllllllllllllllll×Ψ*(aω)eiωbkdω=∫-∞∞∫-∞∞H(a,t)a2π∫-∞∞Ψ*(aω)e-iω(t-bk)dωlllllllllllllllllllllllll×eiateiatdZaa2π∫-∞∞e-iatdt.
It can be found that the item (a/2π)∫-∞∞Ψ*(aω)e-iω(t-bk)dω in (13) is the inverse Fourier transform of the wavelet function at time instant bk. That is,
(14)a2π∫-∞+∞Ψ*(aω)e-iω(t-bk)dω=1aψ-t-bka.
Substituting (14) into (13), the following equation can be derived:
(15)WTf(a,bk)=∫-∞∞∫-∞∞H(a,t)1aψ-t-bkaeiateiatdZalllllllllllllllll×e-iatdt.
Define a new nonstationary stochastic process p(t), whose R-S integration is
(16)pt=∫-∞+∞Da,teiatdZa,
where D(a,t) represents the non-uniformly modulating function of p(t), which is
(17)Da,t=Ha,t1aψ-t-bkaeiat.
Then, (15) can be simplified as
(18)WTfa,bk=∫-∞+∞pte-iatdt.
Therefore, the wavelet coefficients, at a given instant bk, WTf(a,bk), can be interpreted as the Fourier transform of p(t), which is a new nonstationary stochastic process defined by the non-uniformly modulating function D(a,t). The D(a,t) is determined by (17), which indicates that D(a,t) is equal to the wavelet function in time domain multiplied by the non-uniformly modulating function H(a,t) of the original nonstationary stochastic process f(t). According to (6), the EPSD of the nonstationary stochastic process p(t) is
(19)Sppa,t=D(a,t)2Sgg(a)=H(a,t)21aψ-t-bka2Sgg(a)=1aψ-t-bka2Sffa.
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:
(20)WTf(a,bk)2=∫-∞+∞Spp(a,t)dt=∫-∞+∞1aψbk-ta2Sffadt.
Considering several samples of the process f(t), 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
(21)EWTfa,bk2=∫-∞+∞1aψbk-ta2Sff(a,t)dt.
The general solution to (21) takes the following form:
(22)Sffa,t=∑k=1nγka1aψbj-ta2j=1,2,…n,
where γk(a) is the frequency-dependent coefficient. Substituting (22) into (21), the following equation can be made:
(23)EWTfa,bk2=∑k=1nγk(a)∫-∞+∞1aψbk-ta21aψbj-ta2dt=∑k=1nγkaδk,j.To sum up, the EPSD of the processf(t) 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 (24) 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: (1) implementing wavelet transform to obtain wavelet coefficients with selected scales and different time instants; (2) establishing the frequency-dependent coefficients by solving (23) for all scales; (3) assembling the frequency-dependent coefficients and time-dependent wavelet functions to calculate the EPSD of the process f(t). 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.


## 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 ψa,b(t) originate from the “mother” wavelet ψ(t) by
(1)ψa,bt=1aψt-baa,b∈R,a>0,
where a and b represent the scale parameter and shift parameter, respectively.The continuous wavelet transform (CWT) of a finite energy stochastic processf(t) is defined as
(2)WTfa,b=ft,ψa,bt=1a∫-∞+∞ftψ*t-badt,
where WTf(a,b) 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 parametera and the shift parameter b are discretized as
(3)ar=a0rr=1,2,…mbk=k-1b0k=1,2,…n,
where a0 and b0 are constants, m represents the total number of discrete scales, and n denotes the total number of time intervals retained in the discrete process f(t).


## 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 processg(t), its spectral decomposition equation can be expressed as
(4)gt=∫-∞∞eiωtdZω,
where i represents the imaginary unit and Z(ω) 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 f(t) can be represented by the following Riemann-Stieltjes (R-S) integration:
(5)ft=∫-∞+∞Hω,teiωtdZω,
where H(ω,t) denotes the slowly varying time- and frequency-dependent modulating function, namely, non-uniformly modulating function.Then, the EPSD of the nonstationary stochastic processf(t) is defined as
(6)Sffω,t=Hω,t2Sggω,
where Sff(ω,t) is the EPSD of the process f(t) and Sgg(ω) is the time-independent power spectral density of g(t). It can be seen from (6) that, by using the modulating function, the time-independent power spectral density Sgg(ω) 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 functionh(t), and (5) and (6) are simplified to
(7)ft=∫-∞+∞Hω,teiωtdZω=htgt,Sffω,t=Hω,t2Sggω=ht2Sggω.


## 2.3. Wavelet-Based EPSD Estimation

According to the definition of the Fourier transform, the Fourier transform off(t) is
(8)Fω=∫-∞∞fte-iωtdt,
where F(ω) represents the Fourier transform of f(t). Substituting (5) into (8) results in
(9)Fω=∫-∞∞∫-∞∞Hω,teiωtdZωe-iωtdt.
Employing the definition of CWT in frequency domain, the CWT of f(t) is
(10)WTf(a,b)=a2π∫-∞∞∫-∞∞∫-∞∞Hω,teiωtdZωe-iωtdtllllllllllllllllllllllll×Ψ*(aω)eiωbdω.
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 a have one-to-one corresponding relation with the frequency ω in the wavelet transform. Thus, (5) and (6) can be approximately rewritten in the form
(11)ft≈∫-∞∞Ha,teiatdZa,Sffa,t≈Ha,t2Sgga.
Then, the following equation can be obtained as follows:(12)WTf(a,b)≈a2π∫-∞+∞∫-∞+∞∫-∞+∞H(a,t)eiatdZ(a)e-iωtdtlllllllllllllllllllllll×Ψ*(aω)eiωbdω.It is worth noting that (12) becomes more accurate if the change rate of the modulation function H(a,t) is smaller than that of the wavelet functions in frequency domain. At a given time instant b=bk, the wavelet transform of the process f(t) is
(13)WTf(a,bk)≈a2π∫-∞+∞∫-∞+∞∫-∞+∞H(a,t)eiatdZ(a)e-iωtdtlllllllllllllllllllllll×Ψ*(aω)eiωbkdω=∫-∞∞∫-∞∞H(a,t)a2π∫-∞∞Ψ*(aω)e-iω(t-bk)dωlllllllllllllllllllllllll×eiateiatdZaa2π∫-∞∞e-iatdt.
It can be found that the item (a/2π)∫-∞∞Ψ*(aω)e-iω(t-bk)dω in (13) is the inverse Fourier transform of the wavelet function at time instant bk. That is,
(14)a2π∫-∞+∞Ψ*(aω)e-iω(t-bk)dω=1aψ-t-bka.
Substituting (14) into (13), the following equation can be derived:
(15)WTf(a,bk)=∫-∞∞∫-∞∞H(a,t)1aψ-t-bkaeiateiatdZalllllllllllllllll×e-iatdt.
Define a new nonstationary stochastic process p(t), whose R-S integration is
(16)pt=∫-∞+∞Da,teiatdZa,
where D(a,t) represents the non-uniformly modulating function of p(t), which is
(17)Da,t=Ha,t1aψ-t-bkaeiat.
Then, (15) can be simplified as
(18)WTfa,bk=∫-∞+∞pte-iatdt.
Therefore, the wavelet coefficients, at a given instant bk, WTf(a,bk), can be interpreted as the Fourier transform of p(t), which is a new nonstationary stochastic process defined by the non-uniformly modulating function D(a,t). The D(a,t) is determined by (17), which indicates that D(a,t) is equal to the wavelet function in time domain multiplied by the non-uniformly modulating function H(a,t) of the original nonstationary stochastic process f(t). According to (6), the EPSD of the nonstationary stochastic process p(t) is
(19)Sppa,t=D(a,t)2Sgg(a)=H(a,t)21aψ-t-bka2Sgg(a)=1aψ-t-bka2Sffa.
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:
(20)WTf(a,bk)2=∫-∞+∞Spp(a,t)dt=∫-∞+∞1aψbk-ta2Sffadt.
Considering several samples of the process f(t), 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
(21)EWTfa,bk2=∫-∞+∞1aψbk-ta2Sff(a,t)dt.
The general solution to (21) takes the following form:
(22)Sffa,t=∑k=1nγka1aψbj-ta2j=1,2,…n,
where γk(a) is the frequency-dependent coefficient. Substituting (22) into (21), the following equation can be made:
(23)EWTfa,bk2=∑k=1nγk(a)∫-∞+∞1aψbk-ta21aψbj-ta2dt=∑k=1nγkaδk,j.To sum up, the EPSD of the processf(t) 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 (24) 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: (1) implementing wavelet transform to obtain wavelet coefficients with selected scales and different time instants; (2) establishing the frequency-dependent coefficients by solving (23) for all scales; (3) assembling the frequency-dependent coefficients and time-dependent wavelet functions to calculate the EPSD of the process f(t). 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
(24)ψt=e-t2/2eiωt,Ψω=2πe-(ω-ω0)2/2,
where ω0 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 g(t). The process g(t) 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:
(25)ht=exp⁡-t-α2β,
where α and β are problem specific parameters. In this example, the selected α and β are 200 and 18000, respectively; then the sketch of h(t) is shown in Figure 1.Figure 1
Uniformly modulating functionh(t).The simulated wind duration and the sampling frequency are set as 450 s and 4 Hz, respectively. 500 samples of the processf(t) are synthesized and the average values of wavelet coefficients are calculated. Then, the EPSD of simulated typhoon processes is estimated by (22). The estimated EPSD is compared with the theoretical solution to evaluate the reliability of the proposed method. In order to clearly display the results, the three-dimensional EPSD is sliced at different time instants and different frequency points, as shown in Figure 2, which are named as the time slice and the frequency slice in this paper, respectively. In the figure, the unit of frequency is converted to Hz. The comparison between the estimated results and the theoretical solution indicates that the estimated EPSD has a satisfactory agreement with the theoretical solution in both time and frequency domains although the estimated EPSD may deviate slightly from the theoretical solution at some time instants and frequency points. The effectiveness of the proposed wavelet-based EPSD estimation method is validated. The EPSD extends the two-dimensional power spectral density with energy and frequency to three-dimensional evolutionary power spectral density with energy, frequency, and time. One can investigate the energy distribution of a signal at any time instant and frequency point, which is significant for accurately predicting structural dynamic responses under typhoon processes.Comparison of estimated EPSD of uniformly modulated typhoon process with theoretical solution.
(a)
t
=
125 s
(b)
t
=
250 s
(c)
ω
=
0.05 Hz
(d)
ω
=
0.46 HzUsing the marginal features, the mean-square values of time-independent power spectral density of the nonstationary typhoon processf(t) and the instantaneous mean-square values of the nonstationary typhoon process f(t) are obtained. Figure 3 displays the comparison between the estimated marginal features and the exact solution. An excellent agreement can be found between the two results in the figure, which again demonstrates that the proposed method is applicable in estimating the EPSD of nonstationary typhoon processes modulated by the uniformly modulating function.Comparison of marginal features of estimated EPSD with theoretical solution.
(a)
Marginal feature in time domain
(b)
Marginal feature in frequency domain


### 3.3. Non-Uniformly Modulated Nonstationary Typhoon Process

Suppose the non-uniformly modulated nonstationary typhoon processf(t) is formulated as
(26)ft=h1tg1t+h2tg2t,
where h1(t) and h2(t) are two uniformly modulating functions and g1(t) and g2(t) are two stationary typhoon processes described by the Simiu spectra. h1(t) and h2(t) are formulated by
(27)h1t=exp⁡-t-α12β1,h2t=exp⁡-t-α22β2,
where α1, α2, β1, and β2 are problem specific parameters. Here, α1, α2, β1, and β2 are set as 100, 300, 9000, and 9000, respectively. Figure 4 displays the values of h1(t) and h2(t) changing with time. Then, the EPSD of the non-uniformly modulated nonstationary wind typhoon process f(t) is
(28)Sffθ,t=h1t2Sg1g1θ+h2t2Sg2g2θ.Figure 4
Uniformly modulating functionsh1(t) and h2(t).Using the same approach in Section3.2, the non-uniformly modulated nonstationary typhoon process f(t) is synthesized and the EPSD of f(t) is estimated. The comparison of time slices and frequency slices between the estimated results and the exact solutions is illustrated in Figure 5. In most cases, as displayed in Figures 5(a), 5(b), and 5(d), the estimated results agree well with the exact solutions. Since it is impossible to narrow the window of Morlet function to a frequency point, there are some derivations in the time slices, as shown in Figure 5(c). But the proposed method is successful to capture the peak values of the EPSD, which may induce intensive structural responses; therefore, the estimated EPSD can be used in practical application and the big errors in the valley are acceptable in structural dynamic analysis. The proposed wavelet-based method is effective to estimate nonstationary typhoon processes modulated by the non-uniformly modulating function.Comparison of estimated EPSD of non-uniformly modulated typhoon process with theoretical solution.
(a)
t
=
10 s
(b)
t
=
250 s
(c)
ω
=
0.06 Hz
(d)
ω
=
0.24 HzThe marginal features of the EPSD of non-uniformly modulated nonstationary typhoon processes are also calculated. The comparison results between the estimated values and the theoretical solutions are shown in Figure6. The reasonable agreement can also be found, which validates the capability of the proposed method in capturing the time-frequency properties of nonstationary typhoon processes.Comparison of marginal features of estimated EPSD with theoretical solution.
(a)
Marginal feature in time domain
(b)
Marginal feature in frequency domain


## 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
(24)ψt=e-t2/2eiωt,Ψω=2πe-(ω-ω0)2/2,
where ω0 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 g(t). The process g(t) 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:
(25)ht=exp⁡-t-α2β,
where α and β are problem specific parameters. In this example, the selected α and β are 200 and 18000, respectively; then the sketch of h(t) is shown in Figure 1.Figure 1
Uniformly modulating functionh(t).The simulated wind duration and the sampling frequency are set as 450 s and 4 Hz, respectively. 500 samples of the processf(t) are synthesized and the average values of wavelet coefficients are calculated. Then, the EPSD of simulated typhoon processes is estimated by (22). The estimated EPSD is compared with the theoretical solution to evaluate the reliability of the proposed method. In order to clearly display the results, the three-dimensional EPSD is sliced at different time instants and different frequency points, as shown in Figure 2, which are named as the time slice and the frequency slice in this paper, respectively. In the figure, the unit of frequency is converted to Hz. The comparison between the estimated results and the theoretical solution indicates that the estimated EPSD has a satisfactory agreement with the theoretical solution in both time and frequency domains although the estimated EPSD may deviate slightly from the theoretical solution at some time instants and frequency points. The effectiveness of the proposed wavelet-based EPSD estimation method is validated. The EPSD extends the two-dimensional power spectral density with energy and frequency to three-dimensional evolutionary power spectral density with energy, frequency, and time. One can investigate the energy distribution of a signal at any time instant and frequency point, which is significant for accurately predicting structural dynamic responses under typhoon processes.Comparison of estimated EPSD of uniformly modulated typhoon process with theoretical solution.
(a)
t
=
125 s
(b)
t
=
250 s
(c)
ω
=
0.05 Hz
(d)
ω
=
0.46 HzUsing the marginal features, the mean-square values of time-independent power spectral density of the nonstationary typhoon processf(t) and the instantaneous mean-square values of the nonstationary typhoon process f(t) are obtained. Figure 3 displays the comparison between the estimated marginal features and the exact solution. An excellent agreement can be found between the two results in the figure, which again demonstrates that the proposed method is applicable in estimating the EPSD of nonstationary typhoon processes modulated by the uniformly modulating function.Comparison of marginal features of estimated EPSD with theoretical solution.
(a)
Marginal feature in time domain
(b)
Marginal feature in frequency domain


## 3.3. Non-Uniformly Modulated Nonstationary Typhoon Process

Suppose the non-uniformly modulated nonstationary typhoon processf(t) is formulated as
(26)ft=h1tg1t+h2tg2t,
where h1(t) and h2(t) are two uniformly modulating functions and g1(t) and g2(t) are two stationary typhoon processes described by the Simiu spectra. h1(t) and h2(t) are formulated by
(27)h1t=exp⁡-t-α12β1,h2t=exp⁡-t-α22β2,
where α1, α2, β1, and β2 are problem specific parameters. Here, α1, α2, β1, and β2 are set as 100, 300, 9000, and 9000, respectively. Figure 4 displays the values of h1(t) and h2(t) changing with time. Then, the EPSD of the non-uniformly modulated nonstationary wind typhoon process f(t) is
(28)Sffθ,t=h1t2Sg1g1θ+h2t2Sg2g2θ.Figure 4
Uniformly modulating functionsh1(t) and h2(t).Using the same approach in Section3.2, the non-uniformly modulated nonstationary typhoon process f(t) is synthesized and the EPSD of f(t) is estimated. The comparison of time slices and frequency slices between the estimated results and the exact solutions is illustrated in Figure 5. In most cases, as displayed in Figures 5(a), 5(b), and 5(d), the estimated results agree well with the exact solutions. Since it is impossible to narrow the window of Morlet function to a frequency point, there are some derivations in the time slices, as shown in Figure 5(c). But the proposed method is successful to capture the peak values of the EPSD, which may induce intensive structural responses; therefore, the estimated EPSD can be used in practical application and the big errors in the valley are acceptable in structural dynamic analysis. The proposed wavelet-based method is effective to estimate nonstationary typhoon processes modulated by the non-uniformly modulating function.Comparison of estimated EPSD of non-uniformly modulated typhoon process with theoretical solution.
(a)
t
=
10 s
(b)
t
=
250 s
(c)
ω
=
0.06 Hz
(d)
ω
=
0.24 HzThe marginal features of the EPSD of non-uniformly modulated nonstationary typhoon processes are also calculated. The comparison results between the estimated values and the theoretical solutions are shown in Figure6. The reasonable agreement can also be found, which validates the capability of the proposed method in capturing the time-frequency properties of nonstationary typhoon processes.Comparison of marginal features of estimated EPSD with theoretical solution.
(a)
Marginal feature in time domain
(b)
Marginal feature in frequency domain


## 4. Conclusions

Comprehensively understanding the energy distribution of typhoon processes in both time and frequency domains is significant for reasonably predicting typhoon-induced structural responses. The nonstationary behaviors are found in a large number of field measured typhoon processes, which is different from conventional perspective. The EPSD is a powerful index to describe the time-frequency properties of nonstationary typhoon processes, which belong to the class of nonstationary oscillatory processes. A wavelet-based EPSD estimation method is proposed in this paper. Some conclusions are drawn as follows.(1)
Using the definitions of the continuous wavelet transform and the EPSD, the closed-form expressions used to compute the EPSD are rigorously deduced. The EPSD of nonstationary oscillatory processes can be represented as the weighted sum of the square moduli of wavelet functions in time domain. The weighted coefficients are frequency-dependent and determined by wavelet coefficients and the overlapping area of two wavelets at different time instants. So, the EPSD can be estimated only by wavelet functions in time domain and wavelet coefficients.(2)
The proposed wavelet-based EPSD estimation method is verified by two numerical examples, which employ the stationary typhoon processes modulated by the uniformly modulating function and the non-uniformly modulating function. The results indicate that the accuracy of the EPSD estimated by the proposed wavelet-based method is adequate for predicting structural typhoon-induced responses. The proposed wavelet-based method is applicable in the EPSD estimation of nonstationary typhoon processes and can be extended to other nonstationary oscillatory processes.The nonstationary features of signals are important but rather complex. This paper is only a preliminary research for nonstationary typhoon processes. The investigation of nonstationary typhoon processes and nonstationary wind-induced vibrations may be a promising topic in future works.



---

*Source: 870420-2015-03-03.xml*