數值積分公式及龍貝格(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
- 樸素貝葉斯分類和預測演算法的原理及實現演算法
- 高通驍龍675跑分多少?驍龍675效能跑分與規格引數詳解
- java實現有理數分數相關計算(淺度拷貝)案例Java
- MATLAB 檢驗資料正態分佈及程式碼實現Matlab
- # 數值計算:三角形積分
- 基於粒子群演算法的分組揹包MATLAB實現演算法Matlab
- 使用Matlab實現AHP演算法Matlab演算法
- 數值分析Python實現系列—— 一、拉格朗日插值法Python
- matlab如何給未知數及包含未知數的函式賦值Matlab函式賦值
- 全概率公式、貝葉斯公式公式
- 牛頓插值 C++ 和 Matlab實現C++Matlab
- 貝葉斯實現文字分類C++實現文字分類C++
- 【數值框架】戰鬥公式:傷害公式選擇框架公式
- FFT演算法實現與分析MATLABFFT演算法Matlab
- NSGA2 演算法Matlab實現演算法Matlab
- 求取眾數及重數的演算法實現(Java)演算法Java
- 樸素貝葉斯實現文件分類
- 一分鐘瞭解“Matlab統計數值頻率和個數tabulate”Matlab
- 排隊論演算法的matlab實現演算法Matlab
- kmeans聚類演算法matlab實現聚類演算法Matlab
- 蒙特卡羅演算法的matlab實現演算法Matlab
- Word帶數學公式釋出部落格公式
- 龍叔python-直接賦值,深拷貝,淺拷貝的簡單解析Python賦值
- [MatLab]學習筆記2:MatLab數值資料Matlab筆記
- 神經網路 | 基於MATLAB 深度學習工具實現簡單的數字分類問題(卷積神經網路)神經網路Matlab深度學習卷積
- 淺談全概率公式和貝葉斯公式公式
- MATLAB實現頻數表——hist的使用Matlab
- 灰度影像直方圖均衡化公式及實現直方圖公式
- 機器學習之貝葉斯分類(python實現)機器學習Python
- 聚類演算法——DBSCAN演算法原理及公式聚類演算法公式
- 《機器學習實戰》基於樸素貝葉斯分類演算法構建文字分類器的Python實現機器學習演算法文字分類Python