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

1490

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

связи с этим, при решении практических задач необходимо оценить полосу частот исследуемого сигнала, в которой сосредоточена основная часть энергии сигнала. Под основной частью энергии понимают обычно величину не менее чем 0.9 полной мощности сигнала, а полосу частот, в которой сосредоточено 90% мощности сигнала, называют эффективной шириной спектра сигнала.

Определим эффективную ширину спектра для трех типов сигнала, рассмотренных в п. 2.3, а именно: прямоугольный импульс, экспоненциальный импульс и функция Гаусса [11, 12].

2.4.1. Эффективная ширина спектра прямоугольного импульса

Напомним прежде всего, что при анализе спектральной плотности прямоугольного импульса, определенной выражением (2.10а)

ω =τ sinωτ 2 x( )

ωτ 2

было рассмотрено соотношение между длительностью импульса τ и полушириной главного лепестка его спектральной плотности. При этом полушириной главного лепестка был назван интервал угловых частот 0 <ω 2πτ . Представляется чрезвычайно интересным определить часть энергии прямоугольного импульса, содержащейся в данной полосе частот, поскольку именно для этой полосы частот был установлен принцип неопределенности, характеризуемый выражением f τ =1. Учитывая, что амплитуда импульса была принята единичной, нетрудно установить, что полная энергия P прямоугольного импульса длительностью τ определяется площадью импульса, т.к.

62

τ 2

P = 12dt =τ.

τ 2

Поскольку, в соответствии с теоремой Парсеваля, полная энергия определяется как

2

 

1

&

 

2

 

P =

 

 

 

x

(t)dt =

 

 

 

dω,

2π

x(ω)

 

 

−∞

 

 

−∞

 

 

 

 

а для симметричного энергетического спектра это выражение может быть записано в виде

 

 

 

P =

1

 

 

X& (ω)

 

2 dω ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

то энергию, заключенную в полосе частот ω [0, ωВ], где

ωВ

верхняя

частота, можно вычислить, используя соотношение

 

 

 

 

 

 

 

1

 

ω

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Pω =

 

 

 

В

 

X& (ω)

 

2dω.

 

 

 

 

 

 

 

 

 

 

 

 

π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

Тогда, используя спектральную функцию (2.10а), запишем

 

 

 

 

τ2 ωВ sin2 (ωτ / 2)

 

 

 

 

Pω = π

0

 

 

 

 

dω.

 

 

(2.16а)

 

 

 

(ωτ / 2)2

 

 

Используя замену переменных ωτ

 

 

2 = x (dω = 2x τ),

перепишем выраже-

ние (2.16а) в виде

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ωВτ 2

2

 

x dx =τη(ωВτ 2),

 

 

 

Pω =τ π2

 

0

sinx2

 

 

 

 

(2.16б)

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

 

 

 

 

 

 

 

 

2

 

 

ωВτ 2

2

 

 

 

 

 

 

η(ωВτ / 2) =

 

 

 

0

sinx2 xdx

 

 

(2.17)

π

 

 

 

 

 

определяет часть энергии

сигнала

заключенную

в

полосе

частот

ω [0, ωВ].

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

63

Интеграл, входящий в выражение (2.17), после интегрирования по час-

 

=sin

2

x / x

2

, du =(2x

2

sin x cos x 2xsin

2

x)/ x

4

 

 

 

 

 

приво-

тям u

 

 

 

 

 

 

 

, du = dx, v = x

дится к виду

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 ωВτ 2 sin2 x

 

 

 

2

sin2

x

 

ωВτ 2

 

ωВτ 2 2sin xcos x

 

 

 

ωВτ 2 sin2 x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x2

 

dx =

 

 

x

 

 

 

 

 

 

x

 

 

dx + 2

x2

dx ,

 

π

 

 

 

π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

откуда

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

0

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 ωВτ 2 sin2 x

 

 

 

2

 

 

ωВτ 2 sin 2x

 

sin2 x

 

ωВτ 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x2

 

dx

=

 

 

 

 

 

x

dx

 

x

 

 

 

.

 

 

 

(2.18)

 

 

 

 

π

 

 

 

π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Интеграл в

 

 

правой части

 

выражения

(2.18) после

замены

переменных

2x =t ( x =t

 

 

2, dx = dt

2 )

может быть записан как

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z sin tdt =Si(z),

 

 

 

 

 

 

 

 

 

(2.19)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где z =ωВτ . Функция (2.19) представляет собой одну из форм интеграль-

ного синуса и подробно табулирована (см. например таблицы А. Анго [3] ). Используя (2.19), перепишем выражение (2.18), определяющее часть энергии, заключенную в полосе [0,ωВ] , в окончательном виде

 

2 ωВτ

2 sin2 x

dx =

2

 

sin2 (ωВτ 2)

 

 

 

 

 

 

2

 

Si(ωВτ)

 

 

.

(2.20)

 

π 0

x

π

ωВτ 2

 

 

 

 

 

 

 

 

Численный расчет, выполненный с интервалом

0,1π

для переменной

ωВτ2 , дает результат, представленный на рис 2.9. Из этого рисунка вид-

но, что в точке ωВτ

2 =π

( f τ =1), т.е. в полосе частот от 0 до fВ =1 τ

сосредоточено

90%

всей

энергии импульса. В связи с этим формула

f =1 τ дает

простейшее

выражение для выбора полосы пропускания

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

64

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

Рис. 2.9. К определению эффективной ширины спектра

2.4.2. Эффективная ширина спектра экспоненциального импульса

Спектральное разложение

экспоненциального импульса вида

α = exp(αt ), где α =1 τ, t 0 было найдено в п. 2.3.2 в виде

X&

1

.

(ω) =

 

α + jω

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

X&

(ω)

 

2

=

1

.

(2.21)

 

 

 

 

α2 +ω2

 

 

 

 

 

 

 

 

Определим теперь эффективную ширину спектра (2.21), т.е. полосу частотой [ω, ωВ ], в которой сосредоточенно 90% энергии импульса. Для определения эффективной ширины спектра воспользуемся соотношением

0.9

dω

 

1

ωВ

dω

 

π

 

=

 

 

 

.

(2.22)

α2 +ω2

 

π

α2 +ω2

 

0

 

 

 

 

0

 

 

 

Тогда

65

0,9 arctg

ω

 

=

1

arctg

ω

 

ωВ

 

 

 

 

 

 

 

 

πα

α

 

0

πα

 

α

 

0

 

 

 

Подставляя пределы интегрирования, получим arctg ωαВ = 0,45π ,

откуда верхняя частота эффективной ширины спектра определяется, как

ωВ = tg(0,45π)α = 6,28/τ .

Учитывая, что ωВ = 2π fВ , нетрудно увидеть, что верхняя частота эф-

фективной ширины спектра экспоненциального есть

fВ = ωπВ =τ1

2

Таким образом, 90% энергии экспоненциального импульса сосредоточено в полосе частот от 0 до fВ =1τ . Данный результат совпадает с опре-

делением эффективной полосы частот прямоугольного импульса (см. под-

раздел 2.4.1).

2.4.3. Эффективная ширина спектра гауссова импульса

Спектральная плотность гауссова импульса x(t) = exp(t2 2τ2 ) опреде-

лена выражением (2.14)

X (ω) = 2πexp ω2τ2 .

2

Полная энергия гауссова импульса может быть найдена в виде

P = x2 (t)dt =

exp (-t2 /τ2 )dt

(2.24)

−∞

−∞

 

Используя замену переменных t /τ = x , dt = dxτ , преобразующую выражение (2.24) в интеграл Эйлера – Пуассона, определим полную энергию импульса как

66

P =τ π.

(2.25)

Найдем теперь эффективную ширину спектра гауссова импульса, в пределах которой сосредоточено 90% энергии импульса. Для этого проинтегрируем энергетический спектр гауссова импульса

X& (ω) 2 = 2πτ2exp(ω2τ2 )

от нуля до некоторой частоты ωВ , и найдем эту частоту, исходя из условия

90 – процентного содержания энергии импульса в полосе [0, ωВ]

Определим энергию гауссова импульса в полосе [0, ωВ] как

 

1

ω

 

 

 

 

 

 

 

ω

 

Pω =

В

2πτ2exp(ω2τ2 )dω = 2τ2 В exp(ω2τ2 )dω .

(2.26)

π

 

 

0

 

 

 

 

 

 

 

0

 

Преобразуем соотношение (2.26), используя интеграл вероятностей [9]

 

 

 

 

2

z

 

(

 

2

)

 

 

 

 

Φ(z) =

 

exp

x

 

(2.27)

 

 

 

π

 

 

 

dx ,

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

который табулирован во многих справочниках, например [13].

Для преобразования выражения (2.26) используем замену переменных ωτ = r , dω = dr / dτ и перепишем это выражение как

ω

 

ωВτ

 

(

 

)

2

z=ωВτ

 

(

 

)

 

 

x2

 

x2

P

= 2τ

0

exp

 

dr =τ π

π

0

exp

 

dr =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=τ

πΦ(z) = PΦ(z) .

 

 

 

(2.28)

Из выражения (2.28) следует, что отношение энергии Pω , заключенной в полосе частот [0, ωВ], к полной энергии гауссова импульса P равно

Pω P = Φ(z =ωВτ).

(2.29)

Для 90 – процентного содержания энергии импульса в полосе [0, ωВ]

левая часть выражения (2.29) равна 0.9, что и позволяет, используя таблицы интеграла вероятностей, определить величину аргумента Z =ωВτ 1,16.

67

Учитывая, что τ представляет собой половину длительности гауссова импульса на уровне, 1/ e 0,606 т.е. τ = 0,5t , видим, что π fВt =1,16

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

ввиде fв 037 / t .

2.5.Основные свойства преобразования Фурье

(теоремы о спектрах)

Рассмотрим несколько общих теорем о спектрах, основанных на свойствах преобразования Фурье [1, 3, 7, 8, 11].

2.5.1. Теорема сложения (теорема линейности)

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

Преобразование Фурье взвешенной суммы (разности) сигналов равно взвешенной сумме (разности) спектров этих сигналов:

N

fi

 

N

(2.30)

∫ ∑ai

(t) exp(jωt )dt = ai Si (ω) .

−∞ i=1

 

 

i=1

 

Здесь fi (t) и Si (ω) есть функции времени и отвечающие им спектры, а ai

есть некоторые весовые коэффициенты.

Сложение комплексных спектральных плотностей производится по правилам сложения комплексных чисел

N

 

 

N

2

 

N

 

2

0,5

 

 

ai S&i

(ω)

=

ai

Re S&i (ω)

+

ai

Im S&i

(ω)

.

 

 

 

i=1

 

 

i=1

 

 

 

 

i=1

 

 

 

 

 

 

 

 

 

 

 

 

Спектр амплитуд и спектр фаз суммы спектров определяется как

68

N

N

 

 

N

 

ai S&i (ω) = ai Re S&i (ω) + jai Im S&i (ω) ,

i=1

i=1

 

i=1

 

 

 

 

N

 

 

N

 

 

ai Im S&i (ω)

 

arg ai S&i (ω)

= arctg

i=1

 

.

N

&

i=1

 

 

ai

 

 

 

 

Re Si (ω)

 

 

 

i=1

 

 

2.5.2. Теорема смещения (теорема задержки)

Рассмотрим влияние временного смещения (задержки) сигнала на его спектр. Анализ проведем с использованием конкретного примера.

Пусть имеют место два сигнала в виде прямоугольных импульсов длительностью τ , обладающих единичной амплитудой. Однако их расположение на оси времени различаются: один из сигналов расположен симметрично относительно начала координат, а второй смещен относительно первого на величину τ / 2 (см. рис 2.10, а, б).

x(t )

τ / 2

0

t

τ / 2

Рис. 2.10а. Прямоугольный импульс x(t )

0

τ

t

 

Рис. 2.10б. Смещенный прямоугольный импульс

69

Спектр несмещенного сигнала определяется выражением (2.10а)

X& (ω) =

x(t)exp(jωt )dt =τ sin(ωτ 2) .

−∞

ωτ 2

Модуль X& (ω) и фаза ϕ(ω) этого спектра изображены на рис 2.11, а, б.

Рис. 2.11а

ϕ(ω)

π

ω

π

Рис. 2.11б

ϕ(ω)

ω

Рис. 2.11в

Рис. 2.11. Амплитудный и фазовый спектры смещенной функции

70

Рассмотрим теперь выражение для комплексного спектра сигнала, отличающегося от первоначального сигнала произвольным смещением (запаздыванием) на время t

X&

1(ω) = x(t − ∆t)exp(jωt )dt .

(2.31)

 

−∞

 

Используя замену переменных t − ∆t =t1 , t =t1 + ∆t , dt = dt1 , перепишем

(2.31) в виде

X&1(ω) = x(t1)exp[jω(t1 + ∆t)]dt1 =

−∞

= exp(jωt ) x(t1)exp(jωt1 )dt = X& (ω)exp(jωt ), (2.32)

−∞

где X& (ω) есть спектр несмещенной функции. Таким образом, теорема о спектре смещенной функции может быть сформулировать следующим образом:

Если функция x(t) , обладающая спектральной плотностью X& (ω) ,

смещена на время t , то спектр смещенной функции определяется произ-

ведением спектра несмещенной функции X& (ω) и фазового множителя exp(jωt) :

X&1(ω) = X& (ω)exp(jωt) . (2.33)

Для конкретного примера смещенного прямоугольного импульса (см. рис. 2.10, б), получим:

X&1(ω) = τ0 exp(jωt )dt =− j1ω exp(jωt )τ0 =

= − j1ω exp(jωt )1 = j1ω[1cosωτ + j sinωτ].

Используя выражения sin2 α2 = 0,5(1cosα) , sin 2α = 2cosαsinα пре-

образуем последнее соотношение к виду

71

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]