Московский физико-технический институт
Опубликован: 25.10.2007 | Доступ: свободный | Студентов: 3791 / 1086 | Оценка: 4.50 / 4.33 | Длительность: 24:00:00
ISBN: 978-5-9556-0065-9
Специальности: Программист, Математик

Лекция 10: Численные методы решения жестких систем обыкновенных дифференциальных уравнений

9.2. Сингулярно - возмущенные задачи

Рассмотрим простейшую нелинейную жесткую систему А.Н.Тихонова [9.5] (сингулярно - возмущенная задача Коши с малым параметром при производной):

\begin{gather*}
\varepsilon \dot {u} = f(u, v), \\  
\dot {v} = g(u, v), \\  
0 < \varepsilon  \ll 1, 
\end{gather*}

или

\begin{gather*}
\dot {u} = \frac{1}{\varepsilon } f(u, v), \\  
\dot {v} = g(u, v), 
\end{gather*}

с начальными условиями

u(0) = u0, v(0) = v0.

Из характеристического уравнения находим:

\begin{gather*}
\det  \left| \begin{array}{cc}
{\frac{1}{\varepsilon } f^{\prime}_u - \lambda } & {\frac{1}{\varepsilon }f^{\prime}_v}\\ 
{g^{\prime}_u} & {g^{\prime}_v - \lambda }  \\
\end{array} \right| = \lambda ^2 - \lambda  \left({\frac{1}{\varepsilon } f^{\prime}_u + g^{\prime}_u}\right) + \frac{1}{\varepsilon } (f^{\prime}_u g^{\prime}_v - g^{\prime}_u f^{\prime}_v), \\  
\lambda_1 = \frac{1}{\varepsilon } f^{\prime}_u + O(1), \quad \lambda_2 = O(1), 
\end{gather*}

т.е. \lambda _{1} — жесткая, а \lambda _{2} — мягкая часть спектра.

Первое, что хотелось бы сделать, — упростить задачу, положив \varepsilon  \approx  0. В этом случае система приобретает вид:

\begin{gather*}
f(u, v) = 0, \\  
\dot {v} = g(u, v), \\  
u(0) = u_0,  v(0) = v_0
\end{gather*}

и описывает траекторию, проходящую по кривой f(u, v) = 0. Эта кривая играет существенную роль в решении исходной системы ОДУ. Однако полностью поведение решения упрощенная (невозмущенная) система описывать не будет, поскольку начальные данные не обязательно должны лежать на этой кривой.

Для того чтобы понять эту ситуацию, рассмотрим фазовый портрет конкретной системы ОДУ вида:

\begin{gather*}
\varepsilon \dot {u} = v - \frac{1}{3} u^3 + u, \\  
\dot {v} = - u, \\  
u(0) = u_0,  v(0) = v_0 .
\end{gather*}

Эта система эквивалентна нелинейному уравнению второго порядка — уравнению Ван - дер - Поля. Приведенный выше вид иногда называется представлением Льенара. Подробнее о свойствах рассматриваемой задачи можно прочитать в [9.6], [9.7], [9.8], [9.9].


Рис. 9.4.

Кривая f(u, v) = 0 для этого случая изображена на рис. 9.4. Она делит плоскость {u, v} на две части: f > 0 и f < 0 ; вдали от кривой, поле скоростей {du, dv} направлено почти горизонтально влево (вправо) в зависимости от знака f. На самой кривой выделяются две устойчивые ветви AB и CD, где f'u < 0, и неустойчивая ветвь BD, на которой f'u > 0. Опишем теперь качественно поведение траектории рассматриваемой системы ОДУ, состоящей из следующих участков.