День добрый!
Этот вопрос скорее к автору статьи, но все же ...
Вот статья 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) ?
Пасиба
Эллиптический фильтр
- Бахурин Сергей
- Администратор
- Сообщения: 1116
- Зарегистрирован: 05 окт 2010, 19:55
- Контактная информация:
Re: Эллиптический фильтр
В ближайшее время сделаю матлаб файл для расчета эллиптического фильтра для проверки
Re: Эллиптический фильтр
Гы-гы ... Качнул таки Матлаб, а он оказывается не умеет считать э.ф. Якоби для комплексных аргументов ...
Если пользоваться приведением sn(u*j,k)=-j*(sn(u,k')/cn(u,k')) то опять таки получается мимо ...
Что-ж бум думать и/или ждать Ваше решение в матлабе
Если пользоваться приведением sn(u*j,k)=-j*(sn(u,k')/cn(u,k')) то опять таки получается мимо ...
Что-ж бум думать и/или ждать Ваше решение в матлабе
- Бахурин Сергей
- Администратор
- Сообщения: 1116
- Зарегистрирован: 05 окт 2010, 19:55
- Контактная информация:
Re: Эллиптический фильтр
В приложении пример матлаб кода для расчета эллиптического фильтра согласно статье.
Re: Эллиптический фильтр
День добрый! Спасибо за программу !
Наконец разобрался (с определенной степенью) ... Теперь по порядку:
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) и следующих из него
Наконец разобрался (с определенной степенью) ... Теперь по порядку:
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) и следующих из него
- Бахурин Сергей
- Администратор
- Сообщения: 1116
- Зарегистрирован: 05 окт 2010, 19:55
- Контактная информация:
Re: Эллиптический фильтр
да так и есть. опечатки в статье тоже есть. поправлю.sne(u,k) и asne(u,k) в предложенном файле в действительности (если говорить о самих эллиптических функциях)
вычисляются для sn(u*K(k),k) и sne(u*K(k),k)
Re: Эллиптический фильтр
День добрый!
Еще одна незначительная отпечатка
В формуле (25) там где выписано значение сигма1 надо 0.05674 а не 0.005674
В дальнейших расчетах все верно, - используется 0.05674
Еще одна незначительная отпечатка
В формуле (25) там где выписано значение сигма1 надо 0.05674 а не 0.005674
В дальнейших расчетах все верно, - используется 0.05674
Re: Эллиптический фильтр
Тьфу ты, пардон, это в (28) написано 0.005674 а надо 0.05674