Простое умножение матриц. Неверный результат

Dopp_shift = (-10E3 : 200 : 10E3)';
k = (2 * pi) .* Dopp_shift;
u = 1i.*(0:1:length(etal)-1);
w = k * u;

clc;
clear;
Fs = 24.552E6;      % Частота дискретизации в Гц
T_track = 1E-3; % Длительность периода в секундах
T_samples = T_track * Fs;    % Длительность периода в отсчетах
NeedSamples = T_samples * 5;
 
FileName = '.\\1575420__24_MHz__0.dat';
NumOfShiftedSamples = 0; % Считываем фрагмент с самого начала записи
HeadByteSize = 0;
DataType = 'int16';
NumOfChannels = 1;
ChanNum = 0;
Signal = ReadSignalFromFile(FileName, HeadByteSize, ...
    NumOfChannels, ChanNum, NumOfShiftedSamples, NeedSamples, ...
    DataType);
%FsNew = 2.046E6;
etal = repmat( GenCACode(2, 1), 2, 1 ); 
etal = reshape( etal, 1, 2 * length(etal) ); %считали и продублировали
etal = etal .* 2 - 1; %-1,1
Signal = resample( Signal, 1, 12 );
Signal = repmat(Signal, 101, 1);
Dopp_shift = (-10000 : 200 : 10000)';
k = (2 * pi) .* Dopp_shift * (0:1:length(etal)-1);
shift = exp(complex(0, k));
etal = repmat(etal, length(Dopp_shift), 1);
etal = etal .* shift;
A = zeros(length(Dopp_shift), length(Signal)*2-1);
for i=1:1:length(Dopp_shift)
    A(i,:) = xcorr( etal(i,:), conj(Signal(i,:)) );
end
R = abs(A);