數值積分公式及龍貝格(Romberg)演算法實現matlab
一、數值積分方法得基本思想
二、復化求積公式
三、基於復化梯形求積公式的高精度求積演算法——Richardson外推法和Romberg演算法
四、Romberg演算法的matlab程式:
function Romberg
% 龍貝格(Romberg數值求解公式)
% inputs:
% -fun:積分函式控制程式碼
% -a/b:積分上下限
% -tol:積分誤差
% Outputs:
% -R:4階精度Romberg積分值
% -k:迭代次數
% -T:整個迭代過程
%
% Example
% fun=@(x)x^5;
a = 0;
b = 1;
epsilon = 1e-6;
tol = epsilon;
fun = @(x) x^5;
k=0; % 迭代次數
n=1; % 區間劃分個數
h=b-a;
T=h/2*(fun(a)+fun(b));
err=1;
while err>=tol
k=k+1;
h=h/2;
tmp=0;
for i=1:n
tmp=tmp+fun(a+(2*i-1)*h);
end
T(k+1,1)=T(k)/2+h*tmp;
for j=1:k
T(k+1,j+1)=T(k+1,j)+(T(k+1,j)-T(k,j))/(4^j-1);
end
n=n*2;
err=abs(T(k+1,k+1)-T(k,k));
end
R=T(k+1,k+1)
end
可以直接執行的,只出來積分的結果。
相關文章
- 玩轉matlab之一維 gauss 數值積分公式及matlab原始碼Matlab公式原始碼
- 玩轉 matlab 之二維 gauss 數值積分公式使用及 matlab 原始碼(1)-常量區間Matlab公式原始碼
- MATLAB學習筆記:數值積分Matlab筆記
- 數值計算:高斯-勒朗德積分公式公式
- 實驗三 數值積分(android)Android
- 差分進化演算法介紹及matlab實現演算法Matlab
- 貝塞爾曲線原理、推導及Matlab實現Matlab
- 【數值計算方法】數值積分&微分
- 蟻群演算法原理及Matlab實現演算法Matlab
- 樸素貝葉斯分類和預測演算法的原理及實現演算法
- matlab對不定積分和定積分的計算Matlab
- 基於有限體積法和交錯網格的SIMPLE演算法推導及實現演算法
- 奇異值分解以及matlab實現Matlab
- java實現有理數分數相關計算(淺度拷貝)案例Java
- 高通驍龍675跑分多少?驍龍675效能跑分與規格引數詳解
- MATLAB 檢驗資料正態分佈及程式碼實現Matlab
- # 數值計算:三角形積分
- 數值分析Python實現系列—— 一、拉格朗日插值法Python
- 全概率公式、貝葉斯公式公式
- 基於粒子群演算法的分組揹包MATLAB實現演算法Matlab
- 使用Matlab實現AHP演算法Matlab演算法
- matlab如何給未知數及包含未知數的函式賦值Matlab函式賦值
- 貝葉斯實現文字分類C++實現文字分類C++
- 牛頓插值 C++ 和 Matlab實現C++Matlab
- 貝葉斯公式公式
- 【數值框架】戰鬥公式:傷害公式選擇框架公式
- 求取眾數及重數的演算法實現(Java)演算法Java
- FFT演算法實現與分析MATLABFFT演算法Matlab
- NSGA2 演算法Matlab實現演算法Matlab
- 二維或三維的分佈積分方法(格林公式)7公式
- 樸素貝葉斯實現文件分類
- 一分鐘瞭解“Matlab統計數值頻率和個數tabulate”Matlab
- 定積分之換元積分法公式(或稱變數替換法)公式變數
- 排隊論演算法的matlab實現演算法Matlab
- kmeans聚類演算法matlab實現聚類演算法Matlab
- 蒙特卡羅演算法的matlab實現演算法Matlab
- 龍叔python-直接賦值,深拷貝,淺拷貝的簡單解析Python賦值
- Word帶數學公式釋出部落格公式