最小二乘法拟合曲线c语言
时间: 2023-10-01 22:10:28 浏览: 112
c语言-最小二乘法曲线拟合
最小二乘法拟合曲线是一种常见的数学方法,可以用于数据拟合、回归分析等领域。以下是使用C语言实现最小二乘法拟合曲线的示例代码:
```c
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#define MAX_N 100
double x[MAX_N]; // 存储数据点的横坐标
double y[MAX_N]; // 存储数据点的纵坐标
int main() {
int n; // 数据点个数
double sum_x = 0, sum_y = 0, sum_xx = 0, sum_xy = 0;
double a, b, r; // 拟合直线的斜率、截距和相关系数
// 读入数据点
printf("请输入数据点的个数:");
scanf("%d", &n);
printf("请输入数据点的横坐标和纵坐标,用空格分隔:\n");
for (int i = 0; i < n; i++) {
scanf("%lf %lf", &x[i], &y[i]);
sum_x += x[i];
sum_y += y[i];
sum_xx += x[i] * x[i];
sum_xy += x[i] * y[i];
}
// 计算拟合直线的斜率和截距
a = (n * sum_xy - sum_x * sum_y) / (n * sum_xx - sum_x * sum_x);
b = (sum_y - a * sum_x) / n;
// 计算相关系数
double sum_yy = 0, sum_y2 = 0;
for (int i = 0; i < n; i++) {
sum_yy += y[i] * y[i];
sum_y2 += (y[i] - b) * (y[i] - b);
}
r = (n * sum_xy - sum_x * sum_y) / sqrt((n * sum_xx - sum_x * sum_x) * (n * sum_yy - sum_y * sum_y));
// 输出结果
printf("拟合直线的方程为:y = %.2f x + %.2f\n", a, b);
printf("相关系数为:r = %.2f\n", r);
return 0;
}
```
在上面的代码中,我们首先读入数据点的横坐标和纵坐标,并计算出一些用于最小二乘法拟合的参数(如sum_x、sum_y、sum_xx、sum_xy等)。
然后,我们使用这些参数来计算拟合直线的斜率和截距,并使用它们来计算相关系数r。最后,我们将结果输出到屏幕上。
需要注意的是,上面的代码中只给出了一种简单的最小二乘法拟合曲线的实现方法,实际应用中可能需要根据具体情况进行改进和优化。
阅读全文