Метод вариации произвольных постоянных¶

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

¶История

Метод был разработан французским математиком Жозефом Луи Лагранжем в XVIII веке. В 1774–1775 годах Лагранж опубликовал работы, в которых применил этот подход для решения задач небесной механики, в частности, для описания возмущений орбит планет под действием гравитационных сил других небесных тел. Изначально метод назывался «методом вариации орбитальных элементов» или «методом вариации постоянных». Впоследствии он был обобщён на дифференциальные уравнения общего вида и стал одним из фундаментальных инструментов в математической физике, теории колебаний и других областях. В русскоязычной математической литературе метод часто называют методом Лагранжа, а в западной — методом вариации параметров (variation of parameters).

¶Суть метода

Метод вариации произвольных постоянных применяется к линейным неоднородным дифференциальным уравнениям вида:

\[ L[y] = f(x) \]

где \(L\) — линейный дифференциальный оператор, а \(f(x)\) — известная функция (неоднородность). Решение ищется в виде:

\[ y(x) = y_h(x) + y_p(x) \]

где \(y_h(x)\) — общее решение однородного уравнения \(L[y] = 0\), а \(y_p(x)\) — частное решение неоднородного уравнения.

¶Для обыкновенных дифференциальных уравнений (ОДУ)

Рассмотрим линейное неоднородное ОДУ n-го порядка:

\[ y^{(n)} + a_{n-1}(x) y^{(n-1)} + \dots + a_1(x) y' + a_0(x) y = f(x) \]

Пусть известна фундаментальная система решений однородного уравнения \(\{y_1(x), y_2(x), \dots, y_n(x)\}\). Тогда общее решение однородного уравнения записывается как:

\[ y_h(x) = C_1 y_1(x) + C_2 y_2(x) + \dots + C_n y_n(x) \]

где \(C_1, C_2, \dots, C_n\) — произвольные постоянные.

Метод вариации предполагает, что частное решение неоднородного уравнения имеет ту же форму, но с функциями \(C_i(x)\) вместо постоянных:

\[ y_p(x) = C_1(x) y_1(x) + C_2(x) y_2(x) + \dots + C_n(x) y_n(x) \]

Для нахождения \(C_i(x)\) составляется система линейных алгебраических уравнений относительно их производных \(C_i'(x)\):

\[ \begin{cases} C_1'(x) y_1(x) + C_2'(x) y_2(x) + \dots + C_n'(x) y_n(x) = 0 \\ C_1'(x) y_1'(x) + C_2'(x) y_2'(x) + \dots + C_n'(x) y_n'(x) = 0 \\ \vdots \\ C_1'(x) y_1^{(n-1)}(x) + C_2'(x) y_2^{(n-1)}(x) + \dots + C_n'(x) y_n^{(n-1)}(x) = f(x) \end{cases} \]

Эта система имеет единственное решение, так как определитель матрицы коэффициентов (определитель Вронского) для фундаментальной системы решений отличен от нуля. Интегрируя найденные \(C_i'(x)\), получают \(C_i(x)\) (с точностью до произвольных постоянных, которые можно принять равными нулю, так как они дают вклад в общее решение однородного уравнения). Подстановка \(C_i(x)\) в выражение для \(y_p(x)\) даёт искомое частное решение.

¶Для систем линейных дифференциальных уравнений

Метод обобщается на системы линейных неоднородных ОДУ первого порядка. Пусть дана система:

\[ \mathbf{y}' = A(x) \mathbf{y} + \mathbf{f}(x) \]

где \(\mathbf{y}\) — вектор-функция, \(A(x)\) — матрица коэффициентов, \(\mathbf{f}(x)\) — вектор неоднородности. Пусть \(\Phi(x)\) — фундаментальная матрица решений однородной системы (матрица, столбцы которой образуют фундаментальную систему решений). Тогда общее решение однородной системы: \(\mathbf{y}_h = \Phi(x) \mathbf{C}\), где \(\mathbf{C}\) — вектор постоянных. Частное решение ищется в виде \(\mathbf{y}_p = \Phi(x) \mathbf{C}(x)\). Подстановка в исходную систему приводит к уравнению:

\[ \Phi'(x) \mathbf{C}(x) + \Phi(x) \mathbf{C}'(x) = A(x) \Phi(x) \mathbf{C}(x) + \mathbf{f}(x) \]

Учитывая, что \(\Phi'(x) = A(x) \Phi(x)\), получаем:

\[ \Phi(x) \mathbf{C}'(x) = \mathbf{f}(x) \]

Отсюда \(\mathbf{C}'(x) = \Phi^{-1}(x) \mathbf{f}(x)\). Интегрируя, находим \(\mathbf{C}(x)\), а затем \(\mathbf{y}_p(x)\).

¶Для дифференциальных уравнений в частных производных (ДУЧП)

Метод вариации произвольных постоянных применим и к некоторым классам линейных неоднородных ДУЧП, например, к уравнениям первого порядка. В этом случае «постоянные» заменяются на функции от других независимых переменных. Однако для ДУЧП метод менее универсален, чем для ОДУ, и часто требует дополнительных условий (например, использования метода характеристик).

¶Пример применения

Рассмотрим линейное неоднородное ОДУ второго порядка:

\[ y'' + y = \frac{1}{\cos x} \]

Однородное уравнение \(y'' + y = 0\) имеет фундаментальную систему решений: \(y_1 = \cos x\), \(y_2 = \sin x\). Общее решение однородного: \(y_h = C_1 \cos x + C_2 \sin x\).

Ищем частное решение в виде:

\[ y_p = C_1(x) \cos x + C_2(x) \sin x \]

Составляем систему для производных:

\[ \begin{cases} C_1'(x) \cos x + C_2'(x) \sin x = 0 \\ -C_1'(x) \sin x + C_2'(x) \cos x = \frac{1}{\cos x} \end{cases} \]

Решая систему (например, методом Крамера), находим:

\[ C_1'(x) = -\frac{\sin x}{\cos x} = -\tan x, \quad C_2'(x) = 1 \]

Интегрируем:

\[ C_1(x) = -\int \tan x \, dx = \ln|\cos x| + \text{const}, \quad C_2(x) = x + \text{const} \]

Выбирая const = 0 для частного решения, получаем:

\[ y_p = \cos x \cdot \ln|\cos x| + x \sin x \]

Общее решение исходного уравнения:

\[ y = C_1 \cos x + C_2 \sin x + \cos x \cdot \ln|\cos x| + x \sin x \]

¶Сравнение с другими методами

Метод вариации произвольных постоянных является наиболее общим методом нахождения частного решения линейных неоднородных ОДУ, в отличие от метода подбора (метода неопределённых коэффициентов), который применим только для неоднородностей специального вида (многочлены, экспоненты, синусы, косинусы и их комбинации). Метод вариации не требует угадывания формы частного решения и работает для любой непрерывной функции \(f(x)\), при условии, что известно общее решение однородного уравнения. Однако он может быть более громоздким в вычислительном плане, так как требует интегрирования.

¶Применение

Метод широко используется в различных областях науки и техники:

• Теория колебаний: для анализа вынужденных колебаний механических и электрических систем под действием произвольных внешних сил.
• Небесная механика: для расчёта возмущений орбит планет и спутников.
• Теория управления: для решения задач оптимального управления и фильтрации (например, в фильтре Кальмана).
• Математическая физика: для решения неоднородных уравнений теплопроводности, волнового уравнения и уравнений Лапласа (через метод функций Грина, который тесно связан с вариацией постоянных).
• Квантовая механика: в теории возмущений для стационарных и нестационарных задач.

¶Интересные факты

• Метод вариации произвольных постоянных является частным случаем более общего подхода — метода вариации параметров, который применяется не только в дифференциальных уравнениях, но и в вариационном исчислении.
• В русскоязычной литературе метод часто называют «методом Лагранжа», хотя Лагранж разработал его для задач небесной механики, а не для дифференциальных уравнений в их современном виде.
• Для уравнений первого порядка метод вариации сводится к простому интегрированию: если \(y' + p(x)y = q(x)\), то общее решение находится по формуле \(y = e^{-\int p dx} \left( \int q e^{\int p dx} dx + C \right)\), что эквивалентно методу вариации постоянной.

¶Источники

• Эльсгольц Л. Э. «Дифференциальные уравнения и вариационное исчисление». — М.: Наука, 1969.
• Понтрягин Л. С. «Обыкновенные дифференциальные уравнения». — М.: Наука, 1974.
• Арнольд В. И. «Обыкновенные дифференциальные уравнения». — М.: МЦНМО, 2012.
• Тихонов А. Н., Васильева А. Б., Свешников А. Г. «Дифференциальные уравнения». — М.: Физматлит, 2005.