GS
отдаю в функцию взаимокорреляции два набора точек
Size: a a a
2018 July 30
GS
и получаю, что нужно второй (красный) график сдвинуть вправо на 17 * step
GS
график функции взаимокорреляции
GS
и вот, как она мне говорит их надо сдвинуть
g
[R, lag] = xcorr(y2, y(lagDiff:end));
[~,I] = max(abs(R));
lagDiff = lag(I)
plot(lag, R)
plot(x, y, x2 - lagDiff * (right - left) / count, y2)
lagDiff = 17
Почему
max(abs(R))?GS
честно - не знаю. Смотрю пример из матЛаба по xcorr - Delay between two correlated signals - http://www.mathworks.com/help/signal/ref/xcorr.html
GS
и у них всё работает, естественно )
GS
и графики там куда "мощнее" моего синтетического
GS
Почему
max(abs(R))?и, кстати, конкретно в этом случае - без разницы
GS
просто максимум совпадает с макс(абс)
g
Сдвинул, наверное, неправильно. Сейчас, почитаю доку.
GS
давай я свои соображения расскажу
GS
функция эта ничего не знает об абсциссах последовательностей
GS
и, соответственно, результат её работы - это относительный сдвиг одной последовательности от другой по абсциссам одной из них
GS
вот, что я законспектировал, когда читал документацию -
Когда мы найдём максимум взаимокорреляции, то получим некий сдвиг.
Это сдвиг, по формуле, значит, насколько сдвинут по индексам
первый массив.
Если этот сдвиг отрицательный, то при расчёте в формуле
по отрицательным индексам используются нули. Иными словами,
первый график будет дополняться слева нулями (x[i + k]; 0).
Ещё раз иными словами - график сдвигается вправо, а слева
"приклеивается" прямая y = 0.
Если сдвиг положительный: первый график (x) дополняется нулями
справа. Иными словами, сдвигается влево.
GS
что-то я сейчас перечитываю эти строки и осознаю, что вообще картина должна быть другой
GS
мне надо для y2 последовательности взять абсциссы y, со сдвигом на -lagDiff * step
GS
совсем запутался (
g
@geshep Да, это занятно. Проверь второй по величине положительный максимум.
GS
всё таки я ошибся