Просмотр исходного текста страницы Моделирование работы модулярного КИХ фильтра, исследование эффективности экстраполяционных методов для коррекции ошибок

Цель настоящей работы состояла в том, чтобы промоделировать модулярный КИХ фильтр с постоянными коэффициентами, базируясь на идее "Теоремы о Свертке" и ее аналога в конечном поле Галуа. Моделировался фильтр низких частот с опцией обнаружения одиночной ошибки, а также процесс экстраполяции результата в случае обнаружения ошибки.

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

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

<center><math>s(n)=a*b=\sum_{m=0}^{n}a(m)*b(n-m)</math>,     <math>n=0</math> ... <math>N+M-2</math></center>

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

Последовательная схема характеризуется малым числом вычислительных блоков, интенсивным обменом с памятью и низкой производительностью. В крайнем проявлении эта схема представляет собой умножитель с накоплением и управляющее устройство, которое обеспечивает загрузку нужных коэффициентов из памяти. В этом случае, для нахождения одного выходного отсчета требуется <math>N</math> тактов. Этот метод реализуется программным методом на сигнальных процессорах или компьютерах общего назначения.

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

[[Файл:800px-FIR Filter.png]]

Преимущества данной архитектуры - это ее быстродействие и возможность работы в реальном времени. К минусам можно отнести значительное увеличение аппаратурных затрат.

Кроме реализаций во временной области, возможна также реализация в частотной области. Базисом для этого является так называемая '''Теорема о Свертке'''. Спектр циклической свертки есть произведение спектров сворачиваемых сигналов: <math>S(k)=A(k)*B(k)</math>. Где <math>A(k)</math> и <math>B(k)</math> - спектры сворачиваемых сигналов, <math>S(k)</math> - спектр циклической свертки двух сигналов.

<center> [[Файл:Conv Theorem.JPG]] </center>

Таким образом, вместо того чтобы реализовывать свертку непосредственно по формуле, можно перевести сигналы в частотную область с помощью БПФ, там почленно перемножить, и перевести обратно. Результат будет соответствовать циклической свертке [статья в ростове]. Дополняя нулями последовательности, вычисление циклической свертки станет эквивалентно вычислению искомой линейной свертки. Преимуществом этого метода является сокращение операций. Действительно, используя быстрые схемы преобразования Фурье, мы тратим лишь <math>n*lg(n)</math> операций, что дает выигрых при больших значениях <math>n</math>. 

Однако, существуют и недостатки. Зачастую в фильтрах подразумевается входной сигнал бесконечной длительности, который фильтруется некоторым ограниченным набором коэффициентов. Параллельная схема отлично справляется с этой задачей в "онлайн" режиме. С методами на основе БПФ дело обстоит сложнее. Во первых не возможно создать БПФ на бесконечное число точек, а во вторых, даже если мы создадим очень большой преобразователь, то мы потеряем "онлайн" режим, так как для того чтобы БПФ начал выдавать значения, нам необходимо чтобы все значения были туда загружены. Чтобы решить данную проблему были введены различные "overlap" методики. Они позволяют обрабатывая блоки малого размера, совмещать их таким образом, чтобы на выходе мы получили необходимую свертку теоретически бесконечной последовательности отсчетов [та презенташка с оверлап методиками]. 

Резюмируя вышесказанное, что касается методов фильтрации в частотной области, можно отметить что в принципе этот метод можно также реализовывать как в последовательном так и в параллельном виде. Для параллельного метода можно использовать различные конвейерные методики реализации БПФ (типа Radix-2 SPDF). Последовательные методы реализуются программно. Судя по всему, чаще всего эта методика используется именно в последовательной реализации в виде программ на микроконтроллерах и компьютерах общего назначения. Условно классифицировать методы реализации КИХ фильтров можно следующим образом:


<center> [[Файл:Методы.JPG]] </center>

Для реализации фильтров, а в общем и других устройств ЦОС, эффективно использовать аппарат модулярной арифметики. Можно использовать модулярную арифметику при реализации любого КИХ фильтра, реализованного любым методом. Однако, наибольший интерес представляет использование модулярной арифметики при параллельной реализации в частотной области. Существует несколько предпосылок для этого вывода, основная из которых заключается в том, что для конечных полей существует аналог теоремы о свертке. Точнее говоря, для конечных полей существует преобразование, аналогичное преобразованию Фурье для которого выполняется теорема о свертке в конечном поле, и для которого соблюдаются все симметрические свойства позволяющие использовать любой быстрый алгоритм для его вычисления. Формула упомянутого преобразования выглядит следующим образом:

<center><math> A_k=|\sum_{n=0}^{2^t-1}|a_n*\gamma^{nk}|_p|_p </math></center>

Здесь <math>A_k</math> - теоретико-числовой спектр сигнала, <math>a_n</math> - сам сигнал, <math>p</math> - характеристика конечного поля, <math>\gamma</math> - примитивный корень степени <math>2^t</math> [Виноградов, Теория чисел]. Если число <math>p</math> имеет вид <math>p=q*2^t+1</math>, то примитивный корень степени <math>2^t</math> - всегда существует и равен <math>\gamma = \omega^q</math>.

Таким образом, упомянутое преобразованиек является аналогом дискретного преобразования Фурье тригонометрического базиса, в котором поворачивающие множители <math>e^{ij\frac{2\pi}{N}}</math> заменяются на  <math>\gamma^{ij}</math>, а операции умножения и сложения комплексных чисел – на соответствующие операции в <math>GF(p)</math>. 

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

== Моделирование работы модулярного КИХ фильтра ==
=== Краткое описание эксперимента ===
Была создана модель модулярного КИХ фильтра нижних частот порядка 16. КИХ фильтр реализовывался в модулярном виде по принципу теоремы о свертке и методами блочной обработки(overlap-add). В качестве базисного набора оснований использовался набор:

97=6∙16+1
     
193= 12∙16+1 

257=16∙16+1

В качестве дополнительного модуля для контроля ошибок использовался модуль:

353=22∙16+1

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

=== Исходные данные ===

За основу был взят произвольный сигнал на 840 точек.

[[Файл:Firmatlabtest1.JPG]] [[Файл:Firmatlabtest2.JPG]]

<p style="text-align:right;">''Основной сигнал	'''и''' Сигнал после воздействия внешних факторов''</p>

=== Проверка влияния масштабирующего множителя ===

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

Увеличение коэффициентов    
* позволяют уменьшить ошибки представления
* увеличивают аппаратуру и уменьшают производительность схемы

Уменьшение коэффициентов        
* позволяют выбирать меньшие основания, соответственно ускоряя и уменьшая схему
* результат становиться менее точным

[[Файл:Firmatlabtest3.JPG]] [[Файл:Firmatlabtest4.JPG]]

<p style="text-align:right;">''Сигналы сглаженные фильтром с высоким(2048) '''и''' низким(16) разрешением ''</p>

[[Файл:Firmatlabtest31.JPG]] [[Файл:Firmatlabtest41.JPG]]

<p style="text-align:right;">''Сигналы сглаженные фильтром с высоким(2048) '''и''' низким(16) разрешением (укрупненный) ''</p>

=== Проверка характеристик обнаружения ошибки и уменьшение её влияния интерполяционными методами ===

Следующим этапом была проверка характеристик обнаружения ошибки и уменьшение её влияния интерполяционными методами. В процессе разработки модели был получен важный вывод: характер искажения выходных данных существенно зависит от места возникновения ошибки. В том случае, если ошибка возникает в любом из каналов на этапах от прямого преобразователя до обратного БПФ включительно, на выходе искажается серия данных, длиной равной длине блока обработки. В случае возникновения ошибки в конце тракта обработки – искажается одно единственное значение.

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

[[Файл:Firmatlabtest5.JPG]] [[Файл:Firmatlabtest6.JPG]]

<p style="text-align:right;">''График функции с ошибкой в конце канала обработки с вероятностью возникновения 12,5% '''и''' График, сглаженный фильтром с интерполяционным сглаживанием возможных сбоев.''</p>

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

[[Файл:Firmatlabtest7.JPG]] [[Файл:Firmatlabtest8.JPG]]

<p style="text-align:right;">''График функции с ошибкой в начале канала обработки с вероятностью возникновения 0,25% '''и''' График, сглаженный фильтром с интерполяционным сглаживанием возможных сбоев.''</p>