資源簡介
% Integral function f from a0 to b0
% tol0: Total Error limit
% flag:0-->Trapezoid ; 1-->Simpson
代碼片段和文件信息
function?[intPointNumber]?=?AdaptiveIntegral(?fa0b0tol0flag?)
%?Integral?function?f?from?a0?to?b0
%?tol0:?Total?Error?limit
%?flag:0-->Trapezoid?;?1-->Simpson
int=0;n=1;a(1)=a0;b(1)=b0;tol(1)=tol0;app(1)=trapezoid_or_Simpson(fa(1)b(1)flag);PointNumber=0;
while?n>0???????????????????????????????????%n是當前列表結束的位置
????c=(a(n)+b(n))/2;oldapp=app(n);
????app(n)=trapezoid_or_Simpson(fa(n)cflag);app(n+1)=trapezoid_or_Simpson(fcb(n)flag);
????if?abs(oldapp-(app(n)+app(n+1)))<3*tol(n)
????????int=int+app(n)+ap
- 上一篇:Matlab單像素成像算法比較
- 下一篇:ID3決策樹-MATLAB程序
評論
共有 條評論