Matlab 使用Romberg积分法进行求解

a = -1; %下界
b = 4; %上界
epsilon = 1e-6; %精确度
tol = epsilon;
fun = @(x)(1/(1+x^2)); %句柄函数
k=0; % 迭代次数
n=1; % 区间划分个数
h=b-a;
T=h/2*(fun(a)+fun(b));
err=1;
while err>=tolk=k+1;h=h/2;tmp=0;for i=1:ntmp=tmp+fun(a+(2*i-1)*h);endT(k+1,1)=T(k)/2+h*tmp;for j=1:kT(k+1,j+1)=T(k+1,j)+(T(k+1,j)-T(k,j))/(4^j-1);endn=n*2;err=abs(T(k+1,k+1)-T(k,k));
end
R=T(k+1,k+1)

看完记得点赞，（づ￣3￣）づ╭❤～

Matlab 使用Romberg积分法进行求解相关推荐

Romberg积分法计算定积分（Matlab程序）
%Romberg积分法计算定积分 %参考教材:<数值分析>李乃成,梅立泉,科学出版社 %<计算方法教程>第二版凌永祥,陈明逵 clear;clc;close all; for ...
数值分析复化求积matlab,数值分析：复化公式求解定积分
复化公式求解定积分用熟悉的计算机语言编程上机完成 (1)用复化梯形公式.复化Simpson公式和复化Cotes公式计算积分,自己设置不同精度要求,对结果进行比较分析. (2)用Romberg积分法计 ...
fick第二定律matlab求解,(最新整理)实验二MATLAB数值计算常微分方程(组)的求解
<(最新整理)实验二MATLAB数值计算常微分方程(组)的求解>由会员分享,可在线阅读,更多相关<(最新整理)实验二MATLAB数值计算常微分方程(组)的求解(19页珍藏版)> ...
matlab绘制圆，且求解两个圆的交点坐标
前言这里简单记录下对于matlab绘制圆,以及求解两圆交点的坐标的实现(包括C++). 一绘制圆 circle.m文件: function [] = circle(x, y, r, color) ...
Algorithm之PrA：PrA之nLP非线性规划算法+Matlab 优化工具箱的GUI求解非线性规划
Algorithm之PrA:PrA之nLP非线性规划算法+Matlab 优化工具箱的GUI求解非线性规划目录 PrA之nLP非线性规划算法操作图文教程 PrA之nLP非线性规划算法 (1).编写M ...
pagerank算法实现matlab,Matlab 入门及PageRank算法求解.ppt
Matlab 入门及PageRank算法求解矩阵运算 + 矩阵加 - 矩阵减 * 矩阵乘 / 矩阵左除 \ 矩阵右除 ^ 矩阵幂维数相同才能加减:方阵才能求幂. 注意左乘,右乘. a\b?求a* ...
matlab行星运动轨迹仿真动画,利用Matlab可视化功能实现微分方程求解行星运动轨迹...
利用Matlab可视化功能实现微分方程求解行星运动轨迹 1.背景在物理学璀璨的发展史上,物理学家花了很长时间研究我们头顶浩瀚的星空,试图探究星星的运行模式,以及地球自身的运动模式.其中不乏像亚里士多 ...
matlab使用杂谈4-偏微分方程求解之pdede函数使用
matlab使用杂谈4-偏微分方程求解之pdede函数使用偏微分方程求解偏微分方程的数值方法 Matlab解偏微分方程 pdepe()函数 pdepe函数使用示例 PDE方程求解格式 PDE方程初 ...
【优化算法】基于matlab量子粒子群算法求解单目标优化问题【含Matlab源码 2203期】
⛄一.获取代码方式获取代码方式1: 完整代码已上传我的资源:[优化算法]基于matlab量子粒子群算法求解单目标优化问题[含Matlab源码 2203期] 点击上面蓝色字体,直接付费下载,即可. 获 ...

Matlab 使用Romberg积分法进行求解

Matlab 使用Romberg积分法进行求解相关推荐

最新文章

热门文章