I think you're asking two questions:
a) Why is the indexing broken?
b) What is a better way to do this?
I'll ask you my own: Are you a C programmer? =) There seems to be a confusion over 0-based and 1-based indexing in this code.
For example:
Time_delay(d+1)
This will overrun on the last element, and I'm sure you don't mean to add 1 to the index.
Next:
rem(d,8)+1
Is this the desired behaviour?
The sequence will go:
2, 3, 4, 5, 6, 7, 8, 1, 2, 3, ...
Did you mean:
rem(d-1,8)+1
Now, the matrix out. Presumably you want to be putting a single energy value into it after processing each row (ie every 8 times around). So you would end up with a 1x10000 vector. At the moment you are just overwriting it each time.
It looks to me like this code would be a lot clearer if you used nested for-loops over the columns and rows. Perhaps like so:
% Initialise output vector with same number of columns as time delay.
out = zeros(1,size(Time_delay,2));
for col = 1:size(Time_delay,2)
% Zero the accumulator and copy first (non-delayed) row.
mic_out = zeros(size(mic_out));
mic_out(1,:) = mic_in(1,:);
% Accumulate delayed signals for all other rows
for row = 2:size(Time_delay,1)
mic_out(row,:) = del( mic_in(row,:), Time_delay(row,col) );
end
% Store energy value for current set of delays
out(col) = sum( sum(mic_out,1).^2, 2 );
end
I suspect that the calculation for out might not be quite the same as you described. But maybe the above code makes it more obvious where your mistakes are, so you can fix them more easily.
Cheers, -g-
