ФНЧ Баттерворта

Все что касается фильтрации
nikolausi
Сообщения: 5
Зарегистрирован: 22 дек 2016, 21:04

ФНЧ Баттерворта

Сообщение nikolausi »

приветствую
Пытаюсь рассчитать ФНЧ Баттерворта по этой примере
http://www.dsplib.ru/content/filters/bu ... terex.html
Возникло проблемма с формулой (12). не могу понять, как из (11) получить (12) и откуда взялос а0=32,82? Кто нибуд может помочь?
Заранее благодарю

Аватара пользователя
Бахурин Сергей
Администратор
Сообщения: 1114
Зарегистрирован: 05 окт 2010, 19:55
Контактная информация:

Re: ФНЧ Баттерворта

Сообщение Бахурин Сергей »

(12) из 11 получается путем раскрытия скобок и приведения подобных относительно $z^{-1}$. После этого получится коэффициент a0 в знаменателе а0=32,82 = 14.832+12.074+4.914+1, который выносится за скобки учитывается в числителе

nikolausi
Сообщения: 5
Зарегистрирован: 22 дек 2016, 21:04

Re: ФНЧ Баттерворта

Сообщение nikolausi »

Бахурин Сергей
Спасибо за ответ
Вы не подскажете, где нибудь можно найти программу в матлабе для этой примери?

Аватара пользователя
Бахурин Сергей
Администратор
Сообщения: 1114
Зарегистрирован: 05 окт 2010, 19:55
Контактная информация:

Re: ФНЧ Баттерворта

Сообщение Бахурин Сергей »

А функция butter чем не устраивает?

nikolausi
Сообщения: 5
Зарегистрирован: 22 дек 2016, 21:04

Re: ФНЧ Баттерворта

Сообщение nikolausi »

Наверное устраивает, только я не смог использовать функцию butter для параметров
Rp=1;
Rs=30;
Ws=0.6*pi;
Wp=0.2*pi;

чтобы принять коэффициенты фильтра b0, b1, b2, b3, a0, a1 a2, a3

Аватара пользователя
Бахурин Сергей
Администратор
Сообщения: 1114
Зарегистрирован: 05 окт 2010, 19:55
Контактная информация:

Re: ФНЧ Баттерворта

Сообщение Бахурин Сергей »

Butter не поддерживает параметр Rp Rs просто задается порядок и wp

Аватара пользователя
Бахурин Сергей
Администратор
Сообщения: 1114
Зарегистрирован: 05 окт 2010, 19:55
Контактная информация:

Re: ФНЧ Баттерворта

Сообщение Бахурин Сергей »

Rp при этом -3 dB

nikolausi
Сообщения: 5
Зарегистрирован: 22 дек 2016, 21:04

Re: ФНЧ Баттерворта

Сообщение nikolausi »

пробовал так:

Код: Выделить всё

wp=0.2*pi;
n=3;
[b,a] = butter(n,wp)

получил:
b =
0.2871 0.8614 0.8614 0.2871
a =
1.0000 0.7418 0.4798 0.0756

а в примере:
b=0.03047, 0.0914, 0.0914, 0.03047
a=1, -1.483, 0.9296, -0.2033

не подскажете где проблема?

Аватара пользователя
Бахурин Сергей
Администратор
Сообщения: 1114
Зарегистрирован: 05 окт 2010, 19:55
Контактная информация:

Re: ФНЧ Баттерворта

Сообщение Бахурин Сергей »

Матлаб функция не поддерживает rp параметр.

nikolausi
Сообщения: 5
Зарегистрирован: 22 дек 2016, 21:04

Re: ФНЧ Баттерворта

Сообщение nikolausi »

rp параметр
а что это такое?
вы можете подсказать, есть ли в матлабе функция, или можно найти где нибудь код для матлаба, чтобы можно било получить тот результат что в примере http://www.dsplib.ru/content/filters/bu ... terex.html

Ответить