数值分析matlab编程之非线性方程求根.docx

下载文档 降价啦

2
0
约1.28千字
约 4页
2017-09-05 发布于浙江
举报
版权申诉
保障服务

数值分析matlab编程之非线性方程求根.docx

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

数值分析matlab编程之非线性方程求根

1.二分法function [ x,k ] = bisect( a,b,f,emg )digits(8);if nargin 4 emg=1e-5;endmax1=round((log(b-a)-log(emg))/log(2));N=zeros(max1,2);k=0;while abs(b-a)emg fa=feval(f,a); fab=feval(f,(a+b)/2);if fab==0 x=(a+b)/2;return;elseif fa*fab0 b=(a+b)/2;else a=(a+b)/2;end N(k+1,1)=a; N(k+1,2)=b; k=k+1;endx=vpa((a+b)/2);vpa(N)2.牛顿法function [ x ] = Newton( fun2,x0,err ) digits(10);[f0,df]=feval(fun2,x0);if nargin3; err=1e-6;endif df==0; error(df=0 无解);endd=f0/df;i=1;while abs(d)err; x1=x0-d; [f0,df]=feval(fun2,x1);if df==0; error(df=0 无解); end x0=x1; d=f0/df; fprintf(x(%d)=%.10f err(%d)=%.10f\n,i,x1,i,d); i=i+1;endx=x1;3.不动点迭代function [ p0,i,err,p] = NMPiteration( f,p0,tol,max1 )P(1)=p0;for i=1:max1 P(i+1)=feval(f,P(i));%format compact err=abs(P(i+1)-P(i)); fprintf(err(%d)=%f\n,i,err); p=P(i);if(errtol)break;endif i==max1 disp(‘达到最大迭代次数，退出迭代);endendx=P4.Steffensen迭代function [ p0,i,err,p] = Steffensen( f,p0,tol,max1 )P(1)=p0;for i=1:max1 y=feval(f,P(i)); z=feval(f,y); P(i+1)=P(i)-(y-P(i)).^2/(z-2*y+P(i));%format compact err=abs(P(i+1)-P(i)); fprintf(err(%d)=%f\n,i,err); p=P(i);if(errtol)break;endif i==max1 disp(达到最大迭代次数，退出迭代);endendx=P