(1.4)--2023数值实验——练习3.pptVIP

三角分解（Doolittle分解）算法clear,clcn=5;A=randn(n)+n*eye(n);B=A;fork=1:n-1fori=k+1:nA(i,k)=A(i,k)/A(k,k);forj=k+1:nA(i,j)=A(i,j)-A(i,k)*A(k,j);endendendL=tril(A,-1)+eye(n),U=triu(A),E=B-L*U第三章：线性方程组的数值求解

三角分解（Doolittle分解）算法（简化）clear,clcn=5;A=randn(n)+0.5*eye(n);B=A;fork=1:n-1A(k+1:n,k)=A(k+1:n,k)/A(k,k);A(k+1:n,k+1:n)=A(k+1:n,k+1:n)-A(k+1:n,k)*A(k,k+1:n);endL=tril(A,-1)+eye(n),U=triu(A),E=B-L*U

列主元三角分解（Doolittle分解）算法clear,clcn=5;A=randn(n)+0.05*eye(n);B=A;P=eye(n);fork=1:n-1[~,q]=max(abs(A(k:n,k)));q=q+k-1;A([k,q],:)=A([q,k],:);P([k,q],:)=P([q,k],:);A(k+1:n,k)=A(k+1:n,k)/A(k,k);A(k+1:n,k+1:n)=A(k+1:n,k+1:n)-A(k+1:n,k)*A(k,k+1:n);endL=tril(A,-1)+eye(n),U=triu(A),P,E=P*B-L*U

高斯消去法解线性方程组clear,clcn=5;A=randn(n)+0.5*eye(n);B=A;b=randn(n,1);A1=[A,b];fork=1:n-1A1(k+1:n,k)=A1(k+1:n,k)/A1(k,k);A1(k+1:n,k+1:n+1)=A1(k+1:n,k+1:n+1)-A1(k+1:n,k)*A1(k,k+1:n+1);endx=zeros(n,1);fori=n:-1:2x(i)=A1(i,n+1)/A1(i,i);A1(1:i-1,n+1)=A1(1:i-1,n+1)-x(i)*A1(1:i-1,i);endx(1)=A1(1,n+1)/A1(1,1),err=b-B*x

三角分解法解线性方程组clear,clcn=5;A=randn(n)+0.5*eye(n);B=A;b=randn(n,1);c=b;fork=1:n-1A(k+1:n,k)=A(k+1:n,k)/A(k,k);A(k+1:n,k+1:n)=A(k+1:n,k+1:n)-A(k+1:n,k)*A(k,k+1:n);endL=tril(A,-1)+eye(n);U=triu(A);y=zeros(n,1);x=y;fori=1:n-1y(i)=b(i)/L(i,i);b(i+1:n)=b(i+1:n)-y(i)*L(i+1:n,i);endy(n)=b(n)/L(n,n);fori=n:-1:2x(i)=y(i)/U(i,i);y(1:i-1)=y(1:i-1)-x(i)*U(1:i-1,i);endx(1)=y(1)/U(1,1),err=c-B*x

列主元高斯消去法解线性方程组clear,clcn=5;A=randn(n)+0.5*eye(n);B=A;b=randn(n,1);A1=[A,b];fork=1:n-1[~,q]=max(abs(A1(k:n,k)));q=q+k-1;A1([k,q],:)=A1([q,k],:);A1(k+1:n,k)=A1(k+1:n,k)/A1(k,k);A1(k+1:n,k+1:n+1)=A1(k+1:n,k+1:n+1)-A1(k+1:n,k)*A1(k,k+1:n+1);endx=zeros(n,1);fori=n:-1:2x(i)=A1(i,n+1)/A1(i,i);A1(1:i-1,n+1)=A1(1:i-1,n+1)-x(i)*A1(1:i-1,i);endx(1)=A1(1,n+1)/A1(1,1),err=b-B*x

列主元三角分解法解线性方程组clear,clcn=5;A=randn(n)+0.5*eye(n);B=A;b=randn(n,1);c=b;P=eye(n);fork=1:n-1[~,q]=max(abs(A(k:n,k)));q=