Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
BORZENKOV_MATLAB ОДУ.pdf
Скачиваний:
129
Добавлен:
11.05.2015
Размер:
1.91 Mб
Скачать

ГЛАВА 7. ЛИНЕЙНЫЕ УРАВНЕНИЯ ВТОРОГО ПОРЯДКА

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

§1. Линейное уравнение второго порядка

Рассмотрим уравнение следующего вида:

 

 

 

a0(x)y a1(x)y a2(x)y (x), 0 x l,

 

 

 

 

 

(1.1)

в котором коэффициенты a0 (x),a1(x),a2(x) определены и

непрерывны на

отрезке 0 x l, причем производная

a0 (x)

также непрерывна и существует

постоянная

a 0

такая,

что

a0 (x) a.

Функция

 

(x)

предполагается

непрерывной на отрезке 0 x l.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Покажем, что это уравнение можно привести к следующему виду:

 

 

 

 

 

 

 

 

 

d

 

dy

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L[y]

 

 

p(x)

 

 

q(x)y f (x).

 

 

 

 

 

(1.2)

 

 

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

Для этого обе части уравнения (1.1) умножим на (x). Чтобы полученное

уравнение

(x)a0 (x)y (x)a1(x)y (x)a2(x)y (x) (x),

0 x l

можно

было

записать в

виде (1.2),

следует потребовать

выполнения

 

условия

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x a a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(x) exp(

 

 

1

0

 

dx)

 

 

a1 , что,

в свою очередь, выполняется при

 

 

a

 

,

( a0 )

 

 

 

x0

 

 

 

0 x l.

 

 

 

 

 

 

 

(x)

 

 

 

 

 

 

 

0

 

 

 

Выбрав

таким

 

образом

 

и положив

(x) (x)a0(x),

q(x) (x)a2(x), приводим уравнение (1.1) к виду (1.2).

 

 

 

 

 

 

 

 

 

При этом важно, что функция

 

(x) определена на отрезке

[0,l],

непрерывно дифференцируема и удовлетворяет условию (x) 0

0, где

 

0

некоторая положительная постоянная. Уравнение (1.2) обладает рядом интересных свойств. Отметим некоторые из них.

Пусть y(x) и z(x) – решения следующих уравнений:

L[y] f (x), L[z] g(x).

Умножая первое из них на z(x), а второе – на y(x) и вычитая почленно полученные результаты, будем иметь

r(x)L[y(x)] y(x)L[z(

Так как

d z(x)dx

d

dx

 

 

 

d

 

 

dy(x)

 

 

 

d

 

 

x)] z(x)

 

 

 

p(x)

 

 

 

y(x)

 

 

 

p(x)

 

 

 

dx

 

 

 

 

 

 

dx

 

 

 

 

 

 

dx

z(x) f (x) y(x)g(x).

 

 

 

 

 

 

 

dy(x)

 

 

d

 

 

dz(x)

 

 

p(x)

 

 

 

y(x)

 

 

p(x)

 

 

 

 

 

dx

 

 

 

 

 

dx

 

 

 

 

 

 

 

dx

 

 

 

 

 

 

 

 

 

dy(x)

 

 

 

 

dz(x)

 

 

 

 

p(x) z(x)

 

 

 

 

y(x)

 

 

 

,

 

 

 

 

 

 

dx

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dz(x) dx

(1.3)

95

то из равенства (1.3) следует, что

 

 

 

 

d

 

dy(x)

 

dz(x)

 

 

 

p(x) z(x)

 

y(x)

 

 

z(x) f (x) y(x)g(x).

 

 

dx

dx

 

dx

 

 

 

 

Это соотношение называется тождеством Лагранжа. Его часто переписывают в операторной форме

z(x) L y(x) y(x) L z(x) z(x) f (x) y(x) g(x).

(1.4).

Интегрируя тождество (1.4), получаем формулу Грина:

l

 

 

dy(x)

(z(x)L y(x) y(x)L z(x) )dx p(x) z(x)

 

dx

0

 

 

l

z(x) f (x) y(x)g(x) dx.

0

y(x)dz(x) dx

x l

x 0

Из (1.4) следует, что если y(x) и z(x) – решения однородного уравнения

 

 

 

 

 

L[y] 0,

(1.5)

 

dy(x)

 

dz(x)

 

 

 

 

то p(x) z(x)

 

y(x)

 

 

c, где постоянная с не является произвольной,

dx

dx

 

 

 

 

 

 

 

а зависит от выбора решений

y(x) и z(x). Отсюда находим, что определитель

Вронского решений y(x) и z(x)

имеет вид

 

 

 

 

 

 

(y,z)

c

 

 

 

 

 

 

 

p(x) .

(1.6)

 

 

 

 

 

 

Замечание. Из соотношения (1.6) следует, что если известно одно решение

y1(x) уравнения (1.5),

то

любое другое

его

решение y(x)

удовлетворяет

уравнению (y (x),y)

c

 

. Это соотношение представляет собой линейное

 

 

1

p(x)

 

 

 

 

 

 

 

 

 

 

dy

 

dy1(x)

 

 

c

 

неоднородное уравнение первого порядка

y (x)

 

y

 

, зависящее

dx

 

 

 

 

 

 

 

1

 

dx

p(x)

от произвольной постоянной с. Его общее решение можно получить методом вариации произвольной постоянной. В итоге общее решение уравнения (1.5)

 

 

 

 

 

 

x

ds

 

 

 

можно получить в виде y(x) y1(x) c1

c

 

 

.

 

2

 

 

 

 

 

 

 

 

x0

p(s)y1

()s

 

§2. Краевая задача. Функция Грина

 

 

 

 

 

 

Рассмотрим следующую задачу. Требуется найти решение y(x)

уравнения

 

d

 

dy

 

 

 

на отрезке 0,l ,

 

L[y]

 

p(x)

 

q(x)y f (x), непрерывное

которое

 

 

 

dx

dx

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a y(0) a y (0) 0,

 

 

(2.1)

 

 

 

 

0

1

 

 

 

 

 

 

 

 

0 y(l) 1y (l) 0,

 

 

 

где a0,a1, 0, 1 – заданные постоянные, такие, что a02 a12 0 и 02 12 0.

96

Поставленная задача называется краевой задачей.

Если краевые условия неоднородны, т.е. имеют вид

 

 

 

 

,

a y(0) a y (0)

 

0

1

 

1

(2.2)

0 y(l) 1 y (l) 2,

где 1 и 2 – некоторые постоянные, то задачу можно свести к такой же задаче, но с однородными условиями (2.1) следующим образом. Сначала находим функцию u(x) такую, чтобы она удовлетворяла условиям (2.2). Обычно ее можно построить в виде полинома u(x) ax b. Затем в уравнении (1.2) делаем

замену y y1

u(x).

 

 

 

 

 

 

 

 

 

 

 

В результате относительно неизвестной y1

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

 

L[y]

d

 

p(x)

dy

 

q(x)y f (x),

f1(x)

f (x)

d

p(x)

du

 

q(x)u

(2.3)

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

dx

dx

 

 

 

dx

dx

 

 

с однородными граничными условиями (2.1). Поэтому в дальнейшем будем рассматривать только задачу (1.2), (2.1). Эта задача решается с помощью функции Грина, которая определяется следующим образом.

Определение. Функцией Грина будем называть функцию G(x,s), 0 x, s l , удовлетворяющую следующим условиям:

1.G(x,s) непрерывна на x и s при 0 x, s l .

2.G(x,s) как функция переменной x удовлетворяет однородному уравнению

 

d

dy

 

 

0 x l

 

 

 

p(x)

 

 

q(x)y 0,

(2.4)

 

 

 

 

dx

dx

 

 

 

 

 

при любом фиксированном s 0,l и условиям (2.1).

 

3. Первая производная Gx (x,s) имеет разрыв при

 

x s, величина которого

определяется соотношением

 

 

 

1

 

 

 

Gx (s 0,s) Gx

(s 0,s)

 

, 0 s l .

(2.5)

 

 

 

 

 

 

 

 

p(s)

 

Из определения функции Грина еще не следует ее существование для каждой краевой задачи (2.1), (2.4). Докажем, что краевая задача (2.1), (2.4) имеет функцию Грина, если эта задача имеет только тривиальное решение в классе дважды дифференцируемых функций.

2.1. Краевая задача для неоднородного уравнения. Перейдем к решению краевой задачи (1.2), (2.1). Напомним об основных предположениях, при которых рассматривается эта задача. Они состоят в следующем.

1.

Функции p(x), p (x), q(x)

и

f (x) непрерывны при 0 x l.

2.

Существует постоянная p0

(

p0 0) такая, что p(x) p0 при 0 x l.

3. Однородная краевая задача (2.1), (2.4) имеет только тривиальное решение в классе дважды дифференцируемых функций.

97

При выполнении этих условий задача (2.1), (2.4) однозначно разрешима. Если же эти условия не выполняются, то, как показано ниже, подобное утверждение относительно краевой задачи (2.1), (2.4) не всегда верно.

Теорема 1. При выполнении указанных выше предположений 1–3 функция

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y(x) G(x,s) f (s)ds

 

 

(2.6)

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

является решением краевой задачи (1.2), (2.1).

 

 

 

 

 

 

Из (2.6) находим, что

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y (x) G(x,s) f (s)ds.

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

Имея в виду представление функции Грина, запишем:

 

 

 

 

 

 

 

x

 

 

 

l

 

 

 

 

l

 

 

 

 

 

 

 

(x,s)f (s)ds Gx (x,s) f (s)ds Gx(x,s) f (s)ds.

 

 

 

y (x) Gx

 

 

 

 

 

0

 

 

 

x

 

 

 

 

0

 

 

 

Отсюда следует, что

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

l

 

 

 

 

 

(x,x 0) f (x) Gx(x,x 0) f (x)

 

 

 

 

 

 

 

 

 

 

 

y (x) Gxx(x,s) f (s)ds Gxx(x,s) f (s)ds Gx

 

 

0

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Gxx (x,s) f (s)ds Gx (x,x 0) Gx (x,x 0) f (x).

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

Так как, по определению, функция Грина удовлетворяет условию (2.5), то

 

 

 

 

 

 

 

l

 

 

 

 

f (x)

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y (x)

Gxx (x,s)f (s)ds

p(x)

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

L y(x)

 

 

 

 

 

 

l

p(x)G''xx (x,s)

 

 

Следовательно,

 

 

 

 

 

 

p(x)y

(x) p

(x)y

(x) q(x)

 

 

 

 

q(x)G(x,s) f (s)ds

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

f (x)

.

Так

как

выражение

в

квадратных

 

p (x)Gx(x,s)

 

 

 

 

 

 

 

 

 

 

скобках тождественно равно нулю (по определению G(x,s)

– решение

однородного уравнения L[y] 0), то отсюда находим, что

 

 

l

LG(x,s) f (s)ds f (x),

0

что и требовалось доказать.

Рассмотрим вопрос о единственности решения краевой задачи (1.2), (2.1).

Теорема 2.

Если однородная краевая

задача (1.5), (2.1) имеет только

тривиальное

решение

y(x) 0

в

классе

дважды

непрерывно

дифференцируемых функций, то задача (1.2), (2.1) имеет единственное решение.

Доказательство получается методом от противного. Если предположить,

что y1(x) и y2 (x),

y1(x) y2 (x), – два решения задачи, то

функция

y(x) y1(x) y2 (x), не

равная тождественно нулю, является

решением

однородной задачи (1.5), (2.1), а это противоречит условию теоремы.

 

98

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

Теорема 3. Необходимым условием разрешимости краевой задачи (1.2), (2.1) является условие

l

 

(x) f (x)dx 0,

(2.7)

0

 

где (x) – произвольное решение соответствующей однородной задачи (1.5), (2.1).

Доказательство. Пусть y(x) – решение неоднородной краевой задачи (1.2), (2.1), а (x) – решение однородной задачи (1.5), (2.1). Применим формулу Грина к этим функциям:

l

 

 

dy(x)

 

d (x)

( (x)L[y)(x)] y(x)L[ (x)])dx

p(x) (x)

 

y(x)

 

dx

dx

0

 

 

 

l

x l

x 0

(x) f (x) dx.

(2.8)

0

 

Докажем, что выражение, стоящее в середине цепочки равенств, равно нулю. В самом деле, из того, что (x) и y(x) удовлетворяют первому

граничному условию из (2.1), имеем равенства

a0 (0) a1 (0) 0; a0 y(0) a1y (0) 0,

в которых a0 и a1 не могут быть одновременно равны нулю (см. (2.1)). Для определенности предположим, что a1 0. Тогда

 

dy(x)

d (x)

 

 

1

 

a (0)y(0) a

y(0) (0) 0.

(x)

 

y(x)

 

 

 

 

 

 

 

 

 

 

dx

 

 

a

 

dx

 

x 0

 

0

 

0

 

Аналогично доказывается, что

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dy(x)

 

 

 

 

d (x)

0.

 

 

 

(x)

 

 

 

y(x)

 

 

 

 

 

 

dx

 

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

x l

 

 

Следовательно, равенства (2.8) принимают вид

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

l

 

 

 

( (x)L[y(x)] y(x)L[ (x)])dx f (x) (x)dx 0,

 

0

 

 

 

 

 

 

 

 

 

 

 

0

 

 

что и требовалось доказать.

 

 

 

Пример. Рассмотрим краевую задачу

 

 

 

 

y y f (x); 0 x ; y(0) y( ) 0.

(2.9)

Соответствующее однородное уравнение

y y 0

имеет общее

решение

y c1 cosx c2 sin x

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

существует

однопараметрическое

семейство решений

y1(x,c) csin x, удовлетворяющих краевым

условиям

99

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