Продолжение примера

Продолжение примера.

8. Вычисляем значения коэффициентов G_m, b_m и c_m:

- N=4: g = 1.637, G₁ = 0.203, G₂ = 0.149, b₁ = 0.681, b₂ = 0.501, c₁ = 0.492, c₂ = 0.098.

- N=5: g = 1.698, G₁ = 0.203, G₂ = 0.151, b₁ = 0.763, b₂ = 0.568, c₁ = 0.574, c₂ = 0.171.

9. Подставляем вычисленные коэффициенты в выражения (10.1.20, 10.1.24) и вычисляем значения передаточных функций при z = exp(-jwDt). Графики полученных функций приведены на рис. 10.1.3. На рис. 10.1.4

Во временной области фильтрация выполняется последовательной сверткой входного сигнала с операторами ячеек фильтра:

y_k = x_k ③ {h₀(i)} ③ h₁(i) ③ … ③ h_М(i), i = 0,1,2.

Уравнение рекурсивной фильтрации для m-го оператора фильтра:

y_k = G_m (x_k+2x_k-1+x_k-2) + b_m y_k-1 - c_m y_k-2. (10.1.25)

Уравнение рекурсивной фильтрации для дополнительного h₀(i) линейного оператора фильтра при нечетном N:

y₀ = (x_k+x_k-1)/(g+1) + y_k-1·(g-1)/(g+1) (10.1.26)

Рис. 10.1.4.

10. Каждый оператор фильтра имеет определенную передаточную функцию, что можно видеть на рис. 10.1.4. Порядок последовательной свертки сигнала с операторами фильтра значения не имеет, но с учетом разрядности ячеек памяти звено H₁(f) целесообразно реализовать за H₂(f).

11. Для оценки длительности импульсной реакции фильтра подаем на вход фильтра импульс Кронекера на отсчете k = 3, и начинаем фильтрацию со второго отсчета (что обеспечивает начальные условия фильтрации на точках k=0 и k=1). Коэффициент усиления дисперсии шумов (сумма квадратов значений импульсного отклика) равен 0.341 при N=5, и 0.278 при N=4.

Дата добавления: 2014-01-03