135865 (Цифровая обработка сигналов), страница 2

После замены переменной в (1.5) получим

X(jkw₁) = x(nT) ,

x(nT) = X(jkw₁) .

Отсюда после подстановки w₁ = w_д/N, T = 2p/w_д формулы ДПФ принимают окончательный вид

X(jkw₁) = x(nT) - прямое ДПФ ,

x(nT) = X(jkw₁) - обратное ДПФ (1.10)

Сигнал с ограниченным спектром имеет, строго говоря, бесконечную протяженность во времени и, соответственно бесконечное число отсчетов и непрерывный спектр. Спектр останется непрерывным, если число отсчетов сигнала ограничить конечным числом N. Формулы (1.10) в этом случае будут выражать связь между N отсчетами дискретного сигнала и N отсчетами его непрерывного спектра, который можно полностью восстановить по его отсчетам.

Пример. Определить отсчеты спектра сигнала на Рис. 1.5, а.

Здесь N = 2 поэтому X(jkw₁) = x(nT) e^-jpkn следовательно

X(j0w₁) = x(nT)e^-j0 = x(0T) + x(1T) = a + b

X(j1w₁) = x(nT)e^-jpn = x(0T) e^-j0 + x(1T) e^-jp = a - b

график отсчетов спектра приведен на Рис. 1.5, б, где w₁ = w_д/N = 0,5w_д.

Сигнал с конечным числом отсчетов N имеет спектр, который повторяет с конечной погрешностью спектр сигнала с бесконечным числом отсчетов : спектры совпадают на отсчетных частотах kw₁ и отличаются на других частотах. Отличие спектров тем меньше, чем больше N. В самом деле, реальные сигналы обладают конечной энергией и, следовательно, начиная с некоторого номера отсчета остальными номерами можно пренебречь ввиду их малости, что не окажет заметного влияния на спектр сигнала.

Пример. Осуществить дискретизацию экспоненциального импульса X(t) = Ae^-at = 1 e^-10t и сравнить спектры исходного и дискретного сигналов.

Решение.

График сигнала X(t) представлен на Рис. 1.8

Пусть T = 0,02с. В этом случае плавным соединением отсчетов сигнала (штриховая линия на Рис. 1.8) сигнал восстанавливается удовлетворительно хотя заметны искажения в окрестности точки t = 0, поэтому ошибки наложения будут некоторым образом влиять на спектральные характеристики.

Пусть t_u = 0,4с. В этом случае

N = t_u/T = 20.

Расчет спектра по формуле прямого ДПФ в точке w = 0 (k = 0) запишется так

X(j0w₁) = 1,0 + 0,8187 + 0,6703 + 05488 + 0,4493 + 0,368 + 0,3012 + 0,2466 + 0,2019 + 0,1653 + 0,1353 + 0,1108 + 0,09072 + 0,07427 + 0,06081 + 0,04979 + 0,04076 + 0,03337 + 0,02732 + 0,02237 = 5,41

Истинное значение спектра в точке w = 0 можно определить зная спектр аналогового экспоненциального импульса

X_a(jw) = , следовательно X_a(j0) = = 0,1.

чтобы сравнить спектры дискретного и непрерывного сигналов, дискретный спектр необходимо денормировать умножением на T, так как формулы Фурье для дискретных сигналов применяются в нормированном виде. Поэтому

X(jow₁) = 5,41 T = 5,42Ч0,02 = 0,1082.

Таким образом совпадение спектров X_a(jw) и X(jw) в точке w = 0 вполне удовлетворительное. Некоторая неточность объясняется влиянием ошибок наложения.

Уместно заметить, что выбор шага дискретизации достаточно контролировать в точках максимальной крутизны исходной функции X(t). В рассмотренном примере такой точкой является момент времени t = 0.

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

1. Дискретные цепи.

1. Разностное уравнение и дискретная цепь.

Непрерывный сигнал на входе линейной системы x(t) и соответствующий сигнал y(t) на выходе связаны дифференциальным уравнением. Замена непрерывной переменной t на дискретную переменную nT приводит к замене дифференциального уравнения разностным уравнением. Каноническая форма разностного уравнения общего вида, учитывающая в явном виде наличие в системе как прямых, так и обратных связей, запишется так

y(nT) = a_m x(nT - mT) + y(nT - ), (2.1)

где (M + 1) - число прямых связей,

Z - число обратных связей,

m, , n - целые положительные числа.

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

Система работа которой описывается разностными уравнениями, является дискретной так как она способна воздействовать только на отсчеты сигнала. Дискретная система и дискретная цепь осуществляет, согласно (2.1) следующие операции над дискретными сигналами.

1. Сдвиг (запаздывание) на целое число интервалов T

1. Умножение на некоторый коэффициент a_m или b

1. Сложение сигналов.

Перечисленные операции образуют полный базис, в котором можно реализовать заданное воздействие на сигнал.

Набору операций базиса соответствует набор типов элементов дискретной цепи : элементы памяти (задержки), умножители и сумматоры.

Каноническая схема дискретной цепи общего вида, соответствующая разностному уравнению (2.1), приведена на Рис. 2.1.

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

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

Пример. Составить разностное уравнение цепи, схема которой приведена на Рис. 2.2, а.

Решение.

Здесь имеется три пути прохождения сигнала от входа до выхода цепи, по которым сигналы проходят и затем складываются в сумматоре. Поэтому разностное уравнение имеет вид

y(nT) = 0,5 x(nT) - 0,7 x(nT - T) + 0,35 x(nT - 2T).

Пример. Определить y(nT) (Рис. 2.2, б), если x(nT) = {1,0 ; 0,5}.

Решение.

Разностное уравнение цепи y(nT) = 0,5 x(nT - T) + 0,1 x(nT) численное решение разностного уравнения :

n=0; y (0T) = 0,5 x(-T) + 0,1 x(0T) = 0,1;

n=1; y (1T) = 0,5 x(0T) + 0,1 x(1T) = 0,55;

n=2; y (2T) = 0,5 x(1T) + 0,1 x(2T) = 0,25;

n=3; y (3T) = 0,5 x(2T) + 0,1 x(3T) = 0.

Следовательно y(nT) = {0,1; 0,55; 0,25}.

Графики сигналов x(nT) и y(nT) приведены на рис (2.3,а,б).

Пример. Определить сигнал на выходе цепи (рис 2.2,в), если y(nT)={0,1; 0,1}.

Решение.

Цепь содержит обратную связь (ОС), поэтому сигнал на выходе цепи формируется как сигнал со стороны входа, так и со стороны выхода.

y(nT) = 0,4 x(nT-T) - 0,08 y(nT-T)

n=0 y(0T) = 0,4 x(-T) - 0,08 y(-T) = 0

n=1 y(1T) = 0,4 x(0T) - 0,08 y(0T) = 0,4

n=2 y(0T) = 0,4 x(1T) - 0,08 y(1T) = 0,368 и т.д. ...

Следовательно y(nT) = {0; 0,4; 0,368; ...}.

В данном случае за счет циркуляции сигнала по цепи ОС выходной сигнал состоит из бесконечного числа отсчетов.

Дискретная цепь, содержащая ОС, называется рекурсивной. Дискретная цепь без ОС называется нерекурсивной.

2.2 Передаточная функция дискретной цепи.

Замена сигналов в разностном уравнении (2.1) на Z - изображения этих сигналов

,

приводит к алгебраизации разностного уравнения

.

Алгебраизация осуществляется применением теорем линейности и запаздывания.

Переход в область Z - изображений позволяет ввести понятие передаточной функции дискретной цепи H(Z), которая определяется как отношение Z - изображения сигнала на выходе цепи к Z - изображению сигнала на входе цепи. Поэтому, учитывая алгебраическую форму разностного уравнения общего вида, можно записать общий вид передаточной функции дискретной цепи

. (2.3)

Отсюда, в частности, для нерекурсивной цепи

. (2.4)

Если нерекурсивная цепь состоит всего из одного элемента запаздывания, то ,

что находит своё отражение в обозначении элементов памяти на схемах дискретных цепей.

Передаточная функция конкретной цепи формируется по передаточным функциям её элементов согласно общих правил линейных цепей. В частности, для цепи содержащей ОС применяется известная формула

, (2.5)

где - передаточная функция цепи

прямого прохождения сигнала,

- предаточная функция цепи ОС.

Пример. Оперделить передаточную функцию цепи на рис. (2.4,а).

Решение.

, где , .

Пример. Определить передаточную функцию на рис.(2.4,б).

Решение.

,

где - передаточная функция рекурсивной части схемы,

- передаточная функция нерекурсивной части цепи.

По известной передаточной функции можно легко определить разностное уравнение цепи.

Пример. Составить разностное уравнение цепи на рис.(2.2,в).

Решение.

Здесь .

Поэтому .

Отсюда .

Следовательно переходя к оригиналам: y(nT)= 0,4 x(nT-T) - 0,08 y(nT-T).

2.3 Общие свойства передаточной функции.

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

z = x + jy.