'

Спектральный анализ Быстрая свертка Фильтрация

Понравилась презентация – покажи это...





Слайд 0

Спектральный анализ Быстрая свертка Фильтрация Лектор: Лукин Алексей Сергеевич


Слайд 1

План Спектральный анализ Мощность сигнала, теорема Парсеваля Спектрограммы, усреднение во времени Частотно-временное разрешение Быстрая свертка Свертка в частотной области Разбиение ядра Фильтрация Основные понятия Фильтры для изображений, их частотные характеристики Применения: эквалайзеры, реверберация


Слайд 2

Мощность и амплитуда Мощность пропорциональна квадрату амплитуды RMS (среднеквадратичное зачение) сигнала: Окна при вычислении RMS Вычисление RMS в реальном времени: (пропорционально амплитуде) (экспоненциальное окно)


Слайд 3

Теорема Парсеваля Энергия сигнала равна энергии спектра (возможно, с точностью до множителя, зависящего от нормализации в формулах преобразования Фурье) Применение теоремы для оценки энергии различных частотных составляющих сигнала Пример: вычисление THD (total harmonic distortion)


Слайд 4

Спектральный анализ Как вычислить и отобразить спектр сигнала? Взять нужный отрезок сигнала длины 2m; если нужный отрезок короче – дополнить его нулями. Если нужно – умножить сигнал на весовое окно, плавно спадающее к краям. Обычно это улучшает свойства спектра. Вычислить FFT. Перевести комплексные коэффициенты в полярную форму: получить амплитуды. Отобразить график зависимости амплитуды от частоты. Примеры весовых окон


Слайд 5

Свойства ДПФ Линейность Комплексное ДПФ от вещественного сигнала обладает симметрией относительно коэффициента k = N/2. Таким образом, только половина коэффициентов являются значимыми. Периодическое продолжение сигнала Циклический сдвиг отсчётов сигнала не изменяет амплитудного спектра, а меняет только фазовый


Слайд 6

Размытие спектра Как выглядит спектр синусоиды, если её частота совпадает/не совпадает с базисными частотами ДПФ? нет окна (прямоугольное окно)


Слайд 7

Весовые окна Весовые окна уменьшают эффект размытия спектра, подавляя боковые лепестки за счёт расширения главного лепестка окно Хэмминга (Hamming)


Слайд 8

Весовые окна Весовые окна уменьшают эффект размытия спектра, подавляя боковые лепестки за счёт расширения главного лепестка окно Хана (Hann)


Слайд 9

Весовые окна Весовые окна уменьшают эффект размытия спектра, подавляя боковые лепестки за счёт расширения главного лепестка окно Блэкмана (Blackman)


Слайд 10

Весовые окна Весовые окна уменьшают эффект размытия спектра, подавляя боковые лепестки за счёт расширения главного лепестка окно Кайзера (Kaiser), ? = 15


Слайд 11

Весовые окна Прямоугольное (нет окна) Hamming Blackman Kaiser Формулы и картинки: http://en.wikipedia.org/wiki/Window_Function


Слайд 12

Спектральный анализ Отображение спектра звука: спектрограмма Спектрограмма – график зависимости амплитуды от частоты Низкие частоты – слева, высокие – справа Часто применяется логарифмический масштаб частот и амплитуд: “log-log-спектрограмма” Временное и частотное разрешение спектрограммы Децибелы: A1 – амплитуда измеряемого сигнала, A0 – амплитуда сигнала, принятого за начало отсчета (0 дБ) Разница на 6 дБ – разница по амплитуде в 2 раза, разница на 12 дБ – разница по амплитуде в 4 раза. Часто за 0 дБ принимается либо самый тихий слышимый звук, либо самый мощный звук, который может воспроизвести аудиоустройство.


Слайд 13

Спектральный анализ Примеры звуков и их спектров Песня (стерео запись) Нота на гитаре


Слайд 14

Спектральный анализ Отображение спектра звука: спектрограмма (сонограмма) Спектрограмма – график зависимости амплитуды от частоты и от времени, показывает изменение спектра во времени Short Time Fourier Transform (STFT)


Слайд 15

Спектральный анализ Отображение спектра звука: спектрограмма (сонограмма) Спектрограмма – график зависимости амплитуды от частоты и от времени, показывает изменение спектра во времени Низкие частоты – снизу, высокие – сверху Время идет справа налево Амплитуда – яркость или цвет Частотное и временное разрешение Short Time Fourier Transform (STFT) Спектрограмма Осциллограмма (форма волны)


Слайд 16

Спектральный анализ Примеры звуков и их спектрограмм Нота на гитаре


Слайд 17

Спектральный анализ Спектрограммы шума Типы шумов: белый (постоянная PSD), розовый (3 дБ/окт.), коричневый (6 дБ/окт.) Вариация и усреднение спектрограмм


Слайд 18

Построение спектрограмм Мел-шкала частот Недостаточное частотное разрешение на НЧ, недостаточное временное разрешение на ВЧ Обычная спектрограмма, окно 512 точек Обычная спектрограмма, окно 4096 точек


Слайд 19

Построение спектрограмм Адаптивное частотно-временное разрешение спектрограмм Критерий оптимальности: максимально компактная локализация энергии (минимальное размытие спектрограммы) Адаптивная спектрограмма, разрешение зависит от частоты Адаптивная спектрограмма, оптимальное разрешение


Слайд 20

Построение спектрограмм Критерий оптимальности: максимальное компактирование (локализация) энергии 6 ms 12 ms 24 ms 48 ms 96 ms лучшее разрешение Здесь ai,r – амплитуды STFT в блоке разрешения r, Sr – мера размытости спектра для данного разрешения r, r0 – лучшее разрешение, для которого размытость минимальна.


Слайд 21

Свертка и фильтрация Основные термины Свертка (convolution), фильтрация (filtering) Фильтр (filter), ядро фильтра (kernel) Импульсная, частотная и фазовая характеристики (impulse, frequency, phase response) Применения фильтрации Анти-алиасинг изображений, нахождение границ Звуковой эквалайзер Моделирование реверберации помещения


Слайд 22

Быстрая свертка Прямое вычисление: M·N умножений (M – размер ядра свертки, N – длина сигнала) Теорема свертки: свертка* во временной области эквивалентна умножению в частотной области, умножение во временной области эквивалентно свертке* в частотной области. Алгоритм быстрой свертки: Вычислить спектры сигнала и ядра свертки (FFT) Перемножить эти спектры Вернуть полученный спектр во временную область (IFFT) Почему это быстрее? Потому что переход в частотную область и обратно быстрый: FFT * Речь идет о т.н. круговой свертке


Слайд 23

Быстрая свертка Как изменяется длина сигнала при свертке? Она увеличивается на длину ядра минус 1 (т.к. каждый входной отсчет превращается в ядро и они складываются с наложением) Значит, если взять сигнал длины N, ядро длины M и произвести свертку через FFT размера N, то результат свертки (длины N+M-1) не поместится в результате IFFT (длины N). Произойдет круговая свертка (заворачивание результата по времени). Следовательно, для предотвращения круговой свертки надо взять размер FFT как минимум N+M-1


Слайд 24

Вычисление свертки Быстрая свертка через FFT Неоптимальна по скорости, когда длины сигнала и ядра сильно различаются Неприменима при потоковой обработке (когда входной сигнал поступает в реальном времени и его длина неизвестна) Секционная свертка Используем свойство линейности и инвариантности к сдвигу: разбиваем сигнал на блоки и делаем быструю свертку блочно Увеличение длины каждого блока на M-1 > сложение результатов с перекрытием (метод OLA: overlap-add) Подробности – см. в методичке


Слайд 25

Разбиение ядра Проблема секционной свертки: сигнал надо накапливать блоками, сравнимыми по размеру с длиной ядра свертки (иначе – неэффективно) > задержки при real-time Partitioned convolution (свертка с разбиением ядра) Цель – уменьшение задержки Недостаток – некоторое увеличение времени обработки Секционная свертка: Разбиение ядра: (раскрываем скобки)


Слайд 26

Частотная характеристика Получение частотной характеристики по ядру Дополнение нулями FFT Берем модули комплексных коэффициентов Свертка = перемножение частотных характеристик * = Перемножение амплитуд = сложение децибелов


Слайд 27

Виды фильтров Идеальные и реальные фильтры, виды искажений FIR (finite impulse response) и IIR (infinite i. r.) фильтры Линейность ФЧХ («линейность фазы»)


Слайд 28

Построение фильтров Простейший НЧ-фильтр – windowed sinc Метод спектральной инверсии* Вычитание фильтра из единичного импульса – вычитает его АЧХ из единичной Инверсия каждого второго отсчета в фильтре – разворачивает АЧХ фильтра по частоте * Применим для фильтров с линейной ФЧХ


Слайд 29

Построение фильтров Метод весового окна Построение фильтра с линейной фазой по произвольной заданной частотной характеристике Частотная характеристика приближается с любым заданным уровнем точности Более сложные методы Чуть меньшее число точек ядра при соблюдении той же частотной характеристики Специальные требования (к фазовой характеристике, к монотонности частотной характеристики и т.п.)


Слайд 30

Построение фильтров Метод весового окна (идея метода – взять обратное FFT от требуемой АЧХ) Дана требуемая АЧХ (частотная характеристика) Произвольно выбираем размер FFT N и берем амплитуды коэффициентов Фурье по АЧХ (N/2+1 точка) Записываем нули в фазы и делаем IFFT Циклический сдвиг ядра Применение весового окна Контроль реальной полученной АЧХ (см. получение АЧХ по ядру фильтра) Если АЧХ недостаточно соответствует требованиям – увеличиваем размер FFT и снова конструируем фильтр Подробности – см. в методичке


Слайд 31

Единичный импульс Простейшее размытие Двумерные фильтры


Слайд 32

Константное размытие 3х3 Константное размытие 5х5 Двумерные фильтры


Слайд 33

Повышение четкости Выделение границ Двумерные фильтры


Слайд 34

Тиснение Пример спектра изображения Двумерные фильтры


Слайд 35

Эквалайзеры Эквалайзер – устройство коррекции тембра сигнала, изменяющее амплитуды его частотных составляющих Изначально применялись для выравнивания АЧХ неидеального звукового тракта Вскоре стали использоваться и творчески, для создания нужных тембров или аккуратного совмещения инструментов в фонограмме


Слайд 36

АЧХ Амплитудно-частотная характеристика (frequency response) Добротность (Q) определяет ширину полосы воздействия


Слайд 37

Виды эквалайзеров По управлению АЧХ: Параметрические: можно выбирать одну из имеющихся форм АЧХ и задавать ее параметры: центральную частоту, коэффициент усиления и добротность Графические: пользователь «рисует» требуемую АЧХ непосредственно на дисплее или с помощью набора регуляторов усиления на различных частотах Параграфические: гибрид предыдущих По принципу действия: Аналоговые: состоят из конденсаторов, катушек индуктивности, операционных усилителей Цифровые: используют FIR или IIR фильтры


Слайд 38

ФЧХ и ГВЗ Фазово-частотная характеристика (phase response) Изменение фазы в зависимости от частоты Групповое время задержки (group delay) Задержка различных частотных составляющих АЧХ ФЧХ ГВЗ


Слайд 39

Импульсная характеристика Конечная (FIR) Легко достичь линейной ФЧХ Бесконечная (IIR) Нелинейная ФЧХ Звон фильтров (ringing) Линейная ФЧХ Нелинейная ФЧХ


Слайд 40

Измерение реверберации Импульсный отклик помещения Реверберация: прямой звук, ранние и поздние отражения


Слайд 41

Измерение реверберации Способы измерения импульсного отклика: Единичный импульс (плохое соотношение сигнал/шум) Звук стартового пистолета / лопающегося воздушного шарика (неровный спектр, маленькая энергия на НЧ) Шум измеряемого зала Используются материалы Дэвида Гризингера (Lexicon) Стартовый пистолет


Слайд 42

Измерение реверберации Способы измерения импульсного отклика: Белый шум, MLS-последовательности (чувствительны к нелинейным искажениям) Скользящий тон (swept sine, «плавающий синус») Деконволюция (обратная свертка) Скользящий тон


Слайд 43

Реверберация Что почитать: Дэвид Гризингер: http://world.std.com/~griesngr/ Ирина Алдошина: http://www.625-net.ru/archive/z0803/aldo.htm


×

HTML:





Ссылка: