Недостаточно памяти с использованием bsxfun MATLAB

Я пытаюсь реализовать сжатие изображений с помощью преобразования Берроуза-Уилера. Рассмотрим одномерную матрицу сканирования пути:

p = [2 5 4 2 3 1 5];

а затем применить преобразование Берроуза-Уилера:

function output = bwtenc(p) 
n = numel(p);
x = zeros(length(p),1);
for i = 1:length(p)
    left_cyclic = mod(bsxfun(@plus, 1:n, (0:n-1).')-1, n) + 1;
    x = p(left_cyclic);
end
[lex ind] = sortrows(x);
output = lex(:,end);
output = uint8(output(:)');
end

И это работает! Но проблема в том, что когда я пытаюсь реализовать 1D-матрицу из Lena.bmp размером 512 * 512, появляется сообщение об ошибке, показывающее, что bsxfun не хватает памяти. Кто-нибудь, пожалуйста, помогите мне.