Главная  Цепи и сигналы 

[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] [32] [33] [34] [35] [36] [37] [38] [39] [40] [41] [42] [43] [44] [45] [46] [47] [48] [49] [50] [51] [52] [53] [54] [55] [56] [57] [58] [59] [60] [61] [62] [63] [64] [65] [66] [67] [68] [69] [70] [71] [72] [73] [74] [75] [76] [77] [78] [79] [80] [81] [82] [83] [84] [85] [86] [87] [88] [89] [90] [91] [92] [93] [94] [95] [96] [97] [98] [99] [100] [101] [102] [103] [104] [105] [106] [107] [108] [109] [110] [111] [112] [113] [114] [115] [116] [117] [118] [119] [120] [121] [122] [123] [124] [125] [126] [127] [128] [129] [130] [131] [132] [133] [134] [135] [136] [137] [138] [139] [140] [141] [142] [143] [144] [145] [146] [147] [148] [149] [150] [151] [152] [153] [154] [155] [156] [157] [158] [159] [160] [161] [162] [ 163 ] [164] [165] [166] [167] [168] [169]

в пределах одного интервала Лш набег фазы значительно меньше л. На частотах со = пАсо, где скачо 9, (со) превьш1ает 4-л, для восстановления истинного аргумента требуется добавить - 2л; соответственно при скачке

- л требуется добавить + 2л. На рис. 16.17, в изображена корректируюш,ая последовательность, добавляемая к последовательности Qsrji (")•

В заключение рассмотрим важный для практики случай входной последовательности {s {т.)}, m О, когда z-преобразование (одностороннее) определяется выражением

S(z)= V ${т)г-", (16.52)

причем полюсы и нули функции S (z) расположены внутри единичной окружности, т. е. радиус сходимости ряда (16.52) < 1.

Подобные последовательности называются минимально-фазо-

в ы м и, по аналогии с системами, передаточная функция которых К (z) имеет особые точки (полюсы и нули) внутри круга z = 1 (см. § 15.9).

Таким образом, модуль и аргумент z-преобразования минимально-фазовой последовательности однозначно связаны, а комплексный логарифм

In S(e) = In 1 S(e) \ + iG (со)

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

CJm)==y{m) = 2x~ Г in 1 S~(e") cos (mwT) d (соГ),

которая отличается от (16.20) только степенью S (eOI в аргументе логарифма.

Из последнего выражения видно, что комплексный кепстр минимально-фазового сигнала s (т.) равен кепстру мощности (т) того же сигнала, и для его получения можно воспользоваться схемой, представленной на рис. 16.11.

Простейшим примером минимально-фазового сигнала является s (/) =

= е-«, fO, с г-преобразованием S {е) - ~ име-

ющим нуль Zo = О и полюс Zp = е~« <: 1 [см. (12.22)].

Сигнал S (/) = Ate.-, рассмотренный в § 16.7, является другим примером минимально-фазового сигнала.

16.11. НЕКОТОРЫЕ СВОЙСТВА КЕПСТРАЛЬНЫХ ПРЕОБРАЗОВАНИЙ

По аналогии с теоремами о спектрах, изложенными в § 2.8, рассмотрим связь между некоторыми преобразованиями исходного сигнала s (t) и преобразованиями кепстра.

Установление этих связей представляет интерес в основном применительно к комплексным кепстрам.

1. СДВИГ СИГНАЛА ВО ВРЕМЕННОЙ ОБЛАСТИ

Пусть исходному сигналу (/) со спектральной плотностью Si (со) соответствует кепстр Csi ("?) При задержке сигнала на tg получим функцию времени (t) = Sj (/ -

- to) со спектральной плотностью (ы) = e~"°Si (w), логарифм которой

In Sj (ы) = - ibito + In Si (to).



Второму слагаемому соответствует кепстр Cj (q) исходного сигнала Si (/), а первое слагаемое в соответствии с (16.39) приводит к следующему кепстру:

с»

Cto{q) = -to- J iae<ld(o. (16.53)

- оо

Учитывая, что функции б (q) соответствует спектральная плотность, равная единице, множитель 1(0 можно рассматривать как спектральную плотность функции

-6(9). Тогда (16.53) определяет функцию б (о), и, следовательно, dq dq

Cjo (?) = -о -J- б (9). (16.54)

Таким образом, кепстр сигнала (/) = (t - tg)

Cs2 (<?)= -0 -J- б {q)-{-C,, (q). (16.55)

Из сопоставления (16.55) и (16.54) вытекает, что кепстр смещенной дельта-функции

{I - tn) равен - б (q). Оперирование производной дельта-функции затрудни-

тельно, однако при обработке сигнала можно исключить участок кепстрального времени в окрестности точки q = 0.

Рассмотрим соотношение между cs2 (п) и С (т) для цифрового сигнала (т) = = Si (т - /По).

Основываясь на методе г-преобразования, получаем

%{е)= е-"»" S;(e"), Insje) =-ш„ ыТ + \п% {е<).

Применяя к этому выражению обратное г-преобразование по формуле (12.28), получаем

Cs,(n)=-imo J ыТе"" d (ыТ) + Csi (т),

где (m) - кепстр сигнала (т), а

Сго {т)=-Шо J хе"" dx.

При т-0 интеграл обращается в нуль, так что С,„ (0)=0. При тф О

я я я

I хе"" djc= I X cos (mx)dx + i j" [x sin (mx) dx = -я -я -я

sin (mn) -тя cos (тя) „ я

= 2i ----;;--- = - 2« - cos (тя

m" m

-2i - cos (тя) m

c,„(m)= 2„

Таким образом, при m ф О

Cs2 (m)=- (-O--l-Cs, (m). m

Как видим, в случае цифрового сигнала кепстр cs2 (т) задержанного сигнала отличается от (т) лишь знакопеременным сигналом.. (-1)""", убывающим с воз-

растанием т; в точке m = О Cjo (0)= О и дельта-функция не возникает.



2. ИЗМЕНЕНИЕ МАСШТАБА ВРЕМЕНИ

Пусть Sa (t) = Si (nt), re > 0. В соответствии с § 2.8, п. 2

Sa (ы) = - Si (w/re) и In (ы) - -\п {п)+\п Sj (ы/ге). re

Кепстр сигнала (t)

Cs2 (?)=-ln (re) -L j e"dm +

- OO

In Si (ш/ге) e"" d(o= -In (re) б {q) + nCsi (nq).

Изменение масштаба времени t приводит к такому же изменению масштаба кепстрального времени q; кроме того, возникает функция б (q).

При дискретизованном сигнале изменение масштаба времени означает изменение шага Т при неизменном числе отсчетов N (что необходимо для сохранения формы сигнала).

Положим Га = rei, re > О, и запишем выражение (12.20) для г-преобразования сигнала s, (t)

N-l N-l

S;(e"0= 2 «2 (m) e-"»"i= 2 % (n) e-»". in = 0 m=0

При сжатии или растяжении исходного сигнала отсчеты функции {i) сохраняют свое значение (при = const). Таким образом,

S;(e-"-.) Jy si (m) е- =Si (e") m = 0

И кепстр

(m) =

In S(e"*0 cos (отгеыГ,) d (rewTi).

Переход от шага дискретизации к = reTi не изменяет структуры кепстра. Изменяется лишь диапазон частот ш, соответствующий одному обходу окружности единичного радиуса на г-плоскости (от -л/Т, л/Т до - л/пТ,, л/пТ,). Соответственно изменяется и масштаб кепстрального времени; интервалы между отсчетами кепстра на оси q будут = пТ.

3. СМЕЩЕНИЕ СПЕКТРА СИГНАЛА Применим преобразование Фурье к произведению "а (0=40 е»,

где S (О - «медленная» (модулирующая) функция со спектральной плотностью S (ы); е"*° - несущее колебание.

Повторяя рассуждения, приведенные в §2.7, п. 3, придем к спектру

с»

Sa (ш)= J S (О е"»е-"d/=S (й)-(0о).

Тогдг

In Sa (o)) = ln S {(О-СОо) и кепстр сигнала а (i)

Ca{q)=~ J In S (ы-(Оо) е""d(0.



[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] [32] [33] [34] [35] [36] [37] [38] [39] [40] [41] [42] [43] [44] [45] [46] [47] [48] [49] [50] [51] [52] [53] [54] [55] [56] [57] [58] [59] [60] [61] [62] [63] [64] [65] [66] [67] [68] [69] [70] [71] [72] [73] [74] [75] [76] [77] [78] [79] [80] [81] [82] [83] [84] [85] [86] [87] [88] [89] [90] [91] [92] [93] [94] [95] [96] [97] [98] [99] [100] [101] [102] [103] [104] [105] [106] [107] [108] [109] [110] [111] [112] [113] [114] [115] [116] [117] [118] [119] [120] [121] [122] [123] [124] [125] [126] [127] [128] [129] [130] [131] [132] [133] [134] [135] [136] [137] [138] [139] [140] [141] [142] [143] [144] [145] [146] [147] [148] [149] [150] [151] [152] [153] [154] [155] [156] [157] [158] [159] [160] [161] [162] [ 163 ] [164] [165] [166] [167] [168] [169]

0.001