数值分析(第五版)计算实习题第五章作业(共12页).docx

数值分析第五章第一题：LU分解法：建立m文件function h1=zhijieLU(A,b) %h1各阶主子式的行列式值 n n=size(A);RA=rank(A);if RA=n disp(请注意：因为A的n阶行列式h1等于零，所以A不能进行LU分解。A的秩RA如下：) RA,h1=det(A); returnendif RA=n for p=1:n h(p)=det(A(1:p,1:p); end h1=h(1:n); for i=1:n if h(1,i)=0 disp(请注意：因为A的r阶主子式等于零，所以A不能进行LU分解。A的秩RA和各阶顺序主子式h1依次如下：) h1;RA return end end if h(1,i)=0 disp(请注意：因为A的r阶主子式都不等于零，