US 12,451,984 B2
6G pervasive channel modeling method suitable for all frequency bands and all scenarios
Chengxiang Wang, Nanjing (CN); and Zhen Lv, Nanjing (CN)
Assigned to Southeast University, Nanjing (CN); and PURPLE MOUNTAIN LABORATORIES, Nanjing (CN)
Appl. No. 18/696,233
Filed by Southeast University, Nanjing (CN); and PURPLE MOUNTAIN LABORATORIES, Nanjing (CN)
PCT Filed Mar. 19, 2023, PCT No. PCT/CN2023/082380
§ 371(c)(1), (2) Date Mar. 27, 2024,
PCT Pub. No. WO2023/169590, PCT Pub. Date Sep. 14, 2023.
Prior Publication US 2024/0297724 A1, Sep. 5, 2024
Int. Cl. H04B 17/391 (2015.01)
CPC H04B 17/3911 (2015.01) 9 Claims
OG exemplary drawing
 
1. A method for transmitting and receiving signals using a transceiver including pervasively modeling 6G channels configured for frequency bands including sub-6 GHz, millimeter wave, terahertz, and optical wireless frequency bands and scenarios including global-coverage scenarios and full-application scenarios, wherein,
a channel matrix of the pervasively modeling 6G channels is represented as:

OG Complex Work Unit Math
where PL, SH, BL, WE, AL denote large-scale fadings, PL denotes a path loss, SH denotes a shadowing, BL denotes a blockage loss, AL denotes an atmospheric gas absorption loss, WE denotes a weather effect loss, Hs denotes a small-scale fading;
the small-scale fading channel matrix Hs is represented as:

OG Complex Work Unit Math
where MT denotes a number of antenna elements in the transmitter antenna array, MR denotes a number of antenna elements in the receiver antenna array, hqp,fc(t, τ) denotes a channel impulse response between a p-th array element ApT in the transmitter antenna array and a q-th array element AqR in the receiver antenna array at the time instant t, which is represented as a superposition of an LoS component hqp,fcLoS(t, τ) and a NLoS component hqp,fcNLoS(t, τ):

OG Complex Work Unit Math
where KR(t) denotes a Rice factor, hqp,fcLoS(t, τ) and hqp,fcNLoS(t, τ) are respectively represented as:

OG Complex Work Unit Math
where {*}T denotes a transposition operation, fc denotes a carrier frequency, and denote antenna patterns of the array element ApT, for vertica larizations at different frequency bands, Fq,fc,v and Fq,fc,H denote antenna patterns of the array element AqR for vertical and horizontal polarizations at different frequency bands, ϕE,mnR(t) denotes an elevation arrival angle corresponding to a m-th sub-path from a 1st array element of the receiver antenna array to a last-bounce cluster of a n-th path proximity to a receiver side at the time instant t, ϕA,mnR(t) denote an azimuth arrival angle corresponding to the m-th sub-path from the 1st array element of the receiver antenna array to the last-bounce cluster of the n-th path proximity to the receiver side at the time instant t, ϕE,mnT(t) denotes an elevation departure angle corresponding to the m-th sub-path from a 1st array element of the transmitter antenna array to a first-bounce cluster of the n-th path proximity to a transmitter side at the time instant t, ϕA,mnT(t) denotes an azimuth arrival angle corresponding to the m-th sub-path from the 1st array element of the transmitter antenna array to the first-bounce cluster of the n-th path proximity to the transmitter side at the time instant t, κmn(t) denotes a cross polarization power ratio, μ denotes a co-polar imbalance, ϕA,LT(t) and ϕE,LT(t) denote an azimuth departure angle and an elevation departure angle corresponding to an LOS path from A1T to A1R at the time instant t, ϕA,LR(t) and ϕE,LT(t) denote an azimuth arrival angle and an elevation arrival angle corresponding to the LoS path from A1T to A1R at the time instant t, θLVV, θLHH, θmnVV, θmnVH, θmnHV and θmnHH are random phases uniformly distributed over (0, 2π],

OG Complex Work Unit Math
denotes a Faraday rotation angle, a unit of fc in which the Faraday rotation angle is calculated herein GHz, Pqp,mnfc(t) denotes a power of an m-th sub-path in a n-th path from A1T to A1R at a NLoS condition, τqpL(t) denotes a delay of the LoS path at the time instant t,

OG Complex Work Unit Math
dqp(t) denotes a vector distance between the transmitter antenna array APT and the receiver antenna array AqR at the time instant t, c denotes a speed of light, τqp,mn(t) denotes a delay of the m-th sub-path in the n-th path from A1T to A1R at the time instant t, Pqp,mn,fc(t) denotes a power of the m-th sub-path in the n-th path from A1T to A1R at the time instant t, δ denotes Dirac function, and τ denotes a time delay, all of above parameters are time-varying parameters, further, steps for generating the channel matrix H are specifically as follows:
S1, setting propagation scenarios and propagation conditions for the transceiver, and determining a carrier frequency, an antenna array type, a layout of the transceiver and a motion trajectory of the transceiver;
S2, generating, according to standard channel models, the path loss, the shadowing, an oxygen absorption, and blockage effect loss;
S3, generating, according to positions and motion conditions of the transceiver, large-scale parameters with spatial consistency for a delay spread (DS) and 4 angle spreads, wherein the large-scale parameters includes SH, the DS, an azimuth spread of arrival (ASA), an azimuth spread of departure (ASD), an elevation spread of arrival (ESA), an elevation spread of departure (ESD), a Rice factor (KR) and a cross-polarization ratio (XPR);
S4, generating scatterers following an ellipsoid Gaussian scattering distribution, calculating, according to geographical location information of the transceiver and the scatterers, delays, angles and powers of clusters, and generating channel coefficients; and
S5, updating, according to movements of the transceiver and birth-death processes of the clusters, the large-scale parameters and the small-scale parameters; and generating new channel coefficients, wherein
a transmitter communicates with a receiver according to the generated channel matrix.