ОТВЕТНИК ВСЕ В ОДНОМ. ФОРМАТ А5
.pdfw1(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