分段二次插值_第1頁
分段二次插值_第2頁
分段二次插值_第3頁
分段二次插值_第4頁
全文預(yù)覽已結(jié)束

下載本文檔

版權(quán)說明：本文檔由用戶提供并上傳，收益歸屬內(nèi)容提供方，若內(nèi)容存在侵權(quán)，請進行舉報或認領(lǐng)

文檔簡介

1、南昌航空大學(xué)實驗報告 2010.12.06課程名稱：計算方法 實驗名稱：用分段二次插值計算班級：09045131 姓名：楊 望一 、 實驗?zāi)康模?（1）學(xué)會使用分段二次插值方法進行計算。 （2）掌握在程序設(shè)計中用循環(huán)的方法實現(xiàn)分段二次插值算法。 (3) 進一步熟練使用C語言編程。二 、實驗內(nèi)容：從函數(shù)表x0.0 0.1 0.195 0.3 0.401 0.5 f(x) 0.39894 0.39695 0.39142 0.38138 0.36812 0.35206出發(fā)，試使用分段二次插值方法分別計算f(0.15),f(0.31),及f(0.47)的近似值。三 、 實驗步驟：按照流程圖很容易可以實

2、現(xiàn)實驗?zāi)康?。實驗流程圖如下圖所示：四 、 實驗結(jié)果：五、程序代碼：#include"stdio.h"#define N 100void main()int n,j,i,I,k;float xN,yN,a,p,Y=0.0,t;printf("please input nn");scanf("%d",&n);printf("please input xn");for(i=0;i<n;i+) scanf("%f",&xi);printf("please input yn

3、");for(i=0;i<n;i+) scanf("%f",&yi);printf("please input an");scanf("%f",&a);for(j=1;j<n-1;j+) p=0.5*(xj+xj+1); if(a<p) I=j; for(k=I-1;k<I+2;k+) for(j=I-1,t=1.0;j<I+2;j+) if(j!=k)t*=(a-xj)/(xk-xj); Y+=yk*t; printf("%f%f",a,Y);break; if(a>0.5*(xn-3+xn-2) I=n-2;for(k=I-1;k<I+2;k+) for(j=I-1,t=1.0;j<I+2;j+) if(j!=k) t*=(a

溫馨提示

  • 1. 本站所有資源如無特殊說明，都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等，如果需要附件，請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙，網(wǎng)頁內(nèi)容里面會有圖紙預(yù)覽，若沒有圖紙預(yù)覽就沒有圖紙。
  • 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
  • 5. 人人文庫網(wǎng)僅提供信息存儲空間，僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護處理，對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯，并不能對任何下載內(nèi)容負責(zé)。
  • 6. 下載文件中如有侵權(quán)或不適當內(nèi)容，請與我們聯(lián)系，我們立即糾正。
  • 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

評論

0/150

提交評論