LU分解法求解线性方程组

最新推荐文章于 2025-05-14 20:33:23 发布

swordboy_fire

最新推荐文章于 2025-05-14 20:33:23 发布

阅读量1w

点赞数 6

CC 4.0 BY-SA版权

分类专栏：计算方法及程序实现文章标签： LU分解法

本文链接：https://blog.youkuaiyun.com/swordboy_fire/article/details/83243962

计算方法及程序实现专栏收录该内容

9 篇文章

订阅专栏

本文介绍了一种求解线性方程组的算法——LU分解法，该方法通过将系数矩阵分解为上三角矩阵U和下三角矩阵L来简化计算过程。文章详细阐述了LU分解的计算步骤，并提供了完整的C++代码实现。

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

LU分解法是求解线性方程组的一种算法。

先将系数矩阵分别转换成上三角矩阵U和下三角矩阵L，其中U(k,j)=a(k,j) - ∑L(k,m)*U(m,j) m<k，L(i,k)=(a(i,k) - ∑L(i,m)*U(m,j))/(u(k,k)) m<k;

求解Ly=b，b为因变量矩阵；

求解Ux=y；

c++代码如下：

#include<stdio.h>

using namespace std;

#define  n  3    //矩阵阶数
int main()
{
	int i, j, k, m;
	float t, a[n][n], u[n][n], l[n][n], b[n], x[n], y[n];
	for (i = 0; i<n; i++)
	for (j = 0; j<n; j++)
	{
		scanf("%f", &t);
		a[i][j] = t;
	}
	for (i = 0; i<n; i++)
		scanf("%f", &b[i]);
	for (k = 0; k<n; k++)
	{
		for (j = k; j<n; j++)
		{
			t = 0;
			for (m = 0; m <= k - 1; m++)
				t = t + l[k][m] * u[m][j];
			u[k][j] = a[k][j] - t;
		}
		for (i = k + 1; i<n; i++)
		{
			t = 0;
			for (m = 0; m <= k - 1; m++)
				t = t + l[i][m] * u[m][k];
			l[i][k] = (a[i][k] - t) / u[k][k];
		}
	}
	for (i = 0; i<n; i++)
	{
		t = 0;
		for (j = 0; j <= i - 1; j++)
			t = t + l[i][j] * y[j];
		y[i] = b[i] - t;
	}
	for (i = n - 1; i >= 0; i--)
	{
		t = 0;
		for (k = i + 1; k<n; k++)
			t = t + u[i][k] * x[k];
		x[i] = (y[i] - t) / u[i][i];
	}
	for (i = 0; i<n; i++)
		printf("%f\n", x[i]);
	return 0;
}