Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

matan (1)

.pdf
Скачиваний:
10
Добавлен:
01.06.2015
Размер:
819.57 Кб
Скачать

L(y) = f(x) (1’)

Если все коэффициенты P1(x), f(x) непрерывна на [a,b], то по теорема существования и единственности для ОДУ n-го порядка существует единственное решение задачи Коши для этого ДУ с любыми НУ в точке x0 принадлежащей [a,b]

Свойство1: пусть y(x) – решение уравнения

y(n) + P1(x)y(n-1) + … + Pn-1(x)y’ + Pn(x)y = f(x) (1)

а y0(x) – решение соответствующего ЛОДУ L(y) = 0 , тогда y0(x) + y(x) – решение уравнения

(1)

Доказательство:

L[y0(x) + y(x)] = L[y0] + L[y(x)] = f(x)

y0(x) +y(x) – решение (1)

Свойство2: пусть y1(x), … , ym(x) решения L(y) = fi(x) , i=1, … , m, тогда

y(x) =∑

iyi(x) – является решением ДУ L(y) = ∑

ifi(x)

Доказательство вытекает из свойства 1.

Свойство3: (структура общего решения ЛНОДУ n-го порядка)

Общее решение на отрезке [a,b] L(y) = f(x) с непрерывными на этом же отрезке коэффициентами Pi(x) и f(x) равно сумме общего решения соответствующего ЛОДУ n-го

порядка ∑

 

 

 

 

 

 

) и какого нибудь частного решения

 

 

y(x) этого ЛНОДУ n-го порядка,

т.е. y(x)= ∑

) +

 

 

(1)

 

 

 

 

 

y(x)

 

 

 

Доказательство: y(x)= ∑

) +

 

 

 

y(x) – решение ЛНОДУ следует из свойства 1.

 

Докажем что это общее решение: любые НУ: y(x0)=y0 , y’(x0) = y’0 ,…, y(n-1)(x0)=

) (2)

Удовлетворим (1) этим начальным условиям:

 

 

 

 

 

 

 

 

 

 

 

 

(x0) + y(x0)=y0 - первое НУ

 

СЛАУ из n уравнений

 

 

 

 

 

(x0) + y’(x0) =y’0 – второе НУ

 

с n неизвестными С12,…,Сn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

)(x0) +

 

(n-1)(x0) =

) - последнее НУ

 

 

 

i

y

 

 

 

Единственное решение, т.к. = (x0) ≠ 0,

т.к. система линейно независима.

14.Теарема о структуре общего решения ЛНОДУ n-го порядка.

(y(x) – пишется с волнистым подчеркиванием сверху)

Свойство1 ЛНОДУ n-го порядка: пусть y(x) – решение уравнения

y(n) + P1(x)y(n-1) + … + Pn-1(x)y’ + Pn(x)y = f(x) (1)

а y0(x) – решение соответствующего ЛОДУ L(y) = 0 , тогда y0(x) + y(x) – решение уравнения

(1)

Доказательство:

L[y0(x) + y(x)] = L[y0] + L[y(x)] = f(x)

y0(x) +y(x) – решение (1)

Свойство2 ЛНОДУ n-го порядка: пусть y1(x), … , ym(x) решения L(y) = fi(x) , i=1, … , m, тогда

y(x) =∑

iyi(x) – является решением ДУ L(y) = ∑

ifi(x)

Доказательство вытекает из свойства 1.

 

Свойство3 ЛНОДУ n-го порядка: (структура общего решения ЛНОДУ n-го порядка)

Общее решение на отрезке [a,b] L(y) = f(x) с непрерывными на этом же отрезке коэффициентами Pi(x) и f(x) равно сумме общего решения соответствующего ЛОДУ n-го

порядка ∑

 

 

 

 

 

 

) и какого нибудь частного решения

 

 

y(x) этого ЛНОДУ n-го порядка,

т.е. y(x)= ∑

) +

 

 

(1)

 

 

 

 

 

y(x)

 

 

 

Доказательство: y(x)= ∑

) +

 

 

 

y(x) – решение ЛНОДУ следует из свойства 1.

 

Докажем что это общее решение: любые НУ: y(x0)=y0 , y’(x0) = y’0 ,…, y(n-1)(x0)=

) (2)

Удовлетворим (1) этим начальным условиям:

 

 

 

 

 

 

 

 

 

 

 

 

(x0) +

y(x0)=y0 - первое НУ

 

СЛАУ из n уравнений

 

 

 

 

 

(x0) + y’(x0) =y’0 – второе НУ

 

с n неизвестными С12,…,Сn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

)(x0) +

 

(n-1)(x0) =

) - последнее НУ

 

 

 

i

y

 

 

 

Единственное решение, т.к. = (x0) ≠ 0,

т.к. система линейно независима.

15.Нахождение частных решений ЛНОДУ n-го порядка по виду правой части.

(Pm, y, Qm – пишутся с волнистым подчеркиванием сверху)

a0y(n) + a1y(n-1) + … + an-1y’ + any =f(x) (1)

1. пусть f(x) = Pm(x) = A0xm + A1xm-1+ … + Am-1x + Am

а) пусть an≠0, 0 не является корнем характеристического уравнения

y(x) = Qm(x) = B0xm + B1xm-1 + … + Bm (2)

B0, B1 , … ,Bm – находятся методом неопределенных коэффициентов:

Выражение (2) подставляется в выражение (1) и приравниваются коэффициенты с одинаковыми степенями:

xm anB0 = A0 => B0 = A0/an

xm-1 anB1 + mB0 = A1 => B1 = (A1 – mB0)/an

и т.д.

б) an=0, это означет что 0 является корнем характеристического уравнения кратности α

y(x) = xα(B0xm + B1xm-1 + … + Bm), Bi находятся аналогично случаю а

2. пусть f(x) = ePx(A0xm + … + Am)

Является ли Р корнем характеристического уравнения?

а) Р является корнем х.у. =>

y(x) = ePx(B0xm + … + Bm)

Bi – определяются методом неопр коэффициентов б) Р является корнем х.у. кратности α

y(x) = xα ePx(B0xm + … + Bm)

Bi – определяются методом неопр коэффициентов

3. пусть f(x) = ePx(Pm(x)cosqx + Qm(x)sinqx), где Pm(x) ,Qm(x) – многочлены из которых одинр степени m, а другой степени не большей m.

Вопрос: является ли p+iq корнем х.у.?

а) p+iq – не является корнем х.у.

y(x) = ePx(Pm(x)cosqx +Qm(x)sinqx), где Pm , Qm – многочлены степени m с неопределенными коэффициентами , которые определяются методом неопределенных коэффициентов

б) р+iq является корнем х.у. кратности α

y(x) = xα ePx(Pm(x)cosqx +Qm(x)sinqx), где Pm , Qm – многочлены степени m с

неопределенными коэффициентами , которые определяются методом неопределенных коэффициентов

16. Решение ЛНОДУ n-го порядка методом вариации постоянных

Пусть имеем ЛНОДУ n-го порядка с некоторыми коэффициентами:

Yn+P1(x)yn-1+…+Pn-1(x)y’+Pn(x)y=f(x) (1)

 

 

 

 

 

 

Пусть найдено общее решение соответствующееЛОДУ n-го порядка:

 

y(x) = ∑

)

 

(2)

 

 

 

 

 

 

 

Ищем решение в виде y(x)= ∑

)

)

 

(3)

 

 

 

Нахождение Ci(x):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

СЛАУ относительно Ci’(x),i=1,n

 

 

 

 

 

 

 

 

)

) + ∑

)

 

),

где ∑

)

 

)

 

)

 

) + ∑

 

)

 

) , где ∑

)

)

 

 

 

.

 

 

 

 

 

 

 

 

 

 

)

 

) + ∑

)

 

) , где ∑

)

)

)

{

 

 

 

 

 

 

 

 

 

 

 

 

имеет единственное решение, т.к.

=W(x0)

0=>находим Ci(x)

 

 

17. Нормальные системы ОДУ 1-го порядка. Общая теория

Нормальные системы ОДУ 1-го порядка это системы вида:

, , , , )

, , , ,

)

)

),

), ,

)

,, , , )}

Начальное условие: x1(t0)=x10 , x2(t0)=x20,…,xn(t0)=xn0

(2)

Определение 1: Решением нормальной системы (1) называется набор функций: ((x1(t),x2(t),…,xn(t)), который удовлетворяет каждому уравнению системы (1)

Задача (1)-(2) – задача Коши.

Теорема: (Существования и единственности задачи Коши)

1) Пусть fi непрерывна в окрестности начального условия, т.е. М0(t0,x10,…,xn0)

2) Пусть

 

 

 

 

 

, где i,j=1,n

ограничены в некоторой окрестности точки M0 => (существует и

 

при том единственное) решение задачи Коши (1)-(2). Система может быть записана в матричном виде:

 

 

 

)

 

) (

 

)

 

 

 

)

 

)

 

 

условие записи в матричном виде

 

 

 

 

 

 

 

 

 

 

 

 

 

(

 

)

 

(

 

)

 

 

 

 

}

x’(t)=F

(4) - матричный вид.

 

 

18. Приведение нормальной системы ОДУ 1-го порядка к одному ОДУ высшего порядка

, , , , )

, , , ,

)

)

,, , , )

Исключим из уравнения (1) все неизвестные функции, кроме одной, для которой и получим уравнение более высокого порядка:

, , , ,

)

)

исключим

Дифференцируем по переменной t:

+ ∑

f(i)

Ф2(t,x1,…,xn) (3)

+ ∑

f(i)

Ф3(t,x1,…,xn) (4)

, , , )

)

, , ,

)

)

Из системы (3)-(4)-(6) находим x2,x3,…,xn, выражая их через t,x’1,x’’2,…,xn и подставляем их в уравнение (6), в результате получим:

, ,

, ,

)

)

19.Приведение ОДУ высшего порядка к нормальной системе ОДУ 1-го порядка.

20.Системы линейных ОДУ 1-го порядка с постоянными коэффициентами. Метод Эйлера для решения этих систем.

Система уравнений вида

(1)

называется неоднородной системой линейных дифференциальных уравнений с постоянными коэффициентами. Будем считать, что являются непрерывными функциями на (a,b).

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

 

 

 

 

,

 

 

(2)

 

 

 

 

называется

однородной.

Вводя

в

рассмотрение

векторы

 

 

и

матрицу

,

уравнения (1),(2) можно представить в векторной форме

 

 

 

(1')

(2')

Матрица

,

(3)

где - координаты линейно независимых решений (векторов)

...........................

векторного уравнения (2'), называется фундаментальной матрицей этого уравнения.

Иногда ее называют матрицей Вронского.

Определитель

,

составленный из частных решений системы (2), называется определителем Вронского.

Для того, чтобы матрица (3), где - частные решения системы уравнений (2), была фундаментальной, необходимо и достаточно, чтобы при . При этом общее решение векторного уравнения (2') представляется в виде ,

где C - произвольный постоянный вектор. Общее же решение уравнения (1') будет

,

где - какой-нибудь вектор, являющийся частным решением уравнения (1').

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

несложных систем.

Метод Л. Эйлера

y=ekx (2)

y'= k ekx

y''=k2 ekx и так далее

yn = kn ekx

a0kn ekx + a1k(n-1) ekx + ... + an-1k ekx + an ekx = 0

a0kn + a1kn-1 + ... + an = 0 - характеристическоеуравнение

Случаи:

1) Все корни характеристического уравнения действительныеи различные k1, k2, ..., kn

ek1x, ek2x, ..., eknx - линейно независимые решения

y(x) = C1 ek1x + C2 ek2x + ... + Cn eknx

2) k1,2 = α +

e(α+iβ)x = eαx * eiβx =eαx( cosβx + sinβx)

ei = cos + isin - формула Эйлера

e(α-iβ)x = eαx( cosβx - (если спросит, должно ли здесь быть что-то, то скажите забыл(а) i дописать)sinβx)

eαx * cosβx, eαx * sinβx - линейно независимые решения

3) к - корень характеристического уравнения кратности α

ekx, x ekx, ... xα-1 * ekx - линейно независимы на любом отрезке [a, b] 4) k = p - iq - комплексный корень кратности α, то есть повторяется α раз k = p - iq - кратности α

epxcosqx, x epxcosqx, ... , xα-1 epxcosqx

epxsinqx, x epxsinqx, ... , xα-1 epxsinqx

α - линейно независимые решения

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]