原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
数值分析课程设计二
实验报告二实验时间：2012.10.16地点：理学院机房报告人：王振宇一． 题目分别利用雅可比迭代法，赛德尔迭代法解线性方程组10x-y=9-x+10y-2z=7-2y+10z=6二．实验目标利用matlab实现雅克比迭代法，赛德尔迭代法，并显示迭代过程及迭代次数。二． 数学原理雅克比迭代法由aii≠0（i=1,2，……n），选取M为A的对角元素部分，即选取M=D（对角阵），A=D-N，由此而得到雅克比迭代法x0（初始向量）x（k+1)=Bx(k)+f(k=0,1，……)其中B=I-D-1A=D-1L+U≡J,f=D-1b.称J为解Ax=b的雅克比迭代法的迭代阵。记x(k)=(x1k,……xik,……x1k)T，由雅克比迭代公式有Dxk+1=L+Uxk+b在此处键入公式。于是解Ax=b的雅克比迭代法的计算公式为x(0)==（x10,……xn0）Txi(k+1)=(bi-j=1j≠inaijxj(k))/aiii=1,2……n（k=0,1，……表示迭代次数）赛德尔迭代法选取分裂矩阵M为A的下三角部分，即选取M=D-L（下三角阵），A=M=N，于是得到解Ax=b的赛德尔迭代法x0（初始向量）x（k+1)=Bx(k)+f(k=0,1，……)其中B=I-D-L-1A=D-L-1U≡G,f=D-L-1b,称G=D-L-1U为解Ax=b的赛德尔迭代法的迭代阵下面给出赛德尔迭代法的分量计算公式，记x(k)=(x1k,……xik,……x1k)T，于是有 （D-L)xk+1=Uxk+b于是解Ax=b的赛德尔迭代法的计算公式为x(0)==（x10,……xn0）Txi(k+1)=(bi-j=1i-1aijxjk+1-j=i+1naijxj(k))/aiii=1,2……n（k=0,1，……表示迭代次数）四．程序设计(请将Matlab代码打印页附在该报告扉页)雅克比迭代法functionx=jacobi(A,b,P,d,n)%A为n维非奇异矩阵b为n维向量P为初值d为误差界n为给定迭代次数N=length(b);fork=1:nforj=1:Nx(j)=(b(j)-A(j,[1:j-1,j+1:N])*P([1:j-1,j+1:N]))/A(j,j);enderr=abs(norm(x-P));P=x;if(errd)break;endendPx=x;k,errendA=[10-10;-110-2;0-210];b=[976];P=[000];x=jacobi(A,b,P,1e-4,20)P=0.9000000000000000.7000000000000000.600000000000000P=0.9700000000000000.9100000000000000.740000000000000P=0.9910000000000000.9450000000000000.782000000000000P=0.9945000000000000.9555000000000000.789000000000000P=0.9955500000000000.9572500000000000.791100000000000P=0.9957250000000000.9577750000000000.791450000000000P=0.9957775000000000.9578625000000000.791555000000000P=0.9957862500000000.9578887500000000.791572500000000P=0.9957888750000000.9578931250000000.791577750000000P=0.9957888750000000.9578931250000000.791577750000000k=9err=7.320775232161408e-006x=0.9957888750000000.9578931250000000.791577750000000赛德尔迭代法functionx=seidel(A,b,P,d,n)%A为n维非奇异矩阵b为n维向量P为初值d为误差界n为给定迭代次数N=length(A);fork=1;nforj=1:Nifj==1x(1)=(b(1)-A(1,2:N)*P(2:N))/A(1,1);elseifj==Nx(N)=(b(N)-A(N,1:N-1)*(x(1:N-1)))/A(N,N);elsex(j)=(b(j)-A(j,1:j-1)*x(1:j-1)-A(j,j+1:N)*P(j+1:N))/A(j,j);endenderr=abs(norm(x-P));P=xif(errd)break;endendx=x;k,errendy=seidel(A,b,P,1e-4,20)P=0.90
文档评论（0）