写一个QR分解的c语言代码

在C语言中，由于标准库一般不包含直接的QR分解函数，我们可以使用一些第三方数学库如`lapack`或者自编码实现。这里是一个简单的示例，使用`gsl`库中的功能进行QR分解，假设你已经安装了`gsl`并包含了相应的头文件：

#include <stdio.h>
#include <gsl/gsl_linalg.h>

int main() {
    const int N = 3; // 矩阵大小
    double a[N][N] = { /* 初始化你的3x3矩阵 */ };

    gsl_matrix_view Q;
    gsl_matrix R;

    // 使用gsl_linalg_qr_decomp进行QR分解
    if (gsl_linalg_qr_decomp(&a[N*N], &Q.matrix, &R)) {
        printf("QR分解失败\n");
        return 1;
    }

    // 输出结果
    for (size_t i = 0; i < N; ++i) {
        for (size_t j = 0; j <= i; ++j) {
            printf("%.6f ", gsl_matrix_get(Q.matrix, i, j));
        }
        printf("\n");
    }
    printf("\nUpper triangular matrix R:\n");
    for (size_t i = 0; i < N; ++i) {
        for (size_t j = 0; j <= i; ++j) {
            printf("%.6f ", gsl_matrix_get(R, i, j));
        }
        printf("\n");
    }

    return 0;
}

这个例子假设你有一个3x3的矩阵，如果想要处理更大规模的矩阵，你需要相应地调整`N`。请注意实际项目中，使用`lapack`或其他专业库通常会更为高效且准确。