MATLAB: Automatic summation of sub matrices

sub matricessumsummation

So I have a 16×200 matrix of 0's and 1's, I have

t = any(Y,1); fc = strfind(t, [0 1]) + 1; lc = strfind(t, [1 0]);

entered giving me the bookends of my sub matrices under fc and lc. How would I enable MATLAB to automatically sum the 1's in all of the 16 columns of whatever widths were dictated by the previous formula, and give me the sum of each sum matrix. This would be assuming that every time I run this there would be a different number of sub matrices in the 16×200 matrix of different sizes.

Thank you.

Best Answer

Sean's solution to your earlier question was fine, and the logic for the above would be much the same as the summing logic that Sean uses.

Related Solutions

MATLAB: Dividing out smaller matrices bracketed by columns of just 0’s

t = any(YourMatrix,1);
fc = strfind(t, [0 1]) + 1;
lc = strfind(t, [1 0]);

Then each pair fc(K), lc(K) gives the columns to sum submatrix #K over.

You may wish to consider making the above logic more robust to take into account the possibility that a submatrix starts at the very beginning of YourMatrix or ends at the very end of YourMatrix.

MATLAB: Extract and store valuable data in vectors

Simple thresholding won't work because in each burst there are dips that you probably want to keep. Luckily for you I solved this same problem for another person in this discussion: http://www.mathworks.com/matlabcentral/answers/229678#answer_185935 :

% Setup - read in the signal and plot it.
s=load('signal.mat')
signal = s.coarse_d;
subplot(3,1,1);
plot(signal, '-', 'Color', [0,.5,0]);
grid on;
% MAIN CODE IS THE FOLLOWING TWO LINES OF CODE.
% Threshold the signal
lowSignal = abs(signal) < 0.1;
% Remove stretches less than 10,000 elements long.
% And invert it to get the high signal areas instead of the quiet areas.
lowSignal = ~bwareaopen(lowSignal, 10000);
% Now we're done.  Plot the results.
subplot(3,1,2);
plot(lowSignal, 'b-', 'LineWidth', 2);
grid on;
% Plot both on the same graph now.
subplot(3,1,3);
plot(signal, '-', 'Color', [0,.5,0]);
hold on;
plot(lowSignal, 'b-', 'LineWidth', 2);
grid on;

As you can see it identifies each "block" but leaves the inside of the block intact.

Best Answer

Related Solutions

MATLAB: Dividing out smaller matrices bracketed by columns of just 0’s

MATLAB: Extract and store valuable data in vectors

Related Question