Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

МАРТ_Часть1_Главы1_5

.pdf
Скачиваний:
0
Добавлен:
21.03.2026
Размер:
1.76 Mб
Скачать

В частотной области это уравнение примет вид (убедитесь в этом само-

~

 

~

~

 

 

 

 

~

~

0 .

стоятельно) s0 1

 

p

s0 d 2s0

g

 

 

 

 

 

 

 

 

 

 

 

 

 

Используем полученный результат для определения сигнала, у которого

произведение временной протяженности на ширину спектра будет мини-

мально. При этом временную протяженность Тэфф определим как

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Tэфф2

t t0 2 s2 t dt

s2 t dt ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а ширину спектра эфф как

 

 

 

 

 

 

 

 

 

 

 

2

 

 

2

 

~

 

2

 

 

~

 

2

 

 

 

 

 

d

 

 

d ,

эфф

0

 

s

 

 

s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где t0 и 0 – числа, характеризующие положение сигнала на оси t и спектра на оси . Без потери общности они могут быть положены равными нулю.

 

 

 

 

 

 

 

 

1

 

 

Кроме того, с учетом введенной выше нормировки s2 t dt

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

~

 

 

2

 

 

 

 

 

 

 

d 1. Ясно, что максимизация

Tэфф эфф или, что то же самое,

s

 

 

 

 

 

 

 

 

 

 

 

T 2

2

 

 

 

эквивалентна минимизации одного из сомножителей при фикси-

эфф

эфф

 

 

 

 

рованном значении другого. Таким образом, мы имеем предыдущую задачу с

р(t) = t

2

и

~

2

. Так как функции

~

2

 

g

 

g

во временной области соот-

ветствует g(t) = – ''(t), где (t) – дельта-функция (проверьте это), то после подстановки g(t– ) в уравнение (5.9) и использования интегральных свойств

 

t 1

1t

2

s0

t 0. Это ли-

производной дельта-функции, получим 2s0

 

нейное дифференциальное уравнение с переменными коэффициентами называется уравнением Эрмита (в квантовой механике его называют уравнением Шредингера).

Его решениями, собственными функциями дифференциального опера-

 

 

 

 

 

 

 

 

2

 

 

 

тора Эрмита, являются эрмитовы сигналы

1

 

 

e t

2H

 

t . Глобальный

 

 

 

 

k

2k k!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

минимум достигается при k = 0 и соответствует гауссовскому импульсу

e

t 2

2

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

81

1,

 

t

 

 

T 2 ,

 

1,

 

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

 

 

 

 

 

 

 

 

мы получим уже решенную задачу об

 

 

 

 

 

 

 

 

 

 

p t

T 2

и q

 

 

 

 

 

 

0,

t

 

0,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

импульсе, обеспечивающем максимальную концентрацию энергии в полосе частот [– , ].

Оператор Гильберта. Одним из основных понятий теории сигналов является комплексная функция s t , образованная из исходного вещественного сигнала s(t) добавлением мнимой составляющей s (t), которая определя-

1 s

ется преобразованием Гильберта исходного сигнала s t t d . Эта

комплексная функция называется аналитическим сигналом, т. е. s t s t js t .

Вещественная и мнимая части аналитического сигнала позволяют определить огибающую (модуль аналитического сигнала), характеризующую закон амплитудной модуляции и фазу (аргумент аналитического сигнала), определяющие закон угловой модуляции следующим образом:

– огибающая S t s2 t s2 t ,

– фаза t arctg s t .

s t

При этом аналитический сигнал s t может быть записан в форме s t S t e j t . Естественно, если мнимую часть аналитического сигнала

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

изводная от фазы, т. е. t d t . dt

Таким образом, чтобы однозначно определить огибающую, фазу и частоту, необходимо и достаточно указать правило построения мнимой части v(t) по исходному сигналу s(t), т. е. указать оператор А, осуществляющий преобразование v(t) = As(t).

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

1. Потребуем, чтобы малым изменениям исходного сигнала s(t) соответствовали бы малые изменения огибающей S(t) и фазы Ф(t). Так как функ-

82

ции, определяющие S(t) и Ф(t), непрерывны, то сформулированное требование означает непрерывность оператора А.

2. Фаза, а следовательно, и мгновенная частота не должны зависеть от

энергии (нормы) сигнала. Это означает, что оператор А должен быть однородным, т. е. A ks t kA s t . В этом случае

t arctg

A ks t

arctg

A s t

.

ks t

 

 

 

 

s t

3. Для гармонического сигнала s(t) = Um cos ( 0t + ) мы должны получить S(t) = Um и Ф(t) = 0t + . Это означает, что А[Umcos( 0t + )] =

=Umsin ( 0t + ).

В[14] доказано, что единственным линейным оператором, удовлетво-

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

Остановимся на свойствах оператора Гильберта, обозначаемого обычно буквой Н. Это – линейный оператор, являющийся сверткой исходного сигнала и функции 1 t . Обратный оператор Н–1 может быть найден путем решения уравнения s (t) H s(t) относительно s(t) переходом в частотную область с помощью преобразования Фурье. Учитывая правило преобразования Фурье для свертки, получим:

~

~

 

1

s ( ) s ( )F

 

,

 

 

 

 

t

 

1

 

 

e j

 

j

 

sin t

 

sin t

 

где F

 

 

 

 

 

dt

 

 

 

dt . Интеграл

 

dt sign , где

 

t

 

 

 

 

t

 

 

 

 

t

 

t

 

 

 

 

 

 

 

 

 

 

 

1, x 0;

sign x 0, x 0; – знаковая функция.

1, x 0

 

1

Таким образом, преобразование Фурье F

 

, которое можно рас-

 

 

t

сматривать как комплексный коэффициент передачи линейного фильтра,

 

 

1

осуществляющего преобразование Гильберта, равно F

 

= j sign . Сле-

 

 

 

t

~

~

 

 

довательно, s ( ) s ( ) j sign .

 

 

83

Амплитудно-частотная характеристика такого фильтра F 1 t 1, а

 

 

 

1

 

2 , 0;

 

 

фазо-частотная представляет собой arg F

 

 

 

 

 

 

 

. Используя

 

 

 

 

 

 

 

t

 

 

2 ,

0

 

~

~

 

 

1

1

~

 

 

 

 

 

данный результат, получим: s ( ) s ( ) F

 

 

 

 

s ( ) j sign , т. е., в

 

 

 

 

 

 

t

 

 

 

 

отличие от прямого преобразования, изменился лишь знак. Следовательно, обратное преобразование Гильберта имеет вид:

s(t) 1 s ( ) d .

t

С учетом частотного представления оператора Гильберта, H 2 H H E , так как sign 2 1.

1

Оператор Гильберта унитарен, так как для его ядра (t ) выполняется условие (5.4), что доказывается с помощью обобщенного равенства Парсеваля:

 

1

1

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

(t 1) (t 2 )

 

 

=

1

 

 

 

 

 

2

 

1

 

j sign e

j 1

j sign * e j 2 d

 

 

2

 

 

 

 

e j(1 2 )d 1 2 .

Исходная функция s(t) и функция, преобразованная по Гильберту, ортогональны, что также доказывается с помощью обобщенного равенства

Парсеваля:

 

1

~

s(t)s (t)dt

 

s ( )

2

 

 

~

*

 

j

 

~

 

2

 

 

 

 

 

 

 

j sign s ( )

d

 

 

s ( )

 

 

sign d 0

,

2

 

 

 

 

 

 

 

 

 

 

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

Норма оператора Гильберта равна 1. Он не меняет энергии преобразуемых сигналов, т. е. Hs(t) 2 s(t) 2 .

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

= s(t) + js (t) = S(t)e j(t) = S(t)e j 0 (t) .

84

cos 0t , счи-

Комплексная функция S(t) S(t)e

~ s ( )

0

0

 

 

 

0

 

 

~

 

 

 

s ( )

 

 

0

0

 

~ S ( )

 

 

 

 

0

 

 

Рис. 5.5

 

 

j(t) , определяющая законы амплитудной S(t) и угловой (t) модуляции, называется комплексной огибающей. Зная комплексную огибающую, можно определить исходный сигнал s(t) =

= Re S(t)e j 0 (t) . Отметим, что комплексная огибающая является видеофункцией (центральная частота спектра равна нулю). Скорость ее изменения определяется шириной спектра исходного сигнала. Для иллюстрации вычислим спектр аналитического сигнала и комплексной огибающей для колокольного (гауссовского) импульса s(t) = Ume t 2

тая его узкополосным, т. е.

1

 

T

, T

2

– период высоко-

 

 

 

 

 

 

0

0

0

 

 

 

 

 

 

частотного заполнения нашего сигна-

ла.

В соответствии с теоремами о спектрах получим:

 

 

 

 

 

 

 

 

 

2

 

 

 

2

 

~

Um

 

 

 

 

 

 

 

 

 

exp

 

0

 

exp

 

0

.

s ( )

 

 

 

 

 

 

4

 

4

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Спектр аналитического сигнала равен F s(t)

 

= Fs(t) + jFs (t). Как было пока-

~

~

 

 

, поэтому

 

 

зано выше, Fs (t) = s ( ) = s ( ) j sign

 

 

 

~

~

 

 

 

 

 

 

 

2s ( ), 0;

 

 

 

 

Fs(t) s

( )

0.

 

 

 

 

 

 

0,

 

 

 

 

 

 

 

 

 

 

 

 

 

j 0 (t)

,

Для спектра комплексной огибающей, с учетом того, что S(t) s(t)e

 

и теоремы о смещении спектра, выполняется соотношение:

 

 

~

 

 

 

 

 

 

 

 

 

 

FS (t) s ( ) . Таким образом, для рассматриваемого сигнала оконча-

 

0

 

 

 

 

 

 

 

 

 

 

тельно получим:

 

 

 

 

 

 

 

 

 

 

 

~

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

FS(t) S ( ) Um

 

exp

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

85

Рисунок 5.5 иллюстрирует сказанное.

Комплексная огибающая упрощает решение задачи о преобразовании узкополосного сигнала s(t) S(t) cos 0t (t) линейной системой с импульсной характеристикой h(t) H (t) cos 0t (t) , где S(t) и H(t) – огибающие сигна-

ла и импульсной характеристики, а (t) и (t) – отклонение фазы от линейного закона.

С помощью интеграла Дюамеля комплексная огибающая выходного

 

 

1

 

 

 

 

 

сигнала Sвых (t) находится как Sвых (t) =

2

 

S ( )H (t ) d , где S (t) и H (t)

 

 

 

 

 

 

комплексные огибающие входного сигнала и импульсной характеристики фильтра. Выходной сигнал находится как sвых(t) Re Sвых(t)e j 0t .

Контрольные вопросы

1.Дайте определение оператора, функционала, функции.

2.Как определяется понятие обратного оператора?

3.Дайте определение ограниченного оператора. Что называется нормой оператора?

4.Какой оператор называется линейным?

5.Как определяется произвольный линейный оператор в Rn?

6.Какой оператор называется тождественным? Приведите примеры тождественных операторов для Rn и L2.

7.Приведите примеры интегральных операторов Фредгольма и Вольтерра. Сформулируйте требования, предъявляемые к ядрам этих операторов.

8.Дайте определение суммы операторов, произведения и функции от оператора.

9.Какой оператор называется разрешающим, или резольвентой, для оператора А?

10.Что такое спектр оператора, дискретный спектр, собственные зна-

чения?

11.Что такое собственные вектора оператора А?

12.Приведите примеры линейных дифференциальных операторов.

13.Дайте определение линейного функционала. Приведите примеры

для Rn и C[a, b].

14.Поясните геометрический смысл линейного функционала.

15.Дайте определение выпуклого множества.

86

16.Что называется выпуклой оболочкой множества А?

17.Дайте определение симплекса.

18.Дайте определение выпуклого функционала.

19.Как можно представить любой линейный функционал в гильбертовом пространстве Н?

20.Как определяется норма линейного функционала?

21.Как определяется оператор А*, сопряженный оператору А?

22.Приведите примеры линейных сопряженных операторов для Rn и L2.

23.Что такое квадратичная форма? Приведите примеры для Rn и L2.

24.Дайте определение самосопряженного оператора. Приведите при-

меры для Rn и L2.

25.Что такое положительно (неотрицательно) определенная квадратичная форма?

26.Как формируются условия для положительно определенной квадратичной формы в Rn?

27.Сформулируйте свойства собственных значений и собственных векторов самосопряженных операторов.

28.Сформулируйте теорему Гильберта–Шмидта.

29.Дайте определение самосопряженного дифференциального оператора второго порядка.

30.Как записывается оператор Штурма–Лиувилля? Приведите пример его использования в задачах радиотехники.

31.Что такое нормальный оператор, унитарный оператор?

32.Как определяется понятие дельта-функции и ее производных?

33.В чем состоит «фильтрующее» свойство дельта-функции и ее производных?

34.Запишите условие, при котором оператор Фредгольма будет уни-

тарным.

35.Как определяется оператор проектирования? Приведите примеры.

36.Дайте определение оператора Фурье. Докажите, что он является унитарным оператором.

37.Сформулируйте основные свойства оператора Фурье (теоремы о спектрах).

38.Какие функции являются собственными для оператора Фурье?

39.Дайте определение текущего и мгновенного спектра.

40.Как определяется понятие кепстр сигнала? В каких задачах целесообразно использовать кепстральное представление сигнала?

87

41.Как определяется оператор Фурье на конечном интервале?

42.Что такое функции с двойной ортогональностью? Каковы их свойства?

43.Сформулируйте результат решения задачи об отыскании формы конечного сигнала, отличного от нуля лишь в промежутке [–T/2, T/2] и имею-

щего максимум энергии в полосе частот [– , ].

44.Сформулируйте принцип решения задачи об отыскании сигнала, обладающей требуемой концентрацией энергии на плоскости время-частота.

45.Дайте определение оператора Гильберта.

46.С помощью оператора Гильберта определите огибающую сигнала, его фазу и мгновенную частоту.

47.Докажите, что оператор Гильберта унитарен.

48.Определите коэффициент передачи фильтра, реализующего преобразование Гильберта.

49.Докажите, что исходный сигнал s(t) и преобразованный по Гильбер-

ту s (t) ортогональны.

50.Дайте определение узкополосного сигнала.

51.Для узкополосного сигнала определите понятие комплексной огиба-

ющей.

88

Соседние файлы в предмете Математический аппарат радиотехники