%%?LU分解
%?函數(shù)功能：輸入一個方陣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（‘出現(xiàn)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??矩陣分解完整版代碼\大作業(yè)要求.txt

?????文件????????681??2017-11-23?20:02??矩陣分解完整版代碼\矩陣分解代碼說明.txt

?????目錄??????????0??2018-01-10?00:27??矩陣分解完整版代碼

-----------?---------??----------?-----??----

?????????????????6018????????????????????8