Romberg积分法(MATLAB)

% fprintf('请输入区间最小值:\n')
% a=input('');
% fprintf('请输入区间最大值:\n')
% b=input('');
% fprintf('请输入N:\n')
% N = input('');
% fprintf('请输入精度:\n')
% ac = input('');
a=0;
b=1;
N=25;
ac=0.000001;
T = [];
S = [];
C = [];
R = [];
T1=0;
T2=0;
S1=0;
S2=0;
C1=0;
C2=0;
R1=0;
R2=0;
tol = 0;
t = 0;
m=1;
n=1;
h =(b-a)/n;
T1 = h*(fun_rom(b)-fun_rom(a))/2;
T(1)=T1;
n=2;
h =(b-a)/n;
% i = 2;
for i=2:N
% while tol<acii=2^(i-1);sum=0;for k=1:iisum = sum + fun_rom(a+(k-1/2)*h);endT2 = T1/2+h*sum/2;T(i) = T2;S2 = (4*T2-T1)/3;S(i)=S2;if m~=1C2 = (16*S2-S1)/15;C(i)=C2;endif m~=2R2 = (64*C2-C1)/63;R(i) = R2;endif m~=3tol = abs(R2-R1);if tol~=0if tol<act = t+1;break;endendendR1=R2;C1=C2;S1=S2;T1=T2;h=h/2;m=m+1;
%         if i>=N
%            break;
%         endi=i+1;
end

函数

function y=fun_rom(x)
% y = exp(x)*(x^2);
% y= exp(x)*sin(x);
% y = 4/(1+x^2);
% y = 1/(x+1);
% y=sin(x)/x;
y=cos(x)/((1-x)^2);
% y=cos(x)/(1-x)^2;
% y=(cos(x)/(1-x))^2;
end

Romberg积分法(MATLAB)相关推荐

Romberg积分法MATLAB实现（附代码、实例、详解）
第一部分:问题分析 (1)实验题目:龙贝格积分算法具体实验要求:用matlab编写龙贝格积分的代码,要求代码实现用户输入了被积函数.积分区间.精度之后,龙贝格积分表(T-数表). (2)实验目的:让 ...
Romberg积分法计算定积分（Matlab程序）
%Romberg积分法计算定积分 %参考教材:<数值分析>李乃成,梅立泉,科学出版社 %<计算方法教程>第二版凌永祥,陈明逵 clear;clc;close all; for ...
Romberg（龙贝格）积分法 | matlab
function s=rombg(a, b, TOL) % a 为积分下限 % b 为积分上限 % TOL 为误差容限 % s 为最后Romberg积分的值 % rombg_f为积分函数n = 1; ...
Matlab 使用Romberg积分法进行求解
a = -1; %下界 b = 4; %上界 epsilon = 1e-6; %精确度 tol = epsilon; fun = @(x)(1/(1+x^2)); %句柄函数 k=0; % 迭代次数 ...
四阶龙格积分法 matlab,matlab控制系统计算机仿真实验-完整版
MALTAB 仿真实验指导书实验一实验题目: 欧拉法&梯形法的MATLAB 实现实验目的: 1.熟练掌握MATLAB 的使用方法 2.牢记欧拉法.梯形法的计算过程 3.熟悉欧拉法.梯形法 ...
精细时程积分法matlab程序,MATLAB内燃机轴系振动计算程序
用于计算内燃机轴系振动的MTALAB程序,用基于高斯勒让德积分的精细时程积分法进行方程组求解,精度非常高.需要输入内燃机载荷. djwy910=max(xyz1(10,0.0*nn+1:nn-0)); ...
matlab 符号积分法,MATLAB积分方法
一相关函数 %符号积分 int(f,v) int(f,v,a,b) %数值积分 trapz(x,y)%梯形法沿列方向求函数Y关于自变量X的积分 cumtrapz(x,y)%梯形法沿列方向求函数Y关于 ...
Richardson外推加速技术（含Romberg详细分析）的Matlab实现
一.Romberg求积公式的演化根据前面一篇博客<基于复化梯度求积的求积步长自适应matlab实现 >,我们知道可以通过对积分区间不断进行二分,然后采用复化梯形公式计算得到新的求积结果. ...
Romberg（龙贝格）积分法
变长梯形法方法简单,但精度低而且收敛速度慢.如果在保留这种方法使用的逐步二分思想的基础上,将用误差事后估计法所求得的误差作为积分近似值的补偿值,以便进一步提高精度,同时又引入加速收敛的技术,则就可以得 ...

Romberg积分法(MATLAB)

Romberg积分法(MATLAB)相关推荐

最新文章

热门文章