(一)专题实验(Newton-Cotes积分公式)
1、编写[a,b]上梯形积分公式、Simpson积分公式。
2、利用自己编写的程序计算定积分,计算一下数值解和精确解之间差的绝对值。
梯形积分:
function T=TX_int(f,a,b)
T=(b-a)/2*(f(a)+f(b));
>> TX_int(@(x)cos(x),0,pi/4)
ans =文章来源:https://www.toymoban.com/news/detail-445919.html
0.6704
function T=TX_int(f,a,b)
T=(b-a)/2*(f(a)+f(b));
a=abs(sqrt(2)/2-T);
disp(a)
>> TX_int(@(x)cos(x),0,pi/4)
0.0367
ans =
0.6704
Simpson积分;
function T=TX_int(f,a,b)
T=(b-a)/6*(f(a)+4*f((a+b)/2)+f(b));
>> TX_int(@(x)cos(x),0,pi/4)
ans =
0.7072
function T=TX_int(f,a,b)
T=(b-a)/6*(f(a)+4*f((a+b)/2)+f(b));
d=abs(sqrt(2)/2-T);
disp(d);
>> TX_int(@(x)cos(x),0,pi/4)
9.5166e-05
ans =
0.7072文章来源地址https://www.toymoban.com/news/detail-445919.html
到了这里,关于【matlab】数值积分公式的程序实现的文章就介绍完了。如果您还想了解更多内容,请在右上角搜索TOY模板网以前的文章或继续浏览下面的相关文章,希望大家以后多多支持TOY模板网!