How to vectorize this code? element wise to save time

MIDO

2020 4 19

1 个回答

更新时间：2021 8 20

1 次查看（30 天）

Follow Question

信息

此问题已关闭。请重新打开它进行编辑或回答。

显示更早的评论

在 MATLAB Online 中打开

0 个投票

% 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更早的评论

Follow Question

回答（1 个）

Vimal Rathod 2020-4-26

在 MATLAB Online 中打开

0 个投票

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