Generation of Non-gaussian Wide-sense Stationary Random Processes with Desired Psds and Pdfs

This paper describes a new method to generate discrete signals with arbitrary power spectral density (PSD) and first order probability density function (PDF) without any limitation on PDFs and PSDs. The first approximation has been achieved by using a nonlinear transform function. At the second stage the desired PDF was approximated by a number of symmetric PDFs with defined variance. Each one provides a part of energy from total signal with different ratios of remained desired PSD. These symmetric PDFs defined by sinusoidal components with random amplitude, frequency and phase variables. Both analytic results and examples areincluded. The proposed scheme has been proved to be useful in simulations involving non-Gaussian processes with specific PSDs and PDFs.


Introduction
Generating and representing of a Gaussian variable with desired PSD is an easy job.A usual method is to use simple linear time invariant filter with a Gaussian random variable as an input.The problem of generating correlated non-Gaussian random series has been of great interest during the last three decades in connection with the simulation of processes in communication systems such as radar, sonar and speech.Furthermore, it can be applicable in recent military approaches.Groups of stealth aircrafts generate random noise with K distribution to be hidden under clutter noise from active radars.The generation of non-Gaussian signals has given rise to a renewed interest in the mechanical vibration fields and defense industries for some reasons.The realization of many surface transportation and wave environments are non-Gaussian, also the development of shaker control systems that can replicate long time histories are non-Gaussian.Many original current shaker control systems, for generating random vibration tests generate only Gaussian random noise.However, waveform replication techniques now allow the reproduction of any waveforms whose characteristics are within the bounds of a shaker [1,2].Most of the approaches suggested involve the non-linear transformation of colored Gaussian random series or use linear transformation of non-Gaussian no correlated time series.All these methods require a complicated optimization procedure to achieve the de-sired properties [3,4].The zero memory non-linear ZMNL in [3] is chosen so that the desired distribution is exactly realized and the digital filter is designed so that the specific auto-covariance is closely approximated.In this method, the autocorrelation of sequence at the output of the filter may be negative hence the method can fail because in this case we should calculate several hundreds or even thousands of unknowns, on the other hand ergodicity has not been mentioned in this paper.In addition, the method is too complex to be implemented but it can be helpful in this paper for first approximation with some promotion.In [5], a correlated Gaussian random process is multiplied by a modulating sequence.The main problem is in determining the modulating sequence, which can be difficult if not impossible.Furthermore, the generated sequence is not ergodic.In [6], a nonlinearity is used and due to the nonlinearities in these schemes, it seems that the results cannot be extended to the general cases such as multidimensional and multichannel random processes.The method suggested in [7], is based on [3] and decreases some problems and produces wide-sense stationary and ergodic sequences asymptotically or as the number of generated data goes to infinity.In [8] proposed a method that decouples the problem into two separate ones for some signal characteristics such as even and infinitely divisible first-order PDF, but it can generate every continuous PSDs.The method can be limited because the Hankel transform of every PDFs cannot be calculated easily.The restriction of an even PDF limits us to generate vast major PDFs such as Rayleigh, Nakagami, flicker and square Gaussian noises.The assumption of an infinitely divisible PDF may be restrictive too.However, the strategy mentioned in [8] is a constructive tool to be used under some modifications.In this paper, we divide desired PDF with a number of symmetric and infinitely divisible ones that these PDFs depend on the first estimation.In contradiction to [8] there is not any restriction on the PDF.This paper consists of five section, in section 2 a basic function which has been derived from [8] is considered and extended to more general form and the nonlinear transform function is described for closed form first approximation.In section 3, with some assumption we prove that this function has ergodicity properties in both means and variances and the way to compute a nonlinear transform function has been described.In section 4, a complex signal with asymmetric PDF and desired PSD are generated.The conclusion has been given in section 5.

Definition and Properties
We propose a new method to generate arbitrary continuous PSDs and first-order PDFs.In this method, one can generate a discrete signal by dividing the desired PDF into a number of symmetric PDFs.In general form we define   where and show real and imaginary parts of respectively that are independent from each other.In (1) F F L  and are random variables with the following definitions In all sections of this article, we assume , x L and , y L each have equal distributions.By considering the above-mentioned relations the following are shown in appendix The mean of every PDFs is a deterministic value and it forces the PSD to have a  function in zero.Therefore, it can be subtracted from The random process is the wide sense stationary with the following autocorrelation function where , Actually, the power spectral density is equal to Now we can use the nonlinear transform function in order to achieve the first approximation of the desired spectrum and the first order statistical probability density function.This spectrum can be converged exactly on desired one by using where   S P f is the desired power spectrum.

3) The PDF of i
A and can be written as where are the two characteristic functions for symmetric PDFs with proper variance and zero mean, for a Gaussian and we can solve the above relation by using Hankel transformation of similarly, we have In order to use the nonlinear transform function, consider the following relation between two random variables of x and t It is understood that in Equations ( 16) and ( 17), , t t stands for   and , the     2 y g t  are the density functions which are dependent on x and variables.In this paper all the relations and calculations on , variables can be modified to , variables.

y L t L t
If is not monotonically increasing, the function must be broken into parts, and each part is handled separately.For the purposes of this article x g t  will be restricted to its derivation function.Therefore, the sign of d d x t can be negative or positive.The conversion of a realization of to t x using a nonlinear function will always produce a PSD for x that is different from the PSD of .The effect is to add harmonics to all the Fourier components of .This will make the transformed data appear "whiter" than the original data [9,10].However, most of the spectral information is contained in the zero crossings which are preserved, and if the nonlinear-ity is not too great, the spectral change is usually accept-able and the PSD of t t X will be near the PSD of .The practical utility of these results is that the PSD may be approximated from (9) by using nonlinear transform function on n and is converged to exact form by choosing the PDF of frequencies

Ergodicity Properties, Nonlinear Function Synthesizing
The main problem in applications of stochastic processes is the estimation of various statistical parameters in terms of real and imaginary parts.It is clear that the basic Equation ( 1) contains two mathematical terms, the first one is  which its ergodicity of mean and variance has been proven in [1].The second term is , where g t under the limitation of Equation (10).The total power spectrum can be divided as follows are the power spectrum of real and imaginary parts respectively.Therefore, we can assume Both strategies on imaginary and real parts are like each other.Equation (18) can be simplified as and Now we propose a computational method to generate a proper function of   where we select µ as large as possible, considering small positive value for t  , we have Where  is a large number which µ   , Figure 1 shows nonlinear function synthesizing to generate uniform PDF in [-0.5, 0.5] interval.It is clear that generated function guarantees the Equations ( 20) and (21). where

Some Examples of Generating
Asymmetric and PDF with Desired PSD  

One Dimensional Real Signal
Therefore, the frequency PDF can be calculated as the following relation Consider the representation of a WSS Rayleigh random process with the variance of , mean given by , and a given continuous PSD.To represent and generate this process we need determining a nonlinear transform function to convert Gaussian distribution of to Rayleigh PDF, the PDFs of the amplitude and frequency can be calculated after the first approximation of desired PSD. Figure 2 shows a nonlinear transform function to convert Gaussian with shaped spectrum to Rayleigh one.Figure 3 depicts the Rayleigh PDF, the blue curve in Figure 4 shows the first approximated power spectral density and the red dashed line shows desired spectrum.We define the following error function The average of 10,000 realizations of a periodogram along with the true PSD (indicated by the black curve) is shown in Figure 4.It is interesting to note that the generated PDF has less accuracy than what is calculated before adding sinusoidal term and it is depicted in Figure 5. First approximation play an important role in having a more accurate results.bitrary power spectral density, where x and are independent from each other and represent the real and complex parts of respectively.Both y S x and have desired PDF and every one generates a part of total power spectrum which depends on its variance.Therefore , where is the PDF of a complex signal and ,

 
represent the PDF of real and imaginary parts of the signal.Figure 6 shows a joint PDF of p y x and .If we consider y Section 3 we can achieve to proper nonlinear transform function.Figure 7 shows generated PDF and Figure 8 depicts imaginary and real parts of this discrete signal.
The power spectrum can be converged to its correct response by considering proper PDF for both F and F  values as follows Generated and desired power spectral density also has been shown in Figure 9. where

Conclusion
A new method has been presented which can generate a complex WSS random process with a desired PSD and a given first-order PDF.One of the advantages of this method is that we have no limitation on PDF to be infi-After applying the algorithm that was mentioned in    nitely divisible or symmetric.Also with simple computational method we can calculate proper transform function to generate desired PDF from Gaussian PDF with shaped spectrum.Furthermore, this computational method guarantees generated spectrum to be under the wanted spectrum for the first approximation.The error function is reduced to zero by the means of sinusoidal component with proper frequency and amplitude PDF, but with small deviation from exact desired PDF.In other words, there is some tradeoff between exact PDF and PSD.The ergodicity in the mean and variance has been proven under certain conditions.We have avoided bottleneck calculation of Hankel transform function that [8] faces directly to it.Because of negligible deviation from exact PDF it can be used in any practical system.

Appendix
the above relation can be simplified as A F i i i

Figure 2 .
Figure 2. Nonlinear transform function converting non zero mean Gaussian PDF to Rayleigh random process.

Figure 3 .
Figure 3. Generated Rayleigh random process after applying g(t) on Gaussian PDF.

Figure 5 .
Figure 5. First order probability density function of Z[n], red line depicts Rayleigh PDF.

Figure 6 .
Figure 6.Joint PDF of x and y.

Figure 7 .
Figure 7. Generated joint PDF of Re{Z} and Im{Z}.

Figure 8 .
Figure 8. Discrete imaginary and real parts of generated signal.
t has an ergodic properties and