how to plot trapezoidal rule

Question

alp onur karacay 2021-12-29

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1619100-how-to-plot-trapezoidal-rule

回答： Pratyush Roy 2022-1-7

Hello i coded trapezoidal rule but how can i plot it ?here is my code

t = 0:0.2:2*pi;
f = @(x) 3*cos(t*1000);
a = 0;
b = 2*pi;
n  = length(t)-1; 
h =(b-a)/n;
for i=1:1:n-1
    sum = sum + f(a+i*h);
end
result = h/2*(f(a)+f(b)+2*sum);
fprintf('%f',result);

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Torsten 2021-12-29

Your variable "result" is just a single number, namely the approximated area under the function f in the interval [0:2*pi].

So what do you want to plot ? How the area develops from 0 to 2*pi ?

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Pratyush Roy 2022-1-7

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1619100-how-to-plot-trapezoidal-rule#answer_869910

在 MATLAB Online 中打开

Hi,

One can use the polyshape function to create polygons (trapezium in this case) and use hold on and off to show them in a single figure. The following code might be helpful to understand how this works:

h = 0.2;
t = 0:h:2*pi;
f = @(x) 3*cos(x*1000);
a = 0;
b = 2*pi;
n  = length(t); 
sum1 = 0;
plot(t,f(t));
hold on;
for i=0:1:n-1
    sum1 = sum1 + f(a+i*h);
    if (i>=0)
        plot(polyshape([a+i*h,a+i*h,a+(i+1)*h,a+(i+1)*h],[0,f(a+i*h),f(a+(i+1)*h),0]),'FaceColor','red');
        hold on;
    end    
end
hold off;
result = h/2*(f(a)+f(b)+2*sum1);