太阳集团官网1385设Ax=b其中A=DLU为非奇异矩阵且对角阵D也非奇异则当迭代矩阵J的谱半径ρ(J)时雅克比迭代法收

当前位置:太阳集团www.1385.com > 太阳集团官网1385 > 太阳集团官网1385设Ax=b其中A=DLU为非奇异矩阵且对角阵D也非奇异则当迭代矩阵J的谱半径ρ(J)时雅克比迭代法收
作者: 太阳集团www.1385.com|来源: http://www.sifanghuixin.com|栏目:太阳集团官网1385

文章关键词:太阳集团www.1385.com,逐次超松弛

  登录成功,如需使用密码登录,请先进入【个人中心】-【账号管理】-【设置密码】完成设置

  简介:本文档为《数值分析课程实验报告-雅可比和逐次超松弛迭代法doc》,可适用于综合领域

  *若权利人发现爱问平台上用户上传内容侵犯了其作品的信息网络传播权等合法权益时,请按照平台侵权处理要求书面通知爱问!

  关于数值分析课程实验报告-雅可比和逐次超松弛迭代法.doc文档,爱问共享资料拥有内容丰富的相关文档,站内每天千位行业名人共享最新资料。

  数值分析课程实验报告雅可比和逐次超松弛迭代法《数值分析》课程实验报告用雅可比和逐次超松弛迭代法求解方程组算法名称用雅可比和逐次超松弛迭代法求解方程组学科专业机械工程作者姓名程习康作者学号作者班级机电院研究生级C班中南大学二〇一五年十二月《数值分析》课程实验报告《数值分析》课程实验报告用雅可比和逐次超松弛迭代法求解方程组实验名称成绩一、问题背景在自然科学和工程技术中有很多问题的解决常常归结为解线性方程组例如电学中的网络问题化学中的配平方程式模型问题船体数学放样中建立三次样条函数问题用最小二乘法求实验数据的曲线拟合问题非线性方程组求解问题用差分法或者有限元法解常微分方程、偏微分方程边值问题等都导致求解线性方程组。在实践中通常采用数值方法来讨论线性方程组Ax=b的解其中迭代法是一种重要方法。二、数学模型用雅可比迭代法解方程组x,,,,,,,x,,,,,,,,,,,,,,,x用逐次超松弛迭代法解方程组(取ω=ep=e,itmax=),xxx,,,,,,xxx,,,,x,xx,,精确到小数点后四位。三、算法描述雅可比迭代法:首先将方程组中的系数矩阵A分解成三部分即:A=LDU如图所示其中D为对角阵L为下三角矩阵U为上三角矩阵。之后确定迭代格式X^(k)=B*X^(k)f(这里^表示的是上标括号内数字即迭代次数)如图所示其中B称为迭代矩阵雅克比迭代法中一般记为J。(k=,)再选取初始迭代向量X^()开始逐次迭代。设Ax=b其中A=DLU为非奇异矩阵且对角阵D也非奇异则当迭代矩阵J的谱半径ρ(J)时雅克比迭代法收敛。逐次超松弛迭代法:逐次超松弛(SuccessiveOverRelaxation)迭代法简称SOR迭代法它是在GS法基础上为提高收敛速度采用加权平均而得到的新算法设解方程的GS法记为()再由与加权平均得《数值分析》课程实验报告这里ω,称为松弛参数将()代入则得()该法称为SOR迭代法ω,称为松弛因子当ω=时()式即为高斯赛德尔迭代法简记GS法将()写成矩阵形式则得即于是得SOR迭代的矩阵表示()其中()分解后有。四、主要代码雅可比迭代法建立函数文件:functionx,k,index=Jacobi(A,b,ep,itmax)ifnarginitmax=endifnarginep=eendn=length(A)k=x=zeros(n,)y=zeros(n,)index=whilek=itmaxfori=:nifabs(A(i,i))eindex=returnendy(i)=(b(i)A(i,:n)*x(:n)A(i,i)*x(i))A(i,i)endifnorm(yx,inf)epbreak《数值分析》课程实验报告endk=kx=yend逐次超松弛迭代法建立函数文件:functionx,k,index=SOR(A,b,w,ep,itmax)ifnarginitmax=endifnarginep=eendn=length(A)k=w=x=zeros(n,)y=zeros(n,)index=whilek=itmaxfori=:nifabs(A(i,i))eindex=returnendifi==y(i)=(w)*x(i)w*(b(i)A(i,i:n)*x(i:n))A(i,i)elseifi==ny(i)=(w)*x(i)w*(b(i)A(i,:i)*y(:i))A(i,i)elsey(i)=(w)*x(i)w*(b(i)A(i,:i)*y(:i)A(i,i:n)*x(i:n))A(i,i)endendifnorm(yx,inf)epbreakendk=kx=yEnd五、太阳集团官网1385太阳集团官网1385实验结果及分析雅可比迭代法在命令窗口输入:A=b=ep=eitmax=x,k,index=Jacobi(A,b,ep,itmax)得到结果如下:x=k=《数值分析》课程实验报告index=可知方程组的解为,,,迭代次数为由index=表示计算成功。雅克比迭代法的优点明显计算公式简单每迭代一次只需计算一次矩阵和向量的乘法且计算过程中原始矩阵A始终不变比较容易并行计算。然而这种迭代方式收敛速度较慢而且占据的存储空间较大所以工程中一般不直接用雅克比迭代法而用其改进方法。逐次超松弛迭代法在命令窗口输入:A=b=w=ep=eitmax=x,k,index=SOR(A,b,w,ep,itmax)得到结果如下:x=k=index=可知方程组的解为,,,迭代次数为由index=可知计算成功。可知逐次超松弛迭代法与雅可比迭代法相比,收敛速度较快。由逐次超松弛迭代法求出的方程组的数值解与该方程组的精确解十分接近,离散化后线性方程组的逐次超松弛迭代法的精确性较高。逐次超松弛迭代法可以广泛地应用于实际。该算法不仅可以用来求解高阶稀疏线性方程组,还可以用来求解热传导问题这样可以大大减少计算量和计算机的内存储量,从而提高计算效率。

  1620年,著名的“五月花”号船满载被迫害的清教徒到达美洲,但随即他们遭遇寒冬,在印第安人的帮助下,移民们度过困难,为感谢印第安人的真诚帮助,太阳集团官网1385便诞生了感恩节。

网友评论

我的2016年度评论盘点
还没有评论,快来抢沙发吧!