Propagation graph



Propagation graphs are a mathematical modelling method for radio propagation channels. A propagation graph is a signal flow graph in which vertices represent transmitters, receivers or scatterers. Edges in the graph model propagation conditions between vertices. Propagation graph models were initially developed by Troels Pedersen, et al. for multipath propagation in scenarios with multiple scattering, such as indoor radio propagation. It has later been applied in many other scenarios.

Mathematical definition
A propagation graph is a simple directed graph $$\mathcal G = (\mathcal V, \mathcal E)$$ with vertex set $$\mathcal V$$ and edge set $$\mathcal E$$.

The vertices models objects in the propagation scenario. The vertex set $$\mathcal V$$ is split into three disjoint sets as $$\mathcal V = \mathcal V_t \cup \mathcal V_r \cup\mathcal V_s$$ where $$\mathcal V_t $$ is the set of transmitters, $$\mathcal V_r$$ is the set of receivers and $$\mathcal V_s $$ is the set of objects named "scatterers".

The edge set $$\mathcal E $$ models the propagation models propagation conditions between vertices. Since $$\mathcal G$$ is assumed simple, $$\mathcal E \subset \mathcal V^2$$ and an edge may be identified by a pair of vertices as $$e = (v,v')$$ An edge $$e = (v,v')$$ is included in $$\mathcal E$$ if a signal emitted by vertex $$v$$ can propagate to $$v'$$. In a propagation graph, transmitters cannot have incoming edges and receivers cannot have outgoing edges.

Two propagation rules are assumed
 * A vertex sums the signals impinging via its ingoing edges and remits a scaled version it via the outgoing edges.
 * Each edge $$e=(v,v')$$ transfers the signal from $$v$$ to $$v'$$ scaled by a transfer function.

The definition of the vertex gain scaling and the edge transfer functions can be adapted to accommodate particular scenarios and should be defined in order to use the model in simulations. A variety of such definitions have been considered for different propagation graph models in the published literature.

The edge transfer functions (in the Fourier domain) can be grouped into transfer matrices as where $$f$$ is the frequency variable.
 * $$\mathbf D(f)$$ the direct propagation from transmitters to receivers
 * $$\mathbf T(f)$$ transmitters to scatterers
 * $$\mathbf R(f)$$ scatterers to receivers
 * $$\mathbf B(f)$$ scatterers to scatterers,

Denoting the Fourier transform of the transmitted signal by $$\mathbf X(f)$$, the received signal reads in the frequency domain $$\mathbf Y (f) = \mathbf  D(f) \mathbf X (f) + \mathbf R (f)\mathbf T (f) \mathbf X (f) +  \mathbf R (f)\mathbf B(f) \mathbf T (f) \mathbf X (f) +\mathbf R (f)\mathbf B^2(f) \mathbf T (f) \mathbf X (f) + \cdots$$

Transfer function
The transfer function $$\mathbf H(f)$$ of a propagation graph forms an infinite series $$ \begin{align} \mathbf H(f) &= \mathbf D(f)+ \mathbf R (f)[ \mathbf I+ \mathbf B(f) + \mathbf B(f)^{2} + \cdots ] \mathbf T (f)\\ &= \mathbf D(f)+ \mathbf R (f) \sum_{k=0}^\infty \mathbf B(f)^k \mathbf T(f) \end{align} $$ The transfer function is a Neumann series of operators. Alternatively, it can be viewed pointwise in frequency as a geometric series of matrices. This observation yields a closed form expression for the transfer function as $$\mathbf H(f) = \mathbf D(f) + \mathbf R(f) [\mathbf I - \mathbf B(f)]^{-1} \mathbf T(f),\qquad \rho(\mathbf B(f))<1 $$ where $$\mathbf I$$ denotes the identity matrix and $$\rho(\cdot)$$ is the spectral radius of the matrix given as argument. The transfer function account for propagation paths irrespective of the number of 'bounces'. The series is similar to the Born series from multiple scattering theory.

The impulse responses $$\mathbf h(\tau)$$ are obtained by inverse Fourier transform of $$\mathbf H(f)$$

Partial transfer function
Closed form expressions are available for partial sums, i.e. by considering only some of the terms in the transfer function. The partial transfer function for signal components propagation via at least $$K$$ and at most $$L$$ interactions is defined as where $$\mathbf H_k(f) = \begin{cases} \mathbf D(f),& k=0\\ \mathbf R(f) \mathbf B^{k-1}(f) \mathbf T(f), & k = 1,2,3,\ldots \end{cases} $$ Here $$k$$ denotes the number of interactions or the bouncing order.



The partial transfer function is then $$\mathbf H_{K:L}(f) = \begin{cases} \mathbf D(f) + \mathbf R(f) [\mathbf I-\mathbf B^L(f)] \cdot [\mathbf I-\mathbf B(f)]^{-1} \cdot \mathbf T(f), & K = 0\\ \mathbf R(f) [\mathbf B^{K-1}(f)-\mathbf B^L(f)] \cdot [\mathbf I-\mathbf B(f)]^{-1} \cdot \mathbf T(f), & \text{otherwise}.\\ \end{cases} $$ Special cases:
 * $$\mathbf H_{0:\infty}(f) = \mathbf H(f) $$: Full transfer function.
 * $$\mathbf H_{1:\infty}(f) = \mathbf R(f) [\mathbf I-\mathbf B(f)]^{-1} \mathbf T(f) $$: Inderect term only.
 * $$\mathbf H_{0:L}(f)$$: Only terms with $$L$$ or fewer bounces are kept ($$L$$-bounce truncation).
 * $$\mathbf H_{L+1:\infty}(f)$$: Error term due to an $$L$$-bounce truncation.

One application of partial transfer functions is in hybrid models, where propagation graphs are employed to model part of the response (usually the higher-order interactions).

The partial impulse responses $$\mathbf h_{K:L}(\tau) $$ are obtained from $$\mathbf H_{K:L}(f) $$ by the inverse Fourier transform.

Propagation graph models
The propagation graph methodology have been applied in various settings to create radio channel models. Such a model is referred to as a propagation graph model. Such models have been derived for scenarios including
 * Unipolarized inroom channels. The initial propagation graph models  were derived for unipolarized inroom channels.
 * In a polarimetric propagation graph model is developed for the inroom propagation scenario.
 * The propagation graph framework has been extended in to time-variant scenarios (such as the vehicle-to-vehicle). For terrestrial communications, where relative velocity of objects are limited, the channel may be assumed quasi-static and the static model may be applied at each time step.
 * In a number of works including   propagation graphs have been integrated into ray-tracing models to enable simulation of reverberation phenomena. Such models are referred to as hybrid models.
 * Complex environments including outdoor-to-indoor cases. can be studied by taking advantage of the special structure of propagation graphs for these scenarios. Computation methods for obtaining responses for very complex environments have been developed in
 * The graph model methodology has been used to make spatially consistent MIMO channel models.
 * Several propagation graph models have been published for high-speed train communications.

Calibration of propagation graph models
To calibrate a propagation graph model, its parameters should be set to reasonable values. Different approaches can be taken. Certain parameters can be derived from simplified geometry of the room. In particular, reverberation time can be computed via room electromagnetics. Alternatively, the parameters can ben set according to measurement data using inference techniques such as method of moments (statistics), approximate Bayesian computation., or deep neural networks

Related radio channel model types
The method of propagation graph modeling is related to other methods. Noticeably,
 * Multiple scattering theory
 * Radiosity
 * Ray tracing
 * Geometry-based stochastic channel models (GBSCM)