python高斯噪声函数,Python高斯拟合模拟高斯噪声数据

I need to interpolate data coming from an instrument using a gaussian fit. To this end I thought about using the curve_fit function from scipy.

Since I'd like to test this functionality on fake data before trying it on the instrument I wrote the following code to generate noisy gaussian data and to fit it:

from scipy.optimize import curve_fit

import numpy

import pylab

# Create a gaussian function

def gaussian(x, a, b, c):

val = a * numpy.exp(-(x - b)**2 / (2*c**2))

return val

# Generate fake data.

zMinEntry = 80.0*1E-06

zMaxEntry = 180.0*1E-06

zStepEntry = 0.2*1E-06

x = numpy.arange(zMinEntry,

zMaxEntry,

zStepEntry,

dtype = numpy.float64)

n = len(x)

meanY = zMinEntry + (zMaxEntry - zMinEntry)/2

sigm