Array prealocation - Extra array of zeros in 3D array

Question

Christopher McCausland 2021-11-2

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1577520-array-prealocation-extra-array-of-zeros-in-3d-array

评论： Christopher McCausland 2021-11-2

Hi,

I have a 13 channel signal recorded over several hours. I wish to take all channels, split them into 10 second windows with a 5 second overlap and finally concatinate all channels into a 3D array (Buffer length(n), number of buffers, channels).

I have written the function below which works as intended, except my preallocate method results in an 'extra channel' of zeros, 14 rather than 13 with the 1st array containing zeros. What I think would work is windows3D = zeros(size(windows)), however I cannot call that as 'window' is generated after the pre-allocation call.

I could just go into the array after the fact and 'delete' the additonal array of zeros however this feels like cheating, is there a better way to do this?

Kind regards,

Christopher

function [windows3D] = windowGen(eegDataNorm)
% 10 second window, 50% overlap 
EEG = eegDataNorm; % Short label
Fs = 200; % Sampling frequency (Hz)
windows3D = zeros(2000,5147); %Preallocate
for i = 1:size(eegDataNorm)
    windows = buffer(eegDataNorm(i,:),10*Fs,5*Fs);
    windows3D = cat(3,windows3D,windows);
end
end

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

请先登录，再进行评论。

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

Answer 1

Stephen23 2021-11-2

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1577520-array-prealocation-extra-array-of-zeros-in-3d-array#answer_822390

编辑：Stephen23 2021-11-2

在 MATLAB Online 中打开

windows3D = nan(2000,5147,0);

or use actual preallocation:

N = size(eegDataNorm,1);
windows3D = zeros(2000,5147,N);
for k = 1:N
    windows3D(:,:,k) = buffer(..)
end

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Christopher McCausland 2021-11-2

在 MATLAB Online 中打开

Hi Stephen,

I just wanted to say thank you for you answer, I think I may have been looking at my code for too long... not sure how I didn't see this myself! For anyone else here's what the final fuction looks like. I also added a slight modifcation to allow diffrent lenght data to be used which works nicely.

function [windows3D] = windowGen(eegDataNorm)
% 10 second window, 50% overlap 
EEG = eegDataNorm; % Short label
Fs = 200; % Sampling frequency (Hz)
C = size(eegDataNorm,1);
[N,P] = size(buffer(eegDataNorm(1,:),10*Fs,5*Fs));
windows3D = zeros(N,P,C); %Prealocate
for i = 1:C
    windows3D(:,:,i) = buffer(eegDataNorm(i,:),10*Fs,5*Fs);
end
end