最小二乘法拟合一条直线(C语言代码) #include
#define N 9 //N为要拟合的数据的个数
float X[9] = {1,1.477,1.778,2,2.176,2.301,2.398,2.477,2.638};
float Y[9] = {7.55,8.02,8.41,8.64,8.75,8.79,8.81,8.83,8.87};
float K=0; //拟合直线的斜率
float R=0; //拟合直线的截距
float x_sum_average=0; //数组 X[N] 个元素求和 并求平均值
float y_sum_average=0; //数组 Y[N] 个元素求和 并求平均值
float x_square_sum=0; //数组 X[N] 个个元素的平均值
float x_multiply_y=0; //数组 X[N]和Y[N]对应元素的乘机
float Squre_sum(float c[N]) ;
float Sum_Average(float d[N]);
float X_Y_By(float m[N],float n[N]);
float Squre_sum(float c[N]);
void Line_Fit(void);
void Line_Fit(void)
{
x_sum_average= Sum_Average(X);
y_sum_average= Sum_Average(Y);
x_square_sum = Squre_sum(X);
x_multiply_y = X_Y_By(X,Y);
K = ( x_multiply_y - N * x_sum_average * y_sum_average)/( x_square_sum - N * x_sum_average*x_sum_average );
R = y_sum_average - K * x_sum_average;
printf("K = %f\n",K);
printf("R = %f\n",R);
}
float Sum_Average(float d[N])
{
unsigned int i=0;
float z=0;