資源描述:
《數(shù)值分析計(jì)算方法》由會(huì)員上傳分享,免費(fèi)在線閱讀,更多相關(guān)內(nèi)容在行業(yè)資料-天天文庫(kù)。
1、word格式《計(jì)算方法》實(shí)驗(yàn)內(nèi)容一.實(shí)驗(yàn)一:用兩種不同的順序計(jì)算,分析其誤差的變化。1.實(shí)驗(yàn)?zāi)康模和ㄟ^(guò)正序反序兩種不同的順序求和,比較不同算法的誤差;了解在計(jì)算機(jī)中大數(shù)吃小數(shù)的現(xiàn)象,以后盡量避免;體會(huì)單精度和雙精度數(shù)據(jù)的差別。2.算法描述:累加和s=0;正序求和:對(duì)于n=1,2,3,......,10000s+=1.0/(n*n);反序求和:對(duì)于n=10000,9999,9998,.....,1s+=1.0/(n*n);3.源程序:#雙精度型##includecvoidmain(){doubles=0;intn;for(n=1;n<=1000
2、0;n++)s+=1.0/(n*n);printf("正序求和結(jié)果是:%lf",s);s=0;for(n=10000;n>=1;n--)s+=1.0/(n*n);printf("反序求和結(jié)果是:%lf",s);}#單精度型##includevoidmain(){floats=0;intn;for(n=1;n<=10000;n++)s+=1.0/(n*n);printf("正序求和結(jié)果是:%f",s);s=0;for(n=10000;n>=1;n--)s+=1.0/(n*n);printf("反序求和結(jié)果是:%f",s);}..
3、..word格式1.運(yùn)行結(jié)果:雙精度型運(yùn)行結(jié)果:?jiǎn)尉刃瓦\(yùn)行結(jié)果:2.對(duì)算法的理解與分析:舍入誤差在計(jì)算機(jī)中會(huì)引起熟知的不穩(wěn)定,算法不同,肯結(jié)果也會(huì)不同,因此選取穩(wěn)定的算法很重要。選取雙精度型數(shù)據(jù)正反序求和時(shí)結(jié)果一致,但選用單精度型數(shù)據(jù)時(shí),求和結(jié)果不一致,明顯正序求和結(jié)果有誤差,所以第一個(gè)算法較為穩(wěn)定可靠。二.實(shí)驗(yàn)二:1、拉格朗日插值按下列數(shù)據(jù)x-3.0-1.01.02.03.0y1.01.52.02.01.0作二次插值,并求x=-2,x=0,x=2.75時(shí)的函數(shù)近似值2牛頓插值按下列數(shù)據(jù)x0.300.420.500.580.660.72y1.044031.0
4、84621.118031.156031.198171.23223作五次插值,并求x=0.46,x=0.55,x=0.60時(shí)的函數(shù)近似值.1.實(shí)驗(yàn)?zāi)康模和ㄟ^(guò)拉格朗日插值和牛頓插值的實(shí)例,了解兩種求解方法,并分析各自的優(yōu)缺點(diǎn)。2.算法描述:3.源程序:拉格朗日插值:#include#definek2voidmain()....word格式{doubleL,Y,a;doublex[3];doubley[3];for(intp=0;p<=2;p++)scanf("%lf%lf",&x[p],&y[p]);for(intq=0;q<=2;q++){Y=
5、0;scanf("%lf",&a);for(intj=0;j<=k;j++){L=1;for(inti=0;i<=k;i++)if(i!=j)L*=(a-x[i])/(x[j]-x[i]);Y+=L*y[j];}printf("x=%lf的函數(shù)值是:%lf",a,Y);}}牛頓插值:#includevoidmain(){inti,j;doublea[6][6],y,s,x;for(j=0;j<=1;j++)for(i=0;i<6;i++)scanf("%lf",&a[i][j]);for(j=2;j<=5;j++)for(i=j-1;i<
6、=5;i++)a[i][j]=(a[i][j-1]-a[i-1][j-1])/(a[i][0]-a[i-j+1][0]);for(intk=0;k<3;k++){scanf("%lf",&x);y=0;for(i=0;i<5;i++){s=a[i][i+1];for(j=i-1;j>=0;j--)s*=x-a[j][0];y+=s;}....word格式printf("結(jié)果是:%lf",y);}4.運(yùn)行結(jié)果:拉格朗日運(yùn)行結(jié)果:牛頓插值運(yùn)行結(jié)果:5.對(duì)算法的理解與分析:(1)拉格朗日插值:該公式是對(duì)一系列點(diǎn)加權(quán)求和。內(nèi)插通常優(yōu)于外推。選擇的區(qū)間包括x,插值效
7、果越好,高次插值通常優(yōu)于低次插值,但并不是意味著插值次數(shù)越高越好。優(yōu)點(diǎn):編程容易實(shí)現(xiàn)。缺點(diǎn):如果發(fā)現(xiàn)當(dāng)前的插值方法不夠精確,就要增加插值節(jié)點(diǎn)的個(gè)數(shù)。拉格朗日基函數(shù)每一個(gè)都要重新計(jì)算,效率低。(2)牛頓插值:優(yōu)點(diǎn):如果當(dāng)前插值方法不穩(wěn)定,需要增加插值節(jié)點(diǎn)個(gè)數(shù),只需要計(jì)算新家節(jié)點(diǎn)所增加的差商即可,之前的計(jì)算結(jié)果可以被復(fù)用,比拉格朗日插值節(jié)省計(jì)算量,效率高,精度高,更為穩(wěn)定。....word格式三.實(shí)驗(yàn)三:分別用復(fù)化梯形公式和復(fù)化辛卜生公式計(jì)算f(x)=sin(x)/x的積分,并與準(zhǔn)確值比較判斷精度。1.實(shí)驗(yàn)?zāi)康模和ㄟ^(guò)實(shí)例體會(huì)各種算法的精度。熟練掌握復(fù)化梯形,復(fù)化辛
8、普森,復(fù)化柯特斯求積方法的程序。2.算