newton法快速弦截法的实验报告.doc

下载文档 降价啦

40
0
约3.24千字
约 5页
2016-12-31 发布于未知
举报
版权申诉
保障服务

newton法快速弦截法的实验报告.doc

1、本文档共5页，可阅读全部内容。
2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。

Newton 快速弦截法的实验报告一：实验题目：newton法，快速弦截法。二：实验目的：通过用MATLAB语言计算newton法，弦截法对两种方法的结果进行分析。三：实验内容：a：了解MATLAB语言的用法b：用这两种方法求线性方程组四：实验方法与步骤： Newton法 Function [x,k]=Mendnewton(f,x0,emg) [f1,d1]=feval(f,x0); k=1; x(1)=x0; x(2)=x(1)-f1/d1; while abs(f1)emg u=1; k=k+1; [f1,d1]=feval(f,x(k)); X(k+1)=x(k)-u*f1/d1; [f2,d2]=feval(f,x(k+1)); While abs(f2)abs(f1) U=u/2; X(k+1)=x(k)-u*f1/d1; [f2,d2]=feval(f,x(k+1)); end end 例：用牛顿下山法求方程f(x)= 的根，使精度达到10 初值分别选取为：（1）x0=-1;(2)x0=2.0; 编写函数名func3.m Function [f,d]=func3(x) F=sqrt(x^2+1)-tan(x); D1=’sqrt(x^2+1)-tan(x)’; D=subs(diff(d1)); 在命令窗口输入：f=@func3;[x,k]=Mendnewton(f,x0,10^-6);若选初值为x0=-1;运行结果为：迭代次数k x值 f(1)=feval(f,x(k))值 1 -7.069047932971935e-001 2.078789280010764e+000 2 1.942400972108479e-001 8.219695728408301e-001 3 1.163518073303871e+000 -7.838374932606165e-001 4 1.023918977930554e+000 -2.113030290935414e-001 5 9.530711345686330-001 -2.606996588743926e-002 6 9.416925081385333e-001 -5.085688425774393e-004 7 9.414616152761416e-001 -2.012113482496858e-007 8 9.414615238528302e-001 -3.153033389935445e-014 若选初值为x0=2.0;运行结果如下：迭代次数k x值 f1(kfeval(f,x(k))值 1 2.905969917234289e+000 3.313289980588459e+000 2 3.829942435553551e+000 3.135774879468060e+000 3 4.382754035040099e+000 1.572413838570136e+000 4 4.474505813415593e+000 4.607393526441488e-001 5 4.501556126032599e+000 -6.131570643391715e-002 6 4.498750820792893e+000 -8.285614610334946e-004 7 4.498711866735406e+000 -1.555631969907267e-007 8 4.498711859418998e+000 -1.154631945610163e-014 快速弦截法： function [x,k]=Fast_chord(f,x1,x2,emg); k=1; y1=feval(f,x1); y2=feval(f,x2); x(k)=x2-(x2