最小二乘法拟合任意次曲线(C#).docVIP

最小二乘法拟合任意次曲线（C#） 说明： 代码较为简洁没有过多的说明，如有不明白之处可查阅相关最小二乘法计算步骤资料和求解线性方程组的资料。另外该方法只能实现二元N次拟合，多元方程不适用。 以下是最小二乘法类的实现： public class MatrixEquation { private double[,] gaussMatrix; private int coe; public MatrixEquation() { } public MatrixEquation(double[] arrX, double[] arrY, int n) { coe = n; gaussMatrix= GetGauss(GetXPowSum(arrX, n), GetXPowYSum(arrX, arrY, n), n); } public double[,] GetGaussMatrix() { return gaussMatrix; } public double[] GetResult() { return ComputeGauss(gaussMatrix, coe); } /// summary /// 计算获取x散点的幂次和数组 /// /summary /// param name=arrXx散点序列/param /// param name=n函数拟合次数/param /// returns/returns protected double[] GetXPowSum(double[] arrX, int n) { int m = arrX.Length;//X散点的个数 double[] xPow = new double[2 * n + 1]; //存储X散点的幂次值 for (int i = 0; i xPow.Length; i++) { if (i == 0) { xPow[i] = m; } else { //计算x的i次方和 double max = 0; for (int j = 0; j m; j++) { if (arrX[j] == 0) max = max + 1; else max = max + Math.Pow(arrX[j], i); } xPow[i] =Math.Round( max,4); } } return xPow; } /// summary /// 计算获取xy的幂次和序列 /// /summary /// param name=arrXx散点序列/param /// param name=arrYy散点序列/param /// param name=n拟合曲线次数/param /// returns/returns protected double[] GetXPowYSum(double[] arrX, double[] arrY, int n) { int m = arrX.Length;//X散点的个数