§ 1.4. Замены переменных

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

1.4.1. Линейное неоднородное уравнение. Пусть в уравнении

xў = a(t)x + b(t)

(1)

функции a, b: I ® R непрерывны.

(2)

В § 1.2 для соответствующего однородного уравнения

xў = a(t)x

(3)

найдено общее решение

x = Ft0(t)C,

(4)

где

Ft0(t) = expт t t0 a(s) ds.

(5)

Теперь мы будем искать решение линейного неоднородного уравнения (1) в том же виде (4), только считая C не константой, а функцией от t. Такой способ решения называется методом вариации произвольной постоянной. Формулу (4) можно рассматривать как замену старой неизвестной функции x на новую C. Подставим (4) в (1):

Fўt0 + Ft0Cў = a(t)Ft0C + b(t).

(6)

Функция Ft0 есть одно из решений уравнения (3), поэтому

Fўt0 = a(t)Ft0.

Отсюда и из (6) получаем:

Cў = F –1 t0 (t)b(t)

и, следовательно,

C = т t t0 F –1 t0 (s)b(s)ds + C1.

Вместе с (4) это дает:

x = Ft0(t)C + Ft0(t)т t t0 F –1 t0 (s)b(s)ds

(7)

(мы переобозначили произвольную постоянную).

Логическая схема проведенного рассуждения такова: из данного уравнения (1) и уравнения замены переменной (4) мы вывели как следствие уравнение (7): (1) Щ (4) Ю (7). Мы будем говорить в таком случае, что (7) следует из (1) с учетом замены (4), и писать (1) Ю(4) (7). Далее, нетрудно видеть, что и наоборот, (7) Ю(4) (1), т. е. (7) Щ (4) Ю (1). Итак, (1) Ы(4) (7), причем уравнения (1) и (7) не содержат новой неизвестной функции C. Нас интересует вопрос: можно ли в такой ситуации утверждать, что (1) Ы (7) в обычном смысле? Оказывается, да, если уравнение замены удовлетворяет дополнительному условию допустимости (которое, как мы увидим ниже, для (4) выполнено).

1.4.2. Утверждение о замене переменных. Рассматриваются уравнения

F(t, J(m)x, C) = 0,

(8)

F(t, J(m)x, C) = 0,

(9)

и уравнение замены переменных

Z(t, J(p)x, J(q)y, C) = 0,

(10)

которое наряду с x содержит новую неизвестную функцию y. Предполагается, что (8) Ю(10) (9), т. е. (8) Щ (10) Ю (9). Предполагается, далее, что замена переменных допустима для уравнения (1), т. е. любое решение x = j(t) уравнения (8) можно дополнить некоторой функцией y = y(t) до решения (j, y) уравнения (10).

Утверждается, что тогда (8) Ю (9) в обычном смысле.

Д о к а з а т е л ь с т в о.  Пусть x = j(t) — решение уравнения (8). Дополним его до решения x = j(t), y = y(t) уравнения (10). Поскольку (8) Щ (10) Ю (9), полученная пара есть решение уравнения (9). Но поскольку в (9) y не входит, это означает, что j — решение уравнения (9).

1.4.3. Утверждение о линейном неоднородном уравнении. При выполнении условия (2) общее решение уравнения (1) задается формулой (7), где функция Ft0(t) определена равенством (5).

Д о к а з а т е л ь с т в о.  Основная часть доказательства уже проведена в п. 1.4.2. В соответствии с утверждением о замене переменных нам остается проверить, что замена (4) допустима для (1) и для (7). Но она, очевидно, допустима для любого уравнения с неизвестной x, так как Ft0(t) № 0 и из (4) C легко находится по x: C = Ft0–1(t)x.

Формула (7) играет важную роль в теории обыкновенных дифференциальных уравнений. Константа C в ней имеет простой смысл: C = x₀ = x(t₀) — в этом легко убедиться непосредственной подстановкой t = t₀ в (7). В связи с этим саму формулу часто записывают в виде

x = Ft0(t)x0 + Ft0(t)т t t0 F –1 t0 (s)b(s)ds.

(11)

Если a(t) є a не зависит от t, то

Ft0(t) = ea(t – t0)

и из (11) мы получаем еще одну часто используемую формулу:

x = ea(t–t0)x0 + т t t0 ea(t–s)b(s) ds.

(12)

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

1.4.4. Пример на понижение порядка. Уравнение

G(x(k), x(k+1), ..., x(m)) = 0   (k і 0)

(13)

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

t, x, xў, ..., x(k), x(k+1), ..., x(m).

В таком случае часто рекомендуется замена

x(k) = q, x(k+1) = p

(очевидно, допустимая для любого уравнения порядка m і k + 1), которая превращает (13) в уравнение порядка m – k + 1 с двумя неизвестными функциями:

G(q, p, pў, ..., p(m – k + 1)) = 0.

(14)

При его решении можно пользоваться очевидным соотношением

gў = p.

(15)

Если затем удается преобразовать (14) к виду

g(t, q, p, C) = 0,

то для x получается уравнение порядка k + 1:

g(t, x(k), x(k + 1), C) = 0.

Например, уравнение гармонического осциллятора (см. п. 1.3.10)

xўў + w2x = 0

имеет вид (13) с k = 0. Замена

x = q, xў = p

приводит его к виду

dp dt + w2q = 0.

(16)

Если учесть (15), то при p № 0 получим dt = dq/p, поэтому (16) перепишется в форме

p·dp + w2q·dq = 0.

Это уравнение с разделенными переменными легко решается:

p2 + w2q2 = C.

Вернувшись к старым обозначениям, получим для x уравнение первого порядка:

(xў)2 + w2x2 = C,

которое можно далее разрешить относительно xў и затем проинтегрировать как уравнение с разделяющимися переменными. Нетрудно показать, что общее решение уравнения гармонического осциллятора можно записать в виде

x = C1cos(wt + C2) (C1 і 0, C2 О [0,2p)).

(17)

Уравнение гармонического осциллятора относится к важному классу линейных уравнений с постоянными коэффициентами, для которых позже будет описан (см. § 3.5) универсальный алгебраический алгоритм отыскания общего решения.

1.4.5. Сведение уравнения m-го порядка к системе m уравнений первого порядка. Рассмотрим уравнение m-го порядка

f(t, x, xў, ..., x(m–1), x(m)) = 0.

(18)

Это уравнение может быть сведено к системе m дифференциальных уравнений первого порядка с помощью замены

y1 = x, y2 = xў, ... ym-1 = x(m–2), ym = x(m–1).

(19)

Формально продифференцировав уравнения (19) и воспользовавшись уравнениями (18) и (19), получим следующую систему уравнений

yў1 = y2, yў2 = y3, ... yўm-1 = ym, f(t, y1, y2, ..., ym, yўm)= 0.

(20)

Покажем, что уравнения (18) и (20) эквивалентны с учетом замены (19). Действительно, пусть x = j(t), y₁ = y₁(t), ..., y_m = y_m(t) — решение системы (18) Щ (19). Проделав, уже не формально (дифференцировать теперь можно), описанные при выводе системы (20) действия, получим, что y₁ = y₁(t), ..., y_m = y_m(t) — решение системы (20). Импликация (20) Щ (19) Ю (18) доказывается аналогично. Остается заметить, что замена (19) очевидно допустима для уравнения (18) и (20).

1.4.6. Некоторые рекомендуемые замены. В таблице описываются некоторые из наиболее распространенных типов замен переменных. Как правило, их применение (т. е. переход от (8) к (9) с учетом (10) и обоснование допустимости замены) не вызывает затруднений.

1.4.7. Контрольные вопросы

1.4.7.1. Приведите пример недопустимой замены переменных.

1.4.7.2. Транзитивно ли отношение следования с учетом замены (относительно одной и той же замены З), т. е. вытекает ли из утверждения (У1 ЮЗ У2) Щ (У2 ЮЗ У3) утверждение У3 ЮЗ У3?

1.4.7.3. Правильно ли утверждение: если (У1 ЮЗ1 У2) Щ (У2 ЮЗ2 У3), то найдется замена  З3 такая, что У3 ЮЗ3 У3?

1.4.7.4. Сведите уравнение колебаний математического маятника xўў + w²sin x = 0 к нормальной системе дифференциальных уравнений первого порядка.

1.4.7.5. Если два решения уравнения (1) совпадают в одной точке, то они совпадают на общей области определения?

1.4.7.6. Докажите (17).

1.4.8. Задачи

1.4.8.1. Докажите, что если в линейном уравнении

xў = ax + b(t)

a — отрицательная константа, а b : R ® R — непрерывная функция, имеющая конечный предел на бесконечности: lim_t®Ґb(t) = b₀, то при t ® Ґ каждое решение этого уравнения стремится к b₀/a.

1.4.8.2. Докажите, что если в уравнении txў = ax + b(t) константа a положительна, а предел b₀ = lim_t®0b(t) существует, то все решения этого уравнения при t ® 0 имеют один и тот же конечный предел. Найдите его.

1.4.8.3. Пусть в уравнении (1) a(t) Ј c < 0, а lim_{t ® 0}b(t) = 0. Докажите, что все решения этого уравнения стремятся к нулю при t ® + Ґ.

1.4.8.4. Докажите, что если в уравнении (1) a(t) Ј c < 0, то разность любых двух его решения при t ® + Ґ (а если a(t) і c > 0, то при t ® -Ґ) стремится к нулю.

1.4.8.5. Покажите, что уравнение xў + x = b(t), где |b(t)| Ј M при всех t О R имеет единственное ограниченное на R решение. Докажите, что это решения является периодической функцией, если b периодическая.

1.4.8.6. Докажите, что если два различных решения уравнения (1) являются T-периодическими, то любое решение этого уравнения T-периодично.

1.4.8.7. Докажите, что если два различных решения уравнения (1) ограничены на всей оси, то таковыми являются все его решения.

1.4.8.8. Покажите, что уравнение xў = x·sin t + 1 не имеет периодических решений.

1.4.8.9. Пусть x = j_i(t) (i = 1, 2) — решение уравнения (1), в котором b(t) = b_i(t) (b_i: [0, +Ґ) ® R — непрерывные функции). Докажите, что для любого e > 0 найдется такое d > 0, что если

|j1(0) – j2(0)| < d  и  |b1(t) – b2(t)| < d при всех t і 0,

то

|j1(t) - j2(t)| < e при всех t О D(j1) З D (j2).

1.4.8.10. Покажите, что уравнение Рикатти

xў = a(t)x + x2 + c(t)

(21)

с непрерывными периодическими коэффициентами не может иметь более двух периодических решений.

1.4.8.11. Покажите, что любые четыре решения уравнения Рикатти (21) связаны соотношением

j4(t) – j2(t) j4(t) – j1(t) / j3(t) – j2(t) j3(t) – j1(t) = const.

1.4.8.12. Покажите, что если все интегральные кривые уравнения xў = f(t, x) подобны с центром подобия в начале координат (t, x)-плоскости, то это уравнение однородно: f(at, ax) є f(t, x).

1.4.8.13. Докажите, что если x = j(t) — ненулевое решение линейного уравнения m-го порядка

x(m) + am–1x(m–1) + ... + a1(t)xў + a0(t)x = 0,

то после замены x = yj(t) это уравнение переходит в уравнение, не содержащее первой производной, порядок которого, как известно, может быть понижен на единицу.

File based on translation from T_EX by T_TH, version 2.32.
Created 2 Jan 2002, 16:43.
Last modified 8 Apr 2002.