Цифровые методы обработки сигналов



Скачать 236.22 Kb.
Дата30.04.2016
Размер236.22 Kb.
Цифровые методы обработки сигналов

Введение

Сущность цифровой обработки состоит в том, что физический сигнал (напряже­ние, ток и т. д.) преобразуется в последовательность чисел, которая затем подвер­нется математическим преобразованиям в вычислительном устройстве. Трансформированный цифровой сигнал (последовательность чисел) при необходимости может быть преобразована обратно в напряжение или ток.

Далее будут рассмотрены принципы математического описания и ана­лиза дискретных сигналов. Но прежде всего следует пояснить некоторые терми­нологические тонкости.
Аналоговые, дискретные и цифровые сигналы

Исходный физический сигнал является непрерывной функцией времени. Такие сигналы, определенные во все моменты времени, называют аналоговыми, последовательность чисел, представляющая сигнал при цифровой обработке, является дискретным рядом и не может полностью соответствовать аналоговому сигналу. Числа, составляющие последовательность, являются колебаниями сигнала в отдельные (дискретные) моменты времени и называются отсчетами сигнала. Отсчеты берутся через равные про­межутки времени Т, называемые п



ериодом дискретизации (или интервалом , ша­лом дискретизации). Величина, обратная периоду дискретизации, называется частотой дискретизации: fд = 1/Т. Соответст­вующая ей круговая частота определяется следующим образом: д = fд.

Ясно, что в общем случае представление сигнала набором дискретных отсчетов приводит к потере информации, так как мы ничего не знаем о поведении сигнала в промежутках между отсчетами. Однако, как будет показано далее, существует класс аналоговых сигналов, для которых такой потери информации не происходит, и которые могут быть точно восста­новлены по значениям своих дискретных отсчетов.


Процесс преобразования аналогового сигнала в последовательность отсчетов на­зывается дискретизацией, а результат такого преобразования — дис­кретным сигналом.

При обработке сигнала в вычислительных устройствах его отсчеты представля­ются в виде двоичных чисел, имеющих ограниченное число разрядов. Вследст­вие этого отсчеты могут принимать лишь конечное множество значений и, сле­довательно, при представлении сигнала неизбежно происходит его округление. Процесс преобразования отсчетов сигнала в числа называется квантованием по уровню , а возникающие при этом ошибки округления — ошибками (или шумами) квантования.

Сигнал, дискретный во времени, но не квантованный по уровню, называется дис­кретным сигналом. Сигнал, дискретный во времени и квантован­ный по уровню, называют цифровым сигналом. Сигналы, квантованные по уровню, но непрерывные во времени, на практике встречаются редко. Раз­ницу между аналоговыми, дискретными и цифровыми сигналами иллюстрирует рис. 1.



Рис. 1. Аналоговый (слева), дискретный (в центре) и цифровой (справа) сигналы

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

В дальнейшем речь пойдет о дискретных сигналах и методах их обра­ботки, поскольку эффекты, связанные с квантованием по уровню, в большинстве случаев не будут приниматься во внимание.
Аналого-цифровое и цифро-аналоговое преобразование

Обобщенная структура системы цифровой обработки сигналов приведена на рис. 2. На вход поступает аналоговый сигнал sBX(t). Его временная дискретизация и квантование по уровню производятся в аналого-цифровом преобразователе (АЦП). Выходным сигналом АЦП является последовательность чисел, поступающая в цифровой процессор ЦП, выполняющий требуемую обработку. Процессор осуществляет различные математические операции над входными отсчетами. Результатом работы процессора является новая последовательность чисел, представляющих собой отсчеты выходного сигнала. Аналоговый выходной сигнал sBbIX(t) восстанавливается по последовательности чисел с помощью цифро-аналогового преобразователя ЦАП. Напряжение на выходе ЦАП имеет ступенчатую форму при необходимости, оно может быть преобразовано в плавно меняющийся выходной сигнал с помощью сглаживающего фильтра Ф.





Рис. 2. Структурная схема системы цифровой обработки сигналов

Частота Найквиста

Гармонический сигнал может быть адекватно представлен дискретными отсчетами, если его частота не превышает половины частоты дискретизации (эта частота называется частотой Найквиста - fN = fд /2 = 1/(2Т); n = д/2 = π/Т). Происхождение этого ограничения поясняет рис. 3. В зави­симости от соотношения между частотой дискретизируемого гармонического сигнала и частотой Найквиста возможны три случая:



  1. Если частота гармонического сигнала меньше частоты Найквиста, дискретные отсчеты позволяют правильно восстановить аналоговый сигнал (рис. 3, а).

  2. Если частота гармонического сигнала равна частоте Найквиста, то дискрет­ные отсчеты позволяют восстановить аналоговый гармонический сигнал с той же частотой, но амплитуда и фаза восстановленного сигнала (он показан пунктирной линией) могут быть искажены (рис. 3, б). В худшем случае все дискретные отсчеты синусоиды могут оказаться равными нулю.

3. Если частота гармонического сигнала больше частоты Найквиста, восстанов­ленный по дискретным отсчетам аналоговый сигнал (как и в предыдущем слу­чае, он показан пунктирной линией) будет также гармоническим, но с иной частотой (рис. 3, в). Данный эффект носит название появления ложных час­тот.



Рис. 3. Дискретизация гармонических сигналов с разной частотой
Ряд Фурье

Разложению в ряд Фурье могут подвергаться периодические сигналы. При этом они представляются в виде суммы гармонических функций либо комплексных экспонент с частотами, образующими арифметическую прогрессию. Для того чтобы такое разложение существовало, фрагмент сигнала длительностью в один период должен удовлетворять условиям Дирихле:

1. не должно быть разрывов второго рода (с уходящими в бесконечность ветвя­ми функции);

2. число разрывов первого рода (скачков) должно быть конечным;

3. число экстремумов должно быть конечным (в качестве примера функции, ко­торая на конечном интервале имеет бесконечное число экстремумов, можно привести sin(l/x) в окрестности нуля).

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



Синусно-косинусная форма

В этом варианте ряд Фурье имеет следующий вид:



Здесь 1 = 2π/Т — круговая частота, соответствующая периоду Т повторения сиг­нала. Входящие в формулу кратные ей частоты k1 называются гар­мониками, гармоники нумеруются в соответствии с индексом k, частота k = k1 называется k - й гармоникой сигнала. Коэффициенты ряда аk и bk рассчитываются по формулам:





Константа a0 рассчитывается по общей формуле для ak. Ради этой общности и введена несколько странная на первый взгляд форма записи постоянного слагае­мого (с делением на два). Само же это слагаемое представляет собой среднее значение сигнала на периоде





Вещественная форма

Некоторое неудобство синусно-косинусной формы ряда Фурье состоит в том, что для каждого значения индекса суммирования k (то есть для каждой гармони­ки с частотой k1 ) в формуле фигурируют два слагаемых — синус и косинус. Воспользовавшись формулами тригонометрических преобразований, сумму этих двух слагаемых можно трансформировать в косинус той же частоты с иной ам­плитудой и некоторой начальной фазой:



(1)

Если s(t) является четной функцией, фазы φk могут принимать только значе­ния 0 и π, а если s(t) — функция нечетная, то возможные значения для фазы рав­ны ±π/2.



Комплексная форма

Данная форма представления ряда Фурье является, пожалуй, наиболее употре­бимой в радиотехнике. Она получается из вещественной формы представлением косинуса в виде полусуммы комплексных экспонент (такое представление выте­кает из формулы Эйлера ):




Применив данное преобразование к вещественной форме ряда Фурье, получим суммы комплексных экспонент с положительными и отрицательными показате­лями:

А теперь будем трактовать экспоненты со знаком «минус» в показателе как чле­ны ряда с отрицательными номерами. В рамках этого же общего подхода посто­янное слагаемое а0/2 станет членом ряда с нулевым номером. В результате получится комплексная форма записи ряда Фурье:



(2)

Комплексные коэффициенты ряда связаны с амплитудами Ak и фазами φk, фигу­рирующими в вещественной форме записи ряда Фурье (1), следующими не­сложными соотношениями:



, , .

Несложно выглядят и формулы связи с коэффициентами аk и bk синусно-косинусной формы ряда Фурье



, , .

Отсюда сразу же следует и формула непосредственного расчета коэффициентов Ck ряда Фурье в комплексной форме:



(3)

Если s(t) является четной функцией, коэффициенты ряда Ck будут чисто веще­ственными, а если s(t) — функция нечетная, коэффициенты ряда окажутся чис­то мнимыми. Совокупность амплитуд гармоник ряда Фурье часто называют амплитудным спектром, а совокупность их фаз — фазовым спектром.



Преобразование Фурье

Преобразование Фурье является инструментом спектрально­го анализа непериодических сигналов. Для наглядной иллюстрации перехода от ряда Фурье к преобразованию Фурье часто используется не вполне строгий математически, но зато понятный подход. Представим себе периодическую последовательность импульсов произвольно­го вида и сформируем ряд Фурье для нее. Затем, не меняя формы одиночных импульсов, увеличим период их повторения (заполнив промежутки нулевым значением) и снова рассчитаем коэффициенты ряда Фурье. Формула (3) для расчета коэффициентов ряда показывает, что нам придется вычислить тот же самый интеграл, но для более тесно расположенных частот k = k1. Изменение пределов интегрирования не играет роли — ведь на добавившемся между им­пульсами пространстве сигнал имеет нулевое значение. Единственное дополни­тельное изменение будет состоять в уменьшении общего уровня гармоник из-за деления результата интегрирования на увеличившийся период Т.

На рис. 4 описанные изменения иллюстрируются на примере двукратного увеличения периода следования прямоугольных импульсов. Обратите внимание на то, что горизонтальная ось спектральных графиков проградуирована в значени­ях частот, а не номеров гармоник.




Рис. 4. Изменение спектра последовательности импульсов при двукратном увеличении периода их следования

Итак, с ростом периода следования импульсов гармоники располагаются ближе друг к другу по частоте, а общий уровень спектральных составляющих становит­ся все меньше. При этом вид вычисляемого интеграла (3) не меняется.

Наконец, если устремить период к бесконечности (превратив тем самым перио­дическую последовательность в одиночный импульс), гармоники спектра будут плотно занимать всю частотную ось, а их амплитуды упадут до нуля (станут бес­конечно малыми). Однако взаимное соотношение между уровнями гармоник ос­тается неизменным и определяется все тем же интегралом (3). Поэтому при спектральном анализе непериодических сигналов формула для расчета коэффи­циентов комплексного ряда Фурье модифицируется следующим образом:

1. частота перестает быть дискретно меняющейся и становится непрерывным параметром преобразования, (то есть k1 в формуле (3) заменяется на );

2. удаляется множитель 1/Т;

3. результатом вычислений вместо нумерованных коэффициентов ряда Сk явля­ется функция частоты S() — спектральная функция сигнала s(t). Иногда ее называют также спектральной плотностью.

В результате перечисленных модификаций формула (3) превращается в фор­мулу прямого преобразования Фурье:

(3а)

В формуле самого ряда Фурье суммирование заменяется интегрированием (перед интегралом появляется деление на 2). Получающееся выражение называется обратным преобразование Фурье:



(3б)

Дискретное преобразование Фурье
При дискретизации аналогового сигнала его спектр становиться периодическим с периодом повторения, равным частоте дискретизации. Рассмотрим теперь, что представляет собой спектр дискретного периодического сигнала. Итак, пусть последовательность отсчетов x(k) является периодической с перио­дом N:

х (k + N) = х (k) для любого k.

Такая последовательность полностью описывается конечным набором чисел, в ка­честве которого можно взять произвольный фрагмент длиной N, например x(k) { k = 0, 1,..., N-1}. Традиционным представлением является поставленный в соответствие этой последовательности сигнал из смещенных по времени дельта - функций:

(4)

также, разумеется, будет периодическим с минимальным периодом N Т.

где

Так как сигнал (4) является дискретным, его спектр должен быть периоди­ческим с периодом 2/Т. Так как этот сигнал является также и периодическим, его спектр, должен быть дискретным с расстоянием между гармониками, рав­ным 2/(NT).

Итак, периодический дискретный сигнал имеет периодический дискретный спектр, который также описывается конечным набором из N чисел.

Рассмотрим процедуру вычисления спектра периодического дискретного сигна­ла. Так как сигнал периодический, будем раскладывать его в ряд Фурье. Коэф­фициенты Х(n) этого ряда, согласно общей формуле (3), равны







. (5)

Таким образом, формула для вычисления комплексных амплитуд гармоник пред­ставляет собой линейную комбинацию отсчетов сигнала.

В выражении (5) реальный масштаб времени фигурирует только в множите­ле 1/Т перед оператором суммирования. При рассмотрении дискретных последо­вательностей обычно оперируют номерами отсчетов и спектральных гармоник без привязки к действительному масштабу времени и частоты. Поэтому множитель 1/Т из (5) удаляют, то есть считают частоту дискретизации равной единице. Удаляют обычно и множитель 1/N. Получившиеся выражение называется дискретным преобразованием Фурье (ДПФ):
(6)

Существует и обратное дискретное преобразование Фурье. Переход от дискретного спектра к временным отсчетам сигнала выражается следующей формулой:


(7)

Это выражение отличается от формулы прямого ДПФ (6) лишь знаком в показателе комплексной экспоненты и наличием множителя 1/N перед оператором суммирования.



Практическое задание №1.

Создать программу на языке СИ++, реализующую алгоритм дискретного преобразования Фурье (6). Осуществить обработку данной программой различных импульсов и сигналов, представленных в приложении 1 с целью получения их амплитудных и фазовых спектров.



  1. Программа должна выводить амплитудный и фазовый спектр произвольного импульса (сигнала) на экран ПК или в текстовый файл.

  2. Выводить амплитудный и фазовый спектр на экран ПК в виде графиков, используя графические возможности языка СИ++.

  3. Разработать программу расчета амплитудного спектра по формуле (3а) с использованием численных методов (например метод Симсона )

Примечание.

Сигнал х(k) представлять в виде набора N отчетов следующих с одинаковым шагом по времени Т. Таким образом суммарная длительность сигнала составит NT. Для программы в формуле (6) необходимо разделить действительную и мнимую части и получить следующие выражения (такое представление выте­кает из формулы Эйлера ):





где n = 0, 1, 2….., N-1

Амплитудный спектр сигнала получается в виде набора отсчетов на частотной оси по формуле:

для n = 0, 1, 2….., N-1

Фазовый спектр соответственно:



для n = 0, 1, 2….., N-1

Пример представления аналогового прямоугольного сигнала длительностью Т1=NТ на временной оси в дискретном виде :




Т

0 Т1 0 1 ….. N-1



Практическое задание №2.

Составить аналогичную программу на языке СИ++ для обратного дискретного преобразования Фурье (7).



Алгоритм быстрого преобразования Фурье

Для вычисления одного коэффициента ДПФ по формуле (6) необходимо вы­полнить N комплексных умножений и сложений. Таким образом, расчет всего ДПФ, содержащего N коэффициентов, потребует N2 пар операций «умножение-сложение». Число операций возрастает пропорционально квадрату размерности ДПФ. Однако, если N не является простым числом и может быть разложено на множители, процесс вычислений можно ускорить, разделив анализируемый набор отсчетов на части, вычислив их ДПФ и объединив результаты. Такие способы вы­числения ДПФ называются быстрым преобразованием Фурье (БПФ) и повсеместно используются на практике.

При реализации БПФ возможно несколько вариантов организации вычислении в зависимости от способа деления последовательности отсчетов на части (проре­живание по времени либо по частоте) и от того, на сколько фрагментов произво­дится разбиение последовательности на каждом шаге (основание БПФ).

БПФ с прореживанием по времени

Рассмотрим идею БПФ с прореживанием по времени на примере деления набора отсчетов пополам.

Итак, пусть N — четное число. Выделим в формуле (6) два слагаемых, соответ­ствующих элементам исходной последовательности с четными и нечетными но­мерами:

(8)

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


(9)

Две суммы в (9) представляют собой ДПФ последовательностей {у(т)} (от­счеты с четными номерами) и {z(m)} (отсчеты с нечетными номерами). Каждое из этих ДПФ имеет размерность N/2. Таким образом,



(10)

где и — ДПФ соответственно последовательностей отсчетов с четными и нечетными номерами.

Так как ДПФ размерности N/2 дает лишь N/2 спектральных коэффициентов, непосредственно использовать формулу (10) можно только при 0 < n < N/2. Для остальных n (N/2 < n < N) следует воспользоваться периодичностью спектра дискретного сигнала (и, соответственно, периодичностью результатов ДПФ):

,

Следовательно при n ≥ N/2 формула (10) представляется в виде:



(11)

Процесс вычисления 8-точечного ДПФ путем разбиения его на два 4-точечных ДПФ иллюстрируется на рис. 5.



Рис. 5. Вычисление 8-точечного ДПФ с помощью двух 4-точечных ДПФ

Блоки, выполняющие на рис. 5. объединение результатов двух ДПФ, требуют дополнительных комментариев. Каждый такой блок имеет два входных и два выходных сигнала. Один из входных сигналов умножается на комплексную экспоненту , после чего суммируется со вторым входным сигналом и вычитается из него, формируя таким образом два выходных сигнала. Это соответствует реа­лизации формул (10) и (11). Данная операция получила название «бабочки». Расшифровка ее структуры представлена на рис. 6.

Рис. 6. Условное обозначение «бабочки» БПФ с прореживанием по времени (слева) и ее структурная схема (справа)

Оценим количество операций, необходимое для вычисления ДПФ указанным спо­собом. Каждое из двух ДПФ половинной размерности требует N2/4 операций. Кроме того, при вычислении окончательных результатов каждый спектральным коэффициент Z(n) умножается на экспоненциальный комплексный множитель. Это добавляет еще N/2 операций. Итого получается 2N2/4 + N/2 = N(N + l)/2, что почти вдвое меньше, чем при вычислении ДПФ прямым способом.

Если N/2 тоже является четным числом (то есть если N делится на 4), можно продолжить описанную процедуру, выразив результат через четыре ДПФ раз­мерности N/4. Это позволяет еще больше сократить число требуемых вычисли­тельных операций.

Наибольшая степень ускорения вычислений может быть достигнута при N = 2k, в этом случае деление последовательностей на две части можно продолжать до тех пор, пока не получатся двухэлементные последовательности, ДПФ которых рассчитывается вообще без использования операций умножения (достаточно вычислить сумму и разность двух отсчетов). Число требуемых при этом пар опе­раций «умножение — сложение» можно оценить как Nlog2(N). Таким образом, вычислительные затраты по сравнению с непосредственным использованием фор­мулы (6) уменьшаются в N/log2(N) раз.

Практическое задание №3.

Сделать программу на языке СИ++, реализующую алгоритм вычисления 8-точечного ДПФ с помощью двух 4-точечных ДПФ изображенный на рис. 5. Осуществить обработку данной программой различных импульсов и сигналов, представленных в приложении 1 с целью получения их амплитудных и фазовых спектров.



  1. Программа должна выводить амплитудный и фазовый спектр произвольного импульса (сигнала) на экран ПК или в текстовый файл.

  2. Выводить амплитудный и фазовый спектр на экран ПК в виде графиков, используя графические возможности языка СИ++.


БПФ с прореживанием по частоте

Формулы прямого и обратного ДПФ (6) и (7) отличаются только знаком в показателе экспоненты и множителем перед суммой. Поэтому можно получить еще один вариант алгоритма БПФ, выполнив преобразования, показанные на схеме рис. 5, в обратном порядке. Этот способ вычислений называется прореживанием по частоте. Покажем, как получить описание этого метода на основе формулы прямого ДПФ (6).

Разделим исходную последовательность {x(k)} на две следующие друг за другом половины (как и в предыдущем случае, N должно быть четным числом):


Из второй суммы можно выделить множитель



Этот множитель равен 1 или - 1 в зависимости от четности номера вычисляемо­го спектрального отсчета n, поэтому дальше рассматриваем четные и нечетные n по отдельности. После выделения множителя ±1 комплексные экспоненты в обе­их суммах становятся одинаковыми, поэтому выносим их за скобки, объединяя, две суммы:



, (12)

. (13)

Фигурирующие здесь суммы представляют собой ДПФ суммы и разности поло­вин исходной последовательности, при этом разность перед вычислением ДПФ умножается на комплексные экспоненты ехр(-j2m/N). Каждое из двух исполь­зуемых здесь ДПФ имеет размерность N/2.

Итак, при прореживании по частоте вычисления организуются следующим об­разом:

1. Из исходной последовательности {x(k)} длиной N получаются две последова­тельности {у(m)} и {z(m)} длиной N/2 согласно следующим формулам:



,

.

2. ДПФ последовательности {у(m)} дает спектральные отсчеты с четными номерами, ДПФ последовательности {z(m)} - с нечетными:



,

.

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

Процесс вычисления 8-точечного ДПФ путем разбиения его на два 4-точечных ДПФ с прореживанием по частоте показан на рис. 7.


Рис. 7. Вычисление 8-точечного ДПФ с помощью двух 4-точечных ДПФ путем прореживания по частоте

Поскольку комплексный экспоненциальный множитель в данном алгоритме применяется к результату вычитания двух сигналов, «бабочка» БПФ с прорежива­нием по частоте имеет несколько иную структурную схему (рис. 8.).


Рис. 8. Условное обозначение «бабочки» БПФ с прореживанием по частоте (слева и ее структурная схема (справа)



Практическое задание №4.

Разработать программу на языке СИ++, реализующую алгоритм вычисления 8-точечного ДПФ с помощью двух 4-точечных ДПФ путем прореживания по частоте изображенный на рис. 7. Осуществить обработку данной программой различных импульсов и сигналов, представленных в приложении 1 с целью получения их амплитудных и фазовых спектров. Отсчеты амплитудных и фазовых спектров выводить в файл или на экран ПК в текстовом и графическом режиме.



Корреляционный анализ

Корреляционный анализ наряду со спектральным играет большую роль в теории сигналов. Говоря кратко, его смысл состоит в количественном измерении степе­ни сходства различных сигналов. Для этого служат корреляционные функции.



Корреляционная функция

Корреляционная функция (КФ) детерминированного сигнала с конечной энергией представляет собой интеграл (в бесконечных пределах) от произведения двух копий сигнала, сдвинутых друг относительно друга на время τ:



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

1. Значение КФ при τ = 0 равно энергии сигнала, то есть интегралу от его квад­рата:

2. КФ является четной функцией своего аргумента τ:

3. Значение КФ при τ = 0 является максимально возможным значением:

4. С ростом абсолютного значения τ КФ сигнала с конечной энергией затухает:


5. Если сигнал s(t) не содержит особенностей в виде дельта - функций, его КФ не


может иметь разрывов, то есть обязана быть непрерывной функцией.

График КФ прямоугольного импульса показан на рис. 9.





Рис. 9. Корреляционная функция прямоугольного импульса

В случае периодического сигнала (и вообще любого сигнала с бесконечной энер­гией) воспользоваться приведенным определением не удастся. Поэтому КФ пе­риодического сигнала с периодом Т вычисляют, усредняя произведение сдвину­тых копий в пределах одного периода:



Набор свойств такой КФ несколько меняется.

1. Значение при τ = 0 равно не энергии, а средней мощности анализируемого сигнала.

2. Свойство четности сохраняется.

3. Значение КФ при τ = 0 по-прежнему является максимально возможным:
4. КФ периодического сигнала является периодической функцией с тем же периодом, что и сам сигнал:

5. Если сигнал не содержит дельта - функций, его КФ будет непрерывной функцией.



Взаимная корреляционная функция

Если КФ показывает степень сходства между сдвинутыми копиями одного и того же сигнала, то взаимная корреляционная функция (ВКФ) позволяет измерить аналогичную величи­ну для сдвинутых экземпляров двух разных сигналов.

Общий вид формулы КФ сохраняется, но под интегралом стоит произведение двух разных сигналов, один из которых задержан на время τ:

График ВКФ прямоугольного и треугольного импульса показан на рис. 10.


Рис. 10. ВКФ прямоугольного и треугольного импульсов



Автокорреляционная функция стационарной случайной

последовательности отсчетов

Итак, пусть есть последовательность отсчетов {x(k)}, где каждый отсчет является случайной величиной. Величину x(k) можно представить как x(nТ), где Т – временной интервал между отсчетами.

Если известна реализация случайной последовательности, состоящая из N- отсчетов, то оценкой математического ожидания является выборочное среднее:

Оценка автокорреляционной функции определяется как



(14)

Автокорреляционная функция служит мерой корреляции между отсчетами случайной последовательности.



Практическое задание №5.

Создать программу на языке СИ++, по расчету оценки автокорреляционной функции с помощью формулы (14) для произвольной последовательности отсчетов. Определить данной программой автокорреляцию различных импульсов и сигналов, представленных в приложении 1. Полученные, в результате расчетов отсчеты автокорреляционной функции выводить в файл или на экран ПК в текстовом и графическом режиме.



Линейная свертка

Линейная свертка двух бесконечных дискретных последовательностей {x1(k)} и {x2(k)} опре­деляется следующим образом:



Для двух конечных последовательностей по N1 и N2 отсчетов соответственно линейная свертка будет определяться соотношением:



, где k=0,…., N1+N2-2 (15)

Круговая свертка
Если мы рассматриваем периодические последовательности, то и суммирова­ние при вычислении свертки таких последовательностей следует производить по одному периоду. Такую операцию называют круговой сверткой'.

В этой формуле выражение означает взятие по модулю N, то есть вы­числение остатка от деления на N.



Практическое задание №6.

Разработать программу на языке СИ++, по расчету линейной свертки с помощью формулы (15) для двух разных конечных последовательностей по N1 и N2 отсчетов. В качестве последовательностей взять импульсы и сигналы различной формы из приложения 1, представленные дискретным рядом, то есть в виде последовательности отсчетов. Отсчеты линейной свертки выводить в файл или на экран ПК в текстовом и графическом режиме.


Контрольные вопросы.

  1. Что такое частота Найквиста? Какая необходима частота дискретизации исходного сигнала, чтобы однозначно восстановить его по отсчетам дискретизированного сигнала?

  2. Что такое АЦП и ЦАП?

  3. Какие существуют формы представления ряда Фурье?

  4. Дать определение непрерывного прямого и обратного преобразования Фурье.

  5. Записать и объяснить формулу для прямого дискретного преобразования Фурье (ДПФ).

  6. Как определяется амплитудный и фазовый спектр дискретизированного сигнала?

  7. Пояснить схему работы быстрого дискретного преобразования Фурье с прореживанием по частоте?

  8. Пояснить схему работы быстрого дискретного преобразования Фурье с прореживанием по времени?

  9. Дать понятие автокорреляционной функции сигнала. Какие свойства у корреляционной функции сигнала?

  10. Поясните физический смысл взаимной корреляционной функции.

  11. Что такое линейная и круговая свертка сигнала?

ЛИТЕРАТУРА

1. Сергиенко А. Б. Цифровая обработка сигналов — СПб.: Питер, 2003, 608с.

2. Гольденберг Л.М., Матюшкин Б.Д., Поляк М.Н. Цифровая обработка сигналов — М.: Радио и связь, 1985, 312 с.

3. Голд Б., Рэйдер Ч. Цифровая обработка сигналов / Пер. с англ.; Под ред.


А. М. Трахтмана. — М.: Сов. радио, 1973, 368 с.

4. Рабинер Л., Гоулд Б. Теория и применение цифровой обработки сигналов


Пер. с англ.; Под ред. Ю. И. Александрова. — М.: Мир, 1978.

5. Марпл-мл. С. Л. Цифровой спектральный анализ и его приложения / Пер. с


англ. - М.: Мир, 1990, 584с.

Приложение 1.

Примечание.

Все импульсы и сигналы различной формы должны быть представлены дискретным рядом, то есть в виде последовательности отсчетов:




  1. Прямоугольный импульс:





  1. Прямоугольный импульс, задержанный во времени:







  1. Несимметричный треугольный импульс:





4. Симметричный треугольный импульс:




  1. Односторонний экспоненциальный импульс:







  1. Двухсторонний экспоненциальный импульс:







  1. Гауссов импульс:


8. Сигнал вида sin(х)/х:





.

9. Сигнал вида sin(х):


10. Сигнал вида cos(х):


11. Меандр – последовательность прямоугольных импульсов со скважностью равной двум:


n = 0, 1, 2…. N.

12. Пилообразный сигнал:

13. Последовательность треугольных импульсов:


14 Гауссов радиоимпульс:






База данных защищена авторским правом ©bezogr.ru 2016
обратиться к администрации

    Главная страница