9512.net
甜梦文库
当前位置:首页 >> 数学 >>

Gauss全主元消去法


Gaussian elimination with complete pivoting
To keep the Gaussian elimination method smooth operation,We must ensure that the pivot element a
( k ?1) kk

?0

and a

( k ?1) kk

can not too small in each

operation.So the selection of appropriate pivot elements in every step of elimination.,the maximum absolute value or the major elements as the pivot element, this improved Gaussian elimination method named Gauss pivot element elimination. There have 3 methods to select pivot element,Gaussian elimination with partial pivoting 、 scaled partial pivoting、complete pivoting.Through access to some data,I the Gaussian elimination with complete pivoting briefly. Complete entries a
( k ?1) ij

introduce

pivoting
(i ? k , j ? k )

at

the

k

th

step

searches

all

the

,to find the entry to the largest magnitude.Both row

and column interchanges are performed to bring this entry to the pivot position. The specific steps:
? (0) ? (0) ? ? ? ( 0) ? ( 0) , a1n xn a1,n?1 ?a11 x1 a12 x2 (0) (0) ( 0) ( 0) ? ?a21 x1 ? a22 x2 ? ? ? a2 n xn ? a2,n ?1 , ? ?? ? (0) (0) ( 0) ( 0) ? ?an1 x1 ? an 2 x2 ? ? ? ann xn ? an ,n ?1

The

k th

elimination,find all a
( k ?1)

( k ?1) ij

(i=k, k+1,,n; j=k,k+1 ? n),the

a ??

( k ?1)

is

the pivot element. a??

? maxaij

( k ?1)

Both row and column interchanges are performed to bring this entry to the pivot position.Move

a ??

( k ?1)

to

(k,k).Until

the

n-1 step, the

original equation into the same equations on triangular solution:
? ( 0) ? (0) ? ( 0) ? ? ? (0) ? ( 0) , a1n xn a1,n?1 ?a11 x1 a12 x2 a13 x3 ( 1 ) ( 1 ) (1) (1) ? ? a23 x3 ? ? ? a2 n xn ? a2,n ?1 , a x 22 2 ? ( 2) ( 2) ( 2) ? ? a33 x3 ? ? ? a3n xn ? a3,n?1 , ? ? ? ( n ?1) ( n ?1) ? ? an ,n ?1 , a x nn n ? ?

Backward substitution: Backward substitution is the solution of triangular systems .If a then
1 ( n ?1) ? ? xn ? ( n ?1) a n ,n ?1 , an,n ? ? n ? ( i ?1) ? ? 1 ? (i ?1) ? , i ? n ? i, n ? 2,? ,1. ? ? x a a x ( i ? 1 ) i , n ?1 ij j? ? i j ? i ? 1 ? aii ? ?
( n ?1) nn

?0



Algorithm:
function x=Gauss(A,b) B=[A,b]; n=length(A); x=zeros(n,2); for i=1:n x(i,1)=i; %x 第一列数字 i 代表 x(i) end for i=1:n-1 [m,p]=max(abs(B(i:n,i:n)));%求绝对值最大系数所在行号和列号 [l,q]=max(m); if m==0; disp('A 奇异') return

end if p(q)+i-1>i %行变换 pm=p(q)+i-1; C(i:n+1)=B(i,i:n+1); B(i,i:n+1)=B(pm,i:n+1); B(pm,i:n+1)=C(i:n+1); end if q+i-1>i %列变换 q=q+i-1; D(1:n)=B(1:n,q); B(1:n,q)=B(1:n,i); B(1:n,i)=D(1:n); E(1:2)=x(q,:);%调换解的顺序 x(q,:)=x(i,:); x(i,:)=E(1:2); end for k=i+1:n %消元 am=B(k,i)/B(i,i); B(k,i:n+1)=B(k,i:n+1)-am*B(i,i:n+1); end end x(n,2)=B(n,n+1)/B(n,n); for i=n-1:-1:1 %回代 xx(i+1:n)=x(i+1:n,2); x(i,2)=(B(i,n+1)-B(i,i+1:n)*xx(i+1:n)')/B(i,i); end for i=1:n-1 % 调换解的顺序 [u,r]=min(x(i:n,1)); if r+i-1>i r=r+i-1; xf(1:2)=x(i,1:2); x(i,1:2)=x(r,1:2); x(r,1:2)=xf(1:2); end end end For example P89,Exercises2

>> A=[1.19 2.11 -100 1;14.2 -0.122 12.2 -1;0 100 -99.9 1;15.3 0.110 -13.1 -1];b=[1.12;3.44;2.15;4.16]; x=Gauss(A,b)

x= 1.0000 2.0000 3.0000 4.0000 0.1768 0.0127 -0.0207 -1.1826



更多相关文章:
Gauss完全主元消去法.doc
Gauss完全主元消去法 - 2011-2012(1)专业课程实践论文 ( )专业课程实践论文 Gauss 完全主元消去法 韩璐,0718180109,R 数学 07-1 一、算法理论 ...
Gauss全主元消去法.doc
Gauss全主元消去法 - Gaussian elimination with
高斯主元消去法_图文.ppt
高斯主元消去法 - 列主元高斯消去法 高斯-约当消去法... 高斯主元消去法_数学_自然科学_专业资料。列主元高斯消去法 高斯-约当消去法 §4 高斯选主元素...
Matlab程序Gauss主元消去法.pdf
Matlab程序Gauss主元消去法 - .Gauss主元消去法求解线性方
Matlab程序Gauss主元消去法.doc
Matlab程序Gauss主元消去法 - .Gauss主元消去法求解线性方
5.2 Gauss主元消去法_图文.ppt
5.2 Gauss主元消去法_数学_自然科学_专业资料。第二章 解线性方程组的直接法§ 5.2 列主元Gauss消去法 Gauss主元消去法 一、Gauss主元消去法的引入 ...
高斯主元消去法0425121523.doc
高斯主元消去法0425121523 - 数值分析大作业 --——(高斯主元消去法求解线性方程组) 高斯主元消去法求解线性方程组) 高斯主元消去法求解线性方程组 ...
计算方法-6.1Gauss消去法_图文.ppt
除了列主元消去法,求解线性方程组常用的还 有全主元消去法。 2016/8/14 30 四、 Gauss-Jordan消去法 前面所述的消去法均要进行两个过程 , 即消元 过程和...
高斯主元消去法解线性方程组.pdf
线性方程组的数值解法高斯主元消去法解线性方程组的 MATLAB 实现 班级:MATH 20XX 学号:xxxxxxxxxx 姓名:LI Wen 高斯主元消去法解线性方程组的 MATLAB ...
Matlab中Gauss全主元消元法求解线性方程组的程序.doc
Matlab中Gauss全主元消元法求解线性方程组的程序_数学_自然科学_专业资料 暂无评价|0人阅读|0次下载|举报文档 Matlab中Gauss全主元消元法求解线性方程组的程序_...
Gauss主元消去法.ppt
Gauss主元消去法 - 10 算例:线性方程组的Gauss主元消去法求解
实验二:Gauss主元消去法.doc
实验二:Gauss主元消去法 - 实验二:Gauss主元消去法 程序 1:Gauss主元消去法 A=input('请输入线性方程组的增广矩阵 A='); n=length(A)-1; ...
matlab与全主元消去法编写程序.doc
matlab与全主元消去法编写程序 - function [x,qa]=gaussq(a,b)% 输出想 x 为解,qa 为全主元变换后的 a 矩阵 d=[a,b]; RA=rank(a);RD...
解线性代数方程组的直接法之GAUSS主元消去法及其C 编程....pdf
解线性代数方程组的直接法 Gauss 主元消去法最近学习数值分析,为了能够边学
Gauss主元消去法.doc
Gauss主元消去法 - 贵州师范大学数学与计算机科学学院学生实验报告 课程名
全主元消去法.doc
全主元消去法 - 全主元消去法 1. 通过编制全主元消去法的 MATLAB 程序
Gauss消去法和列主元消去法.doc
MATLAB 与科学计算 实验报告 1 Gauss主元消去法学号 实验题目 班级 姓名 指导教师 评分 易昆南 Gauss主元消去法 1、设计(实习)目的运用高斯主元消去法...
数值分析列主元高斯消去法.doc
再利用计算机中的C语言对高斯主元消去法 可以快速得到线性方程组的解,由简单的线性方程组可以推广到一般n阶线性方程组,这对 如何利用高斯主元消去法解决实际...
高斯主元消去法.doc
高斯主元消去法_数学_自然科学_专业资料。北京化工大学 数学 1403 徐彤 学号:2014016544 高斯主元消去法程序 function x = rowGauss(A,b) [N,N]=size(A...
41%20Gauss消去法pdf_图文.pdf
41%20Gauss消去法pdf - 第四章方程组的直接解法 4.1 Gauss消去法 4.1.1 Gauss消去法的计算过程 4.1.2 矩阵的三角分解 4.1.3 主元素消去法 4...
更多相关标签:

All rights reserved Powered by 甜梦文库 9512.net

copyright ©right 2010-2021。
甜梦文库内容来自网络,如有侵犯请联系客服。zhit325@126.com|网站地图