matlab中龙贝格积分,龙贝格积分-matlab通用程序

数值分析课程的大作业，教材《数值分析》李乃成.梅立泉

clear

clc

format long

f=input(‘请输入原函数f=’,’s’);

a=input(‘积分下限a=’);

b=input(‘积分上限b=’);

eps1=input(‘精度eps1=’);

T(1)=double((b-a)/2*(limit(sym(f),findsym(sym(f)),a)+limit(sym(f),findsym(sym(f)),b)));

for k=2:4

sum1=0;

for i=1:2^(k-2)

sum1=sum1+subs(sym(f),findsym(sym(f)),(a+(2*i-1)*(b-a)/2^(k-1)));

end

T(k)=1/2*T(k-1)+(b-a)/(2^(k-1))*sum1;

end

for k=1:3

S(k)=T(k+1)+1/(4-1)*(T(k+1)-T(k));

end

for k=1:2

C(k)=S(k+1)+1/(4^2-1)*(S(k+1)-S(k));

end

R(1)=C(2)+1/(4^3-1)*(C(2)-C(1));

k=3;

while 1

T(1)=T(2);

T(2)=T(3);

T(3)=T(4);

sum2=0;

for i=1:2^k

sum2=sum2+subs(sym(f),findsym(sym(f)),(a+(2*i-1)*(b-a)/2^(k+1)));

end

T(4)=1/2*T(4)+(b-a)/2^(k+1)*sum2;

S(1)=S(2);

S(2)=S(3);

S(3)=T(4)+1/(4-1)*(T(4)-T(3));

C(1)=C(2);

C(2)=S(3)+1/(4^2-1)*(S(3)-S(2));

R(2)=C(2)+1/(4^3-1)*(C(2)-C(1));

if abs(R(2)-R(1))

break;

end

R(1)=R(2);

end

vpa(R(2),9)

matlab中龙贝格积分,龙贝格积分-matlab通用程序相关推荐

matlab显示循环里进度,关于matlab中的waitbar(用于显示程序运行的进度)
关于matlab中的waitbar(用于显示程序运行的进度) 应用1: h = waitbar(x,'message'); %其中x必须为0到1之间的数,message为显示的信息,其实这个x大于1也 ...
matlab中特殊字符6,matlab中特殊字符
matlab特殊字符_工学_高等教育_教育专区.本文说明了matlab中如何输入特殊字符,如希腊字母字符映射表 C:\\WINDOWS\\system32\\charmap.exe M 文件里是可以输 ...
matlab中插入希腊,matlab中希腊字母打法
如何在matlab中输入希腊字母_数学_自然科学_专业资料.matlab数据可视化处理,绘图技巧如何在 matlab 中输入希腊字母 matlab 中用转义符来输入希腊字母的方法上标用 ..... ...
matlab excel导入矩阵,如何将excel中数据导入到matlab中!Excel导入matlab生成矩阵
如何将excel中数据导入matlab并存入数组中使用Matlab中的Excel Link就可以了. 如果使用Excel 2003版本,操作如下: 1 打开Excel. 2 在菜单栏中选择→加载项→ ...
MATLAB中hough反变换,matlab中hough变换
hough 变换提取直线(Matlab 实现) Hough 变换提取直线一.实验目标实现用 Hough 变换检测直线的算法二.实验内容 1.读入图像选取有较多直线及部分曲线以作对比的...... j= ...
matlab中的routte,matlab外文翻译外文文献英文文献MATALAB混合仿真平台控制算法的概述...
<matlab外文翻译外文文献英文文献MATALAB混合仿真平台控制算法的概述>由会员分享,可在线阅读,更多相关<matlab外文翻译外文文献英文文献MATALAB混合仿真平台控制算 ...
matlab中real函数,Matlab中del2()函数学习笔记
在L=del2(U)表达式中,Matlab帮助文件的解释是:li,j=(ui+1,j+ui-1,j+ui,j+1+ui,j-1)/4 - ui,j 但帮助文件没有解释边缘点如何处理.根据数学中的原理, ...
matlab中的图像,MATLAB中图像的基本操作
MATLAB中图像的基本操作 1.读取.显示图片 MATLAB中提供了immread()与imshow()函数读取和显示图片.其中读取函数imread()原型: imread: A= imread(f ...
quad8是matlab中调用那个,Matlab 数值积分
MATLAB数值积分与微分 8.1 数值积分 8.1.1 数值积分基本原理求解定积分的数值方法多种多样,如简单的梯形法.辛普生(Simpson)法.牛顿-柯特斯(Newton-Cotes)法等都是经 ...
matlab中的控制语句,MATLAB控制语句
目的:研究控制结构(用于, 是否, 切换, 中断, 继续, 输入/输出功能, 读取和存储数据). If:If评估逻辑表达式并根据表达式的值执行一组语句. If语句的语法 if expression 1 ...

matlab中龙贝格积分,龙贝格积分-matlab通用程序

matlab中龙贝格积分,龙贝格积分-matlab通用程序相关推荐

最新文章

热门文章