python添加高斯白噪声及其原理
生活随笔
收集整理的這篇文章主要介紹了
python添加高斯白噪声及其原理
小編覺得挺不錯的,現在分享給大家,幫大家做個參考.
python添加高斯白噪聲及其原理
python添加高斯白噪聲及其原理
1 原理公式
Ps、PnPs、PnPs、Pn分別表示信號和噪聲的有效功率,信噪比 (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添加高斯白噪声及其原理的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: 音频功放的种类和基本原理
- 下一篇: python脚本使用hive正则表达式斜