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

Методы математической физики. Интегральные уравнения (90

.pdf
Скачиваний:
16
Добавлен:
15.11.2022
Размер:
362.61 Кб
Скачать

Федеральное агентство по образованию Российской Федерации

Государственное образовательное учреждение высшего профессионального образования

Ивановский государственный химико-технологический университет

МЕТОДЫ МАТЕМАТИЧЕСКОЙ ФИЗИКИ Интегральные уравнения

Методические указания

Составитель: Г.А. Зуева

Иваново 2006

Составитель: Г.А. Зуева

УДК 613.19

Методы математической физики. Интегральные уравнения: Методические указания / Иван. гос. хим.-технол. ун-т; Сост. Г.А. Зуева. - Иваново, 2006. – 32 с.

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

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

Рецензент доктор химических наук, профессор В.В. Рыбкин

(ГОУ ВПО «Ивановский государственный химико-технологический университет»)

2

Вопросы и задачи для самопроверки

1.Что такое неопределенный интеграл?

2.Напишите формулу интегрирования по частям.

3.Напишите формулу Ньютона-Лейбница.

4.В чем состоит геометрический смысл определенного интеграла?

5.Чему равна производная определенного интеграла по переменному

верхнему пределу?

6.Что такое несобственный интеграл?

7.Как определяются пространства R2 , R3 , Rn ?

11

8.Найдите интеграл ∫∫x2 y2dxdy .

00

9.Найдите общее решение для уравнения

1)у-6у+9у=е;

2)у+3+2у=cos2x+2sin2x.

10.Найдите интегралы вида:

 

1) x n dx;

2) x cos(nx)dx ;

3) xsin(nx)dx ;

4) cos 2 (nx)dx;

5) sin2 (nx)dx;

6) cos( nx ) sin( mx )dx ;

те же интегралы найдите в пределах

 

8) от 0 до π,

9) от -π до π,

10) от 0 до l.

l

11) найдите cos( πn / l )dx , n=1,2,…

0

11.Что такое скалярное произведение? Какими свойствами оно обладает? Какие векторы называются ортогональными?

12.Какой может быть область сходимости степенного ряда?

13.Найдите область сходимости ряда

xn

n2

.

n =1

+1

14.Сформулируйте задачу Коши для дифференциального уравнения первого порядка.

15.Сформулируйте теорему единственности решения задачи Коши для дифференциального уравнения первого порядка.

3

1. Интегральные уравнения. Основные понятия и определения. Классификация уравнений

Определение: Интегральными уравнениями называются уравнения, содер-

жащие неизвестную функцию под знаком интеграла.

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

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

Пример линейного уравнения

1

e

x

 

ϕ(x) e xyϕ( y)dy = x

 

, 0 x 1.

x

0

 

 

 

 

Пример нелинейного уравнения

1xyϕ( y)

ϕ(x) 01 +ϕ2 ( y)dy = f (x), 0 x 1.

Решить интегральное уравнение – значит найти такую функцию, которая обращает данное уравнение в верное тождество.

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

Уравнения первого и второго рода с постоянными пределами интегрирования называются уравнениями Фредгольма, а уравнения с переменным верхним пределом называются уравнениями Вольтерра.

Таким образом, линейное интегральное уравнение Фредгольма первого рода имеет вид:

b

K (x, y)ϕ( y)dy = f (x), a x b.

a

Линейное интегральное уравнение Фредгольма второго рода имеет вид

b

ϕ(x) K(x, y)ϕ( y)dy = f (x), a x b.

a

Линейное интегральное уравнение Вольтерра первого рода имеет вид

x

K (x, y)ϕ( y)dy = f (x), a x.

a

Линейное интегральное уравнение Вольтерра второго рода имеет вид

4

x

ϕ(x) K(x, y)ϕ( y)dy = f (x), a x.

a

Здесь

ϕ(x) – неизвестная искомая функция;

f(x) – заданная непрерывная функция, называемая свободным членом интегрального уравнения;

K(x,y) – заданная непрерывная функция, называемая ядром интегрального уравнения.

Уравнения второго рода иногда записывают с параметром λ так:

b

 

 

ϕ(x) λK (x, y)ϕ( y)dy = f (x),

a x b.

(1)

a

 

 

или

 

 

x

 

 

ϕ(x) λK(x, y)ϕ( y)dy = f (x),

a x.

(2)

a

Тогда уравнения (1) или (2) представляют собой не одно уравнение, а семейство уравнений, зависящее от числового параметра λ.

Если f(x) = 0, то интегральное уравнение называется однородным, в противном случае оно называется неоднородным.

Уравнения Вольтерра более просты, чем уравнения Фредгольма; уравнения второго рода более просты, чем уравнения первого рода. Уравнение Вольтерра можно при некоторых ограничениях рассматривать как частный случай уравнения Фредгольма.

2. Метод последовательных приближений решения интегральных уравнений

2.1. Построение решения уравнения Фредгольма второго рода при малых значениях параметра методом последовательных приближений.

Будем рассматривать интегральное уравнение Фредгольма второго рода

b

 

ϕ(x) λK (x, y)ϕ( y)dy = f (x), a x b,

(3)

a

где ядро K(x,y) и правая часть f(x) – заданные непрерывные функции переменных x, y, ϕ(x) - искомая функция, λ - параметр.

Вслучае, когда параметр λ мал, т.е. удовлетворяет условию

λ< M1 ,

где М – положительное число, такое, что

5

b b

 

 

 

 

1

 

b

 

2

 

2

 

 

 

K (x, y)

 

dxdy

 

M или

 

K(x, y)

 

dy M , a x b ,

 

 

 

 

∫∫

 

 

 

 

 

 

 

a a

 

 

 

 

 

 

 

 

a

решение ϕ(x) уравнения (3) существует и его можно построить методом по-

следовательных приближений.

Суть метода: функция ϕ(x) ищется в виде предела последовательности

ϕ0 (x) = f (x),

b

 

ϕn (x) = f (x) +λK(x, y)ϕn 1( y)dy, n =1,2,...

(4)

a

т.е.

ϕ(x) = lim ϕn (x)

n→∞

и других решений уравнение (3) не имеет.

Методом последовательных приближений найти решение уравнения Фредгольма 2-го рода, предварительно убедившись, что выполнено условие

 

 

 

 

b b

 

 

 

1

 

 

 

 

1

 

 

2

 

2

 

 

 

 

 

 

∫∫

 

K(x, y)

 

dxdy

 

 

 

 

 

 

 

 

 

λ

<

M

=

 

 

 

 

 

:

(5)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a a

 

 

 

 

 

 

 

b

1. ϕ(x) xyϕ( y)dy = 2x.

 

 

a

 

 

Здесь f(x) = 2х, K(x,y) = xy, λ= 1. Имеем

Р е ш е н и е.

b b

 

 

 

2

 

 

 

 

 

 

 

b b

 

 

 

 

 

 

K (x, y)

 

dxdy

 

 

= M

2

 

 

 

2

y

2

dxdy =1/ 9.

 

 

 

 

 

∫∫

 

 

 

 

 

 

 

= ∫∫x

 

 

a a

 

 

 

 

 

 

 

 

 

 

 

a a

 

 

 

 

Значит M = 1/3,

 

λ

 

=1 <

 

1

= 3. Следовательно, условие (5) выполнено.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1/ 3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Построим последовательность приближенных решений (4)

ϕ0 (x) = 2x,

 

 

b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ϕn (x) = 2x +1 xyϕn 1( y)dy, n =1,2,...

т.е.

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

ϕ1(x) = 2x +1 xy2 ydy = 2x + 2x

;

 

3

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6

 

 

b

 

 

 

 

 

2

 

 

 

 

 

 

 

1

 

1

 

 

 

 

 

 

 

ϕ2

(x) = 2x +1 xy(2 y + y

)dy = 2x + 2x(

+

);

 

 

 

 

 

3

3

9

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

 

 

 

 

 

 

2

 

 

2

 

 

 

 

 

 

1

 

1

 

1

 

ϕ

3

(x) = 2x +1 xy(2 y + y(

+

 

))dy = 2x + 2x(

+

+

);

 

9

 

3

9

27

 

a

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

1

 

1

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

ϕn

(x) = 2x + 2x(

+

+

+...

+

 

);

 

 

 

 

 

 

 

 

 

 

3

9

27

 

3n1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

… .

ϕ(x) = lim ϕn (x) .

n→∞

Т.к. сумма убывающей геометрической прогрессии

 

1

 

+

1

+

1

 

+... =

 

 

b1

=

 

 

1/ 3

 

=

1

, то

 

 

 

 

 

 

 

 

 

1 q

1

1/ 3

 

 

 

3 9 27

 

 

 

 

 

 

 

2

 

 

 

 

ϕ

(x) = 2x + 2x ·

1

 

= 3x.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2.

ϕ(x) +

 

 

 

cos2

yϕ( y)dy =1.

 

 

 

 

 

 

 

π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ответ: ϕ(x) =

.

 

 

 

 

 

 

 

 

 

 

 

 

 

π

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3.

y(x)

 

 

 

 

sin x y(t)tdt = 2 sin x.

 

 

 

 

2π

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ответ: y(x) = 4sin x .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

4.

y(x) π (1 x) sin 2πt

y(t)dt =

(1 x).

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ответ: y(x) = 1 x.

1 1

5. ϕ(x) ϕ( y)dy = sinπx. 2 0

Ответ: ϕ(x) =sin πx + π2 .

2.2. Построение решения уравнения Вольтерра второго рода методом последовательных приближений.

Будем рассматривать интегральное уравнение Вольтерра второго рода

7

x

 

ϕ(x) λK(x, y)ϕ( y)dy = f (x), a x.

(6)

a

 

Методом последовательных приближений функция ϕ(x) ищется в виде

 

предела последовательности

 

ϕ0 (x) = f (x),

 

x

 

ϕn (x) = f (x) +λK(x, y)ϕn 1( y)dy, n =1,2,...,

(7)

a

 

т.е.

 

ϕ(x) = lim ϕn (x) .

(8)

n→∞

 

Интегральное уравнение Вольтерра (6) при требовании непрерывности его ядра K (x, y) и правой части f (x) имеет единственное решение (8) для каждого

конечного значения параметра λ.

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

Обычно полагают ϕ0 (x) = f (x), однако это вовсе не обязательно: удачный

выбор ненулевого приближения часто позволяет ускорить сходимость последовательности ϕn (x) к точному решению.

Методом последовательных приближений найти решение уравнения Вольтерра 2-го рода:

x

6. ϕ(x) =1 (x y)ϕ( y)dy.

0

Ре ш е н и е. Положим

ϕ0 (x) =1.

Тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

x

2

 

 

 

x 2

 

 

 

 

 

 

 

 

 

 

 

 

 

ϕ

1

(x) =1

(x y)dy =1

 

 

 

 

=1

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

2

 

2!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

y 2

 

 

 

 

 

 

 

x 2

 

x 4

 

 

 

 

x2

 

 

 

x4

 

ϕ

2

(x) =1

(x y)(1

 

)dy =1

 

 

 

+

 

=

1

 

 

 

+

 

 

,

2!

 

2

 

24

2!

 

4!

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

y 2

 

 

 

y 4

 

 

 

 

 

x 2

 

 

 

x 4

 

 

x6

 

 

ϕ

3

(x) =1

(x y)(1

 

 

+

 

 

 

)dy =1

 

 

+

 

 

 

 

 

 

.

 

 

2!

 

4!

2!

 

4!

 

6!

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

… .

8

Для n-ого приближения

 

 

 

 

 

 

 

 

n x2n

 

 

x2k

 

 

 

x2

 

x

4

 

 

 

 

x6

 

 

 

 

 

 

n

k

 

ϕn (x) =

1

 

 

 

+

 

 

 

 

 

 

 

+... +... +

 

(1)

 

 

 

= (1)

 

 

 

.

2!

 

4!

 

6!

 

 

(2n)!

 

(2k)!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k=0

 

 

Тогда решение уравнения

2k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1)k

x

 

 

 

 

 

 

 

 

 

ϕ(x) =

lim ϕn (x) =

 

 

 

 

= cos x.

 

 

 

 

 

 

(2k)!

 

 

 

 

 

n→∞

 

 

 

 

 

k =0

 

 

 

 

 

 

 

 

 

 

 

 

П р о в е р к а: Подставим ϕ(x) = cos x в исходное уравнение

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ϕ(x) =1 (x y)ϕ( y)dy.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

cos x =1 x cos xdy +y cos xdy.

 

 

 

 

 

 

 

 

(9)

 

0

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Находим первый интеграл

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x cos xdy = x sin y

 

= x sin x .

 

 

 

 

 

 

 

 

 

 

 

(10)

0

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Находим второй интеграл

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

u = x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dv = cos ydy

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y cos xdy =

= y sin y

= x sin y +cos x 1.

 

 

 

(11)

0

 

 

v = sin y

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

du = dy

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Подставим (10), (11) в (9)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

cos x =1 x sin x + x sin x +cos1,

 

 

 

 

 

 

 

 

 

cos x = cos x - верно.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ответ: ϕ(x) = cos x.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

7. ϕ(x) =1 + ϕ( y)dy.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ответ: ϕ(x) = exp x.

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8. y(x) =1 x2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ xy(t)dt.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a) y0 (x) =1 x2 ;

 

 

 

 

 

 

 

 

 

б) y0 (x) =1.

 

 

 

 

Ответ: y(x) =1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

9. y(x) = x2 + x y(t)dt.

0

9

a) y

0

(x) =1;

б) y0 (x) =

x

2

+ x.

2

 

 

 

 

Ответ: y(x) = x.

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

10. ϕ(x) =1 + yϕ( y)dy,

ϕ0 (x) =1.

 

 

 

 

 

0

 

 

 

 

Ответ: ϕ(x) = exp(x2 / 2).

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

11. ϕ(x) = x (x y)ϕ( y)dy,

ϕ0 (x) = 0 .

 

 

 

0

 

 

 

 

Ответ: ϕ(x) = sin x .

 

 

x

 

 

 

 

 

 

 

 

 

 

 

12. ϕ(x) =1 (x y)ϕ( y)dy,

ϕ0 (x) = 0 .

 

 

 

0

 

 

 

 

Ответ: ϕ(x) = ch x .

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

13. ϕ(x) = 2 x + 2 x y ϕ( y)dy,

ϕ0 (x) = 0 .

 

 

 

0

 

 

 

 

Ответ: ϕ(x) = (2e) x .

 

 

 

 

 

 

 

3.Интегральные уравнения с вырожденными ядрами

3.1Решение уравнений Фредгольма 2-го рода с вырожденным ядром.

Ядро K(x,y) называется вырожденным, если оно имеет вид

 

n

 

 

K(x, y) = pi (x)qi ( y) .

 

(12)

i =1

 

 

Соответствующее интегральное уравнение (1)

 

b

n

 

 

 

pi (x)qi ( y)

 

(13)

ϕ(x) λ

ϕ( y)dy = f (x)

a i =1

 

 

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

Перепишем уравнение (13) в виде

n

ϕ(x) λ ci pi (x) = f (x), (14) i =1

где неизвестные ci определяются через искомое решение ϕ(x) равенствами

b

 

сi = qi ( y)ϕ( y)dy, i =1,2,..., n

(15)

a

10

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