信道模型：Rayleigh、Rician、卫星→地面：Shadowed-Rician

比较

	C.Loo	Corazza	Lutz
	非静止轨道卫星	非静止轨道卫星	静止轨道卫星
	乡村	所有	所有
	模拟卫星的直升机	实际卫星	实际卫星

1. C.Loo模型和Corazza模型适用于描述非静止轨道卫星信道的传播特性
   Lutz模型适用于描述静止轨道卫星信道的传播特性；
2. C.Loo模型只适用于乡村信道环境，
   而Corazza模型和Lutz模型均适用所有的卫星移动通信信道环境（公路、乡村、郊区和城市）；
3. C.Loo模型拟合参数时依据的实测数据来自于模拟卫星的直升机所发射的信号；
   而Corazza模型和Lutz模型则采用实际卫星所发射信号的实测数据，对卫星移动通信信道传输特性的反映更为真实。

C. Loo模型：直射阴影，多径不阴影

A statistical model for a land mobile satellite link
Chun Loo, “A statistical model for a land mobile satellite link,” in IEEE Transactions on Vehicular Technology, vol. 34, no. 3, pp. 122-127, Aug. 1985, doi: 10.1109/T-VT.1985.24048.

假设接收信号中只有直射信号分量受到阴影遮蔽的作用而多径信号分量不受阴影遮蔽的作用，因此该模型又称为部分阴影信道模型。
C. Loo模型接收信号包络的概率密度函数
$f_r\left(r\right)={\int }_0^{\infty }{f}_r\left(r|z\right)f_z\left(z\right)dz=\frac{r}{{\text{b}}_{\text{0}}\sqrt{2\pi d_0}}{\int }_0^{\infty }\frac{1}{z}\exp \left[-\frac{{\left(\ln z-\mu \right)}^2}{2d_0}-\frac{\left(r^2+z^2\right)}{2b_0}\right]\text{dz }$
z表示直射波，
b0是平均散射多径功率，
μ和d0分别是lnz的均值和方差。

Corazza模型：直射和多径都阴影

假设接收信号中的直射信号分量和多径信号分量同时都受到阴影遮蔽的作用，因此该模型又称为全阴影信道模型
接收信号包络的概率密度函数为
$f_r\left(r\right)={\int }_0^{\infty }{f}_r\left(r|s\right)f_s\left(s\right)\text{ = }\frac{\text{2}\left(\text{k + 1}\right)r}{h\sigma s\sqrt{2\pi }}\exp \left(-k\right){\int }_0^{\infty }\exp \left[-\frac{\left(k+1\right)r^2}{s^2}-\frac{{\left(\ln s-\mu \right)}^2}{2{\left(h\sigma \right)}^2}\right]I_0\left(\frac{2r\sqrt{k\left(k+1\right)}}{s}\right)ds$
S(t)表示阴影慢衰落效应，
h=(ln10)/20，
μ和σ2为lns的均值和方差，
k为Rice因子，
r为接收信号包络。
$$\begin{gathered} \text{K}\left(\alpha \right)\text{ = }{\text{k}}_{\text{0}}\text{ + }{\text{k}}_{\text{1}}\alpha \text{ + }{\text{k}}_{\text{2}}{\alpha }^{\text{2}} \\ \mu \left(\alpha \right)\text{ = }{\mu }_{\text{0}}\text{ + }{\mu }_{\text{1}}\alpha \text{ + }{\mu }_{\text{2}}{\alpha }^{\text{2}}\text{ + }{\mu }_{\text{3}}{\alpha }^{\text{3}} \\ \sigma \left(\alpha \right)\text{ = }{\sigma }_{\text{0}}\text{ + }{\sigma }_{\text{1}}\alpha \end{gathered}$$

Lutz模型：好坏2个状态

E. Lutz, D. Cygan, M. Dippold, F. Dolainsky and W. Papke, “The land mobile satellite communication channel-recording, statistics, and channel model,” in IEEE Transactions on Vehicular Technology, vol. 40, no. 2, pp. 375-386, May 1991, doi: 10.1109/25.289418.

将卫星移动通信信道分为两个状态——“好状态”和“坏状态”。
当信道为“好状态”时接收信号中的直射信号分量和多径信号分量均不受阴影遮蔽的作用；
当信道为“坏状态”时接收信号中只有受到阴影遮蔽作用的多径信号分量且不存在直射信号分量，
因此该模型又称为两状态信道模型。

好状态下接收信号功率S的归一化概率密度函数为
$f_{s\_rician}\left(s\right)=c\exp \left[-c\left(s+1\right)\right]I_0\left(2c\sqrt{s}\right)$
C为归一化因子，c = 1/(2σ2)

坏状态下接收信号功率S的归一化概率密度函数为
$$\begin{gathered} f_{s\_rayl\_LN}\left(\text{s}\right)\text{ = }{\int }_{\text{0}}^{\infty }{\text{f}}_{\text{s}}\left(s|s_0\right)f_{s0}\left(s_0\right)\text{d}{\text{s}}_{\text{0}} \\ \text{ = }\frac{\text{10}}{\sigma \text{ln10}\sqrt{\text{2}\pi }}{\int }_{\text{0}}^{\infty }\frac{\text{1}}{{{\text{s}}_{\text{0}}}^{\text{2}}}\text{exp}\left[\text{ - }\frac{\text{s}}{{\text{s}}_{\text{0}}}-\frac{{\left(10{\log }_{10}{s}_0-\mu \right)}^2}{2{\sigma }^2}\right] \end{gathered}$$
S0 为短时平均接收功率

总的接收信号功率S的概率密度函数为
$f_s\left(\text{s}\right)\text{ = }\left(\text{1 - A}\right){\text{f}}_{\text{s\_rician}}\left(\text{s}\right)\text{ + A}{\text{f}}_{\text{s\_rayl\_LN}}\left(\text{s}\right)$

Rayleigh and Rician 信道生成

2.4.2 Rayleigh and Rician fading
(2.54)
Fundamentals of Wireless Communication

Rayleigh channel (NLoS), Rician channel (LoS)
https://blog.csdn.net/weixin_41181312/article/details/123516627

bandwidth
$B=180kHz$

Noise power spectral density
$N_0=10^{\frac{-170}{10}}mW/Hz$
$N_{0,dBm}=-170dBm/Hz$
Noise power
$P_N=N_{0}B\left(mW\right)=10^{\frac{-170}{10}}\times 180000\left(mW\right)$
$P_{N,dBm}=10{\log }_{10}\left(P_N\right)=N_{0,dBm}+10{\log }_{10}\left(B\right)\left(dBm\right)=-170+10{\log }_{10}\left(180000\right)\left(dBm\right)$

pathloss
$L\left(d,f\right)={\left(\frac{4\pi df}{c}\right)}^2$
$L_{dB}\left(d,f\right)=10{\log }_{10}\left(L\left(d,f\right)\right)\left(dB\right)$

Further advancements for E-UTRA physical layer aspects (Release 9). 3GPP TS 36.814, Mar. 2010.
3GPP TS 36.814
the 3GPP propagation environment [52, Table B.1.2.1-1].
LOS
lossdB=28 + 22*log10(dist_) + 20*log10(freq_/10^9);
NLOS
lossdB=22.7 + 36.7*log10(dist_) + 26*log10(freq_/10^9);

attenuation
$\alpha \left(d,f\right)=\sqrt{L{\left(d,f\right)}^{-1}/P_N}=\sqrt{{10}^{\frac{-L_{dB}\left(d,f\right)}{10}}/{10}^{\frac{P_{N,dBm}}{10}}}=\sqrt{{10}^{\frac{-P_{N,dBm}-L_{dB}\left(d,f\right)}{10}}}$

1 M-antenna transmitter → 1 single-antenna receiver
NLOS
${\mathbf{h}}_{NLOS}=\sqrt{1/2}\left(\mathbf{c}+j\mathbf{e}\right)\in {\mathbb{C}}^{M\times 1}$
均值为0，方差为var的复高斯分布
sqrt(var/2)*(randn(1,K) +j*randn(1,K))
randn产生均值为0，方差为1的正态分布随机数
Rayleigh
${\mathbf{h}}_{Rayleigh}\left(d,f\right)=\alpha \left(d,f\right){\mathbf{h}}_{NLOS}\in {\mathbb{C}}^{M\times 1}$
Rician
${\mathbf{h}}_{Rician}\left(d,f\right)=\alpha \left(d,f\right)\left(\sqrt{\frac{\epsilon }{1+\epsilon }}\mathbf{a}\left(\theta \right)+\sqrt{\frac{1}{1+\epsilon }}{\mathbf{h}}_{NLOS}\right)\in {\mathbb{C}}^{M\times 1}$

Shadowed-Rician 直射径 散射径

优点
计算量小

Although some mathematical models, such as those of Loo, Barts–Stutzman, and Karasawa et al.,have been developed to describe the satellite channel, the shadowed-Rician model proposed in [19] is a popular model, which provides a significantly less computational burden than other channel models.

A. Abdi, W. C. Lau, M. . -S. Alouini and M. Kaveh, “A new simple model for land mobile satellite channels: first- and second-order statistics,” in IEEE Transactions on Wireless Communications, vol. 2, no. 3, pp. 519-528, May 2003, doi: 10.1109/TWC.2003.811182.

复包络
the lowpass-equivalent complex envelope of the stationary narrowband shadowed Rice single model can be written as
$\Re (t)=A(t)\exp [j\alpha (t)]+Z(t)\exp (j{\zeta }_0)$
$\alpha (t)$ is the stationary random phase process with uniform distribution over $[0,2\pi )$, 随机，多天线咋办？
$\zeta 0$ is the deterministic phase of the LOS component. 随机选一个就行，多天线咋办？
The independent stationary random processes $A(t)$ and $Z(t)$, which are also independent of $\alpha (t)$, are the amplitudes of the scatter and the LOS components, following Rayleigh and Nakagami distributions, respectively 根据分布及其参数生成
$p_A(a)=\frac{a}{b_0}\exp \left(\frac{-a^2}{2b_0}\right),a\ge 0$
$p_Z(z)=\frac{2m^m}{\Gamma (m){\Omega }^m}{z}^{2m-1}\exp \left(\frac{-m{z}^2}{\Omega }\right),z\ge 0$
$\Gamma (.)$ is the gamma function,
$2b_0=E[A^2]$ is the average power of the scatter component,
$m=(E[Z^2]{)}^2/\mathrm{Var}[Z^2]\ge 0$ is the Nakagami parameter with $Var[.]$ as the variance,
$\Omega =E[Z^2]$ is the average power of the LOS component.

包络
Let us define the envelope as $R(t)=\left|\Re (t)\right|$
the shadowed Rice PDF for the signal envelope in an LMS channel can be written as
$p_R(r)=E_Z\left[\frac{r}{b_0}\exp \left(-\frac{r^2+Z^2}{2b_0}\right)I_0\left(\frac{Zr}{b_0}\right)\right],r\ge 0$
$EZ[.]$ is the expectation with respect to $Z$
$In(.)$ is the $n$ th-order modified Bessel function of the first kind.
$p_R(r)={\left(\frac{2b_{0}m}{2b_{0}m+\Omega }\right)}^m\frac{r}{b_0}\exp \left(-\frac{r^2}{2b_0}\right){}_1{F}_1\left(m,1,\frac{\Omega r^2}{2b_0(2b_{0}m+\Omega )}\right),r\ge 0$
where ${}_1{F}_1(.,.,.)$ is the confluent hypergeometric function [28].
For $m=0$, (3) simplifies to the Rayleigh PDF in (1), i.e., $(r/b_0)\exp (-r^2/2b_0)$,
for $m=\infty$, it reduces to the Rice PDF $(r/b_0)\exp (-(r^2+\Omega )/2b_0)I_0(\sqrt{\Omega }r/b_0)$.

包络的平方
$S(t)=R^2(t)$ is the instantaneous power, and its PDF can be derived from (3)as
$p_S(s)={\left(\frac{2b_{0}m}{2b_{0}m+\Omega }\right)}^m\frac{1}{2b_0}\exp \left(-\frac{s}{2b_0}\right){}_1{F}_1\left(m,1,\frac{\Omega s}{2b_0(2b_{0}m+\Omega )}\right),s\ge 0$

the probability density function (PDF)
$f_{|h{|}^2}(x)=\alpha \exp (-\beta x{)}_1{F}_1(m;1;\delta x)$
$f_{{|h|}^2}\left(x\right)=\alpha e^{-\beta x}{}_1{F}_1\left(m;1;\delta x\right)$
$\alpha =\frac{1}{2b}{\left(\frac{2bm}{2bm+\Omega }\right)}^m$
$\beta =\frac{1}{2b}$
$\delta =\frac{\Omega }{2b(2bm+\Omega )}$

总结

the confluent hypergeometric function of the first kind [31, Eq. (9.14.1)],
https://blog.csdn.net/gsgbgxp/article/details/126136164
https://www.mathworks.com/help/symbolic/hypergeom.html#bt1nkmw-2
${}_1{F}_1(m;1;\delta x)=e^{\delta x}{\sum }_{k=0}^{m-1}\frac{{(-1)}^k{(1-m)}_k{(\delta x)}^k}{{(k!)}^2}$
the pochhammer symbol
(x)K = x (x + 1) (x + 2) · · · (x − k + 1)

The SAT is equipped with $N_t$ -antennas
$\mathbf{h}=\sqrt{b(\phi )}\tilde{\mathbf{h}}\in {\mathbb{C}}^{N_t\times 1}$
$b(\phi )={\left(\frac{J_1(u)}{2u}+36\frac{J_3(u)}{u^3}\right)}^2$
$u=2.07123\frac{\sin \phi }{\sin {\phi }_{3\mathrm{dB}}}$
$\tilde{\mathbf{h}}=A\exp (j\mathbf{\psi })+Z\exp (j\mathbf{\varphi })$

$G_s(\theta )g=G_s(\theta )|\bar{g}/D{|}^2$
$g=|\bar{g}/D{|}^2$
$\bar{g}=A\exp \left(j\varphi \right)+Z\exp \left(j\phi \right)$
$G_s\left(\theta \right)=G_{tx}{\left(\frac{J_1\left(u\right)}{2u}+36\frac{J_3\left(u\right)}{u^3}\right)}^2$

${\mathbf{h}}_{SJ}=\sqrt{L\frac{G_{t,SJ}{G}_{r,SJ}}{K_{B}TB}}{\mathbf{g}}_{SJ}=\frac{C}{4\pi f{d}_{SJ}}\sqrt{\frac{G_{t,SJ}{G}_{r,SJ}}{K_{B}TB}}{\mathbf{g}}_{SJ}$

$y=\sqrt{L{G}_{s}G(\phi )}hx+n$
$L={\left(\frac{c}{4\pi f_{c}d}\right)}_{}^2$
$G(\phi )=G_j{\left(\frac{J_1\left(u\right)}{2u}+36\frac{J_3\left(u\right)}{u^3}\right)}_{}^2$
$f_{{|h|}^2}\left(x\right)=\alpha e^{-\beta x}{}_1{F}_1\left(m;1;\delta x\right)$

多天线 2018LWC北京理工大学

B. Li, Z. Fei, X. Xu and Z. Chu, “Resource Allocations for Secure Cognitive Satellite-Terrestrial Networks,” in IEEE Wireless Communications Letters, vol. 7, no. 1, pp. 78-81, Feb. 2018, doi: 10.1109/LWC.2017.2755014.
Resource Allocations for Secure Cognitive Satellite-Terrestrial Networks

The SAT is equipped with $N_t$ -antennas
$\mathbf{h}=\sqrt{b(\phi )}\tilde{\mathbf{h}}\in {\mathbb{C}}^{N_t\times 1}$
$b(\phi )={\left(\frac{J_1(u)}{2u}+36\frac{J_3(u)}{u^3}\right)}^2$
$u=2.07123\frac{\sin \phi }{\sin {\phi }_{3\mathrm{dB}}}$
$\tilde{\mathbf{h}}=A\exp (j\mathbf{\psi })+Z\exp (j\mathbf{\varphi })$
where $b(\phi )$ is the corresponding beam gain factor, which is determined by their location.
$\phi$ is the angle between the corresponding receiver and the beam center, 怎么确定？？
${\phi }_{3\mathrm{dB}}=0.4/180\ast \pi$ is the 3-dB angle.
$J_1(\cdot )$ and $J_3(\cdot )$ represent the first-kind Bessel function of order 1 and 3.
$\tilde{\mathbf{h}}\in {\mathbb{C}}^{N_t\times 1}$ denotes the channel fading vector from SAT to the receiver, which include the scattering and the line-of-sight (LOS) components.
$\mathbf{\psi }\in [0,2\pi )$ denotes the stationary random phase 随机，多天线咋办？
$\mathbf{\varphi }$ denotes the deterministic phase of the LOS component. 随机选一个就行，多天线咋办？
$A$ and $Z$ are the amplitudes of the scattering and the LOS components. 根据分布及其参数生成
The beam angles from SAT to PU, to Eve and to SU are set as 0.01°, 0.4° and 0.8°, respectively

多天线 2020TVT

H. Zhang, C. Jiang, J. Wang, L. Wang, Y. Ren and L. Hanzo, “Multicast Beamforming Optimization in Cloud-Based Heterogeneous Terrestrial and Satellite Networks,” in IEEE Transactions on Vehicular Technology, vol. 69, no. 2, pp. 1766-1776, Feb. 2020, doi: 10.1109/TVT.2019.2959933.

多天线 2020TCOMM空军工程大学

Q. Huang, M. Lin, J. -B. Wang, T. A. Tsiftsis and J. Wang, “Energy Efficient Beamforming Schemes for Satellite-Aerial-Terrestrial Networks,” in IEEE Transactions on Communications, vol. 68, no. 6, pp. 3863-3875, June 2020, doi: 10.1109/TCOMM.2020.2978044.

Joint CoMP Transmission for UAV-Aided Cognitive Satellite Terrestrial Networks

For the SAT channel links, the channel coefficients from SAT to SU/PU follow shadowed-Rician (SR) model [13], [33], which can be written as
$\bar{g}=A\exp \left(j\varphi \right)+Z\exp \left(j\phi \right)$

Denote the transmit power of SAT as $Ps$ , which is assumed to be invariant over each time slot.
In addition, the transmit power of BS and UAV within time slot $n$ are respectively denoted as $pb[n]$ and $pu[n]$ .
Therefore, the achievable rate of SU (bps/Hz) within time slot $n$ can be expressed as
$R_m\left[n\right]={\log }_2\left(1+\frac{G_u{p}_u\left[n\right]h_m\left[n\right]+G_b{p}_b\left[n\right]f_m}{G_s(\theta )P_s{g}_m+{\sigma }_m^2}\right)$
where $g_m$ represents the channel gain from SAT to SU and can be calculated as $g=|\bar{g}/D{|}^2$ , where
$D$ denotes the distance between SAT and ground user.
${\sigma }_m^2$ denotes the power of additive Gaussian white noise at SU,
$G_u$ and $G_b$ represent the UAV and BS directional antenna gain,
$G_s(\theta )$ denotes the SAT’s transmit antenna gain, which is determined by the angle between SAT and SU, and can be approximately represented as [34]
$G_s\left(\theta \right)=G_{tx}{\left(\frac{J_1\left(u\right)}{2u}+36\frac{J_3\left(u\right)}{u^3}\right)}^2$

Secure Transmission in Cognitive Satellite Terrestrial Networks

K. An, M. Lin, J. Ouyang and W. -P. Zhu, “Secure Transmission in Cognitive Satellite Terrestrial Networks,” in IEEE Journal on Selected Areas in Communications, vol. 34, no. 11, pp. 3025-3037, Nov. 2016, doi: 10.1109/JSAC.2016.2615261.

the SAT, PU, SU and EVE are all equipped with a single antenna
the terrestrial BS has N antennas

the overall channel for satellite link can be modeled as
$g=\sqrt{b\left(\phi \right)}\bar{g}$

Given a receiver’s position within the satellite spot beam coverage area, the beam gain factor can be approximated as [1], [2]
$b\left(\phi \right)={\left(\frac{J_1\left(u\right)}{2u}+36\frac{J_3\left(u\right)}{u^3}\right)}^2$
$u=2.07123\frac{\sin \phi }{\sin {\phi }_{3dB}}$
$\phi$ is the angle between the location of the corresponding receiver and the beam center with respect to the satellite, and ${\phi }_{3dB}$ is the 3-dB angle.

the channel fading coefficient $\bar{g}$ is described as
$\bar{g}=A\exp \left(j\psi \right)+Z\exp \left(j\zeta \right)$

$A$ is the amplitudes of the scattering components, which are independent stationary random processes following Rayleigh distributions.
$p_A(a)=\frac{a}{b_0}\exp \left(\frac{-a^2}{2b_0}\right),a\ge 0$
$\psi$ is a stationary random phase vector with elements uniformly distributed over $[0,2\pi )$ ,

$Z$ is the amplitudes of the LOS components, which are independent stationary random processes following Nakagami-$m$ distributions.
$p_Z(z)=\frac{2m^m}{\Gamma (m){\Omega }^m}{z}^{2m-1}\exp \left(\frac{-m{z}^2}{\Omega }\right),z\ge 0$
$\zeta$ the deterministic phase vector of the LOS component.

The SR fading distribution
$\bar{g}\sim \mathrm{SR}\left(\Omega ,b,m\right)$
$\Omega$ being the average power of the LoS component,
$2b$ the average power of the multipath component, and
$m$ the Nakagami-$m$ parameter corresponding to the fading severity.

仿真参数
Omega=0.835
b=0.126
m=10
${\phi }_{3\mathrm{dB}}=0.4/180\ast \pi$

$X=\bar{\gamma }b\left(\phi \right){\left|\bar{g}\right|}^2$
$f_X(x)=\frac{\alpha }{b(\phi )\bar{\gamma }}\exp \left(-\frac{\beta }{b(\phi )\bar{\gamma }}x\right){}_1{F}_1\left(m;1;\frac{\delta }{b(\phi )\bar{\gamma }}x\right)$
the probability density function (PDF)
$f_{|h{|}^2}(x)=\alpha \exp (-\beta x{)}_1{F}_1(m;1;\delta x)$
$\alpha =\frac{1}{2b}{\left(\frac{2bm}{2bm+\Omega }\right)}^m$
$\beta =\frac{1}{2b}$
$\delta =\frac{\Omega }{2b(2bm+\Omega )}$

Outage Analysis of Cognitive Hybrid Satellite-Terrestrial Networks With Hardware Impairments and Multi-Primary Users

$M$ terrestrial primary users (single antenna)
secondary satellite (single antenna) → relay ($N$ antennas)→ destination (single antenna)

对噪声归一化了？？
${\mathbf{h}}_{SJ}=F_{SJ}{\mathbf{g}}_{SJ}=\sqrt{L\frac{G_{t,SJ}{G}_{r,SJ}}{K_{B}TB}}{\mathbf{g}}_{SJ}=\frac{C}{4\pi f{d}_{SJ}}\sqrt{\frac{G_{t,SJ}{G}_{r,SJ}}{K_{B}TB}}{\mathbf{g}}_{SJ}$

${\mathbf{g}}_{SJ}$ is the channel coefficient vector, which is often assumed to undergo SR fading [6]–[8],

$F_{SJ}$ a scaling parameter including many practical effects, such as free space loss (FSL) and antenna pattern, which is given by
$F_{SJ}=\frac{C}{4\pi f{d}_{SJ}}\sqrt{\frac{G_{t,SJ}{G}_{r,SJ}}{K_{B}TB}}$ with
$C$ being the light speed,
$f$ the carrier frequency,
$d_{SJ}$ the distance between the satellite and the user,

$K_B=1.38\times 10^{-23}J/K$ the Boltzman constant,
$T$ the receiver noise temperature,
$B$ the carrier bandwidth.

$G_{r,SJ}$ denotes the receive gain,
$G_{t,SJ}$ the satellite beam gain, which can be approximately written as
$G_{t,SJ}=G_{\max }{\left(\frac{J_1(u)}{2u}+36\frac{J_3(u)}{u^3}\right)}^2$ [20],
${\mathrm{G}}_{\max }$ being the maximal beam gain
$u=2.07123\frac{\sin \phi }{\sin {\phi }_{3\mathrm{dB}}}$ ,
$\phi$ is the angle between the location of the corresponding receiver and the beam center with respect to the satellite, and ${\phi }_{3\mathrm{dB}}$ is the 3 – dB angle.

The system parameters are given as the GEO, $f=2$ GHz, ${\phi }_{3\mathrm{dB}}=0.8^{\circ }$ , $G_{\max }=48dB$ , $G_{r,SJ}=4dB$ , $B=15$ MHz, and $T=300^{\circ }\text{K}$ .
Without loss of generality, we set $\delta 2R=\delta 2D=1$ and in all plots we denote $\gamma ¯¯¯SR=\gamma ¯¯¯SP=\gamma ¯¯¯RD=\gamma ¯¯¯RP=\gamma ¯¯¯$ , $kSR=KRD=kSP=kRP=k=kSP=kRP=k$ and $mU=1$ , $bU=0.063$ , $\Omega U=0.0007$ .

The Application of Power-Domain Non-Orthogonal Multiple Access in Satellite Communication Networks

an entire link budget of satellite$\to$ User $j$ ($j\in 1,2,\dots ,M,$ ),
including propagation loss, beam gain, channel statistical property, and receive antenna gain,
can be given by
$Q_j=L_j{G}_j{G}_s({\phi }_j){\hat{h}}_j^2$

Lj : Propagation loss which is mainly caused by attenuation or free space loss.
It is a function of distance $d$ from the satellite to the User $j$ , the light speed $c$ , and the frequency $fc$ of the transmission.
Gj : The received antenna gain at the User $j$ .
$Gs(\phi j)$ : Received beam gain at the User $j$ which is related to the angle between User $j$ and the beam center with respect to the satellite $\phi j$ and the maximum beam gain at the on-board boresight $GSmax$ [10].
${\hat{h}}_j$ : The estimated channel coefficient of link satellite$\to$ User $j$ can be given by ${\hat{h}}_j=h_j+e_j$ with $hj$ and $ej$ denoting the perfect channel coefficient and the estimated channel error caused by user mobility and/or path loss, respectively [11].
Without loss of generality, we assume that the channel power gain $|h_j{|}^2$ follows a Shadowed-Rice (SR) model [12], which is flexible to fit data and evaluate performance for satellite propagation environments by characterizing a wide range of elevation angle $\theta j$ , and has been widely applied in various frequency bands such as the UHF-band, L-band, S-band, and Ka-band.

Performance Analysis of NOMA-Based Land Mobile Satellite Networks

All nodes in the proposed model are also assumed to equip with a single antenna

The satellite channel model including beam gain, fading model, and free space loss (FSL) is described in the following.
$G_j({\phi }_j)=G_j{\left(\frac{J_1\left(u_j\right)}{2u_j}+36\frac{J_3\left(u_j\right)}{u_j^3}\right)}_{}^2$
$f_{{|h_j|}^2}\left(x\right)={\alpha }_j{e}^{-{\beta }_{j}x}{}_1{F}_1\left(m_j;1;{\delta }_{j}x\right)$
$L_j={\left(\frac{c}{4\pi f_c{d}_j}\right)}_{}^2$
$y_j=\sqrt{L_j{G}_s{G}_j({\phi }_j)}{h}_{j}x+n_j$

Specifically, the channel parameters depending on the referred shadowing scenario for satellite links are given in Table I [27].
Moreover, we set the carrier frequency to be 1.6 GHz,
$\phi p=0.1\circ$ , $\phi q=0.6\circ$ , ${\phi }_{p\text{3dB}}={\phi }_{q\text{3dB}}=0.4^{\circ }$ , and
$G_p=G_q=3.5$ dBi, $Gs=24.3$ dBi,
$\xi =2$ , and $Pint=50$ W [7], [40].
The label (LS/HS) denotes the link shadowing severity of User-$p$ / User-$q$ .

On the Performance of Cognitive Satellite-Terrestrial Networks

a multi-beam satellite
BSs (equipped with multiple antennas)

the beam gain can be approximated as [3]
$G_{ii}={\mathcal{L}}_{\max }{G}_{s,i}{G}_{r,i}{\left(\frac{J_1(x)}{2x}+36\frac{J_3(x)}{x^3}\right)}^2$
where $Lmax$ is the free space loss [24]

The overall channel gain between the j th beam and i th user of the satellite can be given as
$h_{pp}^{ij}=h_{pp}^j{G_{ij}\left({\varphi }_{ij}\right)}^{1/2},i,j=1,\dots ,K.$

The received signal at i th beam user can be formulated as
$y_i=\sqrt{P_{si}}{G}_{ii}{h}_{pp}^i{x}_p^i+\sum _{j\in {\Phi }_U,j\ne i}\sqrt{P_{sj}}{G}_{ij}{h}_{pp}^i{x}_p^j+{\mathcal{I}}_{\text{BS}}+{\omega }_i$

测试代码：已知信道包络平方的概率密度函数，生成信道包络平方

close all
clear

parameter_Shadowed_Rician=[0.158,19.4,1.29; % Infrequent light shadowing
    0.063,0.739,8.97*10^(-4); % Frequent heavy shadowing
    0.251,5.21,0.278; % Overall results
    0.126,10.1,0.835];% Average shadowing
shadowing_scenario = 4;
bb=parameter_Shadowed_Rician(shadowing_scenario,1);
mm=parameter_Shadowed_Rician(shadowing_scenario,2);
omega_=parameter_Shadowed_Rician(shadowing_scenario,3);
prob_density_f=@(xx) prob_density_Shadowed_Rician_power(bb,mm,omega_,xx);

range_begin=0;
range_end=5;
hh=linspace(range_begin,range_end,50); 
ff=prob_density_f(hh);
intergra_prob_density=trapz(hh,ff);  %计算整个区间概率密度的积分
ff=ff/intergra_prob_density;         %归一化概率密度

figure 
plot(hh,ff);%根据公式画概率密度曲线
legend('prob density');
xlabel('The square of the channel envelope');
% ylabel('prob density');
hold on;
grid on

num_rand=1000; %需要随机数的个数
rand_from_prob_density = gene_rand_from_prob_density(prob_density_f,...
    range_begin,range_end,intergra_prob_density,num_rand);

% 统计检验随机数列是否符合分布
num_interval=100; %直方图的区间个数
[counts,centers]=hist(rand_from_prob_density,num_interval);    %统计不同区间出现的个数
range_interval=(range_end-range_begin)/num_interval;        %区间大小
counts=counts/num_rand/range_interval;         %根据统计结果计算概率密度
bar(centers,counts,1);   %根据统计结果画概率密度直方图
% title(['生成了',num2str(num_rand),'个随机数，并分了',num2str(num_interval),'个区间画直方图'])
title(['Generate ',num2str(num_rand),' random numbers, and divide ',num2str(num_interval),' intervals to draw histogram'])


function rand_from_prob_density = gene_rand_from_prob_density(prob_density_f,...
    range_begin,range_end,intergra_prob_density,num_rand)
% 舍选法 (Acceptance-Rejection Method)
rand_from_prob_density = zeros(num_rand,1);
nn=1;
while nn<=num_rand
    %生成[range_begin,range_end]均匀分布随机数
    xx=rand(1)*(range_end-range_begin)+range_begin;
    %生成[0,1]均匀分布随机数
    rand_prob=rand(1);  
    %如果随机数r小于密度函数值f(t)，接纳该t并加入序列a中
    if rand_prob<=prob_density_f(xx)/intergra_prob_density     
        rand_from_prob_density(nn) = xx;
        nn = nn+1;
    end
end
end

function prob_den_Shadowed_Rician_power = prob_density_Shadowed_Rician_power(bb,mm,omega_,hsquare)
alpha = (2*bb*mm/(2*bb*mm+omega_))^mm / (2*bb);
beta = 1 / (2*bb);
delta = omega_/(2*bb*(2*bb*mm+omega_));
% https://www.mathworks.com/help/symbolic/hypergeom.html#bt1nkmw-2
prob_den_Shadowed_Rician_power = alpha .* exp(-beta*hsquare) .* hypergeom(mm,1,delta.*hsquare);
end

the generalized-K model

K. P. Peppas, " Accurate closed-form approximations to generalised- $K$ sum distributions and applications in the performance analysis of equal-gain combining receivers ", IET Commun., vol. 5, no. 7, pp. 982-989, May 2011.
I. S. Ansari, S. Al-Ahmadi, F. Yilmaz, M.-S. Alouini, and H. Yanikomeroglu, “A new formula for the BER of binary modulations with dual-branch selection over generalized-K composite fading channels,” IEEE Trans. Commun., vol. 59, no. 10, pp. 2654–2658, Oct. 2011.

Energy Efficient Adaptive Transmissions in Integrated Satellite-Terrestrial Networks With SER Constraints

Y. Ruan, Y. Li, C. -X. Wang, R. Zhang and H. Zhang, “Power Allocation in Cognitive Satellite-Vehicular Networks From Energy-Spectral Efficiency Tradeoff Perspective,” in IEEE Transactions on Cognitive Communications and Networking, vol. 5, no. 2, pp. 318-329, June 2019, doi: 10.1109/TCCN.2019.2905199.

Ka频段 自由空间损耗 雨衰 波束增益

Physical Layer Security in Multibeam Satellite Systems
Ka频段

Z. Lin, M. Lin, J. -B. Wang, T. de Cola and J. Wang, “Joint Beamforming and Power Allocation for Satellite-Terrestrial Integrated Networks With Non-Orthogonal Multiple Access,” in IEEE Journal of Selected Topics in Signal Processing, vol. 13, no. 3, pp. 657-670, June 2019, doi: 10.1109/JSTSP.2019.2899731. 不太一样自由空间损耗 雨衰 波束增益 阵列控制向量
mmWave

M. Lin, C. Yin, Z. Lin, J. -B. Wang, T. de Cola and J. Ouyang, “Combined Beamforming with NOMA for Cognitive Satellite Terrestrial Networks,” ICC 2019 - 2019 IEEE International Conference on Communications (ICC), 2019, pp. 1-6, doi: 10.1109/ICC.2019.8761139.
Ka频段

Joint Beamforming for Secure Communication in Cognitive Satellite Terrestrial Networks
mmWave

the free space loss (FSL)
$C_L={\left(\frac{\lambda }{4\pi d}\right)}^2$

the rain attenuation fading vector $\mathbf{r}\in C^{N_s\times 1}$ between the satellite antenna array and the ground user
$\mathbf{r}={\xi }^{-\frac{1}{2}}{e}^{-j\mathbf{p}}$
$\mathbf{p}$ is the $N_s\times 1$ sized phase vector with its components uniformly distributed over $[0,2\pi )$ 随机
According to ITU-R Recommendation P.1853 [39], the power gain due to rain attenuation in dB, ${\xi }_{\mathrm{dB}}=20{\log }_{10}\left(\xi \right)$ , commonly follows a lognormal distribution, namely, $\ln \left({\xi }_{\mathrm{dB}}\right)\sim \mathcal{C}\mathcal{N}\left(\mu ,{\sigma }^2\right)$ , 随机
where $\mu$ and $\sigma$ , both expressed in dB, are the lognormal location and scale parameter, respectively.
$\mu =-3.125$
$\sigma =1.591$

the antenna gain from the $m$ -th on-board beam to the user
$b_m=b_{\max }{\left(\frac{J_1\left(u_m\right)}{2u_m}+36\frac{J_3\left(u_m\right)}{u_m^3}\right)}^2$
$b_{\max }$ is the maximal satellite antenna gain of the $m$ -th beam,
$u_m=2.07123\sin {\varphi }_m/\sin {\left({\varphi }_{3dB}\right)}_m$ , 怎么确定？？
$J_1\left(\cdot \right)$ and $J_3\left(\cdot \right)$ the first-kind Bessel functions of order 1 and 3, respectively.

the overall satellite channel
$\mathbf{f}=\sqrt{C_L}\mathbf{r}\odot {\mathbf{b}}^{\frac{1}{2}}$
$\mathbf{b}={\left[b_1,\cdots ,b_{N_s}\right]}^T\in C^{N_s\times 1}$ is the beam gain vector

f=20GHz
$b_{\max }=52dB$
${\phi }_{3\mathrm{dB}}=0.4/180\ast \pi$
噪声带宽=50MHz

L个散射体的阵列响应

Adaptive scheduling for millimeter wave multi-beam satellite communication systems

Secure Satellite-Terrestrial Transmission Over Incumbent Terrestrial Networks via Cooperative Beamforming
J. Du, C. Jiang, H. Zhang, X. Wang, Y. Ren and M. Debbah, “Secure Satellite-Terrestrial Transmission Over Incumbent Terrestrial Networks via Cooperative Beamforming,” in IEEE Journal on Selected Areas in Communications, vol. 36, no. 7, pp. 1367-1382, July 2018, doi: 10.1109/JSAC.2018.2824623.
mmWave

the channel vector ${\mathbf{h}}_n\in {\mathbb{C}}^{N_s\times 1}$ of FSS terminal $n$ ( n ∈ N ≜ { 1 , 2 , ⋯   , N } n \in \mathsf {\mathcal {N}} \triangleq \left \{{ 1,2, \cdots,N }\right \} n∈N≜{1,2,⋯,N} )
${\mathbf{h}}_n=\sqrt{\frac{N_s}{L}}{\sum }_{l=1}^L{\delta }_{n,l}\mathbf{\alpha }\left({\theta }_{n,l}\right),\forall n\in \mathcal{N}$
where $L$ is the number of scatters,
${\delta }_{n,l}$ and ${\theta }_{n,l}$ are the complex gain and normalized direction of the LOS path for FSS $n$ , respectively.
In addition, ${\delta }_{n,l}^2\sim \mathcal{C}\mathcal{N}\left(0,{\sigma }_0^2\right)$ is independent identically distributed (i.i.d.) complex Gaussian distribution with zero-mean, and covariance ${\sigma }_0^2=1$ , which indicates Rician factor, and 随机
${\theta }_{n,l}\sim U\left[-1,1\right]$ is i.i.d. uniformly distributed. 随机
L=2 [36]

Moreover, when a uniform linear array (ULA) is adopted, the normalized array response $\alpha (\theta )$ is given by
$\mathbf{\alpha }\left(\theta \right)=\frac{1}{\sqrt{N_s}}{\left[1,e^{-j\frac{2\pi }{\lambda }d\sin \left(\phi \right)},\cdots ,e^{-j\frac{2\pi }{\lambda }\left(N_s-1\right)d\sin \left(\phi \right)}\right]}^T$
normalized direction $\theta n$ is related to the physical azimuth angle of departure (AoD) of $\phi \in [-\pi /2,\pi /2]$ as $\theta =(2d/\lambda )sin(\phi )$ , 怎么确定？？
where $d$ is the antenna spacing (i.e., the distance between the two adjacent antennas), and
$\lambda$ is the carrier wavelength.
In this work, we assume the critically sampled environment, i.e, $d/\lambda =0.5$ , considering that the normalized AoD is the sine function of the actual AoD.