下载此文档

列主高斯消去法.doc


文档分类:IT计算机 | 页数:约9页 举报非法文档有奖
1/9
下载提示
  • 1.该资料是网友上传的,本站提供全文预览,预览什么样,下载就什么样。
  • 2.下载该文档所得收入归上传者、原创者。
  • 3.下载的文档,不会出现我们的网址水印。
1/9 下载此文档
文档列表 文档介绍
课程名称:数值分析项目名称:线性方程组的直接解法姓名:专业:班级:(1)学号:20104161003同组成员实验准备注:1、实验准备部分包括实验环境准备和实验所需知识点准备。2、若是单人单组实验,同组成员填无。:线性方程组的直接解法用列主高斯消去法解下面线性方程组输出解向量x实验过程记录注:实验过程记录要包含实验步骤,页码不够可自行添加。:%列主元Gauss校区法解线性方程组A=[4111;13-11;1-120;1102;];%系数矩阵b=[]';%n维向量y=inv(A)*b%matlab的计算结果n=length(b);%方程个数nx=zeros(n,1);%未知向量%-------------消去-----------fork=1:n-1Auk=A(k:n,k);[m,u]=max(abs(Auk));u=u+k-1%u为最大元所在的列%------交换最大的行和当前行的值-------forj=k:ntemp=A(u,j);A(u,j)=A(k,j);A(k,j)=temp;endtemp=b(k);b(k)=b(u);b(u)=temp;%ifA(k,k)==0;%error('Error');%endfori=k+1:n%A(i,k)=A(i,k)/A(k,k);Aik=A(i,k)/A(k,k)forj=k:nA(i,j)=A(i,j)-Aik*A(k,j);endAb(i)=b(i)-Aik*b(k)endend%-------------回代-----------x(n)=b(n)/A(n,n)fork=n-1:-1:1S=b(k);forj=k+1:nS=S-A(k,j)*x(j);endx(k)=S/A(k,k)endx%程序的计算结果error=abs(x-ones(n,1))%误差%误差小于直接进行高斯消去的计算误差实验操作成绩(百分制)__________实验指导教师签字:__________实验小结:注:1、实验小结应包含实验所需知识点和实验方法的总结,实验心得体会等。2、分组实验需包含同组讨论内容。y=--=1Aik==--=-==---=--==----=--=2Aik=-=--=-

列主高斯消去法 来自淘豆网www.taodocs.com转载请标明出处.

相关文档 更多>>
非法内容举报中心
文档信息
  • 页数9
  • 收藏数0 收藏
  • 顶次数0
  • 上传人drp539601
  • 文件大小37 KB
  • 时间2019-11-19