数值分析(任玉杰)第二章习题.docVIP

第二章线性方程组的解法 习题3.3 高斯消元法主程序 function [RA,RB,n,X]=gaus(A,b) B=[A b]; n=length(b);RA=rank(A); RB=rank(B);zhica=RB-RA; if zhica0, disp(‘请注意：因为RA?=B，所以此方程组无解.’) return end ifRA=RB ifRA==n dispf请注意：因为RA=RB=n，所以此方程组有唯一解 X=zeros(n, 1 );C=zeros(l，n+1); for p=l:n-l for k=p+l:n m=B(k，p)/B(p，p); B(k，p:n+1 )= B(k，p:n+1 )-m*B(p，p:n+1); end end b=B(l :n，n+ l);A=B(l:n,l:n);X(n)=b(n)/A(n，n); for q=n-l:-l:l X(q)=(b(q)-sum(A(q，q+l :n)*X(q+l :n)))/A(q，q); end else disp(‘请注意：因为RA=RB0，所以此方程组有无穷多解 end end 1. (1) A=[4 8 4 0;l 5 4-3;l 4 7 2;1 3 1 -7];b=[8;-4;10;-4]; [RA，RB，n，X]=gaus(A，b) 请注意：因为RA=RB=n,所以此方程组有唯一解. RA = 4 RB = 4 n = 4 X = 8.399999999999999 -4.600000000000000 2.800000000000000 0.200000000000000 (2) A=[l 1 2 3;1 3 6 1;3-1 -3 15;1 -5-10 12];b=[l;3;3;5]; [RA，RB，n，X]=gaus(A，b) 请注意：因为RA=RB=n,所以此方程组有唯一解. RA = RB = x = -13.333333333333334 -1.000000000000000 2.666666666666667 3.333333333333334 习题3.4 1.主程序 function hl=zhjLU(A) [n n]=size(A);RA=rank(A); if RA?=n disp(请注意：因为A的n阶行列式hl等于零，所以A不能进行LU分解.A 的秩 RA 如下：*),RA,hl=det(A); return end if RA==n for p=l:n h(p)=det(A(l:p，l:p)); end hl=h(l:n); for i=l:n if h(l,i)==0 dispf请注意：因为A的r阶主子式等于零，所以A不能进行LU 分解.A的秩RA和各阶主子式值好依次如下：’)，hl;RA return end end if h（l，i）~=O disp（*请注意：因为A的各阶主子式都不等于零，所以A能进行LU分 解.A的秩RA和各阶顺序主子式值hl依次如下： for j=l:n U（l，j）=A（l，j）; end for k=2:n for i=2:n for j=2:n L（l，l）=l;L（i，i）=l; ifi〉j L（ 1，1 ）= 1 ;L（2,1 ）=A（2,1）/U（ 1，1 ）;L（i，1 ）=A（i，1 ）/U（ 1，1）; L（i，k）=（A（i，k）-L（i，1 :k-1 ）*U（ 1 :k-1，k））/U（k，k）; else U（kj）=A（kj）-L（k,l:k- l）*U（l:k-lj）; end end end end hl;RA,U,L end end （1）. A=[2 4 -6;1 5 3;1 3 2];hl=zhjLU（A） 请注意：因为A的各阶主子式都不等于零，所以A能进行LU分解.A的秩RA 和各阶顺序主子式值hl依次如下： RA 二 2.000000000000000 4.000000000000000 -6.000000000000000 0 5.000000000000000 6.000000000000000 0 0 3.800000000000000 1.000000000000000 0 0 0.500000000000000 1.000000000000000 0 0.500000000000000 0.200000000000000 1.000000000000000 hl = 2 6 18 . A=[l 1 6;-l 2 9;1 -2 3];hl=zhjLU(A) 请注意：因为A的各阶主子式都不等于零，所以A能进行LU分解.A的秩RA 和各阶顺序主子式值hl依次如下： RA = 3 1.00000000000000000L = 1.00000000