java直方图算法_使用pymc与MCMC拟合两个正态分布（直方图）？

我正在尝试使用CCD上的光谱仪检测线条轮廓 . 为了便于考虑，我已经包含了一个演示，如果解决了，它与我实际想要解决的演示非常相似 .

import pymc as mc

import numpy as np

import pylab as pl

def GaussFunc(x, amplitude, centroid, sigma):

return amplitude * np.exp(-0.5 * ((x - centroid) / sigma)**2)

wavelength = np.arange(5000, 5050, 0.02)

# Profile 1

centroid_one = 5025.0

sigma_one = 2.2

height_one = 0.8

profile1 = GaussFunc(wavelength, height_one, centroid_one, sigma_one, )

# Profile 2

centroid_two = 5027.0

sigma_two = 1.2

height_two = 0.5

profile2 = GaussFunc(wavelength, height_two, centroid_two, sigma_two, )

# Measured values

noise = np.random.normal(0.0, 0.02, len(wavelength))

combined = profile1 + profile2 + noise

# If you want to plot what this looks like