Тема нерекурсивные частотные цифровые фильтры недостаточно овладеть премудростью, нужно уметь пользоваться ею (стр. 3 из 5)
N = 2p/Dp. (7.3.2)
Для примера на рис. 7.3.1 значение N принято равным 200, при этом крутизна переходной зоны увеличилась (тонкая кривая H'(w), N=200), создавая запас на последующее сглаживание весовой функцией.
Выбор весовых функций целесообразно осуществлять по допустимой величине осцилляций усиления сигнала в полосе подавления, т.е. по относительному значению амплитуды первого выброса на передаточных характеристиках весовых функций. Для выбранной весовой функции (с учетом числа ее членов по (7.3.2)) производится расчет весовых коэффициентов pn, после чего устанавливаются окончательные значения оператора фильтра:
hn = pn h(n). (7.3.3)
 Рис. 7.3.2. Полосовая фильтрация (вверху – входной сигнал, внизу – выходной). |
Подстановкой коэффициентов (7.3.3) в (7.3.1) рекомендуется произвести построение полученной передаточной характеристики фильтра и непосредственно по ней оценить пригодность фильтра для поставленных задач. Это наглядно видно на рис. 7.3.1, где для нашего примера была применена весовая функция Гаусса. Передаточная функция Hp(w) имеет практически такую же крутизну, как и функция H'(w) при N=100 и практически плоскую вершину в интервале спектра сигнала. Качество работы фильтра для сигнала, приведенного на рис. 7.2.1, можно видеть на рис. 7.3.2.
При необходимости более точной оценки полученной передаточной функции можно рекомендовать увеличение ее частотного разрешения в 2-4 раза перед выполнением преобразования Фурье, что можно выполнить путем увеличения размеров оператора hn дополнением нулями.
Основные весовые функции. Ниже в таблицах приведены формулы и основные спектральные характеристики наиболее распространенных весовых окон. Носители весовых функций, в принципе, являются неограниченными и при использовании в качестве весовых окон действуют только в пределах окна и обнуляются за его пределами. Для упрощения записи формулы приводятся в аналитической форме с временным окном 2t, симметричным относительно нуля (0
t). При переходе к дискретной форме окно 2t заменяется окном 2N+1, а значения t – дискретами t = nDt. Большинство весовых функций на границах окна (n = N) принимают нулевые или близкие к нулевым значения. Последнее исключается, если принять 2t = (2N+3)Dt, при этом близкие к нулю значения перемещаются за границы окна.Основные весовые функции.
| Временное окно | Весовая функция | Фурье-образ |
| Естественное (П) | П(t) = 1, |t|£t; П(t) = 0, |t|>t | П(w) = 2t sinc[wt] |
| Бартлетта (D) | b(t) = 1-|t|/t | B(w) = t sinc2(wt/2). |
| Хеннинга, Ганна | p(t) = 0.5[1+cos(pt/t)] | 0.5П(w)+0.25П(w+p/t)+0.25П(w-p/t) |
| Хемминга | p(t) = 0.54+0.46 cos(pt/t) | 0.54П(w)+0.23П(w+p/t)+0.23П(w-p/t) |
| Карре (2-е окно) | p(t) = b(t) sinc(pt/t) | t·B(w)*П(w), П(w) = 1 при |w|<p/t |
| Лапласа-Гаусса | p(t) = exp[-b2(t/t)2/2] | [(t/b)  exp(-t2w2/(2b2))] ③ П(w) |
| Кайзера-Бесселя | p(t) =  , Jo[x] =  [(x/2)k/k!]2 | Вычисляется преобразованием Фурье. Jo[x] - модифицированная функция Бесселя нулевого порядка |
Характеристики спектров весовых функций.
| Параметры | Ед. изм. | П- Окно | Барт- летт | Лан-цош | Хен- нинг | Хемминг | Кар- ре | Лаплас | Кайзер |
| Амплитуда: Главный пик 1-й выброс(-) 2-й выброс(+) Ширина Гл. пика Положения: 1-й нуль 1-й выброс 2-й нуль 2-й выброс | t %Гл.п. - “ - wt/2p wt/2p wt/2p wt/2p wt/2p | 2 0.217 0.128 0.60 0.50 0.72 1.00 1.22 | 1 - 0.047 0.89 1.00 - - 1.44 | 1.18 0.048 0.020 0.87 0.82 1.00 1.29 1.50 | 1 0.027 0.0084 1.00 1.00 1.19 1.50 1.72 | 1.08 0.0062 0.0016 0.91 1.00 1.09 1.30 1.41 | 0.77 - - 1.12 - - - - | 0.83 0.0016 0.0014 1.12 1.74 1.91 2.10 2.34 | 0.82 .00045 .00028 1.15 1.52 1.59 1.74 1.88 |
Весовая функция Кайзера. Наибольшее распространение при расчетах частотных НЦФ получила весовая функция Кайзера:
p(n) =
.Это объясняется тем, что параметры функции Кайзера могут устанавливаться непосредственно по техническим требованиям к передаточным функциям проектируемых фильтров – допустимой ширине переходной зоны Dp и значению коэффициента шума фильтра d (максимальным значениям осцилляций передаточной функции в единицах коэффициента передачи в полосе пропускания).
Кайзером установлено, что для заданного значения d произведение количества членов оператора НЦФ на ширину переходной зоны является величиной постоянной. Оно получило название D-фактора:
D = N·Dp/p.
С другой стороны, установлены следующие эмпирические соотношения между D-фактором и параметром b функции Кайзера:
D = (А-7.95)/14.36 при А>21.
= 0.9222 при А<21.
b = 0.1102(A-8.7) при А>50.
= 0 при А<21.
= 0.5842(A-21)0.4 + 0.07886(A-21), 21<А<50.
где: А = -20 log d - затухание в децибелах.
Приведенные выражения позволяют по заданному значению коэффициента шума d определить параметр b функции Кайзера, а через D-фактор число членов фильтра:
N = pD/Dp.
При проектировании полосовых фильтров проверка передаточной функции полученного оператора НЦФ исходному заданию по значению коэффициента шума является обязательной. Это объясняется тем, что поскольку полоса пропускания полосового фильтра ограничена двумя скачками, на передаточной характеристике возникают два центра осцилляций, при этом наложение осцилляций может как уменьшить, так и увеличить амплитуду суммарных осцилляций. Если за счет наложения произойдет увеличение амплитуды осцилляций, то расчет НЦФ следует повторить с уменьшением исходного значения d.
Пример расчета полосового фильтра.
Произвести расчет ПФ при следующих исходных параметрах: wн = 0.3p, wв = 0.6p, Dp = 0.1p, d = 0.02.
1. А = -20 log d. А = 34.
2. N = p (A-7.95)/(14.36 Dp). N = 18.
3. b = 0.5842(A-21)0.4 +0.07886(A-21). b = 2.62.
4. hо = (wв-wн)/p. hо = 0.3
5. h(n) = (sin nwв-sin nwн)/(np). h(n)= 0.04521, -0.24490, -0.09515, ... , 0.02721.
6. pn= Jo{b
} / Jo{b}. pn = 1.00, 0.997, 0.9882, .......7. Оператор фильтра: hn = pn h(n), n = 0, 1, 2,..., N. h-n = hn. hn = 0.3000, 0.04508, -0.2420, ....
8. Проверка по формуле: H(w) =
hn cos nw, 0 £ w £ p.Для оценки формы передаточной функции количество точек спектра в интервале 0-p достаточно задать равным 2N, т.е. с шагом Dw £ p/36.
Влияние конечной разрядности на цифровые фильтры должно быть минимальным и не создавать на их частотных характеристиках дополнительных неравномерностей и отклонения от заданной формы. С чисто практической точки зрения ограничение разрядности коэффициентов фильтра в целях повышения производительности вычислений лучше всего (и проще всего) выполнять непосредственно сравнением частотных характеристик с изменением разрядности от большей к меньшей. Следует учитывать, что ограничение разрядности может по разному сказываться на неравномерности фильтра в полосе пропускания и степени затухания сигналов в полосе подавления.
Ошибки отклонения e(w) частотной характеристики относительно заданной при проектировании кроме разрядности коэффициентов В в битах зависит также от размеров N оператора фильтра и в первом приближении может оцениваться по формулам: