Разработка, производство и продажа радиоэлектронной аппаратуры
|
Карта сайта
|
Пишите нам
|
В избранное
Требуется программист в Зеленограде
- обработка данных с датчиков; ColdFire; 40 тыс.
e-mail:
jobsmp@pochta.ru
Телесистемы
|
Электроника
|
Конференция «Микроконтроллеры и их применение»
Нужен быстрый (!!) алгоритм нахождения комплексных корней полинома вида (+):
Отправлено
argus98
(81.22.205.230)
16 июля 2012, г. 10:05
В ответ на:
Значит, вопросы должны усложниться
отправлено
Крок
16 июля 2012, г. 00:30
z^N + a[N-2]*z^(N-2) + a[N-3]*z^(N-3) + a[N-4]*z^(N-4) + ... + a[1]*z + a[0]
N - большие (10, 20, ..., 100500)
a[n] - действительные
|a[n]| всегда <1
действительных корней нет
кратных корней тоже нет
Ньютона, Монте-Карло, подбор - не предлагать
Составить ответ
|
Вернуться на конференцию
Ответы
У Марпла в главе про метод Прони приводится программа и ссылки на алгоритм факторизации полинома (с комплексными коэффициентами) не знаю насколько быстрый, но, наверное, не самый худший.
-
GF(6)
(
пустое
, 18.07.2012, 15:41:45
193.33.145.67
)
дпф возьми от последовательности an...a0. Потом минимумы от его модуля. Они тебе некое приближение к мнимой части корней дадут:-)
-
-=ВН=-
(
пустое
, 17.07.2012, 15:36:53
178.249.67.15
)
дпф возьми от последовательности an...a0. Потом минимумы от его модуля. Они тебе некое приближение к мнимой части корней дадут:-)
-
-=ВН=-
(
пустое
, 17.07.2012, 15:36:52
178.249.67.15
)
А что будет если модули несколько корней будут ближе чем 2/N ?
-
GF(6)
(
пустое
, 18.07.2012, 15:23:50
193.33.145.67
)
А с каким окном брать?
-
GF(6)
(
пустое
, 18.07.2012, 15:15:6
193.33.145.67
)
попробуем.. Я правда больше склоняюсь к разложению a[n] по полиномам Лагранжа.... или зря потрачу время?
-
argus98
(
пустое
, 17.07.2012, 16:48:25
81.22.205.230
)
...Мюллера
-
GF(6)
(
пустое
, 16.07.2012, 16:03:59
193.33.145.67
)
Старик Кардано спёкся после N=3 :-) Кстати, если Z=2, не боитесь неожиданных результатов?
-
Крок
(54 байт, 16.07.2012, 12:05:59
62.118.59.14
)
а чего ее бояться? нули ких фильтра имеют полное право быть снаружи единичной окружности, в том числе нуль могет быть и в z=2 :-)
-
-=ВН=-
(
пустое
, 17.07.2012, 15:40:27
178.249.67.15
)
Там пристутсвует сотая степень. А я с детства боюсь чисел, больших, чем 2^16
-
Крок
(
пустое
, 17.07.2012, 19:44:38
188.123.248.109
)
зы Z не может =2 по определению
-
argus98
(
пустое
, 16.07.2012, 20:05:30
87.117.189.102
)
Даже 1.5 - много. Или Зы всегда меньше 1?
-
Крок
(
пустое
, 16.07.2012, 20:23:38
188.123.248.109
)
Ньютона я знаю. Хочется знать то, чего не знаю
-
argus98
(
пустое
, 16.07.2012, 20:03:35
87.117.189.102
)
Отправка ответа
Имя*:
Пароль:
E-mail:
Тема*:
Сообщение:
Ссылка на URL:
URL изображения:
если вы незарегистрированный на форуме пользователь, то
для успешного добавления сообщения заполните поле, как указано ниже:
введите число 85:
Перейти к списку ответов
|
Конференция
|
Раздел "Электроника"
|
Главная страница
|
Карта сайта
Web
telesys.ru