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

Physics_(i)

.PDF
Скачиваний:
9
Добавлен:
19.04.2015
Размер:
122.05 Кб
Скачать

Вычислим её непосредственно (А.Ф.П.) и с помощью ”те-

лескопирования сумм”(Щепин).

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

eftb +ist − kseftb !

(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(ω) = √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

 

2e83 fbs2

−1/2

eftb +ist

 

eftb

s

fb

 

 

 

 

4

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(12)

 

 

(s = 2πfc, где fc - обозначение из Матлаба).

26

О численных квадратурах.

 

 

Z

 

t2e−t

/2dt = 1

1

 

+

2

 

 

 

 

−∞

 

Если вычислять этот интеграл численно, например, с помо-

щью матлабовской QUAD, то

1

+20000

2

 

 

Z−20000

t2e−t

/2dt = 0

(должен быть ≈ 1). (Матлабовский int - символьные вычи-

сления, поэтому вычисляет правильно)

27

Целая функция Целая функция - функция, голоморфная во всей комплекс-

ной плоскости. Типичным примером целой функции может служить многочлен или экспонента, а также суммы, про-

изведения и суперпозиции этих функций. Ряд Тейлора це-

лой функции сходится во всей плоскости комплексного пере-

менного. Логарифм, квадратный корень не являются целыми функциями.

Отметим, что целая функция может иметь особенность (в

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

должна быть постоянной (это свойство может быть исполь-

зовано для элегантного доказательства основной теоремы ал-

гебры).

Целая функция, имеющая на бесконечности полюс, долж-

на быть многочленом. Таким образом, все целые функции,

не являющиеся многочленами (в частности, тождественно

28

постоянными) имеют на бесконечности существенно особую точку. Такие функции называются трансцендентными целы-

ми функциями.

Малая теорема Пикара значительно усиливает теорему Лиувилля: не равная тождественно постоянной целая функ-

ция принимает все комплексные значения, кроме, возможно,

одного. Примером является экспоненциальная функция, при-

нимающая в качестве значений все комплексные числа, кро-

ме нуля.

Дж. Литлвуд в одной из своих книг указывает сигма-

функцию Вейерштрасса в качестве ”типичного”примера це-

лой функции.

29

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