瑞利噪声python

瑞利噪声通常用于模拟无线信道信号传输的衰落。在Python中，我们可以使用NumPy中的random模块来生成瑞利分布的随机变量序列。下面是一个简单的例子展示了如何使用Python生成瑞利噪声：

import numpy as np
import matplotlib.pyplot as plt

# 设置噪声参数
N = 20000  # 采样点数
fs = 2000  # 采样频率
fc = 20    # 信号频率
A = 1     # 信号幅度
sigma = 0.2 # 噪声方差

# 生成瑞利噪声
rayleigh_noise = np.random.rayleigh(scale=1, size=N)*1e-3

# 生成正弦信号
t = np.arange(N) / float(fs)
signal = A * np.sin(2.0 * np.pi * fc * t)

# 加入噪声
noisy_signal = signal + rayleigh_noise

# 绘制图像
plt.subplot(3, 1, 1)
plt.plot(t, signal)
plt.title('Signal')
plt.xlabel('Time(s)')
plt.ylabel('Amplitude')

plt.subplot(3, 1, 2)
plt.plot(t, rayleigh_noise)
plt.title('Rayleigh Noise')
plt.xlabel('Time(s)')
plt.ylabel('Amplitude')

plt.subplot(3, 1, 3)
plt.plot(t, noisy_signal)
plt.title('Noisy Signal')
plt.xlabel('Time(s)')
plt.ylabel('Amplitude')

plt.show()

上述代码中，我们使用NumPy中的random模块生成了一个瑞利分布的随机变量序列(rayleigh_noise)，然后将其加到一个正弦信号(signal)上，得到一个带有瑞利噪声的信号(noisy_signal)。最后，我们使用Matplotlib库绘制了信号和噪声的波形图以及带噪声的信号波形图。