Creating matrix nested for loop

I need to write a matrix where the first column has the varying mass of an ODE, the second to n-th column all the y(t) values of that ODE. I tried to solve this with a nested for loop but somehow fail so hard eventhough I am sure it is pretty easy.
for row_index = 1:10
for col_index = 1:length(time)
Matrix(row_index,column_index) = [mass(row_index)];
end
end
Somehow something after "[mass(row_index)]" is missing or do you know a way to make it more elegant?
Thanks in advance!

1 个评论

@Detox: Please do not cross-post a question. If it is really important to post a question in different forums, please add at least links to the other forum. Otherwise the voluntary helpers waste time with writing answers, which have been given elsewhere already. Thanks.

请先登录,再进行评论。

回答(1 个)

Chaya N
Chaya N 2016-10-22
编辑:Chaya N 2016-10-22
If you already have all these data available (in separate vectors, I assume) then you do not need any loops.
For example, if you had
x = [all the ODE mass values];
y1 = [values for y(1)]; % at time = 1
y2 = [values for y(2)]; % at time = 2
.
.
.
yt = [values for y(t)]; % at some time 't'
you could do:
matrix = [x y1 y2 ... yt]; if all of these are column vectors and of the same length
or
matrix = [x; y1; y2; ...; yt]'; if all of these are row vectors of the same length
NOTE: Please observe the differences in syntax between the two versions!

4 个评论

Detox
Detox 2016-10-22
编辑:Detox 2016-10-22
Thanks for your answer. But I want to solve this Problem with a nested loop if possible. Any ideas on How i could Achieve that?
Alright! Here is a simple illustration. Assume that your ODE mass values are represented by 'x', your y(t) functions are represented by the anonymous function 'y', time values are contained in 't' and the final matrix is called 'm'.
x = 1:10;
y = @(x,t) x*t^2; % just an example
t = 1:10;
m = [x' zeros(numel(x),numel(t))]; % preallocate and fill in first column with 'x'
for r = 1:numel(x) % row counter
for c = 1:numel(t) % column counter
m(r,c+1) = y(x(r),t(c)); % note the 'c+1' column subscript for the output matrix
end
end
If you still have trouble, please also include your equation for y(t) when you post here.
Function is:
function ydot = Exp_01(t,y,m)
ydot = zeros(size(y));
ydot(1) = y(2);
ydot(2) = (-y(2)-100*y(1))./m;

请先登录,再进行评论。

类别

帮助中心File Exchange 中查找有关 Programming 的更多信息

标签

提问:

2016-10-22

评论:

Jan
2016-10-24

Community Treasure Hunt

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

Start Hunting!

Translated by