3. Свойства симметрии
Если периодическая последовательность xp(n) с периодом в N отсчетов действительная, то ее ДПФ Xp(k) удовлетворяет условиям симметрии:
Re[Xp(k)] = Re[Xp(N – k)], Im[Xp(k)] = –Im[Xp(N – k)],
Xp(k) = Xp(N – k) , (5.38)
arg Xp(k) = –arg Xp(N – k).
Аналогичные равенства справедливы и для конечной действительной последовательности x(n), имеющей N-точечное ДПФ Xp(k). Если ввести дополнительное условие симметрии последовательности xp(n), т. е. считать, что xp(n) = xp(N – n), то окажется, что Xp(k) может быть только действительной.
Чаще всего приходится иметь дело с действительными последовательностями, поэтому, вычислив одно ДПФ, можно получить ДПФ двух последовательностей, используя свойства симметрии (5.38). Рассмотрим действительные периодические последовательности xp(n) и yp(n) с периодами в N отсчетов каждая и N-точечными ДПФ Xp(k) и Yp(k) соответственно. ДПФ комплексной последовательности
zp(n) = xp(n) + j yp(n)
равно
Zp(k) = [xp(n) + j yp(n)] ;
Zp(k) = Xp(n) + j Yp(n). (19)
Выделяя действительную и мнимую части равенства (19), получим
Re[Zp(k)] = Re[Xp(k)] – Im[Yp(k)];
Im[Zp(k)] = Im[Xp(k)] + Re[Yp(k)].
Действительные части Xp(n) и Yp(n) симметричны, а мнимые — антисимметричны, поэтому их легко разделить, используя операции сложения и вычитания:
Re[Xp(k)] = {Re[Zp(k)] + Re[Zp(k)]}/2;
Im[Yp(k)] = {Re[Zp(k)] – Re[Zp(k)]}/2;
Re[Yp(k)] = {Im[Zp(k)] + Im [Zp(k)]}/2;
Im[Xp(k)] = {Im[Zp(k)] – Im [Zp(k)]}/2.
Итак, вычисляя одно N-точечное ДПФ, удается преобразовать сразу две действительные последовательности длиной по N отсчетов. Если эти последовательности еще и симметричные, то число операций, необходимых для получения их ДПФ, можно сократить еще больше.
3. Спектральный анализ в точках z-плоскости
Рис. 9. Фильтр для скользящего спектрального анализа:
блоки с обозначением z–1 – элементы задержки;
величины, равные степеням z1, – коэффициенты умножителей
Спектральный анализ можно рассматривать как задачу вычисления z-преобразования модифицированного сигнала в некоторой области на z-плоскости. Спектральные составляющие сигнала х(п) можно измерять в любой точке z1 на z-плоскости
Sn(z1) = , (20)
где N — число отсчетов, по которым находится оценка спектра.
Во многих приложениях, например, когда спектр сигнала меняется, приходится измерять Sn(z1) для последовательных значений п, т. е. значения S0(z1), S1(z1), S2(z1) и т. д. Такой способ измерений называют скользящим спектральным измерением; оно обеспечивается за счет смещения на один отсчет вперед временнόго окна (содержащего N отсчетов) и повторения измерения. Из формулы (5.40) видно, что скользящее спектральное измерение в одной точке z = z1 эквивалентно преобразованию фильтром с импульсной характеристикой вида
h(n) = , 0 n N – 1, (21)
По формуле (5.40) составляется схема вычисления прямой свертки, обеспечивающая спектральные измерения – рис. 5.9.
По выражениям для двух последовательных спектральных измерений Sn – 1(z1) и Sn(z1), можно получить рекуррентную формулу:
Sn (z1) = Sn – 1(z1) + x(n) – x(n – N) (22)
– рис.10.
Рис. 10. Рекурентный метод скользящего спектрального анализа
Лекция 6
Do'stlaringiz bilan baham: |