约当消去法与高斯消去法.docVIP

约当消去法和高斯消去法的应用 题目： 求解方程组： 误差分析范围为e=0.0001. 二、引言 线性方程组的解法大致分为直接法与迭代法两大类。迭代法的突出优点是算法简单，编程容易。但迭代法的缺点是它要求迭代过程具有收敛性，发散的迭代过程是没有实用价值的。 前面我们通过对迭代法的学习，已经大致掌握了有关迭代法的基本思想，此处我们将学习解决线性方程组的另一种方法——直接法 这种算法思想简单并且容易掌握，结果精确，这虽然是一种古老的算法，但用在现代计算机上仍十分有效，所以这依然是一种十分有效的解决线性方程组的一种方法，下面我们一起来学习这种算法！ 三、算法 直接法中我们接触到的有约当消去法和高斯消去法，下面我们将逐一介绍两种算法： 约当消去法： 考察一般形式的线性方程组 ,i=1，2，3，…，n （1） 第一步，先把方程（1）中的第一个方程的 的系数化为1，并从其余的方程中消去，得到方程（1）1 第二步，先把方程中的第二个方程的的系数化为1，并从其余的方程中消去，得到方程（1）2 第K步，先把方程中的第k个方程的的系数化为1，并从其的方程中消去，得到方程（1）K 第n步，先把方程中的第n个方程的系数化为1，并从其与方程中消去。这样经过n步的消去就可以得到原方程组的解。 2、高斯消去法 高斯消去法包括两个过程：消去过程和回代过程 （1）消去过程： 第一步，把方程中的第一个方程的x（1）（1）1 Program fjs3 implicit none integer::k,i,j integer,parameter::m=3,n=3 real::a(m,n)=(/10,-1,-1,-1,10,-1,-2,-2,5/) real::b(m)=(/7.2,8.3,4.2/) do k=1,n do j=k+1,n a(k,j)=a(k,j)/a(k,k) end do b(k)=b(k)/a(k,k) do i=1,n if(i/=k) then do j=k+1,n a(i,j)=a(i,j)-a(i,k)*a(k,j) end do b(i)=b(i)-a(i,k)*b(k) end if end do end do write(*,*) b(1),b(2),b(3) stop end 1.100000 1.200000 1.300000 Press any key to continue 高斯消去法 program fjs3 implicit none integer::k,i,j integer,parameter::m=3,n=3 real::a(m,n)=(/10,-1,-1,-1,10,-1,-2,-2,5/) real::b(m)=(/7.2,8.3,4.2/) do k=1,n do j=k+1,n a(k,j)=a(k,j)/a(k,k) end do b(k)=b(k)/a(k,k) do i=k+1,n do j=k+1,n a(i,j)=a(i,j)-a(i,k)*a(k,j) end do end do do i=k+1,n b(i)=b(i)-a(i,k)*b(k) end do end do do i=n-1,1,-1 do j=i+1,n