Пример реализации

В качестве примера использования DSP сравним аналоговый и цифровой фильтры низкой частоты (ФНЧ), каждый с частотой среза 1 кГц.

Цифровой фильтр реализован в виде типовой цифровой системы, показанной на рис. 1.4. Обратите внимание, что в диаграмме принято несколько неявных допущений. Во -первых, чтобы точно обработать сигнал, принимается, что тракт АЦП/ЦАП обладает достаточными значениями частоты дискретизации, разрешающей способности и динамического диапазона. Во -вторых, для того, чтобы закончить все свои вычисления в пределах интервала дискретизации (1/f_s), устройство ЦОС должно иметь достаточное быстродействие. В -третьих, на входе АЦП и выходе ЦАП сохраняется потребность в аналоговых фильтрах ограничения и восстановления спектра сигнала (anti-aliasing filter и anti -imaging filter), хотя требования к их производительности невелики. Приняв эти допущения, можно сравнить цифровой и аналоговый фильтры.

Рис. 13.4 Структурная схема цифрового фильтра.

Требуемая частота среза обоих фильтров - 1 кГц. Аналоговое преобразование реализуется фильтром Чебышева первого рода шестого порядка (характеризуется наличием пульсаций коэффициента передачи в полосе пропускания и отсутствием пульсаций вне полосы пропускания). Его характеристики представлены на рис. 1.5. На практике этот фильтр может быть представлен тремя фильтрами второго порядка, каждый из которых построен на операционном усилителе и нескольких резисторах и конденсаторах. С помощью современных систем автоматизированного проектирования (САПР) фильтров создать фильтр шестого порядка достаточно просто, но чтобы удовлетворить техническим требованиям по неравномерности характеристики 0,5 дБ, требуется точный подбор компонентов.

Представленный же на рис 1.4 цифровой КИХ-фильтр со 129 коэффициентами имеет неравномерность характеристики всего 0,002 дБ в полосе пропускания, линейную фазовую характеристику и намного более крутой спад. На практике такие характеристики невозможно реализовать с использованием аналоговых методов. Другое очевидное преимущество схемы состоит в том, что цифровой фильтр не требует подбора компонентов и не подвержен дрейфу параметров, так как частота синхронизации фильтра стабилизирована кварцевым резонатором. Фильтр со 129 коэффициентами требует 129 операций умножения с накоплением (MAC) для вычисления выходного отсчёта. Эти вычисления должны быть закончены в пределах интервала дискретизации 1/fs, чтобы обеспечить работу в реальном масштабе времени. В этом примере частота дискретизации равна 10 кГц, поэтому для обработки достаточно 100 мкс, если не требуется производить существенных дополнительных вычислений. Семейство DSP ADSP-21xx может закончить весь процесс умножения с накоплением (и другие функции, необходимые для реализации фильтра) за один командный цикл. Поэтому фильтр со 129 коэффициентами требует быстродействия более 129/100 мкс = 1,3 миллиона операций с секунду (MIPS). Существующие DSP имеют намного большее быстродействие и, таким образом, не являются ограничивающим фактором для этих приложений. Быстродействие серии 16-разрядных ADSP-218x с фиксированной точкой достигает 75MIPS. На рис. 1.6 приведен ассемблерный код, реализующий фильтр на DSP процессорах семейства ADSP-21xx. Обратите внимание, что фактические строки исполняемого кода помечены стрелками; остальное - это комментарии.

Рис. 13.5 Амплитудно-частотные характеристики фильтров

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

Программа на ассемблере: fir фильтр для adsp-21xx

.MODULE fir_sub;

{ Подпрограмма КИХ фильтра

Параметры вызова подпрограммы

I0 --> Наиболее старые данные в линии задержки

I4 --> Начало таблицы коэффициентов фильтра

L0 = Длина фильтра (N)

L4 = Длина фильтра (N)

M1,M5 = 1

CNTR = Длина фильтра - 1 (N-1)

Возвращаемые значения

MR1 = Результат суммирования (округлённый и ограниченный)

I0 --> Наиболее старые данные в линии задержки

I4 --> Начало таблицы коэффициентов фильтра

Изменяемые регистры

MX0,MY0,MR

Время работы

(N - 1) + 6 cycles = N + 5 cycles

Все коэффициенты записаны в формате 1.15 }

.ENTRY fir;

fir: MR=0, MX0=DM(I0,M1), MY0=PM(I4,M5)

CNTR = N-1;

DO convolution UNTIL CE;

convolution: MR=MR+MX0*MY0(SS), MX0=DM(I0,M1), MY0=PM(I4,M5);

MR=MR+MX0*MY0(RND);

IF MV SAT MR;

RTS;

.ENDMOD;