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

Математика для инженеров(теория)I том

.pdf
Скачиваний:
128
Добавлен:
18.02.2016
Размер:
4.09 Mб
Скачать

Запишем его так: (λ − 2)(λ2 +1)2 = 0 . Отсюда получаем λ1 = 2 – однократный корень и λ2,3,4,5 = ±i пара двукрат- ных мнимых корней. Общее решение есть

y= C1e2x + (C2 + C3x)cos x + (C4 + C5x)sin x .

Вслучае, когда характеристическое уравнение (3) имеет кратные

вещественные и комплексные корни, нужно объединить случаи 2 и 4.

§ 4. Линейные неоднородные уравнения n-го порядка

Линейное неоднородное уравнение n-го порядка имеет

вид

yn + p (x)y(n−1)

+ p

2

(x)y(n−2) + ... + p

n

(x) y = f (x)

(1)

1

 

 

 

 

 

 

или

L[ y] = f (x) ,

 

 

(2)

 

 

 

 

 

где pi (x) , i =

 

и

f (x) непрерывны на I = (a;b) .

 

1,n

 

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

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

(1) представляется в виде суммы какого-нибудь частного решения y * этого уравнения и общего решения y соот-

ветствующего однородного уравнения, т.е. y = y * + y .

Доказательство. Из условия теоремы 1 вытекает, что L[ y*] = f (x) , а L[ y] = 0 . Тогда, в силу линейности оператора,

имеем L[ y * + y] = L[y*] + L[y] = f (x) + 0 = f (x) . Это

значит, что

функция

y = y * + y

есть решение уравнения (2).

 

 

 

Для доказательства того, что функция

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

y = y * + y = y * +åCi yi ,

(3)

 

 

 

 

i=1

 

 

 

где

Ci , i =

 

 

произвольные постоянные,

yi ,i =

 

, –

1,n

1,n

фундаментальная система решений соответствующего од- нородного уравнения, необходимо показать, что любое решение уравнения (2) можно получить из формулы (3) соответствующим выбором постоянных C1,C2 ,...,Cn . Это

543

можно сделать по аналогии с доказательством теоремы

2.1.

Упражнение 1. Доказать, что функция (3) является общим решением уравнения (1).

Для отыскания частного решения линейного неоднородного уравнения (1) применяется метод вариации произвольных постоян- ных, который заключается в следующем. Пусть известна фундамен- тальная система y1, y2 ,..., yn решений соответствующего однородного

уравнения. Тогда общее решение неоднородного уравнения (1) будем

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

 

 

 

 

n

 

 

 

 

 

y = åCi (x)yi ,

(4)

 

 

 

 

i=1

 

 

но с переменными коэффициентами. Для определения

Ci (x) , i =

 

,

которых n, а уравнение одно, требуется на-

1,n

ложить (n −1)

условие при вычислении производных функ-

ции y, а именно:

n

 

n

 

 

 

 

 

 

 

 

y′ = åCi′(x)yi +

åCi (x) yi′,

 

 

 

 

i=

 

i=1

 

 

 

 

14243

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

n

 

n

 

 

 

y′′ = åCi′(x) yi′ + åCi (x) yi′′ ,

 

 

 

 

i=

 

i=1

 

 

 

 

14243

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

........................................,

 

 

 

 

14243

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

n

 

n

 

 

 

y(n)

= åCi′(x)yi(n−1) + åCi (x) yi(n) .

 

 

 

 

i=1

i=1

Подставляя эти производные в уравнение (1), полу-

чим еще одно уравнение для нахождения Ci′(x) :

 

 

 

 

 

n

 

 

 

 

Ln[ y] = Ln[ y] +

åCi′(x) yi(n−1) = f (x) .

 

 

 

{

i=1

 

 

 

 

 

 

 

 

 

 

0

 

 

 

544

Будем иметь следующую систему линейных алгеб-

раических уравнений

 

 

для

определения

функций Ci′(x) ,

i =

 

:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1,n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Cy + Cy

2

 

+ ... + Cy

n

= 0,

 

 

 

 

1 1

2

 

 

 

n

 

 

 

 

 

 

 

 

 

Cy

+ Cy

 

+ ... + Cy

= 0,

 

 

 

 

1 1

2

2

 

 

 

n

n

 

 

 

 

 

 

(5)

........................................,

 

 

 

 

 

 

 

 

 

Cy(n−1) + C

y(n−1)

+ ... + Cy(n−1)

=

f (x).

1 1

 

 

2

2

 

 

 

 

n n

 

 

 

Решение системы (5) можно найти, например, по формулам

Крамера (§1.7):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ci′ =

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

, i =1,n ,

(6)

 

 

 

W[y1,..., yn ]

 

 

 

 

 

 

 

 

 

 

причем определителем этой системы является определитель Вронского.

Интегрируя (6), получаем функции Ci (x),i =1,n .

Пример 1. Решить уравнение y′′′ − 1x y′′ + y′ − 1x y = x , если

известно, что x,cos x,sin x фундаментальная система ре- шений соответствующего однородного дифференциально- го уравнения.

Решение. Согласно методу вариации произвольных постоянных, общее решение заданного уравнения будем искать в виде (4), т.е.

y = C1(x)x + C2 (x)cos x + C3 (x)sin x .

Вычислим определители:

 

x

cos x

sin x

 

 

 

 

W[x,cos x,sin x] =

1

−sin x

cos x

 

= x ,

 

0

−cos x

−sin x

 

 

 

0

cos x

sin x

 

 

x

0

sin x

 

1 =

0

−sin x

cos x

= x ,

2 =

1

0

cos x

= −x(xcos x − sin x) ,

 

 

x

−cos x

−sin x

 

 

0

x

−sin x

 

 

 

x

cos x

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3 =

1

−sin x

0

 

= x(−xsin x − cos x) .

 

 

 

 

0

−cos x

x

 

 

 

 

 

 

 

545

 

По

формулам

(6) вычислим

Ci′(x) : C1′(x) =1,

C

(x) = -xcos x + sin x , C

(x) = -xsin x - cos x .

Тогда

2

 

3

 

 

 

 

C1(x) = x + C1 ,

 

 

C2 (x) = -òxcos xdx - cos x + C2 = -xsin x - 2cos x + C2 ,

 

C3 (x) = -òxsin xdx - sin x + C3 = xcos x - 2sin x + C3.

 

Окончательно получим:

 

 

y = (x + C1)x + (-xsin x - 2cos x + C2 )cos x + (x cos x - 2sin x + C3)sin x =

 

2

 

 

 

 

= {x - 2 + C x + C2 cos x + C3 sin x.

 

 

y

14444244443

 

 

 

y

 

§ 5. Линейные неоднородные дифференциальные

уравнения n-го порядка с постоянными коэффициентами

 

Рассмотрим линейное неоднородное уравнение n-го

п

о

р

я

д

к

а

 

 

y(n) + a1y(n−1) + ... + an−1y¢ + an y = f (x) ,

 

(1)

где a1,a2 ,...,an постоянные вещественные коэффици-

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

частного решения и общего решения соответствующего ему однородного уравнения. Общее решение однородного уравнения находится методом Эйлера, описанным в § 4. Таким образом, интегрирование (1) сводится к отысканию его частного решения y *. В общем случае это можно осуществить ме-

тодом вариации произвольных постоянных.

Пример 1. Решить уравнение y′′′ + y= tgx .

 

Решение.

Составим

характеристическое

уравнение

λ3 + λ = 0 .

Имеем

λ(λ2 +1) = 0 Þ λ

= 0, λ = ±i ;

 

 

1

2,3

y = C1 + C2 cos x + C3 sin x

общее решение соответствующего

однородного уравнения.

 

Для определения функций Ci′(x),i =1,2,3, вычислим оп-

ределитель Вронского

и дополнительные определители

546

Di ,

 

i =1,2,3

и подставим их значения в формулу (4.6). По-

лучим:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

cos x

sin x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

W[1,cos x,sin x] = D =

 

0

-sin x

cos x

=1,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

-cos x

-sin x

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

cos x

sin x

 

 

= tg x Þ C1¢(x) = tg x ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D1 =

0

 

 

 

-sin x

cos x

 

 

 

 

 

 

 

 

 

tg x -cos x

-sin x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

0

 

sin x

 

 

 

 

 

 

 

 

¢

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D2 =

0

0

 

cos x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= -sin x Þ C2 (x) = -sin x ,

 

 

 

 

 

 

 

 

 

0

tg x

-sin x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

cos x

 

 

0

 

 

 

= -tg x sin x Þ C3¢ (x) = - tg x sin x .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D3 =

0

-sin x

 

0

 

 

 

 

 

 

 

 

 

 

 

0

-cos x

 

tg x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

После

 

 

 

 

 

 

 

 

 

интегрирования

получим:

C1(x) = òtg xdx = -ln | cos x | +C1 ,

C2 (x) = -òsin xdx = cos x + C2 ,

C (x) = -

sin

2 x

dx =

 

æcos x -

 

 

1

ö dx = sin x + +ln

 

tgæ

x

- π ö

 

+ C .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

ò cos x

 

ò

ç

 

 

 

 

÷

 

 

ç

 

÷

 

3

 

 

 

 

è

 

 

 

 

 

cos x ø

 

 

è 2

4 ø

 

 

Тогда

 

 

 

 

 

 

 

 

y = C1(x) + C2 (x)cos x + C3 (x)sin x = −ln | cos x | +1+

+ln

 

æ x

-

π

ö

 

sin x

+ C1 + C2 cos x + C3 sin x .

 

 

 

 

 

 

 

 

 

tgç

 

4

÷

 

 

 

 

2

 

 

 

 

 

è

 

 

ø

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В некоторых случаях для определения частного ре- шения уравнения (1) применяется метод подбора.

Общий вид правой части f (x) уравнения (1), при котором можно

применить метод подбора (метод неопределенных коэф- фициентов) следующий:

f (x) = eα x[P (x)cos β x + Q (x)sin β x] ,

(2)

l

m

 

где Pl (x) и Qm (x) –

многочлены степени l и

m соот-

ветственно.

В этом случае частное решение y * уравнения (1) ищется в

виде

y *(x) = x

s α x

%

%

(3)

e

[Pk (x)cos β x + Qk (x)sin β x] ,

547

где

%

%

k = max(l,m) – многочлены

k -ой степе-

Pk (x),Qk (x) ,

ни от

x

с неопределенными коэффициентами, а s крат-

ность

корня

λ = α + iβ

характеристического

уравнения,

причем,

если

α + iβ

не

является корнем характеристиче-

ского уравнения, то s = 0 . Для удобства приведем таблицу видов частных решений для различных правых частей (2) уравнения (1).

Здесь первые три вида правых частей являются част- ными случаями вида IV.

 

 

 

 

 

 

 

Таблица 1

 

Правая

 

 

Корни

 

Виды

 

часть диф-

характеристического

 

частного реше-

 

ференциального

 

уравнения

 

 

ния

 

 

уравнения

 

 

 

 

 

 

 

 

 

 

1.

Число

0

 

 

 

 

не является

кор-

Pm (x)

 

 

нем

характери-

 

 

 

 

 

 

 

%

 

 

стического

 

 

 

 

 

Pm (x)

уравнения

 

 

 

 

 

 

2.

Число 0

 

 

 

 

 

 

 

 

 

корень

характе-

x Pm (x)

 

 

ристического

 

 

 

 

 

 

 

 

s

%

 

 

уравнения

крат-

 

 

 

 

ности s

 

 

 

 

 

 

 

1.

Число

α

 

 

 

 

не

является кор-

 

 

 

 

нем

характеристи-

%

α x

 

 

Pm

(x)e

 

 

ческого

уравне-

 

 

I

P (x)eα x

ния

 

 

 

 

 

 

m

 

2.

Число

α

 

 

 

 

 

 

 

 

является корнем

 

 

 

 

характеристиче-

s %

α x

 

 

x Pm (x)e

 

 

ского

уравнения

 

 

 

 

кратности s

 

 

 

 

 

 

 

 

Продолжение

таблицы 1

548

 

 

 

 

 

1.

Число

+iβ

 

 

%

(x)cos β x +

 

 

 

 

не

является

кор-

 

 

Pk

 

 

 

 

 

 

 

%

 

 

 

 

нем

характеристи-

 

 

+Qk (x)sin β x,

 

 

 

 

 

 

k = max(l, m)

 

 

 

 

ческого

 

уравне-

 

 

 

P (x)cos β x +

 

 

 

 

 

 

 

+Q (x)sin β

 

ния

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

II

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2.

Число

+iβ

 

 

x

(Pk (x)cos β

 

 

 

 

 

 

 

 

 

 

является

корнем

 

 

 

 

 

 

характеристиче-

 

 

s

%

 

 

 

 

 

 

+Qk (x)sin β x)

 

 

 

 

 

 

 

 

 

 

 

 

%

 

 

 

 

ского

уравнения

 

 

 

 

 

 

 

 

кратности s

 

 

 

 

 

 

 

 

 

 

1.

Число

α + iβ

 

 

%

 

 

 

 

не является корнем

 

 

(Pk (x)cos β x +

 

 

 

 

 

 

 

%

 

 

 

 

характеристическо-

 

 

+Qk (x)sin β x)

 

 

 

 

 

 

 

 

 

eα x (P (x)cos

 

го уравнения

 

 

 

 

 

 

 

 

2.

 

Число

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V

+Qm (x)sin β

 

α + iβ

является

 

 

 

 

 

 

 

 

 

s

 

 

 

 

 

корнем

характери-

 

 

%

 

 

 

 

 

 

x

(Pk (x)cos β

 

 

 

 

стического

урав-

 

 

+Qk (x)sin β x)

 

 

 

 

 

 

 

 

 

 

 

 

%

 

 

 

 

нения

кратности

 

 

 

 

 

 

 

 

s

 

 

 

 

 

 

 

 

 

Пример

2.

 

Найти

общее

решение

уравнения

y′′′ − y′′ + y′ − y = x2 + x .

 

 

 

 

 

 

 

 

 

 

Решение.

 

 

Характеристическое

 

уравнение

λ3 − λ2 + λ −1 = 0

имеет различные корни λ

=1, λ

= −i, λ = i .

 

 

 

 

 

 

 

 

1

2

 

3

Поэтому общее решение соответствующего однородного уравнения y будет:

y = C1ex + C2 cos x + C3 sin x .

Так как число 0 не является корнем характеристиче- ского уравнения, то частное решение исходного уравнения y * надо

искать в виде (табл. 1, случай I):

y* = A1x2 + A2 x + A3 ,

549

где A1, A2 , A3

неизвестные пока коэффициенты. Под-

ставляя выражение для y *

в данное уравнение, получаем:

-A x2

+ (2A - A )x + (A - 2A - A ) = x2

+ x .

1

1

2

2

1

3

 

Отсюда

 

 

 

 

 

 

 

ìA = -1,

 

 

 

 

ï

1

 

 

 

 

 

í2A1 - A2 =1,

 

 

 

ïA - 2A - A = 0.

 

 

 

î

2

1

3

 

 

 

Решая эту систему,

найдем A1 = -1, A2 = -3, A3 = -1.

Сле-

довательно, частное решение будет

y* = -x2 - 3x -1,

и общее реше-

ние исходного уравнения имеет вид

 

 

 

 

 

 

 

 

 

 

 

y(x) = C1ex + C2 cos x + C3 sin x - x2 - 3x -1.

 

 

Пример

3. Найти общее

решение

уравнения

y¢¢ + 2y¢ + 5y = ex cos 2x .

 

 

 

 

 

 

 

 

 

 

 

 

Решение.

 

 

Характеристическое

 

 

уравнение

λ2 + + 5 = 0

имеет

корни

λ

= -1± 2i

 

и,

 

 

значит,

 

 

 

 

 

 

 

 

 

1,2

 

 

 

 

 

 

 

 

y = e x (C cos2x + C

2

sin 2x) . Так как число α + iβ = −1+ 2i явля-

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ется

простым

корнем

характеристического

уравнения, то

y *

надо искать

в

виде

(табл.

1

случай

IV):

y* = xex (Acos2x + Bsin 2x) .

Тогда:

 

 

 

 

 

 

 

 

 

 

 

(y*)¢ = ex[(A - Ax + 2Bx)cos2x + (B - Bx - 2Ax)sin 2x],

 

 

(y*)¢¢ = ex[(-2A - 3Ax + 4B - 4Bx)cos2x + (-2B - 3Bx - 4A + 4Ax)sin 2x].

 

Подставляя выражения для функции

y *

и ее произ-

водных в исходное уравнение и сокращая на ex , будем

иметь: −4Asin 2x + 4B cos2x = cos 2x . Отсюда

A = 0, B =

 

1

и,

зна-

4

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

чит,

y* =

xex sin 2x .

 

 

 

 

 

 

 

 

 

 

 

4

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

y(x) = ex (C1 cos2x + C2 sin 2x) + 14 xex sin 2x .

550

§ 6. Линейная краевая задача

10 . Постановка краевых задач. В инженерной прак-

тике часто приходится решать дифференциальное уравне-

ние

Ln[ y] = f (x) ,

(1)

отыскивая его решение, которое на заданном отрезке [a;b] удовлетворяет некоторым дополнительным условиям.

Примером такой задачи служит задача Коши. Характерной особенностью задачи Коши является то обстоятельство, что значение искомой функции y = y(x) и ее производных

до (n −1) -го порядка включительно задаются в одной точке x = x0 данного отрезка. Однако, при решении некоторых

физических и технических задач появляется необходи- мость нахождения решения y = y(x) линейного дифферен-

циального уравнения n-го порядка, если начальные усло- вия заданы на концах некоторого отрезка. Эти условия на- зывают краевыми условиями, и они состоят в том, что на обоих концах некоторого отрезка [a;b] задаются значения искомого

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

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

Краевая задача называется линейной, если дифферен- циальное уравнение и краевые условия линейны. Для про- стоты предположим, что в краевые условия входят две

концевые абсциссы

 

x1 = a

и

x2 = b . Такие краевые условия

называют двухточечными.

 

 

 

 

 

 

Краевые условия называют линейными, если они

имеют вид

 

 

 

 

 

 

 

 

 

Rν (y) = γν

(ν =1, 2, ..., n) ,

 

(2)

n−1

é

 

(k)

(a) + βkν y

(k)

ù

αkν , βkν ν

за-

где Rν (y) = å

ëαkν y

 

 

(b)û ,

k=0

 

 

 

 

 

 

 

 

 

данные

 

n−1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

постоянные, причем

å(kν | + | βkν |) ¹ 0 .

 

 

k=0

551

Так, например, для уравнения второго порядка a0 (x)y′′ + a1(x)y+ a2 (x) y = f (x) , a x b ,

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

 

α1y(a) +α2 y′(a) = γ1 ,

 

β1y(b) + β2 y′(b) = γ 2 .

 

Отметим, что краевая задача не всегда разрешима, а если разре-

шима, то не всегда единственным образом.

 

 

 

 

Пример 1. Найти решение уравнения y′′ + y = 0 , удов-

летворяющее краевым условиям

y(0) = 0 ,

y(a) = y0 .

 

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

вид

y = C1 cos x + C2 sin x .

Условие

y(0) = 0

выполняется, если

C1 = 0 , при этом y = C2 sin x .

 

 

 

 

 

 

 

Если a ¹ π n , где n целое число, то из второго крае-

вого

условия находим

y = C

2

sin a , C

=

y0

. Вследствие

 

 

 

0

 

2

 

sin a

 

 

 

 

 

 

 

этого, существует единственное решение данной краевой

задачи y =

y0

sin x .

 

 

 

 

 

 

sin a

 

 

Если

a = π n , то

из второго

краевого условия имеем

C2 sinπ n = y0 . Значит,

если y0 = 0 ,

то краевая задача имеет

бесконечное множество решений y = C2 sin x , где C2 может принимать любые значения, а при y0 ¹ 0 краевая задача не имеет ре-

шений.

Двухточечные краевые задачи возникают в различного рода технических и физических задачах, в частности, находят широкое применение при расчете балок.

Пример 2. Найти максимальный прогиб консольной балки длины l, на-

груженной равномерно распределенной нагрузкой q (рис 1.)

 

ис. 1

Р

 

 

Решение.

Дифференциаль-

 

 

ное

уравнение

изгиба балки

имеет вид

 

 

 

 

 

 

 

 

 

EI

d 4 y

= q ,

(3)

 

 

 

dx4

 

 

 

 

 

 

552