杜利特尔分解Doolittle转化为克洛特分解Crout_解线性方程组的直接解法

最新推荐文章于 2021-10-08 11:06:01 发布

原创

最新推荐文章于 2021-10-08 11:06:01 发布

· 3.3k 阅读

3 ·

版权

该博客介绍了如何将杜利特尔(Doolittle)分解转换为克洛特(Crout)分解，这两种方法都是求解线性方程组的直接解法。实验内容包括数据文件和结果展示。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

杜利特尔分解Doolittle转化为克洛特分解Crout_解线性方程组的直接解法

标签：计算方法实验

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

const int maxn = 15;

int main(){
    double a[maxn][maxn], b[maxn], y[maxn], x[maxn], l[maxn][maxn], u[maxn][maxn], d[maxn];
    int i, j, k, r, n, sum;

    freopen("lu.txt", "r", stdin);
    scanf("%d", &n);
    for(int i = 1; i <= n; i++){
        for(int j = 1; j <= n; j++)  scanf("%lf", &a[i][j]);
        scanf("%lf", &b[i]);
    }

    for(i = 1; i <= n; i++)  u[1][i] = a[1][i];  //U的第一行元素
    for(i = 1; i <= n; i++)  l[i][1] = a[i][1] / u[1][1];  //L的第一列元素

    for(k = 2; k <= n; k++){
        for(j = k; j <= n; j++){  //计算U的第k行元素