Добрый день. 1. Если говорить о стандартное задаче обнаружения сигнала, то существует схема оптимального алгоритма, на рисунке пример. Т.е. в качестве фильтра используется сам сигнал. Снимок.PNG Я промоделировал грубо схему. На графиках сверху входной сигнала, снизу сигнал поступающий на решающее ус...

Есть обсуждение об оценки длительности импульса с шумом. https://dsp.stackexchange.com/questions/20143/duration-of-unknown-rectangular-pulse-with-additive-white-gaussian-noise Границы импульса выделяют по порогу (можно по порогу производной от сигнала). Найдя границы импульсов, можно определить и вс...

В силу некоррелированности E_{s}=\int s^{2} dt = \int \left ( x + n \right ) ^ {2} dt = \int (x^2 + n^2 + 2xn) dt = E_x + E_n + 2B_{xn} = \left \langle B_{xn} = 0 \right \rangle = E_x + E_n Тогда SNR = \frac{E_x}{E_n} = \frac{E_s-E_n}{E_n} Что то не так в рассуждении? Нам ведь известны необходимые п...

А разве все так сложно?
Если шум белый, то сигнал некоррелирован с шумом и взаимная энергия ~= 0,
тогда энергия смеси ~= энергия сигнала + энергия шума.

Т.е ОСШ ~= (энергия смеси - энергия шума) / энергия шума

Изменил расчет невалидной части передискретизированного сигнала: invalid_out_length = fix((up_rate*(phase_length-1) - (ord/2)) / down_rate); if rem((up_rate*(phase_length-1) - (ord/2)), down_rate) invalid_out_length = invalid_out_length + 1; end Теперь идея такая - считаю длину нулевой вставки после...

фильтр с частотой среза 1/5 то и генерируется фильтр у которого каждый 5 отсчёт импульсной характеристики 0

Не обязательно же, можно сформировать и другой фильтр с такой частотой среза, у которого не будет такого количества нулей в ИХ. Эффективность реализации от этого не снизится.

resample(x, p, q, n) генерирует фильтр, а дальше вызывается % ---- HERE'S THE CALL TO UPFIRDN ---------------------------- y = upfirdn(x,h,p,q); которой не важны частотные характеристики фильтра. Порядок выбирается достаточно произвольно, в справке написано, что увеличение n дает лучшую точность и б...

up_rate = 1001; down_rate = 1000; Эти коэффициенты как пример неправильной работы алгоритма, в данном случае не важна эффективность, просто хотелось бы разобраться с полифазной передискретизацией. Однако отмечу, что единственной проблемой является хранение и доступ к элементам фильтра такого порядка...

Сигнала input длины input_length фильтруется полифазным фильтром polyphase_filter для передискретизации в up_rate/down_rate раз. Добиваю input в начале и конце нулями по phase_length-1 штук (phase_length = polyphase_filter / up_rate). Т.е получаю на выходе фильтра отсчеты с первого "ненулевого". Хоч...

Я потихоньку пытаюсь сделать передискретизацию в DSPL c КИХ фильтрами и полифазной фильтрацией.
Но есть трудности с отладкой без среды. Как вы отлаживаете? Если через gdb, то надо менять make фалы?