雅克比迭代法 实验二用雅克比迭代法解方程组.doc

雅克比迭代法 实验二用雅克比迭代法解方程组 数值计算方法实验 实验2 用雅克比迭代法解方程组 一、 实验目的 1、 学习使用雅克比迭代法求解方程组，加深对雅克比迭代法的认识。 2、 在误差允许的范围内，在循环次数较少的情况下，求解方程组的稳定解。 3、 进一步加深对Matlab的学习。 二、 … 2010年-2012年双向转诊管理实施方案 为促进社区卫生事业全面,协调,可持续发展,做到医疗卫生机构优势互补,资源共享,为广大人 民群众提供优质,便捷,廉价的医疗卫生服务,解决人民群众”看病难,看病贵”的问题,特制定本 实施方案. 一,双向转诊原则… 西南大学荣昌校区 《动物生物化学》复习思考题 （供动物医学、动物科学专业用） 动物生理生化教研室 第一章：绪论 一、思考题： 1、生物化学的研究内容是什么？概述生物化学的发展阶段。 2、生物化学和畜牧兽医专业之间有何联系？ 二、名词解释 生物化学 静… 数值计算方法实验 实验2 用雅克比迭代法解方程组 一、 实验目的 1、 学习使用雅克比迭代法求解方程组，加深对雅克比迭代法的认识。 2、 在误差允许的范围内，在循环次数较少的情况下，求解方程组的稳定解。 3、 进一步加深对Matlab的学习。 二、 实验题目 用雅克比迭代法解方程组： 三、 实验原理 设有n元线性方程组AX=b, 其中系数矩阵A的对角元素 0（i=1,2,3,…,n） ,可得到以下等价方程 从方程组的第i个方程中可以解出组： (( ) , ) , ………………………………………………， . 记D为A的对角线矩阵，L为A的下三角阵，U为A的上三角阵， 即 则有： X=-(L+U)X+ b; b; 令 B=-(L+U)， d= 则有：X=BX+d; 其中B成为迭代矩阵 上式即为雅克比迭代法的迭代公式。 取定一个 以后，便可得 , , k=0,1,2,…, }.若 ,则 就 , 再往下迭代得到：如此反复迭代，一般的有：由此便得到一个向量序列{ 是方程组的解。反之：若就不能用雅克比迭代法计算。 可将上述迭代过程改为如下形式： 不存在，即{}不收敛， (i=1,2,3,…,n) , 四、 实验内容 按照上述迭代法可得到如下的算法： 设方程组为：AX=b, 其中系数矩阵A的对角元素 0（i=1,2,3,…,n） max_iter为迭代数容许的最大值，eps为容许误差： (0) ,...,1. 取初始向量x??x1(0),x(0) xn?,令k=0; 2 2. 对i=1,2,…,n,计算： ??n k? xik?1?1?. ??baiijxj??j?1?aii?j?i?? 3. 如果?xik?1?xik i?1 n ?eps,则输出xi k?1 ,结束。否则，执行下一步； 4. 如果max_iter=50,则该迭代法不收敛，终止程序；否则转2. 建立一个M文件，程序如下： function [x,iter]=jacobi1(A,b,x_0) A=[5 2 1;2 8 -3;1 -3 -6]; b=[8 21 1]; eps=10e-6; x_0=[0 0 0]; n=3; x=x_0; max_iter=100; y=[0 0 0]; for k=1:n if(A(k,k)==0) return end end iter=1; while iterfor i=1:n for j=1:n temp=temp+A(i,j)*x(j); end temp=(b(i)-temp)/A(i,i); y(i)=x(i)+temp; end x=y; iter=iter+1; if abs(temp)return end end 五、 实验结果 ans = 1.0000 2.0000 -1.0000 六、 实验结果分析 该方程组用雅克比迭代法求解是收敛的。 七、 实验出现的问题及体会 1、用M函数文件创建的函数可以输出两个参量。 数值计算方法实验 实验2 用雅克比迭代法解方程组 一、 实验目的 1、 学习使用雅克比迭代法求解方程组，加深对雅克比迭代法的认识。 2、 在误差允许的范围内，在循环次数较少的情况下，求解方程组的稳定解。 3、 进一步加深对Matlab的学习。 二、 … 数值计算方法实验 实验2 用雅克比迭代法解方程组 一、 实验目的 1、 学习使用雅克比迭代法求解方程组，加深对雅克比迭代法的认识。 2、 在误差允许的范围内，在循环次数较少的情况下，求解方程组的稳定解。 3、 进一步加深对Matlab的学习。 二、 … 数值计算方法实验 实验2 用雅克比迭代法解方