- About this Journal ·
- Abstracting and Indexing ·
- Aims and Scope ·
- Annual Issues ·
- Article Processing Charges ·
- Articles in Press ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents
International Journal of Antennas and Propagation
Volume 2013 (2013), Article ID 103461, 8 pages
A TDL Based Non-WSSUS Vehicle-to-Vehicle Channel Model
1State Key Laboratory of Rail Traffic Control and Safety, Beijing Jiaotong University, No. 3 Shang Yuan Cun, Haidian District, Beijing 100044, China
2School of Electronics Engineering and Computer Science, Peking University, No. 5 Yiheyuan Road, Haidian District, Beijing 100871, China
3School of Electronic and Information Engineering, Beijing Jiaotong University, No. 3 Shang Yuan Cun, Haidian District, Beijing 100044, China
Received 9 July 2013; Accepted 18 October 2013
Academic Editor: César Briso Rodríguez
Copyright © 2013 Yan Li 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.
This paper proposes a non-wide-sense-stationary-uncorrelated scattering (WSSUS) channel model for vehicle-to-vehicle (V2V) communication systems. The proposed channel model is based on the tapped-delay line (TDL) structure and considers the correlation between taps both in amplitude and phase. Using the relationship between the correlation coefficients of complex Gaussian, Weibull, and Uniform random variables (RVs), the amplitude and the phase of taps with different delays are modeled as correlated RVs to reflect the non-WSSUS properties of V2V channels. The effectiveness of the proposed channel model and simulation method is validated by the measurements in different scenarios.
Vehicle-to-vehicle (V2V) communication system is a promising part of Intelligent Transportation Systems (ITS) , which can provide various of information services, such as the accident notification and the real-time road traffic reports, to improve the safety and the efficiency of the transportation system. The V2V communication system is a type of ad hoc communication network; different cars can exchange information with each other to provide a comfort and “infotainment” environment for driver. The IEEE 802.11p is a new V2V communication standard in the 5.9 GHz Unlicensed National Information Infrastructure (UNII) band, which is an approved amendment to the IEEE 802.11a standard by adding some extended functions to adapt to vehicle communication environment. IEEE 802.11p has been admitted as the communication standard for dedicated short range communications (DSRC) of the ITS in USA .
Channel modeling is the foundation for design and simulation of communication systems, whose accuracy plays a key role in optimizing and improving the performance of V2V communication systems . The propagation characteristics of V2V channels in different scenarios, such as in the urban, the suburban, and the highway, are quite different from those of traditional cellular channels. The characteristics of V2V channels contain (i) both the transmitter and the receiver are moving, (ii) the significant scatterers/reflectors are also changing continuously, (iii) the omnidirectional antennas for both Tx and Rx are at low heights, and (iv) the environment changes dynamically as vehicles move. In traditional channel modeling, the wide-sense stationary uncorrelated scattering (WSSUS) assumption has been widely used to describe random linear time-varying cellular channels . In V2V channels, due to the dynamically changing environment, the statistics of channel changes over moderate time scales (non-WSS); due to the several taps interacting with one-and-the-same object, such as buildings on the road side, the taps at different delays show correlated fading property (non-US). Therefore, the WSSUS assumption is not valid for accurate description of V2V radio communication channels. This specific characteristics of V2V channels are not adequately reproduced by standard TDL models, which leads us to characterize the non-WSSUS properties of V2V channels by developing a non-WSSUS channel model.
To fill the aforementioned gap, this paper proposes a new TDL based non-WSSUS V2V channel model that for the first time considers the amplitude and phase correlations between different taps. Using the relationship between the correlation coefficients of complex Gaussian, Weibull, and Uniform random variables (RVs), correlated taps are generated with arbitrary amplitude and phase correlation coefficients, arbitrary fading parameters, and arbitrary energies. In this model, the correlated taps represent the correlation at different delays, which is caused by the several taps interacting with one-and-the-same object. The model multiplies the switching function with the first-order two-state Markov chains to represent the non-WSS properties, which is caused by the short-time averaged power associated with a reflector fluctuates and the delays of the channel taps drift as the distances between the transmitter and the receiver.
The remainder of this paper is organized as follows. In Section 2, we review the related works of V2V channel models. In Section 3, we review the TDL models and the characteristics of V2V channels. Section 4 proposes a new non-WSSUS V2V channel model, whose taps at different delays are correlated in both the amplitude and phase. Section 5 validates the accuracy of the proposed channel model and evaluates the performances of OFDM systems under the condition of V2V channels.
2. Related Work
There are many literatures focusing on developing accurate models for V2V channels. V2V channel models can be classified as two classes, the geometry-based deterministic models (GBDMs) and the stochastic models, while the stochastic models can be further categorized as the geometry-based stochastic models (GBSMs) and the nongeometrical stochastic models (NGSMs) . The GBDMs characterize typical channel parameters in a deterministic approach. A GBDM based on ray-tracing methods for V2V systems is the pioneering work by Maurer et al. [6, 7]; this model aims at developing a site-specific and actual physical radio propagation process for a given environment. However, the model cannot easily be generalized to more scenarios . The GBSMs are used to simulate the channel by predefined stochastic distribution of effective scatters around the transmitter and the receiver. The GBSMs can be further divided into the regular-shaped GBSMs and the irregular-shaped GBSMs. The simplest regular-shaped GBSM is the two-ring model, which is proposed by Akki and Haber . The two-ring and ellipse model constructs the received signal as a sum of the line-of-sight, single, and double bounced rays with different energies . The classic two-ring model assumes that all scatters are placed in the horizontal plane, while in reality the scatters and reflectors are in different heights. Then the two-ring model is extended to the three-dimensional model, which places the scatters in the two-concentric cylinder around the transmitter and the receiver . In the irregular-shaped GBSMs model, mobile discrete scatters, static discrete scatters, and diffuse scatters are randomly generated . Due to the static nature of the geometry, the GBSMs can generalize to more scenarios; however, the model cannot represent the non-WSSUS properties of the V2V channels. NGSMs determine the model parameters in a completely statistic manner with a certain delay and Doppler shift. In NGSMs, the impulse response of the V2V channel is modeled with taps at certain delays; hence the channel is modeled based on the tapped-delay line (TDL) structure . The sum-of-sinusoids (SoS) simulation method can be used to simulate the V2V channel models . Furthermore, the wideband channel model and channel statistics are discussed in [14, 15]. Among all kinds of channel models, NGSMs are widely used to mimic non-WSSUS for V2V channels due to its low complexity and acceptable accuracy.
IEEE 802.11p employs the 6-tap and 12-tap NGSMs [16, 17] as basic channel model; each tap consists of several paths, and each path has different types of Doppler spectrum (flat shape, round shape, classic 3 dB shape, and classic 6 dB shape). This allows almost arbitrary Doppler spectrum to be synthesized for each tap. However, these models are based on the WSSUS assumption, since they assume that the different taps in TDL structure are independent and uncorrelated to each other. David and coworkers describe the non-WSS property of channel by modeling multipath component persistence via Markov chains . Over moderate time scales, the scatters and reflectors will be “Death and Birth” caused by dynamic blockage or obstacles nearby; the model describes the non-WSS phenomenon by random on/off process. Furthermore, the model generates multivariate Weibull variables with arbitrary fading parameters, arbitrary average fading powers, and arbitrary amplitude correlations to represent the non-WSSUS properties, while none of the above models considers the complex correlations between different taps, including both the amplitude correlations and the phase correlations. To accurately describe the channel characteristics of V2V communication systems, the V2V channel should be modeled as non-WSSUS process.
3. Channel Models
In this section, we first review the traditional TDL based models and discuss the non-WSS and non-US properties of V2V channels. Then the TDL based V2V channel model is proposed to represent the non-WSSUS properties.
3.1. TDL Based Models
The TDL based models are originated from the assumption that the resolvable taps of channels can be approximated as discrete taps with different delays . The models are used to represent the WSSUS linear randomly time-varying channel. The channel impulse response (CIR) of the TDL models can describe completely the channel characteristics, which can be denoted as where represents the index of taps, represents the th resolved amplitude, the exponential term represents the phase of the th tap, and represents the express delay of the th tap. is the radian carrier frequency; , where is the Doppler frequency induced by the motion of the Tx and Rx, denotes the relative velocity, is the aggregate phase angle of all components arriving in the th tap, and is the speed of light. The TDL model is illustrated in Figure 1.
In V2V channels, due to the dynamic environment and large obstacles, the channels are non-WSSUS. Traditional TDL models assume that the taps are independent, while this assumption cannot be hold in V2V channels.
3.2. Non-WSSUS Model
In V2V channels, the transmitter, receiver, reflectors, and scatters around are mobile, so the propagation environment changes dynamically. The significant characteristics of V2V channels are non-WSS and non-US. Non-WSS property describes the channel stochastic characteristics changed with time scale, which represents the time domain characteristics. While non-US property reflects the impacts of correlation on different paths, which represents the delay domain characteristics. In V2V channels, the scatterers and the reflectors will be “Death and Birth” over moderate time scales, which can be modeled by the multipath component persistence as a random on/off process. It can be described as the first-order two-state Markov chains  where is the state transition matrix, is the transition probability of transition from state to state , the state represents the existence of taps, and state represents the disappearance of taps. The represents the steady-state matrix and represents the steady-state probability of state , which satisfies , , and .
The state decision interval of “Death and Birth” state is defined as channel coherence time , where is the maximum Doppler frequency. Compared to traditional TDL models, a switching function is proposed to reflect the non-WSS property.
In V2V channels, the taps at different delays are correlated, due to the several taps interacting with one-and-the same object, such as buildings on the road side. As shown in Figure 2, the correlated complex Gaussian RVs can be generated with any desired correlation coefficients.
The process to generate correlated complex Gaussian RVs with correlation from pairs of uncorrelated Gaussian RVs is as follows. (i) Generate uncorrelated complex Gaussian RVs, which are denoted as . (ii) Determine the coloring matrix by Cholesky decomposition of the correlation matrix , where is the Hermitian transpose of . (iii) Generate correlated complex Gaussian RVs using . The channel correlation property can be realized with correlated taps.
According to the wireless propagation schemes between the transmitter and the receiver in V2V channels, the receivers receive the signal primarily from a particular direction through a narrow beamwidth. On the other hand, the scatters around may lead to different impacts of signal propagation toward the receiver. In the general environment of nonisotropic scattering, which corresponds to directional signal reception, the receivers receive the signal only from particular directions . The uniform distribution is not suitable for the angle distribution of the taps in V2V channels, while the Laplace distribution is suitable for the description of the signal from a particular direction to the receiver. In this paper, the Laplace distribution is used to represent the angle distribution, which is an extension of the classic Uniform distribution model where is the location parameter and which is the scale parameter referred to as the diversity.
In V2V channels, the Weibull distribution is a better fit over all empirical data than the Nakagami distribution . The taps amplitude statistics can be modeled as the flexible Weibull distribution where is the fading parameter. When , the Weibull distribution can be transformed to the Rayleigh distribution. With being increased, the channel quality becomes better; when is large enough, it means that there exits line-of-sight components. The average power is given by where is the expected value, is the gamma function, and is the power parameter. Rayleigh distribution can be transformed to Weibull distribution through the variable transformation .
3.3. Proposed Channel Model
In this section, a TDL based non-WSSUS V2V channel model is proposed. The channel model is defined as CIR , where the amplitude follows the Weibull distribution, the phase follows the Laplace distribution, and is the switching function. As shown in Figure 3, the model is used to generate correlated taps with arbitrary amplitude and phase correlation coefficients, fading parameters, and energy parameters.
The complex Gaussian correlation coefficients can be derived with the desired Weibull correlation coefficients , Uniform correlation coefficients , and energy parameters . Define as complex Gaussian RVs, where the amplitude follows Rayleigh distribution, and the phase follows Uniform distribution. The element is the th row and th column of complex Gaussian correlation coefficients matrix , which is defined as where and are complex Gaussian RVs and denotes variance.
Since the complex random variables are generated from the independent Gaussian random variables, we can prove that the amplitude and phase of the complex random variables are independent on each other; then the first term of the numerator in (6) is given by
The Rayleigh RVs are and , respectively; using the transformation between Rayleigh RVs and Weibull RVs, the expected value is where is the joint pdf of the bivatiate Weibull distribution, which is given by  the following: where is the zero-order modified Bessel function of the first kind. This joint bivariate Weibull distribution assumes that the correlation coefficients of the Rayleigh and Weibull RVs are identical. Even for , the joint pdf can be given by (9).
The Uniform distribution RVs are and ; then the second term of (7) is given as follows:
To calculate the joint distribution, a general family of bivariate joint distribution is proposed by Morgenstern’s family found in  with specified marginal distributions and as where is the association parameter between and . Then we can derive the bivariate uniform distribution which belongs to the Morgenstern’s family.
Consider the following: where and are uniformly distributed RVs in the interval . The correlation coefficients for the bivariate uniform distribution under study are given by ; that is, the correlation coefficients cannot exceed .
Substite the general expression for the joint Weibull distribution (9) into (8) to get the expectation of . Use the classic expression of the joint Uniform distribution (12) into (10) to get the exponential expectation . Then calculating (7), we can derive the element of the complex Gaussian correlation coefficients matrix where is the traditional Gaussian hypergeometric function. It is ideal to realise that the function is related to the Weibull RVs power and the correlation coefficients and .
As shown in Figure 2, generate correlated complex Gaussian RVs with the complex Gaussian correlation coefficients matrix and the RVs of uncorrelated complex Gaussian . The coloring matrix can be generated by the Cholesky decomposition ; then the correlated complex Gaussian RVs are derived by .
Separate the correlated complex Gaussian RVs into the amplitude parts and the phase parts, where the amplitude parts follow Rayleigh distribution and the phase parts follow Uniform distribution.
Using the convenient transformation between Weibull and Rayleigh RVs, correlated Weibull RVs can be derived from the Rayleigh RVs.
The method for generating correlated Laplace RVs is to derive from the correlated Uniform RVs, because the convenient transformation from Uniform RVs to Laplace RVs. Assuming that are uniformly distributed in the interval , then we can find the Laplace distribution with parameters and as follows: where is function.
Multiply the switching function , which is decided by the first-order two-state Markov chains illustrated in (2).
Then a TDL based non-WSSUS channel model is proposed by , which considers the correlation between taps both in amplitude and phase.
4. Simulation Results and Analysis
In this section, we validate the accuracy of the TDL based non-WSSUS V2V channel model and evaluate the performances of the OFDM system under the proposed channel model condition. Table 1 provides the channel parameters for 10 MHz bandwidth V2V channels . The measurement scenarios are classified as the urban-antenna inside car (UIC), the urban-antenna outside car (UOC), the Small City (S), the open-area low traffic density (OLT), and the open-area high traffic density (OHT), which are conducted in five cities in different time and under different traffic conditions. The sum of all tap energies multiplied by their steady-state probabilities equals to unity, and the energy is concentrated in the first tap. The Weibull shape factors represent the fading property of the Weibull distribution. For the tap persistence process, we use the method of modeling an on/off process by Markov chains. is the probability of transition from the state to the state , where the states and describe the existence of taps, means the tap is not existing, means the tap is existing, and is the number of the tap. represents the probability of state and is the probability of state . The statistics of V2V channels change over time; the taps may be disappeared by the obstacles; then we use the random on/off process by Markov chains to describe the “Death and Birth” process.
Tables 2, 3, 4, 5 and 6, describe the matrices of the tap correlation coefficients for the scenarios: UIC, UOC, Small City, OHT, and OLT. Since the matrices of the correlation coefficients are symmetric about the diagonal, the tables are simplified as follows. The upper triangular part is phase correlation coefficients, which is generated randomly. The lower triangular part is amplitude correlation coefficients, which is measured from the scenarios: UIC, UOC, Small City, OHT, and OLT.
The procedure to develop a TDL based V2V channel model which considers the non-WSSUS channel characteristics is given as follows. (i) The correlated taps are generated with the amplitude and phase correlation coefficients and average energies, which represents the correlation at different delays (non-US). (ii) The amplitude and phase are separated from complex correlated taps; then the correlated Weibull RVs can be transformed from Rayleigh RVs and the correlated Laplace RVs can be transformed from Uniform RVs. (iii) Multiply the switching function with the first-order two-state Markov chains; the random on/off process represents the changing of fading process (non-WSS). The taps of the TDL based V2V channel model are correlated, where the amplitude of the generating taps follows Weibull distribution and the phase follows Laplace distribution. Then we calculate that the correlation coefficients matrix of the generated Weibull RVs is and the correlation coefficients matrix of the generated Uniform RVs is
A commonly used method to evaluate the performance is the mean-square error () between the estimated and the desired data
To evaluate the accuracy of the TDL based V2V channel model, we employ the UOC channel parameters to generate non-WSSUS channel model. The desired correlation coefficients matrix is in Table 3, and the generated correlation coefficients matrix is and . The between the desired correlation coefficients matrix and the generated correlation coefficients matrix can be used to show the performance of the TDL based non-WSSUS V2V channel model. The in matrix form is where the matrix is the of the amplitude correlation coefficients matrix and the matrix is the of the phase correlation coefficients matrix. The specifies that the generated taps have a good approximation to the desired correlated taps.
The accuracy of the proposed channel model with respect to measurements is compared with the basis of the simulated power delay profiles (PDPs). In Figure 4, we compare the PDPs of the proposed channel model with measurements in scenarios: UIC, UOC, Small City, OHT, and OLT. The spread delay in UIC scenario is more severe than others, because the UIC scenario is antenna inside the car, which leads to more obstacles in channel propagation. In OLT scenario, which is low traffic density, the scattering and reflection caused by moving cars are less than other scenarios, the spread delay is short, and the energies are centered in the first tap. The simulation results of PDPs are consistent with the measurements, which validates that the proposed channel model can accurately represent the non-WSSUS properties of V2V channels.
Furthermore we evaluate the performances of the OFDM systems under the proposed channel model condition. In simulation, the data rate is 20 Mb/s, the subcarrier number is , the cyclic prefix number is , and channel bandwidth is under 10 MHz. The simulation of the TDL based non-WSSUS channel model in OFDM system is always compared on the basis of bit error ratio (BER). Comparing the BER of simulating five scenarios and AWGN channels in Figure 5.
The presented results have been calculated with a SNR of −5 and 30 dB; it is clear that OLT scenario channel condition is the best and UIC scenario channel condition is the worst. The reason is that OLT is open-area low traffic density scenario; the obstacles such as moving cars are less than other scenarios; UIC is antenna inside the car; the propagation paths are all through the transmitter and receiver cars, which causes severe fading and strong reflections. Due to the antenna inside the car, the channel condition of the UIC scenario presents a large spread delay and achieves the highest BER. The low traffic scenario presents a short delay spread and no further multipath contributions thus achieving a small BER.
In this paper, a TDL based channel model that characterized the non-WSSUS properties of V2V communications is proposed. The taps with different delays are correlated both in amplitude and phase. Utilizing the relationship between the correlation coefficients of complex Gaussian, Weibull, and Uniform RVs, the correlated taps are generated with arbitrary amplitude and phase correlation coefficients, fading parameters, and energy parameters. Using the to evaluate the accuracy of the proposed channel model, the generated correlation coefficients are close to the desired correlation coefficients. The simulation results of the V2V channel PDPs show that the proposed channel model can simulate the V2V channel accurately and effectively.
This work is supported by the National Natural Science Foundation of China under Grant 61222105, the Fundamental Research Funds for the Central Universities under Grant 2013YJS005, Beijing Municipal Natural Science Foundation under Grant 4112048, Project of State Key Lab under Grant no. RCS2012ZT013, key Project for Railway Ministry of China under Grant 2012X008-A, interdiscipline cooperation projects of the New-Star of Science and Technology supported by Beijing Metropolis under Grant no. xxhz201201, and the Key Grant Project of Chinese Ministry of Education under Grant 313006.
- ITS Project, 2007, http://www.its.dot.gov/index.htm.
- IEEE P802.11p/D2.01, “Standard for wireless local area networks providing wireless communication,” in Proceedings of the 5th European Congress Exhibition Intelligent Transportation System Services, Hannover, Germany, June 2005.
- J. Wang, S. Jia, and J. Song, “Signal vector based detection scheme for spatial modulation,” IEEE Communications Letters, vol. 16, no. 1, pp. 19–21, 2012.
- P. A. Bello, “Characterization of randomly time-variant linear channels,” IEEE Transactions on Communications, vol. 11, no. 4, pp. 360–393, 1963.
- C.-X. Wang, X. Cheng, and D. Laurenson, “Vehicle-to-vehicle channel modeling and measurements: recent advances and future challenges,” IEEE Communications Magazine, vol. 47, no. 11, pp. 96–103, 2009.
- J. Maurer, T. Fügen, T. Schäfer, and W. Wiesbeck, “A new inter-vehicle communications (IVC) channel model,” in Proceedings of the IEEE 60th Vehicular Technology Conference (VTC '04), pp. 9–13, September 2004.
- J. Maurer, T. Fgen, M. Porebska, T. Zwick, and W. Wisebeck, “A ray-optical channel model for mobile to mobile communications,” in Proceedings of the 4th MCM COST, vol. 2100, pp. 6–8, 2008.
- J. Wang, S. Jia, and J. Song, “Generalised spatial modulation system with multiple active transmit antennas and low complexity detection scheme,” IEEE Transactions on Wireless Communications, vol. 11, no. 4, pp. 1605–1615, 2012.
- A. S. Akki and F. Haber, “A statistic model for mobile-to-mobile land communication channel,” IEEE Transactions on Vehicular Technology, vol. 35, no. 1, pp. 2–7, 1986.
- X. Cheng, C.-X. Wang, D. I. Laurenson, S. Salous, and A. V. Vasilakos, “An adaptive geometry-based stochastic model for non-isotropic MIMO mobile-to-mobile channels,” IEEE Transactions on Wireless Communications, vol. 8, no. 9, pp. 4824–4835, 2009.
- A. G. Zajić and G. L. Stüber, “Three-dimensional modeling and simulation of wideband MIMO mobile-to-mobile channels,” IEEE Transactions on Wireless Communications, vol. 8, no. 3, pp. 1260–1275, 2009.
- J. Karedal, F. Tufvesson, N. Czink et al., “A geometry-based stochastic MIMO model for vehicle-to-vehicle communications,” IEEE Transactions on Wireless Communications, vol. 8, no. 7, pp. 3646–3657, 2009.
- X. Cheng, Q. Yao, C. X. Wang et al., “An improved parameter computation method for a MIMO V2V rayleigh fading channel simulator under non-isotropic scattering environments,” IEEE Communications Letters, vol. 17, no. 2, pp. 265–2268, 2013.
- X. Cheng, C. X. Wang, B. Ai, and H. Aggoune, “Envelope level crossing rate and average fade duration of non-isotropic vehicle-to-vehicle ricean fading channels,” IEEE Transactions on Intelligent Transportation, no. 99, pp. 1–11.
- X. Cheng, Q. Yao, M. Wen, C. X. Wang, L. Song, and B. Jiao, “Wideband channel modeling and ICI cancellation for vehicle-to-vehicle communication systems,” IEEE Journal on Selected Areas in Communications, vol. 31, no. 9.
- G. Acosta and M. A. Ingram, “Model development for the wideband expressway vehicle-to-vehicle 2.4 GHz channel,” in Proceedings of the IEEE Wireless Communications and Networking Conference (WCNC '06), pp. 1283–1288, April 2006.
- G. Acosta-Marum and M. A. Ingram, “Six time- and frequency- selective empirical channel models for vehicular wireless LANs,” IEEE Vehicular Technology Magazine, vol. 2, no. 4, pp. 4–11, 2007.
- I. Sen and D. W. Matolak, “Vehicle-vehicle channel models for the 5-GHz band,” IEEE Transactions on Intelligent Transportation Systems, vol. 9, no. 2, pp. 235–245, 2008.
- E. Kreyszig, Advanced Engineering Mathematics, Wiley, New York, NY, USA, 8th edition, 1988.
- J. Fühl, J.-P. Rossi, and E. Bonek, “High-resolution 3-d direction-of-arrival determination for urban mobile radio,” IEEE Transactions on Antennas and Propagation, vol. 45, no. 4, pp. 672–682, 1997.
- N. C. Sagias and G. K. Karagiannidis, “Gaussian class multivariate Weibull distributions: theory and applications in fading channels,” IEEE Transactions on Information Theory, vol. 51, no. 10, pp. 3608–3619, 2005.
- N. L. Johnson and S. Kotz, Distributions in Statistics: Continuous Multivariate Distributions, Wiley, New York, NY, USA, 1972.