matlab实现复化Newton-Cotes公式求积分的程序应用和代码.pdf
版权申诉
126 浏览量
2021-09-14
15:50:44
上传
评论
收藏 113KB PDF 举报
执行函数为 mymulNewtonCotes.m
1、使用方法：
Step1：在 MATLAB 命令窗口输入被积函数
t e dt
。
0
1
t
2
2
输入应为：ft=@(t)t.*exp(t^2/2)。
Step2：执行函数。输入形式为 mymulNewtonCotes(ft,a,b,m,n)；
其中ft—被积函数，此体重ft=@(t)t.*exp(t^2/2)，已经在第一步赋值；
a—积分下限，本题中为0；
b—积分上限，本题中为1；
m—将区间[a,b]等分的子区间数量，本题可选为10；
n—采用的Newton-Cotes公式的阶数，必须满足n<8，否则积分没法保
证稳定性。
当n=1时，即为复化梯形公式；n=2时，即为复化复化辛普森公式。
所以，分别输入mymulNewtonCotes(ft,0,1,10,1)和
mymulNewtonCotes(ft,0,1,10,2)就可以得到两种方法的积分计算结果。
2、计算结果
而根据积分运算，可得：
1
t
2
x
2
x
t e dt e d ( ) e dx e
0
0
0
2
1 1
t
2
2
t
2
2
1
2
0
e e
0
1.6487 1 0.6487
1
2
说明复化梯形和复化辛普森公式计算出的结果基本一致，与实际结果相符。
资源评论
