Trapezoidal rule in a for loop

8 次查看（过去 30 天）
Mathias Pettersen
Mathias Pettersen 2019-9-23
评论： James Tursa 2019-9-24
d
  2 个评论
awezmm
awezmm 2019-9-23
Are you getting an error or the values are just not coming out right
James Tursa
James Tursa 2019-9-24
Mathias, it is considered rude to delete your question after it has received answers. This also can render this thread useless to other readers who might be having similar problems.

请先登录，再进行评论。

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

采纳的回答

James Tursa
James Tursa 2019-9-23
编辑：James Tursa 2019-9-24
You made a good start, but there are several problems with your summation. You don't accumulate into your trapez variable, you don't use your index i in your function evaluation, you have your endpoints evaluated multiple times inside your loop, etc.
It would probably be best for you to code the logic up exactly as it appears in the formula. So, first do the summation. Then add the endpoints. Then multiply by (b-a)/n. Also, why not make your life simpler and use the same index variable k that is used in the formula? So, the outline would be this:
trapez = 0; % Initialize your summation to 0
for k = 1: n - 1 % use k as the index so it looks like the formula
trapez = trapez + calculateFormula(_____); % you fill in the value here based on formula
end
trapez = _____ + trapez + _____; % you add in the end points here
trapez = _____ * trapez; % you fill in the multiplying factor here
Try to fill in the blanks above to see if you can get it to work. If you have more problems, come back and ask for more help with the specific problems you are having.
  2 个评论
Mathias Pettersen
Mathias Pettersen 2019-9-24
thanks so much. Does this look right?
a = -3;
b = 3;
n = 10;
trapez = 0;
for k = 1: n - 1
trapez = trapez + sum(calculateFormula(a + (k.*(b-a)./n)) + calculateFormula(b)./2);
end
trapez = trapez + calculateFormula(a)/2;
trapez = trapez * ((b-a)/n)
James Tursa
James Tursa 2019-9-24
编辑：James Tursa 2019-9-24
No, it doesn't. But you are getting closer. You don't need the sum( ) function inside your loop since you are in fact doing the summation manually. Also, the f(b)/2 is one of the endpoints that gets added after the loop (look at your formula closely and you will see that it is not intended to be inside the summation).

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Pulse and Transition Metrics 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by