资源简介
对矩阵进行分解,包括LU分解,QR 分解,Householder分解,Givens分解,代码采用matlab编写,包括清楚的注释,下载后即可直接运行。
代码片段和文件信息
%% LU分解
% 函数功能:输入一个方阵A,对A进行LU分解,得到上三角矩阵U,下三角矩阵L,伴随矩阵置换矩阵。
% 输入:方阵A
% 输出:下三角矩阵L,上三角矩阵U,置换矩阵P,L*U结果。
function [LUP] = LU_fc(A)
a=A;
N = size(A1);
M = size(A2);
if((N-M) ~= 0)
error(‘不是方阵,无法进行LU分解!‘);
end
if(det(A) == 0)
error(‘是奇异矩阵,无法进行LU分解!‘);
end
L = zeros(N);
P = eye(N);
%% LU分解:PA=LU
for i=1:N-1
for j=i+1:N
if (a(ii)==0)
disp(‘出现0主元,需要对矩阵进行置换变换‘);
for k = N:-1:i+1
if (a(ki)~=0)
L([i k]:)=L([k i]:);
a([i k]:)=a([k i]:);
P([i k]:)=P([k i]:);
break;
end
end
end
x=a(ji)/a(ii);
L(ji) = x;
a(j:)=a(j:)-x*a(i:);
end
end
U=a;
L=L+eye(N);
属性 大小 日期 时间 名称
----------- --------- ---------- ----- ----
文件 962 2017-11-23 18:41 矩阵分解完整版代码\LU_fc.m
文件 1854 2017-11-23 20:01 矩阵分解完整版代码\Matrix_Factorization.m
文件 730 2017-11-23 19:56 矩阵分解完整版代码\QR_Givens.m
文件 727 2017-11-23 18:59 矩阵分解完整版代码\QR_Householder.m
文件 636 2017-11-23 19:58 矩阵分解完整版代码\QR_Schmidt.m
文件 428 2017-11-23 12:57 矩阵分解完整版代码\大作业要求.txt
文件 681 2017-11-23 20:02 矩阵分解完整版代码\矩阵分解代码说明.txt
目录 0 2018-01-10 00:27 矩阵分解完整版代码
----------- --------- ---------- ----- ----
6018 8
评论
共有 条评论