如何在python中拟合三个高斯峰值？

import numpy as np
import matplotlib.pyplot as plt
from scipy.optimize import curve_fit

def gaussian(x, a1, b1, c1, a2, b2, c2, a3, b3, c3):
    return (a1 * np.exp(-(x - b1)**2 / (2 * c1**2)) +
            a2 * np.exp(-(x - b2)**2 / (2 * c2**2)) +
            a3 * np.exp(-(x - b3)**2 / (2 * c3**2)))

x = np.linspace(0, 10, 100)
y = gaussian(x, 1, 5, 1, 2, 7, 0.5, 0.5, 9, 0.3) + np.random.normal(0, 0.2, 100)

popt, pcov = curve_fit(gaussian, x, y)

plt.plot(x, y, 'b-', label='data')
plt.plot(x, gaussian(x, *popt), 'r-', label='fit')
plt.legend()
plt.show()