File:Fading corraza vatalaro.png

import numpy as np
import matplotlib.pyplot as plt

def K_alfa(alfa):
    K0 = 2.731
    K1 = -1.074e-1
    K2 = 2.2774e-3
    K = K0 + K1*alfa + K2*(alfa**2)
    return K

def mu_alfa(alfa):
    mu0 = -2.331
    mu1 = 1.142e-1
    mu2 = -1.939e-3
    mu3 = 1.094e-5
    mu = mu0 + mu1*alfa + mu2*(alfa**2) + mu3*(alfa**3)
    return mu

def sigma_alfa(alfa):
    sigma0 = 4.5
    sigma1 = -0.05
    sigma = sigma0 + sigma1*alfa
    return sigma

alfa = np.arange(20, 80) # elevation angles
K = K_alfa(alfa) # Rician factor
mu = mu_alfa(alfa) # lognormal mean
sigma = sigma_alfa(alfa) # lognormal variance

fig, axs = plt.subplots(nrows=1, ncols=2,\
                constrained_layout=True, figsize=(10, 5), dpi=150)

axs[0].plot(alfa, K, label='K')
axs[0].plot(alfa, sigma, label='$\sigma$')
axs[0].set_xlabel('elevation angles (degrees)')
axs[0].grid()
axs[0].legend()


axs[1].plot(alfa, mu, label='$\mu$', color='r')
axs[1].set_xlabel('elevation angles (degrees)')
axs[1].grid()
axs[1].legend()

plt.show()