Сделать полосовой цифровой фильтр несложно. Я только не понял данные идут сплошным потоком или нет? Если фильтровать каждый фрагмент 500 отсчетов то это не даст результата. Надо фильтровать поток с ацп без разрыва. Ещё вопрос зачем 1 мгц частоты ацп если сигнал 8 кГц?

Если через преобразование Фурье, то частотное разрешение 1/T, T - время наблюдения. Так что за 0.1 сек частотное разрешение будет 10Гц. Существуют параметрические методы оценки: Прони, Писаренко и MUSIC. Они дают более высокое частотное разрешение для процесса если известно что это гармоническое кол...

Выяснилось, что я языками не владею. Покажите спектр который подается на iffy

Проведите код чтобы посмотреть

Не совсем понятно между паразитными и полезными сигналами есть разница? Как понять это тренд надо убрать а это сигнал оставить, если все одинаково близкотк нулевой частоте.

Если однократно (по одной длинной реализации данных), то https://habr.com/ru/post/448848/
Если надо постоянно удалять (данные поступают непрерывно), то так или иначе задача сводится к построению фильтра верхних частот. Вопрос в близосте информационной компоненты к нулевой частоте.

Посмотрите здесь: http://www.dsplib.ru/content/goertzelmo ... elmod.html
Там для одной спектральной компоненты. Если надо N, то параллельно N таких блоков

Если необходимо соблюсти идентичность фильтров, то я бы строил по следующей схеме: transform.png 1) Спектр исходного сигнала (L+R)/2 - синий (L-R)/2 - красный. 2) делаем частотный перенос (L+R)/2 на нулевую частоту (обозначим этот сигнал как Z(t) запомним его, заметим что он комплексный I + jQ) 3) Ф...

Да похоже на то lear all; close all; clc; pkg load signal i = randi([0,1], 512, 1); M = 128; s = []; b = []; for k = 1:length(i) if (i(k) > 0.5) s = [s, [ones(1, M/2), -ones(1, M/2)]]; else s = [s, [-ones(1, M/2), ones(1, M/2)]]; end b = [b, i(k)*ones(1,M)]; end h = fir1(64, 0.02); s = filter(h,1,s)...

Операция смещения по частоте является нелинейной, и просто поменять ее с фильтрацией местами нельзя (менять можно только линейные операторы). Пусть имеется s(t) входной сигнал на некоторой частоте, который мы переносим по частоте в ноль и получаем z(t) = s(t) \exp(-j\omega_0 t) = \underbrace{s(t)\co...