Читайте также:
|
|
Пример 1
Задано разностное уравнение:
y(n)=a0x(n)+a1x(n-1)+a2x(n-2)+a3x(n-3) (63)
1 Записать передаточную функцию H(z) цепи.
Для уравнения (63) найдем его z-преобразование:
Y(z)=a0z(z)+a1z-1x(z)+a2z-2x(z)+a3z-3x(z) (64)
По определению, z – преобразование функции передачи равно отношению z – преобразования отклика цепи Y(z) к z – преобразованию входного воздействия x(z):
(65)
2 Записать АЧХ дискретной цепи (63), полагая a0=a1=a2=a3=1.
Пусть входная последовательность представляет собой отсчеты комплексной экспоненты:
(66)
Выходная последовательность y(n) равна дискретной свертке входной последовательности x(n) и дискретной импульсной характеристики ЛДС h(n):
= (67)
Здесь H() – комплексная частотная характеристика дискретной цепи.
(68)
Допустим, что входной сигнал – отсчеты гармонической функции амплитуды A, частоты ω0, начальной фазы :
(69)
(70)
(71)
Выходной сигнал y1(n) на x1(n) равен:
(72)
Выходной сигнал y2(n) на x2(n) равен:
(73)
Отсчеты выходной последовательности y(n) равны сумме y1(n) и y2(n)
(74)
Учитывая, что и - комплексно-сопряженные функции, можно записать:
(75)
Пусть - частота Найквиста
- интервал дискретизации, выбранный в соответствии с теоремой Котельникова для сигнала с ограниченным спектром.
(76)
где - круговая частота дискретизации.
Обозначим нормированную частоту
произведение:
(77)
Выразим из (75) Tg и подставим (77):
(78)
Определим комплексную частотную характеристику ЛДС как функцию частоты :
(79)
Для разностного уровня (63) при заданных исходных данных комплексная частотная характеристика дискретной цепи равна
(80)
где h(n) – дискретная импульсная характеристика.
(81)
| |||
Рисунок 3 – График дискретной импульсной
характеристики ЛДС, N =4
(82)
Формула (82) получена на основе вычисления суммы геометрической прогрессии.
| |||
Рисунок 4 – График АЧХ ЛДС (пример 2)
Пример реализации фильтра Баттерворта 3го порядка
Аналоговый прототип Фильтра Баттерворта 3го порядка имеет квадрат модуля комплексного коэффициента передачи:
|H(jΩ)|2 = = H(jΩ)·H(-jΩ) (1)
Произведение операторного коэффициента передачи H(S) на зеркальную функцию H(-S) будет иметь вид:
H(S) H(-S) = (2)
Pk=jΩk=jej , k= 1,2,3,4,5,6
|
Рисунок 1 – Расположение полюсов аналогового фильтра
на комплексной плоскости
Первые три полюса P1, P2, P3 расположены в левой полуплоскости. В соответствии с условиями физической реализуемости, они принадлежат передаточной функции H(p).
Полюса P4, P5, P6 расположены в правой полуплоскости и принадлежат зеркальной функции H(-P).
Следовательно:
H(S)= , (3)
где P1=- ;
P2=-1;
P3=-
Подставим значения, получим формулу для операторного коэффициента передачи в следующем виде:
H(S)= (4)
Используем билинейное преобразование для перехода от операторного коэффициента передачи к Z-преобразованию функции передачи фильтра:
P=2fg (5)
где fg – частота дискретизации.
Пусть частота среза фильтра определяется равенством:
fc=0,2 fg (6)
S=jΩ=j =j =j (7)
S= (8)
на частоте функция передачи будет иметь вид:
H(P)= = (9)
Подставив в (9) (5) и (6) получим:
H(Z)= (10)
X=
После упрощения (10) будет иметь вид:
H(Z)=
Этой функции соответствует структурная схема рисунок 2.
X(k) Y(k)
Рисунок 2 – Структурная схема цифрового фильтра
Баттерворта 3го порядка с fc=02fg
Проверим с помощью пакета MATLAB 6.5, что эта структурная схема действительно соответствует фильтру Баттерворта.
>> в = [ ] % вектор коэффициентов числителя Z - преобразования функции передачи;
>> а = [ ] % вектор коэффициента знаменателя той же функции;
>> f gz(в,а).
Пример:
Фильтр Баттерворта 3го порядка с частотой среза fc=0,2fg или fN= , fc=0,4fN,
где fc – частота среза ФНЧ (граница ПП);
fg = - частота дискретизации;
fN – частота Нейквиста
имеет Z-преобразование передаточной функции, полученное с помощью билинейного преобразования:
H(Z)= =
При расчете по этому оператору по умолчанию используются нормированные значения частот, измеряемые в радианах на отсчет. При такой нормировке частота дискретизации равна 2π, а частота Найквиста (максимальная частота спектра сообщения) равна π. При этом число частотных точек равно 512 на интервале 0/π с постоянным частотным шагом. [с. 218, 1]
Варианты заданий
Вариант 1
Дата добавления: 2015-11-03; просмотров: 57 | Нарушение авторских прав
<== предыдущая страница | | | следующая страница ==> |
Линейной дискретной системы | | | Дискретная цепь описывается разностным уравнением |