Отправлено Mega 10 февраля 2004 г. 09:40

Нужна теоретико-практическая помощь. В матлабе генерю несколько синусоид различной частоты с различными начальными фазами. Делаю FFT, затем нахожу через ANGLE фазы. Для тех синусоид которые не вписываются в бины, а болтаются между ними, измеренные фазы отличаются от истинных значений, причем иногда очень сильно. Как вооще можно оценивать фазу синусоиды частота которой не кратна бинам, и соответственно имеет отклик на двух бинах. У меня получалось так примерно: начальная частота f с фазой 0.7, после FFT дала отклик на бинах n и n+1 (ну и на других, но они достаточно слабые), и каждый этот бин имеет фазу 0.4. Т.е ошибка получилась 0.3. При тех же самых условиях берешь FFT большего размера, ошибка не уменьшается, а по разному - то больше, то меньше. Сигнал/шум ломовой, так как все чисто в модели.
Буду благодарен любой подсказке, но особенно подкрепленным конкретикой.

• Ответ: Всем спасибо! Все получилось. — Mega   (11.02.2004 09:54, пустое)
• забавная штука - но все повторяется в этой жизни - если не сложно дайте определение фазы которую вы хотите вычислить? ;-) — net_net   (10.02.2004 14:45, пустое)
  • Ответ: В вопрос не въехал... — Mega   (10.02.2004 23:51, 93 байт)
  • где-то это уже было :) arg(signal(0)) — SM   (10.02.2004 15:50, пустое)
    • у человека то явно не так ;-) поскольку он из сигнала выделет нечто и потом говорит это фаза - и она должна совпадать с фазой другого - и тут совсем не arg(sin(x)) ж-) — net_net   (10.02.2004 22:40, пустое)
      • Судя по его вопросу там все понятно (+) — SМ   (11.02.2004 01:28, 341 байт)
        • Ответ: — net_net   (11.02.2004 10:51, 327 байт)
          • А про "фазу всего сигнала" вроде ничего и не было. — SM   (11.02.2004 11:12, пустое)
            • я про определение - говорю. причем тут всего сигнала? — net_net   (11.02.2004 23:01, пустое)
• Ответ: — GroundCtrl   (10.02.2004 10:31, 778 байт)
• Ответ: — GroundCtrl   (10.02.2004 10:31, 778 байт)
• Если Вы знаете, что у вас 3 частоты F1, F2, F3 - то вам никакой FFT не нужен. делаете корреляцию с sin(F1) & cos(F1) - получаете проекции на x и y - по ним находите угол вектора, повторяется для F2 и F3 — AntZ   (10.02.2004 10:07, пустое)
  • Ответ: Количество частот изначально неизвестно, номиналы частот неизвестны и оцениваются по ходу дела ... — Mega   (10.02.2004 10:35, 167 байт)
    • тогда с помошью FFT находите грубо, затем делаете интерполяцию, затем находите экстремумы, затем каждый экстремум коррелируете с sin(Fext) и cos(Fext) — AntZ   (10.02.2004 10:39, пустое)
      • Интерполяцию спектра, я имел в виду — AntZ   (10.02.2004 10:40, пустое)
  • и какая точность данного метода - если частоты отличаются на 10-ки герц - будет работать? — 99600036   (10.02.2004 10:13, пустое)
    • Корреляция c sin(2piF) и cos(2piF) иногда называется Single Point FFT, ну есть еще алгоритм Герцеля, который Вас может сильно выручить. FFT на N точек - есть не что иное, как N корреляций c cos(2*pi*N*(Fs/N)) sin(2*pi*N*(Fs/N)) — AntZ   (10.02.2004 10:54, пустое)
      • по поводу Герцеля (+) — 99600036   (10.02.2004 14:10, 215 байт)
    • Будет. Если достаточно долго коррелировать. — SМ   (10.02.2004 10:25, пустое)