Структура общего решения линейного неоднородного ДУ n-го порядка.



Рассмотрим линейное неоднородное дифференциальное уравнение

y(n) + an-1(x)y(n- 1) + ... + a1(x)y' + a0(x)y = f(x).

Общим решением этого уравнения на отрезке [a;b] называется функция y = Φ(x, C1,..., Cn ), зависящая от n произвольных постоянныхC1,..., Cn и удовлетворяющая следующим условиям :

− при любых допустимых значениях постоянных C1,..., Cn функция y = Φ(x, C1,..., Cn ) является решением уравнения на [a; b] ;

− какова бы ни была начальная точка (x0, y0, y1,0 ,..., yn − 1,0 ) , x0∈ [a;b] , существуют такие значения C1 =C10 , ..., Cn = Cn0 , что функция y = Φ(x, C10 , ..., Cn0) удовлетворяет начальным условиям y(x0) = y0, y '(x0) = y1,0 ,..., y(n − 1) (x0) = yn− 1,0 .

 

Справедливо следующее утверждение (теорема о структуре общего решения линейного неоднородного уравнения).

Общее решение лин.неоднородного ду есть y(x,C1,..., Cn) = y0(x) + y*(x), где y0 – решение лин.однрород. ду соотв.данному, y*(x) – частное решение неоднородного ду.

{ Если все коэффициенты уравнения линейного однородного дифференциального уравнениния непрерывны на отрезке [a;b] , а функцииy1(x), y2(x),..., yn(x) образуют фундаментальную систему решений соответствующего однородного уравнения, то общее решение неоднородного уравнения имеет вид

y(x,C1,..., Cn) = C1 y1(x) + C2 y2(x) + ... + Cn yn(x) + y*(x),

где C1,...,Cn — произвольные постоянные, y*(x) — частное решение неоднородного уравнения. }

Пример.

Рассмотрим линейное неоднородное дифференциальное уравнение с непрерывными на (e, ∞) коэффициентами и непрерывной правой частью

y’’-y = 2-x^2

y’’-y = 0, k^2 – 1 = 0;

k1,2 = +-1, y1 = e^x, y2 = e^-x.

 

частным решением неодродного является функция

y*(x) = A*x^2+B*x + C

y*(x) = x^2

Фундаментальную систему решений соответствующего однородного уравнения образуют функции y1(x) = e^x , y2(x) = e^-x.

Общее решение неоднородного уравнения имеет вид

Y(C1,C2, x) = C1*y1 + C2*y2 + y*(x)

Метод Лагранжа

Задача состоит в вычислении какого–либо частного решения линейного неоднородного дифференциального уравнения с непрерывными коэффициентами и непрерывной правой частью.

Рассмотрим линейное неоднородное дифференциальное уравнение

y(n) + an-1(x)y(n- 1) + ... + a1(x)y' + a0(x)y = f(x).

с непрерывными на [a; b] коэффициентами и непрерывной правой частью.

Предположим, что известна фундаментальная система y1(x), y2(x),..., yn(x) решений соответствующего однородного уравнения y(n) + an-1(x)y(n- 1) + ... + a1(x)y' + a0(x)y = 0.

Будем искать частное решение неоднородного уравнения в виде

y*(x) = C1(x) y1(x) + C2(x) y2(x) + ... + Cn(x) yn(x) ,

где C1(x), C2(x) , ... , Cn(x) — неизвестные, n раз дифференцируемые на [a; b] функции. Их называют варьируемые постоянные общего решения однородного уравнения.

Справедливо следующее утверждение.

Пусть y1(x), y2(x),..., yn(x) — фундаментальная система решений однородного уравнения y(n) + an-1(x)y(n- 1) + ... + a1(x)y' + a0(x)y = 0 с непрерывными на отрезке [a; b] коэффициентами. Если правая часть f(x) неоднородного уравнения y(n) + an-1(x)y(n- 1) + ... + a1(x)y' +a0(x)y = f(x) непрерывна на [a; b], то его частное решение можно искать в виде

y*(x) = y(x,C1,..., Cn) = C1(x) y1(x) + C2(x) y2(x) + ... + Cn(x) yn(x) .

Неизвестные функции C1(x), C2(x) , ... , Cn(x) находятся из системы

Такой метод отыскания частного решения неоднородного уравнения называется методом вариации произвольных постоянных илиметодом Лагранжа.

Подробнее (кому нужно)

Посмотрим как можно найти методом Лагранжа частное решение для уравнения 2–го порядка y'' + a1(x)y' + a0(x)y = f(x) с нерерывными коэффициентами и непрерывной правой частью.

Предположим, что известна фундаментальная система y1(x), y2(x) решений соответствующего однородного уравнения y'' + a1(x)y' + a0(x)y =0.

Будем искать частное решение неоднородного уравнения в виде

y*(x) = C1(x) y1(x) + C2(x) y2(x),

где C1(x), C2(x) — такие неизвестные,дважды дифференцируемые на [a; b] функции.

Для того чтобы подставить функцию y*(x) в исходное уравнение, найдём сначала первую производную y*(x):

(y*(x))' = C1'(x) y1(x) + C1(x) y1'(x) + C2'(x) y2(x) + C2(x) y2'(x).

Будем искать C1(x), C2(x) такими, чтобы C1'(x) y1(x) + C2'(x) y2(x) = 0 и, следовательно,

(y*(x))' = C1(x) y1'(x) + C2(x) y2'(x).

Тогда

(y*(x))'' = C1'(x) y1'(x) + C1(x) y1''(x) + C2'(x) y2'(x) + C2(x) y2''(x).

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

y'' + a1(x)y' + a0(x)yC1'(x) y1'(x) + C1(x) y1''(x) + C2'(x) y2'(x) + C2(x) y2''(x) + a1(x)(C1(x) y1'(x) + C2(x) y2'(x)) + a0(x)(C1(x) y1(x) + C2(x)y2(x))= f(x).

После простых преобразований имеем:

C1(x)(y1'' + a1(x)y1' + a0(x)y1) + C2(x)(y2'' + a1(x)y2' + a0(x)y2) + C1'(x) y1'(x) + C2'(x) y2'(x) = f(x).

Но поскольку y1(x), y2(x) — решения однородного уравнения y'' + a1(x)y' + a0(x)y = 0, то

C1'(x) y1'(x) + C2'(x) y2'(x) = f(x).

Для неизвестных функций C1(x), C2(x) получили систему линейных дифференциальных уравнений второго порядка:

Определитель этой линейной относительно C1'(x), C2'(x) системы — это отличный от нуля на [a; b] вронскиан фундаментальной системы решений.

Следовательно, система имеет единственное решение, которое можно выписать в явном виде (имеем два дифференциальных уравнения первого порядка):

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

Неизвестные варьируемые постоянные найдены — найдено частное решение линейного неоднородного уравнения второго порядка.

Пример

Найдём методом Лагранжа (методом вариации произвольных постоянных) частное решение линейного неоднородного дифференциального уравнения второго порядка

с непрерывными на (e, ∞) коэффициентами и непрерывной правой частью.

Фундаментальную систему решений соответствующего однородного уравнения образуют функции y1(x) = ln x , y2(x) = x.

Будем искать частное решение неоднородного уравнения в виде y*(x) = C1(x) lnx + C2(x) x:

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

Для неизвестных функций C1(x), C2(x) получили систему линейных дифференциальных уравнений второго порядка:

Окуда имеем:

И тогда частным решением исходного уравнения второго порядка является функция y*(x) = C1(x) lnx + C2(x) x:


Дата добавления: 2018-10-27; просмотров: 578; Мы поможем в написании вашей работы!

Поделиться с друзьями:






Мы поможем в написании ваших работ!