Эллиптический фильтр

[_usa_))]

День добрый!

Этот вопрос скорее к автору статьи, но все же ...
Вот статья http://www.dsplib.ru/content/filters/ch6/ch6.html все расписано и разжевано, пытаюсь
повторить расчет (кодю прям по алгоритму) все замечательно но спотыкаюсь об (24), что есть (11) и (16)
получаемое значение упорно другое (-0.01828959) что не верно конечно ...

В связи с этим пока что два вопроса:
1. Действительно ли k1 не пересчитывается (вместо 0.0028615 -> 0.999996)
(у меня при значении k1=0.99999786586 дает таки результат в 0.18181
правда при k1=0.999996 дает уже 0.02396)
2. Нет ли опечатки в (24), (11), (16) ?


Пасиба

[Бахурин Сергей]

В ближайшее время сделаю матлаб файл для расчета эллиптического фильтра для проверки

[_usa_))]

Гы-гы ... Качнул таки Матлаб, а он оказывается не умеет считать э.ф. Якоби для комплексных аргументов ...
Если пользоваться приведением sn(u*j,k)=-j*(sn(u,k')/cn(u,k')) то опять таки получается мимо ...
Что-ж бум думать и/или ждать Ваше решение в матлабе

[Бахурин Сергей]

В приложении пример матлаб кода для расчета эллиптического фильтра согласно статье.

ellip_example.zip

(5.93 КБ) 710 скачиваний

[_usa_))]

День добрый! Спасибо за программу !

Наконец разобрался (с определенной степенью) ... Теперь по порядку:

1. sn^-1 для меня все же соответствует 1/sn, оказывается это все-таки инверсный sn по первому аргументу u.
не знаю, может так принято обозначать инверсные функции, но для меня это было загадкой

2. расчет всех выражений включающих эллиптические функции Якоби в статье приводится следуя алгоритму
http://www.dsplib.ru/content/ellipfunc/ellipfunc.html формулы (19) и (17), там вроде все тоже ОК [я правда
пользовался выкладками из других источников (справочник по эллиптическим функциям и т.д.)]
т.е. в качестве первого аргумента используется нормированный на K(k) параметр: вычисляемые
sne(u,k) и asne(u,k) в предложенном файле в действительности (если говорить о самих эллиптических функциях)
вычисляются для sn(u*K(k),k) и sne(u*K(k),k)

3. и самое интересное ... ПО МОЕМУ ВСЕ-ТАКИ в (24), (16) и (11) допущена опечатка, -
есть: v0=(-j/n)*sn(j*K(k1)/ep,k1)
а должно быть: v0=(-j/n)*sn(j/ep,k1)/K1
что соответствует формуле (65) в http://www.ece.rutgers.edu/~orfanidi/ece521/notes.pdf

вот вроде все ... теперь мои вычисления совпадают с предложенной программой на MatLab'е ...
этот путь наверное все же менее затратен чем самостоятельный вывод выражения (24) и следующих из него

[Бахурин Сергей]

sne(u,k) и asne(u,k) в предложенном файле в действительности (если говорить о самих эллиптических функциях)
вычисляются для sn(u*K(k),k) и sne(u*K(k),k)

да так и есть. опечатки в статье тоже есть. поправлю.

[_usa_))]

День добрый!

Еще одна незначительная отпечатка
В формуле (25) там где выписано значение сигма1 надо 0.05674 а не 0.005674
В дальнейших расчетах все верно, - используется 0.05674

[_usa_))]

Тьфу ты, пардон, это в (28) написано 0.005674 а надо 0.05674