高斯消去法上机实验报告.docxVIP

高斯消去法上机实验报告 　　?2?????????+??????=??　　??????+??????+??????=??????+??????=??　　一、课题名称:高斯消去法　　二、引言　　为了节省计算量，对约当消去法进行改进　　三、算法　　1.先将方程中x1的系数化为1，并从方程组1的其余方程中消去x1,得x1?+=　　2?x3=2　　2+=　　2.再将方程2中的x2系数化为1，并从方程3中消去x2　　x1?+=　　2?=　　3=?6　　3.将方程3回代到方程2中，然后再回代到方程1最终得到所求的解x1=9,x2=?1,x3=?6,　　高斯消去法分为消元过程和回代过程两个环节　　1.消元过程对k?1,2,?,n?1　　(1)选主元，找　　(2)如果ik??k,k?1,?,n?使得aik,k?maxaikk?i?naik,k?0，则矩阵A奇异，程序结束；否则执行(3)。　　a?aikj(3)如果ik?k，则交换第k行与第ik行对应元素位置，kj，　　j?k,?,n?1。　　(4)消元，对i?k,?,n，计算lik?aik/akk,对j?k?1,?,n?1，计算aij?aij?likakj.　　2.回代过程　　(1)若ann?0,则矩阵奇异，程序结束；否则执行(2)。　　n??xi??ai,n?1??aijxj?/aiix?an,n?1/ann;j?i?1??(2)n对i?n?1,?,2,1，计算　　四、程序设计　　programmain　　implicitnone　　integer::k,j,i　　real::s　　integer,parameter::m=3,n=3　　real::a(m,n)=(/10,-1,-1,-1,10,-1,-2,-2,5/)　　real::b(m)=(/,,/)　　dok=1,n　　doj=k+1,n　　a(k,j)=a(k,j)/a(k,k)　　enddo　　b(k)=b(k)/a(k,k)　　doi=k+1,n　　doj=k+1,n　　a(i,j)=a(i,j)-a(i,k)*a(k,j)　　enddo　　b(i)=b(i)-a(i,k)*b(k)　　enddo　　enddo　　doi=n-1,1,-1　　s=0　　doj=i+1,n　　s=s+a(i,j)*b(j)　　enddo　　b(i)=b(i)-s　　enddo　　write(*,*)b(1),b(2),b(3)　　end　　五、结果及讨论分析　　Pressanykeytocontinue　　优点：这种算法是对约当消去法的一种改进，这种改进明显的减少了计算量。当n充分大是，高斯消去法的消元过程约需3可以忽略不计，因此，它比约当消去法大约节省了50%的计算量。　　缺点：由于添加了回代过程，算法的结构较约当法稍复杂。在用高斯消去法求解一般的线性方程时，即使akk不为0，但如果其绝对值很小，舍入误差的影响也会严重的损失精度，实际计算时要尽量避免这类情况的发生。　　n3　　深圳大学实验报告　　课程名称：　　实验项目名称：　　学院：计算机与软件学院　　华　　科技大数值分析实验报告　　学号姓名类别硕士　　XX年5月6日　　中学　　实验　　实验要求：　　根据教材实验做出相应改编：分别使用Gauss消元、列选主元。全选主元的方法求解线性方程组，分别比较三种消元方法的结果和算法的区别，并说明主元的选取在Gauss消元的中的作用。问题提出：　　Gauss消去法是我们在线性代数中已经熟悉的。但由于计算机的数值运算是在一个有限的浮点数集合上进行的，如何才能确保Gauss消去法作为数值算法的稳定性呢？Gauss消去法从理论算法到数值算法，其关键是主元的选择。主元的选择从数学理论上看起来平凡，它却是数值分析中十分典型的问题。一般来说书本上采用的列选主元的办法对其线性方程组进行求解的，那么我们是否可以选择一种行列都选取主元消去的办法来减小相应的误差呢？全主元消元法和列主元消元　　法一样都是由高斯消元法演变而来。只不过选取主元的范围有所加大。全选主元相对于列选主元的更加复杂化了，因为在运算的过程中导致了元的位置发生了变化，这样我们就不得不追踪每个元的位置。本次实验就几个问题进行了matlab实验分析，比较几种计算方法的优劣性。　　实验内容：　　考虑线性方程组　　Ax?b,A?Rn?n,b?Rn　　编制一个程序：分别能进行Gauss消去、列选主元Gauss消去、全选主元Gauss消去法进行解线性方程组。对三种算法所得到的结果进行比较，分析三种计算方法的准确性。　　具体内容：　　?61??7?????861???