1 лр / myCF

function mKF = myCF(M1, M2) % функция для получения корелляционной функции
    M1inv = M1(length(M1):-1:1); % разворачиваем М1 (ну короче я мыслю, как на бумаге)
    t = zeros(1, length(M1)); % временная ячейка t - память триггеров
    mKFsize = length(M2)*2-1;  % количество итераций для получения КФ
    mKF = zeros(1, mKFsize); % создание заготовки под КФ
    for j = 1:(mKFsize)  
        for i = (length(M1)):-1:2 %сдвигаем в триггерах элементы на 1 вперед, 
                                  %а на место первого итерации идут с последнего до второго элементов
            t(i)=t(i-1);
            if j>length(M2) %если нам уже сигнал закончился, то добавляем 0
                t(1)=0;
            else
                t(1)= M2(j); %если нет - то берем дальше
            end
        end
        mKF(j) = sum(M1inv.*t);  %перемножаем и закидываем в житый элемент КФ
    end
end