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

ОТВЕТНИК ВСЕ В ОДНОМ. ФОРМАТ А5

.pdf
Скачиваний:
13
Добавлен:
08.02.2015
Размер:
981.9 Кб
Скачать

w1(x) = 0

Теперь найдем V (y).

 

 

dV

=

dV

 

 

dy

 

=

1

V 0(y)

 

 

dx

dy dx

 

 

 

 

 

"

 

 

 

 

d2V

=

d2V

 

dy

2 =

1

V 00(y)

 

 

 

 

 

Подставляем разложение для

 

в исходное уравнение, домножив

 

dx2

 

dy2

dx

"2

 

V

åãî íà ":

v000(y)+"v100(y)+(1+"y)v00 (y)+"(1+"y)v10 (y) "v0(y) "2v1(y)+: : : = 0

Выпишем коэффициенты при степенях ":

"0 : v000 + v00 = 0

"1 : v100 + yv00 + v10 v0 = 0

Решим сначала первое уравнение: v000 + v00 = 0

v0(y) = C2e y

Учтем, что v0(y) ! 0; v1(y) ! 0; y ! 1.

Из граничного условия U(0) = 0 следует, что v0(0) = w0(0) = 12 = C2, v1(0) = 0.

v0(y) = 12e y

Теперь решим второе уравнение. v00 (y) = 12 e y. v100 + v10 = 12e y(1 + y)

Ищем частное решение в виде v1p = (Ay2 + By)e y.

v10 p = (2Ay + B)e y e y(Ay2 + By) = e y(2Ay + B Ay2 By)

v100p = e y(2Ay + B Ay2 By) + e y(2A 2Ay B) =

= e y(2A 4Ay 2B + Ay2 + By)

41

e y(2A 4Ay 2B+Ay2+By)+e y(2Ay+B Ay2 By) = e y

1

(1+y)

 

2

 

 

 

 

 

 

 

 

( B + 2A =

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2Ay =

 

 

 

 

 

1 y;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A =

 

1

; B = 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

1

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v1p = e

 

1

 

 

 

 

y

 

 

+ y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v

= e y

 

 

y2

+

y

+

C

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Èç

v1

(0) = 0

следует,1

÷òî

 

 

4

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C3 = 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v1(y) = e y

 

1

y2 + y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

1

 

 

x

 

 

 

 

 

 

 

 

 

 

x

 

 

1 x2

 

 

 

x

 

 

 

 

U

 

 

(1 + x)

 

 

e "

 

+ "e "

 

 

 

 

+

 

 

 

 

 

 

 

 

2

2

 

4

"2

"

 

 

Ответ.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U

1

(1 + x) + x

1

e

x

+ e

x 1 x2

 

 

 

 

 

 

 

 

 

"

"

 

 

 

 

 

 

 

 

2

2

4

"

 

 

42

Задача 7.2.

Известно, что гамильтонова система

(

p0 = Hx; p; x 2 R3;

x0 = Hp; H(x; p) 2 C1(R6);

имеет два первых интеграла: f1 = p3x2, f2 = p22 + x3. Следует ли отсюда, что функция f3 = 2p2p3 + x2 также является первым интегралом для этой системы?

Решение.

Функция f является первым интегралом системы, если fH; fg =

0.

Вычислим ff1; f2g:

ff1; f2g =

@f1 @f2

 

@f1 @f2

+

@f1 @f2

 

@f1 @f2

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

@p3 @x3

@x3 @p3

@p2

@x2

@x2

@p2

 

 

 

 

 

 

 

 

 

 

 

 

@f1 @f2

 

@f1 @f2

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

= x2

p32p2

= f3

 

 

 

 

 

 

@p1

@x1

@x1

@p1

Получили, что ff1; f2g = f3.

fH; f3g = fH; ff1; f2gg

Воспользуемся тождеством Якоби:

fH; ff1; f2gg + ff1; ff2; Hgg + ff2; fH; f1gg = 0

Из условия задачи следует, что fH; f1g = fH; f2g = 0. Получаем fH; f3g = 0.

43

8.2 Решить начально-краевую задачу

8@2u = @2u

<

> @t2 @x2

ujt=0 = 2 cos 52x ; @u@t jt=0 = 0

>

:@u@x jx=0 = 0; ujt= = 0

Решение:

Решение уравнения будем искать в виде u = T (t)X(x).

T 00X = T X00

TT00 = XX00 =

X00 + X = 0

T 00 + T = 0

Ищем нетривиальное решение задачи. 1) < 0

X00 + X =

0

 

 

 

k2 = ; k = p

 

 

 

p

 

 

p

Xxjx=0 = C1 C2 = 0 ) C1 = C2

p p p p

Xj =0 = C1e + C1e = C1e (e2 + C1) = 0 )

C1 = 0

X = 0 получили тривиальное решение 2) = 0

X00 = 0

X= C1x + C2

Xxjx=0 = C1 = 0

X= C2; Xjx= = C2 = 0

X= 0 получили тривиальное решение

3) > 0

X00 + X = 0 p

k2 = ; kp= i p

X = C cos x + C sin x

1 p p 2 p p

X = C1 p sin x + C2 cos x

Xxjx=0 = C2p= 0 ) C2 = 0

Xjxp= = C1 cos = 0

cos = 0 p

= 2 + n; n = 0; 1; 2; :::

44

p

 

=

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n =

 

 

1

+ n 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

2

 

 

 

 

 

2 +

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( n

= 1 + n

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X

 

 

 

 

 

2

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= cos

 

 

 

 

 

 

 

 

n x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u =

1

 

Tn(t) cos

 

21

+ n

x

 

 

1

 

 

 

 

 

 

 

1nP

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8n=0 Tn00Xn

 

= n=0 TnXn00 = n=0 nTnXn

 

 

 

 

P

 

 

 

 

 

 

 

 

 

 

 

 

P

 

 

 

 

5x

 

 

 

P

 

 

 

 

 

>

 

 

Tn(0)Xn = 2 cos

 

2

 

 

 

 

 

 

 

 

 

 

 

 

<

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

n

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n00P

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

>

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n=0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

=

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Tn = C3 cos p

 

 

t + C4 sin p

 

 

t

 

 

 

 

 

n

n

 

 

 

 

 

T

0 =

 

 

C p

 

 

 

 

sin p

 

t + C

 

 

p

 

 

cos p

 

 

 

 

n

 

4

 

n

 

t

 

n

 

 

3

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

n

 

T

0

 

 

 

= p

 

C

= 0

!

C = 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

njt=0

 

 

 

 

 

 

n

 

 

14

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n(0) cos

 

2

+ n

 

 

= 2 cos 2

 

 

Tn = C3 cos

 

 

 

2 + n

t

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

5x

 

 

P= 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

n=0

n

Коэффициенты при остальных n равны нулю.

T2(0) = 2

T2 = C3 cos 52 t

T

2jt=0

= C3 = 2

 

 

 

 

 

 

 

 

5

 

 

 

 

 

 

T2 = 2 cos

2 t

 

 

 

 

 

 

Ответ:

u(x; t) = 2 cos

 

2t cos

 

2x

 

 

 

 

 

 

 

 

5

 

 

5

 

 

 

 

 

 

 

 

 

 

45

Задача 8.2 Закодировать следующую последовательность сообщений алго-

ритмом Шеннона-Фено

 

1

2

3

4

 

 

 

 

 

P ( )

1

1

1

1

 

2

4

8

8

Построить кодовое дерево и вычислить эффективность данного кода.

Решение:

Метод Шеннона-Фено:

1.Составить список элементарных сообщений в порядке убывания вероятностей их появления

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

3.Каждая группа разбивается на 2 подгруппы и применяется тот же принцип и т.д.

 

1

2

3

4

 

 

 

 

 

P ( )

1

1

1

1

2

4

8

8

 

i = 1

A

B

B

B

 

 

 

 

 

i = 2

-

A

B

B

 

 

 

 

 

i = 3

-

-

A

B

V ( )

A

BA

BBA

BBB

 

 

 

 

 

V ( ) - набор элементов кодовой последовательности. Построим кодовое дерево:

46

Вычислим эффективность кода.

 

 

=

l min

 

TEK

l

 

 

TEK - эффективность текущего кодирования

0 < TEK 6 1; åñëè TEK = 1, то текущее кодирование является оптимальным.

l - средняя длина кодового слова

P

l = P ( )l( )

l( ) - длина кодового слова

l min = H log D

D - основание кода

D = 2; log D = 1

H - средняя энтропия Больцмана-Шеннона (входного сообщения)

47

H = PP ( ) log P1( )

l = 12 1 + 14 2 + 18 3 + 18 3 = 74

lmin = 12 log 2 + 14 log 4 + 18 log 8 + 18 log 8 = 74

TEK = 74 47 = 1 ) построили оптимальный код. Теорема. Средняя длина кодового слова l меньшая, чем

является недостижимой ни при каком кодировании

l > H log D

H

log D

Определение. Текущее кодирование - поэлементное кодирование, т.е. каждому элементу ставится во взаимно однозначное соответствие кодовое слово.

Определение. Текущий код l = lmin называется оптимальным текущим кодом.

48

Задача 9.1.

Пусть A и B эрмитовы n n матрицы. Для A известны ее

собственные значения 1; : : : ; n и собственные векторы Y1; : : : ; Yn. Вычислить главный член асимптотики при " ! 0 величины

(et(A+"B)Yj; Yk)

при условии, что j 6= k.

Решение.

Воспользуемся рядом Ньютона:

 

 

2

1

3

"2

 

 

 

 

 

f(A + "B) = f(A) + " B f(A; A) +

:

 

 

 

 

Разложим функцию et(A+"B) в ряд Ньютона:

O

 

 

 

 

 

et(A+"B) = etA + " B

1

d t exp [t A

+t(1 ) A] + O "2

=

 

 

 

Z

 

 

 

2

 

1

 

3

 

 

 

 

 

 

 

0

1

 

 

 

 

 

 

 

= etA + "t Z

d e(1 )tABe tA + O "2

 

(19)

 

 

0

 

 

 

 

 

 

 

Тогда

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

(et(A+"B)Yj; Yk) = (etAYj; Yk) + ("t Z d e(1 )tABe tAYj; Yk) + O "2

 

 

 

0

 

 

 

 

 

 

 

Сначала вычислим (etAYj; Yk). Воспользуемся следующей форму-

E x = (x; Y )Y

(ò.ê. = ).

 

 

n

 

 

jP

ëîé. Òàê

êàê

матрица

A эрмитова, то

f(A)

=

f( j)E j ,

 

 

 

j 6

 

 

 

 

=1

j

j j

k

 

 

 

 

 

n

 

 

n

 

n

 

 

 

X

 

X

Xj

 

 

f(A)Yi =

 

f( j)E j Yi =

 

f( j)(Yi; Yj)Yj =

f( j) ijYj = f( i)Yi

 

j=1

 

 

j=1

=1

 

 

49

Применим эту формулу.

(etAYj; Yk) = (et j Yj; Yk) = et j (Yj; Yk) = et j ij

Учитывая условие задачи j 6= k получаем, что (etAYj; Yk) = 0.

 

 

 

1

 

 

 

 

 

 

Теперь рассмотрим ("t R0

d e(1 )tABe tAYj; Yk).

1

 

 

 

 

1

 

 

 

 

("t Z0

d e(1 )tABe tAYj; Yk) = ("t Z0

d e(1 )tABe t j Yj; Yk) =

 

 

1

 

 

 

 

1

 

 

= ("t Z0

d e t j e(1 )tABYj; Yk) = "t Z0

e t j (e(1 )tABYj; Yk)d =

 

1

 

 

 

 

1

 

 

= "t Z0

e t j (BYj; e(1 )tAYk)d = "t Z0

e t j (BYj; e(1 )t k Yk)d =

 

1

 

 

 

 

 

1

 

 

= "t Z0

e t j e(1 )t k (BYj; Yk)d = "t Z0

e (t j t k)d et k (BYj; Yk) =

 

 

 

= "t

et( j k) 1

et k (BYj

; Yk) = "

et j et k

(BYj; Yk)

 

 

 

 

t( j k)

 

 

( j k)

Ответ.

Получили формулу:

(et(A+"B)Yj; Yk) = "et j et k (BYj; Yk) + O "2 :j k

50