Цифровая обработка сигналов (стр. 1 из 11)
ВВЕДЕНИЕ
В ЦИФРОВУЮ ОБРАБОТКУ
СИГНАЛОВ
Содержание.
1. Дискретные сигналы
1.1. Дискретизация непрерывных сигналов
1.2. Связь спектров дискретных и непрерывных сигналов
1.3. Преобразование Фурье и Лапласа для дискретных сигналов
1.4. Z - преобразование
1.5. Основные теоремы Z - преобразования
1.6. Дискретное преобразование Фурье
2. Дискретные цепи
2.1. Разностное уравнение и дискретная цепь
2.2. Передаточная функция дискретной цепи
2.3. Общие свойства передаточной функции
2.4. Частотные характеристики
2.5. Импульсная характеристика. Свертка.
2.6. Круговая свертка
2.7. Энергия дискретного сигнала. Корреляция и энергетический спектр
2.8. Расчет энергии сигнала в дискретной цепи
2.9. Секционирование
3. Цифровые фильтры
3.1. Цифровая система обработки сигналов
3.2. Расчет не рекурсивных ЦФ общего вида
3.3. Схема и характеристики фильтров с линейной фазой
3.4. Общие свойства фильтров с линейной фазой
3.5. Расчет ЦФ с линейной фазой. Метод взвешивания.
3.6. Метод частотной выборки
3.7. Расчет рекурсивных фильтров. Метод билинейного преобразования
4. Эффекты конечной разрядности и их учет.
4.1. Шум квантования и шумовая модель
4.2. Расчет шумов квантования
4.3. Влияние структуры ЦФ на шум квантования
4.4. Квантование коэффициентов. Расчет разрядности.
4.5. Чувствительность
4.6. Масштабирование сигнала в цепи
4.7. Динамический диапазон ЦФ
4.8. Предельные циклы
5. Восстановление непрерывного сигнала
5.1. Характеристики ЦАП
5.2. Погрешности восстановления
Литература
Обсуждены основные положения теории дискретных сигналов и способы их обработки. Рассмотрены особенности цифровой реализации дискретных систем. Изложены методы расчета цифровых фильтров, получившие наибольшее распространение.
Эффекты конечной разрядности ЦФ и их учет рассмотрены применительно к системам с фиксированной запятой. Погрешности дискретизации и восстановления обсуждены на уровне необходимом для понимания вопроса.
Для технических факультетов.
1. Дискретные сигналы.
1.1 Дискретизация непрерывных сигналов.
Обработка сигналов на цифровых ЭВМ начинается с замены непрерывного сигнала X(t) на дискретную последовательность, для которой применяются такие обозначения
x(nT) , x(n) , xn , {x0 ; x1 ; x2 ; … } .
Дискретизация осуществляется электронным ключом (ЭК) через равные интервалы времени T (Рис. 1.1).
Дискретная последовательность аппроксимирует исходный сигнал X(t) в виде решетчатой функции X(nT). Частота переключения электронного ключа fд и шаг дискретизации T связаны формулой
fд = 1 / T . (1.1)
Дискретная последовательность или дискретный сигнал выражается через исходный непрерывный (аналоговый) сигнал следующим образом
x(nT) = x(t)
d(t - nT) , (1.2)где d(t) - дискретная d - функция (Рис. 1.2, а),
d(t - nT) - последовательность d - функций (Рис. 1.2, б). 
Погрешность, возникающую при замене аналогового сигнала дискретным сигналом, удобно оценить сравнивая спектры этих сигналов.
1.2. Связь спектров дискретного и непрерывного сигналов.
Исходное выражение для спектра дискретного сигнала с учетом (1.2) запишется следующим образом
X(jw) =
x(nT) e-jwt dt = x(t) 
d(t - nT) e-jwt dt .Периодическую последовательность d - функций здесь можно разложить в ряд Фурье
d(t - nT) = 
,где с учетом формулы связи спектров периодического и непериодического сигналов
, поскольку Fd(jw) = 1После замены в исходном выражении периодической последовательности d - функций ее разложением в ряд Фурье получим
X(jw) =
x(t)( 
) e-jwt dt = 
x(t) 
e-jwt dt .Учитывая здесь теорему смещения спектров, т.е. :
если f(t) ® F(jw), то f(t)
® F[j(w±w0)] ,последнее равенство можно представить в виде формулы, выражающей связь спектров дискретного X(jw) и аналогового Xa(jw) сигналов
X(jw) =
Xa[j(w - 
)] . (1.3)На основании формулы (1.3) с учетом поясняющих рисунков 1.3, а, б можно сделать следующие выводы :
1. Спектр дискретного сигнала состоит из суммы спектров исходного непрерывного сигнала, сдвинутых друг относительно друга по оси частот на величину равную частоте дискретизации wд
2. Спектры аналогового и дискретного сигналов совпадают в диапазоне частот [-0,5wд ; 0,5wд], если удовлетворяется неравенство
wвЈ 0,5wд , (1.4)
где wв - верхняя частота спектра аналогового сигнала.
Равенство в (1.4) соответствует утверждению теоремы Котельникова о минимальной частоте wд.
1. Смежные спектры Xa(jw) в (1.3) частично перекрываются, если условие (1.4) не выполняется (Рис 1.3, б). В этом случае спектр дискретного сигнала искажается по отношению к спектру аналогового сигнала. Эти искажения являются неустранимыми и называются ошибками наложения.
2. Аналоговый сигнал можно восстановить полностью по дискретному сигналу с помощью ФНЧ, частота среза которого wс = 0,5wд. Это утверждение основано но совпадении спектров дискретного сигнала на выходе ФНЧ и непрерывного сигнала. Сигнал восстанавливается без искажений, если выполняется условие (1.4). в противном случае сигнал восстанавливается с искажениями, обусловленными ошибками наложения.
Выбор частоты дискретизации осуществляется в соответствии с (1.4). если частота wв не известна, то выбор из wд определяется расчетом по формуле (1.1), в которой интервал T выбирается приближенно с таким расчетом, чтобы аналоговый сигнал восстанавливался без заметных искажений плавным соединением отсчетов дискретного сигнала.
1.3 Преобразование Фурье и Лапласа для дискретных сигналов.
Для дискретных сигналов формулы Фурье и Лапласа представляется возможным упростить. Действительно, поскольку
то после перехода к дискретной переменной пара преобразований Фурье принимает вид
Здесь применяются формулы одностороннего преобразования Фурье, так как начало отсчета совмещается с началом действия дискретного сигнала.
Формулы Фурье для дискретных сигналов применяются в нормированном виде, поэтому после замены X(nT) ® X(nT) / T преобразование Фурье принимает окончательный вид
(1.5)Формулы Лапласа для дискретных сигналов получаются на основании (1.5) после обобщения частоты на всю плоскость комплексного переменного, то есть jw® P = d + jw
(1.6)1.4. Z - преобразование.
Эффективность частотного анализа дискретных сигналов существенно возрастает, если заменить преобразование Лапласа Z - преобразованием. В этом случае изображение сигнала X(p), которое представляет собой трансцендентную функцию переменной P = d + jw, заменяется Z - изображением сигнала X(Z), которое является рациональной функцией переменной Z = x + jy.
Формулы Z - преобразования получаются из формулы Лапласа (1.6) заменой переменных
epT = Z . (1.7)
Подстановка (1.7) и ее производной
dZ / dp = TepT
в (1.6) приводит к формулам прямого и обратного Z - преобразования
(1.8)Точки на мнимой оси комплексного переменного p = d +jw, то есть точки p = jw, определяют реально частотные характеристики сигнала. Мнимой оси соответствует на плоскости Z единичная окружность, так как в этом случае согласно (1.7)
Z = ejwT =
(1.9)Поэтому непрерывному росту переменной на мнимой оси плоскости p = d + jw, соответствует многократный обход единичной окружности на плоскости z = x + jy (Рис. 1.4). Этим фактом объясняется, в частности, то обстоятельство, что интегрирование в формуле обратного z - преобразования (1.8) осуществляется вдоль единичной окружности плоскости z взамен интегрирования вдоль прямой параллельной мнимой плоскости p.