I keep getting 'Index exceeds array bounds.' my code doesn't like r(i)

2 次查看（过去 30 天）
Jessica Gerberdolan
Jessica Gerberdolan 2018-11-13
评论： Sudhaunsh Deshpande 2021-4-28
%function [int_ft] = gerberdolan_jessica_hw09_int(t, f)
%INPUT:
%t = vector of t values at which the function is provided
%f = vector of values of the function at the points contained in t
%OUTPUT:
%int_ft = value of integral
%start this function as a script with the following test senario
t = [0 pi/4 pi/2 3*pi/4 pi]
f = cos(t)+5
for i = 1:(length(t)-1)
t = t(i);
r(i) = f(i).*(t(i+1)-t(i))
%((t(1+1)-t(i)) is delta x
%f(t(i)) is the base rantagle for the retangular rule
%r is rectangle
tr(i) = (f(t(i))-f(t(i+1)))*(t(1+1)-t(i))*1/2
%tr is the triangle on top
%(f(t(i))-f(t(i+1))) is the height difference
%(t(1+1)-t(i)) is delta x
%1/2 is to make it a triangle
a_of_trap(i) = r(i) + tr(i)
%area of trapizoid is the rectangle plus the triangle
end
int_ft = sum(a_of_trap)

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

采纳的回答

Adam Danz
Adam Danz 2018-11-13
编辑：Adam Danz 2018-11-13
Initially in your code, t is a vector.
t = [0 pi/4 pi/2 3*pi/4 pi]
t =
0 0.7854 1.5708 2.3562 3.1416
Then within your i loop, you overwrite t as a scalar.
t = t(i)
t =
0
Immediately after this, you try to grab the 2nd element of t even though t now only has 1 element. This is why the index (i+1) exceed array bounds. Because at this point, t is only 1 element long.
r(i) = f(i).*(t(i+1)-t(i))

更多回答（1 个）

Guillaume
Guillaume 2018-11-13
"my code doesn't like r(i)"
There's no liking involved. Your code does exactly what you tell it to. If you tell it to do something impossible it's not going to work.
I haven't looked further than
t = t(i);
to know that your code is never going to work. That line replace the whole vector t by a scalar, such that indexing t after that is always going to be an error.
Note that you can easily debug your own code which would be the best way for you to understand where you're going wrong. Step through your code one line at a time in the debugger. See what happens to your variables after each statement and check that it conforms to your expectations.
  1 个评论
Jessica Gerberdolan
Jessica Gerberdolan 2018-11-13
I appreciate you help! I have been attempting to de-bug my code for over thirty minutes now and I figured I needed an extra pair of eyes

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Programming 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by