實(shí)驗(yàn)1拉格朗日插值與牛頓插值

《實(shí)驗(yàn)1拉格朗日插值與牛頓插值》由會(huì)員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在教育資源-天天文庫。

1、西華大學(xué)計(jì)算機(jī)系學(xué)生上機(jī)實(shí)踐報(bào)告西華數(shù)學(xué)與計(jì)算機(jī)學(xué)院上機(jī)實(shí)踐報(bào)告課程名稱：計(jì)算方法A年級(jí)：2010級(jí)上機(jī)實(shí)踐成績(jī)：指導(dǎo)教師：嚴(yán)常龍姓名：上機(jī)實(shí)踐名稱：拉格朗日插值和牛頓插值法學(xué)號(hào)：上機(jī)實(shí)踐日期：上機(jī)實(shí)踐編號(hào)：1上機(jī)實(shí)踐時(shí)間：一、目的1．通過本實(shí)驗(yàn)加深對(duì)拉格朗日插值和牛頓插值法構(gòu)造過程的理解；2．能對(duì)上述兩種插值法提出正確的算法描述編程實(shí)現(xiàn)。二、內(nèi)容與設(shè)計(jì)思想自選插值問題，編制一個(gè)程序，分別用拉格朗日插值法和牛頓插值法求解某點(diǎn)的函數(shù)近似值。（從課件或教材習(xí)題中選題）已知y=f(x)的數(shù)據(jù)表如下，求t=0.63處的函數(shù)值z(mì)=f(t)。ixiyiixiyi10.100.90483760.570.56

2、552520.150.86070870.700.49658530.250.77880180.850.42741540.400.67032090.930.39455450.500.606531101.000.367879三、使用環(huán)境操作系統(tǒng)：軟件環(huán)境：四、核心代碼及調(diào)試過程1.拉格朗日插值代碼#includemain(){floatx[8],y[8];floatt,z=0.0,s;intn=0,i=0,tp=0,start,end,j;printf("請(qǐng)輸入節(jié)點(diǎn)控制數(shù)n");printf("");scanf("%d",&n);printf("請(qǐng)輸入節(jié)點(diǎn)t的值");print

3、f("");第5頁共5頁西華大學(xué)計(jì)算機(jī)系學(xué)生上機(jī)實(shí)踐報(bào)告scanf("%f",&t);//節(jié)點(diǎn)賦值for(i=0;in-1){end=n-1;start=end-7;}//用所選插值點(diǎn)做拉格朗日插

4、值多項(xiàng)式運(yùn)算for(i=start;i#definemax20main(){doublex[max],y[max];//插值點(diǎn)xi以及對(duì)應(yīng)函數(shù)yi的函數(shù)值doublet,fw;//t為待求的參數(shù)點(diǎn)fw為對(duì)應(yīng)要求的值doubleft[max

5、][max];//表示各階差商第二維為階數(shù)intn=0,i=0,j=0;doubleb[max];/*printf("請(qǐng)輸入節(jié)點(diǎn)控制數(shù)n");printf("");scanf("%d",&n);printf("請(qǐng)輸入節(jié)點(diǎn)t的值");printf("");scanf("%f",&t);//節(jié)點(diǎn)賦值for(i=0;i

6、;y[1]=0.860708;x[2]=0.25;y[2]=0.778801;x[3]=0.40;y[3]=0.670320;x[4]=0.50;y[4]=0.606531;x[5]=0.57;y[5]=0.565525;x[6]=0.70;y[6]=0.496585;x[7]=0.85;y[7]=0.427415;第5頁共5頁西華大學(xué)計(jì)算機(jī)系學(xué)生上機(jī)實(shí)踐報(bào)告x[8]=0.93;y[8]=0.394554;x[9]=1.00;y[9]=0.367879;for(i=0;i

7、t[i][0]=y[i];//計(jì)算0階差商}for(i=1;i=1;i--){b[i-1]=ft[i-1][i-1]+b[i]*(t-x[i-1]);}fw=b