- •Министерство образования и науки Российской Федерации
- •2.2. Энергия и мощность сигнала ………………………………………………...11
- •1. ОСНОВНЫЕ ПОНЯТИЯ ЦИФРОВОЙ ОБРАБОТКИ СИГНАЛОВ
- •1.1. Понятие о первичной и вторичной обработке сигналов
- •1.2. Технические средства комплекса обработки сигналов
- •2. ПОНЯТИЕ СИГНАЛОВ. ВИДЫ СИГНАЛОВ
- •2.2. Энергия и мощность сигнала
- •2.3. Представление периодических сигналов в частотной области
- •2.4. Представление в частотной области непериодических сигналов
- •2.4.1. Введение в теорию ортогональных преобразований
- •2.5. Свойства преобразования Фурье
- •2.6. Интегральное преобразование Хартли
- •2.7. Случайные сигналы
- •2.7.1.Модели случайных процессов
- •Числовые характеристики
- •Примеры случайных процессов с различными законами распределения
- •3. КОРРЕЛЯЦИОННЫЙ АНАЛИЗ СИГНАЛОВ
- •3.1. Корреляционная функция (КФ):
- •3.2. Взаимная корреляционная функция
- •3.3. Взаимный спектр сигналов
- •3.4. Корреляционные функции случайных процессов
- •3.4.1. Стационарные и эргодические случайные процессы
- •3.5. Спектральные характеристики случайных процессов
- •3.5.1. Теорема Винера-Хинчина
- •3.6. Комплексная огибающая сигнала
- •4. ПЕРЕХОД ОТ АНАЛОГОВЫХ СИГНАЛОВ К ЦИФРОВЫМ
- •4.1. Дискретизация сигналов
- •Изменение частоты дискретизации. При решение различных задач обработки сигналов достаточно часто требуется изменение частоты дискретизации сигнала.
- •5.1. Линейные и нелинейные преобразования
- •Импульсная характеристика
- •Переходная характеристика
- •Комплексный коэффициент передачи (передаточная функция) системы:
- •Коэффициент передачи по мощности:
- •Взаимный спектр входного и выходного сигналов
- •Взаимная корреляция между входом и выходом
- •Корреляционная функция
- •Дисперсия на выходе:
- •5.3. Циклическая свертка и корреляция
- •5.5. Двумерная апериодическая свертка и корреляция
- •5.6 Нерекурсивные и рекурсивные фильтры
- •5.7. Метод синхронного или когерентного накопления
- •5.8. Адаптивные фильтры.
- •Рис.5.5. Адаптивный фильтр
- •5.8.1. Фильтр Винера-Хопфа.
- •5.10. Фильтр Калмана.
- •F=kBNX
- •6.1. Дискретное преобразование Фурье
- •6.2. Дискретное преобразование Хартли
- •6.3. Двумерные дискретные преобразования Фурье и Хартли
- •6.4. Ортогональные преобразования в диадных базисах
- •6.7. Выполнение фильтрации в частотной области
- •1) Фильтры нижних частот (ФНЧ) low-pass filter
- •2) Фильтры верхних частот (ФВЧ) hight-pass filter
- •3) Полосовые фильтры (ПФ) band-pass filter
- •4) Режекторные фильтры (ПФ) band-stop filter
- •Фильтр Баттерворта:
- •Фильтр Чебышева 1-го рода:
- •7. ВЕЙВЛЕТ ПРЕОБРАЗОВАНИЯ ИЛИ РАЗЛОЖЕНИЕ ПО ВСПЛЕСКАМ
- •7.1. Понятие о Wavelet-преобразованиях. Преобразование Хаара
- •7.2. Вейвлеты
- •7.2.1. Непрерывные вейвлет преобразования
- •7.2.2. Частотный подход к вейвлет преобразованиям
- •7.2.4. Дискретное вейвлет-преобразование
- •7.2.4.1. Условия полного восстановления сигнала
- •7.2.5. Пакеты вейвлетов (алгоритм одиночного дерева)
- •7.2.6. Целочисленное вейвлет-преобразование
- •Вейвлет-преобразование Лэйзи. Вейвлет-преобразование Лэйзи заключается в простом разбиении входного сигнала на четную и нечетную части. На этапах декомпозиции и реконструкции используются одни и те же формулы:
- •7.3. Применение вейвлет-преобразований для сжатия изображения
- •8. БЫСТРЫЕ АЛГОРИТМЫ ОРТОГОНАЛЬНЫХ ПРЕОБРАЗОВАНИЙ
- •8.1. Вычислительная сложность ДПФ и способы её сокращения
- •8.2. Запись алгоритма БПФ в векторно-матричной форме
- •8.3. Представление алгоритма БПФ в виде рекурсивных соотношений
- •8.4. Алгоритмы БПФ с прореживанием по времени и по частоте
- •Рис.8.4. Граф БПФ
- •8.6. Вычислительная сложность алгоритмов БПФ
- •8.7. Выполнение БПФ для случаев
- •8.8. Быстрое преобразование Хартли
- •Рис.8.7. Граф базовой операции БПХ, где
- •8.9. Быстрое преобразование Адамара
- •8.10. Выбор метода вычисления свертки / корреляции
- •9.1. Ранговая фильтрация
- •Рис. 9.1. Гистограмма распределения элементов по уровням
- •Гистограммный алгоритм ранговой фильтрации для окна размером М х М может быть представлен в следующем виде [16,21]:
- •9.2. Взвешенная ранговая фильтрация
- •9.3. Скользящая эквализация гистограмм
- •9.4. Преобразование гистограмм распределения
- •Рис.9.3. Глобальная эквализация гистограмм
- •ЛИТЕРАТУРА
11
а) Потенциальный сигнал (сигнал постоянного уровеня, например, «0» → логический «0», 5В → логическая «1»). Допускает только амплитудную модуляцию.
Б) Гармонический сигнал:
|
|
|
S(t) = Acos(ωt +ϕ) |
|
||
A |
– амплитуда; |
ω = |
2π |
= 2π f |
– круговая частота; ϕ |
– начальная |
|
|
|
T |
|
|
|
задержка или фаза. Допускает амплитудную, фазовую или частотную модуляцию или их одновременную комбинацию.
В) Дельта-функция δ(t) или единичный импульс
t0, t ≠ 0
δ( ) = ∞, t = 0
∞∫δ(t)dt =1
−∞
Фильтрующее свойство:
∞∫ f (t)δ(t −t0 )dt =f (t0 )
−∞
∞
comb(t) = ∑δ(t −tk ) – гребенчатая функция
k=−∞
г) Функция Хевисайда (функция включения) или едичный скачок
0, t < 0
σ(t) = 12, t = 0
1, t > 0
д) |
Прямоугольный импульс: |
|
|
|
|
|
|
|
|
|
|
|||
|
S(t) = A(σ(t) −σ(t −t0 )) |
|||||||||||||
|
rect |
( |
t |
T ) |
1, |
|
|
t |
|
|
|
≤T |
2 |
|
|
|
|
||||||||||||
|
= |
|
|
|
|
|
|
|
|
|
||||
|
|
t |
>T |
|
||||||||||
|
|
|
0, |
|
|
2 |
||||||||
|
|
|
|
|
|
|
|
Сигналы могут описываться действительной функцией или, существенно реже, комплексной функцией:
S(t) = A e− jω0t
2.2. Энергия и мощность сигнала
Энергия сигнала:
E = T∫S2 (t)dt |
(B2 C) |
(2.1) |
0
Мгновенная мощность (instantaneous power):
P(t) = S2 (t)
Средняя мощность (average power):
12
P = |
1 |
T S2 (t)dt |
(B2 ) |
(2.2) |
|
T |
|||||
ср |
∫0 |
|
|
||
|
|
|
|
Если энергия сигнала бесконечна как, например, у периодических сигналов, то:
Pср = Tlim→∞ T1 T∫2 S2 (t)dt
−T2
Среднеквадратическое (действующее) значение сигнала (root mean square;
RMS)
|
|
|
|
|
|
|
|
|
|
|
|
σ |
|
= |
|
= lim |
1 |
T 2 |
S2 (t)dt |
(2.3) |
|||
S |
P |
||||||||||
|
|||||||||||
|
|
ср |
|
T →∞ T −T∫ |
|
|
|
||||
|
|
|
|
|
2 |
|
|
2.3. Представление периодических сигналов в частотной области
Периодические сигналы могут быть описаны в виде суммы (или суперпозиции) гармонических составляющих или гармоник, каждая из которых имеет определённую частоту, амплитуду и начальную фазу. Конкретный набор таких составляющих будет определяться видом сигнала S(t). Для того, чтобы
такое представление можно было бы осуществить, фрагмент сигнала на периоде T должен удовлетворять условиям Дирихле, т.е.:
1)не должно быть разрывов II рода;
2)число разрывов I-го рода (или скачков) должно быть конечным;
3)число экстремумов должно быть конечным.
При соблюдении этих требований периодический сигнал S(t) может быть представлен в виде ряда Фурье:
|
A0 |
∞ |
|
|
S(t) = |
+ ∑ak cos(kω1t) +bk sin(kω1t), |
(2.4) |
||
2 |
||||
|
k=1 |
|
где ω1 = 2Tπ – круговая частота или период повторения сигнала;
A0 = 1 T∫2 S(t)dt – постоянная составляющая сигнала;
2 T −T 2
ak = 2 T∫2 S(t)cos(kω1t)dt ;
T −T2
bk = 2 T∫2 S(t)sin(kω1t)dt ;
T −T2
kω1 – k -я частотная составляющая сигнала или k -я гармоника.
Заметим, что пределы интегрирования могут быть и другими, например, от 0 до T – важно, чтобы охватывался бы лишь весь период сигнала S(t).
Если S(t) – чётная функция (это значит, что S(t) ≡ S(−t) ), то все bk ≡1 и, наоборот, если S(t) – нечётная функция (S(t) = −S(−t) ), то все ak ≡ 0 .
13
Такое разложение можно записать в тригонометрической форме:
|
|
|
|
A0 |
∞ |
|
|
|
|
S(t) = |
+ ∑Ak cos(kω1t +ϕk ) , |
(2.5) |
|
|
|
|
|
|||
|
|
|
2 |
k=1 |
|
|
где A = |
a2 |
+ b2 |
– амплитуда k -й гармоники; |
|
||
k |
k |
k |
|
|
|
|
ϕk = arctg(bk |
ak ) – начальная фаза. |
|
||||
При этом если S(t) – чётная, |
то ϕk [0;π] и если |
S(t) – нечётная, то |
||||
ϕk [±π 2]. |
|
|
|
|
|
|
Множество амплитуд гармоники называют амплитудным спектром, а множество фаз – фазовым спектром.
Если S(t) – действительная функция, то:
A−k = Ak
ϕ−k = −ϕk .
Пример: Пусть исходный сигнал представляет собой последовательность прямоугольных импульсов с постоянным периодом Т (рис. 2.1).
|
S(t) |
|
∆t |
−∆t 2 |
∆t 2 |
|
∆t |
|
T |
Рис.2.1. Последовательность прямоугольных импульсов
Т.к. такой сигнал чётный, то надо определить
|
2 ∆t 2 |
2πk |
|
2A |
|
2π∆t |
|||
ak = |
|
|
Acos |
|
t dt = |
|
sin |
|
|
T −∆∫t 2 |
T |
πk |
T |
||||||
|
|
|
|
q – скважность; q =T∆t ;
q−1 = ∆Tt – коэффициент заполнения (duty cycle);
A0 = AT∆t = qA
|
A |
∞ |
2A |
|
πk |
2πk |
|
|||
S(t) = |
|
+ ∑ |
|
sin |
|
cos |
|
t |
||
q |
πk |
q |
T |
|||||||
|
k=1 |
|
|
|
|
Амплитудный спектр подобного сигнала показан на рис.2.2..