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

Lesn3_Integraly

.pdf
Скачиваний:
9
Добавлен:
11.05.2015
Размер:
3.63 Mб
Скачать

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

n

 

 

 

c*jY j является решением системы уравнений (3). Это вытека-

j 1

 

 

 

ет из теоремы 1.

 

 

2. Пусть

U(x) некоторое решение системы уравнений (3).

Доказывается,

что найдутся числа c*,c*,...,c* R, для которых

 

1 2

n

 

 

 

n

 

на интервале (a; b) выполняется равенство U ( x) c*jY j .

j 1

Построение фундаментальной системы решений рассмотрим для систем уравнений частного вида.

2.Системы линейных однородных уравнений

спостоянными коэффициентами

Вместо функциональной матрицы P(x) будем рассматривать числовую матрицу

a

a

11

12

a21

a22

A

 

 

 

... ...

 

an2

an1

...

a

 

 

 

1n

 

...

a2n

(5)

...

...

.

 

 

...

 

 

 

ann

 

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

Y A Y

(6)

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

удем искать решение системы уравнений (6) в виде

x ekx

 

x

 

 

 

1

 

 

1

 

 

x

 

ekx

x

 

ekx x ekx ,

Y

 

2

 

 

2

 

 

........

...

 

 

 

 

 

 

 

 

 

 

xnekx

xn

 

140

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

где x1, x2 ,...,xn и k – неизвестные числа.

Продифференцируем вектор-функцию: Y kx ekx . Подставим вектор-функции Y и Y в систему уравнений (6):

kx ekx Ax ekx . Отсюда следует равенство

 

 

Ax kx .

(7)

Равенство (7) означает, что искомый вектор

x является

собственным вектором матрицы A, а искомое число

k - соответ-

ствующим x собственным числом этой матрицы.

 

Собственные числа находятся как корни характеристиче-

ского уравнения

 

 

A kE

 

0 ,

(8)

 

 

 

 

 

 

 

которое является алгебраическим уравнением порядка n. Собственные векторы x находятся как решения системы

линейных однородных алгебраических уравнений

A kE x

 

.

 

0

(9)

зависимости от корней характеристического уравнения рассмотрим три случая.

1. Уравнение (8) имеет n различных действительных кор-

ней k1, k2 ,...,kn .

Этим собственным числам соответствуют n линейно независимых cобственных векторов x1, x2 ,...,xn . Данным векторам

соответствует n различных решений Y1,Y2 ,...,Yn . Они образуют фундаментальную систему решений системы уравнений (3).

Пример 1.

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

 

2 y1

y2 ,

 

y1

(*)

y

3y

4 y

 

2

.

2

1

 

 

Имеем систему линейных однородных уравнений с постоянными коэффициентами. Решения ФСР будем искать в виде Y x ekx .

141

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

 

 

 

 

 

1.

 

 

 

2

1

 

 

 

Запишем матрицу системы уравнений A

 

 

.

 

 

 

 

 

 

 

3

4

 

 

 

 

 

 

 

 

 

 

 

2.

Составим характеристическое уравнение

 

2 k

1

 

0 .

 

 

 

 

3

 

4 k

 

3.

Решим его:

 

 

 

 

 

 

 

 

 

(2 k)(4 k) 3 0;

 

 

 

 

 

 

k 2 6k 5 0 ;

 

 

 

 

 

 

k1 1,

k2 5 . Корни действительны и различны.

 

 

4.

Запишем систему для нахождения собственных векторов x :

(2 k)x1 x2 0,3x1 (4 k)x2 0.

5. Найдем собственные векторы, отвечающие собственным зна-

чениям k1 1,

k2

5 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x x 0,

 

 

 

x

 

x ,

 

 

 

 

 

1

k 1 :

 

1

2

 

x x 0 ,

 

 

2

1

,

 

x

 

 

.

1

 

 

 

 

 

 

1

 

2

 

 

t.

 

 

 

1

 

 

 

 

 

 

3x1 3x2 0.

 

 

 

x1

 

 

 

 

 

 

1

 

 

5 :

3x x 0,

3x x 0 ,

x

 

3x ,

 

 

 

 

 

1

 

k

 

3x

1

 

2

 

 

2

1

,

x

 

 

.

 

2

 

x

2

0.

 

1

2

x t.

 

 

 

2

 

3

 

 

 

 

 

 

1

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

6. Запишем решения фундаментальной системы:

 

 

 

 

 

 

 

 

 

 

Y

 

1

 

1

e5x .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ex , Y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

1

2

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

7. Согласно равенству (4) общее решение системы (*) имеет вид: Y C1Y1 C2Y2 . Запишем его в развернутом матричном виде:

y

 

C

 

1

1

 

 

ex

 

 

1

 

y2

 

 

 

1

 

 

1

 

 

C ex C

2

e5x

 

C

 

 

1

 

 

.

 

 

e5x

C ex 3C

 

2

3

 

 

e5x

 

 

 

 

 

1

 

2

 

 

Обычная запись общего решения системы (*):

y C ex C

e5x ,

1

1

2

 

 

y

2

C ex 3C

e5x.

 

 

1

2

 

Рассмотрим теперь другой случай.

142

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

2. Характеристическое уравнение (8) имеет n различных корней, как действительных k, так и комплексных i .

В этом случае каждому действительному корню k

по-

прежнему ставится в соответствие решение Y x ekx .

Оно

включается в фундаментальную систему решений системы уравнений (3).

С парой комплексно сопряженных корней i поступа-

ем следующим образом.

Подставив i в систему (9) для нахождения координат

собственных векторов, получим систему с комплексными коэффициентами. Ее решение тоже будет комплексным:

 

u1 v1i

 

 

u1

 

 

v1

 

 

 

 

 

 

 

 

 

 

 

 

 

x

u2

v2i

 

u2

 

 

v2

 

(10)

 

 

 

 

 

 

i u vi .

1

..........

 

 

...

 

 

...

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

un vni

 

un

 

vn

 

Тогда числу i будет соответствовать собственный

вектор x2 u vi .

Оказывается, в этом случае в фундаментальную систему решений системы уравнений (3) можно включить следующие две вектор-функции:

Y e x

(

 

 

 

sin x v cos x),

 

u

 

1

 

 

 

 

 

 

 

 

 

(11)

Y e x (

 

 

 

cos x v sin x).

u

 

2

 

 

 

 

 

 

 

 

 

 

Рассмотрев таким образом все корни характеристического

уравнения (8), получим

 

 

 

фундаментальную систему решений

системы уравнений (3).

 

 

 

 

 

 

 

 

 

Пример 2.

 

 

 

 

 

 

 

 

 

 

Решим систему уравнений

 

 

 

 

 

 

 

 

y1 5y2 ,

 

y1

 

(*)

y

 

 

2 y

y

 

.

 

 

2

 

 

2

1

 

 

 

1. Матрица системы

 

 

 

 

 

1

5

 

 

A

 

.

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

2

 

 

 

143

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

 

 

2.

Характеристическое уравнение

 

1 k

5

 

0 .

 

 

 

2

1 k

 

3.

Решаем его:

 

 

 

 

 

 

k 2 9 0 ;

 

 

 

 

 

 

k1,2 3i . Корни комплексны и различны.

 

4. Запишем систему для нахождения собственных векторов:

 

 

(1 k)x1 5x2 0,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2x1 ( 1 k)x2

0.

 

 

 

 

 

5. Найдем собственный вектор, отвечающий собственному зна-

чению k1 3i .

 

 

 

 

 

 

(1 3i)x1 5x2 0,

 

 

x

2

1 (1 3i)x ,

(1 3i)x1 5x2 0 ,

 

 

5

1

 

( 1 3i)x2 0.

 

 

 

 

2x1

 

 

x

 

t.

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

5 0i

 

Взяв x

 

5 , получим собственный вектор

x

 

 

.

 

1

 

 

1

 

 

 

 

 

 

 

 

 

1 3i

 

Представим его так: x

 

5

 

0

 

 

 

 

 

i u vi .

1

 

 

 

 

3

 

 

1

 

 

 

6. Согласно равенству (11) запишем решения фундаментальной системы:

 

5

 

 

0

 

 

;

Y

5

 

 

 

0

sin 3x .

Y

sin 3x

cos 3x

 

 

cos 3x

 

 

 

1

1

 

 

3

 

 

 

2

1

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

7. Общее решение имеет вид

 

Y C1Y1 C2Y2 .

 

 

 

 

Обычная запись общего решения системы (*):

 

 

 

 

 

 

y

C 5sin3x C

2

5cos 3x,

 

 

 

 

 

 

 

 

 

1

1

 

 

 

 

 

 

 

 

 

 

 

 

 

y2

C1(sin 3x 3cos 3x) C2 (cos 3x 3sin3x).

 

 

 

3. Среди корней характеристического уравнения (8) есть кратные корни. В этом случае построение фундаментальной системы решений является более сложным. Перейдем к исследованию систем линейных неоднородных

уравнений.

144

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

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

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

Y P(x) Y Q(x) .

(1)

тесно связано с решением соответствующей системы линейных однородных дифференциальных уравнений

Y P(x) Y .

(2)

Теорема 1. (О структуре общего решения). Пусть Y1,Y2 ,...,Yn -

фундаментальная система решений системы уравнений

(2) и Yчн - частное решение системы уравнений (1). Тогда общее решение системы уравнений (1) имеет вид:

n

 

Yон Yчн C jY j

(3)

j 1

Доказательство опускаем.

Для решения системы уравнений (1) используются метод вариации произвольных постоянных и метод подбора частного решения.

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

Алгоритм решения системы линейных неоднородных уравнений

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

1. Находим общее решение системы уравнений (2):

n

Yоо C jY j .

j 1

145

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

2. Записываем форму общего решения системы (1):

n

 

Yон Ci (x) Yi .

(4)

i1

3.Записываем систему уравнений для нахождения производных неизвестных функций C j ( x) :

n

 

C j ( x) Y j Q( x) .

(5)

j1

4.Находим производные функций C j ( x) .

5.Интегрированием находим сами функции C j ( x) .

6.Подставляя найденные функции C j ( x) в равенство (4), записываем общее решение системы уравнений (1).

Пример 1.

Решим систему уравнений

y

2 y y

2

3e5x ,

 

1

1

 

(*)

 

 

3y 4 y

 

y

2

e5x.

 

2

1

 

 

1. Запишем общее решение однородной системы:

 

1

 

 

1

 

 

 

 

Y C

 

ex C

2

 

e5x . (См. пример 1 из §3).

 

1

 

 

 

3

 

 

 

 

 

1

 

 

 

 

 

2. Общее решение неоднородной системы будем искать в виде:

 

 

 

1

 

 

 

 

1

 

(**)

Y C (x)

ex

C (x)

e5x .

1

 

 

 

 

 

2

3

 

 

 

 

 

1

 

 

 

 

 

3. Запишем систему уравнений (5) для нахождения производных

неизвестных функций

 

C1(x),

 

C2 (x) :

 

 

1

 

 

1

 

2

1

 

3

 

 

 

 

 

 

 

 

C (x)

e x C

(x)

e5x

 

e5x .

 

 

 

1

 

 

 

3

1

В развернутом виде:

 

 

 

 

 

 

C

(x)ex

C (x)e5x 3e5x ,

 

 

 

1

 

 

2

 

 

 

 

 

 

 

(x)ex

3C

(x)e5x

e5x .

 

 

C

 

 

 

1

 

 

2

 

 

 

 

 

146

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

4. Складывая уравнения,

получаем

C (x) 1 . Подставляя эту

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

производную в первое уравнение, находим

C (x) 2e4x .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

5. Интегрируя производные, находим сами функции:

C (x)

1

 

e4x C

;

C (x) x C .

 

 

 

 

 

 

1

 

 

 

 

 

2

 

 

1

2

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6. Подставляя найденные функции в равенство (**), получаем

общее решение исходной системы уравнений:

 

 

y1

 

 

1

 

 

 

1

 

 

1

 

 

 

 

(

 

e4x C )

ex (x C ) e5x .

 

 

2

 

 

 

 

 

 

1

 

 

2

3

 

y2

 

 

 

 

 

 

 

 

1

 

 

В обычной записи:

 

 

 

 

 

 

 

y

1

e5x xe5x C ex

C e5x ,

 

 

 

 

 

1

 

2

 

 

 

 

 

 

 

 

 

1

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

2

 

1

e5x 3xe5x

C ex 3C e5x .

 

 

 

 

 

2

 

 

 

 

 

 

 

 

1

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

то касается метода подбора частного решения по правой части специального вида, ограничимся одним замечанием.

Для нахождения частного решения системы уравнений с постоянными коэффициентами

Y A Y Q(x)

(6)

методом подбора используются утверждения, аналогичные теоремам 2 - 4 (§4 предыдущей главы) для одного линейного неоднородного уравнения с постоянными коэффициентами.

На этом мы закончим исследование обыкновенных дифференциальных уравнений и систем дифференциальных уравнений. Обратимся к вопросам интегрирования функций нескольких переменных.

147

Глава VI. Интегрирование функций нескольких переменных

Лекция 17

Глава VI

Интегрирование функций нескольких переменных

До сих пор мы занимались интегрированием функций одной переменной. Теперь будем исследовать функции нескольких переменных. Начнем с функций двух переменных.

§1. Понятие двойного интеграла

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

1. Задача об объеме пространственной области

Рассмотрим пространственную область, которая определяется графиком функции аналогично криволинейной трапеции.

 

 

z

 

 

 

 

 

Пусть

функция

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z = f(x, y)

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

на на замкнутой

 

 

 

 

 

 

 

 

 

 

 

 

 

S

 

 

ограниченной

 

 

 

 

 

 

 

 

 

 

 

 

 

 

области D R 2.

 

 

 

 

 

 

 

 

Ее

графиком

 

 

 

O

 

 

 

 

является

неко-

 

 

 

 

D

 

y

торая

 

поверх-

 

 

 

 

 

 

 

x

 

 

 

 

 

 

ность

S

в трех-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

мерном

геомет-

рическом пространстве R 3, определяемая уравнением z = f(x, y).

148

§1. Понятие двойного интеграла.

Определение 1. Пусть пространственная область V ограничена плоскостью z = 0, поверхностью z = f(x, y) и цилиндрической поверхностью, направляющая которой совпадает с границей области D, а образующие параллельны оси Oz. Тогда область V называется цилиндроидом функции f(x, y) на области D.

Задача.

На области D R 2 задана положительная непрерывная функция f(x, y). Вычислить объем V ее цилиндроида с основанием D.

Решение. Предположим сначала, что функция f(x, y) является постоянной на области D: f(x, y) = c. Тогда задача решается достаточно просто: V = c SD.

Пусть теперь функция f(x, y) не является постоянной на области D. Если бы область D была достаточно малой, то в силу непрерывности функции f(x, y) ее значения на этой области мало отличались друг от друга. Учитывая это, постараемся решение задачи свести к решению предыдущей. Проделаем следующие операции.

1. Разобьем область D на n малых областей D1, D2,…, Dn, которые не будут иметь общих внутренних точек. Площадь ма-

лой области Di

обозначим через

si.

 

 

 

 

 

 

Цилиндроид с основанием

D разобьется на n малых ци-

линдроидов

с

ос-

 

 

 

 

 

 

 

нованиями

 

Di.

 

z

z = f(x, y)

 

 

 

f (xi,yi)

 

 

Объем

V цилин-

 

 

 

 

 

 

 

 

 

 

дроида

будет

ра-

 

 

 

 

 

 

 

вен

сумме

объе-

 

 

 

 

 

 

 

мов Vi всех малых

 

 

 

 

 

 

 

 

 

 

 

 

 

 

цилиндроидов.

 

 

O

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2.

Рассмот-

 

 

 

 

 

 

y

рим

малый

 

ци-

 

 

 

 

 

 

 

x

 

Mi(xi,yi)

 

 

 

 

 

 

линдроид с

осно-

 

 

Di

 

 

 

 

 

 

 

ванием

Di

 

(на

 

 

 

 

 

 

 

рисунке

Di

-

прямоугольник).

Возьмем произвольную точку

149

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