数值分析编程及运行结果(高斯顺序消元法).pdf

数值分析编程及运行结果(高斯顺序消元法) 高斯消元法 1.程序: clear format rat A=input(输入增广矩阵 A=) [m,n]=size(A); for i=1:(m-1) numb=int2str(i); 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 数值分析编程及运行结果(高斯顺序消元法) . 2.运行结果： . 数值分析编程及运行结果(高斯顺序消元法) . 高斯选列主元消元法 1. 程序: clear format rat A=input(输入增广矩阵 A=) [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 2.运行结果： . 数值分析编程及运行结果(高斯顺序消元法) . . 数值分析编程及运行结果(高斯顺序消元法) . 追赶法 1. 程序： function [x,L,U]=zhuiganfa(a,b,c,f) a=input( 输入矩阵-1 对角元素 a=); b=input(输入矩 阵对角元素 b=); c=input(输入矩阵+1 对角 元素 c=); f=input(输入增广矩阵最后一列元 素 f=); n=length(b); % 对 A 进行分解 u(1)=b(1); for i=2:n