實驗一題目一、已知函式表如下:
試用三次lagrange插值多項式求x=0.5635時的函式值。
解:lagrange插值法公式:
程式框圖如下:
本題中n=4
分析設計思想:
設計思路如上面的流程圖所示,用我們學過的c語言程式設計實現。
程式清單:
#include
#include
#include
float lagrange(float *x,float *y,float xx,int n) /*拉格朗日插值演算法*/
free(a);
return yy;
}main()
printf("\n");
for(i=0;i<=n-1;i++)
printf("\n");
printf("input xx:");
scanf("%f",&xx);
yy=lagrange(x,y,xx,n);
printf("x=%f,y=%f\n",xx,yy);
getch();
}程式執行結果:
按照試驗題目中要求輸入xi和yi的資料,得到如下結果:
經驗證,執行結果正確。
除錯過程分析:
需要使用陣列語句定義xi和yi,
當要定義素組的長度時,標頭檔案中要有
#include
題目二:已知函式表如下:
造出差商表,並用newton插值多項式求nn(0.596)和nn(0.895)。
解:程式流程圖如下:
分析設計思想:
程式清單:
#include
#include
#include
void difference(float *x,float *y,int n)
return;
}main()
printf("\n");
for(i=0;i<=n-1;i++)
printf("\n");
difference(x,(float *)y,n);
printf("input xx:");
scanf("%f",&xx);
yy=y[20];
for(i=n-1;i>=0;i--) yy=yy*(xx-x[i])+y[i];
printf("newtoninter(%f)=%f",xx,yy);
getch();
}程式執行結果:
經驗證執行結果正確。
數值分析實驗報告
數值分析 課程設計實驗報告模板 常微分方程資料值解 用龍格 庫塔法分析lorenz方程的特性 一 考慮著名的lorenz方程 其中s,r,b為變化區域內有一定限制的實引數,該方程形式簡單,表面上看並無驚人之處,但由該方程揭示出的許多現象,促使 混沌 成為數學研究的嶄新領域,在實際應用中產生了巨大的影...
數值分析實驗報告
實驗專案 a 列主元高斯消去法 b 矩陣直接三角分解法 c 雅可比迭代法 d 高斯 賽德爾迭代法 e sor迭代法 f 冪法求最大特徵值 實驗名稱 列主元高斯消去法 實驗目的 熟悉了解高斯消去法的原理,並知道如何用高斯消去法求方程組的解 實驗所用的軟體 win tc 版本 1.9.0.0 例項用高斯...
數值分析實驗報告
課題一 解線性方程組的直接方法 1.實驗目的 1 通過該課題的實驗,體會模組化結構程式設計方法的優點 2 運用所學的計算方法,解決各類線性方程組的直接演算法 3 提高分析和解決問題的能力,做到學以致用 4 通過三對角形線性方程組的解法,體會稀疏線性方程組解法的特點。2.實驗過程 實驗 include...