许可数值.docVIP

3、解：Doolittle分解 A== Crout分解 A== 4、解：= = =\(\)= 解：= =\(\)= 解：A= ||A||1=7，||A||2=9.5364，||A||=7， 25、%将矩阵A的LU分解 function LU(A) [N N]=size(A); R=rank(A); if R~=N disp(A不能进行LU分解) return end if R==N for i=1:N h(i)=det(A(1:i,1:i)); if h(1,i)==0 disp(A不能进行LU分解) return end end for j=1:N U(1,j)=A(1,j); end for i=2:N L(1,1)=1;L(i,i)=1; L(i,1)=A(i,1)/U(1,1); end for k=2:N for i=2:N for j=2:N if ij L(i,1)=A(i,1)/U(1,1); L(i,k)=(A(i,k)-L(i,k-1)*U(1:k-1,k))/U(k,k); else U(k,j)=A(k,j)-L(k,1:k-1)*U(1:k-1,j); end end end end L,U, End %改进平方根法分解矩阵A function gcholesky(A) [N N]=size(A); R=rank(A); if R~=N disp(A不能进行LU分解) return end if A~=A disp(A不能进行LU分解) return end if R==N for i=1:N h(i)=det(A(1:i,1:i)); if h(1,i)==0 disp(A不能进行LU分解) return end end for j=1:N U(1,j)=A(1,j); end for i=2:N L(1,1)=1;L(i,i)=1; L(i,1)=U(i,1)/U(1,1); end for k=2:N for i=2:N for j=2:N if ij L(i,k)=U(k,i)/U(k,k); else U(k,j)=A(k,j)-L(k,1:k-1)*U(1:k-1,j); end end end end L,U, end %追赶法 function x=chasing(A,b) n=length(b); if rank(A)~=n disp(A不能进行分解) return end for i=1:n if A(i,i)==0 disp(错误) return end end d=ones(n,1); a=ones(n-1,1); c=ones(n-1); for i=1:n-1 a(i,1)=A(i+1,i); c(i,1)=A(i,i+1); d(i,1)=A(i,i); end d(n,1)=A(n,n);z for i=2:n d(i,1)=d(i,1)-(a(i-1,1)/d(i-1,1))*c(i-1,1); b(i,1)=d(i,1)-(a(i-1,1)/d(i-1,1))*b(i-1,1); end x(n,1)=b(n,1)/d(n,1); for i=n-1:-1:1 x(i,1)=(b(i,1)-c(i,1)*x(i+1,1))/d(i,1); end x=x End  取节点，、，； 则 （2）取节点，、，、，； 则 取节点，，，， 则 =1/3 =0.05174 =0.5001 7、 则=3.128 8、