Передискретизация Опции

[Ivan55]

Есть комлексный сигнал, ширина полосы сигнала 20 кГц, т.е. от -10 кГц до +10 кГц. Частота дискретизации сигнала 41.25 кГц, необходимо передискретизовать его в 48 кГц. Для этого необходимо проинтерполировать исходный сигнал с коэффициентом 64 и продицимировать с коэффициентом 55. Для простоты реализации это можно разбить на два каскада т.е. 8 вверх 5 вниз 8 вверх и 11 вниз. Вопрос следующий, какие должны быть параметры интерполирующего ФНЧ фильтра для каждого каскада (т.е. ширина полосы пропускания, задержания, порядок фильтра)?

[STAR_IK]

У вашего фильтра слишком широкая полоса пропускания. Им разве что в 2 раза можно продецимировать или проинтерполировать, причем сигнал с полосой не больше +/-0,37 от частоты Найквиста (Fn). Фильтр должен иметь полосу пропускания такую, чтоб пролазил спектр сигнала, и обеспечивал нужное подавление на той частоте, где будет находиться копия сигнала (интерполяция) ну или то, что завернется на спектр сигнала (децимация). Т.е. в вашем случае для интерполяции в 8 раз полоса пропускания должна быть 10 кГц, а подавление начиная с частоты 41,25-10 = 31,25 кГц и выше. Замечу что частота дискретизации при расчете фильтра уже 330 кГц, а не 41.25. К тому же т. к. частота начала заграждения 31,25 кГц меньше частоты 8/5 Fn т.е. 33 кГц, то для децецимации в 5 раз дополнительная фильтрация не нужна. Во втором каскаде можно тоже обойтись одним фильтром, но подавление должно начинаться уже с частоты не выше 48 кГц-10 кГц=38 кГц, при частоте дискретизации (8/5 fs1)*8 = 528 кГц. И зачем такое большое подавление в 153 дБ??? По моему это перебор, хватит и 70 - 90 дБ. Кстати, использование fdatool в матлабе гораздо упрощает процедуру расчета фильтров. И не забывайте про полифазную реализацию фильтров, которая уменьшит вычислительные затраты

[Ivan55]

У вашего фильтра слишком широкая полоса пропускания. Им разве что в 2 раза можно продецимировать или проинтерполировать, причем сигнал с полосой не больше +/-0,37 от частоты Найквиста (Fn). Фильтр должен иметь полосу пропускания такую, чтоб пролазил спектр сигнала, и обеспечивал нужное подавление на той частоте, где будет находиться копия сигнала (интерполяция) ну или то, что завернется на спектр сигнала (децимация). Т.е. в вашем случае для интерполяции в 8 раз полоса пропускания должна быть 10 кГц, а подавление начиная с частоты 41,25-10 = 31,25 кГц и выше. Замечу что частота дискретизации при расчете фильтра уже 330 кГц, а не 41.25. К тому же т. к. частота начала заграждения 31,25 кГц меньше частоты 8/5 Fn т.е. 33 кГц, то для децецимации в 5 раз дополнительная фильтрация не нужна. Во втором каскаде можно тоже обойтись одним фильтром, но подавление должно начинаться уже с частоты не выше 48 кГц-10 кГц=38 кГц, при частоте дискретизации (8/5 fs1)*8 = 528 кГц. И зачем такое большое подавление в 153 дБ??? По моему это перебор, хватит и 70 - 90 дБ. Кстати, использование fdatool в матлабе гораздо упрощает процедуру расчета фильтров. И не забывайте про полифазную реализацию фильтров, которая уменьшит вычислительные затраты

Сделал так как вы сказали, но результата особого не дало, на 8 дБ сигнал ухудшается, когда делаю встроенной функцией матлаб resampling все работает идеально, странно как то... Но не падает даже на 1 дБ

Сообщение отредактировал Ivan55 - Feb 13 2013, 14:41

[Make_Pic]

Сделал так как вы сказали, но результата особого не дало, на 8 дБ сигнал ухудшается, когда делаю встроенной функцией матлаб resampling все работает идеально, странно как то... Но не падает даже на 1 дБ

Уже наступали на эти грабли - в матлаб точность расчетов выше, а вы не учитываете переполнение в своем железе/программе.

[Ivan55]

Уже наступали на эти грабли - в матлаб точность расчетов выше, а вы не учитываете переполнение в своем железе/программе.

Я делаю в матлабе без железа