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

1490-spectry_i_analiz

.pdf
Скачиваний:
119
Добавлен:
11.05.2015
Размер:
3.68 Mб
Скачать

 

x0

ω

 

 

x0

 

 

 

 

x0ωВ

sinωВt

 

x(t )=

В

exp[jωt]dω =

exp[jωt]

 

ωωВ В

=

(5.3)

 

 

 

 

 

2π

 

2 jπt

 

 

π ωВt

 

 

ω

В

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и называется функцией отсчетов [25]. Происхождение этого термина и смысл функции отсчетов будет рассмотрен ниже. Иногда функцию (5.3) называют идеальным низкочастотным сигналом, поскольку эту функцию можно рассматривать как выходную реакцию фильтра нижних частот, обладающего прямоугольной амплитудно частотной характеристикой, на входной сигнал в виде δ функции, имеющей равномерную, бесконечно протяженную спектральную плотность.

График функции отсчетов приведен на рис. 5.3.

Рис 5.3. Идеальный низкочастотный сигнал (функция отсчетов)

Функция отсчетов принимает нулевые значения при пересечении оси времени в точках t =n / 2 fВ (n=1,2…N), где fВ - верхняя частота спектра

X (ω)

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

 

0,

ω < −ωB ,

 

X1

 

 

(5.4)

(ω)= x0 exp{jωt0}, ωB ω ωB ,

 

 

ω >ωB

 

 

0,

 

202

Функция отсчетов, отвечающая спектральной плотности (5.4) определяется как:

 

 

)

 

x0

ωB

 

 

 

]

 

 

 

 

 

 

 

 

x0

ωB

 

 

 

 

x

t

=

exp

jωt

exp

jωt

dω=

 

exp j

t t

ω dω=

 

 

 

 

 

1

(

 

2π

 

[

0

[

]

 

 

 

2π

(

0 )

 

 

 

 

 

 

ω

 

 

 

 

 

 

 

 

 

 

 

ω

 

 

 

 

 

 

 

 

 

B

 

 

 

 

 

 

 

 

 

 

 

 

 

B

 

 

 

 

 

=

 

 

x0

 

exp[j(t t0 )ω]

 

ωωB

 

=

x0ωB

sinωB (t t0 ) ,

(5.5)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 jπ(t t0 )

 

B

 

 

π ωB (t t0 )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

т.е. представляет собой функцию x(t), смещенную по оси времени на величину t0 , что соответствует теореме задержки для преобразования Фурье.

Пример 2. Высокочастотная функция отсчетов (идеальный полосовой сигнал).

Рассмотрим сигнал, спектральная плотность, которого ограничена двумя частотными интервалами шириной 2ωB , каждая с центрами на частотах

±ω0 (см. рис. 5.4а)

 

 

 

 

 

 

x , ω ω

B

ω ≤ −ω +ω

,

 

 

0

0

0 B

 

 

X 2

(ω)= x0 , ω0 ωB ω ω0 +ωB ,

 

(5.6)

 

 

для всех других значений

ω.

 

0

X (ω)

x0

 

0

ω0 (ω0 +ωBω)

(ω0 +ωB ) ω0 (ω0 ωB ) (ω0 ωB )

Рис. 5.4а. Спектральная плотность идеального полосового сигнала

Функцию времени x2 (t ), отвечающую спектральной плотности (5.6),

найдем, используя обратное преобразование Фурье и учитывая при этом четность функции X (ω):

203

 

 

 

 

(t )=

x0

(ω0 +ωB )

 

 

x0

 

 

 

ω0 +ωB

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

ω

cosωt dω =

sinωt

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

π (ω

 

 

 

πt

 

 

 

 

 

 

 

 

 

 

 

)

 

 

 

ω ω

B

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

B

 

 

 

 

 

 

 

0

 

 

=

x0

sin (ω

 

 

+ω )t sin

(ω ω

 

)t

=

2x0ωB

sinωBt

cosω t.

(5.7)

 

 

 

 

 

 

πt

0

B

 

 

 

0

B

 

 

 

π

 

ωBt

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Выражение (5.7) представляет собой гармоническое колебание cosω0t ,

умноженное на функцию отсчетов sinωBt /ωBt . С точки зрения теории сиг-

налов функция x2 (t ) представляет собой балансно модулированный сиг-

нал, где несущим является сигнал cosω0t , а модулирующим – функция от-

счетов, которая осциллирует значительно медленнее модулирующего сигнала. Особенности балансной модуляции будут подробно рассмотрены в главе 6, посвященной анализу спектров модулированных сигналов. Функ-

цию x2 (t ) будем называть высокочастотной функцией отсчетов. Возмож-

ный способ получения такого сигнала заключается в подаче δ функции, обладающей равномерным, бесконечно протяженным спектром на вход идеального полосового фильтра с прямоугольной амплитудно частотной характеристикой в интервале (ω0 ωB , ω0 +ωB ). Высокочастотная функ-

ция отсчетов изображена на рис. 5.4б.

1

0,5

0

-0,5

-1

Рис. 5.4б. Высокочастотная функция отсчетов

На рисунке ясно видны характерные скачки фазы на 180° в точках, где огибающая sinωBt /ωBt принимает нулевое значение.

204

5.3 Ортогональные сигналы с ограниченным спектром

Рассмотренные выше способы разложения периодических и непериодических сигналов в ряд и интеграл Фурье по ортогональным тригонометрическим функциям не являются единственными. Существуют также другие классы ортогональных функций и полиномов, которые могут быть использованы для разложения функций времени. Один из таких возможных классов ортогональных функций может быть построен на основе использования функций с ограниченным спектром [1, 7]. Для этого необходимо рассмотреть метод ортогонализации этих функций. Для решения задачи зададим две функции:

x(t )=

x ω

sinω

B , x1 (t )= x(t t0 )=

x ω

 

sinωB (t t0 )

0 B

 

0 B

 

 

π ωBt

π ωB (t t0 )

 

 

и, следуя изложенному в главе 1, запишем условие ортогональности этих функций в виде равенства нулю их скалярного произведения

 

 

 

 

2 2

 

sinωB (t t0 )

 

 

x(t ) x(t t

 

)

=

x0ωB

−∞

sinωBt

dt = 0.

(5.8)

 

 

 

 

0

 

 

π2

ωBt ωB (t t0 )

 

Для вычисления интеграла (5.8), который представляет собой свертку двух функций времени, произведем преобразование Фурье над правой и левой частями выражения (5.8):

 

(

 

)

 

(

 

0 )

[

 

]

0

[

 

0

]

 

 

t

x

t t

jωt

jωt

,

(5.9)

 

x

 

 

 

exp

 

 

dt = x2 exp

 

 

 

−∞

поскольку преобразование Фурье свертки есть произведение преобразований Фурье свертываемых функций, а спектральные плотности свертываемых функций определяются соотношениями (5.2) и (5.4).

Применяя теперь к выражению (5.9) обратное преобразование Фурье,

найдем скалярное произведение функций x(t ) и x1 (t ):

205

 

 

 

x2 ωB

exp{jωt0

}exp{jωt}dω =

 

 

 

 

 

0

ω

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2π

 

 

 

 

 

 

 

 

 

 

B

 

 

 

 

 

 

=

x02

ωb

exp{jω(t t0 )}dω =

x02ωB

sinωBt

,

(5.10)

 

 

 

2π ω

 

 

 

 

π ωBt

 

 

 

 

b

 

 

 

 

 

 

где t =t t0 означает временной сдвиг между функциями отсчета x1 (t ) и x(t ). Из выражения (5.10) следует, что скалярное произведение обращает-

ся в нуль и, соответственно, две функции отсчета будут ортогональны, если взаимный сдвиг этих функций во времени удовлетворяет условию

ωBt = nπ, (n = ±1,±2,...),

откуда

 

 

 

 

 

 

t = nπ =

 

n

 

(n = ±1,±2,...)

(5.11)

 

2 fB

ωB

 

 

 

 

 

Минимально возможный сдвиг t , при котором функции отсчета орто-

гональны, имеет место для n = ±1

 

 

 

 

 

 

 

t =

 

1

.

 

 

 

 

 

 

 

 

2 fB

 

Таким образом, используя бесконечную последовательность функций отсчета, в которой каждая последующая функция сдвинута относительно предыдущей на величину t =1/ 2 fB , можно построить ортогональную

систему функций, которая будет служить в качестве базиса для разложения

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

та которого есть ωB . Взаимное расположение функций отсчета, образую-

щих ортогональную систему, показано на рис. 5.5.

5.4Теорема отсчетов во временном представлении

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

206

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

относящиеся к теории интерполяции функций с ограниченным спектром,

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

лее строгой теорией, а также будет дана краткая справка об истории во-

проса.

Рис. 5.5.

Как в задачах математического анализа и его приложений, так и в зада-

чах радиоэлектроники и радиофизики возникает проблема интерполирова-

ния функции, т.е. задача определения искомой функции x(t ) для всех зна-

чений t , заключенных в некотором интервале a t b , если известны зна-

чения x(t1 ), x(t2 ),... x(tn ) в ряде точек t1, t2 ,... tn [19]. Задача, поставленная

таким образом, может иметь большое (вплоть до бесконечности) количест-

во решений. Однако ограничение класса функций x(t ) может привести к

существованию только единственного решения. Кроме того, при ограни-

чении свойств класса функций, из которых выбирается интерполирующая функция, задача становится более определенной.

207

Рассмотрим непрерывный сигнал x(t ), стремящийся к нулю при t →∞

(что обеспечивает его интегрируемость) и будем полагать, что этот сигнал

имеет ограниченный спектр. Для сигналов этого типа докажем следующую

теорему [7, 25].

Теорема. Если функция x(t ) не содержит составляющих с частотой

выше fB герц, то она полностью определяется последовательностью её значений в точках, отстоящих на расстояние 1/ 2 fB секунд друг от друга.

Данная формулировка представляет собой теорему отсчетов во временном представлении.

Приступая к доказательству теоремы, запишем прежде представление

функции x(t ) через её

спектральную плотность, учитывая,

что спектр

функции заключен в интервале (ωB , ωB )

 

x(t )=

1

2π f

B

X (ω)exp[jωt]dω.

(5.11)

 

 

 

 

 

 

2π 2π fB

 

 

Будем рассматривать значения t , удовлетворяющие условию:

 

 

t =

n

 

 

, (n = ±1, ± 2,...) ,

(5.12)

 

2 fB

 

 

 

 

 

где n – целые числа. Для этих моментов времени выражение (5.11) принимает вид:

 

n

 

 

1

2π fB

 

 

X (ω)exp j (n / 2 fB )ω dω.

 

x

 

=

(5.13)

 

2π

 

2 fB

 

2πfB

 

Из выражения (5.13) следует, что совокупность значений x(n2 fB ) пол-

ностью определяет спектральную плотность X (ω), которая, в свою оче-

редь, однозначно определяет функцию

x(t )=

1

2π fB

X (ω)exp j (n / 2 fB )ω dω.

2π

 

 

 

 

 

2π fB

 

 

208

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

дающая ограниченным спектром X (ω) с верхней частотой fB , которая в точках отсчета t = n / 2 fB принимает заданные значения x(n / 2 fB ).

Теперь необходимо найти аналитическое выражение для функции, заданной ее значениями в точках отсчета.

Используя периодическое продолжение спектральной плотности X (ω),

заданной на частотном интервале 2π fB ω 2π fB (рис.5.6),

Рис. 5.6. Периодическое продолжение ограниченной спектральной плотности

разложим эту функцию в ряд Фурье:

 

 

2πn

 

 

 

 

 

 

 

 

n

 

 

 

 

X (ω)= ∑

 

 

 

 

 

 

 

 

 

 

 

С&n exp j

 

 

 

 

 

ω =

C&

exp j

 

 

 

ω

, (5.14)

4π f

 

 

2π f

 

n=−∞

 

 

B

 

n=−∞

n

 

 

 

 

B

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где комплексные коэффициенты Фурье определяется выражением

 

 

 

1

 

2π fB

 

 

 

 

n

 

 

 

 

 

 

 

 

 

C&n =

 

 

 

 

x(ω)exp

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j

 

 

 

 

ω dω,

 

 

 

(5.15)

 

4π f

 

 

 

 

2 f

 

 

 

 

 

 

 

 

b 2π f

B

 

 

 

 

 

B

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а величина 4π fB

есть период периодически продолженной спектральной

плотности. Сравнивая выражения (5.13) и (5.15), можно записать

 

209

C&n =

1

x(n / 2 fB ),

(5.16)

 

 

2 fB

 

 

откуда следует, что коэффициенты C&n

разложения периодически продол-

женной спектральной плотности x(ω)

в ряд Фурье есть ни что иное, как

отсчеты функции x(t ), взятые в моменты времени t = n / 2 fB и взвешенные с весом 1/ 2 fB . Тогда ряд Фурье (5.14) можно записать в виде

 

1

 

 

π

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

X (ω)=

 

 

x

 

 

exp j

 

 

 

ω .

2 f

 

2 f

 

2 f

 

 

B

n=−∞

 

B

 

 

B

 

 

 

 

 

 

 

 

Из выражений (5.11) и (5.17) следует, что

 

 

 

 

 

 

 

 

 

1

 

 

2π f

 

 

 

 

 

 

 

 

 

 

 

 

 

x(t )=

 

 

 

B

x(ω)exp( jωt )dω =

 

 

 

 

 

2π

 

 

 

 

 

 

 

 

 

 

 

 

 

2π fB

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

2π fB

1

 

 

 

 

n

 

 

 

 

n

 

 

 

 

=

 

 

 

 

 

x

 

exp

jω

exp{jωt}dω =

2π

 

 

 

 

 

 

2 fB

 

 

 

2πfB

2 fB n=−∞

 

 

 

 

 

 

2 fB

 

 

 

 

 

 

1

 

 

 

 

 

 

n

 

 

 

2π fB

 

 

n

 

 

 

 

 

 

−∞

 

 

 

 

 

 

 

 

 

 

=

 

 

x

 

 

 

 

 

 

exp jω t

 

 

dω.

 

 

4π f

B

 

2 f

B

 

2 f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2π f

B

 

 

 

B

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Вычислив интеграл

2π fB

 

 

 

 

 

 

 

 

exp jω(t n / 2 f

B

)

 

2π fB

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(

t

n / 2 f

B )

 

 

 

 

 

 

 

 

 

=

 

 

exp jω

 

dω =

j (t n / 2 fB )

 

 

 

 

 

2π f

B

 

 

 

 

 

 

 

 

 

 

 

2π fB

 

 

 

 

 

(t n / 2 f

 

 

 

 

 

 

 

 

 

 

exp j2π f

B

B

) exp j2π f

B

(t n / 2 f

B

)

=

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

j

(t n / 2 fB )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= 2sin 2π fb (t n / 2 fB ) t n / 2 fB

и, подставляя результат в выражение (5.18), получим

 

n

 

2sin

2π fB (t n / 2 fb )

x(t )= x

 

 

 

 

 

.

 

 

 

 

n=−∞

 

2 fB

 

 

(t n / 2 fb )

(5.17)

(5.18)

(5.19)

210

Соотношение (5.19) решает проблему интерполирования функции x(t ),

т.е. задачу определения этой функции всюду по её значениям в точках отсчета n / 2 fB .

Разложение (5.19) для класса целых функций конечной степени было известно математикам как разложение по формуле Лагранжа очень давно, и работы Уиттекера и Гельфонда в этой области выполнены в 1915, 1935 и 1937 годах [8, 16, 25]. В этих, чисто математических исследованиях, были установлены причины, порождающие замечательное свойство функций с ограниченным спектром – свойство однозначности восстановления значений функций на всей оси времени по её значениям в некоторой периодической последовательности временных отсчетов. Прием, использованный выше для доказательства теоремы отсчетов, не использовал теорию интерполяции целых функций, поэтому его использование для получения более общих результатов невозможно. Однако для решения технических задач полученный результат вполне достаточен. Идея использования теоремы отсчетов в задачах связи была развита Найквистом и Хартли в 1924 1928 годах. В связи с этим величину 1/ 2 fB в западной литературе называют интервалом Найквиста. В дальнейшем теорема отсчетов была независимо доказана Котельниковым в 1933 году, а также исследована Шенноном в

1949 году [7, 8, 16, 25].

Вернемся теперь к выводам, следующим из теоремы отсчетов. Было установлено, что для неискаженной передачи функции с ограниченным спектром достаточно передавать ее значения в моменты времени, следующие с интервалом t =1/ 2 fB . Выборочные значения могут быть получены мето-

дом стробирования с использованием схемы электронного ключа, открывающегося на короткое время под воздействием импульсов стробирования, следующих с интервалом 1/ 2 fb (см. рис. 5.7)

211