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

Все в одном

.pdf
Скачиваний:
39
Добавлен:
13.02.2015
Размер:
23.06 Mб
Скачать

4. Машина Тьюринга

Опр: Алфавит - конечное непустое множество. Опр: Буквы – элементы алфавита.

Опр: Словом длины n над алфавитом А называется отображение U : [1..n]N A , где

U(i) i-я буква слова.

Опр: Равенство слов – это равенство отображений.

Опр: Пусть А-алф-т и ΛÏА. Тогда АÈ{Λ} алф-т с пустым символом.

Опр: Упорядоченной склейкой слов u и v длин n и m над А называется слово длины

n m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u(i) 1 i n

 

u v над А, задающееся по правилу: (u v)(i)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v(i 1) n 1 i m n

Склейка не коммутативна. Пример:

 

 

 

 

Пусть A:{|,*}, u=||*, v=** => u v || ***,

v u **|| * - не коммутативно.

Опр: Сдвигом на h Z называется отображение h

: Z Z ( h )(i) i h . Ясно, что

 

0

id

Z

,

h1

 

h2

 

h2

 

h1

 

h1 h2

, (

h

) 1

h

.

 

 

На ∞ ленте введем

0

1 отношение : u v h Z :

u h v .

 

Опр: Машина Тьюринга: ∞ лента, разбитая на ячейки (для записи слов в A { } ), СЗУ-устройство(счит-е запис-е устр-во), к-рое на каждом шаге работы машины умеет счит-ть содерж-ое одной ячейки и записывать вместо него букву из A { } . СЗУ

может сдвигаться на –1, 0, +1 ячейку. Работает в соответствии с программой. На такте известно, какую ячейку обозревает СЗУ, в каком состоянии находится машина. УУ обращается к соответствующей клетке Пт, а там всегда тройка (a, qi, s), и на место прочитанного знака записывается буква a, машина переходит в состояние qi и делает сдвиг s. Работа while q0.

МТ имеет 3 алфавита: внешний - A { } , сдвигов - S { 1,0, 1} и внутренний - Q {q0 ,..., qn } . q0 заключительное состояние, q1…qn рабочие состояния.

Опр: Программа для МТ-это отобр-ие (Пт): (( A { }) (Q \ {q0 })) ( A { }) Q S , т.е.

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

Пример: машина – сумматор (в унарной сист-е счисл-я). Сложить все | на ленте в формате:

 

^

^

|

 

|

 

+

|

|

|

^

^

 

 

 

 

 

 

q1

 

 

 

 

 

и стать на начало ответа.

 

 

 

 

 

 

 

 

 

 

 

 

 

Внешний алф-т машины A { } {, , }

 

 

 

 

 

Решение – программа:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q1

 

 

 

q2

 

 

q3

 

 

|

 

| q1 +1

^ q3 -1

| q3 -1

 

 

+

 

| q1 +1

 

 

 

 

 

 

 

 

 

^

 

^ q2

-1

 

 

 

 

^ q0

 

 

 

 

 

 

 

+1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

def
(T2 T1 )(u) T2 (T1 (u))

Пусть Т1 и Т2 – машины с общим внешним алфавитом A { } , но Q1 Q2 . Опр: Композицией машин Т1 и Т2 называется машина Т2◦Т1, действующая по

правилу: . Теорема: Композиция всегда .

Пусть Ti машина имеет программу:

 

 

qi 1

qi ni

 

 

A

Ti

 

Тогда машина Т2◦Т1

имеет программу:

 

 

 

q11

q1 ni

q2 1

q2 ni

A

ПT’1

ПT2

 

 

 

, где ПT’1 было получено из ПT1

заменой всех q10

на q21, т.е. Пт1

передает управление

Пт2.

 

 

 

 

 

Нач-м сост-ем Т2◦Т1 явл-ся начальное сост-е 1й машины, а заключит-ым – заключит-ое

состояние 2й машины.

 

 

 

 

 

Опр: МТ с полулентами бывают с левой и правой полулентой. Символ -

неподвижный ограничитель на таких лентах. За него нельзя выйти, нельзя затереть,

возможно лишь чтение. Ячейка после информации заполнена спецсимволом . Т.е.

рабочая зона МТ – диапазон [.. ] – для правой полуленты и [ ..] – для левой

полуленты:

 

 

 

 

 

 

 

 

 

правая полулента

нерабочая зона

информация

 

 

рабочая зона

 

 

 

 

 

 

 

 

 

 

 

 

 

левая полулента

информация

 

нерабочая зона

 

 

рабочая зона

 

 

 

 

 

 

 

 

Спецсимвол можно перемещать, т.е. рабочая зона может менять размеры.

def

Опр: Объединение Мт называется Мт (T1 T2)(uv) T1(u)T2(v)

Теорема: Объединение Мт всегда .

Построим T1 T2 в виде композиции: T1 T2 = T T2 T T1, где Т1– аналог Т1 на левой полуленте,

T2 - аналог T2 на правой полуленте,

T - машина(умеет переходить через ограничитель ), применимая к uv: начиная с 1й буквы u она останавливается на 1й букве v, оставляя на ленте слово uv.

T - отправляясь от 1й буквы v она останавливается на 1й букве u. 2

Очевидно, машины T и T сущ-т.

 

 

u1st

u

u

ulast

 

v1st

 

v

v

vlast

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

После Т1

 

 

 

 

 

 

u1st

u

u

u

u

ulast

 

v1st

 

v

v

vlast

 

 

 

⌂ ↔ ↔ ↔ ↔ ↔

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

После T

 

 

 

 

 

 

u1st

u

u

u

u

ulast

 

v1st

 

v

v

vlast

 

 

 

→ → → → → → → ⌂

 

 

 

 

 

 

 

 

 

 

 

После T2

 

 

 

 

 

 

u1st

u

u

u

u

ulast

 

v1st

 

v

v

v

v

v

vlast

 

 

 

 

 

 

 

 

⌂ ↔ ↔ ↔ ↔ ↔ ↔

 

 

 

 

 

 

После T

 

 

 

 

 

 

u1st

u

u

u

u

ulast

 

v1st

 

v

v

v

v

v

vlast

← ← ← ← ← ← ←

 

 

 

 

 

 

Пусть Т1 и Т2 – МТ с алф-ом A { } .

Опр: Машина – предикат Р – машина с внешнии алфавитом A {0,1}, применимая к

u A : P(u) {0,1} .

Опр: Разветвлением машин Т1 и Т2, управляемым предикатом Р называется

1

0

 

 

 

 

 

 

 

 

 

1

 

 

 

0

T2)(u)

 

T1 (u) if

 

P(u) 1

машина T1 p

T2, которая работает по правилу: (T1

p

 

 

 

P(u) 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T2 (u) if

 

Теорема: Разветвление машин всегда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

1

 

 

Построим разветвление в виде композиции

 

П

T1

 

 

(P E) K , где П

T1

- машина-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

1

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T2

 

 

 

 

 

 

 

 

T2

 

переключатель, P-предикат, E-тождественная машина, К1 – копирующая машина:

К1 (u)=uu, P E K1 - такое, что получаем либо 1u либо 0u.

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Программа П

T1

выглядит так:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T2 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q

q

q

q11

 

 

q1 ni

q2 1

 

q2 ni

 

 

 

 

 

A

 

 

 

 

 

 

 

 

Пт1

 

 

 

 

 

 

Пт2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

^ q11

^ q21

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+1

+1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

^ q2 +1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

^ q3 +1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

Опр: Итерация машин: пусть есть Мт Т и машина-предикат Р. Работа происходит так: применяем Т к слову пока P(u)=1. Иначе – стоп. Обозначается Т.

Теорема: Итерация всегда .

Составим T1 p 0Е. Чтобы получить Т, нужно преобразовать все клетки вида

(a, q0T, s) в (a, начальное состояние разветвления, s).

4

5. Графы

Опр: Граф – тройка g(X,U,f), где X – множество вершин (непусто), U – множество дуг графа, f :U X X – отображение инцидентности(т.е. правило, к-рое каждой дуге сопост-т упорядоченную пару вершин).

Если g(X,U,f) – граф, u U и f(u)=(x,y) => x – начало дуги, y – конец дуги. Если u v, f(u)= f(v), то это параллельные дуги.

Если u v, f(u)=(x,y), f(v)=(y,x), то это противоположные дуги.

Если f(u)=(x,x), то u – петля с вершиной x.

Опр: Графы g1(X1,U1,f1) и g2(X2,U2,f2) изоморфны (g1 g2), если

 

 

biect

 

biect

 

,

: X1

X 2 ,

:U1

U 2 :

( pi f2 )(u) ( pi f1 )(u) , причем

p1 (x, y) x

p2 (x, y) y .

 

 

Опр: Геометрический граф – граф, реализованный в геометрическим виде, у которого множество вершин есть точки пространства Rn, дуги – отрезки непрерывных кривых, заданных параметрически так, что начальное значение параметра соответствует началу дуги, конечное – концу.

Пример: в R

2

x x(t),

t [0,T ]. Вершина графа, явл-ся началом дуги, (x(0),y(0));

 

 

 

 

y y(t),

 

(x(0),y(0)) – конец дуги.

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

Геометрическая реализация графа – изоморфный данному геометрический граф. Теорема: Любой граф имеет правильную геометрическую реализацию в R3.

Приведем алгоритм построения правильной геометрической реализации в R3

словесно: фиксируем в R3 прямую l, вершине ставим в соответствие точку на l. дуге графа в правильной реализации ставим в соответствие полуплоскость, проходящую через l. Если u – не петля, т.е. f(u)=(x,y) x y , то в соответствующей полуплоскости изображаем полуокружность, имеющую диаметр отрезка [x,y]. Если v – петля в точке z, то изображаем единичную окружность, касающуюся точки z. Её обход произвольный. Так как полуплоскости не имеют общих точек кроме l, а на l

расположены только вершины, то реализация правильная.

Опр: Цепь длины n на g(X,U,f) - отображение :[1..n]N U : у дуги (i) одна вершина общая с (i 1) , а другая – с (i 1) .

Пусть есть граф g(X,U,f),x X вершина. Рассм-м множ-во

Sx= y X (x y) ( цепь_из _ x _ в _ y) .

Опр: Компонентой связности вершины x называется подграф, порожденный множеством Sx.

Опр: Число связности графа g - характеристика c(g) = количеству различных компонент связности графа.

Опр: Дуга u U называется мостом графа g(X,U,f), если c(g) c(g

) , где g

- граф без

 

G X ,U /{u}, f

 

 

.

{u}

{u}

 

дуги u,т.е. g

 

U /{u}

 

 

 

 

 

 

 

{u}

 

 

 

 

 

 

Теорема: (о мостах) Если u – мост g c(g) 1 c(g

)

 

 

 

 

 

 

{u}

 

 

 

1

Т.к. при удалении моста «разваливается» компонента связности, которой мост принадлежит, и остальные компоненты связности остаются без изменений, то нам

 

?

 

 

достаточно д-ть, что если u мост и c(g) 1 c(g

) 2 .

 

 

{u}

 

Допустим противное: G и мост u на G: c(g

) 3 . Тогда вершины x, y, z X такие,

 

{u}

 

 

что на g

они принадлежат разным компонентам связности. А на исход. графе g x,y,z

{u}

 

 

 

принадлежат одной компоненте связности (т.к. g связен) => цепь x y на g и цепь

x z на g. По лемме о простых цепей эти цепи м. считать простыми. Обозн: ηxy – простая цепь, соединяющая x и y и ηxz – простая цепь, соединяющая x и z.

Так как после удаления дуги u x, y, z оказались в разных компонентах связности => ηxy и ηxz разорвались => они обе содержали u. Так как обе цепи – простые => проходили через u по 1 разу. Возможен один из четырех случаев:

1)

 

 

 

 

 

 

 

 

 

2)

x◦

◦y

ηxy

x◦

 

◦y

ηxy

 

 

u

 

 

 

u

 

 

 

 

x◦

 

◦z

ηxz

x◦

 

◦z

ηxz

 

 

 

 

3)

 

 

 

 

 

 

 

 

 

4)

x◦

◦y

ηxy

x◦

 

◦y

ηxy

 

 

u

 

 

 

u

 

 

 

 

x◦

 

◦z

ηxz

x◦

 

◦z

ηxz

 

 

 

 

Часть цепей не пострадали при удалении u => они g

.

 

 

 

 

 

 

 

 

 

{u}

 

 

 

 

 

1) из

 

и

 

склеивается цепь, соединяющая y и z на графе

g

, что

 

 

 

 

 

 

 

 

 

{u}

 

противоречит тому, что y и z в разных компонентах связности.

2)аналогично

3)из непострадавших кусков цепей получим цепь, соединяющая y и x на графе что противоречит тому, что y и x разным компонентам связности

4)аналогично

Итак, во всех 4х случаях получили противоречие => если c(g) 1 c(g

) 2 .

{u}

 

g

,

{u}

 

Опр: Вершина x X точка сочленения графа g, если c(g) c(g

) , где g

– граф без

 

{x}

{ x}

 

вершины x, т.е. g

G X /{x},... - граф, полученный из g, удалением вершины x и

{x}

 

 

 

приходящих/уходящищ из нее дуг.

2

1. Уравнения в дифференциалах (УВД)

Они имеют вид P(x,y)dx+Q(x,y)dy=0, P и Q – определены и непрер-ны в некоторой обл. R2. Область задания дифура – там, где P и Q определены и непрерывны. Решение УВД – параметризованный образ гладкой кривой Г: x=x(t), y=y(t), t0<t<t1 (x,y непрер-но диф-мы на [t0,t1]), подстановка которых в уравнение обращает его в верное

тождество: P(x(t), y(t))

dx(t)

Q(x(t), y(t))

dy(t)

0 . Предполагается, что

dt

dt

 

 

 

 

 

 

 

 

2

2

(t) 0

t (t0 , t1 ) .

 

 

 

x

 

(t) y

 

 

 

Пример: xdx+ydy=0. Решение – единичная окружность. Параметризуем её: x(t):=sin(t),

* *

y(t):=cos(t), t [0,2π]. Проверка: sin t (sint) cost (cost) sin t cost cost sin t 0

Уравнение в полных дифференциалах

(1) M(x,y)dx+N(x,y)dy=0 – ур-ие в полн. диф-х, DÎR2, если такая непрер-но диф-мая функция F(x,y), что левая часть уравнения является её дифференциалом, т.е.

M(x,y)dx+N(x,y)dy=dF(x,y) или M (x, y) dFdx (x, y) N (x, y) dFdy (x, y) .

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

(x0,y0).

Теорема: (о ! решения ЗК в ПД)

Пусть в некоторой области D уравнение (1) является уравнением в полных

дифференциалах. Функции M, N определены и непрерывны в области D, (x0,y0) D – не является особой (т.е. M (x0 , y0 ) N (x0 , y0 ) 0 ) и выполнено условие:

(2) M 2 (x, y)dx N 2 (x, y)dy 0

в некоторой окрестности точки (x0,y0) ! решение (1), проходящее через (x0,y0).

По условию теоремы N (x0 , y0 ) 0

N отлично от 0 в некоторой окр-ти (x0 , y0 ) .

Перейдем к ЗК:

 

 

 

 

dy(x)

 

M (x, y)

 

 

 

 

(3,4)

 

 

. Из (3) (5)

M (x, y) N (x, y)

dy

0 . Это уравнение в полных

dx

N (x, y)

 

dx

 

y(x0 ) y0

 

 

 

 

 

 

 

 

дифференциалах M (x, y)

dF

(x, y)

N (x, y)

dF

(x, y) подставим в (5), получим

dx

dy

 

 

 

 

 

 

 

 

 

 

 

 

F

 

F dy

0

 

 

(F (x, y(x))) 0,

F (x, y(x)) C (6)

 

 

 

 

x

y dx

x

 

 

 

 

 

 

 

 

 

 

С учётом (4) C=F(x0,y0), а из (6) Ф(х,у(х))=F(x,y)-F(x0,y0), Ф(х,у(х))=0 (7). Применим теорему о неявной функции:

Ф(х,у) определена и непрерывна в некоторой окр-ти (x0,y0)

Ф(x0,y0) =F(x0,y0)-F(x0,y0)=0

Ф

 

F

M (x, y),

Ф

 

F

N (x, y) - определены и непрерывны в окр-ти (x0,y0).

x

 

x

 

y

 

y

 

Ф

N (x, y) N (x, y)

 

Ф

M (x, y) M (x, y)

y

 

 

 

 

 

x

 

1

Ф

(x0,y0)0

Ф

(x0,y0)0

y

 

x

 

окрестность х0, на которой определена функция у=у(х) такая, что Ф(х,у(х))=0х интервалу, у(х0)=у0, эта функция !, непрерывно дифференцируема.

Условия теоремы выполнены, значит теорема справедлива: интервал, содержащий (x0,y0), на котором определено решение (7) решение УвПД.

Аналогично M (x0 , y0 ) 0 => х=х(у).

 

 

 

Признаки уравнений в полных дифференциалах

Опр: Область D на плоскости называется односвязной если замкнутая, без самопересечений кривая D ограничивает область, целиком лежащую в D.

Теорема: (признак уравнения в полных дифференциалах).

Чтобы уравнение M(x,y)dx+N(x,y)dy=0 (1) было уравнением в полных дифференциалах в области D необходимо (а в случае односвязности области D и

достаточно), чтобы M N . Предполагается,чтоM, N(они

y x

непр.диф),

M ,

N

определены и непрерывны в D.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F

(x, y)

M (x, y)

 

 

 

2 F

 

(x, y)

M

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Если (1) – УвПД F(x,y) такая, что: x

 

 

y x

y

. Если

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F

(x, y)

N (x, y)

 

 

 

2 F

(x, y)

N

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

x y

x

 

смешанные производные непрерывны порядок дифференцирования не влияет:

 

N

 

2 F

 

 

2 F

 

 

M

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x y

 

y x

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Докажем лишь частный случай: D – прямоугольник. Пусть

M

N

, приведем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

x

 

 

алгоритм построения F: В обл. D выберем тчк (x0,y0). Так как

F

 

M (x, y)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

проинтегрировав:

 

F (x, y) F (x0 , y) M (s, y)ds (3), где F (x0 , y)

 

- неизвестная. Подставим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

получившееся в F

 

 

d

 

 

x

 

M (s, y)ds N (x, y) . Так как M N , то под

N (x, y) :

F (x0 , y)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

dy

 

 

x

y

 

 

 

 

 

 

 

 

 

y

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

N (s, y)

. Поэтому

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

 

 

 

 

 

x

 

 

 

 

Ньютон Лейбниц d

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F (x0 , y)

 

 

N (s, y)ds

 

 

 

F (x0 , y) N (x, y) N (x0 , y) . При приведении N(x,y)

 

dy

s

 

dy

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

получаем

 

F (x0 , y) N (x0 , y) F (x0 , y) F (x0 , y0 ) N (x0 , v)dv . Подставим в (3)

dy

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

x

 

 

 

F (x, y) F (x0 , y) M (s, y)ds , получим: F (x, y) F (x0 , y0 ) N (x0 , v)dv M (s, y)ds -

 

 

 

 

 

 

 

 

 

 

 

 

 

x0

 

 

 

 

 

 

 

 

 

y0

 

 

 

 

 

 

 

x0

 

 

определена с точностью до произвольной const.

2

Замечание

Теорема доказана для прямоугольной области D. Любую область можно приблизить прямоугольниками.

Если D – неодносвязна, то, поочередно рассматривая её односвязные подобласти, получим УвПД.

Интегрирующий множитель

Рассмотрим M(x,y)dx+N(x,y)dy=0 (1) – уравнение в дифференциалах. Пусть M N .

y x

Будем искать непрер.диф-ю функцию =(х,у) такую, что умножение на неё даст УвПД. То есть (х,у)M(x,y)dx+(х,у)N(x,y)dy=0 (2) – УвПД.

Опр: Интегрирующий множитель – это такая (х,у).

Опр: Уравнение интегрирующего множителя:

 

M

M

 

 

N

N

(3) при

 

 

 

y

 

y

 

x

 

x

 

( M )

 

( N ) .

 

 

 

 

 

 

 

 

y

 

x

 

 

 

 

 

 

 

 

Ур-ие (3) – диф.ур. в частных произв-х. Оно линейно.

Уравнение имеет решений. Если (x,у)-инт-ий множ-ль и с=const, то c (x,у) -инт-ий

множ-ль; если 1(x,у), 2(x,у)-инт-ие множ-ли и с12=const, то c1 1(x,у)+ c2 2(x,у) -инт- ий множ-ль.

Некоторые частные случаи:

а) = (x): 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d (x)

 

 

M

N

 

 

 

M 0 ,

N (

M

N )

 

1

 

 

 

 

y

 

 

x

Правая часть

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

y

x

 

 

 

 

x

 

(x) dx

 

N

 

 

 

 

зависит только от х (это необходимое условие ) => это уравнение с

 

 

 

разделяющимися переменными.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

б) =(у) аналогично.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример: 1. (x-y2)dx + 2xydy = 0. M(x,y)=x2-y2 , N(x,y)=2xy. M

N . Ищем (х):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

x

 

 

 

 

 

 

 

d (x)

 

M

N

2 y 2 y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

y

x

 

2

 

(x)

 

c

. Здесь =R2\{x=0}, т.к. х=0 – решение

 

 

 

 

 

 

 

 

 

 

 

 

(x) dx

 

N

2xy

 

x

 

 

 

 

 

x 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

уравнения. Умножим уравнение на :

(

1

 

 

y 2

 

)dx

2 y

dy 0 . Теперь M

N , это

 

 

x 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

x

 

 

 

y

 

 

x

уравнение в полных дифференциалах. Тогда F(x,y)=

2 y

dy (x)

y2

(x) ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

x

 

 

 

 

 

 

 

y 2

 

ln

 

 

(x) ln

x

C . Обобщенный интеграл исходного уравнения – это

 

x

.

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

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

Случай простых корней.

a0z(n)+a1z(n-1)+...+an-1z’+anz=0 (1). z(t): --> , ai , a00. Это однородное уравнение.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Обозначим

 

p

d

 

, p

2

d

 

,..., pn

d

 

. Тогда L(p):=a0pn+a1pn-1+...+an – диф-ый

 

dt

dt

2

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

многочлен=> (1) можем представить в виде L(p)z=0. (1’).

 

 

 

 

 

 

Лемма (1).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L(p)e t=L( )e t.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

p k e t

d

 

e t k e t L( p)e k an k p k e t

an k k e t L( )e t

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k 0

 

 

 

 

k 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Следствие:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

* - корень L( )=0-харак-ого ур-я

e *t-решение (1).

 

L(p)e *t=L(*)e *t=0

Лемма (2).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L( p)z 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z(t0 )

 

 

начальных данных (t0,z0,z0’,...,z0(n-1)) n+1 решение З.К. z'(t0 ) z'0

,! и

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

...

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z (n 1)

(t

0

) z (n 1)

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

определено на всей .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Сведем систему к НСДУ. Так как а0 0, то поделим на него: z(n)=b1z(n-1)+ b2z(n-

2)

+...+b z, где b

 

ak

. Заменим u :=z

(n-1)

. Система примет вид:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

k

 

 

a0

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u

u

2

,u

2

u

,...,u

n 1

u

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

. Решение , ! и определено на .

 

 

 

 

 

 

 

 

un

 

b1un

... bn u1

 

 

 

 

u (t

0

) z

, u

2

(t

0

) z'

0

,...,u

n

(t

0

) z (n 1)

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Теорема:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

все корни характеристического полинома L( ) простые общее решение

уравнения L(p)z=0 (1) имеет вид:

z(t) c e

1t

c

e

2t

... c

e

nt (6).

1..n

-различные корни

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

2

 

 

n

 

 

 

 

 

 

 

L( )=0, ci .

Реш-е явл-ся общим в том смысле, что решение уравнения (1) однозначно представимо в виде (6), причем const c1..n опред-ны ! образом.

1)Покажем, что фиксированного набора c1..n (6) – решение (1). k – корни L( )=0e k t -решение (1) – по следствию из Л.1 k=1..n. С др. стороны по свойству ЛСДУ л. комбинация e k t является решением.

2)Пусть (t) настоящее решение L(p)z=0. Будем считать, что (t) определено на

(Л.2). Докажем, что (t) можно представить в виде (6). (t) удовлетворяет З.К.:

 

L( p)z 0

 

.

 

0 , z'(0) '0 ,...,z (n 1)

(0) (n 1)

z(0)

0

1