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

advertisement
Спектральный анализ сигналов
1
Спектральный анализ сигналов
1. Непрерывное преобразование Фурье
Спектральное представление сигналов основано на разложении функций в ряд.
0
t
T
Рис.1. Непрерывная периодически повторяющаяся функция (сигнал) U(t).
Периодически повторяющаяся функция U(t) (далее - сигнал) любой формы с
периодом повторения T на бесконечном интервале времени может быть представлена
бесконечной суммой элементарных тригонометрических функций с надлежащим образом
подобранными параметрами: амплитудой Ak и начальной фазой φk:
a

U (t )  0   A cos(k t   )
1
k
2 k 1 k
(1)
или в форме:
a

U (t )  0   (ak cos(k1t )  bk sin( k1t ))
2 k 1
(2)
Преобразование выражения (1) в форму (2) основано на представлении
гармонического колебания в виде двух квадратурных составляющих: косинусной и синусной
с нулевыми начальными фазами.
Это разложение (1) или (2) периодической функции (сигнала) в бесконечный ряд
тригонометрических функций называется рядом Фурье, а функции ряда называются
гармониками.
Частоты гармоник кратны основной частоте ω1:
1  2  
T
f1  1
T
- круговая частота (рад/сек)
- циклическая частота (герц)
Коэффициенты (амплитуды гармоник) ряда (2) могут быть вычислены:
T
2 2
a 
U (t )dt
0 T T
2
T
2 2
a   U (t ) cos(k t )dt
k TT
1
2
T
2 2
b   U (t ) sin( k t )dt
k TT
1
2
(3)
Спектральный анализ сигналов
2
Амплитуды Ak и фаза φk гармоники ряда (1) связаны с коэффициентами ak и bk ряда
(2) связаны соотношениями:
b
A  a2  b2
k
k
k
  arctg k
k
a
k
Ряд Фурье обычно принято представлять в комплексной форме.
Для преобразования выражения (2) в комплексную форму следует воспользоваться
формулами Эйлера (представление тригонометрических функций экспоненциальными):
e j  e  j
cos( ) 
2
e j  e  j
sin(  ) 
2j
Тогда разложение U(t) в ряд может быть представлено:

a
U (t )  0   (a
2 k 1 k
e
jk t
 jk t
1 e
1
2
 j b
k
e
jk t
 jk t
1 e
1
)
2
или в форме:
a
 a  jbk jk t
 a  jbk  jk t
1   k
1
U (t )  0   k
e
e
2 k 1
2
2
k 1
(4)
Обозначим коэффициенты ряда (4)
a  jb
k
C  k
k
2
C
a  jb
k
 k
k
2
С учетом соотношений (3) коэффициенты Ck и C-k вычисляются:
T
a  jb
1 2
k
k
C 

 U (t )(cos(k1t )  j sin( k1t )dt
k
2
T T
2
T
a  jb
2
1
k 
C  k
U (t )(cos(k t )  j sin( k t )dt

1
k
1
2
T T
2
Так как
cos(α) ± jsin(α) = e±jα
то соотношения (5) в комплексной форме будут иметь вид
 jk t
1 2
1 dt
C 
U
(
t
)
e

k
T T
T
2
 jk t
1 2
1 dt
C 
U (t )e

k
T T
T
2
(5)
Спектральный анализ сигналов
Легко заметить, что значения C-k отличаются от Ck лишь знаком показателя
экспоненты.
Если ввести отрицательные значения k и учесть, что
a0
 C0
2
то разложение в ряд U(t) (4) можно представить в следующей форме:
jk t  
jk t
jk t

1
1
1
U (t )   C e
  C e
  C e
k
k
k
k 0
k  1
k  

И окончательно запишем формулы преобразования Фурье для периодически
повторяющихся сигналов с периодом повторения T:
 jk t
1 2
1 dt
C 
U
(
t
)
e
k T  T
T
- прямое преобразование Фурье
(6)
- обратное преобразование Фурье
(7)
2
U (t ) 


k  
C e
k
jk t
1
где Сk комплексные гармоники – спектр периодически повторяющегося сигнала.
Еще раз обратим внимание на свойства спектра периодически повторяющегося
сигнала на бесконечном интервале времени:
- спектр является дискретным;
- гармоники спектра кратны основной частоте f = 1/T (ω = 2π/T);
- число гармоник бесконечно;
- математически спектр содержит как реальные – положительные по частоте
гармоники, так и отрицательные гармоники (k – отрицательные значения).
Спектр (а точнее спектральная плотность) не периодического, одиночного, сигнала
может быть получен путем предельного перехода при T  
2
   d  0
T
k  
1
1 
C  S ( j )
k
- интервал между соседними гармониками
ω - непрерывная частота
- спектральная плотность
и преобразование Фурье для одиночного сигнала выполняется по формулам:

S ( j )   U (t )e
 jt
dt
(8)

U (t ) 
1 
jt
d
 S ( j  )e
2  
(9)
3
Спектральный анализ сигналов
4
Огибающая дискретного спектра периодического сигнала полностью совпадает со
спектральной плотностью одиночного сигнала. Поэтому для получения спектра
периодического сигнала с периодом T достаточно по (8) вычислить спектральную плотность
для одиночного сигнала и взять дискретные его значения через f = 1/T.
Пример:
Спектральная плотность прямоугольного импульса длительностью τ и амплитудой Um

2
S ( j )   U  e
m

 jt

2
dt  U m  e

2
 jt
U
dt  m e  jt
 j
2

2

Пределы
интегрирования
2
После подстановки пределов интегрирования и с учетом формул Эйлера получим
S ( j ) 
2 U

j
m e

2 e
2j
j

2
 U m  
sin(

)
2

2
1
   

 2 
sin
 
2
0.5
0

Рис.2. Модуль функции sin(x)/x (амплитудный спектр прямоугольного импульса)
Как следует из приведенного графика, функция sin(x)/x имеет лепестковый характер и
принимает значения, равные 0 при x = 2nπ n = 1, 2, 3…
или для спектра при значениях частоты f = 1/τ, 2/τ, 3/τ…
Уровень лепестков относительно главного составляет:
- первый - 0.217
- второй - 0.13
- третий - 0.09
- четвертый – 0.07
т.е., хотя в главном лепестке спектра и сосредоточена основная мощность прямоугольного
импульса (более 90%), все же спектр (спектральная плотность) его достаточно медленно
убывает с частотой.
Спектральный анализ сигналов
5
2. Дискретное преобразование Фурье (ДПФ)
Дискретное преобразование Фурье основано на дискретизации непрерывных
сигналов.
Пусть спектр сигнала U(t), заданного на интервале T, ограничен верхней частотой FВ.
Тогда в соответствии с теоремой Котельникова (теоремой отсчетов Винера) такой
сигнал может быть представлен своими дискретными значениями Un (n=0,1…N-1),
взятыми через интервалы времени
Δt=1/2FВ
Число дискретных значений – отсчетов сигнала будет равно N=T/Δt.
Тогда при условии дискретизации сигнала соотношение для расчета гармоник спектра
 jk t
1 2
1 dt
C 
U (t )e

k T T
T
2
следует записать в виде суммы, заменив
dt  t
t  N  t
T  N  t
1  2   
1
N  t
- текущее время
N 1
1
1
ck 
 U n  exp(  j  k  (2   
)  n  t )  t
N  t n 0
N  t
окончательно дискретное преобразование Фурье принимает вид
1 N 1
1
ck   U n  exp(  j  2     k  n)
N n 0
N
(10)
Рассмотрим некоторые свойства дискретного преобразования Фурье.
1. Сравним гармоники при k=0 и k=N; при k=1 и k=N+1 т.е отстоящих на N:
1 N 1
c0    U n
N n 0
1 N 1
c N   U n  exp(  j  (2    n)
N n 0
Т.к. индексы отсчетов n –целочисленные значения
то exp(-j2πn)=cos(2πn) – j*sin(2πn)=1
и cN=c0
1 N 1
1
c1   U n  exp(  j  2     n)
N n 0
N
c N 1
1 N 1
1
  U n  exp(  j  2     ( N  1)  n) 
N n 0
N
1 N 1
1
  U n  exp(  j  2     n)  exp(  j  2    n)
N n 0
N
То есть c1=cN+1
Таким же образом можно показать, что c2=cn+2 и т.д.
Спектральный анализ сигналов
6
Т.е. при дискретном преобразовании Фурье число рассчитываемых гармоник равно
числу отсчетов N. (Далее значения гармоник повторяются).
2. Рассмотрим гармоники c1 и cN-1
гармоники cN/2
c N 1
c2 и cN-2, т.е. симметричные относительно
1 N 1
1
  U n  exp(  j  2     ( N  1)  n) 
N n 0
N
1 N 1
1
  U n  exp(  j  2     n)  exp(  j  2    n)
N n 0
N

1 N 1
1
 U n  exp(  j  2     n)  c1*
N n 0
N
Также и
c N 2  c2*
и т.д
Т.е. гармоники, симметричные относительно cN/2, являются комплексно сопряженными.
Обратное преобразование Фурье выполняется по формуле:
N 1
U n   ck  exp(  j  2   
k 0
1
 k  n)
N
(11)
Для выполнения дискретного преобразования Фурье (прямого и обратного) при
прямом использовании соотношений (10) (11) потребуется N*N операций комплексных
умножений и сложений (операция умножения требует больших временных или
аппаратурных затрат по сравнению с операций суммирования).
Для уменьшения временных затрат для выполнения дискретного преобразования
Фурье разработаны различные алгоритмы быстрого преобразования Фурье, основанные на
прореживании данных по времени или по частоте. Для БПФ число отсчетов сигнала N
должно быть степенью числа 2. Тогда необходимое число операций для полного
преобразования Фурье составит N*log2N.
В пакете Mathcad имеются строенные функции для ДПФ:
c:=FFT(U)
U:=IFFT(c)
FFT - прямое быстрое преобразование Фурье (БПФ)
Преобразуемая переменная U должна быть представлена вектором вещественных
чисел размерностью N=2m m>2. Результат преобразования FFT – вектор комплексных чисел
размерностью 1+2m-1 (иными словами FFT вычисляет только часть возможных гармоник ck
k:=0…N/2 без комплексно сопряженных гармоник – см. свойство 2)
IFFT - обратное быстрое преобразование Фурье
Здесь c – вектор комплексных чисел (гармоники) размерностью 1+2m-1 m>2.
IFFT –возвращает вектор вещественных чисел размерностью N=2m
Многомерное преобразование Фурье
G:=CFFT(A)
A:=ICFFT(G)
Данные преобразования применимы как к векторам, так и к матрицам комплексных
чисел. Как прямое преобразование CFFT, так и обратное ICFFT возвращает вектор или
матрицу той же размерности, что и преобразуемый вектор или матрица. При этом не
накладывается ограничений типа равенства размерности строк или колонок степени числа 2.
Спектральный анализ сигналов
7
Комплексное преобразование Фурье CFFT может быть применено и вектору
вещественных чисел U вместо FFT, если, к примеру, не удается обеспечить требование
равенства N=2m. Однако, если над гармониками затем выполняется некое преобразование,
например, моделируется прохождения сигнала через частотно зависимое устройство в виде:
S ВЫХ ( j)  S ВХ ( j)  K ( j)
то в силу комплексной сопряженности гармоник, симметричных относительно cN/2, обратное
преобразование ICFFT уже не обеспечит получение временной функции (вида сигнала) на
выходе.
Ниже приведены амплитудный и фазовый спектры сигнала (прямоугольного
импульса), полученные с использованием функций FFT и CFFT.
arg c k
ck
k
k
а)
N
2
ck
k
N
2
arg c k
k
б)
Рис.2. Амплитудный и фазовый спектры, полученные с помощью FFT (а) и CFFT (б).
Кроме того, в Mathcad предусмотрены встроенные функции
c:=fft(U)
U:=ifft(c)
которые выполняются по формулам
1 N 1
1
ck 
 U n  exp(  j  2     k  n)
N
N n 0
Un 
1 N 1
1
  ck  exp(  j  2     k  n)
N
N k 0
Как видно из приведенных формул, гармоники спектра, рассчитанные по fft, по
амплитуде в 1/√N будут больше и комплексно сопряжены по сравнению с гармониками,
рассчитанными по FFT.
Спектральный анализ сигналов
8
(Соответственно, следует пользоваться парами FFT – IFFT или fft - ifft)
Аналогично действуют и функции комплексного преобразования Фурье cfft – icfft.
Примечание:
При описании формирования сигналов использованы обозначения переменных
T:=
t:=0…T-1
и сигнал рекомендуется сразу задать в виде вектора Ut, т.е. сигнал фактически представлен в
виде дискретных отсчетов для последующего дискретного преобразования Фурье. Иными
словами, значению числа отсчетов N, использованному в описании выше соответствует
значение T, а индексам n значения t.
3. Примеры расчета спектров некоторых сигналов
3.1. Прямоугольный импульс.
0.1
0.2
c1 k
c2 k
0
0
10
20
0
30
0
10
20
k
k
а)
б)
30
0.4
c3 k
0.2
0
0
10
20
30
k
в)
Рис.3. Амплитудные спектры прямоугольного импульса.
а – длительность импульса T/10
б – длительность T/5
в – длительность T/2 (меандр).
Как уже отмечалось выше, огибающая спектра описывается функцией |sin(x)/x|.
Спектр периодической последовательности импульсов с τ=T/2 (меандра) не содержит
четных гармоник.
3.2. Прямоугольный радиоимпульс.
0.05
ck
0
80
100
120
k
Рис.4. Спектр радиоимпульса длительностью T/10 и частотой заполнения 100/T.
Спектр радиоимпульса симметричен относительно частоты заполнения и каждая
боковая полоса соответствует спектру видеоимпульса.
Спектральный анализ сигналов
3.3. Импульс с экспоненциальными фронтами.
0.2
1
Ut
ck
0
0
500
0
1000
0
10
20
t
30
k
Рис.5. Импульс с экспоненциальными фронтами (длительностью T/5) и его спектр.
При уменьшении крутизны фронтов импульса уменьшается уровень боковых
лепестков спектра – следующих за главным.
3.4. Колоколообразный (гауссов) видеоимпульс
1
0.1
Ut
ck
0
0
500
0
1000
0
10
t
20
k
Рис.6. Гауссов импульс длительностью T/10 на уровне 0.5 и его спектр.
Огибающая спектра также описывается гауссовой зависимостью. Спектр гауссового
импульса является самым компактным (сравни со спектром прямоугольного импульса
длительностью T/10 на рис.3).
3.5. Случайный телеграфный сигнал с фазовой модуляцией
0.1
c_фм k
0
200
300
400
500
600
700
800
k
Рис.7. Спектр случайного телеграфного сигнала с ФМ (число элементов 64).
9
Download