Saving vectors from a loop into a matrix

Question

KJones 2016-11-12

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/311840-saving-vectors-from-a-loop-into-a-matrix

评论： Changoleon 2016-11-13

I have this code that I made and I would like to save my final vector as a row in a matrix as I would like to plot the averages. I would like this matrix to consist of each row representing a run at a different H value where the H values change from R>0.1 up until R>0.9. How can I do this automatically because I don't want to change my matrix A so it would need to be done in one run?

A=createsrandomnetwork(100,5); %Network Matrix
I=zeros(100,1);
positions=[1,1];
I(positions)=1;
v=zeros(1,100);        %Vector which proportion value will be saved into
n=numel(I);            %n=number of array elements
v(1)=nnz(I)/n;         %At time 0 the proportion infected (nnz=number of non zero elements in I)
for i=1:100
a = rand(100);
R = triu(a) + triu(a,1)';
H=R>0.5;     %Matrix with % chance of passing on infection
P=A.*H;      %New Network Matrix including the probability
P(logical(eye(size(P)))) = 0; %makes sure the diagonal is 0 then people can recover
I=P*I;      %New Network Matrix times the Vector of Infected
v(i)=nnz(I)/n        %Proportion infected at time i
end

Thank you

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Changoleon 2016-11-13

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/311840-saving-vectors-from-a-loop-into-a-matrix#answer_243078

编辑：Walter Roberson 2016-11-13

在 MATLAB Online 中打开

Hi. This is a very simplified solution, I believe you should be able to modify this code.

h=5;
for i=1:h
x(i) = i * 10;
s(i,:) = [i,x(i)];
end
Sina

2 个评论
显示无隐藏无

KJones 2016-11-13

在 MATLAB Online 中打开

Thank you, I have managed to get it so it saves in the different columns each time but it creates a different matrix.

A=createsrandomnetwork(20,3);
for k= 1:9 
v=zeros(1,20);        %Vector which proportion value will be saved into
n=numel(I);
I=zeros(20,1);
positions=[1,1];
I(positions)=1;
a = rand(20);
R = triu(a) + triu(a,1)'; 
x = k;
H=R>x/10;
for i=1:20  
P=A.*H;
P(logical(eye(size(P)))) = 0; %makes sure the diagonal is 0 then people can recover
I=P*I;
v(i)=nnz(I)/n  
end 
s=zeros(20,9);
s(:,k)=v
end