0 个投票

Hi, All!
I have 3 "for" loops but i confuse when combining one and another. I need vector u, depend on x(j)-x(i) and x(i,t)-x(i,t-1)
this is the best function i've tried before. So glad to know that anyone can correct and help me :)
for t=1:1000
for i=1:5
for j=1:5
u(i,j)= a1*adj2(i,j)*(x(j)- x(i))+a2*(x(i,t)-x(i,t-1));
end
end
u=sum(u,2)
end

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

tmarske
tmarske 2020-3-16
Please use auto-indentation when nesting multiple loops like this - it makes things much easier to read.
What are the sizes of a1, adj2, a2 and x?
You are indexing the array 'x' in two different ways on the same line:
(x(j)- x(i))
and
(x(i,t)-x(i,t-1))
Both of these will return output if x is a 2-d array, but I suspect it's not what you meant to do.
Bob Thompson
Bob Thompson 2020-3-16
Also, what exactly is the problem you're experiencing? Are you getting an error message, or are you just not getting the output you would expect. If the output is different, what does it look like, and how were you expecting it to be different?
To add to tmarske's point about the different forms of indexing, if x is a 2D array then x(j) - x(i) will return a reduced vector, while x(i,t) - x(i,t-1) will return a single element. This is not inherently a problem, depending on the size of a1, adj2, and a2, but MATLAB will not let you perform a '+' operation on two components of different sizes like that.
Nur Amalina
Nur Amalina 2020-3-17
First, thanks for the reply. a1 and a2 are only a constant value while u and x are 2D array. row by n. column by t. in this case it will be 5x1000. adj2 is 5x5 matrix.
I am getting an error message "Index in position 2 exceeds array bounds (must not exceed 1)." hmm maybe the problem like Bob said. "x(j) - x(i) will return a reduced vector, while x(i,t) - x(i,t-1) will return a single element" any advice what should i do to make it be the same size?

请先登录，再进行评论。

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

 采纳的回答

Sriram Tadavarty
Sriram Tadavarty 2020-3-17
编辑：Sriram Tadavarty 2020-3-17

0 个投票

Hi Nur,
If the code is updated to run from time t set from 2 to 1000, this should go away with the error you see.
The subtraction of x(j) and x(i) would only do it for the first five elements all the time. I think this is not the intention. It should be x(j,t) and x(i,t)
With these changes the code looks as such:
for t=2:1000
for i=1:5
for j=1:5
u(i,j,t)= a1*adj2(i,j)*(x(j,t)- x(i,t))+a2*(x(i,t)-x(i,t-1));
end
end
u=sum(u,2)
end
u = permute(u,[1 3 2]);
Hope this helps.
Regards,
Sriram

6 个评论
显示 4更早的评论 隐藏 4更早的评论

Nur Amalina
Nur Amalina 2020-3-17
Hi Sriram, ive tried the code but the error message still show up :(
Sriram Tadavarty
Sriram Tadavarty 2020-3-17
Hi Nur, I do see the outputs for the code. Did you change the variable t?
Nur Amalina
Nur Amalina 2020-3-17
yes, i did. what is the output did you get? is it u(1:5,1:1000)?
Sriram Tadavarty
Sriram Tadavarty 2020-3-17
Should the output be of size 5 x 1000, then you need to even use the third variable t.
Like,
u(i,j,t)= a1*adj2(i,j)*(x(j,t)- x(i,t))+a2*(x(i,t)-x(i,t-1));
% At the end of all the loops
u = permute(u,[1 3 2]);
Let me know what is intended in u?
Nur Amalina
Nur Amalina 2020-3-17
problem solve Sriram. Thank you very much!
Sriram Tadavarty
Sriram Tadavarty 2020-3-17
That's nice to hear Nur. Do accept the answer if it helped.

请先登录，再进行评论。

更多回答（0 个）

类别

在 帮助中心 和 File Exchange 中查找有关 Loops and Conditional Statements 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by