/// <summary> /// 对一组点通过最小二乘法进行线性回归 /// </summary> /// <param name="parray"></param> public static void linearregression(point[] parray) { //点数不能小于2 if (parray.length < 2) { console.writeline( "点的数量小于2,无法进行线性回归" ); return ; } //求出横纵坐标的平均值 double averagex = 0, averagey = 0; foreach (point p in parray) { averagex += p.x; averagey += p.y; } averagex /= parray.length; averagey /= parray.length; //经验回归系数的分子与分母 double numerator = 0; double denominator = 0; foreach (point p in parray) { numerator += (p.x - averagex) * (p.y - averagey); denominator += (p.x - averagex) * (p.x - averagex); } //回归系数b(regression coefficient) double rcb = numerator / denominator; //回归系数a double rca = averagey - rcb * averagex; console.writeline( "回归系数a: " + rca.tostring( "0.0000" )); console.writeline( "回归系数b: " + rcb.tostring( "0.0000" )); console.writeline( string .format( "方程为: y = {0} + {1} * x" , rca.tostring( "0.0000" ), rcb.tostring( "0.0000" ))); //剩余平方和与回归平方和 double residualss = 0; //(residual sum of squares) double regressionss = 0; //(regression sum of squares) foreach (point p in parray) { residualss += (p.y - rca - rcb * p.x) * (p.y - rca - rcb * p.x); regressionss += (rca + rcb * p.x - averagey) * (rca + rcb * p.x - averagey); } console.writeline( "剩余平方和: " + residualss.tostring( "0.0000" )); console.writeline( "回归平方和: " + regressionss.tostring( "0.0000" )); } |