用c语言实现最小二乘法系统辨识方法

最小二乘法是一种常用的系统辨识方法，可以用于估计线性系统的参数。下面是用C语言实现最小二乘法系统辨识方法的示例代码：

#include <stdio.h>
#include <stdlib.h>
#include <math.h>

#define MAX_POINTS 1000

// 用于存储数据点的结构体
struct point {
    double x;
    double y;
};

int main() {
    int n; // 数据点数量
    struct point data[MAX_POINTS]; // 存储数据点的数组

    // 读入数据点
    printf("请输入数据点的数量：");
    scanf("%d", &n);
    printf("请输入数据点，每行一个点，格式为：x y\n");
    for (int i = 0; i < n; i++) {
        scanf("%lf %lf", &data[i].x, &data[i].y);
    }

    // 计算最小二乘法的参数
    double sum_x = 0, sum_y = 0, sum_xy = 0, sum_xx = 0;
    for (int i = 0; i < n; i++) {
        sum_x += data[i].x;
        sum_y += data[i].y;
        sum_xy += data[i].x * data[i].y;
        sum_xx += data[i].x * data[i].x;
    }
    double mean_x = sum_x / n;
    double mean_y = sum_y / n;
    double a = (sum_xy - n * mean_x * mean_y) / (sum_xx - n * mean_x * mean_x);
    double b = mean_y - a * mean_x;

    // 输出结果
    printf("y = %fx + %f\n", a, b);

    return 0;
}

这段代码读入了数据点，然后根据最小二乘法的公式计算出了线性拟合的参数，最后输出结果。需要注意的是，这段代码只适用于线性拟合，如果需要进行非线性拟合，需要使用其他方法。