高斯滤波器实验

高斯平滑滤波，可以有效的去除数据中的高斯噪声。可以用简单的实验来验证。
1. 首先生成一组正弦离散数据
2. 生成一组噪声
3. 把噪声叠加到数据中，生成噪声数据
4. 在噪声数据上用高斯滤波器

 

python 代码如下：

import matplotlib.pyplot as plt
import numpy as np
import scipy.ndimage as ndimage

x=np.linspace(0,10,50)
y=5*np.sin(2*np.pi*x/2.0)
n=np.random.normal(0,1,50)
r=y + n
r1=ndimage.gaussian_filter1d(r,1)

plt.figure()
plt.subplot(211)
plt.plot(x,y,'bo')
plt.plot(x,r1,'ro')
plt.subplot(212)
plt.plot(x,r,'bo')
plt.plot(x,r1,'ro')