數(shù)值分析上機(jī)

《數(shù)值分析上機(jī)》由會(huì)員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在工程資料-天天文庫(kù)。

1、1＞分別用不動(dòng)點(diǎn)迭代與Newton法求解方程x-e2x+5=0的正根與負(fù)根。解答：(1)開(kāi)始把原方程變形為x二"2x+5;用不動(dòng)點(diǎn)迭代法發(fā)現(xiàn)不收斂pO=2.389056098930650p0二1.138797166964409e+02pO二8.216088789469673e+98后來(lái)變形為x二ln(x+5)/2i=l;N=20;T=0?0001;p0=l;disp(，不動(dòng)點(diǎn)迭代法求正根,);fori=l:Np=(log(pO+5))/2;if(abs(p-pO)

2、endzuoyel不動(dòng)點(diǎn)迭代法求正根0.886320947105208.⑵對(duì)x-e2x+5=0求導(dǎo)yJ：L?2"2x;i=l;N=20;T=0.0001;p0=-l;formatlongdisp(*牛頓法求負(fù)根T;fori=l:Np=p0-(pO-exp(2*p0)+5)/(l-2*exp(2*p0));if(abs(p-pO)

3、0/使用下面的各種方法求位^[0.1,1]內(nèi)的解，精確到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：應(yīng)用Newton法求f(x)的零點(diǎn)，e=10A-6這里f(x)=x-sinx再用求重根的兩種方法求f(x)的零點(diǎn)。i=l;N=100;T=0.000001;p0=l;disp(fA￡3[U-;for