METHOD FOR DESIGNING TIME-DOMAIN NON-STATIONARY V2V MIMO COMMUNICATION CHANNEL EMULATOR

Abstract

A method for designing a time-domain non-stationary V2V MIMO communication channel emulator includes determining basic parameters for the V2V MIMO communication channel; generating a V2V 2D time-domain non-stationary communication channel environment, by using a MATLAB, that is, the numbers of the scatterers and the positions of the scatterers and the like; importing parameters generated in the previous step into a hardware simulation platform to calculate communication channel parameters for clusters, such as an angle distribution and a power distribution, writing a Verilog code for running, and eventually calculating to obtain a channel impulse response of the time-domain non-stationary V2V MIMO communication channel; and comparing with a statistical characteristic of a theoretical communication channel model, and designing an appropriate hardware diagram of a communication channel emulator. The method supports the simulation of time-domain non-stationary V2V MIMO communication channel, filling the gap in the field of communication channel emulators.

Claims

1. A method for designing a time-domain non-stationary V2V MIMO communication channel emulator, wherein the method comprises following steps: Step S1, determining basic parameters for the V2V MIMO communication channel; Step S2, generating, by using a MATLAB, a V2V 2D time-domain non-stationary communication channel environment, specifically including: a number of scatterers and positions of the scatterers, a random phase of a NLoS path, an angle spread of the NLoS path, a sine function lookup table, and an arctangent function lookup table; Step S3, importing parameters generated in Step S2 into a hardware simulation platform to calculate communication channel parameters for clusters, wherein the communication channel parameters include an angle distribution and an amplitude distribution, writing a Verilog code for running, and eventually calculating to obtain a channel impulse response of the time-domain non-stationary V2V MIMO communication channel; and Step S4, comparing with a statistical characteristic of a theoretical communication channel model, and designing an appropriate hardware diagram of a communication channel emulator.

2. The method for designing the time-domain non-stationary V2V MIMO communication channel emulator according to claim 1, wherein in Step S1, in a geometry-based stochastic model for the time-domain non-stationary V2V MIMO communication channel, the basic parameters for the V2V MIMO communication channel include a number of a simulation time point, a simulation time interval, a position of a transmitter, a position of a receiver, a velocity of the transmitter, a velocity of the receiver, a Rice factor, an angle spread coefficient, a total link power, a LoS path, a power ratio of a SB path, a power ratio of a DB path, a number of initial scatterers, a velocity of the initial scatterers, a number of sub-paths in each cluster, a generation rate for scatterers, a disappearance rate for the scatterers, a motion ratio of the scatterers, a relative coordinate of a receiver antenna and a relative coordinate of a transmitter antenna.

3. The method for designing the time-domain non-stationary V2V MIMO communication channel emulator according to claim 1, wherein Step S2 specifically includes: Step S201, generating the numbers of the scatterers and the positions of the scatterers, and expressing an average survival probability P.sub.remain of effective scatterers within a time period t as $P_{remain}=\exp \left\{-{}_R\left[P_m\left(.Math.{\stackrel{.fwdarw.}{v}}_{S_1}.Math.+.Math.{\stackrel{.fwdarw.}{v}}_{S_2}.Math.\right)t+.Math.{\stackrel{.fwdarw.}{v}}_{Rx}-{\stackrel{.fwdarw.}{v}}_{Tx}.Math.t\right]\right\},$ where .sub.R denotes a generation rate of the scatterers, P.sub.m denotes a motion percentage, {right arrow over (v)}.sub.S.sub.1, and {right arrow over (v)}.sub.S.sub.2 denote an average velocity of scatterers S.sub.1 and an average velocity of scatterers S.sub.2, respectively, {right arrow over (v)}.sub.Tx and {right arrow over (v)}.sub.Rx denote a velocity of the transmitter and a velocity of the receiver, respectively; denoting a probability function P(N.sub.new(t)=k) of a number of newly generated scatterers as $P\left(N_{new}\left(t\right)=k\right)=\frac{{}^k}{k!}{e}^{-},k=0,1,.Math.,\mathrm{and}$ $=\frac{{}_G}{{}_R}\left(1-P_{remain}\right),$ where denotes a parameter for a Poisson distribution probability distribution function; .sub.G denotes a disappearance rate of the scatterers, k denotes a parameter; then, an average number of the scatterers E{N(t)} is: $E\left\{N\left(t\right)\right\}=\frac{{}_G}{{}_R},$ eventually calculating to obtain a number of the scatterers within each time period t, and generating, through a MATLAB, the numbers of the scatterers within each time period t, and sequentially storing the number in a file with a coe suffix; uniformly and randomly generating, in a rectangular area around the receiver and the transmitter, position coordinates of the scatterers; and generating, through the MATLAB, the position coordinates of the scatterers, and sequentially storing the position coordinates in the file with a coe suffix; Step S202, generating the random phase of a NLoS path, wherein the random phase of the NLoS path follows a uniform distribution in [,), and generating, through the MATLAB, the random phase of the NLoS path, and sequentially storing the random phase in the file with a coe suffix; Step S203, generating the angle spread of the NLoS path, and expressing a calculation formula as ${}_{AoD}^m\left(t\right)=AS\left({}_{AoD}\right)Y_{AoD}^m,$ ${}_{AoA}^m\left(t\right)=AS\left({}_{AoA}\right)Y_{AoA}^m,$ where .sub.AoD.sup.m(t) and .sub.AoA.sup.m(t) denote an angle spread of an AoD and an angle spread of an AoA, respectively, AS(.sub.AoD) and AS(.sub.AoA) denote an angle spread coefficient for the AoD and an angle spread coefficient for the AoA, respectively, and Y.sub.AoD.sup.m and Y.sub.AoA.sup.m follow a standard normal distribution N(0,1), respectively, and generating, through the MATLAB, the angle spread of the AoD and the angle spread of the AoA, and sequentially storing the angle spread of the AoD and the angle spread of the AoA in the file with a coe suffix; Step S204, generating the sine function lookup table, uniformly sampling 65536 points within one cycle to be symmetric about a y-axis, amplifying function values to 4096 times, generating, through the MATLAB, the sine function lookup table, and sequentially storing the table in the file with a coe suffix; and Step S205, generating the arctangent function lookup table, uniformly sampling 524288 points around an origin point with an interval $\frac{1}{65536}$ between each sampling point to be symmetric about the y-axis, amplifying function values to $\frac{2^{15}}{}$ times, generating, through the MATLAB, the arctangent function lookup table, and sequentially storing the table in the file with a coe suffix.

4. The method for designing the time-domain non-stationary V2V MIMO communication channel emulator according to claim 1, wherein in Step S3, the angle distributions of the clusters are expressed as .sub.AoD.sup.LoS(t), .sub.AoA.sup.LoS(t), .sub.AoD.sup.i,n(t)and .sub.AoA.sup.i,n(t), the amplitude distribution is expressed as H.sub.qp.sup.LoS(t), H.sub.qp.sup.SB.sup.i(t) and H.sub.qp.sup.DB(t); wherein an angle of the clusters is determined by the coordinate of the transmitter, the coordinate of the receiver, and the coordinates of the scatterers; the amplitude is determined by the Rice factor, the total power, a proportion of a SB ray and a DB ray to a total scattering power on the NLoS paths, a number of an i-th class scatterers at a time instant t, and the number of sub-paths in each cluster of the NLoS path, Step S3 specifically includes: Step S301, generating the AoD of all paths and the AoA of all paths, and expressing calculation formulas of the AoD and the AoA as ${}_{AoD}^{LoS}\left(t\right)=\arctan \frac{R{x}^y\left(t\right)-T{x}^y\left(t\right)}{R{x}^x\left(t\right)-T{x}^x\left(t\right)},$ ${}_{\mathrm{AoA}}^{LoS}\left(t\right)=\arctan \frac{T{x}^y\left(t\right)-R{x}^y\left(t\right)}{T{x}^x\left(t\right)-R{x}^x\left(t\right)},$ ${}_{\mathrm{AoD}}^{i,n}\left(t\right)=\arctan \frac{S_{i,n}^y\left(t\right)-{\mathrm{Tx}}^y\left(t\right)}{S_{i,n}^x\left(t\right)-{\mathrm{Tx}}^x\left(t\right)},$ ${}_{AoA}^{i,n}\left(t\right)=\arctan \frac{S_{i,n}^y\left(t\right)-R{x}^y\left(t\right)}{S_{i,n}^x\left(t\right)-R{x}^x\left(t\right)},$ where S.sub.i,n.sup.x(t), Tx.sup.x(t) and Rx.sup.x(t) denote a horizontal coordinate of a n-th scatterer in an i-th class(i=1,2,3), a horizontal coordinate of the transmitter, and a horizontal coordinate of the receiver, respectively, and S.sub.i,n.sup.y(t), Tx.sup.y(t) and Rx.sup.y(t) denote a vertical coordinate of the n-th scatterer in the i-th class, a vertical coordinate of the transmitter, and a vertical coordinate of the receiver, respectively; Step S302, generating a sub-path angle of the NLoS path, and expressing calculation formulas of the sub-path angle as ${}_{AoD}^{i,n,m}\left(t\right)={}_{AoD}^m\left(t\right)+{}_{\mathrm{AoD}}^{i,n},$ ${}_{\mathrm{AoA}}^{i,n,m}\left(t\right)={}_{AoA}^m\left(t\right)+{}_{\mathrm{AoA}}^{i,n},$ where .sub.AoD.sup.i,n,m(t) and .sub.AoA.sup.i,n,m(t) denote an AoD of an m-th sub-path and an AoA of an m-th sub-path, respectively, .sub.AoD.sup.m(t) and .sub.AoA.sup.m(t) denote an angle spread value for the AoD of the m-th sub-path and an angle spread value for the AoA of the m-th sub-path, respectively, and .sub.AoD.sup.i,n and .sub.AoA.sup.i,ndenote an average value for AoDs of clusters passing though the n-th scatterer in the i-th class and an average value for AoAs of the clusters passing though the n-th scatterer in the i-th class, respectively; Step S303, generating a time delay value, and expressing calculation formulas of the time delay value as: ${}^{LoS}\left(t\right)=\frac{\sqrt{{\left[R{x}^x\left(t\right)-T{x}^x\left(t\right)\right]}^2+{\left[R{x}^y\left(t\right)-T{x}^y\left(t\right)\right]}^2}}{c}$ ${}_n^{{\mathrm{SB}}_i}\left(t\right)=\frac{\sqrt{\begin{array}{c}{\left[S_{i,n}^x\left(t\right)-{\mathrm{Tx}}^x\left(t\right)\right]}^2+{\left[S_{i,n}^y\left(t\right)-{\mathrm{Tx}}^y\left(t\right)\right]}^2+\\ {\left[S_{i,n}^x\left(t\right)-{\mathrm{Rx}}^x\left(t\right)\right]}^2+{\left[S_{i,n}^y\left(t\right)-{\mathrm{Rx}}^y\left(t\right)\right]}^2\end{array}}}{c}$ ${}_{n_1,n_2}^{\mathrm{DB}}\left(t\right)=\frac{\sqrt{\begin{array}{c}{\left[S_{1,n}^x\left(t\right)-{\mathrm{Tx}}^x\left(t\right)\right]}^2+{\left[S_{1,n}^y\left(t\right)-{\mathrm{Tx}}^y\left(t\right)\right]}^2+{\left[S_{1,n}^x\left(t\right)-S_{2,n}^x\left(t\right)\right]}^2+\\ {\left[S_{1,n}^y\left(t\right)-S_{2,n}^y\left(t\right)\right]}^2+{\left[S_{2,n}^x\left(t\right)-{\mathrm{Rx}}^x\left(t\right)\right]}^2+{\left[S_{2,n}^y\left(t\right)-{\mathrm{Rx}}^y\left(t\right)\right]}^2\end{array}}}{c}$ where, .sup.LoS(t) denotes a time delay value for a LoS path, .sub.n.sup.SB.sup.i(t) denotes a time delay value for a SB path cluster passing through the n-th scatterer in the i-th class, .sub.n.sub.1.sub.,n.sub.2.sup.DB(t) denotes a time delay value for a DB path cluster passing through a n.sub.1-th scatterer in the first class and a n.sub.2-th scatterer in the second class, and c denotes a velocity of light; Step S304, generating Doppler phase values of .sup.LoS(t), .sub.n,m.sup.SB.sup.i(t) and .sub.n.sub.1.sub.,n.sub.2.sub.,m.sup.DB(t), and expressing calculation formulas of the Doppler phase values as $f^{LoS}\left(t\right)=f_c\left(\frac{{\stackrel{.fwdarw.}{v}}_{Tx}^T{\stackrel{}{r}}_{Tx}^{LoS}\left(t\right)}{c}+\frac{{\stackrel{.fwdarw.}{v}}_{Rx}^T{\stackrel{}{r}}_{Rx}^{LoS}\left(t\right)}{c}\right),$ where f.sup.LoS(t) denotes a Doppler frequency of the LoS path, f.sub.c denotes a carrier frequency; and {circumflex over (r)}.sub.Tx.sup.LoS(t)=(cos.sub.AoD.sup.LoS(t), sin.sub.AoD.sup.LoS(t)) ${\widehat{r}}_{\mathrm{Rx}}^{\mathrm{LoS}}\left(t\right)=\left(\cos {}_{\mathrm{AoA}}^{\mathrm{LoS}}\left(t\right),\sin {}_{\mathrm{AoA}}^{\mathrm{LoS}}\left(t\right)\right)$ ${}^{\mathrm{LoS}}\left(t\right)={}_{t_0}^{t}2f^{\mathrm{LoS}}\left(t^{}\right){\mathrm{dt}}^{},$ where {circumflex over (r)}.sub.Tx.sup.LoS(t) and {circumflex over (r)}.sub.Rx.sup.LoS(t) denote a unit vector for the AoD of the LoS path and a unit vector for the AoA of the LoS path, respectively, and .sup.LoS(t) denotes a Doppler phase of the LoS path; and $f_{n,m}^{{\mathrm{SB}}_i}\left(t\right)=f_c\left(\frac{{\stackrel{.fwdarw.}{v}}_{n_i,\mathrm{Tx}}^T{\widehat{r}}_{\mathrm{Tx},n,m}^{\mathrm{SB}}\left(t\right)}{c}+\frac{{\stackrel{.fwdarw.}{v}}_{n_i,\mathrm{Rx}}^T{\widehat{r}}_{\mathrm{Tx},n,m}^{\mathrm{SB}}\left(t\right)}{c}\right),$ where f.sub.n,m.sup.SB.sup.i(t) denotes a Doppler frequency of an m-th sub-path of the SB path passing through the n-th scatterer in the i-th class, and ${\stackrel{.fwdarw.}{v}}_{n_i,\mathrm{Rx}}={\stackrel{.fwdarw.}{v}}_{\mathrm{Rx}}-{\stackrel{.fwdarw.}{v}}_{n_i}$ ${\widehat{r}}_{\mathrm{Tx},n,m}^{{\mathrm{SB}}_i}\left(t\right)=\left(\cos {}_{\mathrm{AoD}}^{i,n,m}\left(t\right),\sin {}_{\mathrm{AoD}}^{i,n,m}\left(t\right)\right)$ ${\widehat{r}}_{\mathrm{Rx},n,m}^{{\mathrm{SB}}_i}\left(t\right)=\left(\cos {}_{\mathrm{AoA}}^{i,n,m}\left(t\right),\sin {}_{\mathrm{AoA}}^{i,n,m}\left(t\right)\right)$ ${}_{n,m}^{{\mathrm{SB}}_i}\left(t\right)={}_{t_0}^{t}2f_{n,m}^{{\mathrm{SB}}_i}\left(t^{}\right){\mathrm{dt}}^{},$ where {right arrow over (v)}.sub.n.sub.i.sub.,Tx and {right arrow over (v)}.sub.n.sub.i.sub.,Rx denote a relative velocity of the transmitter with a n-th scatterer and a relative velocity of the receiver with the n-th scatterer, respectively; {circumflex over (r)}.sub.Tx,n,m.sup.SB.sup.i(t) and {circumflex over (r)}.sub.Rx,n,m.sup.SB.sup.i(t) denote a unit vector for an AoD of the m-th sub-path of the SB path passing through the n-th scatterer in the i-th class, and a unit vector for an AoA of the m-th sub-path of the SB path passing through the n-th scatterer in the i-th class, respectively, and .sub.n,m.sup.SB.sup.i(t) denotes a Doppler phase of the m-th sub-path of the SB path passing through the n-th scatterer in the i-th class; $f_{n_1,n_2,m}^{\mathrm{DB}}\left(t\right)=f_c\left(\frac{{\stackrel{.fwdarw.}{v}}_{n_1,\mathrm{Tx}}^T{\widehat{r}}_{\mathrm{Tx},n_1,m}^{\mathrm{DB}}\left(t\right)}{c}+\frac{{\stackrel{.fwdarw.}{v}}_{n_2,\mathrm{Rx}}^T{\widehat{r}}_{\mathrm{Rx},n_2,m}^{\mathrm{DB}}\left(t\right)}{c}\right),$ where f.sub.n.sub.1.sub.,n.sub.2.sub.,m.sup.DB(t) denotes a Doppler frequency of a m-th sub-path of the DB path passing through a n.sub.1-th scatterer in the first class and a n.sub.2-th scatterer in the second class; and ${\stackrel{.fwdarw.}{v}}_{n_1,\mathrm{Tx}}={\stackrel{.fwdarw.}{v}}_{\mathrm{Tx}}-{\stackrel{.fwdarw.}{v}}_{n_1}$ ${\stackrel{.fwdarw.}{v}}_{n_2,\mathrm{Rx}}={\stackrel{.fwdarw.}{v}}_{\mathrm{Rx}}-{\stackrel{.fwdarw.}{v}}_{n_2}$ ${\widehat{r}}_{\mathrm{Tx},n_1,m}^{\mathrm{DB}}\left(t\right)=\left(\cos {}_{\mathrm{AoD}}^{1,n_1,m}\left(t\right),\sin {}_{\mathrm{AoD}}^{1,n_1,m}\left(t\right)\right)$ ${\widehat{r}}_{\mathrm{Rx},n_2,m}^{\mathrm{DB}}\left(t\right)=\left(\cos {}_{\mathrm{AoA}}^{2,n_2,m}\left(t\right),\sin {}_{\mathrm{AoA}}^{2,n_2,m}\left(t\right)\right)$ ${}_{n_1,n_2,m}^{\mathrm{DB}}\left(t\right)={}_{t_0}^{t}2f_{n_1,n_2,m}^{\mathrm{DB}}\left(t^{}\right){\mathrm{dt}}^{},$ where {circumflex over (r)}.sub.Tx,n.sub.1.sub.,m.sup.DB(t), and {circumflex over (r)}.sub.Rx,n.sub.2.sub.,m.sup.DB(t) denote a unit vector for an AoD of the m-th sub-path of a cluster of the DB path passing through the n.sub.1-th scatterer in the first class and the n.sub.2-th scatterer in the second class, and a unit vector for an AoA of the m-th sub-path of the cluster of the DB path passing through the n.sub.1-th scatterer in the first class and the n.sub.2-th scatterer in the second class, respectively, and .sub.n.sub.1.sub.,n.sub.2.sub.,m.sup.DB(t)denotes a Doppler phase of the m-th sub-path of the cluster of the DB path passing through the n.sub.1-th scatterer in the first class and the n.sub.2-th scatterer in the second class; Step S305, generating an antenna phase value, and expressing a calculation formula of the antenna phase value as: a relative coordinate {right arrow over (d)}.sub.Tx of a transmitter antenna being: ${\stackrel{.fwdarw.}{d}}_{\mathrm{Tx}}={\left[d_{\mathrm{Tx}}^z,d_{\mathrm{Tx}}^y\right]}^T,$ where d.sub.Tx.sup.x denotes a horizontal ordinate of {right arrow over (d)}.sub.Tx and dry denotes a vertical ordinate of {right arrow over (d)}.sub.Tx; a phase difference .sub.Tx(t) of the transmitter antenna being: ${}_{\mathrm{Tx}}\left(t\right)=2f_c\frac{{\widehat{r}}_{\mathrm{tx}}^T\left(t\right).Math.{\stackrel{.fwdarw.}{d}}_{\mathrm{Tx}}}{c},$ where {circumflex over (r)}.sub.Tx(t) denotes a unit vector for an AoD at the time instant t; a relative coordinate {right arrow over (d)}.sub.Rx of the receiver antenna being: ${\stackrel{.fwdarw.}{d}}_{\mathrm{Rx}}={\left[d_{\mathrm{Rx}}^x,d_{\mathrm{Rx}}^y\right]}^T,$ where d.sub.Rx.sup.x, denotes a horizontal ordinate of {right arrow over (d)}.sub.Rx, d.sub.Rx.sup.y, denotes a vertical ordinate of {right arrow over (d)}.sub.Rx, a phase difference .sub.Rx(t) of the receiver antenna being: ${}_{\mathrm{Rx}}\left(t\right)=2f_c\frac{{\widehat{r}}_{\mathrm{tx}}^T\left(t\right).Math.{\stackrel{.fwdarw.}{d}}_{\mathrm{Rx}}}{c},$ where {circumflex over (r)}.sub.Rx(t) denotes a unit vector for an AoA at the time instant t; and a total antenna phase value (t) being: $\left(t\right)={}_{\mathrm{Tx}}\left(t\right)+{}_{\mathrm{Rx}}\left(t\right);$ Step S306, generating amplitude values, expressing calculation formulas of the amplitude values as: $H^{\mathrm{LoS}}\left(t\right)=\sqrt{\frac{K_{\mathrm{qp}}{P}_{\mathrm{qp}}}{K_{\mathrm{qp}}+1}}$ $H^{{\mathrm{SB}}_i}\left(t\right)=\sqrt{\frac{{}_{{\mathrm{SB}}_i}{P}_{\mathrm{qp}}}{\left(K_{\mathrm{qp}}+1\right)N_i\left(t\right)M}}$ $H_{\mathrm{qp}}^{\mathrm{DB}}\left(t\right)=\sqrt{\frac{{}_{\mathrm{DB}}{P}_{\mathrm{qp}}}{\left(K_{\mathrm{qp}}+1\right)N_1\left(t\right)N_2\left(t\right)M}},$ where, H.sub.qp.sup.LoS(t) denotes a channel impulse response amplitude value for a LoS path between a q-th receiving antenna and a p-th transmitting antenna, H.sub.qp.sup.SB.sup.i(t) denotes a channel impulse response amplitude value for a sub-path of a SB path passing through the scatterer in the i-th class located between the q-th receiving antenna and the p-th transmitting antenna, H.sub.qp.sup.DB(t) denotes a channel impulse response amplitude value for the sub-path of the DB path between the q-th receiving antenna and the p-th transmitting antenna, K.sub.qp denotes a Rice factor of a p-q link, P.sub.qp denotes a total power of the p-q link, .sub.SB.sub.i and .sub.SB.sub.i denote the proportion of the total scattering power of the SB ray and the DB ray on the NLoS path, respectively, N.sub.i(t) denotes a number of the scatterer in the i-th class at the time instant t, and M denotes a number of the sub-paths in each cluster of the NLoS path; and Step S307, generating channel impulse response, and expressing a calculation formula of the CIR as $h_{\mathrm{qp}}^{\mathrm{LoS}}\left(t,\right)=H^{\mathrm{LoS}}\left(t\right)e^{-j2f_c{}^{\mathrm{LoS}}\left(t\right)}{e}^{j{}^{\mathrm{Los}}\left(t\right)}{e}^{j{}^{\mathrm{LoS}}\left(t\right)}\left(-{}^{\mathrm{LoS}}\left(t\right)\right)$ $h_{\mathrm{qp}}^{\mathrm{SB}}\left(t,\right)=\underset{i=1}{\overset{3}{.Math.}}\underset{n_i=1}{\overset{N_i\left(t\right)}{.Math.}}{H}^{{\mathrm{SB}}_i}\left(t\right)e^{-j2f_c{}_n^{{\mathrm{SB}}_i}\left(t\right)}{e}^{j{}_{n,m}^{{\mathrm{SB}}_i}\left(t\right)}{e}^{j{}_{n,m}^{{\mathrm{SB}}_i}\left(t\right)}{e}^{j{}_{n,m}^{{\mathrm{SB}}_i}\left(t\right)}\left(-{}_n^{{\mathrm{SB}}_i}\left(t\right)\right)$ $h_{\mathrm{qp}}^{\mathrm{DB}}\left(t,\right)=\underset{n_1,n_2=1}{\overset{N_1\left(t\right),N_2\left(t\right)}{.Math.}}{H}_{n_1,n_2,m}^{\mathrm{HB}}\left(t\right)e^{-j2f_c{}_{n_1,n_2}^{\mathrm{DB}}\left(t\right)}{e}^{j{}_{n_1,n_2,m}^{\mathrm{DB}}\left(t\right)}{e}^{j{}_{n_1,n_2,m}^{\mathrm{DB}}\left(t\right)}{e}^{j{}_{n_1,n_2,m}^{\mathrm{DB}}\left(t\right)}\left(-{}_{n_1,n_2}^{{\mathrm{DB}}_i}\left(t\right)\right)$ $h_{\mathrm{qp}}\left(t,\right)=h_{\mathrm{qp}}^{\mathrm{LoS}}\left(t,\right)+h_{\mathrm{qp}}^{\mathrm{SB}}\left(t,\right)+h_{\mathrm{qp}}^{\mathrm{DB}}\left(t,\right),$ where h.sub.qp.sup.LoS(t,) denotes a channel impulse response for the LoS path between the q-th receiving antenna and the p-th transmitting antenna, h.sub.qp.sup.SB(t,) denotes the channel impulse response for the SB path between the q-th receiving antenna and the p-th transmitting antenna, h.sub.qp.sup.DB(t,) denotes a channel impulse response for the DB path between the q-th receiving antenna and the p-th transmitting antenna, f.sub.c denotes a carrier center frequency, denotes a time delay, denotes a Doppler phase, denotes an antenna phase difference, .sub.n,m.sup.SB.sup.i(t) and .sub.n.sub.1.sub.,n.sub.2.sub.,m.sup.DB(t) denote a random phase of a m-th SB sub-path passing through the n-th scatterer in the i-th class and a random phase of a m-th SB sub-path passing through the n.sub.1-th scatterer in the first class and a n.sub.2-th scatterer in the second class, respectively; and importing from the file with a coe suffix and exporting a calculated channel impulse response to a txt file.

5. The method for designing the time-domain non-stationary V2V MIMO communication channel emulator according to claim 1, wherein in Step S4, formulas for the statistical characteristic specifically include: Step S401, expressing a calculation formula of a time autocorrelation function TACF curve as: $r_{\mathrm{qp},q^{}{p}^{}}=E\left\{{h_{\mathrm{qp}}\left(t,\right)\left[h_{q^{}{p}^{}}\left(t,\right)\right]}^{*}\right\},$ where h.sub.qp(t,) denotes a channel impulse response between a q-th receiving antenna and a p-th transmitting antenna in a case where a time period is t and a time delay is , h.sub.qp(t,) denotes a channel impulse response between a q-th receiving antenna and a p-th transmitting antenna in a case where a time period is t and a time delay is ;(.Math.)* denotes a conjugate complex of(.Math.); Step S402, expressing a calculation formula of a spatial cross-correlation function SCCF curve as $r_{\mathrm{qp}}\left(t\right)=E\left\{{h_{\mathrm{qp}}\left(t,\right)\left[h_{\mathrm{qp}}\left(t+t,\right)\right]}^{*}\right\},$ and Step S403, expressing a calculation formula of a delay power spectral density(PSD)curve as $\left(,{}_{\mathrm{AoA}}^{i,n}\right)={.Math.h_{\mathrm{qp}}\left(t,\right).Math.}^2{.Math.}_{{}_{\mathrm{AoA}}^{i,n}}.$

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0042] FIG. 1 illustrates a flow chart diagram of a method for designing a time-domain non-stationary V2V MIMO communication channel emulator provided in Embodiment 1 of the present disclosure.

[0043] FIG. 2 illustrates a schematic diagram of a geometric stochasticity-based model for the time-domain non-stationary V2V MIMO communication channel provided in Embodiment 1 of the present disclosure.

[0044] FIG. 3 illustrates a schematic diagram of a hardware block diagram of a time-domain non-stationary V2V MIMO communication channel emulator provided in Embodiment 1 of the present disclosure.

[0045] FIG. 4 illustrates a schematic diagram of a TD sum-of-frequency-modulation(SoFM) module in a hardware block diagram of a time-domain non-stationary V2V MIMO communication channel emulator provided in the Embodiment 1 of the present disclosure.

[0046] FIG. 5 illustrates a schematic diagram of TACF verification results provided in the Embodiment 1 of the present disclosure, where f.sub.c=500 MHZ, K.sub.qp=1, .sub.SB.sub.i =0.2, i=1, 2, 3, .sub.DB=0.4, P.sub.qp=30 dBm, M=8, {right arrow over (v)}.sub.Tx=(10,0), {right arrow over (v)}.sub.Rx=(5,0)and {right arrow over (v)}.sub.S.sub.1(2)=(1,0).

[0047] FIG. 6 illustrates a schematic diagram of SCCF verification results provided in the Embodiment 1 of the present disclosure, where f.sub.c=500 MHZ, K.sub.qp=1, .sub.SB.sub.i =0.2, i=1, 2, 3, .sub.DB=0.4, P.sub.qp=30 dBm, M=8, {right arrow over (v)}.sub.Tx=(10,0), {right arrow over (v)}.sub.Rx=(5,0)and {right arrow over (v)}.sub.S.sub.1(2)=(1,0).

[0048] FIG. 7 illustrates a schematic diagram of SCCF simulation results provided in the Embodiment 1 of the present disclosure, where f.sub.c=500 MHZ, K.sub.qp=1, .sub.SB.sub.i =0.2, i=1, 2, 3, .sub.DB=0.4, P.sub.qp=30 dBm, M=8, {right arrow over (v)}.sub.Tx=(10,0), {right arrow over (v)}.sub.Rx=(5,0)and {right arrow over (v)}.sub.S.sub.1(2)=(1,0), .sub.G=0.08/m, .sub.R=0.03/m, and P.sub.m=0.6.

[0049] FIG. 8 illustrates a schematic diagram of simulation results for delayed PSD hardware platform provided in the Embodiment 1 of the present disclosure, where f.sub.c=500 MHZ, K.sub.qp=1, .sub.SB.sub.i =0.2, i=1, 2, 3, .sub.DB=0.4, P.sub.qp=30 dBm, M=8, {right arrow over (v)}.sub.Tx=(10,0), {right arrow over (v)}.sub.Rx=(5,0)and {right arrow over (v)}.sub.S.sub.1(2)=(1,0), .sub.G=0.08/m, .sub.R=0.03/m, and P.sub.m=0.6.

DETAILED DESCRIPTION OF THE EMBODIMENTS

[0050] In order to clarify the objectives, the technical solutions, and the advantages of the embodiments of the present disclosure to be clearer, the technical solutions in the embodiments of present disclosure will be clearly and completely described in conjunction with the accompanying drawings. Obviously, the described embodiments are one part of the embodiments of the present disclosure, not all of them. Based on the embodiments in the present disclosure, all other embodiments obtained by a person skilled in the art without creative labor fall within the protection scope of the present disclosure.

Embodiment 1

[0051] With reference to FIGS. 1 to 8, provided in this embodiment is a method for designing a time-domain non-stationary V2V MIMO communication channel emulator. The method for designing the communication channel emulator specifically includes the following steps.

[0052] In Step S1, basic parameters for the V2V MIMO communication channel are determined.

[0053] Specifically, in this embodiment, c a method for establishing a geometry-based stochastic channel model (GBSM). Firstly, the application scenario is determined as V2V, and the initial coordinates of the transmitter and the receivers are determined as (8,8) and (50,50) at the same time. The carrier center frequency is set as 500 MHz, the antenna spacing distance between the transmitter and the receiver is set as 1 meter, the elevation angle and the horizontal angle of the transmitting antenna array and the receiving antenna array are set as 0 and 0, the motion velocity of the transmitter and the motion velocity of the receiver are set as 10 m/s and 5 m/s, respectively, and the direction of the motion is that the elevation angle is 0 and the horizontal angle is 0. K.sub.qp=1, .sub.SB.sub.i=0.2, i=1,2,3, .sub.DB=0.4, P.sub.qp=30 dBm, and the number of the sub-paths in each cluster is 8. Non isotropic omnidirectional antenna is adopted by both the transmitting antenna and the receiving antenna.

[0054] In Step S2, a V2V 2D time-domain non-stationary communication channel environment is generated by using a MATLAB, and the environment specifically includes the number of scatterers and positions of the scatterers, a random phase of a NLoS path, an angle spread of the NLoS path, a sine function lookup table, and an arctangent function lookup table.

[0055] Specifically, in this embodiment, the V2V communication channel model is adopted in the model, and the schematic diagram of the specific communication channel model is illustrated in FIG. 2.

[0056] A uniform linear array is adopted at the antenna terminal, which can be placed arbitrarily in a 2D space. For simplicity, only one cluster in each of LoS path, SB path, and DB path is shown in the drawings. Three classes of the scatterers are located on a circle with the transmitter as the center, a circle with the receiver as the center, and an ellipse representing the street environment, respectively, and only one scatterer of all three classes of the scatterers is shown in the drawings. {right arrow over (v)}.sub.S.sub.1(2) denotes a motion velocity of scatterer in the 1.sup.st(2.sup.nd)class.

[0057] More specifically, in this embodiment, generating a V2V 2D time-domain non-stationary communication channel environment specifically includes the following steps.

[0058] In Step S201, the numbers of the scatterers and positions of the scatterers are generated.

[0059] An average survival probability of effective scatterers within the time period t is expressed as

[00024]$P_{\mathrm{remain}}=\exp \left\{-{}_R\left[P_m\left(.Math.{\stackrel{.fwdarw.}{v}}_{S_1}.Math.+.Math.{\stackrel{.fwdarw.}{v}}_{S_2}.Math.\right)t+.Math.{\stackrel{.fwdarw.}{v}}_{\mathrm{Rx}}-{\stackrel{.fwdarw.}{v}}_{\mathrm{Tx}}.Math.t\right]\right\},$

where .sub.R denotes a generation rate of the scatterers, P.sub.m denotes a motion percentage, {right arrow over (v)}.sub.S.sub.1 and {right arrow over (v)}.sub.S.sub.2 denote an average velocity of the scatterers S.sub.1 and an average velocity of S.sub.2, respectively. {right arrow over (v)}.sub.Tx and {right arrow over (v)}.sub.Rx denote a velocity of the transmitter and a velocity of the receiver, respectively. A probability function N.sub.new(t) of the number of newly generated scatterers is expressed as

[00025]$P\left(N_{\mathrm{new}}\left(t\right)=k\right)=\frac{{}^k}{k!}{e}^{-},$$k=0,1,.Math.,$$\mathrm{and}$$=\frac{{}_G}{{}_R}\left(1-P_{\mathrm{remain}}\right)$

where denotes a parameter for a Poisson distribution probability distribution function; and .sub.G denotes a disappearance rate of the scatterers. By combining the above two formulas, the average number of scatterers can be obtained:

[00026]$E\left\{N\left(t\right)\right\}=N\left(t-t\right)P_{\mathrm{remain}}+E\left\{N_{\mathrm{new}}\left(t\right)\right\}=\frac{{}_G}{{}_R}.$

[0060] Eventually, the numbers of the scatterers within each time period are obtained by calculation, which are generated through a MATLAB, and are sequentially stored in a file with a coe suffix. Position coordinates of the scatterers are uniformly and randomly generated in a rectangular area around the receiver and the transmitter, which are generated through the MATLAB, and are sequentially stored in the file with a coe suffix.

[0061] In Step S202, the random phase of a NLoS path is generated.

[0062] The random phase of a NLoS path follows a uniform distribution in [,), which are generated through the MATLAB, and are sequentially stored in the file with a coe suffix.

[0063] In Step S203, the angle spread of the NLoS path is generated.

[0064] Calculation formula of the angle spread of the NLoS path is expressed as

[00027]${}_{AoD}^m\left(t\right)=AS\left({}_{AoD}\right)Y_{AoD}^m,$${}_{AoA}^m\left(t\right)=AS\left({}_{AoA}\right)Y_{AoA}^m,$

where AS(.sub.AoD)and AS(.sub.AoA)denote an angle spread coefficient for the AoD and an angle spread coefficient for the AoA, respectively, and Y.sub.AoD.sup.m and Y.sub.AoA.sup.m follow a standard normal distribution N(0,1), respectively, which are generated through the MATLAB, and are sequentially stored in the file with a coe suffix.

[0065] In Step S204, the sine function lookup table is generated.

[0066] 65536 points are uniformly sampled within one cycle to be symmetric about a y-axis, and function values are amplified to 4096 times, which are generated through the MATLAB, and are sequentially stored in the file with a coe suffix.

[0067] In Step S205, the arctangent function lookup table is generated.

[0068] 524288 points are uniformly sampled around an origin point with an interval 1/65536 between each sampling point to be symmetric about the y-axis, and function values are amplified to 2.sup.15/ times, which are generated through the MATLAB, and are sequentially stored in the file with a coe suffix.

[0069] In Step S3, parameters generated in Step S2 are imported into a hardware simulation platform to calculate the communication channel parameters of the clusters, such as an angle distribution and an amplitude distribution, a Verilog code is written for running, and a channel impulse response of the time-domain non-stationary V2V MIMO communication channel is eventually obtained by calculation.

[0070] In this embodiment, Step S3 specifically includes the following steps.

[0071] In Step S301, the AoD of all paths and the AoA of all paths are generated.

[0072] The calculation formulas of the AoD of all paths and the AoA of all paths are expressed as

[00028]${}_{\mathrm{AoD}}^{\mathrm{LoS}}\left(t\right)=\arctan \frac{{\mathrm{Rx}}^y\left(t\right)-{\mathrm{Tx}}^y\left(t\right)}{{\mathrm{Rx}}^x\left(t\right)-{\mathrm{Tx}}^x\left(t\right)},$${}_{\mathrm{AoD}}^{\mathrm{LoS}}\left(t\right)=\arctan \frac{{\mathrm{Tx}}^y\left(t\right)-{\mathrm{Rx}}^y\left(t\right)}{{\mathrm{Tx}}^x\left(t\right)-{\mathrm{Rx}}^x\left(t\right)},$${}_{\mathrm{AoD}}^{i,n}\left(t\right)=\arctan \frac{S_{i,n}^y\left(t\right)-{\mathrm{Tx}}^y\left(t\right)}{S_{i,n}^x\left(t\right)-{\mathrm{Tx}}^x\left(t\right)},$${}_{\mathrm{AoA}}^{i,n}\left(t\right)=\arctan \frac{S_{i,n}^y\left(t\right)-{\mathrm{Rx}}^y\left(t\right)}{S_{i,n}^x\left(t\right)-{\mathrm{Rx}}^x\left(t\right)},$

where S.sub.i,n.sup.x(t), Tx.sup.x(t) and Rx.sup.x(t) denote a horizontal coordinate of the n-th scatterer in the i-th class(i=1,2,3), a horizontal coordinate of the transmitter, and a horizontal coordinate of the receiver, respectively, and S.sub.i,n.sup.y(t), Tx.sup.y(t) and Rx.sup.y(t) denote a vertical coordinate of the n-th scatterer in the i-th class, a vertical coordinate of the transmitter, and a vertical coordinate of the receiver, respectively.

[0073] In Step S302, a sub-path angle of the NLoS path is generated.

[0074] The calculation formulas of a sub-path angle of the NLoS path are expressed as

[00029]${}_{\mathrm{AoD}}^{i,n,m}\left(t\right)={}_{\mathrm{AoD}}^m\left(t\right)+{}_{\mathrm{AoD}}^{i,n},$${}_{\mathrm{AoD}}^{i,n,m}\left(t\right)={}_{\mathrm{AoA}}^m\left(t\right)+{}_{\mathrm{AoA}}^{i,n},$

where .sub.AoD.sup.m(t) and .sub.AoA.sup.m(t) denote an angle spread value for the AoD of the m-th sub-path and an angle spread value for the AoA of the m-th sub-path, respectively, and .sub.AoD.sup.i,n and .sub.AoA.sup.i,n denote an average value for the AoDs of clusters passing through the n-th scatterer in the i-th class and an average value for the AoAs of the clusters passing through the n-th scatterer in the i-th class, respectively.

[0075] In Step S303, a time delay value is generated.

[0076] The calculation formula of a time delay value is expressed as

[00030]${}^{\mathrm{LoS}}\left(t\right)=\frac{\sqrt{{\left[{\mathrm{Rx}}^x\left(t\right)-{\mathrm{Tx}}^x\left(t\right)\right]}^2+{\left[{\mathrm{Rx}}^y\left(t\right)-{\mathrm{Tx}}^y\left(t\right)\right]}^2}}{c}$${}_n^{{\mathrm{SB}}_i}\left(t\right)=\frac{\sqrt{\begin{array}{c}{\left[S_{i,n}^x\left(t\right)-{\mathrm{Tx}}^x\left(t\right)\right]}^2+{\left[S_{\mathrm{in}}^y\left(t\right)-{\mathrm{Tx}}^y\left(t\right)\right]}^2+\\ {\left[S_{i,n}^x\left(t\right)-{\mathrm{Rx}}^x\left(t\right)\right]}^2+{\left[S_{i,n}^y\left(t\right)-{\mathrm{Rx}}^y\left(t\right)\right]}^2\end{array}}}{c}$${}_{n_{1,}{n}_2}^{\mathrm{DB}}\left(t\right)=\frac{\sqrt{\begin{array}{c}{\left[S_{1,n}^x\left(t\right)-{\mathrm{Tx}}^x\left(t\right)\right]}^2+{\left[S_{1,n}^y\left(t\right)-{\mathrm{Tx}}^y\left(t\right)\right]}^2+\\ {\left[S_{1,n}^x\left(t\right)-S_{2,n}^x\left(t\right)\right]}^2+[S_{1,n}^y\left(t\right)-\\ {S_{2,n}^y\left(t\right)]}^2+[S_{2,n}^x\left(t\right)-\\ {{\mathrm{Rx}}^x\left(t\right)]}^2+{\left[S_{2,n}^y\left(t\right)-{\mathrm{Rx}}^y\left(t\right)\right]}^2\end{array}}}{c}$

where, .sup.LoS(t) denotes a time delay value for a LoS path, .sub.n.sup.SB.sup.i(t) denotes a time delay value for a SB path cluster passing through the n-th scatterer in the i-th class, .sub.n.sub.1.sub.,n.sub.2.sup.DB(t) denotes a time delay value for a DB path cluster passing through the n.sub.1-th scatterer in the first class and the n.sub.2-th scatterer in the second class, and c denotes a velocity of light.

[0077] In Step S304, Doppler phase values are generated.

[0078] The calculation formulas of the Doppler phase values are expressed as

[00031]$f^{LoS}\left(t\right)=f_c\left(\frac{{\stackrel{.fwdarw.}{v}}_{Tx}^T{\stackrel{}{r}}_{Tx}^{LoS}\left(t\right)}{c}+\frac{{\stackrel{.fwdarw.}{v}}_{Rx}^T{\stackrel{}{r}}_{Rx}^{LoS}\left(t\right)}{c}\right),$

where f.sup.LoS(t) denotes a Doppler frequency of the LoS path, f.sub.c denotes a carrier frequency; and {circumflex over (r)}.sub.Tx.sup.LoS(t)=(cos.sub.AoD.sup.LoS(t),sin.sub.AoD.sup.LoS(t))

[00032]${\stackrel{}{r}}_{Rx}^{LoS}\left(t\right)=\left(\cos {}_{AoA}^{Los}\left(t\right),\sin {}_{AoA}^{Los}\left(t\right)\right)$${}^{LoS}\left(t\right)={}_{t_0}^{t}2f^{LoS}\left(t^{}\right){\mathrm{dt}}^{},$

where {circumflex over (r)}.sub.Tx.sup.LoS(t) and {circumflex over (r)}.sub.Rx.sup.LoS(t) denote a unit vector for the AoD of the LoS path and a unit vector for the AoA of the LoS path, respectively, and .sup.LoS(t) denotes a Doppler phase of the LoS path; and

[00033]$f_{n,m}^{S{B}_i}\left(t\right)=f_c\left(\frac{{\stackrel{.fwdarw.}{v}}_{n_i,\mathrm{Tx}}^T{\stackrel{}{r}}_{\mathrm{Tx},n,m}^{\mathrm{SB}}\left(t\right)}{c}+\frac{{\stackrel{.fwdarw.}{v}}_{n_i\mathrm{Rx}}^T{\stackrel{}{R}}_{x,n,m}^{\mathrm{SB}}\left(t\right)}{c}\right),$

f.sub.n,m.sup.SB.sup.i(t) denotes a Doppler frequency of the m-th sub-path of the SB path passing through the n-th scatterer in the i-th class, and {right arrow over (v)}.sub.n,.sub.i.sub.,Tx ={right arrow over (v)}.sub.Tx{right arrow over (v)}.sub.n,.sub.i

[00034]${\stackrel{.fwdarw.}{v}}_{n_i,Rx}={\stackrel{.fwdarw.}{v}}_{Rx}-{\stackrel{.fwdarw.}{v}}_{n_i}$${\stackrel{}{r}}_{\mathrm{Tx},n,m}^{S{B}_i}\left(t\right)=\left(\cos {}_{AoD}^{i,n,m}\left(t\right),\sin {}_{AoD}^{i,n,m}\left(t\right)\right)$${\stackrel{}{r}}_{\mathrm{Rx},n,m}^{S{B}_i}\left(t\right)=\left(\cos {}_{\mathrm{AoA}}^{i,n,m}\left(t\right),\sin {}_{\mathrm{AoA}}^{i,n,m}\left(t\right)\right)$${}_{n,m}^{S{B}_i}\left(t\right)={}_{t_0}^{t}2f_{n,m}^{S{B}_i}\left(t^{}\right){\mathrm{dt}}^{},$

where {right arrow over (v)}.sub.n.sub.i.sub.,Tx and {right arrow over (v)}.sub.n.sub.i.sub.,Rx denote a relative velocity of the transmitter with a n-th scatterer and a relative velocity of the receiver with the n-th scatterer, respectively; {circumflex over (r)}.sub.Tx,n,m.sup.SB.sup.i(t) and {circumflex over (r)}.sub.Rx,n,m.sup.SB.sup.i(t) denote a unit vector for an AoD of the m-th sub-path of the SB path passing through the n-th scatterer in the i-th class, and a unit vector for an AoA of the m-th sub-path of the SB path passing through the n-th scatterer in the i-th class, respectively, and .sub.n,m.sup.SB.sup.i(t) denotes a Doppler phase of the m-th sub-path of the SB path passing through the n-th scatterer in the i-th class;

[00035]$f_{n_1,n_2,m}^{\mathrm{DB}}\left(t\right)=f_c\left(\frac{{\stackrel{.fwdarw.}{v}}_{n_i,\mathrm{Tx}}^T{\stackrel{}{r}}_{\mathrm{Tx},n_1,m}^{\mathrm{DB}}\left(t\right)}{c}+\frac{{\stackrel{.fwdarw.}{v}}_{n_i\mathrm{Rx}}^T{\stackrel{}{r}}_{x,n_2,m}^{\mathrm{DB}}\left(t\right)}{c}\right),$

where f.sub.n.sub.1.sub.,n.sub.2.sub.,m.sup.DB(t) denotes a Doppler frequency of the m-th sub-path of the DB path cluster passing through the n.sub.1-th scatterer in the first class and the n.sub.2-th scatterer in the second class; and

[00036]${\stackrel{.fwdarw.}{v}}_{n_1,\mathrm{Tx}}={\stackrel{.fwdarw.}{v}}_{Tx}-{\stackrel{.fwdarw.}{v}}_{n_1}$${\stackrel{.fwdarw.}{v}}_{n_2,Rx}={\stackrel{.fwdarw.}{v}}_{Rx}-{\stackrel{.fwdarw.}{v}}_{n_2}$${\stackrel{}{r}}_{\mathrm{Tx},n_1,m}^{DB}\left(t\right)=\left(\cos {}_{AoD}^{1,n_1,m}\left(t\right),\sin {}_{AoD}^{1,n_1,m}\left(t\right)\right)$${\stackrel{}{r}}_{\mathrm{Rx},n_2,m}^{DB}\left(t\right)=\left(\cos {}_{AoA}^{2,n_2,m}\left(t\right),\sin {}_{AoA}^{2,n_2,m}\left(t\right)\right)$${}_{n_1,n_2,m}^{DB}\left(t\right)={}_{t_0}^{t}2f_{n_{1,}{n}_2,m}^{\mathrm{DB}}\left(t^{}\right){\mathrm{dt}}^{},$

where {circumflex over (r)}.sub.Tx,n.sub.1.sub.,m.sup.DB(t), and {circumflex over (r)}.sub.Rx,n.sub.2.sub.,m.sup.DB(t) denote a unit vector for an AoD of the m-th sub-path of a cluster of the DB path passing through the n.sub.1-th scatterer in the first class and the n.sub.2-th scatterer in the second class, and a unit vector for an AoA of the m-th sub-path of the cluster of the DB path passing through the n.sub.1-th scatterer in the first class and the n.sub.2-th scatterer in the second class, respectively, and .sub.n.sub.1.sub.,n.sub.2.sub.,m.sup.DB(t)denotes a Doppler phase of the m-th sub-path of the cluster of the DB path cluster passing through the n.sub.1-th scatterer in the first class and the n.sub.2-th scatterer in the second class.

[0079] In Step S305, an antenna phase value is generated.

[0080] The calculation formula of the antenna phase value is expressed as follows.

[0081] A relative coordinate {right arrow over (d)}.sub.Tx of a transmitter antenna is:

[00037]${\stackrel{.fwdarw.}{d}}_{Tx}={\left[d_{\mathrm{Tx}}^x,d_{Tx}^y\right]}^T,$

where d.sub.Tx.sup.x denotes a horizontal ordinate of {right arrow over (d)}.sub.Tx and d.sub.Tx.sup.y denotes a vertical ordinate of {right arrow over (d)}.sub.Tx

[0082] A phase difference .sub.Tx(t) of the transmitter antenna is:

[00038]${}_{Tx}\left(t\right)=2f_c\frac{{\hat{r}}_{\mathrm{Tx}}^T\left(t\right).Math.{\stackrel{.fwdarw.}{d}}_{\mathrm{Tx}}}{c},$

where {circumflex over (r)}.sub.Tx(t) denotes a unit vector of an AoD at the time instant t.

[0083] A relative coordinate {right arrow over (d)}.sub.Rx of the receiver antenna i:

[00039]${\stackrel{.fwdarw.}{d}}_{Rx}={\left[d_{\mathrm{Rx}}^x,d_{Rx}^y\right]}^T,$

where dix denotes a horizontal ordinate of {right arrow over (d)}.sub.Rx, d.sub.Rx.sup.y denotes a vertical ordinate of {right arrow over (d)}.sub.Rx. A phase difference .sub.Tx(t) of the receiver antenna is:

[00040]${}_{Rx}\left(t\right)=2f_c\frac{{\hat{r}}_{\mathrm{Rx}}^T\left(t\right).Math.{\stackrel{.fwdarw.}{d}}_{Rx}}{c},$

where {circumflex over (r)}.sub.Rx(t) denotes a unit vector for an AoA at the time instant t.

[0084] The calculation formula of the total antenna phase value (t) is expressed as:

[00041]$\left(t\right)={}_{\mathrm{Tx}}\left(t\right)+{}_{\mathrm{Rx}}\left(t\right).$

[0085] In Step S306, amplitude values are generated.

[0086] The calculation formulas of the amplitude values are expressed as:

[00042]$H^{LoS}\left(t\right)=\sqrt{\frac{K_{qp}{P}_{qp}}{K_{qp}+1}}$$H^{S{B}_i}\left(t\right)=\sqrt{\frac{{}_{S{B}_i}{P}_{qp}}{\left(K_{qp}+1\right)N_i\left(t\right)M}}$$H_{qp}^{DB}\left(t\right)=\sqrt{\frac{{}_{DB}{P}_{qp}}{\left(K_{qp}+1\right)N_1\left(t\right)N_2\left(t\right)M}},$

where, H.sub.qp.sup.LoS(t) denotes channel impulse response amplitude value for a LoS path between the q-th receiving antenna and the p-th transmitting antenna, H.sub.qp.sup.SB.sup.i(t) denotes a channel impulse response amplitude value for a sub-path of a SB path passing through the scatterer in the i-th class located between the q-th receiving antenna and the p-th transmitting antenna, H.sub.qp.sup.DB(t) denotes a channel impulse response amplitude value for the sub-path of the DB path between the q-th receiving antenna and the p-th transmitting antenna, K.sub.qp denotes a Rice factor of a p-q link, P.sub.qp denotes a total power of the p-q link, .sub.SB.sub.i and .sub.SB.sub.i denote the proportion of the total scattering power of the SB ray and the DB ray on the NLoS path, respectively, N.sub.i(t) denotes the number of the scatterer in the i-th class at the time instant t, and M denotes the number of the sub-paths in each cluster of the NLoS path.

[0087] K.sub.qp, .sub.SB.sub.i and .sub.DB in the model are related to the macro cells, the micro cells, and the micro micro cells, as well as a SB path and a DB path. The model proposed by the present disclosure can adapt to various V2V propagation environments by adjusting the model parameters .sub.SB.sub.i, .sub.DB and the Rice factor K.sub.qp. In the macro cell scenarios, due to a relative large distance between the transmitter and the receiver, the DB rays carry more energy than that of the SB rays(the relative large distance between the transmitter and the receiver leads to a greater independence of an AoD and an AoA, that is, .sub.DB>max{.sub.SB.sub.1,.sub.SB.sub.2}>>.sub.SB.sub.3, and the received signal power mainly comes from the SB rays and the DB rays of the dual loop model. Therefore, the Rice factor K.sub.qp and the energy parameter .sub.SB.sub.3 are extremely little, even close to zero. This means that in the macro cell scenario, it can be characterized by using a dual loop model that ignores the LoS radial component. Compared with the macro cell scenario, in the micro cell scenario and the micro-micro cell scenario, VTD significantly affects the channel characteristics. In consideration of the impacts of VTD on channel statistics, the present disclosure distinguishes the mobile cars and the stationary roadside environments(such as buildings, trees, parked cars, etc.)around the transmitter and the receiver. Therefore, the present disclosure simulates a mobile car by using a dual loop model and describes a stationary roadside environment by using an elliptical model. For low VTD, since the LoS path has a high power, the value for K.sub.qp is large. In addition, the received scattering power mainly comes from the signal reflected by the stationary roadside environment described by the scatterers freely located on the ellipse. Since the mobile car represented by the scatterers located on the dual loop is sparse, it is more likely to be a SB rather than a DB, indicating a validity of .sub.SB.sub.3>max{.sub.SB.sub.1,.sub.SB.sub.2}>.sub.DB is. Under a condition of a high VTD, the value for K.sub.qp is smaller than that in a low VTD scenario, and due to the relative large number of the mobile cars, the DB rays of the dual loop model carry more energy than those of the SB rays of the dual loop model and the elliptical model, that is, .sub.DB>max {.sub.SB.sub.1,.sub.SB.sub.2,.sub.SB.sub.3}. Thus, the micro cell scenario and the micro-micro cell scenario in consideration of VTD can be well characterized by utilizing the combined dual loop model and the elliptical model with a LoS radial component.

[0088] In Step S307, channel impulse response is generated.

[0089] The calculation formula of the channel impulse response is expressed as

[00043]$h_{qp}^{LoS}\left(t,\right)=H^{LoS}\left(t\right)e^{-j2f_c{}^{LoS}\left(t\right)}{e}^{j{}^{{\mathrm{LoS}}_{\left(t\right)}}}{e}^{j{}^{LoS}\left(t\right)}\left(-{}^{LoS}\left(t\right)\right)$$h_{qp}^{SB}\left(t,\right)=\underset{i=1}{\overset{3}{.Math.}}\underset{n_i=1}{\overset{N_i\left(t\right)}{.Math.}}{H}^{S{B}_i}\left(t\right)e^{-j2f_c{}_n^{S{B}_i}\left(t\right)}{e}^{j{}_{n,m}^{S{B}_i}\left(t\right)}{e}^{j{}_{n,m}^{S{B}_i}\left(t\right)}{e}^{j{}_{n,m}^{S{B}_i}\left(t\right)}\left(-{}_n^{S{B}_i}\left(t\right)\right)$$h_{qp}^{DB}\left(t,\right)=\underset{n_1,n_2=1}{\overset{N_1\left(t\right),N_2\left(t\right)}{.Math.}}{H}_{n_1,n_2,m}^{\mathrm{DB}}\left(t\right)e^{-j2f_c{}_{n_1,n_2}^{DB}\left(t\right)}{e}^{j{}_{n_1,n_2,m}^{DB}\left(t\right)}{e}^{j{}_{n_1,n_2,m}^{DB}\left(t\right)}{e}^{j{}_{n_1,n_2,m}^{DB}\left(t\right)}\left(-{}_{n_1,n_2}^{DB}\left(t\right)\right)$$h_{qp}\left(t,\right)=h_{qp}^{LoS}\left(t,\right)+h_{qp}^{SB}\left(t,\right)+h_{qp}^{DB}\left(t,\right),$

where h.sub.qp.sup.LoS(t,) denotes a channel impulse response for the LoS path between the q-th receiving antenna and the p-th transmitting antenna, h.sub.qp.sup.SB(t,) denotes the channel impulse response for the SB path between the q-th receiving antenna and the p-th transmitting antenna, h.sub.qp.sup.DB(t,) is a channel impulse response for the DB path between the q-th receiving antenna and the p-th transmitting antenna, f.sub.c denotes a carrier center frequency, denotes a time delay, denotes a Doppler phase, denotes an antenna phase difference, .sub.n,m.sup.SB.sup.i(t) and .sub.n.sub.1.sub.,n.sub.2.sub.,m.sup.DB(t) denote a random phase of the m-th SB sub-path passing through the n-th scatterer in the i-th class and a random phase of a m-th SB sub-path passing through the n.sub.1-th scatterer in the first class and the n.sub.2-th scatterer in the second class, respectively; the channel impulse response is imported from the file with a coe suffix and a calculated channel impulse response is exported to a txt file.

[0090] In Step S4, a comparison is performed with a statistical characteristic of a theoretical communication channel model, and an appropriate hardware diagram of a communication channel emulator is designed.

[0091] The specific hardware block diagram is illustrated in FIGS. 3 and 4. In this embodiment, Step 4 specifically includes the following steps.

[0092] In Step S401, a calculation formula of a time autocorrelation function TACF curve is expressed as:

[00044]$r_{\mathrm{qp},q^{}{p}^{}}=E\left\{{h_{qp}\left(t,\right)\left[h_{q^{}{p}^{}}\left(t,\right)\right]}^{*}\right\},$

where h.sub.qp(t,)denotes a channel impulse response between the q-th receiving antenna and the p-th transmitting antenna in a case where a time period is t and a time delay is , h.sub.qp(t,) denotes a channel impulse response between the q-th receiving antenna and the p-th transmitting antenna in a case where a time period is t and a time delay is ; (.Math.)* denotes a conjugate complex of(.Math.).

[0093] In Step S402, a calculation formula of a SCCF curve is expressed as

[00045]$r_{qp}\left(t\right)=E\left\{{h_{qp}\left(t,\right)\left[h_{qp}\left(t+t,\right)\right]}^{*}\right\}.$

[0094] In Step S403, a calculation formula of a delay PSD curve is expressed as

[00046]${\left(,{}_{AoA}^{i,n}\right)={.Math.h_{qp}\left(t,\right).Math.}^2.Math.}_{{}_{\mathrm{AoA}}^{i,n}}.$

[0095] In order to verify the correctness of the method provided in this embodiment, experiments are conducted, and specifically lie as follows.

[0096] The TACF between the transmitting antenna and the receiving antenna at different time intervals are verified, and the results are as illustrated in FIG. 5. As the time interval increases, the TACF of the NLoS path is a constant, while the TACF of the SB path and DB path decrease with the increase of the time intervals. And the reason of that is that the time changes cause the changes on the angular distributions of sub-path and the Doppler phase.

[0097] The SCCF between the receiving antennas at different antenna intervals are verified, and the results are as illustrated in FIG. 6. As the antenna interval increases, the SCCF of the NLoS path is a constant, while the SCCF of the SB path and DB path decrease with the increase of the antenna intervals. And the reason of that is that the changes in antenna intervals cause changes on the phase differences of the antenna array in the sub-path.

[0098] The delay PSD of channel impulse response at different antenna intervals are verified, and the results are as illustrated in FIG. 7 and FIG. 8. Observing the drawings, the hardware simulation results match well with the MATLAB simulation results. The coherent lines in FIG. 7 and FIG. 8 represent the delay PSD of the channel impulse response of cluster passing through the scatterers in the first class in the LoS path and SB path. This is because the scatterers in the first class is relative closer to the transmitter and relative farther from the receiver. Therefore, even if one part of the scatterers in the first class is generated or disappears, the AoA is basically the same as the angle from the transmitter to the receiver. The other scatter points in FIG. 7 and FIG. 8 indicate that as the time delay increases, the clusters of other NLoS paths are generating or disappearing. It can be seen from the delay PSD drawings, this embodiment has a characteristics of birth-death process of clusters in time-domain, which is an important characteristic of time-domain non-stationary communication channels.

[0099] In summary, provided in the present disclosure is a method for designing a time-domain non-stationary V2V MIMO communication channel emulator, which can describe the fading characteristics of macro cell, micro cell, and micro-micro cell, as well as the birth and death process of the clusters. The statistical characteristics of the simulations have important reference value for the designs of time-domain non-stationary communication channel emulators.

[0100] The unspecified parts in the present disclosure are all the common sense for a person skilled in the art. The preferred specific embodiments of the present disclosure are described in details in above. It should be understood that various amendments and changes can be made by an ordinary person skilled in the art according to the concept of the present disclosure with no creative efforts. Therefore, all technical solutions that can be obtained by a person skilled in the art on a basis of the prior art according to the concept of the present disclosure through the logical analysis, reasoning, or limited experiments should be within the protection scope determined by the claims.