数值计算方法上机答案.docVIP

下载本文档

159
0
约1.63万字
约 23页
2016-12-30 发布于重庆
举报
版权申诉

数值计算方法上机答案.doc

1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
7、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

数值计算方法上机答案

本科实验报告课程名称：实验项目：实验地点：机房专业班级：采矿1206 学号： 2012002896 学生姓名：康明月指导教师： 2014年 7 月 3 日求非线性方程的根。求方程在附近的是根，要求精度满足.(牛顿切线法) f=inline(x-cos(x)); % f(x) df=inline(1+sin(x)); %f(x) n=1; x0=input(x0=); del=input(del=); N=input(N=); fprintf(\n k x(k)); fprintf(\n %2d %f ,0,x0); F0=f(x0); dF0=df(x0); while nN if dF0==0 fprintf(导数为0,迭代无法继续进行.); return; end x1=x0-F0/dF0; F1=f(x1); dF1=df(x1); if ((abs(x1-x0)del) |abs(F1)del) fprintf(\n \n 结果：%f\n,x1); return; end fprintf(\n %2d %f ,n,x1); n=n+1; x0=x1; F0=F1; dF0=dF1; end fprintf(\n\n % d 次迭代后未达到精度要求.\n,N); NewtonIteration x0=1.5 del=1e-4 N=100 k x(k) 0 1.500000 1 0.784472 2 0.739519 结果：0.739085 求方程在附近的是根，求出具有思维有效数字的根近似值..(简单迭代法) clear clc phi=inline((0.8+x^2)^(1/3)); %迭代函数 x0=input(x0=); del=input(del=); N=input(N=); n=1; fprintf(\n %2d %f ,0,x0); while nN x=phi(x0); if abs(x-x0)del fprintf(\n \n 近似解=%f \n,x); return end fprintf(\n %2d %f ,n,x); n=n+1; x0=x; end fprintf(\n \n %f d次迭代后未达到精度要求. \n,N) x0=1 del=1e-4 N=100 0 1.000000 1 1.216440 2 1.316116 3 1.363004 4 1.385180 5 1.395688 6 1.400671 7 1.403034 8 1.404155 9 1.404687 10 1.404939 11 1.405059 近似解=1.405116 100.000000 d次迭代后未达到精度要求. 求解线性方程组（直接法或迭代法） 1、 (列主元素消元法) a=input(a=)%[2,2,1,-3,8;-2,1,-1,-3,1;8,-1,3,8,-1;10,4,4,3,8]; [p,n]=size(a); for w=1:p [x,y]=find(a(w:p,w)==max(max(a(w:p,w)))); q=a(w,:); a(w,:)=a(x,:); a(x,:)=q; end for j=1:(p-1) for i=(j+1):p a(i,:)=a(j,j)/a(i,j).*a(i,:)-a(j,:); end end m=p; while m1 s(m)=a(m,n); j=p; while ((j2) (j=m+1) (jn)) s(m