91
6 ЦИФРОВАЯ ФИЛЬТРАЦИЯ
6.1 Исходные понятия и определения
Дискретный сигнал. Под дискретным сигналом понимают дискретную последовательность значений реального сигнала, поступающих в определенные фиксированные моменты времени ti .
Дискретные системы. Системы, работающие с дискретным сигналом, называют дискретными.
Источники сигнала могут быть дискретны и непрерывны по своей природе возникновения. В случае непрерывного сигнала, его преобразование в дискретную форму производится с помощью стробирования и выборки значений в фиксированные моменты времени. Из соображений удобства обычно используют дискретные сигналы с равными промежутками времени между отсчетами. Интервал времени между последовательными отсчетами называется периодом последовательности T . Величина обратная периоду дискретизации называется частотой дискретизации и определяется
соотношением d |
2 |
/T . |
Дискретную последовательность или сигнал удобно описывать в виде |
||
{x(t0 ), x(t0 T ), x(t0 |
2 |
T ), , x(t0 k T )} или {x0, x1, x2, , xk }. Значение |
отсчета соответствует амплитуде реального сигнала в соответствующий момент времени.
Напомним, что оператор сдвига последовательности yk на один такт вперед (упреждение) записывается в виде E yk yk 1 , а оператор сдвига на
один такт назад (задержка) описывается соотношением E |
1 y |
y |
1 |
. В |
|
k |
k |
|
области изображений действие операторов упреждения и задержки соответствует умножению функции изображения Yz , соответственно, на z и
z 1. Символически это соответствие можно отобразить соотношениями
E y |
y |
z Y ; E 1 |
y |
y |
z |
1 Y . |
k |
k 1 |
z |
|
k k 1 |
|
z |
Часто для представления значения используют дискретные значения амплитуды для удобства дальнейшей обработки. Операцию фиксации уровня называют дискретизацией по амплитуде или квантованием уровней и выполняют ее соответствующими устройствами – компараторами. Число уровней квантования существенно влияет на точность обработки сигнала. Отношение максимально возможного уровня к минимальному уровню, называют динамическим диапазоном.
Часто, чтобы снизить требования к рабочему диапазону амплитуд, значение отсчета представляется набором нулей и единиц двоичного кода, определенной разрядности, с определенной тактовой частотой, позволяющей передать двоичное представление отсчета до поступления, следующего отсчета. При этом повышаются требования к быстродействию устройств обработки сигнала в виде двоичной последовательности импульсов.
|
92 |
|
|
|
Преобразование |
непрерывного |
сигнала |
в |
двоичную |
последовательность, путем дискретизации по времени и амплитуде выполняется аналого-цифровыми преобразователями (АЦП). Обратное преобразование сигналов из двоичной последовательности в непрерывный путем восстановления отсчетов и аппроксимации значений сигнала между отсчетами выполняется цифро-аналоговыми преобразователями (ЦАП).
Цифровые сигналы. Дискретные сигналы, представленные в двоичном коде, называются цифровыми.
Цифровые системы. Под цифровыми системами понимают разновидность дискретных систем обрабатывающих цифровую последовательность в виде двоичных m - разрядных последовательностей импульсов соответствующих дискретным отсчетам входного сигнала.
Дискретные, а особенно цифровые сигналы удобны для обработки с помощью импульсных логических схем, процессоров и ЭЦВМ.
В настоящее время для обработки цифровых сигналов используют цифровые фильтры, которые могут быть реализованы аппаратно и/или программно.
Цифровые фильтры. Цифровые фильтры (ЦФ), это класс систем (устройств) дискретной обработки цифровых сигналов, так называемые линейные стационарные цифровые фильтры. В общем случае, под цифровыми фильтрами будем понимать аппаратную и/или программную реализацию цифровой обработки сигналов.
Цифровые фильтры выполняют операцию частотной фильтрации, преобразуя цифровой поток входного сигнала в выходной сигнал, используя взвешенное суммирование текущих и предыдущих отсчетов сигнала. Современные ЦФ обладают высокой стабильностью параметров, возможностью реализации любой формы АЧХ и ФЧХ, не требуют подстройки и легко реализуются на микропроцессорах и ЭЦВМ. В низкочастотном диапазоне цифровые фильтры способны реализовать полосы пропускания либо заграждения и избирательность, не достижимые другими типами фильтров.
На рисунке 6.1 приведена функциональная схема цифровой обработки аналоговых сигналов.
Цифровая обработка производится следующим образом. Непрерывный входной сигнал x(t) поступает на вход АЦП и, в момент подачи
синхроимпульса от тактового генератора (ТГ), его значение преобразуется в двоичную последовательность импульсов определенной разрядности. Двоичная последовательность фиксированной разрядности передается последовательным либо параллельным кодом на арифметическое устройство (АУ) при поступлении на соответствующие шины уровней напряжений соответствующих разрядам двоичного числа. Далее цифровой процессор АУ, выполняющий циклические операции умножения и сложения по модулю два и операции сдвига разрядов, обрабатывает оцифрованные значения отсчетов входного сигнала. Оперативное запоминающее устройство (ОЗУ) служит для
93
запоминания необходимого числа предшествующих значений входного и выходного сигналов необходимых для реализации алгоритма цифровой фильтрации. С выхода цифрового процессора АУ цифровое представление выходного сигнала при необходимости поступает на ЦАП для перевода в аналоговую форму, в противном случае передается на выход в цифровой форме.
Быстродействие, как основной показатель ЦФ зависит как от сложности алгоритма, разрядности двоичного представления, так и от быстродействия процессора. Частотный диапазон сигналов обрабатываемых ЦФ постоянно увеличивается с ростом тактовой частоты современных микропроцессоров.
Квантование уровней. Как уже отмечалось, значение входного сигнала дискретизируется по амплитуде, то есть квантуется по уровням, таким образом, чтобы минимальное значение сигнала соответствовало младшему разряду двоичного представления
|
N |
2k , |
x |
ik |
|
i |
|
|
k |
0 |
|
где ik 0, 1- весовые коэффициенты, принимающие значения базового
набора двоичного представления; N - разрядность двоичного представления. Максимально представимое целое десятичное число или число уровней сигнала, представимое N разрядным двоичным числом, определяется
формулой M 2N 1. При конечном числе разрядов погрешность, связанная с округлением значения до представимого уровня.
Сигналы, принимающие счетное значение уровней называются квантованными.
Таким образом, квантование сигналов приводит к специфической погрешности их представления называемого шумом квантования. Для снижения шумов квантования следует увеличивать разрядность цифровых процессоров, однако это снижает их быстродействие и скорость цифровой обработки.
94
На практике, в зависимости от решаемых задач, обычно используют 4, 8, 16, 32– х разрядные АЦП, при реализации устройств цифровой обработки сигналов.
6.2 Алгоритм цифровой фильтрации
Алгоритм цифровой фильтрации. Известно, что стационарная аналоговая линейная система преобразует входной сигнал, в соответствии с выражением
y(t) x( ) g(t )d x(t ) g( )d ,
где g(t) - импульсная характеристика стационарной линейной системы.
Как отмечалось, линейный ЦФ представляет собой дискретную систему либо программу, преобразующие входную последовательность числовых отсчетов {xk } в выходную последовательность отсчетов выходного
сигнала {yk }
{x0, x1, x2, } {y0, y1, y2, } ; {xk } {yk }.
Из свойства линейности ЦФ следует, что сумма любого числа входных сигналов умноженных на коэффициент преобразуется в сумму выходных сигналов, умноженных на тот же коэффициент
m |
m |
{x(i)} |
{y(i)} |
k |
k |
i 1 |
i 1 |
Для того, чтобы воспользоваться исходным соотношением для |
|
преобразования входного сигнала в выходной для цифрового фильтра, по аналогии с непрерывными системами вводят понятие импульсной характеристики gk .
Напомним, что импульсной характеристикой дискретной системы
или ЦФ, находящихся в исходном состоянии покоя, называется его
реакция на одиночный единичный - импульс 10 |
{1, 0, 0, 0, } |
10 {1, 0, 0, 0, } {g0, g1, g2, g3, |
} {gk }. |
Под состоянием покоя понимается полное установление реакции на предыдущее воздействие и отсутствие сторонних источников (возмущений).
Линейный ЦФ стационарен, то есть, при смещении входного единичного импульса на любое число интервалов дискретизации по времени, его импульсная характеристика смещается таким же образом
10 {gk }; 10 m {gk m}.
Таким образом, алгоритм цифровой фильтрации при условии линейности и стационарности ЦФ по аналогии с непрерывными системами, можно записать в виде
|
95 |
|
|
|
m |
ym x0 gm x1 gm 1 |
xm g0 |
xk gm k . |
|
k |
0 |
Данное соотношение представляет собой свертку входного сигнала и импульсной характеристики ЦФ. Из выражения следует, что каждый отсчет выходного сигнала представляет собой сумму произведений предшествующих отсчетов входного сигнала на соответствующие отсчеты импульсной характеристики ЦФ.
Таким образом, ЦФ обладает памятью по отношениям к прошлым воздействиям.
Физически реализуемые ЦФ, исходя из принципа причинности, обладают тем свойством, что сигнал на выходе ЦФ равен нулю до подачи входного сигнала. Это означает, что отсчеты импульсной характеристики {g 1, g 2, g 3, } равны нулю и суммирование следует распространять
только на положительные значения индекса k . Это позволяет записать основное соотношение для образования выходного сигнала в виде
|
m |
ym |
xk gm k , |
k |
0 |
где m 0,1, 2, .
Дискретные гармонические последовательности. Как известно, в
теории непрерывных систем особая роль отводится комплексному представлению гармонических сигналов вида
x(t) A e j ( t ) .
Дискретизация по времени такого сигнала приводит к дискретной гармонической последовательности
{x } {A e j ( k T |
)} |
|
k |
|
|
такой, что |
|
|
Re{xk } {A cos( |
k T |
)} , |
где T - период дискретизации по времени; |
- начальный фазовый сдвиг. |
|
Отметим, что эти последовательности представляют дискретизированные гармонические сигналы неоднозначно, так как последовательности не изменятся при замене
2 n |
n |
d , |
|
|
|||
T |
|||
|
|
где n - любое целое число; d - угловая частота дискретизации.
Это означает, что отсчеты сигнала можно выбрать либо за один период гармонического сигнала, либо выбирать по одному из каждого k - го периода, где k 1, 2, 3, . Данное свойство дискретного представления периодических
последовательностей соответствует так называемому стробоскопическому эффекту.
Частотная характеристика ЦФ. Пусть на вход ЦФ поступает неограниченная во времени гармоническая последовательность {xk } , где
| 05_Холера |
| 1 |
| 10.4. Исследование регистров |
| 1112 |
| 12 |
| 1285 |
| 1560 |
| 1568 |
| 1604248853606027 |
| 1673 |