浅谈对矩阵理论数值分析及matab应用的理解 7.docVIP

附录1 Ax=b 利用高斯消元法 clear format rat A=randn(5) b=randn(5,1) [m,n]=size(A); for i=1:(m-1) numb=int2str(i); for j=(i+1):m A(j,:)=A(j,:)-A(i,:)*A(j,i)/A(i,i); end A end x(m)=A(m,n)/A(m,m); for i=(m-1):-1:1 x(i)=(A(i,n)-A(i,i+1:m)*x(i+1:m))/A(i,i); end x 2、利用高斯列主元消元法 clear format rat A=randn(5) b=randn(5,1) [m,n]=size(A); for i=1:(m-1) numb=int2str(i); disp([第,numb,次选列主元后的增广矩阵]) temp=max(abs(A(i:m,i))); [a,b]=find(abs(A(i:m,i))==temp); tempo=A(a(1)+i-1,:); A(a(1)+i-1,:)=A(i,:); A(i,:)=tempo disp([第,numb,次消元后的增广矩阵]) for j=(i+1):m A(j,:)=A(j,:)-A(i,:)*A(j,i)/A(i,i); end A end %回代过程 disp(回代求解) x(m)=A(m,n)/A(m,m); for i=(m-1):-1:1 x(i)=(A(i,n)-A(i,i+1:m)*x(i+1:m))/A(i,i); end x 3、LU分解 A = randn(5); [L,U,P] = lu(A) L = P\L 附录2 n =5; A = randn(n); for i = 1:n A(i,i) = sum(abs(A(i,:))); end b = randn(n, 1); D = diag(diag(A)); L=tril(A,-1) U=triu(A,1) M = D\(L+U); f = D\b; Rho = max(abs(eig(M))) R = 1e-07; switch sign(1 - Rho) case -1 disp(The Jacobian Method is not applicable) otherwise x(:,1) = normrnd(0, 9, n, 1); k = 1; while k = 50*n x(:, k+1) = M*x(:, k) + f; if norm(x(:, k+1) - x(:, k))=R k = k+1; else X = x(:, k+1) IterN = k break end end end Y = A\b ER_1 = norm(Y - X, 1) ER_2 = norm(Y - X, 2) ER_inf = norm(Y - X, inf) 改变谱半径的大小 看收敛速度的变化 n =5; A = randn(n) for i = 1:n A(i,i) = sum(abs(A(i,:))); end b = randn(n, 1) D = diag(diag(A)); L=tril(A,-1); U=triu(A,1); M = D\(L+U); f = D\b; Rho = max(abs(eig(M)))/2 R = 1e-07; switch sign(1 - Rho) case -1 disp(The Jacobian Method is not applicable) otherwise x(:,1) = normrnd(0, 9, n, 1); k = 1; while k = 50*n x(:, k+1) = M*x(:, k) + f;