bet188金宝搏

#这个python脚本用于产生瑞利信道系数
进口numpy np
进口matplotlib。pyplot作为plt

#正弦波形的一代
t = np。linspace (100, 750)
x_volts = 20 * np.sin (t / (2 * np.pi))
x_watts = x_volts * * 2
x_db = 10 * np.log10 (x_watts)

#参数仿真
v = 60 #速度(米/秒)
center_freq = 100 e6 #射频100 MHz
Fs = 2 e5 #采样率0.2 MHz
N = 1000 #总数量的正弦波
π= 3.14
fd = v * center_freq / 3 e8 #多普勒频移(最大)
打印(“多普勒频移(Max):“, fd)
t = np。不等(0,1,1 / Fs) #时间向量。(启动、停止、步骤)
x = np.zeros (len (t))
y = np.zeros (len (t))
因为我在范围(N):
α= (np.random.rand() - 0.5) * 2 *π
φ= (np.random.rand() - 0.5) * 2 *π
x = x + np.random.randn () * np。因为(2 *π* fd * t * np.cos(α)+φ)
y = y + np.random.randn () * np。罪(2 *π* fd * t * np.cos(α)+φ)

z = (1 / np.sqrt (N)) * (x + 1 j * y) #这是通道响应用于卷积传输数据或信号
z_mag = np.abs (z) #用于情节
z_mag_dB = 10 * np.log10 (z_mag) #转换为数据库

#卷积与瑞利衰落信道的正弦波形
y3 = np。卷积(z, x_volts)

#地块
图,轴= plt。次要情节(2,2)
轴(0,0).plot (x_volts)
轴(0,0)。set_title(“纯正弦波信号”)
轴[0,1].plot (z)
轴[0,1]。set_title(“瑞利信道响应”)
轴(1,0).plot (z_mag_dB)
轴(1,0)。set_title(“瑞利信道响应(dB)”)
轴[1].plot (y3)
轴[1]。set_title(“正弦波信号卷积”)
plt.tight_layout ()
plt.show ()