![]() |
|
![]() ![]() ![]() |
|
vagon1 |
|
|||
Новичок Профиль Группа: Участник Сообщений: 9 Регистрация: 18.1.2006 Репутация: нет Всего: нет |
Подскажите пожалуйста, есть однородное дифференциальные уравнения типа
а(58)*y(58)^58+a(57)*y(57)^57+......+a(1)*y(1)^1+a(0)*y(0)=0 или а(16)*y(16)^16+a(16)*y(16)^16+......+a(1)*y(1)^1+a(0)*y(0)=0 можно ли его решить программированием и как коэффициенты при этом уравнении изменяются при определенных условиях Известно, что аналитического решения этих уравнений нет (в природе есть решение диф-го уравнения только пятого порядка) известно что подобные уравнения решались на машинах "Нари-С" ![]() |
|||
|
||||
Cr@$h |
|
||||
![]() Исследователь ![]() ![]() ![]() Профиль Группа: Участник Клуба Сообщений: 1693 Регистрация: 3.4.2005 Где: Санкт-Петербург, Россия Репутация: 4 Всего: 41 |
Насколько я знаю уже для полиномов пятого порядка нельзя найти решение аналитически...
![]() Переведи дифуру в пространтсво состояний, составив систему из 58 дифуров первого порядка, и решай полную проблему собственных значений для матрицы или просто как обычную систему дифур. |
||||
|
|||||
vagon1 |
|
|||
Новичок Профиль Группа: Участник Сообщений: 9 Регистрация: 18.1.2006 Репутация: нет Всего: нет |
то есть для определения корней этого диф-го уравнения а58*y^58+a57*y^57+.....+a1*y^1+a0*y=0
мне необходимо написать систему уравнений вида y1=y; y2=y' ; y3=y'' ..... |
|||
|
||||
Cr@$h |
|
||||
![]() Исследователь ![]() ![]() ![]() Профиль Группа: Участник Клуба Сообщений: 1693 Регистрация: 3.4.2005 Где: Санкт-Петербург, Россия Репутация: 4 Всего: 41 |
Да, если я правильно тебя понял. Должен в итоге получить систему x' = Ax 58 порядка. После замены
матрица получится примерно такого вида:
Дерзай! |
||||
|
|||||
vagon1 |
|
|||
Новичок Профиль Группа: Участник Сообщений: 9 Регистрация: 18.1.2006 Репутация: нет Всего: нет |
Cr@$h
То усть после написания этой матрицы без труда будет отыскать 58 корней уравнения |
|||
|
||||
SoWa |
|
|||
![]() Харекришна ![]() ![]() ![]() ![]() Профиль Группа: Комодератор Сообщений: 2422 Регистрация: 18.10.2004 Репутация: 6 Всего: 74 |
Потрудится придется. Но создание этой матрицы- уже больше чем половина проблеммы решена.
-------------------- Всем добра ![]() |
|||
|
||||
maxim1000 |
|
|||
![]() Эксперт ![]() ![]() ![]() ![]() Профиль Группа: Участник Сообщений: 3334 Регистрация: 11.1.2003 Где: Киев Репутация: 33 Всего: 110 |
в принципе задача сводится либо к поиску собственных чисел матрицы (если переходитьк многомерному уравнению первого порядка), либо к поиску корней полинома, если решать, как обычно...
и та, и та решаются численно... -------------------- qqq |
|||
|
||||
Cr@$h |
|
|||
![]() Исследователь ![]() ![]() ![]() Профиль Группа: Участник Клуба Сообщений: 1693 Регистрация: 3.4.2005 Где: Санкт-Петербург, Россия Репутация: 4 Всего: 41 |
Найдя собственные числа матрицы, можешь записать решение аналитически через экспоненты, ну это ты, наверное, и сам знаешь. Можешь решать задачу Коши численно, как угодно, т.е. собственно как систему дифур. Собственные числа через полиномы стараются не искать -- неустойчивые алгоритмя: при малых изменениях коэффициентов полинома имеем большие изменения его корней. |
|||
|
||||
mvdr |
|
|||
физик ![]() ![]() ![]() Профиль Группа: Участник Сообщений: 1349 Регистрация: 31.12.2004 Где: Волгоград, Россия Репутация: 1 Всего: 42 |
все возможно.
Решение - Вопрос времени -------------------- Появляюсь редко, но часто метко Изображать идиота сложнее, чем изображать умного: полезнее и не каждому дано |
|||
|
||||
220v |
|
|||
![]() Шустрый ![]() Профиль Группа: Участник Сообщений: 57 Регистрация: 15.1.2006 Где: г.Владимир Репутация: нет Всего: 1 |
Есть много способов, например метод Адамса, или Рунге-Кутта. Попробуй набрать в поиске - наверняка найдешь много интересного.
|
|||
|
||||
Cr@$h |
|
|||
![]() Исследователь ![]() ![]() ![]() Профиль Группа: Участник Клуба Сообщений: 1693 Регистрация: 3.4.2005 Где: Санкт-Петербург, Россия Репутация: 4 Всего: 41 |
Да с этим он разобрался. Просто такой метод даст одно решение, а вот, если получить собственные значения, то получим семейство решений. Ему это, вроде, надо было.
Кстати, рекомендую системный метод Ю.В. Ракитского. Задача свелась к x' = Ax. Ее точное решение x = e^A + С. Все. Можно разложить е в ряд и до приемлемой точности взять членов. Всех подробностей уже и не помню. |
|||
|
||||
vagon1 |
|
|||
Новичок Профиль Группа: Участник Сообщений: 9 Регистрация: 18.1.2006 Репутация: нет Всего: нет |
Cr@$h
Спасибо |
|||
|
||||
Alex1984 |
|
|||
![]() Опытный ![]() ![]() Профиль Группа: Участник Сообщений: 299 Регистрация: 6.3.2005 Где: Киев Репутация: нет Всего: 2 |
Уравнения такого порядка можно решить написав программку. Знакомый писал для 11 степени, просто программу немного дописать придеться. Вот и все.
Это не проблемма. Проблемма например найти функцию в 58 мерном пространстве. Тут нужно трудиться (вопрос времени, но писать много и есть вариант просто механическу допустить ошибку) |
|||
|
||||
esperant0 |
|
|||
![]() Опытный ![]() ![]() Профиль Группа: Участник Сообщений: 714 Регистрация: 20.5.2005 Репутация: 4 Всего: 14 |
Вам же сказали, что аналитически выше 5 порядка не решается -------------------- Student->Teacher Assistant ->Research assistant->Microsoft Software Development Engineer Пользователь получил наказание за то, что проигнорировал замечание которое было написано модератором а затем стерто и которое он - пользователь не мог видеть. |
|||
|
||||
Cr@$h |
|
|||
![]() Исследователь ![]() ![]() ![]() Профиль Группа: Участник Клуба Сообщений: 1693 Регистрация: 3.4.2005 Где: Санкт-Петербург, Россия Репутация: 4 Всего: 41 |
||||
|
||||
![]() ![]() ![]() |
Правила форума "Алгоритмы" | |
|
Форум "Алгоритмы" предназначен для обсуждения вопросов, связанных только с алгоритмами и структурами данных, без привязки к конкретному языку программирования и/или программному продукту.
Если Вам понравилась атмосфера форума, заходите к нам чаще! С уважением, maxim1000. |
1 Пользователей читают эту тему (1 Гостей и 0 Скрытых Пользователей) | |
0 Пользователей: | |
« Предыдущая тема | Алгоритмы | Следующая тема » |
|
По вопросам размещения рекламы пишите на vladimir(sobaka)vingrad.ru
Отказ от ответственности Powered by Invision Power Board(R) 1.3 © 2003 IPS, Inc. |