资源简介
用于计算数值积分的龙贝格算法MATLAB程序,文件为.m文件,在MATLAB可直接运行。
代码片段和文件信息
function[] = Romberg(abmMd)
disp(‘请输入分子多项式a,分母多项式b,积分下限m,积分上限M,计算精度‘)
f = poly2sym(a)/poly2sym(b) %用于给用户显示被积函数的形式
disp(‘利用梯形公式计算数值积分的结果‘)
kk = zeros();
kk(11)=1/2*(M-m)/1*(subs(f‘x‘m)+subs(f‘x‘M)); %先存储首项?
for i=1:1:2^40
t=0;
for j=0:1:2^(i-1)-1
v=m+(2*j+1)*(M-m)/(2^i);
vv=polyval(av)/polyval(bv);
t=t+(M-m)/(2^i)*vv;
end
y=1/2*kk(i1)+t; %通项公式计算各项值?????????
kk(i+11)=y %存储其他项?
if(abs(1/3*(kk(i+11)-kk(i1)))<=d) %判断梯形公式值是否达到要求????
disp(‘The?result?is:‘)
kk()
kk(i+11)%梯形值满足要求,输出结果??????????
break;
else
s=(4*kk(i+11)-kk(i1))/(4-1);%构造simpson各项?????????
kk(i+12)=s %存储????????????
if(i+1>=3)
if(i+1>=3 & abs(1/15*(kk(i+12)-kk(i2)))<=d)
kk()
disp(‘The result is:‘)
- 上一篇:MATLAB 贪吃蛇
- 下一篇:用MATLAB实现彩色图像的灰度化处理
相关资源
- matlab_OFDM调制解调(来自剑桥大学)
- Matlab路面裂缝识别69319
- 高灵敏度GPS接收机MATLAB仿真,附捕获
- 基于MATLAB的质点弹道计算与外弹道优
- 阵列天线的matlab仿真
- MATLAB 经典程序源代码大全
- MATLAB小波软阈值去噪代码33473
- 天线阵的波束形成在MATLAB仿真程序及
- 非线性SVM算法-matlab实现
- 《MATLAB 智能算法超级学习手册》-程序
- 组合导航matlab程序
- 读取txt文件内容matlab代码实现
- Matlab实现基于相关的模板匹配程序
- matlab优化工具箱讲解
- 基于MATLAB的快速傅里叶变换
- 光纤传输中的分布傅立叶算法matlab实
- 基于matlab的图像处理源程序
- matlab 椭圆拟合程序
- 算术编码解码matlab源代码
- optical_flow 光流法 matlab 实现程序
- 引导图像滤波器 Matlab实现
- 分形几何中一些经典图形的Matlab画法
- OFDM系统MATLAB仿真代码
- SVM工具箱(matlab中运行)
- 图像小波变换MatLab源代码
- LU分解的MATLAB实现
- 冈萨雷斯数字图像处理matlab版(第三
- 替代数据法的matlab程序
- 用matlab实现的多站定位系统性能仿真
- 通过不同方法进行粗糙集属性约简m
评论
共有 条评论