Помогите c системой дифуров

thunder87l
Сообщений: 2
Зарегистрирован: 01 янв 2010, 21:00

Помогите c системой дифуров

Сообщение thunder87l » 02 янв 2010, 14:03

Есть задание - решить систему методом Эйлера.

Система задана следующим образом:

A=
4 -1 -1 0
-1 4 0 -2
-1 0 4 -1
0 -1 -1 4

начальное условие X(t=0)=(1, 0, 0, 0)

интервал 0<t<5уравнение вида dX/dt + AX = 01) Самое главное - объясните как это выглядит если записать систему в нормальном виде, по-человечески? 80% на то что дальше я как-нибудь соображу как писать код.У меня вот закралась мысль, что это что-то вида:a11*x1 + a12*x2 + a13*x3 + a14*x4 = C1a21*x1 + a22*x2 + a23*x3 + a24*x4 = C2a31*x1 + a32*x2 + a33*x3 + a34*x4 = C3a41*x1 + a42*x2 + a43*x3 + a44*x4 = C4где C1..4 - условие, хотя наверняка я ошибаюсь.2) Если не затруднит - объясните что делать дальше c системой (ибо ещё остаётся 20% на то, что до меня не дойдёт) B смысле - метод Эйлера это же тупо итерация вида:while (условие){new_X=old_X+step_X; //step_X - приращение по оси Хnew_Y=nextE(old_Y,old_X,new_X);old_X=new_X; old_Y=new_Y;}, где nextE(old_Y,old_X,new_X) это функция возвращающая new_Y=old_Y+(new_X-old_X)*F(old_X,old_Y)откуда добыть результат F(X,Y)? и что эта функция есть такое? строка из решаемой системы?3) Если кто-то может кинуть ссылку на реализацию кода на паскале/бейсике/си/си++ - пожалуйста киньте. Лучше конечно без этого, ибо всё одно необходимо будет пояснить суть программы, но если кинуть линк проще чем объяснить - тоже помощь. B крайнем случае - в коде покопаюсь.Спааaасибо всем огромное!
M Для написания формул используйте $$\LaTeX$$, разбором каракуль мы не занимаемся.
A Для написания формул используйте $$\LaTeX$$, разбором каракуль мы не занимаемся.
Последний раз редактировалось thunder87l 29 ноя 2019, 20:59, всего редактировалось 1 раз.
Причина: test

a_l_e_x86
Сообщений: 985
Зарегистрирован: 02 мар 2007, 21:00

Помогите c системой дифуров

Сообщение a_l_e_x86 » 02 янв 2010, 19:06

Если уравнение именно такого вида
$$dX/dt + AX = 0$$
(обычно подразумевается
$$dX/dt = AX $$
)
то система дифуров выглядит так

$$\{ x_1^\prime = -4x_1+x_2+x_3 \\   x_2^\prime = x_1-4x_2+2x_4 \\   x_3^\prime = x_1-4x_3+x_4 \\   x_4^\prime = x_2+x_3-4x_4 \\$$


Если не изменяет память итерации в методе Эйлера проводятся по формулам

$$x^{<n+1>}=x^{<n>}+hAx^{<n>}$$
h - шаг (число), вычисляется в зависимости от требуемой точности либо по формулам для метода Эйлера, либо по принципу Рунге (проводятся вычисления c шагами h h/2 и сравнивается усредненная разница между значениями в совпадающих точках)
Количество итераций равно длине интервала (5) деленное на шаг
в нашем случае
$$X^{<0>}=\(1\\0\\0\\0\)$$
Последний раз редактировалось a_l_e_x86 29 ноя 2019, 20:59, всего редактировалось 1 раз.
Причина: test

thunder87l
Сообщений: 2
Зарегистрирован: 01 янв 2010, 21:00

Помогите c системой дифуров

Сообщение thunder87l » 02 янв 2010, 21:25

OOO! Гранд мерси товарисч!

Я подозревал, что система открывается подобным образом, просто как-то не очень понял куда начальное условие засунуть

Bсё-таки трудно достаточно использовать численные методы, когда системы дифуров вырезали из курсa c понтом: "да вам оно всё равно в рот не упёрлось" И ещё раз спасибо за скорый ответ - eсли до послезавтра реализую в коде, автомат по численным считай - в кармане
Последний раз редактировалось thunder87l 29 ноя 2019, 20:59, всего редактировалось 1 раз.
Причина: test


Вернуться в «Математический анализ»

Кто сейчас на форуме

Количество пользователей, которые сейчас просматривают этот форум: нет зарегистрированных пользователей и 7 гостей