Просмтор этого файла доступен только зарегистрированным пользователям. Но у нас супер быстрая регистрация: достаточно только электронной почты!

Текст из файла (страница 4)

Проведите исследование для различных структур, в том числе и для каскадной. Укажите, для какой структуры разрядность коэффициентов наименьшая.

- 25 -

Если при уменьшении формата каждый раз на единицу характеристики меняются незначительно, используйте более эффективные процедуры, например, метод дихотомии. Так после формата [16 15] задайте формат [8 7]. Если АЧХ укладывается в поле допусков, перейдите к формату [4 3], а если не укладывается, к формату [12 11] , и т.д., пока не будет найден формат с минимально допустимым количеством двоичных разрядов.

При малой разрядности коэффициентов возможна потеря устойчивости фильтра, вследствие того что модули некоторых полюсов из-за округления могут превысить единицу (в особенности тех из них, которые, обладая высокой добротностью, располагаются близко к единичной окружности). Информация об устойчивости или неустойчивости фильтра заносится в раздел Current Filter Information основного окна fdatool (Stable: Yes или Stable: No). Проводя исследование эффектов квантования, следует отмечать, сохраняет ли фильтр устойчивость при последовательном уменьшении разрядности его коэффициентов.

9.2. Расчёт наименьшей разрядности входного сигнала и выходных регистров умножителей цифрового фильтра.

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

D = 20 lg (A_max / A_min) [дБ],

где А_max и A_min – максимальная и минимальная амплитуды входного сигнала, и отношение сигнал/шум квантования на выходе фильтра

R = 10 lg (P_с / P_ш) [дБ],

где P_с – минимальная мощность сигнала (в случае синусоидального сигнала P_с=A²_min/2), P_ш=²_вых – допустимая мощность (дисперсия) шума квантования на выходе фильтра.

Примем, что разрядность входного сигнала и разрядность умножителей цифрового фильтра одинаковы. Кроме того, будем полагать, что входной сигнал масштабирован, и его максимальная амплитуда А_max = 1. Масштабированы пусть будут и коэффициенты фильтра, так что максимальное значение коэффициента передачи в полосе пропускания равно единице. В этом случае можно считать, что

P_ш = A²_max/2 10 ^{– (R+D)/10} . (6)

- 26 -

В соответствии с линейной шумовой моделью фильтра на входе фильтра и на выходе каждого умножителя «точно» представленные отсчёты сигналов суммируются с шумом квантования:

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

²=2 ^{– 2p}/12, (7)

где p – количество разрядов сигнала на выходе умножителя (без учёта знакового). На выходе фильтра каждый источник шума квантования создаёт шум с дисперсией ²_{вых i} =² (g_i[n])², где g_i[n] – импульсная характеристика

ⁿ

части фильтра от i-го источника шума до выхода. Суммирование квадратов

отсчётов импульсной характеристики ведётся для всех номеров n, при которых значения g_i[n] существенны (не являются пренебрежимо малыми).В силу независимости источников шума полная дисперсия шума квантования на выходе фильтра равна сумме дисперсий отдельных источников:

²_вых=²_{вых i} . В результате анализа, основанного на изложенном подходе,

ⁱ

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

1. Прямая и транспонированная структуры.

²_вых=²_вх( (g[n])²+(k+m)(g_рек[n])²), (8)

где k – количество умножителей в обратных связях (с коэффициентами a),

m – количество умножителей в прямых связях (с коэффициентами b),

g_рек[n] – импульсная характеристика рекурсивной части фильтра.

- 27 -

В случае прямой структуры, как видно из её схемы (см. подраздел 6.1), шум всех умножителей проходит только через рекурсивную часть (умножители с коэффициентами a), в то время как входной шум проходит через весь фильтр. То же самое можно сказать и о транспонированной структуре (см. подраздел 6.3). Дисперсия шума умножителей равна дисперсии входного шума ²_вх, поскольку, как указывалось выше, отсчёты сигнала везде представлены одинаковым количеством разрядов p. Следует отметить, что числа k и m необязательно равны количеству коэффициентов a_k и b_k , т.е. числам N и M+1 соответственно, так как некоторые из коэффициентов могут быть нулевыми или равняться единице. В этих случаях умножители не применяются.

1. Каноническая структура.

²_вых=²_вх ((k+1)(g[n])² + m). (9)

Анализ канонической структуры (см. подраздел 6.2) показывает, что входной шум и шум умножителей рекурсивной части (коэффициенты a) проходят через весь фильтр, а шум умножителей с коэффициентами b непосредственно проходит на выход.

1. Каскадная структура со звеньями в виде прямых или транспонированных структур.

²_вых=²_вх c₁c₂…c_L , (10)

где c_i= (g_i[n])² + (k_i+m_i)(g_{рек i} [n])² , i=1,2,…, L, (11)

L – количество каскадов.

Шум квантования, прошедший через первый каскад, характеризуется дисперсией ²_вхc₁. Этот шум является входным для следующего каскада, поэтому дисперсия шума на выходе второго каскада ²_вхc₁c₂ и т.д. Из этого рассуждения становится понятным, каким образом составлено выражение (10).

1. Каскадная структура со звеньями в виде канонических структур.

²_вых=²_вх c₁c₂…c_L ,

где c_i= (k_i+1) (g_i[n])² + m_i , i=1,2,…, L (12)

1. Параллельная структура со звеньями в виде прямых или транспонированных структур.

²_вых=²_вх (c₁ +c₂ +…+c_L ) , (13)

- 28 -

где c_i определяются выражением (11), L – количество параллельно включённых звеньев.

1. Параллельная структура со звеньями в виде канонических структур.

²_вых=²_вх (c₁ +c₂ +…+c_L ) ,

где c_i определяются выражением (12), L – количество параллельно включённых звеньев.

1. Нерекурсивный фильтр.

²_вых=²_вх ((g[n])² + m) =²_вх(b_k² + m) (14)

(см. рис. на с.19) .

Допустимую дисперсию (среднюю мощность ) шума квантования на выходе можно рассчитать по формуле (6), в которой положить А_max=1. Затем из формул (8) – (14) выразить дисперсию входного шума квантования ²_вх, предварительно рассчитав отношение дисперсий для нужной структуры в соответствии с приведёнными выражениями. Далее на основании выражения (7) получаем наименьшее количество двоичных разрядов:

p= int [ 0.5 log₂ (1/(12²_вх)) ] +1, (15)

где int [] – операция взятия целой части. С учётом знакового разряда нужно полученное по формуле (15) значение увеличить ещё на единицу.

При работе в среде MatLab для расчёта наименьшей разрядности сигнала и выходных регистров умножителей цифрового фильтра можно применить программу minubit. Она вызывается следующим образом:

>> minubit (b, a, D, R)

Здесь b и a – векторы коэффициентов передаточной функции фильтра, D – динамический диапазон входного сигнала [дБ], R – допустимое отношение сигнал/ шум квантования на выходе фильтра [дБ].

Программа рассчитывает наименьшее количество разрядов (с учётом знакового) для структур перечисленных выше типов. Если фильтр рекурсивный, то производится расчёт для девяти структур (см. с. 26 – 28). Если фильтр нерекурсивный, то для одной структуры (см. с.28 и рис. на с.19). Результаты расчётов выводятся в командное окно по завершении работы программы. Кроме наименьшей разрядности приводятся также значения дисперсии шума квантования на входе и выходе фильтра. Анализируя полученные результаты, можно выбрать оптимальную структуру, обеспечивающую заданное отношение сигнал/ шум квантования на выходе фильтра при заданном динамическом диапазоне и позволяющую установить самую маленькую разрядность отсчётов сигнала по сравнению с другими структурами.

- 29 -

9.3. Расчёт дисперсии шума квантования на выходе фильтра при заданной разрядности отсчётов сигнала

Поставим теперь задачу несколько иначе. Пусть разрядность входного сигнала, а также разрядность сигналов на выходах умножителей известна. Нужно рассчитать дисперсию шума квантования на выходах различных структур цифрового фильтра, обладающих одной и той же передаточной функцией K(z). Поскольку дисперсия входного шума однозначно определяется количеством разрядов (см. (7)), то она одинакова для всех структур, а так как выражения (8) – (14), связывающие дисперсии шума квантования на входе и выходе различны, получается, что разные структуры будут давать на выходе шум квантования различной средней мощности.

Чтобы произвести расчёт, вызовите программу quanod:

>> quanod (b, a, p)

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

10. Моделирование работы цифрового фильтра

Моделирование работы цифрового фильтра предполагает задание тестового сигнала, использование его отсчётов в качестве входных в алгоритме цифровой фильтрации, нахождение выходного сигнала и сравнение его с входным. Кроме того, полезно рассмотреть спектры входного и выходного сигналов и сопоставить их с частотной характеристикой фильтра.

10.1. Задание тестовых сигналов

Данная процедура осуществляется в рабочей области MatLab. Сигнал задаётся в виде вектора, сопоставленного с вектором моментов времени. Перед вводом непосредственно модели сигнала нужно указать частоту дискретизации и сформировать вектор-столбец моментов времени. Например,

>> Fs = 1e3; t=0:1/Fs:1; t=t’;

- 30 -

В данном случае введена частота дискретизации 1кГц. Сигнал будет задан на интервале времени 1с (1001 отсчёт). Последний оператор означает преобразование вектора-строки в вектор-столбец ( ‘ – операция транспонирования матрицы). Не следует забывать ставить точку с запятой в конце каждого оператора, чтобы подавить вывод значений на экран монитора.

Рассмотрим некоторые из возможных сигналов.

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

Характеристики

Список файлов курсовой работы