matlab中几种求积分的方法

matlab中几种求积分的方法

最近看了几天matlab的基础知识，总结一下求简单积分的方法

举例图形

x = linspace(0, 2*pi, 100);
y = sin(x) + cos(x);
plot(x, y);

line([0, 2*pi], [0, 0]);
set(gca, 'XLim', [0, 2*pi]);
set(gca, 'XTick', 0:pi/2:2*pi);
set(gca, 'FontName', 'symbol');
set(gca, 'XTickLabel', {'0', 'p/2', 'p', '3p/2', '2p'});

Midpoint Rule

fi = f((x1 + x2) / 2)
取两点之间中点的平均值作为该区域的高

h = pi/100;
x = 0:h:2*pi;
mid = (x(1:end-1) + x(2:end)) / 2;
y = sin(x) + cos(x);
s1 = h * sum(y)

Trapezoid Rule

fi = (f(x1) + f(x2)) / 2;
取两点函数值得平均值作为该区域的高

h = pi/100;
x = 0:h:2*pi;
y = sin(x) + cos(x);
yy = (y(1:end-1) + y(2:end))/2;
s2 = h * sum(yy)

除此以外，matlab中的trapz函数可用来调用求其值，如下

h = pi/100;
x = 0:h:2*pi;
y = sin(x) + cos(x);
s3 = h * trapz(y)

1/3 Simpson’s

fi = 3/h*(f0 + 4f1 + f2),利用公式求解即可

h = pi/100;
x = 0:h:2*pi;
y = sin(x) + cos(x);
s4 = h * (y(1) + y(end) + 2*sum(y(3:2:end-2)) + 4*sum(y(2:2:end-1)))

利用integral

f = @(x) sin(x) + cos(x);
s5 = integral(f, 0, 2*pi)

结果显示

对比一下各种方式求出的积分值