欢迎访问 生活随笔!

生活随笔

当前位置: 首页 > 编程语言 > python >内容正文

python

python添加高斯白噪声及其原理

发布时间:2024/3/12 python 39 豆豆
生活随笔 收集整理的这篇文章主要介绍了 python添加高斯白噪声及其原理 小编觉得挺不错的,现在分享给大家,帮大家做个参考.

python添加高斯白噪声及其原理

python添加高斯白噪声及其原理

1 原理公式

Ps、PnPs、PnPsPn分别表示信号和噪声的有效功率,信噪比 (SNR) 是正常信号功率与无信号时噪声信号功率的比值,单位为 dB,信噪比应该越高越好。

下式中x(t)x(t)x(t)n(t)n(t)n(t) 分别为信号和噪声:

高斯白噪声公式如下:

2 程序

import numpy as np import pylab as plt import math def wgn(x, snr): # 输出为高斯白噪声'''程序中用hist()检查噪声是否是高斯分布,psd()检查功率谱密度是否为常数。'''snr = 10**(snr/10.0)xpower = np.sum(x**2)/len(x)npower = xpower / snrreturn np.random.randn(len(x)) * np.sqrt(npower)t = np.arange(0, 1000000) * 0.1 x = np.sin(t) snr=6 n = wgn(x, 6) xn = x+n # 增加了6dBz信噪比噪声的信号 print(10*math.log10(sum(x**2) / sum(n**2))) # 验算信噪比 plt.figure(figsize=(10,7)) plt.subplot(411) plt.title('Gauss Distribution') plt.hist(n, bins=100)plt.subplot(412) plt.psd(n) # psd()检查功率谱密度是否为常数plt.subplot(413) plt.plot(t[0:100],x[0:100]) plt.title('The Original Sin Signal')plt.subplot(414) plt.plot(t[0:100],xn[0:100]) plt.title('The Noisy Sin Signal') plt.tight_layout() plt.show()

以上就是全部内容

总结

以上是生活随笔为你收集整理的python添加高斯白噪声及其原理的全部内容,希望文章能够帮你解决所遇到的问题。

如果觉得生活随笔网站内容还不错,欢迎将生活随笔推荐给好友。