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

ПОСОБИЕ ПО ЛАБАМ СФУ

.pdf
Скачиваний:
63
Добавлен:
09.06.2015
Размер:
3.38 Mб
Скачать

2 T 2

bn = T òs(t)sin nω0tdt .

T 2

В общем случае периодический сигнал содержит постоянную составляющую и бесконечный набор гармонических колебаний гармоник, с частотами ωn = nω0, (n = 0, 1,..) кратными основной частоте ω0=2π/T исследуемой последовательности с известным периодом T. Четный сигнал содержит только косинусоидальные составляющие, а нечетный сигнал синусоидальные.

Каждую гармонику можно описать её амплитудой Аn и начальной фазой ϕn, тогда коэффициенты ряда Фурье принимают вид:

an = An cosϕn ;

 

bn =An sinϕn ;

 

 

 

 

 

 

=

bn

 

A = a

2 + b 2

;

tgϕ

n

;

 

n

n

n

 

 

 

an

 

 

 

 

 

 

 

 

 

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

s(t) = a0 + åAn cos(nω0t −ϕn ) . 2 n=1

Спектральная диаграмма периодического сигнала это графическое изображение коэффициентов ряда Фурье для конкретного сигнала.

Спектральные диаграммы бывают амплитудные и фазовые (рис. 11). Как следует из приведенных графиков, спектры периодических сигналов состоят из отдельных линий, то есть они дискретны.

An

ϕ n

 

1

3

w¤w0

1

3

w¤w0

 

 

 

 

Рис. 11. Амплитудная и фазовая спектрограммы периодического сигнала

1.3.2. Комплексная форма ряда Фурье. Спектральное разложение

периодического сигнала можно провести в системе экспоненциальных базисных функций. Функции этого базиса периодичны с периодом Т и ортонормированны на отрезке времени [-T/2,T/2]. Тогда комплексный ряд Фурье, с учетом нормы комплексного сигнала, характеризуется формулами:

s(t) = åcne jnω0t ,

 

n=−∞

cn =

1

T 2s(t)e jnω0t dt.

 

 

T Tò

 

2

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

ejnω0t = cos nω0t j sin nω0t ;

cosnω0t = ejnω0t + ejnω0t ,

2

j sin nω0t = ejnω0t ejnω0t .

2

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

1.3.3. Спектральное представление непериодических сигналов

Предположим, что имеется одиночный импульсный сигнал s(t) конечной длительности. Мысленно дополним его такими же сигналами, периодически следующими через некоторые интервалы времени Т. В результате получим периодическую последовательность sпер(t), которую можно представить в виде ряда Фурье (23).

Для перехода к одиночному импульсу увеличим до бесконечности период повторения импульсов Т. В этом случае:

1). Частоты nω0 и (n+1)ω0 окажутся сколь угодно близкими и поэтому дискретную переменную nω0 можно заменить непрерывной – ω, то есть текущей частотой.

2). Амплитудные коэффициенты cn станут бесконечно малыми (стремящимися к нулю) из-за наличия Т→∞ в знаменателе в формулах ряда.

3). Рассмотрим интервал частот Δω→0 в окрестностях некоторой частоты ω0. В пределах этого интервала будет содержаться N отдельных пар спектральных составляющих, частоты которых будут отличаться друг от друга сколь угодно мало (N=Δω/ω0=ΔωT/2π).

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

Комплексная амплитуда эквивалентного гармонического сигнала внутри интервала Δω:

Aω0

=

2N

j

t

dt =

ω

j

t

dt .

T

òs(t)e

ω0

π

òs(t)e− ω0

 

 

 

−∞

 

 

 

−∞

 

 

 

 

 

 

 

 

 

 

 

 

 

При этом интеграл:

òs(t)ejωt dt = S( jω)

−∞

называется спектральной плотностью сигнала s(t) или прямым преобразованием Фурье данного сигнала.

С точки зрения физического смысла спектральная плотность S(ω0)=S(2πf0) – масштабный множитель, связывающий малую длину интервала частот f и соответствующюю комплексную амплитуду A гармонического сигнала на центральной частоте f0.

При решении обратной задачи, то есть при нахождении сигнала по его спектральной плотности, необходимо воспользоваться обратным преобразованием Фурье для сигнала s(t):

s(t) = 1 òS(ω)e jωt dω .

2π −∞

Следует отметить основное условие существования спектральной плотности сигнала: для того, чтобы сигналу s(t) можно было бы сопоставить его спектральную плотность S(ω), необходимо, чтобы сигнал был абсолютно интегрируем, то есть, чтобы существовал интеграл:

òs(t) dt < +∞ .

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

Приложение 2

К лабораторной работе №2

Дискретизация информационных сигналов

Задача преобразования непрерывных функций времени в дискретные встречается часто в радиотехнической практике: при организации

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

В зависимости от ограничений, накладываемых на информационный сигнал и его спектральные характеристики, преобразование непрерывного сигнала s(t) в дискретно-непрерывный s(k t) может быть выполнено на основании частотного критерия дискретизации, корреляционного, критерия наибольшего отклонения и др.

Частотный критерий выбора шага дискретизации основан на теореме КотельниковаШеннона в которой отмечено, что непрерывная функция времени s(t), удовлетворяющая условиям Дирихле и не содержащая частот выше ωm, полностью определяется отсчетами мгновенных значений s(k t), отстоящих друг от друга на интервалах t=π /ωm, где k=1, 2, 3 ...

В соответствии с этой теоремой отсчетов, функцию s(t), удовлетворяющую условиям теоремы, можно представить рядом в виде:

s(kDt)sinωm (t-kDt)

 

s(t) = å

(1)

k=−∞

ω (t-kDt) .

 

m

 

Каждое слагаемое в формуле (1) представляет собой произведение функции времени s(t), взятой в дискретные моменты времени (k t), на непрерывную функцию времени, называемую функцией отсчетов. В моменты времени tk=k t все слагаемые, кроме одного, обращаются в ноль. Внутри промежутка t восстанавливаемая функция определяется всеми слагаемыми. Функция s(t) с ограниченным спектром полностью задается счетным множеством мгновенных значений, определенных через равные промежутки времени t – шаг дискретизации.

Теорема КотельниковаШеннона применима к непрерывным функциям s(t), не ограниченных во времени. Только в этом случае функция s(t) может иметь ограниченный спектр частот. Реальные физические процессы имеют начало и конец, а функции, их описывающие, ограничены во времени и, следовательно, имеют неограниченные спектры [2]. Введением понятия реальной ширины спектра как интервала частот, в котором сосредоточена преобладающая часть энергии функции s(t), сохраняет общий смысл теоремы о представлении функции на интервале [0, T] счетным числом значений:

n = ωmT π .

Но значение всех n отсчетов на интервале [0, T] не позволяет полностью определить функцию s(t), т. к. определяющие функцию точки, лежащие вне интервала [0, T], тоже влияют на восстанавливаемый сигнал.

Ограничение числа членов суммы ведет к появлению погрешности

восстановления функции, производимого с помощью ряда КотельниковаШеннона, равную:

n

 

sinωm(t-kDt)

 

ε(t) = s(t) - å

s(kDt)

 

.

(2)

ω (t-kDt)

k=-n

 

m

 

В точках tk=k t погрешность восстановления минимальна и равна нулю. Внутри интервала t погрешность отлична от нуля, причем при приближении к границам интервала [0, T] она возрастает (рис. 1).

Определим величину среднеквадратической погрешности восстановления дискретизированного сигнала таким идеальным фильтром, используя подобие процессов дискретизации и амплитудно-импульсной модуляции (АИМ). Дискретизированный сигнал (рис. 2) часто рассматривается как АИМ-сигнал, для которого несущим колебанием

является последовательность прямоугольных импульсов единичной амплитуды с периодом следования, равным шагу дискретизации, и бесконечно малой длительностью (рис. 3). Амплитудный спектр такого несущего колебания имеет вид, показанный на рис. 4. В свою очередь, спектр АИМ-сигнала (дискретизированного) примет вид, показанный на рис. 5. Если информационный сигнал удовлетворяет требованиям теоремы отсчетов, то

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

 

 

 

 

δ (t)

 

 

 

 

 

 

 

 

t

−3Δt

−2Δ t

− t

0

t

2Δ t

3Δ t

 

 

Рис. 1. График погрешности

 

 

 

s(k

t)

 

 

 

 

 

1

 

 

 

 

 

 

 

0

 

 

 

 

 

 

t

 

t

2Δ t

3Δ t

4Δ t

5Δ t

6Δ t

 

 

Рис. 2.

Дискретизированный сигнал

 

1

u(t)

 

 

 

 

 

 

0

 

 

 

 

 

 

t

 

t

2Δ t

3Δ t

4Δ t

5Δ t

6Δ t

 

 

Рис. 3. Несущее колебание

 

 

 

C k ( ω )

 

 

 

 

 

0

 

 

 

0

 

 

ω

 

ω 0

 

0

 

0

Рис. 4. Амплитудный спектр идеального импульсного АИМ-колебания

 

SАИМ(ω)

 

 

 

 

 

 

 

 

 

 

 

 

ω

 

0

ω 0

0

0

 

0

 

 

Рис. 5. Спектр дискретизированного сигнала

Процесс восстановления информационной функции идеальным фильтром (рис. 6 а) поясняется на рис. 6. Для обеспечения абсолютной

точности восстановления дискретизированного сигнала следует строго выдержать соотношение:

ω0 =

= 2ωm ,

ωc = ωm .

(3)

t

При нарушении последнего условия, как следует из рис. 6, возрастает

погрешность восстановления, т. к. при ωсm в полосу пропускания фильтра непопадают составляющие спектра информационной функции, лежащие в диапазоне [ωс, ωm]; при ωсm в полосу пропускания фильтра попадают составляющие спектра нижней боковой при ω0, лежащие в диапазоне частот

[ωm, ωc].

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

длительность его импульсов τ хотя бы на порядок меньше периода их следования:

Т= t, (скважность q=T/τ = t /τ >10).

Спектр АИМ-сигнала для идеальной информационной функции (идеальная информационная функция удовлетворяет условиям теоремы Котельникова-Шеннона) и реального несущего колебания приведен на рис. 7, период следования импульсов которого Т= t выбран в соответствии с теоремой t=π /ωm, поэтому спектры боковых составляющих не перекрываются. Это позволяет восстановить сигнал (информационную функцию) с абсолютной точностью.

Реальные сигналы имеют неограниченные спектры. Спектр АИМ- сигнала для случая такой информационной функции приведен на рис. 8 и 9 (на рис. 9 выделен фрагмент спектра). Из рис. 9 следует, что сигнал, восстановленный даже идеальным фильтром, отличается от исходного информационного, т. к.:

1) в полосу пропускания фильтра попали только лишь составляющие спектра сигнала в диапазоне частот [-ωc , ωc], в то время как сигнал имеет бесконечный спектр; ограничение частотного диапазона приводит к появлению среднеквадратической погрешности:

σ

2

 

1

S 2 (ω ) d ω

 

 

 

=

 

ò

,

(5)

0

π

 

 

ωc

 

где S2(ω) спектральная плотность мощности информационной функции;

2) в полосу пропускания фильтра попала часть составляющих спектра верхней боковой при (-ω0) и нижней боковой при ω0; появление этих

составляющих увеличивает среднеквадратическую погрешность на величину:

σ

2

 

1 ω

Sáî2

ê ,ω0 (ω)dω,

 

1

=

 

−ωò c

(6)

π

где S2бок,ω0 спектральная плотность мощности нижней боковой при ω0 в диапазоне [- ωc , ωc];

а )

K (ω )

c

 

0

 

 

 

 

б ) K m

K ( ω )

 

 

 

 

0 .

7 0 7 K m

 

 

ω c ω

0 .1 K m

ω c

ω

 

ω c

S А И М (ω )

ω

 

 

 

 

 

 

− ω 0

− ω m

0

m

m

ω 0

 

 

 

m

 

 

 

=ω

 

 

 

c

c

c

Р и с . 6 .

 

 

ω

ω

ω

В о с ста н о вл е н и е и н ф о р м а ц и о н н о й ф у н кц и и

и д е а л ьн ы м ф и л ьтр о м

S А И М (ω )

0

ω 0 2 ω 0 3 ω 0 4 ω 0 5 ω 0 6 ω

ω

0 7 ω 0 8 ω 0

 

-

 

 

2 π /τ u

Р и с . 7 . С п е ктр д и скр е ти зи р о в а н н о го (А И М ) с и гн а л а п р и

пр ям о у го л ьн о й и н ф о р м а ц и о н н о й ф у н кц и и и р е а л ьн о м и м пул ьс н о м н е су щ е м ко л е б а н и и

S А И М (ω )

 

ω 0 2 ω 0 3 ω 0 4 ω 0 5 ω 0 6 ω 0 7 ω 0

ω

0

8 ω 0

 

-

 

 

2 π /τ u

 

Р и с . 8 . С п е ктр д и скр е ти зи р о в а н н о го (А И М )

с и гн а л а п р и

р е а л ьн ы х и н ф о р м а ц и о н н о й ф у н кц и и и и м пу л ьс н о м н е су щ е м ко л е ба н и и

3) аналогично устанавливается составляющая погрешности, обусловленная попаданием части верхней боковой при (-kω0) и части нижней боковой при kω0 в полосу пропускания фильтра.

Суммарная среднеквадратическая погрешность восстановления дискретизированного сигнала может быть определена по формуле:

σ 2

 

σ 2

 

= å

(7)

 

k =

0

k .

Анализ формул (5) и (6) , а также рис. 9 и 10 приводит к выводу, что

σ 02

 

 

2 = 02 .

 

= å σ k2

;

σ

(8)

 

k =1

 

 

 

 

Последнее соотношение позволяет установить относительную среднеквадратическую погрешность δ%2 восстановления дискретизированного сигнала идеальным фильтром в следующем виде:

2

2

 

ò

S2

(ω)dω

 

é

ò

ωc

S

2

(ω)dω

ù

 

 

ωc

 

 

 

ê

 

 

ú

 

δ

= 0

 

 

 

 

 

0

 

 

 

 

= 2

ò0S2

(ω)dω

= 2

ê1-

 

 

ú .

(9)

òS2

(ω)dω

 

 

 

 

 

 

 

ë

 

0

 

 

 

û

 

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

 

 

 

1

 

T

2

 

 

1

 

 

 

R

 

=

 

ò

 

s2 (t )∂ t = C2

4 +

 

å C2

,

(10)

 

T

 

 

2

 

s

 

T 2

0

 

k=1

k

где С0, Сk коэффициенты ряда Фурье.

Первый вариант из (10) (определения средней мощности во временной области) реализуется проще т.к. не требует нахождения бесконечной суммы, но требуется информация о временной функции. Поскольку в полосу пропускания идеального фильтра попадают лишь m первых гармоник сигнала, то мощность сигнала ограниченного в частотной области полосой пропускания фильтра ωс=mΩ , равна:

 

 

1

m

 

Rs (mW) = C02

4 +

 

å Ck2 .

(12)

2

 

 

k=1

 

Рис. 9. Фрагмент спектра АИМ-сигнала

Рис. 10. Спектр реального сигнала и его

ограничение идеальным фильтром

Относительная среднеквадратичная погрешность восстановления дискретизированного периодического сигнала в этом случае будет равна:

 

 

 

 

 

æ

 

ö

δ

2

 

2

ç

Rs (mW ) ÷

=

= 2 ç1-

 

÷ .

 

0

Rs

 

 

 

 

 

è

ø

К лабораторным работам №5, 6

Краткие сведения о расчете спектра при нелинейных преобразованиях сигналов

(13)

Приложение 3

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