Круто! :D Тогда очевидным решением для преобразования ФНЧ -> ФНЧ будет преобразование к другой заданной частоте заграждения, что я и попробовал, в этом случае АЧХ удовлетворяет заданному коридору, но коэффициенты немного отличаются от тех что выдает программа filter, наверное предложенный вариант не...

Тут с инверсным Чебышевым что-то не так, в вашей статье приводится пример расчета фильтра с частотой среза 0.7 рад/cек, причем не производится преобразования к заданной частоте среза. Я рассчитал нормированный ФНЧ инверсного Чебышева исходя из примера в статье и посмотрел в матлабе характеристику с ...

Логично, сижу проверяю, но в то же время не уверен, что не требуются какие-нибудь дополнительные преобразования частоты для инверсного фильтра Впрочем я только рад что принципиально ничего неправильного не делаю, требуется только отладить имеющуюся методику расчета.

Добрался до Чебышева 1го рода, для РФ фильтра применил следующую формулу для частоты ФНЧ и все заработало: \frac {1}{ \cosh \left( \frac{1}{N} \mathrm{acosh} \left(\frac{ \varepsilon_s}{\varepsilon_p} \right) \right) } ; Сейчас делаю Чебышева 2го рода (инверсный), проверил расчет нормированной перед...

Спасибо большое! У меня получились такие же коэффициенты как и у вас, но с небольшим НО. Дело в том что все заработало когда я применил формулу \Omega_0 = \sqrt[N]{\frac{\varepsilon_p}{\varepsilon_s}} а не \Omega_0 = \sqrt[N]{\frac{\varepsilon_s}{\varepsilon_p}} Я так думаю что для Режекторных фильт...

Моя программная реализация выдает неверные коэффициенты для Режекторных фильтров. Для примера решил рассчитать простенький цифровой фильтр 1го порядка Баттерворта со следующими характеристиками: \omega_{pl} = 0.2\pi рад/сек \omega_{sl} = 0.3\pi рад/сек \omega_{sh} = 0.6\pi рад/сек \omega_{ph} = 0.8\...

Да согласен, затупил, с заменой переменных все стало гораздо проще. Кстати, можно реализовать вычисление производной частотной характеристики по аналогии как функции freqz, и использовать полученные мнимые I'(w) и действительные части R'(w) в формуле выше, тогда все будет ещё точнее

Хочется отобразить график групповой задержки. Как правильно это сделать? По определению групповая задержка это минус производная от ФЧХ по частоте. Как самый простой вариант решения задачи это вычислить N точек ФЧХ. Сделать unwrap для фазы, и взять от полученного массива производную скажем по 3ем то...

Спасибо за статью Пример расчета цифрового фильтра нижних частот Баттерворта. Попробовал сверить конечные коэффициенты с тем что выдает программа filter, и получил другой порядок коэффициентов b0 и b3. Возможно это ошибка при подстановке чисел в формулу.

Да, вы правы, мои формулы для диапазона [0...\pi] (что вообще говоря неверно). А для устойчивого фильтра необходимо чтобы его действительная часть была меньше 0, т.е. чтобы коэффициенты затухания в s-плоскости были отрицательные. Поэтому нужен диапазон углов полюсов [-\pi/2 ... \pi/2] . Для этого ди...