C#中的高斯拟合

f(x) = A * exp(-(x - mu)^2 / (2 * sigma^2))

using System;
using System.Collections.Generic;
using System.Linq;

public class GaussianFit
{
    public static (double A, double mu, double sigma) Fit(List<(double x, double y)> data)
    {
        double xMean = data.Average(d => d.x);
        double yMean = data.Average(d => d.y);

        double num = 0;
        double denom = 0;

        foreach (var (x, y) in data)
        {
            num += (x - xMean) * (y - yMean);
            denom += (x - xMean) * (x - xMean);
        }

        double mu = num / denom;
        double sigma = Math.Sqrt(data.Sum(d => (d.x - xMean) * (d.x - xMean)) / denom);
        double A = yMean - mu * sigma;

        return (A, mu, sigma);
    }
}