Добрый день.
Допустим есть два синусоидальных сигнала одинаковой частоты (порядка 30Гц).
После оцифровки в DSP получаем два массива A и B с дискретными значениями этих сигналов.

Теперь вопрос - годится ли следующий алгоритм для вычисления разности фаз этих сигналов?

1. Берем первый массив A, осуществляем над ним преобразование Гильберта, в результате чего получаем еще один массив A2, который будет представлять собой этот же самый сигнал, но сдвинутый на 90 градусов.
2. Берем первое значение из массива A2 и делим его на первое значение из массива A. Из полученного частного находим арктангенс, т.е. : phi_1 = arctg(A2[1]/A[1]). Таким образом получаем мгновенную фазу первого сигнала.
3. Тоже самое проделываем со вторым сигналом и получаем его мгновенную фазу phi_2.
4. Вычисляем phi_1 - phi_2. Это и будет разностью фаз. Для большей точности можно вычислить некоторое количства фаз для разных индексов массива и усреднить их.