Physics_(i)
.PDFВычислим её непосредственно (А.Ф.П.) и с помощью ”те-
лескопирования сумм”(Щепин).
1. Заметим, что
N
X k2 = 1+(2+2)+(3+3+3)+(4+4+4+4)+. . .+(N+N+. . . N).
k=1
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(7) |
|
||
поскольку |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
N |
k = N |
k |
|
n−1 k = (N + 1)N − n(n − 1), |
||||||||||||||||||||||||||||
|
X |
|
|
|
X |
− |
X |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
k=n |
|
|
|
k=1 |
k=1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
то, учитывая (8), имеем |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
N k2 = |
N k+ |
N k+. . . |
N |
|
|
k+N = N−1 |
|
(N + 1)N − n(n − 1) |
+N = |
||||||||||||||||||||||||
=X |
|
||||||||||||||||||||||||||||||||
X |
X |
|
|
X |
|
|
1 |
|
|
|
|
|
|
|
|
X |
|
|
|
|
|
|
2 |
|
|||||||||
k=1 |
k=1 |
|
|
k=2 |
|
|
k N |
− |
|
|
|
|
|
|
|
|
|
n=1 |
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 N−1 n + N. |
||||||||||
|
= (N + 1)N(N − 1) 1 N−1 n2 + |
|
|||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
− |
|
|
|
X |
|
|
|
|
|
|
|
X |
|
|
|
|
||||
|
|
|
|
|
|
|
2 |
|
|
|
|
|
2 |
n=1 |
|
|
|
|
2 n=1 |
|
|
|
|
||||||||||
Обозначим искомую сумму символом P. Имеем |
|
|
|
|
|||||||||||||||||||||||||||||
|
= |
|
(N2 − 1)N |
+ |
N(N |
|
1) |
+ N |
|
|
1 |
+ |
1 |
N2, |
|||||||||||||||||||
|
|
|
|
|
|
|
− 2 |
2 |
|||||||||||||||||||||||||
X |
|
|
|
2 |
|
|
|
|
|
|
4− |
|
|
|
|
|
X |
|
|
|
|||||||||||||
откуда |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
X = |
2 |
N3 |
+ 3N2 + N) |
|
|
N |
( |
N |
+ |
|
1)(2N + |
1) |
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
= |
|
|
|
|
|
|
|
|
|
. |
|
|||||||||||||
|
|
|
|
6 |
|
|
|
|
|
|
|
|
|
|
|
|
6 |
|
|
|
|
||||||||||||
2. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
21
Разность функций:
f(k) = f(k + 1) − f(k).
Сумма разностей:
n
X F (k) = F (2)−F (1)+F (3)−F (2)+. . . F (n+1)−F (n) = F (n+1)−F (1).
k=1
Т.е. для простого суммирования нужно выразить суммируе-
мую функцию f(k) через разность какой-то другой функции: f(k) = F (k).
x = x + 1 − x = 1,
x2 = (x + 1)2 − x2 = 2x + 1,
|
|
|
|
|
|
x3 = (x + 1)3 − x3 = 3x2 + 3x + 1. |
|
|
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
x3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
− x2 = x + 1/3 |
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
3 |
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
x2 |
|
|
|
|
x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
− |
|
|
|
= x + 1/2 − 1/6 = x + 1/3 |
|
|||||||||||||||||||
|
|
|
|
|
2 |
6 |
|
||||||||||||||||||||||||
Следовательно |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
x2 = |
|
x3 |
x2 |
x |
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
− |
|
|
+ |
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
3 |
|
2 |
|
|
6 |
|
|
|
|
|
|
|
|||||||||
n |
= |
x |
3 |
|
|
x |
2 |
+ x |
= (n + 1) |
3 |
− 1 |
|
(n + 1) |
2 |
− 1+ |
(n + 1) − 1 |
|||||||||||||||
x2 |
|
|
|
|
|
|
|
|
|||||||||||||||||||||||
X |
|
|
|
|
− |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
− |
|
|
|
|
|
|
|
x=1 |
|
3 |
|
|
2 |
|
6 |
|
|
|
|
|
3 |
|
|
|
|
|
|
2 |
|
|
|
6 |
22
Окончательно
n |
= |
N(N + 1)(2N + 1). |
|
x2 |
|||
X |
|
|
|
x=1 |
|
6 |
|
23
9К вейвлетам и Фурье
Для вычислений нужен следующий табличный интеграл (Ма-
ричев, Брычков, Прудников):
Z |
|
|
|
|
√ |
|
|
|
|
q2 |
|
+ |
∞(ax2 |
2 |
|
π |
|
|
|||||
|
+2bx+c)e−px |
−qxdx = |
|
p−5/2[4cp2 |
−4pqb+a(q2+2p)] exp( |
|
) |
||||
|
4 |
4p |
|||||||||
|
−∞ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(Rep > 0). |
|
(8) |
|
|
|||
|
Вейвлет Морле (Goupilland, Grossmann, Morlet, 1984). |
|
|
||||||||
|
|
|
Ψs(t) = Csπ−1/4 |
e−ftb +ist − kse−ftb ! |
(9) |
|
|
||||
|
|
|
|
|
2 |
2 |
|
|
|
Здесь введена величина fb для удобства сравнения с матла-
бовским ”Морлетом”(т.е. s = 2πfc), Cs -нормировочная кон-
станта.
Величина ks находится из критерия допустимости (admissibility criterion):
Z ∞
−∞ Ψs(t)dt = 0 (10)
Используя (8) из (10) получаем ks = e−s2 f4b
24
Преобразование Фурье функции (9), если его определять как
Ψˆ s(ω) = √2π Z ∞ |
Ψs(t)e−iωtdt, |
||
1 |
|
|
|
|
|
−∞ |
|
вычисляется с помощью (8): |
|
ˆ (ω) = C π−1/4 e−(s−ω)2 f4b − k e−ω2 f4b .
Ψs s s
Параметр s позволяет взаимно изменять разрешения по времени и по частоте. Малое s -высокое временное разреше-
ние, но возникают вычислительные проблемы. При s > 5 ks < 10−5, поэтому им (ks) часто пренебрегают. (Ма-
лые значения s нет необходимости использовать, если сиг-
нал содержит только медленно изменяющуюся частоту (!) и
амплитудную модуляцию (аудио)).
Если константу Cs вычислить из следующего условия нор-
мировки
Z ∞ |Ψs(t)|2dt = 1,
−∞
то, с помощью (8) имеем
Cs = |
2 |
!1/4 |
1 + e− |
fb |
s2 |
− |
2e− |
83 fbs2 |
|
−1/2 . |
|
2 |
|||||||||||
|
|||||||||||
|
fb |
|
|
|
|
25
При таком значении нормировочной константы Cs пре-
образование Фурье вейвлета ˆ также нормировано
Ψs,fb(ω) (??)
на единицу (что легко проверить с помощью (8)):
Z ∞ |ˆ |2
−∞ Ψs,fb(ω) dω = 1.
Окончательно, имеем:
Вейвлет Морле:
Ψs,fb(t) = |
fb |
! |
1/4 |
1 + e− |
fb |
s2 |
|
2e−83 fbs2 |
−1/2 |
e−ftb +ist |
|
e−ftb − |
s |
fb |
|
|
|
|
4 |
||||||||||||
2 |
|
|
|||||||||||||
2π |
|
|
|||||||||||||
|
|
− |
|
|
2 |
− |
2 |
|
2 |
|
|||||
|
|
|
|
|
|
|
|
|
|
(11) |
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
Преобразование Фурье ВМ:
Ψˆ s,fb(ω) = |
fb |
!1/4 |
1 + e− |
fb |
s2 |
− |
2e− |
83 fbs2 |
|
−1/2 |
|
e−(s−ω)2 |
fb |
− |
kse−(ω2+s2) |
fb |
. |
|
2 |
4 |
4 |
||||||||||||||||
|
||||||||||||||||||
|
2π |
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(12) |
|
|
(s = 2πfc, где fc - обозначение из Матлаба).
26
О численных квадратурах. |
|
||||
|
√2π Z |
|
∞ t2e−t |
/2dt = 1 |
|
1 |
|
+ |
2 |
|
|
|
|
|
−∞ |
|
Если вычислять этот интеграл численно, например, с помо-
щью матлабовской QUAD, то
1 |
+20000 |
2 |
|
|
√ |
|
Z−20000 |
t2e−t |
/2dt = 0 |
2π |
(должен быть ≈ 1). (Матлабовский int - символьные вычи-
сления, поэтому вычисляет правильно)
27
Целая функция Целая функция - функция, голоморфная во всей комплекс-
ной плоскости. Типичным примером целой функции может служить многочлен или экспонента, а также суммы, про-
изведения и суперпозиции этих функций. Ряд Тейлора це-
лой функции сходится во всей плоскости комплексного пере-
менного. Логарифм, квадратный корень не являются целыми функциями.
Отметим, что целая функция может иметь особенность (в
т.ч. даже существенную особенность) в бесконечности. Как следует из теоремы Лиувилля, функция, которая не имеет особых точек на всей расширенной комплексной плоскости,
должна быть постоянной (это свойство может быть исполь-
зовано для элегантного доказательства основной теоремы ал-
гебры).
Целая функция, имеющая на бесконечности полюс, долж-
на быть многочленом. Таким образом, все целые функции,
не являющиеся многочленами (в частности, тождественно
28
постоянными) имеют на бесконечности существенно особую точку. Такие функции называются трансцендентными целы-
ми функциями.
Малая теорема Пикара значительно усиливает теорему Лиувилля: не равная тождественно постоянной целая функ-
ция принимает все комплексные значения, кроме, возможно,
одного. Примером является экспоненциальная функция, при-
нимающая в качестве значений все комплексные числа, кро-
ме нуля.
Дж. Литлвуд в одной из своих книг указывает сигма-
функцию Вейерштрасса в качестве ”типичного”примера це-
лой функции.
29