Newton-Raphson 函数.docxVIP

下载本文档

5
0
约1.84千字
约 5页
2018-02-15 发布于河南
举报
版权申诉

Newton-Raphson 函数.docx

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

Newton-Raphson 函数

function hom[P,iter,err]=newton(‘f’,’JF’,[7.8e-001;4.9e-001;3.7e-001],0.01,0.001,1000);disp(P);disp(iter);disp(err);function Y=f(x,y,z)Y=[x^2+y^2+z^2-1;2*x^2+y^2-4*z;3*x^2-4*y-z^2];function y=JF(x,y,z)f1=x^2+ y^2+z^2-1;f2= 2*x^2+y^2-4*z;f3=2*x^2+y^2-4*z;df1x=diff(sym(f1),x)df1y=diff(sym(f1),y)df1z=diff(sym(f1),z)df2x=diff(sym(f2),x);df2y=diff(sym(f2),y);df2z=diff(sym(f2),z);df3x=diff(sym(f3),x);df3y=diff(sym(f3),y);df3z=diff(sym(f3),z);j=[df1x,df1y,df1z; df2x,df2y,df2z; df3x,df3y,df3z;];y=(j);function [P,iter,err]=newton(F,JF,P,tolp,tolfp,max)%输入P为初始猜测值，输出P则为近似解%JF为相应的Jacobian矩阵%tolp为P的允许误差%tolfp为f(P)的允许误差%max:循环次数Y=f(F,P(1), P(2), P(3));for k=1:maxJ=f(JF, P(1), P(2), P(3));Q=P-inv(J)*Y;Z=f(F,Q(1), Q(2), Q(3));err=norm(Q-P);P=Q;Y=Z;Iter=k;If (errtolp)||(abs(Y)tolfp||abs(Y)0.0001)breakendendpre lang=”matlab” line=”1” file=”test.m”function homework4[P,iter,err]=newton(‘f’,’JF’,[7.8e-001;4.9e-001;3.7e-001],0.01,0.001,1000);disp(P);disp(iter);disp(err);function Y=f(x,y,z)Y=[x^2+y^2+z^2-1;2*x^2+y^2-4*z;3*x^2-4*y=z^2];function y=JF(x,y,z)f1=‘x^2+ y^2+z^2-1’;f2=‘ 2*x^2+y^2-4*z’;f3=‘ 2*x^2+y^2-4*z’;df1x=diff(sym(f1),’x’);df1y=diff(sym(f1),’y’);df1z=diff(sym(f1),’z’);df1x=diff(sym(f1),’x’);df2x=diff(sym(f2),’x’);df2y=diff(sym(f2),’y’);df2z=diff(sym(f2),’z’);df3x=diff(sym(f3),’x’);df3y=diff(sym(f3),’y’);df3z=diff(sym(f3),’z’);j=[df1x,df1y,df1z; df2x,df2y,df2z; df3x,df3y,df3z;];y=(j);function [P,iter,err]=newton(F,JF,P,tolp,tolfp,max)%输入P为初始猜测值，输出P则为近似解%JF为相应的Jacobian矩阵%tolp为P的允许误差%tolfp为f(P)的允许误差%max:循环次数Y=f(F,P(1), P(2), P(3));For k=1:maxJ=f(JF, P(1), P(2), P(3));Q=P-inv(J)*Y;Z=f(F,Q(1), Q(2), Q(3));err=norm(Q-P);P=Q;Y=Z;Iter=k;If (errtolp)||(abs(Y)tolfp||abs(Y)0.0001)breakend牛顿下山法引入因子t ，建立迭代法在迭代时，选择一个t使得T的选取方式按t=1,1/2;1/2^2;1/2^3,…的顺序7eval 和 feval 函数y=sym(‘x’);f=y^2+5*y;x=3;eval(f)