Купить Matlab  |  Mathematica  |  Mathcad  |  Maple  |  Statistica  |  Другие пакеты Поиск по сайту
Internet-класс  |  Примеры  |  Методики  |  Форум  |  Download
https://hub.exponenta.ru/


Курс ОДУ.
Готовые занятия
 
Занятие 7
Теоретический материал Теоретическая справка Примеры Задачи для самостоятельного решения Контрольные вопросы

Решение системы ~ Задача Коши ~ Интегральная кривая ~ Теорема существования и единственности ~ Численное решение задачи Коши

 

Рассмотрим систему дифференциальных уравнений
y'1 = f1(x, y1, y2 , ..., yn),
y'2 = f2 (x, y1, y2 , ..., yn),
..............................
y'n = fn (x, y1, y2 , ..., yn),
где x — независимая переменная, а
y1(x), y2(x), ..., yn(x) — неизвестные функции, n — порядок системы.

Обозначив
t_1.gif (1415 bytes)  t_2.gif (1800 bytes)
запишем систему в векторной форме
`Y '=`F(x,`Y ).

Решением системы называется вектор-функция `Y , которая определена и нерерывно дифференцируема на интервале (a, b) и удовлетворяет системе, т.е. для всех x0О (a, b) справедливо
`Y '(x) =`F(x,`Y (x)).

Задачей Коши (задачей с начальными условиями) называется следующая задача: найти такое решение`Y (x) системы `Y '=`F(x,`Y ), что `Y (x0) =`Y 0,
где x0 — заданное число, а `Y 0 — заданный вектор.

Интегральной кривой системы называется кривая в (n+1) -мерном пространстве Rn+1x,y,  заданная уравнением `Y =`Y (x), где `Y (x) - решение системы.Таким образом, решить задачу Коши — это значит найти интегральную кривую, проходящую через заданную точку пространства Rn+1x,y.

Для нормальных систем обыкновенных дифференциальных уравнений справедлива следующая теорема существования и единственности решения задачи Коши.

Если вектор-функция `F(x,`Y (x)) и ее частные производные по переменным yi , i = 1, 2, ..., n, непрерывны в области G пространства Rn+1x,y, то на некотором интервале (x0 -h, x0+h) существует единственное решение системы
`Y '(x) =`F(x,`Y (x)) ,
удовлетворяющее начальному условию
`Y (x0) =`Y 0,
т.е. через каждую точку области G проходит единственная интегральная кривая системы.

Подробнее геометрическая интерпретация систем обыкновенных дифференциальных уравнений и их решений рассмотрена в разделе, посвященном изучению автономных систем.

 

ПРИМЕР 1. Интегральная кривая для системы дифференциальных уравнений 2-го порядка.

 

Если задачу об отыскании всех решений системы дифференциальных уравнений удается свести к конечному числу алгебраических операций, операций интегрирования и дифференцирования известных функций, то говорят, что система интегрируется в квадратурах.
В приложениях крайне редко встречаются системы, интегрируемые в квадратурах. Поэтому для исследования дифференциальных уравнений широко используются приближенные, численные методы их решения.

Численное решение задачи Коши для системы обыкновенных дифференциальных уравнений состоит в построении таблицы приближенных значений
yi,1, yi,2 , ..., yi,N   компонент yi(xj,) вектора решения `Y (x) в точках
x1, x2 , ..., xN.

Чтобы получить расчетные формулы метода Рунге—Кутты для систем дифференциальных уравнений, достаточно в расчетных формулах для уравнений первого порядка заменить
y, f(x, y), k1, k2, k3 , k4 на
`Y , F(x,`Y ), `k 1, `k 2, `k 3, `k 4.

 

ПРИМЕР 2. Численное решение задачи Коши для системы дифференциальных уравнений 2-го порядка методом Рунге-Кутты. Графики чтсленного решения.

 

Задача Коши для любого дифференциального уравнения n -го порядка 
y(n) = f(x, y' , y' , ..., y(n-1)),
y(x0 )= y0 , y'(x0 )= y0,1, y''(x0 )= y0,2 , ..., y(n-1)(x0 )= y0,n-1
легко сводится к задаче  для системы дифференциальных уравнений порядка
`Y '=`F(x,`Y ), где
t_3.gif (1844 bytes) t_4.gif (1579 bytes)t_5.gif (1330 bytes)

 

ПРИМЕР 3. Сведение задачи Коши для дифференциального уравнения 2-го порядка к задаче Коши для системы  дифференциальных уравнений 2-го порядка.

В начало страницы

Примеры Задачи для самостоятельного решения Контрольные вопросы
| На первую страницу | Поиск | Купить Matlab

Исправляем ошибки: Нашли опечатку? Выделите ее мышкой и нажмите Ctrl+Enter


Copyright © 1993-2024. Компания Softline. Все права защищены.

Дата последнего обновления информации на сайте: 04.03.17
Сайт начал работу 01.09.00