计算方法实验六高荣双.doc

实 验 报 告 课程名称： 计算方法 院 系： 数学科学系 专业班级： 数应1001B 学 号： 1031110107 学生姓名： 高荣双 指导教师： 沈林 开课时间： 2012 至 2013 学年第 一 学期 一、学生撰写要求 按照实验课程培养方案的要求，每门实验课程中的每一个实验项目完成后，每位参加实验的学生均须在实验教师规定的时间内独立完成一份实验报告，不得抄袭，不得缺交。 学生撰写实验报告时应严格按照本实验报告规定的内容和要求填写。字迹工整，文字简练，数据齐全，图表规范，计算正确，分析充分、具体、定量。 二、教师评阅与装订要求 1.实验报告批改要深入细致，批改过程中要发现和纠正学生实验报告中的问题，给出评语和实验报告成绩，签名并注明批改日期。实验报告批改完成后，应采用适当的形式将学生实验报告中存在的问题及时反馈给学生。 2.实验报告成绩用百分制评定，并给出成绩评定的依据或评分标准（附于实验报告成绩登记表后）。对迟交实验报告的学生要酌情扣分，对缺交和抄袭实验报告的学生应及时批评教育，并对该次实验报告的分数以零分处理。对单独设课的实验课程，如学生抄袭或缺交实验报告达该课程全学期实验报告总次数三分之一以上，不得同意其参加本课程的考核。 3.各实验项目的实验报告成绩登记在实验报告成绩登记表中。本学期实验项目全部完成后，给定实验报告综合成绩。 4.实验报告综合成绩应按课程教学大纲规定比例（一般为10-15%）计入实验课总评成绩；实验总评成绩原则上应包括考勤、实验报告、考核（操作、理论）等多方面成绩； 5.实验教师每学期负责对拟存档的学生实验报告按课程、学生收齐并装订，按如下顺序装订成册：实验报告封面、实验报告成绩登记表、实验报告成绩评定依据、实验报告（按教学进度表规定的实验项目顺序排序）。装订时统一靠左侧按“两钉三等分”原则装订。 实验名称 数值积分与数值微分 实验时间 2012年12月13日 学生姓名 高荣双 实验地点 9#405数学实验室 1、实验所用软件 WIN7操作系统、Matlab 2、实验目的 设计程序实现雅可比迭代法，高斯-赛德尔迭代法等算法。掌握方程求根的基本理论。并用实例在计算机上计算。 3、实验内容 1、用雅可比迭代法求解方程组：x1+2*x2-2*x3=1; x1+x2+x3=1; 2*x1+2*x2+x3=1. 2、用高斯-赛德尔迭代法求解方程组：x1+0.4*x2+0.4*x3=1; 0.4*x1+*x2+0.8*x3=2; 0.4*x1+0.8*x2+x3=3. 4、实验方法、步骤 1、function varargout=Jacobiliu(varargin) A=[1,2,-2;1 1 1;2 2 1]; b=[1;1;1]; x0=[0;0;0]; x=Jacobi(A,b,x0) function x=Jacobi(A,b,x0) n=100;%最大迭代次数 ee=0.0001;%精度 n1=length(b); x1=zeros(n1,1); for i=1:n for j=1:n1 s=0; for k=1:n1 if k~=j s=s+A(j,k)*x0(k); end end x1(j)=(b(j)-s)/A(j,j); end if norm(x1-x0)ee break end x0=x1; end x=x1 2、function x=Gauss_Seidel(A,b,x0,tol) if (nargin==2) x0=ones(size(b)); tol=1e-6; elseif (nargin==3) tol=1e-6; else sprintf(USAGE:Gauss_Seidel(A,b,x0,tol)) end D=diag(diag(A)); U=triu(A,1); L=tril(A,-1); G=-(D+L)\U; d1=(D+L)\b; x=G*x0+d1; n=1; while norm(x-x0)=tol