c數(shù)值分析上機(jī)實(shí)驗(yàn)

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

1、第一早舍入誤差和數(shù)值穩(wěn)定性1:例題利用遞推公式計(jì)算定積分1■■0xndx方法一：利用遞歸公式兒=；_5兒一1(門=0/1,2,3?……20)11方法二利用遞推公式=-^--yn(n=0,l,2,3?……20)這兩種方法的思路是：1：逆代初值2：根據(jù)范圍來確定循環(huán)次數(shù)3：根據(jù)逆代函數(shù)確定逆代形似4：如何排版1.1：方法一的C語言實(shí)現(xiàn)：#inelude#includefloaty0=log(6.0)-log(5.0),yl;intn二1;printf(,,y[0]=%-20f"/y0);while(n<=20){yl

2、=1.0/n-5*y0;printf(,Iy[%d]=%-20f,,,n,yl);yO=yl;n++;if(n%3==0)printf(,,");}}5，[0]=0.182322y[31=0.043138^61=0.024428y[91=0.004076y[121=1.619237&[151=-200.772052kl8J=25097.828993y[ll=0.088392y[41=0.034310y[71=0.020719y[10]=0.079618y[13l=-8.019263y[161=1003.922745y[19]=-125489.0

3、87993y[21=0.058039y[51=0.028448y[81=0.021407y[lll=-0.307181y[141=40.167745y[l?]=-5019.554824y[201=627445.4796881.2：方法二的c語言實(shí)現(xiàn):#includefloaty20二(1/105.0+1/126.0)/20)19;intn二20;printf(,,y[20]=%-20f,,,y20);while(n>0){yl9=1.0/(5*n)-y20/5;printf(,,y[%d]=%-20f,,,n-l,yl9);y20

4、=yl9;n—;if(n%3==0)printf(,,");}9[201=0.000873y[171=0.009399y[141=0.011229y[111=0.014071y[8J=0.01883?9(51=0.028468y[2J=0.058039y[191=0.009825^[161=0.009885p[131=0.012040p[10]=0.015368y[?]=0.021233y[4J=0.034306y[ll=0.088392y[18]=0.008561Sf[151=0.010523y[121=0.0129??y[9]=0.0169

5、26y[61=0.024325y[31=0.043139y[0]=0.182322通過兩種計(jì)算方式可以看出第二種的算法收斂2：課后習(xí)題13根據(jù)pn=(10^)*pn.1-pn.2可知道pn的相對誤差限為en=(10^)*en.ren_2n二23…?…根據(jù)書中的例題方法很難比較誤差是否增大，現(xiàn)在經(jīng)過以下遞歸程序可以看出其結(jié)果：#includedoublefact(intn);voidmain(){intn;doubleresult;printf("lnputn:");scanf("%d",&n);result=fact(n);訐(r

6、esult==-l)printf("lnputerror");/*n的值輸入不正確*/elseprintf("wc=%.5f",result);}doublefact(intn){訐(*2)return-1;elseif(n==2)/*終止條件*/return0.005;/*起始誤差為0.005*/elsereturn(10.0/3)*fact(n-l)-fact(n-2);/*調(diào)用遞歸*/Inputn：4wc=3.38389Pressanykeytocontinue結(jié)果表明當(dāng)n=4時(shí)誤羌要大于nJ時(shí)的誤羌，所以不穩(wěn)定。17：編寫程序如下#in

7、cludevoidmain(){floati,j,n;floatsuml=0;floatsum2=0;floatzhi,wclzwc2;for(i=2;i<=100;i++){suml+=l/(i*i-l)；}printf("suml=%fu,suml);for(j=100;j>=2;j-){sum2+=l/(j*j-l)；}printf("sum2=%f"/sum2);zhi=l.O/2*(3.O/2-(l.O/n)-l.O/(n+l));printf(,,zhi=%f',,zhi);wcl=zhi-suml;wc2

8、=zhi-sum2;printf(,,wcl=%f"/wcl);printf(,,wc2=%f",wc2);}三種