>
Одержуємо:
Wn = 0.1000
b = 0.7294 -2.1883 2.1883 -0.7294
a = 1.0000 -2.3741 1.9294 -0.5321
Для порівняння отриманих коефіцієнтів
* У зарубіжній літературі і в MATLAB коефіцієнти відносяться В«aВ» відносяться до отсчетам вихідного сигналу, а В«bВ» до отсчетам вхідного. У вітчизняній літературі ж і в цій роботі навпаки. br/>
Таблиця 4 - Порівняння коефіцієнтів розрахованих вручну і коефіцієнтів отриманих в системі MatLab.
КоэффициентыРассчитанныеMatLab 0,7310,7294 -2,193-2,1883 2,1942,1883 -0,731-0,7294 11,0000 -2,380-2,3741 1,9521,9294 -0,545-0,5321
З таблиці видно, що результати обчислень збігаються з невеликою відносною похибкою.
6. Частотні характеристики цифрового фільтра
Для того щоб побудувати АЧХ і ФЧХ синтезованого фільтра знову скористаємося середовищем MATLAB. Для цього до попередньої програмі додамо наступні рядки:
n = 3; fc = 500; fs = 10000;
Wn = 2 * (fc/fs)
[b, a] = butter (n, Wn, 'high') (b, a, 500, fd)
Т.ч. отримаємо:
В
Рисунок 4 - АЧХ і ФЧХ синтезованого цифрового фільтра Баттерворта верхніх частот третього порядку.
. Імпульсна характеристика цифрового фільтра
Для того щоб побудувати імпульсну характеристику синтезованого фільтра (для 50 перших відліків) використовуємо MATLAB. Для цього до попередньої програмі додамо наступні рядки:
n = 3; fc = 500; fs = 10000;
Wn = 2 * (fc/fs)
[b, a] = butter (n, Wn, 'high')
impz (b, a, 50)
Т.ч. отримаємо:
В
Рисунок 5 - Імпульсна характеристика цифрового синтезованого фільтра Баттерворта верхніх частот третього порядку (50 відліків)
Щоб побудувати БІХ, візьмемо 350 перших відліків. Отримаємо:
n = 3; fc = 500; fs = 10000;
Wn = 2 * (fc/fs)
[b, a] = butter (n, Wn, 'high')
impz (b, a, 350)
Т.ч. отримаємо:
В
Малюнок 6 - Імпульсна характеристика цифрового синтезованого
фільтра Баттерворта верхніх частот третього порядку (350 відліків).
8. Карта нулів і п...
