Спектральный анализ на сверхнизких частотах Опции

[Crowbar]

Допустим, требуется получить разложение спектра частот до 20Гц с точностью 0,01Гц и выше. Каким образом это реализуется, помимо самого простого способа, как поставить частоту отцифровки на 40Гц, выставить кол-во сэмплов на 4000 и ждать больше полутора минут завершения очередного цикла?

[mikalaha]

Допустим, требуется получить разложение спектра частот до 20Гц с точностью 0,01Гц и выше. Каким образом это реализуется, помимо самого простого способа, как поставить частоту отцифровки на 40Гц, выставить кол-во сэмплов на 4000 и ждать больше полутора минут завершения очередного цикла?

Не обязательно. Просто набираете некоторое количество отсчетов (например 80 - 2 секунды записи). Затем производите БПФ требуемой размерности (например, 4096 - точность как раз около 0.01Гц). Недостающие до 4096 отсчеты (4016) заполняете нулями. И выполняете БПФ.

[Alex255]

Не обязательно. Просто набираете некоторое количество отсчетов (например 80 - 2 секунды записи). Затем производите БПФ требуемой размерности (например, 4096 - точность как раз около 0.01Гц). Недостающие до 4096 отсчеты (4016) заполняете нулями. И выполняете БПФ.

То есть Вы хотите получить информацию о процессе периодом 100сек за 2секунды? Увы...

To Crowbar: Да, именно так как Вы сказали

[mikalaha]

То есть Вы хотите получить информацию о процессе периодом 100сек за 2секунды? Увы...

To Crowbar: Да, именно так как Вы сказали

На самом деле все отлично работает.
Например:
сигнал - тон с частотой 10 Гц (период 0.1 сек) - в одной секунде 10 периодов!!
при записи с дискретизацией 40 Гц одной секунды получаем 40 отсчетов, в которых содержится
целых 10 периодов тона.
Если теперь дополнить нулями и взять FFT мы получим этот тон.

Прикреплены картинки
time - сигнал во временной области;
freq - после FFT 4096
freq2 - увеличенная частотная позиция 10 Гц после FFT (что и требовалось)

На самом деле все отлично работает.
Например:
сигнал - тон с частотой 10 Гц (период 0.1 сек) - в одной секунде 10 периодов!!
при записи с дискретизацией 40 Гц одной секунды получаем 40 отсчетов, в которых содержится
целых 10 периодов тона.
Если теперь дополнить нулями и взять FFT мы получим этот тон.

Прикреплены картинки
time - сигнал во временной области;
freq - после FFT 4096
freq2 - увеличенная частотная позиция 10 Гц после FFT (что и требовалось)

То же самое наблюдается и при тоне с частотой например (10.05 Гц). Единственное различие вышло следующее: без прореживания нулями пик после Фурье вышел на частоте 10.0488 Гц, а с прореживанием - 10.0399 Гц

Эскизы прикрепленных изображений