Неравномерные отсчеты и DSP
Jan. 27th, 2011 11:07 am![[personal profile]](https://www.dreamwidth.org/img/silk/identity/user.png){kind=small}
metaclassА вот скажите, как правильно делать FFT, цифровую фильтрацию и тому подобное, если отсчеты сигнала идут неравномерно во времени? Вся теория, которую я когда либо изучал, обычно заточена под постоянный период отсчетов.
В принципе, можно было бы в лоб интерполировать равномерные отсчеты, исходя из неравномерных данных, но может есть еще какие-нибудь вуду-оккультные методы?
В принципе, можно было бы в лоб интерполировать равномерные отсчеты, исходя из неравномерных данных, но может есть еще какие-нибудь вуду-оккультные методы?
![[identity profile]](https://www.dreamwidth.org/img/silk/identity/openid.png){kind=small}
no subject
Date: 2011-01-27 09:16 am (UTC)FFT, скорее всего, не получится в принципе (там вся оптимизация основана на регулярности отсчётов), а вот преобразование Фурье проблем не вызовет. FIR/IIR станут посложнее, чем для равномерных отсчётов (коэффициенты придётся вычислять каждый раз, а не брать из таблицы), но терпимо.
no subject
Date: 2011-01-27 09:30 am (UTC)no subject
Date: 2011-01-27 09:41 am (UTC)no subject
Date: 2011-01-27 09:46 am (UTC)будет как фильтр отсекающий очень высокие частоты - в районе предельных для FFT, которые все равно вряд ли удастся засечь
no subject
Date: 2011-01-27 09:50 am (UTC)no subject
Date: 2011-01-27 09:53 am (UTC)поделить на нольустроить децимацию и тупо отрезать кусок спектра фильтром, и дальше работать с уменьшенной частотой отсчетов - один хрен в высоких частотах мусор.no subject
Date: 2011-01-27 01:59 pm (UTC)no subject
Date: 2011-01-27 02:05 pm (UTC)Если этот фильтр ещё и высокие частоты убирать будет, тогда даже интуитивно будет похож на интерполятор.
Но только, с неравномерными отсчётами, спектр этого фильтра получится ну очень непохож на отсекание по частоте.
no subject
Date: 2011-01-27 09:35 am (UTC)Data": http://www.mpi-hd.mpg.de/astrophysik/HEA/internal/Numerical_Recipes/f13-8.pdf
и еще десять верхних ссылок с гугла. FAST версия тоже имеется.
no subject
Date: 2011-01-27 09:45 am (UTC)http://window.edu.ru/window_catalog/pdf2txt?p_id=28533&p_page=1 но там в основном расчет автокорелляций всяческих.
Все что я для себя оттуда вынес - это мрак и ужос с червями, если нет крайней необходимости, лучше этим не заниматься. Если, например, сетка отсчетов приблизительно равномерна, а имеется просто некоторое дрожание то, может, и фиг с ним?
no subject
Date: 2011-01-27 02:08 pm (UTC)Если базовые всплески имеют компактный носитель, что в данном случае, попросту говоря, означает, что на пределами некой области они обращаются в ноль, или же, для практики достаточно, чтобы были достаточно малы, как какие-нибудь гауссианы.
no subject
Date: 2011-01-27 06:39 pm (UTC)Может имеет смысл тупо свёртывать с (ко)синусоидами? Метод, формально, процессороёмче, но на практике выкидывается интерполяция, то есть, надо аккуратно считать сложность. Или не заморачиваться сложностью, если явно есть запас производительности.
no subject
Date: 2011-01-27 06:50 pm (UTC)no subject
Date: 2011-01-27 06:53 pm (UTC)