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

книжка умняшкина по моцос

.pdf
Скачиваний:
133
Добавлен:
16.04.2013
Размер:
1.55 Mб
Скачать

Министерство образования Российской Федерации

Московский государственный институт электронной техники (технический университет)

______________________________________________________________

С.В. Умняшкин

Математические основы цифровой обработки и кодирования сигналов

Учебное пособие

Утверждено редакционно-издательским советом института

Москва 2004

ББК 32.81 У54

УДК 621.372

Рецензенты: доктор физ.-мат. наук, проф. А.В. Алиев доктор техн. наук, проф. В.П. Дворкович

Умняшкин С.В.

У54 Математические основы цифровой обработки и кодирования сигналов: Учебноепособие. - М.: МИЭТ, 2004. - 176 с.: ил.

ISBN 5-7256-0366-0

Учебное пособие предназначено для студентов, обучающихся по направлениям «Информатика и вычислительная техника» и «Прикладная математика». Материал пособия включает в себя в основном наиболее общие теоретические вопросы, связанные с цифровым представлением сигналов и основами анализа линейных дискретных систем, на базе чего в дальнейшем возможно построение специализированных курсов цифровой обработки сигналов, ориентированных на конкретные специальности и направления. По этой причине учебное пособие может быть также рекомендовано для направлений радиотехнического и телекоммуникационного профиля, так как изложение материала опирается лишь на общеобразовательные математические курсы, читаемые в МИЭТе для всех технических специальностей.

ISBN 5-7256-0366-0

© МИЭТ, 2004

Содержание

Введение

5

Глава 1.

Элементы функционального анализа

6

1.1. Неравенства Гельдера и Минковского

6

1.2.

Линейные нормированные пространства

7

1.3. Анализ в линейных нормированных пространствах

9

1.4.

Банаховы пространства

12

1.5. Пространства со скалярным произведением

13

1.6. Аппроксимация в гильбертовом пространстве

16

1.7. Примеры ортогональных систем в пространстве L2

22

Глава 2. Спектральное представление функций

28

2.1. Тригонометрические ряды Фурье. Интеграл Фурье

28

2.2.

Обобщенное преобразование Фурье

35

2.3.

Принцип неопределенности время-частотного

 

 

представления сигналов

38

2.4.

Энергетический спектр. Спектр мощности

43

Глава 3. Дискретизация и квантование сигналов. Дискретные

 

ортогональные преобразования

46

3.1. Преобразование непрерывных сигналов в дискретные

46

3.2. Дискретизация по критерию наибольшего отклонения

47

3.3. Частотный критерий выбора шага дискретизации

48

3.4.

Спектр дискретного сигнала. Спектральная иллюстрация

 

 

теоремы отчетов

54

3.5.

Дискретное преобразование Фурье (ДПФ)

56

3.6.

Быстрое преобразование Фурье (БПФ)

61

3.7.

Дискретное преобразование Уолша

68

3.8.

Дискретное преобразование Хаара

71

3.9.

Некоторые применения дискретных ортогональных

 

 

преобразований

73

3.10.

Квантование дискретных сигналов

75

Глава 4. Линейные дискретные системы

79

4.1.

Z-преобразование

79

4.2.

Линейные дискретные фильтры (ЛДФ)

84

4.3.

Соединения и структурные схемы фильтров

89

4.4.

Устойчивость ЛДФ

92

4.5.

Частотная характеристика ЛДФ

95

Глава 5. Элементы прикладной теории информации

100

5.1. Дискретный источник сообщений без памяти,

 

 

количество информации

100

5.2.

Эффективное кодирование сообщений

 

 

дискретного источника без памяти

106

5.3.

Условная энтропия

113

5.4. Кодирование сообщений источника с памятью

116

5.5. Неопределенность непрерывного источника сообщений.

 

 

Дифференциальная энтропия

120

Глава 6. Теоретические основы применения ортогональных

 

преобразований для представления дискретных сигналов

123

6.1. Корреляция как мера статистической зависимости данных

123

6.2.

Преобразование Карунена-Лоэва

125

6.3.Эффективность использования дискретных ортогональных преобразований для кодирования коррелированных данных 127

6.4. ДПФ в вещественной форме

132

6.5. Дискретный марковский процесс первого порядка.

 

Дискретное косинусное преобразование

133

6.6. Аппроксимационный подход к выбору преобразований для

 

кодирования дискретных сигналов. Частотная трактовка

138

6.7. Время-частотный анализ. Оконное преобразование Фурье

143

Глава 7. Вейвлет-преобразования и их приложения для обработки

дискретных сигналов

147

7.1.

Кратно-масштабный анализ (КМА)

147

7.2. Проецирование функций на подпространства КМА

152

7.3.

Вычисление вейвлет-преобразований

158

7.4.

Квадратурно-зеркальные фильтры

163

7.5.

Биортогональные вейвлет-преобразования

169

7.6. Применение вейвлет-преобразований для сжатия данных

171

Рекомендуемая литература

174

4

Введение

На сегодняшний день вопросы, связанные с цифровой обработкой и кодированием информации, перестали быть узкоспециальными, основы знаний в этой области требуются большинству специалистов с высшим техническим образованием. В последние десятилетия методы цифровой обработки сигналов (ЦОС) в радиотехнике и радиоэлектронике, системах связи, контроля и управления приобрели большую важность и постепенно вытесняют методы аналоговой обработки. Способствует этому стремительно увеличивающаяся производительность вычислительной техники, сопровождаемая её проникновением во все области человеческой деятельности. К цифровому виду хранения и передачи звука все уже привыкли: сегодня, например, виниловые звуковые диски являются такой же редкостью, как 15 лет назад – цифровые компакт-диски. Никого уже не удивляют и цифровые видеозаписи, которые можно воспроизвести на обычном персональном компьютере.

Важность изучения методов и алгоритмов ЦОС, таким образом, трудно переоценить. Для этого, прежде всего, необходимо понимание тех математических основ теории, на которые они опираются. В главах 1-5 мы рассмотрим общие теоретические вопросы, связанные с представлением (фактически, кодированием) сигналов в цифровом виде,

иосновы теории линейных систем обработки дискретных сигналов. На этой теоретической базе строятся специализированные методы ЦОС, изучение которых предстоит в дальнейшем, в рамках специальных курсов. Материал, приведенный в главах 6 и 7, представляет собой более углубленное теоретическое рассмотрение разделов ЦОС, связанных с эффективным представлением сигналов, и ориентирован в основном на студентов, обучающихся по специальностям «Прикладная математика»

и«Программное обеспечение вычислительной техники и автоматизированных систем».

5

Глава 1. Элементы функционального анализа

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

1.1. Неравенства Гёльдера и Минковского

Не останавливаясь на доказательствах, которые хорошо известны и приведены в многочисленных изданиях, сформулируем следующие важные теоремы.

Теорема 1. (Неравенство Гёльдера для интегралов.) Пусть существуют

 

 

 

 

 

 

 

b

 

 

 

b

 

 

 

 

и отличны от нуля интегралы:

 

u(x)

 

p dx ,

 

v(x)

 

q dx

(пределы интег-

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

a

 

 

 

 

рирования – не обязательно конечные), где p>0, q>0,

1

+

1

=1 . Тогда

p

q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

существует также интеграл

 

u( x)v( x)

 

dx и верна оценка:

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

b

 

 

1 p b

 

 

1 q

 

 

 

 

u(x)v( x)

dx

u(x)

p dx

 

v(x)

q dx .

 

 

 

(1)

a

a

 

 

 

 

 

a

 

 

 

 

 

 

 

 

Теорема 2. (Неравенство Минковского для интегралов.) Пусть для p≥1

 

 

 

 

 

 

 

 

b

 

p dx ,

b

p dx

 

 

 

 

 

 

существуют интегралы

 

u(x)

 

 

v(x)

 

 

(пределы интегрирова-

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

a

 

 

 

 

 

 

 

ния – не обязательно

 

конечные).

Тогда

существует

также

интеграл

b

 

p dx , причём верна оценка:

 

 

 

 

 

 

 

 

u(x) + v(x)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

 

 

 

 

 

 

1 p

 

b

 

 

 

 

 

1 p

 

b

1 p

 

 

 

u(x) + v(x)

 

p dx

u(x)

p dx

 

+

v(x)

p dx

.

(2)

6

 

a

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1.2. Линейные нормированные пространства

Определение. Множество Е элементов произвольной природы называется линейным пространством, если в нём определены операции сложения элементов x+y и умножения элементов на скаляр λ (вещественное или комплексное число) λx, результатом которых является элемент из того же множества Е, причём выполняются следующие аксиомы.

. x E, y E: x+y=y+x (свойство коммутативности).

. x E, y E, z E: (x+y)+z = x+(y+z) (свойство ассоциативности).

. θE, x E: θ+x = x (существование нулевого элемента).

. x E, λ,µ (λ, µ – скаляры): λ(µx) = (λµ)x.

. Специальные виды умножения на скаляр λ=0 и λ=1: 0x=θ, 1x=x.

. x E, λ,µ: (λ+µ)x = λx +µx (свойство дистрибутивности).

. x E, y E, λ: λ(x+y) = λx +λy. (свойство дистрибутивности).

Назовем противоположным элементом для x E такой элемент y E, что x+y=θ. Из аксиом 5° и 6° видим, что y=(-1)x (элемент x, умноженный на число -1). Будем обозначать противоположный элемент: -x.

В курсе линейной алгебры систематически изучались линейные пространства Rn арифметических векторов размерности n. Приведённое выше аксиоматическое определение обобщает понятие линейного пространства Rn на множества произвольной природы. По аналогии элементы любого линейного пространства также будем называть вектора-

ми, а сами линейные пространства - векторными пространствами. В

том же обычном смысле будем понимать термины «линейная зависимость» (независимость) векторов, и «размерность пространства». Напомним, что число n называется размерностью векторного пространства Е (обозначается n=dimE), если в Е найдётся n линейно независимых ненулевых элементов, а любые (n+1) ненулевых элементов пространства Е будут линейно зависимы. Данное выше аксиоматическое определение линейного пространства допускает его бесконечную размерность.

7

Пример 1. Пусть C[a;b] – множество всех функций, непрерывных на отрезке [a;b]. Является ли это множество линейным пространством и если да, то какова его размерность?

◄ Выполнение аксиом 1°–7° очевидно, нулевым элементом является функция f(x)=0, x [a;b]. Покажем, что C[a;b] – бесконечномерное пространство. Выберем из множества C[a;b] n ненулевых элементов – степенных функций yi=xi-1, i=1,2,…,n. Очевидно, что при любом числе n эти элементы являются линейно независимыми, так как выполнение

n

n

равенства λi yi

= λi xi1 = 0 для всех точек отрезка x [a,b] воз-

i=1

i=1

можно лишь в случае λi=0, i=1,…,n. Поскольку число n можно выбрать как угодно большим, то dim(C[a;b])=. ►

Определение. Множество M E векторного пространства Е называется линейным многообразием, если М включает в себя все линейные комбинации всех своих элементов. То есть x, y M , λ, µ :

(λx + µy) M .

Определение. Линейное пространство Е называется нормированным, если каждому элементу x E поставлено в соответствие число, называемое нормой (обозначается x ), и выполняются следующие

аксиомы.

. Невырожденность нормы: x0 , x = 0 x = θ . . Однородность нормы: λx = λ x .

. Неравенство треугольника. x,y E: x + y x + y .

В одном и том же векторном пространстве Е норму можно вводить различными способами.

Пример 2. Рассмотрим векторное пространство C[a;b] из предыдущего примера. Покажите самостоятельно, что приводимые ниже способы вычисления нормы удовлетворяют аксиомам 1°–3°:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

 

 

 

1 p

 

а)

 

x

 

= max

 

x(t)

 

,

б)

 

 

 

x

 

 

 

=

 

x(t)

 

p dt

, где p1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t [a,b]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

Указание. В пункте б) для доказательства аксиомы треугольника воспользоваться неравенством Минковского.

8

Определение. Расстоянием между элементами x, y нормированного векторного пространства Е назовем число ρ( x, y) = x y .

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

. ρ(x,y)=ρ(y,x).

 

. ρ(x,y)=0 x=y.

 

. x,y,z E: ρ(x,y)ρ(x,z)+ρ(y,z)

(неравенство треугольника).

◄ Первое и второе свойства очевидны. Для неравенства треугольника имеем:

ρ( x, y) = x y = (x z) + (z y)

 

x z + z y = ρ(x, z) + ρ( y, z) .

Расстояние между элементами называют также метрикой пространства. Пространство с введённым правилом нахождения расстояния (не обязательно нормированное) со свойствами 1°-3° называется метриче-

ским.

Определения. В метрическом пространстве E открытым шаром

радиуса

r>0

c

центром

x0 E

назовем

множество

Sr (x0 ) = {x E

 

ρ( x, x0 ) < r},

замкнутым

шаром

– множество

 

 

 

r (x0 ) = {x E

 

ρ( x, x0 ) r},

 

сферой

множество

 

S

 

 

σr (x0 ) = {x E

 

 

ρ( x, x0 ) = r}.

Окрестностью точки x0 E будем на-

 

 

зывать открытый шар произвольного радиуса ε, т.е.

Sε (x0 ) . Если

существует некоторый шар конечного радиуса R<, который содержит в себе все элементы множества X E, то множество X называют

ограниченным.

Понятия нормы, расстояния, окрестности являются исходными для построения анализа в линейных нормированных пространствах.

1.3. Анализ в линейных нормированных пространствах

Определение. В линейном нормированном пространстве (ЛНП) Е

элемент y называется пределом последовательности {xk}, если

lim ρ( y, xk ) = 0 . Говорим, что последовательность {xk} сходится к

k→∞

элементу y и обозначаем это как lim xk = y .

k→∞

9

Определение. Точка a из ЛНП E называется предельной точкой множества M E, если в любой окрестности a содержится хотя бы одна точка x M, отличная от a. То есть r > 0 : (Sr (a) \ a)M .

Почему такая точка называется предельной? Ответ дает следующая теорема.

Теорема 3. Для того, чтобы элемент a из ЛНП E был предельной точкой множества M E, необходимо и достаточно существование последова-

тельности {xk}, xk M, xka, сходящейся к a: lim xk = a .

k→∞

Необходимость. Возьмем сходящуюся к нулю числовую последовательность из ненулевых элементов, например, εk=1/k, k=1,2,… Так как a

– предельная точка M,

то, по определению, εk > 0

xk M , xk

 

 

 

a :

 

 

 

xk Sεk (a) . Поскольку

ρ( xk ,a) = xk a

< εk =1 k

и

lim

 

 

 

xk a

 

 

 

= 0 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

{xk}

 

 

k→∞

 

 

 

 

a:

то так

 

 

 

построенная последовательность

сходится

к точке

 

lim xk

 

= a .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k→∞

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Достаточность. Так как

{xk }, xk M , xk

a ,

причем

lim xk = a ,

то,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k→∞

 

 

 

 

 

по определению предела,

lim

 

 

 

xk a

 

 

 

= 0

и ε > 0 N (ε) :

 

n > N (ε)

 

 

 

 

 

 

 

 

xn a

 

 

 

 

 

k→∞

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

< ε . То есть

в

любой ε-окрестности точки

 

содержатся

 

 

 

 

 

элементы xn M, xna, поэтому точка a является предельной для множества M. ►

Определение. Пусть M – подмножество в ЛНП E, а M' – множество

всех предельных точек M. Объединение множеств M = M M называется замыканием множества M. Если M содержит все свои предельные точки, т.е. M' M, то множество M называется замкнутым. Определение. Замкнутое линейное многообразие L в ЛНП E, L E,

назовём подпространством.

Определение. Расстоянием от точки x из ЛНП E до подпростран-

ства L E называется величина ρ( x, L) = inf x u .

u L

Аналогично вводится и расстояние от точки x множества E до произвольного подмножества M E. В общем курсе математического анализа было показано, что для ограниченного снизу числового множества всегда найдётся точная нижняя грань. Поскольку любая норма неотри-

10