Отправлено pvr 24 ноября 2006 г. 17:32

У меня есть периодический процесс с некоторыми искажениями. Мне нужно определить пару десятков его гармоник при условии, что 2048 точек точно не укладываются в N периодов. Как это делается? (О БПФ и окне Блэкмана я знаю).

• Оказалось просто (+) — pvr   (26.11.2006 18:49 217.197.113.41, 396 байт)
  • Амлитуду : — -=ВН=-   (27.11.2006 11:56 217.25.80.230, 238 байт)
  • А какая точность получается? У меня -80 дБ на окне Кайзера. Или даже более того. — st256   (26.11.2006 19:46 217.151.231.218, пустое)
    • Она фактически определяется разрядность АЦП. Когда я моделирую сигнал преобразователем бесконечной разрядности, то вижу разницу разве что в 7-м десятичном разряде. — pvr   (26.11.2006 21:27 217.197.113.41, пустое)
      • А гармоники у Вас далеко отстоят друг от друга? — st256   (27.11.2006 07:33 217.151.231.218, пустое)
        • Примерно каждая 16-я частота БПФ - гармоника. — pvr   (27.11.2006 11:21 81.95.129.82, пустое)
          • Тогда понятно. Я и не знал это свойство окна Блэкмана. — st256   (27.11.2006 14:03 217.151.231.218, 240 байт)
• Частоты гармоник какие? — st256   (25.11.2006 07:53 217.151.231.218, пустое)
• Попробуйте такой способ (+) — Бяка   (24.11.2006 19:08 80.253.9.6, 396 байт)
• Что-то я ни хрена не понял. В чем экзотика-то? БПФ с окном и просится. Причем настойчиво. — -=ВН=-   (24.11.2006 18:55 193.125.71.140, пустое)
  • Делаю БПФ с окном, получаю БПФ'овские гармоники (БПФ считает, что периоды точно укладываются в 2048 точек). А мне нужны гармоники основной частоты. — pvr   (24.11.2006 19:02 217.197.113.41, пустое)
    • Что-то Вы не то всё пишете. — -=ВН=-   (24.11.2006 19:46 193.125.71.140, 1459 байт)
• Интерполяцией не получается. у меня есть две соседние гармоники с примерно одинаковыми амплитудами, а нужная мне гармоника - между ними, и ее амплитуда раза в полтора больше. — pvr   (24.11.2006 18:46 217.197.113.41, пустое)
  • так интерполируйте sinc'ом :) Он Вам и выдаст ткой пик. — SM   (24.11.2006 19:22 213.141.159.26, пустое)
• А если вначале корреляционным методом определить период, а потом поставить гребенку БИХ-фильтров (или один КИХ-фильтр), так, чтобы пики попадали как раз в частоты k/T, где Т - найденный период период, k - целое число. — homekvn   (24.11.2006 18:42 212.185.161.237, пустое)
  • А что почитать? — pvr   (24.11.2006 19:16 217.197.113.41, пустое)
    • Ничего не надо читать. Надо взять для начала Матлаб и (+) — homekvn   (24.11.2006 19:37 212.185.161.237, 1352 байт)
      • Добавлю, что при таком методе, чем длиннее будет выборка, тем меньше будут поганить жизнь посторонние сигналы, которые присутствуют в смеси с Вашим периодическим сигналом. — homekvn   (24.11.2006 19:41 212.185.161.237, пустое)
• Интерполяцией например делается между этими точками. — SM   (24.11.2006 18:37 213.141.159.26, пустое)
  • Когда периоды процесса не укладываются точно в 2048 точек, то гармоники основной частоты размазываются по БПФ'овским гармоникам. И как восстановить пик интерполяцией? — pvr   (24.11.2006 19:06 217.197.113.41, пустое)
    • Ну частоту-то Вы точно знаете (+) — SM   (24.11.2006 19:21 213.141.159.26, 320 байт)