数学建模社区-数学中国
标题:
【原创】LU分解求解线性方程组Matlab代码
[打印本页]
作者:
建不了的模。
时间:
2015-1-9 10:01
标题:
【原创】LU分解求解线性方程组Matlab代码
function [x,L,U]=lusolve(A,b)
% LU分解求解线性方程组Ax=b
%
% 参数说明
% A,b：方阵A满足前n-1阶顺序主子式不为零，b为非零列向量
% x：方程组根
% L,U：L为单位下三三角阵，U的上三角阵(Doolittle分解)，满足A=LU，
%
% LU分解的说明：
% 方阵A满足前n-1阶顺序主子式不为零
% 则A经过行变换的得到一个上三角阵U，而初等行变换相当于左乘一个初等矩阵L，故A可分解为A=LU，
% 其中L或U为下或上三角矩阵，可以证明LU分解不唯一
%
% 如果L为单位下三角阵，U为上三角阵，那么称为Doolittle分解
% 如果L为下三角阵，U为单位上三角阵，那么称为Crout分解
%
% 可以证明Doolittle和Crout分解满足唯一性
%
% LU分解程组的思路：
% Ax=(LU)x=L(Ux)=Ly=b，其中x,y,b都是列向量
% 先求解方程Ly=b，计算出y，然后通过Ux=y求解x
%
% 对A进行Doolittle分解
[L,U]=Doolittle(A);
% 调用trimatrix()求解Ly=b
y=trimatrix(L,b,'L');
% 调用trimatrix求解Ux=y
x=trimatrix(U,y,'U');
作者:
WSHXPY
时间:
2015-1-9 15:25
学习一下……
作者:
WSHXPY
时间:
2015-1-9 15:25
学习一下……
作者:
WSHXPY
时间:
2015-1-9 15:26
学习一下……
作者:
WSHXPY
时间:
2015-1-9 15:26
学习一下……
作者:
WSHXPY
时间:
2015-1-9 15:26
学习一下……
作者:
WSHXPY
时间:
2015-1-9 15:26
学习一下……
作者:
WSHXPY
时间:
2015-1-9 15:26
学习一下……
作者:
WSHXPY
时间:
2015-1-9 15:26
学习一下……
作者:
WSHXPY
时间:
2015-1-9 15:26
学习一下……
作者:
WSHXPY
时间:
2015-1-9 15:27
学习一下……
作者:
WSHXPY
时间:
2015-1-9 15:27
学习一下……
作者:
WSHXPY
时间:
2015-1-9 15:27
学习一下……
作者:
WSHXPY
时间:
2015-1-9 15:27
学习一下……
作者:
WSHXPY
时间:
2015-1-9 15:27
学习一下……
作者:
w785485068
时间:
2015-1-9 23:10
顶一下。。。。。。
作者:
meciny
时间:
2015-5-6 19:58
终于可以回复了。。。
作者:
huang_ray
时间:
2015-10-5 12:34
学习一下
作者:
imagination13
时间:
2015-10-24 14:19
学习一下，这个是紧凑格式吧
作者:
zsedc
时间:
2015-11-6 22:53
学习一下……
作者:
627743701
时间:
2015-12-7 17:17
RE: 【原创】LU分解求解线性方程组Matlab代码
作者:
，——！
时间:
2016-4-6 22:37
kankan~~~~~~~~~~~~~~~~
作者:
，——！
时间:
2016-4-6 22:37
kankan~~~~~~~~~~~~~~~~
作者:
神水那月
时间:
2016-5-28 22:20
用下试试吧w
作者:
FatSheep-SYSU
时间:
2016-10-5 18:01
来学习一下
作者:
SiriusT
时间:
2017-3-19 17:34
学习学习
作者:
Danru
时间:
2017-3-21 20:00
不错，顶一个
作者:
niuzt
时间:
2017-6-25 20:23
学习一下。。。
作者:
niuzt
时间:
2017-6-25 20:24
谢谢分享，学习一下
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5