数值分析上机

《数值分析上机》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、1＞分别用不动点迭代与Newton法求解方程x-e2x+5=0的正根与负根。解答：(1)开始把原方程变形为x二"2x+5;用不动点迭代法发现不收敛pO=2.389056098930650p0二1.138797166964409e+02pO二8.216088789469673e+98后来变形为x二ln(x+5)/2i=l;N=20;T=0・0001;p0=l;disp(，不动点迭代法求正根,);fori=l:Np=(log(pO+5))/2;if(abs(p-pO)

2、endzuoyel不动点迭代法求正根0.886320947105208.⑵对x-e2x+5=0求导yJ：L・2"2x;i=l;N=20;T=0.0001;p0=-l;formatlongdisp(*牛顿法求负根T;fori=l:Np=p0-(pO-exp(2*p0)+5)/(l-2*exp(2*p0));if(abs(p-pO)

3、0/使用下面的各种方法求位^[0.1,1]内的解，精确到10A-4.a•二分法i=l;T=0.0001;a=0=l;N=20;disp('二分法T)f=@(x)600*xA4-550*xA3+200*xA2-20*x-l;FA=f(a);fori=l:Np=(a+b)/2;FP=f(p)；if(FP==0

4、

5、(b-a)/2

6、30;T=0.0001;p0=l;disp(■牛顿法求根T;formatlongf=@(x)600*xA4-550*xA3+200*xA2-20*x-l;s=@(x)2400*xA3-1650*xA2+400*x-20;fori=l:NP=pO-f(pO)/s(pO);if(abs(p-pO)

7、ongf=@(x)600*xA4-550*xA3+200*xA2-20*x-l;qO=f(pO);ql=f(pl)；fori=l:Np=pl-ql*(pl-pO)/(ql-qO)；if(abs(p-pl)

8、3+200*xA2-20*x-l;qO=f(pO);ql=f(pl)；fori=l:Np=pl-ql*(pl-pO)/(ql-qO)；if(abs(p-pl)

9、0*xA4-550*xA3+200*xA2-20*x-l;hl=pl-pO;h2=p2-pl;gl=(f(pl)-f(pO))/hl;g2=(f(p2)-f(pl))/h2;d=(g2-gl)/(h2+hl);fori=l:Nb=g2+h2*d;D=sqrt((匕入2-4*£(p2)*d));if(abs(b-D)(b+D))E=b+D;elseE=b_D;endh=-2*f(p2)/E;p=p2+h;if(h

10、)-f(pO))/hl;g2=(f(p2)-f(pl))/h2;d=(g2-gl)/(h2+hl);endzuoye2emuller法求根0.2359990704136793：应用Newton法求f(x)的零点，e=10A-6这里f(x)=x-sinx再用求重根的两种方法求f(x)的零点。i=l;N=100;T=0.000001;p0=l;disp(fA£3[U-;for