LU矩阵分解法解线性方程组（C++实现）

最新推荐文章于 2024-01-16 00:01:41 发布

原创

最新推荐文章于 2024-01-16 00:01:41 发布 · 1.7k 阅读

26 ·

CC 4.0 BY-SA版权

文章标签：

#线性代数 #矩阵 #算法

本文介绍了一种使用LU矩阵分解法求解线性方程组的方法，并提供了完整的C++代码实现。该方法首先将系数矩阵分解为下三角矩阵L和上三角矩阵U，然后通过前向和后向代换求解未知数。

LU矩阵分解法解线性方程组（C++实现）

#include<cstdio>
#include<iostream>
using namespace std;
void fenjie(double a[],double b[],int n){
   
   
	double aa[n][n];
	for(int i=0;i<n;i++){
   
   
		for(int j=0;j<n;j++){
   
   
			aa[i][j]=a[n*i+j];
		}
	}
   double l[n][n];
   double u[n][n];
   for(int i=1;i<=n;i++){
   
   
   	u[0][i-1]=aa[0][i-1];
   	l[i-1][0]=aa[0][i-1]/u[0][