Расчёт и анализ нерекурсивного цифрового фильтра
Общая характеристика цифрового фильтра, его описание разностным уравнением. Импульсная характеристика и комплексная частотная характеристика, их связь парой преобразований Фурье. Метод частотной выборки и наименьших квадратов, их сравнение и отличия.
| Рубрика | Коммуникации, связь, цифровые приборы и радиоэлектроника |
| Вид | курсовая работа |
| Язык | русский |
| Дата добавления | 22.02.2011 |
| Размер файла | 138,0 K |
Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже
Студенты, аспиранты, молодые ученые, использующие базу знаний в своей учебе и работе, будут вам очень благодарны.
Размещено на http://www.allbest.ru/
1. Краткое математическое описание методов расчёта
1.1. Общие положения
Цифровой фильтр полностью описывается своим разностным уравнением:
(1)
Для нерекурсивного цифрового фильтра и уравнение принимает вид:
(2)
Зная коэффициенты разностного уравнения, можно легко получить выражение для передаточной функции фильтра (для НЦФ):
(3)
Для образа выходного сигнала НЦФ справедливо выражение
, (4)
где - z-преобразования выходного и входного сигналов фильтра.
Зная выражение (4) и учитывая, что z-преобразование функции единичного скачка равно 1, можно получить выражение для z-образа импульсной характеристики :
(5)
Из (5) следует, что отсчеты импульсной характеристики НЦФ численно равны коэффициентам разностного уравнения НЦФ, а сама импульсная характеристика и передаточная функция связаны парой z-преобразований (прямым и обратным).
Заменив в (4) z на , получим комплексную частотную характеристику:
(6)
Импульсная характеристика и комплексная частотная характеристика связаны парой преобразований Фурье:
(7)
(8)
Из комплексной частотной характеристики можно получить выражения для АЧХ и ФЧХ:
(9)
(10)
Во все вышеприведённые формулы входит интервал квантования . Чтобы от него избавиться, частоту обычно нормируют. Это можно сделать с помощью замены:
(11)
Так как интервал определения , то интервал определения . Исходными данными для проектирования фильтра является его АЧХ. Как правило, в зонах неопределённости АЧХ некоторым образом доопределяют с тем, чтобы избежать явления Гиббса («выбросы» характеристики в точках разрыва первого рода - «скачках»). В простейшем случае доопределить АЧХ можно линейным законом. В этом случае АЧХ проектируемого полосового фильтра будет выглядеть таким образом.
Аналитически АЧХ будет записываться в виде:
(12)
При проектировании часто полагают, что ФЧХ фильтра является линейной. В [1] показывается, что в этом случае импульсная характеристика фильтра является либо симметричной (), либо антисимметричной (). Учитывая, что порядок фильтра может быть чётным и нечётным, существует четыре вида ИХ с линейной ФЧХ:
1. N - нечётное, ИХ - симметричная
2. N - чётное, ИХ - симметричная
3. N - нечётное, ИХ - антисимметричная
4. N - чётное, ИХ - антисимметричная
цифровой фильтр выборка частотный
1.2 Метод частотной выборки
Основная идея метода частотной выборки - замену в выражениях (7) и (8) непрерывную частоту дискретизированной. В этом случае выражения (7) и (8) превращаются в пару дискретных преобразований Фурье:
(13)
(14)
Существует 2 метода дискретизации частоты (выражения записаны для нормированной частоты):
(15)
(16)
Выражения (13) и (14) записаны для первого метода дискретизации частоты. По условию задания необходимо использовать второй метод дискретизации частоты, в этом случае выражение (14) приобретает вид:
(17)
Из (17) следует, что для определения импульсной характеристики необходимо знать частотную характеристику. Её можно записать в показательной форме:
(18)
(19)
При чётном N:
(20)
При нечётном N:
(21)
Подставляя вместо , по выражениям (20) и (21) можно найти , а из (17) - .
1.3 Метод наименьших квадратов
При расчете коэффициентов импульсной характеристики используется формула вида:
после чего решается система уравнений:
и находятся коэффициенты Ск.
Далее из найденных Ск можно найти коэффициенты импульсной характеристики:
2. Расчётная часть
2.1 Расчёт методом частотной выборки
2.1.1 Расчёт импульсной характеристики
Расчёт импульсной характеристики для нечётных N осуществлялся по формулам (21) и (17), для чётных - по формулам (20) и (17). Результаты расчёта импульсной характеристики для N=15, 25 и 32 представлены в таблице 1.
Таблица 1. Результаты расчёта импульсной характеристики методом частотной выборки
|
i |
Значение импульсной характеристики |
|||
|
N=15 |
N=25 |
N=32 |
||
|
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 |
0,081 -0,013 0,025 -0,052 -0,303 0,03 0,46 0,03 -0,303 -0,052 0,025 -0,013 0,081 |
0,001497 0,001756 -0,02 -0,007456 -0,007554 0,028 0,061 -0,004905 0,034 -0,048 -0,297 -0,035 0,45 0,035 -0,297 -0,048 0,034 -0,004905 0,061 0,028 -0,007454 -0,007456 -0,02 0,001756 0,001497 |
0,001488 -0,008534 0,008698 -0,000256 0,003711 -0,011 0,015 -0,007875 -0,001266 0,053 0,029 0,0009025 0,04 -0,193 -0,224 0,321 0,321 -0,224 -0,193 0,04 0,0009025 0,029 0,053 0,001266 -0,007875 -0,015 -0,011 -0,003711 -0,000256 0,008698 -0,0008534 0,001488 |
2.1.2 Расчёт АЧХ и ФЧХ
Расчёт АЧХ и ФЧХ осуществлялся по формулам (9) и (10) для 50 значений частоты , взятой с шагом 0,01 (). На рисунках приведены графики рассчитанной АЧХ фильтра.
Для расчёта точности аппроксимации запишем функцию ошибки аппроксимации:
, (32)
В таблице 2 приведены результаты расчёта точности аппроксимации .
Таблица 2. Результаты расчета точности аппроксимации для метода частотной выборки
График функции точности аппроксимации для N=25
Максимальные ошибки аппроксимации (абсолютная погрешность) для трёх значений N приведены в таблице 3:
Абсолютная погрешность аппроксимации АЧХ, рассчитанной методом частотной выборки
|
Абсолютная погрешность аппроксимации АЧХ |
|||
|
N=13 |
N=25 |
N=32 |
|
|
0,125 |
0,082 |
0,049 |
2.2 Расчёт методом наименьших квадратов
2.2.1 Расчёт импульсной характеристики
Результаты расчёта импульсной характеристики для N=13, 25 и 32 представлены в таблице. Учитывая симметрию импульсной характеристики, приведена только половина отсчётов.
Результаты расчёта импульсной характеристики методом наименьших квадратов
|
i |
Значение импульсной характеристики |
|||
|
N=13 |
N=25 |
N=32 |
||
|
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 |
0,055 -0,004049 0,035 -0,042 -0,296 0,03 0,45 |
-0,003929 -0,003499 -0,012 0,008469 -0,008832 -0,026 0,055 0,035 -0,042 -0,296 0,03 0,45 |
0,002208 -0,005211 0,003349 0,003189 -0,003929 -0,003499 -0,012 -0,008469 -0,008832 0,026 0,055 -0,004049 0,035 -0,042 -0,296 0,45 0,45 |
2.2.2 Расчёт АЧХ и ФЧХ
Расчёт АЧХ и ФЧХ осуществлялся по формулам (9) и (10) для 50 значений частоты , взятой с шагом 0,01 ().
Заданная по условию и рассчитанная АЧХ фильтра для N=25 (метод наименьших квадратов)
2.2.3 Расчёт точности аппроксимации
Точность аппроксимации оценивалась по формуле (32). В таблице (5) приведены результаты расчёта
Результаты расчета точности аппроксимации для метода наименьших квадратов
В таблице 6 приведена максимальная (абсолютная) погрешность аппроксимации для различных значений N.
Абсолютная погрешность аппроксимации для метода наименьших квадратов
|
Абсолютная погрешность аппроксимации АЧХ |
|||
|
N=135 |
N=25 |
N=32 |
|
|
0,125 |
0,057 |
0,051 |
2.3 Сравнение методов расчёта
Сравнивая результаты расчётов точности аппроксимации, приведённые в таблицах 2 и 6, можно сделать вывод, что метод наименьших квадратов обеспечивает более точную аппроксимацию при N=25 амплитудно-частотной характеристики по сравнению с методом частотной выборки. С увеличением порядка фильтра N точность аппроксимации увеличивается для обоих методов, но точность метода наименьших квадратов начинает уменьшаться по сравнению с методом частотной выборки.
Заключение
В данной курсовой работе был рассмотрен расчёт нерекурсивного цифрового фильтра двумя методами: методом наименьших квадратов и методом частотной выборки. Результаты расчётов точности аппроксимации для каждого метода позволяют сделать следующие выводы:
· Точность аппроксимации увеличивается с увеличением N (порядка фильтра)
· Метод наименьших квадратов обеспечивает более точную аппроксимацию при средних значениях N.
Размещено на Allbest.ru
Подобные документы
Испытание синтезированного нерекурсивного и рекурсивного цифрового фильтра стандартными и гармоническими сигналами. Расчет реакции фильтра на четырехточечный входной сигнал. Получение системной функции и частотных характеристик цифрового фильтра.
курсовая работа [3,0 M], добавлен 19.05.2015Расчет цифрового и аналогового фильтра-прототипа. Структурные схемы и реализационные характеристики фильтра. Синтез цифрового фильтра в системе программирования MATLAB. Частотные и импульсные характеристики цифрового фильтра, карта его нулей и полюсов.
курсовая работа [564,8 K], добавлен 24.10.2012Алгоритм расчета фильтра во временной и частотной областях при помощи быстрого дискретного преобразования Фурье (БПФ) и обратного быстрого преобразования Фурье (ОБПФ). Расчет выходного сигнала и мощности собственных шумов синтезируемого фильтра.
курсовая работа [679,2 K], добавлен 26.12.2011Расчет характеристик фильтра во временной и частотной областях с помощью быстрого преобразования Фурье, выходного сигнала во временной и частотной областях с помощью обратного быстрого преобразования Фурье; определение мощности собственных шумов фильтра.
курсовая работа [2,8 M], добавлен 28.10.2011Алгоритм, реализующий заданный тип фильтра в частотной области. Спектр входного, выходного сигнала. Спектральная (амплитудно-частотная) характеристика окна. Отклик фильтра на заданный сигнал. Двусторонний экспоненциальный радиоимпульс с несущей частотой.
курсовая работа [318,2 K], добавлен 07.07.2009Определение импульсной характеристики фильтра. Расчет амплитудно- и фазово-частотной характеристик и методами разложения в ряд Фурье, наименьших квадратов и частотной выборки. Построение графиков и оценка точности аппроксимации (абсолютной погрешности).
курсовая работа [677,0 K], добавлен 21.12.2012Представление чисел в дополнительном двоичном коде. Номенклатура арифметических операций в цифровом фильтре. Назначение аналого-цифрового преобразователя. Амплитудно-частотная характеристика рекурсивного фильтра. Составление схемы лабораторного макета.
реферат [215,6 K], добавлен 14.02.2016Структурная схема цифрового фильтра. Расчет устойчивости, построение графиков. Виды свертки дискретных сигналов. Определение выходного сигнала в частотной области с помощью алгоритма "бабочка". Схема шумовой модели фильтра, мощность собственных шумов.
курсовая работа [641,3 K], добавлен 15.10.2013Построение графиков амплитудного и фазового спектров периодического сигнала. Расчет рекурсивного цифрового фильтра, цифрового спектра сигнала с помощью дискретного преобразования Фурье. Оценка спектральной плотности мощности входного и выходного сигнала.
контрольная работа [434,7 K], добавлен 10.05.2013Разработка и описание общего алгоритма функционирования цифрового режекторного фильтра на основе микропроцессорной системы. Обоснование аппаратной части устройства. Отладка программы на языке команд микропроцессора. Расчёт быстродействия и устойчивости.
курсовая работа [266,1 K], добавлен 03.12.2010
