資源描述:
《計(jì)算方法源程序》由會(huì)員上傳分享,免費(fèi)在線閱讀,更多相關(guān)內(nèi)容在工程資料-天天文庫(kù)。
1、計(jì)算方法實(shí)驗(yàn)報(bào)告實(shí)驗(yàn)一舍入誤差與數(shù)值穩(wěn)定性班級(jí):計(jì)算機(jī)cO96姓名:郭睿學(xué)號(hào):088115實(shí)驗(yàn)一舍入誤差與數(shù)值穩(wěn)定性目的與要求:1、通過上機(jī)編程,復(fù)習(xí)鞏固以前所學(xué)程序設(shè)計(jì)語(yǔ)言;2、通過上機(jī)計(jì)算,了解舍入誤差所引起的數(shù)值不穩(wěn)定性。3、通過上機(jī)計(jì)算,了解運(yùn)算次序?qū)τ?jì)算結(jié)果的影響,從而盡量避免大數(shù)吃小數(shù)的現(xiàn)象。實(shí)驗(yàn)內(nèi)容:?通過止反兩個(gè)實(shí)例的計(jì)算,了解利用計(jì)算機(jī)進(jìn)行數(shù)值計(jì)算中舍入誤并所引起的數(shù)值不穩(wěn)定性,深入理解初始小的舍入誤差可能造成誤差積累從而對(duì)計(jì)算結(jié)果的巨大影響。?通過實(shí)際編程,了解運(yùn)算次序?qū)τ?jì)算結(jié)果的影響,了解實(shí)數(shù)運(yùn)算符合的結(jié)合律和分配律在計(jì)算機(jī)里不一定
2、成立。1對(duì)n=0,1,2,...,20計(jì)算定積分算法1利用遞推公式nln6-ln5?0.182322算法2利用遞推公式人n=20,19,...,l注意到110511112011nr訂x%*占注和x%算法一:#include#include#include#includeHstdlib.hHvoidmain(void){floaty_0=log(6.0)-log(5.0),y_l;intn=l;system(”cis”);/*清屏*/printf(”y[0]=%?20f”,y_0);while(l){y_
3、l=1.0/n-5*y_0;printf(Ky[%d]=%-20f,y_l);/*輸出*/if(n>=20)break;y_0=y_l;n++;if(n%3==0)printf(Hn);}getch();/*保持用戶屏幕*/運(yùn)行結(jié)杲:Y[0]=0.182322Y[l]=0.088392Y⑵=0.058039Y[3]=0.043128Y[41=0.034310Y[5]=0.028448Y[6]=0.024428Y[7]=0.020719Y[8]=0.021407Y[9]=0.004076Y[10]=0.079618Y[ll]=-0.307181Y[
4、12]=1.619237Y[13]=-8.019263Y[14]=40.167744Y[15]=-200.772069Y[16]=1003.922729Y[17]=-5019.554688Y[18)=25097.828125Y[19]=-125489.085938Y(20
5、=627445.500000算法二#include#include#include#include"stdlib.hHmain()floaty_0=(l/105.0+l/l26.0)/2,y_l;intn=20;system(ncls
6、");/*清屏*/printf(Hy[201=%-20fy_0);while(l){y_l=l/(5.0*n)-y_0/5.0;printf(Ky[%d]=%-20f-1,y_l);if(n<=l)break;y_O=y_l;n-;if(n%3==0)printf(nn);}getch();運(yùn)行結(jié)杲:Y[20]=0.008730Y[19]=0.008254Y[18]=0.008876Y[17]=0.009336Y[16]=0.009898Y[151=0.010520Y[14]=0.011229Y[13]=0.012040Y[12]=0.0129
7、77Y[ll]=0.014071Y[10]=0.015368Y[9]=0.016926Y⑻=0.018837Y[7]=0.021233Y[6]=0.024325Y[5]=0.028468Y[4]=0.034306Y[3]=0.043139Y(2]=0.058039Y(l]=0.088392Y[0
8、=0.182322說明:從計(jì)算結(jié)果可以看出,算法1是不穩(wěn)定的,而算法2是穩(wěn)定的。計(jì)算方法實(shí)驗(yàn)報(bào)告實(shí)驗(yàn)二拉格朗日插值與牛頓插值班級(jí):計(jì)算機(jī)c096姓名:學(xué)號(hào):088115實(shí)驗(yàn)二拉格朗日插值與牛頓插值一、目的與要求:熟悉拉榕朗日插值多項(xiàng)式和牛頓插值多項(xiàng)式,注意其不同
9、特點(diǎn);二、實(shí)驗(yàn)內(nèi)容:通過拉格朗口插值和牛頓插值多項(xiàng)式的兩個(gè)實(shí)例的計(jì)算,了解兩種求解方法,分析他們的優(yōu)缺點(diǎn)。三、程序與實(shí)例算法1.輸入X,.,y;.(i=0,l,2,...,n),令L(xM)=0;2.對(duì)=0,1,2,.,?n計(jì)算?1/W=nn/=0j右X-Xjxi-XjL,kL“+l,(x)y,程序與實(shí)例例1已知函數(shù)農(nóng)Xj0.561600.562800.564010.56521y,-0.827410.826590.825770.82495用三次拉格朗日多項(xiàng)式求x=0.5635的函數(shù)近似值。拉格朗日插值:#include#include<
10、conio.h>#include#include"stdli