How to vectorize this code? element wise to save time

Question

MIDO 2020-4-19

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/519187-how-to-vectorize-this-code-element-wise-to-save-time

关闭： MATLAB Answer Bot 2021-8-20

% size(matAL) > 32,500,5
matCT=zeros(500,32,32); %Pre-allocate
x=0;
for t=1:500
    x=1+x;
    for i=1:32
        for j=1:32
            xsd = squeeze(matAL(i,t,:) .* conj(matAL(j,t,:)));
            matCT(x,j,i) = abs(mean(exp(1i*angle(xsd))));
        end
    end
end

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

此问题已关闭。

Answer 1

Vimal Rathod 2020-4-26

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/519187-how-to-vectorize-this-code-element-wise-to-save-time#answer_428425

在 MATLAB Online 中打开

Hi,

As all the three for loops have some dependencies it is difficult to vectorize the code.

you can try the following code, it is partially vectorized.

for i = 1:32
for j = 1:32
xsd = squeeze(matAL(i,1:500,:) .* conj(matAL(j,1:500,:)));
matCT(1:500,j,i) = abs(mean(exp(1i*angle(xsd(1:500,:)))'));
end
end