Python生成高斯白噪声white noise

原理

$random(N)\cdot \sqrt{\frac{\sum x_i^2}{N}\cdot 10^{-snr/10}}$

假设序列X长度为N。

import numy as np

def generate_white_noise(X, N, snr)
	noise = np.random.randn(N)
	snr = 10 ** (snr/10)
	power = np.mean(np.square(X))
	npower = power / snr
	noise = noise * np.sqrt(npower)