數(shù)值分析計(jì)算方法

數(shù)值分析計(jì)算方法

ID:47688292

大小:104.50 KB

頁(yè)數(shù):11頁(yè)

時(shí)間:2020-01-24

數(shù)值分析計(jì)算方法_第1頁(yè)
數(shù)值分析計(jì)算方法_第2頁(yè)
數(shù)值分析計(jì)算方法_第3頁(yè)
數(shù)值分析計(jì)算方法_第4頁(yè)
數(shù)值分析計(jì)算方法_第5頁(yè)
資源描述:

《數(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.算

當(dāng)前文檔最多預(yù)覽五頁(yè),下載文檔查看全文

此文檔下載收益歸作者所有

當(dāng)前文檔最多預(yù)覽五頁(yè),下載文檔查看全文
溫馨提示:
1. 部分包含數(shù)學(xué)公式或PPT動(dòng)畫的文件,查看預(yù)覽時(shí)可能會(huì)顯示錯(cuò)亂或異常,文件下載后無(wú)此問(wèn)題,請(qǐng)放心下載。
2. 本文檔由用戶上傳,版權(quán)歸屬用戶,天天文庫(kù)負(fù)責(zé)整理代發(fā)布。如果您對(duì)本文檔版權(quán)有爭(zhēng)議請(qǐng)及時(shí)聯(lián)系客服。
3. 下載前請(qǐng)仔細(xì)閱讀文檔內(nèi)容,確認(rèn)文檔內(nèi)容符合您的需求后進(jìn)行下載,若出現(xiàn)內(nèi)容與標(biāo)題不符可向本站投訴處理。
4. 下載文檔時(shí)可能由于網(wǎng)絡(luò)波動(dòng)等原因無(wú)法下載或下載錯(cuò)誤,付費(fèi)完成后未能成功下載的用戶請(qǐng)聯(lián)系客服處理。