7——matlab——利用逆Broyden秩1拟牛顿法求解非线性方程组.doc

下载文档

132
0
约小于1千字
约 3页
2020-09-08 发布于浙江
举报
版权申诉
保障服务

7——matlab——利用逆Broyden秩1拟牛顿法求解非线性方程组.doc

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

利用逆Broyden秩1拟牛顿法求解非线性方程组：函数文件M-文件： function s=dBroyden2(x,eps) if nargin==1 eps=1e-6; else if nargin1 error return; end end [ff dff]=NewtonSfun(x); f1=ff; b1=inv(dff); p=-b1*f1; x3=x+p; [ff dff]=NewtonSfun(x3); d=ff; q=(d-f1); b1=b1+(p-(dli5_19fun(x)\q))*p*b1*inv(p*(dli5_19fun(x)\q)); while (norm(p)eps) [ff dff]=NewtonSfun(x3); f1=ff; p=-b1*f1; x3=x3+p; [ff dff]=NewtonSfun(x3); q=(ff-f1); b1=b1+(p-b1*q)*p*b1*inv(p*(dli5_19fun(x)\q)); end s=x3; return; 求解方程组的函数值和导数值文件：M-文件： function [ff dff]=NewtonSfun(x) syms x1 x2 syms y1 y2 syms y syms y11 y12 syms y21 y22 dy y1=x1*x1-10*x1+x2*x2+8; y2=x1*x2*x2+x1-10*x2+8; y=[y1 y2]; ff=subs(y,[x1,x2],x); y11=diff(y1,x1); y12=diff(y1,x2); y21=diff(y2,x1); y22=diff(y2,x2); dy=[y11,y12;y21,y22]; dff=subs(dy,[x1,x2],x); 调用命令及结果： dBroyden2([0 0]) ans = 1.000000001115433 1.000000001835155