Creating a Matrix with for loop

I have a fairly large vector named blender (45,000+). I have extracted n elements for which blender is greater than b (irrelevant). Now my difficulty is the following:

I am trying to create a 21 x n matrix with each extracted element of blender plus 10 things before, and the 10 things after, hence the reason why each column is 21.

I have been trying variations of the following:

   element=find(blender >= 120);

   NewMatrix = zeros(21, length(element));
   for i = 1:length(element)
       k = element(i);
       NewMatrix(:,i)= Matrix(blender(k-10:k+10));
   end

then I want to plot one column of the matrix at the time when I hit Enter. This second part I can figure out later, but I would appreciate some help making the Matrix

Thanks

Answer 1

[Cedric]

If I understand well and if you want to keep it simple with an explicit for loop, what you want to do is something like that I guess:

 element=find(blender >= 120);

 NewMatrix = zeros(21, length(element));
 for ii = 1:length(element)
    k = element(ii);
    r = max(1, k-10) : min(length(blender), k+10) ; % Range of indices +/- 10
                                                    % limited to valid ids.
    NewMatrix(1:length(r),ii) = blender(r);         % Store blender on the
                                                    % valid range; next 
                                                    % elements left to 0.
 end

Cheers,

Cedric

[Diego Yanez]

Thank you so much, I don't understand your solution much but it works exactly like I want it to. It didn't have to be a for loop, but that is the first idea that I got.

[Cedric]

You're welcome! About how it works, I just followed your initial idea indeed. The logic is as follows:

• You want to extract blocks from blender at each k +/-10, k going iteratively through all indices returned by find().
• During iteration ii, you want a priori to perform something like: NewMatrix(:,ii) = blender(r) , r being the set/range of indices that cover the block k +/-10. This is almost what you implemented.
• Now there is a small additional difficulty; if k<=10, k-10 is not a valid index (below 1). In other words, when k is too close to element 1 of blender, this "window" that you want to extract must be truncated. The same applies when k is close to the end of blender.
• For this reason, we define r that is this window, truncated when it has to.
• Now when there is truncation, blender(r) is not 21 elements long, so you cannot store it in NewMatrix(:,ii) as the : addresses the full length of NewMatrix along the first dimension, which is 21. So we have to replace : with a vector of indices whose length equals the length of r.

[Diego Yanez]

Any chance you can help with the second part of my question. I want to plot one column of NewMatrix at the time, cycling through the next column when I hit Enter. This is as far as I got:

   for j=1:length(element)
       figure(1)
       plot(NewMatrix(:,j),'bo--')
       grid on
       keyboard
   end

This stops at the first column, but does not keep cycling

[Cedric]

Yes, replace

 keyboard

with

 input( 'Press a key to continue..' )