紧凑存储的杜利特尔分解法Doolittle(LU分解法)_解线性方程组的直接解法

最新推荐文章于 2021-12-04 14:34:45 发布

李典金

最新推荐文章于 2021-12-04 14:34:45 发布

阅读量6k

点赞数 1

CC 4.0 BY-SA版权

分类专栏：计算方法实验文章标签： Doolittle 杜利特尔分解法

本文链接：https://blog.youkuaiyun.com/sunshineacm/article/details/78838721

本文介绍了如何使用紧凑存储的杜利特尔分解法（Doolittle LU 分解）来直接求解线性方程组。这种方法属于计算方法实验的一部分。

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

紧凑存储的杜利特尔分解法Doolittle(LU分解法)_解线性方程组的直接解法

标签：计算方法实验

/*
紧凑存储的杜利特尔分解法Doolittle：
    如果初始矩阵不要求保留的话，可以紧凑存储。
    因为每个a[i][j]用来计算u[i][j]或l[i][j]之后不需要了，u[i][j]或l[i][j]算出后可存入a[i][j]所占的单元；
    y[i]可存入b[i]所占的单元，最后x[i]再取代y[i]存入最初b[i]所占的单元。
    注意用 '////////' 标记部分。
*/
#include <stdio.h>
#include <math.h>

const int maxn = 15;

int main(){
    double a[maxn][maxn], b[maxn];
    int i, j, k, r, n, sum, temp;

    freopen("lu.txt", "r", stdin);
    scanf(