資源描述：
《計(jì)算方法實(shí)驗(yàn)指導(dǎo)書》由會員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在教育資源-天天文庫。
1、計(jì)算方法實(shí)驗(yàn)指導(dǎo)書實(shí)驗(yàn)一：拉格朗日插值算法目的與要求：熟悉拉格朗日插值多項(xiàng)式構(gòu)造類型：必做命名（源程序.cpp及工作區(qū).dsw）：lagrange問題：0.40.550.800.910.410750.578150.888111.026521.17520用四次拉格朗日多項(xiàng)式求的函數(shù)近似值運(yùn)行及結(jié)果顯示：#include#include#defineN4voidmain(){doublex[N+1]={0.4,0.55,0.8,0.9,1};doubley[N+1]={0.41075,0.57815,0.88811,1.02652,1.
2、17520};doublevarx=0.5;intcheckvalid(doublex[],intn);doubleLagrange(doublex[],doubley[],doublevarx,intn);if(checkvalid(x,N)==1)printf("插值結(jié)果： P(%f)=%f",varx,Lagrange(x,y,varx,N));elseprintf("輸入的插值節(jié)點(diǎn)的x值必須互異！");getch();}intcheckvalid(doublex[],intn){…}//檢查插值節(jié)點(diǎn)是否互異，有重節(jié)點(diǎn)返回-1；否則返回1double
3、Lagrange(doublex[],doubley[],doublevarx,intn){…}//參數(shù)說明：x[]，y[]存放插值節(jié)點(diǎn)值，varx變量x的值，n插值多項(xiàng)式的次數(shù)容錯(cuò)測試：對doublex[N+1]={0.4,0.4,0.8,0.9,1};將能提示錯(cuò)誤。實(shí)驗(yàn)二：牛頓插值算法目的與要求：熟悉牛頓插值多項(xiàng)式構(gòu)造，注意與拉氏插值的不同特點(diǎn)16類型：必做命名（源程序.cpp及工作區(qū).dsw）：newton_cz問題：0.40.550.800.910.410750.578150.888111.026521.17520用牛頓插值多項(xiàng)式求運(yùn)行及結(jié)果顯示：#inclu
4、de#include#defineN4staticdoublex[N+1]={0.4,0.55,0.8,0.9,1};staticdoubley[N+1]={0.41075,0.57815,0.88811,1.02652,1.17520};voidmain(){doublevarx,f[N+1][N+1];intcheckvalid(doublex[],intn);voidchashang(doublex[],doubley[],doublef[][N+1]);doublecompvalue(doublet[][N+1],doubl
5、ex[],doublevarx);varx=0.5;if(checkvalid(x,N)==1){chashang(x,y,f);printf("牛頓插值結(jié)果： P(%f)=%f",varx,compvalue(f,x,varx));}else printf("輸入的插值節(jié)點(diǎn)的x值必須互異！");getch();}intcheckvalid(doublex[],intn){}//檢查插值節(jié)點(diǎn)是否互異，有重節(jié)點(diǎn)返回-1；否則返回1voidchashang(doublex[],doubley[],doublef[][N+1]){…}//參數(shù)說明：x[]，y[]
6、存放插值節(jié)點(diǎn)值，varx變量x的值，f[][N+1]返回各階差商值要求：在函數(shù)chashang中完成對各階16差商的計(jì)算及以下三角格式打印.顯示量，本題可采用一維數(shù)組僅保存差商表對角線元素，則算法可進(jìn)一步優(yōu)化，有興趣的同學(xué)可嘗試設(shè)計(jì)算法。doublecompvalue(doublef[][N+1],doublex[],doublevarx){…}.//參數(shù)說明：f[][N+1]傳入各階差商值，varx變量x的值，x[]傳入插值節(jié)點(diǎn)值//返回結(jié)果：插值多項(xiàng)式的值，利用嵌套乘法實(shí)現(xiàn)求值運(yùn)算實(shí)驗(yàn)三：自適應(yīng)梯形求積公式目的與要求：通過實(shí)際計(jì)算體會自適應(yīng)梯形公式的原理命名（源程
7、序.cpp及工作區(qū).dsw）：autotrap注意事項(xiàng)：對于中間計(jì)算結(jié)果不得采用數(shù)組來保存，否則空間浪費(fèi)太嚴(yán)重；類型：必做問題：計(jì)算的近似值，使得誤差（EPS）不超過運(yùn)行及結(jié)果顯示：#include#include#include#defineEPS1e-6voidmain(){doubles;doublef(double);doubleAutoTrap(double(*)(double),double,double,double);s=AutoTrap(f,0.0,1.0,EPS);getch()
