龍貝格 實驗報告

2022-09-23 16:27:09 字數 975 閱讀 8440

學院:數學科學學院專業:

實驗三龍貝格求積分

一、 實驗目的

用演算法實現龍貝積分運算,求解i=∫01x3/2dx的積分二、 實驗原理

三、 演算法設計與程式流程圖

演算法設計:(先假定所求積分二分最大次數次數為20)1、 先求t[k][0]

2、 再由公式求t[i][j]

3、 在求出的同時比較t[k][k]與t[k-1][k-1]的大小,如果二者之差的絕對值小於1e-5,就停止求t[k][k];此時的k就是所求的二分次數,而此時的t[k][k]就是最終的結果

4、 列印出所有的t[i][j];

程式流程圖:ny

n nn

y四、 源程式

#include<>

#include<>

#include<>

#define a 0 //區間下限#define b 1 //區間上限#define n 20 //最大加速次數double f(double x)

void main()

t[k][0] = 0.5*t[k-1][0] + h*s;

}for(k=1;k

if(fabs(t[k][k]-t[k-1][k-1])<1e-5)}

printf("r=%.6lf\n",t[k][k]);

printf("k=%d \n",k);

printf("t:\n");

for(i = 0;i<=k;i++)

}五、 程式執行

六、 結果分析

如上所示的結果與課本中求得的結果完全一樣,表明程式編寫正確,且符合要求,事實上,只要再將所求值的精度設定得更小,則所求的結果將更加準確,最終將無限接近於標準值,由上表也可以看出用龍貝格積分法求函式的積分值在精度比較低的情況下就能求到很準確的值!

七、 實驗體會

通過這次實驗上機,使我更進一步了解了龍貝格法的計算思想,其在精度上很***,收斂較快,是解積分問題的有效方法。

數值分析實驗報告 龍貝格演算法

利用matlab編寫龍貝格演算法,並測試。1龍貝格求積演算法matlab主程式 2呼叫函式解題 1寫出龍貝格求積演算法matlab的程式 function t rbg f,a,b,c t zeros 15,4 t 1,1 b a 2 f a f b for k 2 4 sum 0 for i 1 2...

數值計算方法龍貝格積分報告

實驗目的 1通過本實驗加深對 龍貝格公式 的認識 2通過c 編寫程式求解下列積分 e 1e 6 3預期達到題中結果 4注意實驗給出的程式有兩個子程式為 f f1使用時只需該程式名即可。實驗題目 復合求積公式計算定積分 實驗原理及理論基礎 用梯形公式得到的積分近似值的誤差是,因此,人們希望用這個誤差作...

微格實驗報告

課程名稱現代教育技術 實驗專案微格教學 實驗地點微格實驗室 指導教師張志潔 專業化學教育班級 08化本 1 班 學號 2008364155姓名代甜甜 2011 年 4 月 28 日 學校 湛江師範學院專業 化學教育課題 苯的結構授課時間 8分鐘受訓班級 08化本 1 班受訓人 代甜甜受訓技能 匯入和...