Добавил:
ИВТ (советую зайти в "Несортированное")rnПИН МАГА Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

TsOS_Ekzamen_2021_1

.pdf
Скачиваний:
80
Добавлен:
19.06.2024
Размер:
33.08 Mб
Скачать

Пример 5. Вырезка гармоники Φ0 = π/2. Сдвиг шкалы.

Так как φ0 = π/2, то исходный сигнал уже не косинус, а синус.

На амплитудный спектр это никак не повлияет, а вот весь фазовый спектр будет идентичен фазовому спектру из примера 4, но вдобавок к этому он будет поднят на значение φ0 = π/2.

10. Пример 6. Сигнал типа Sinc. Пример 7. Сигнал типа SInc с заполнением гармоникой. Техническая полоса сигнала ΔfT.

Пример 6.

Вместо того, чтобы брать интеграл, воспользуемся тем обстоятельством, что преобразование Фурье и обратное преобразование Фурье полностью симметричны, то

есть имеет смысл такая замена:

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

При этом заметим, что если в исходном сигнале типа sinc измерить значение на

 

(−

2∆ ; +

2∆ )

 

 

интервале:

 

1

1

, оно будет равным

1

, где

 

это полоса сигнала.

Пример 7.

Тот же сигнал sinc, только внутри заполнен гармоникой с частотой f0.

1)Т.к. внутри есть заполнение гармоникой, то амплитудный спектр будет симметричен относительно нулевой частоты, при этом он сядет на частоту ±f0 (здесь “спектр сядет” означает, что его центры будут находиться на частоте f0).

2)Т.к. sinc огибает внутреннюю гармонику, то можно воспользоваться обстоятельством, что преобразование Фурье и обратное преобразование Фурье полностью симметричны, то есть имеет смысл такая замена:

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

При этом амплитуды у двух прямоугольников одинаковы и вдвое меньше, чем амплитуда у исходного сигнала, т.к. имеется два представления спектра.

Важно: если у сигнала на середине интервала главного лепестка измерить значение,

снова получится

1

, где

полоса спектров.

 

 

 

Техническая полоса сигнала .

Если над идеальным прямоугольным импульсом провести преобразование Фурье, то получится сигнал типа sinc в диапазоне частот (-∞;+∞). Никакой аппарат с таким диапазоном работать не будет.

Тогда считаем, что если аппарат способен работать на полосе частот (0; 21и )

(половинка положительной половины главного лепестка - область, выделенная красным цветом), то он применим.

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

Согласно формуле (1), по ширине он увеличился в 4 раза, согласно формуле (2) его

амплитуда уменьшилась в 2 раза.

В случае сигнала типа pass (косинусоида) возьмем половинку главного

положительного лепестка (область, выделенная красным цветом).

Тогда после обратного преобразования Фурье сигнал во времени будет таким:

По ширине он расплылся только в 2 раза, потому что мы брали половинку главного лепестка (в случае сигнала типа base мы брали четверть главного лепестка), а на

середине высоты он остался таким же по ширине. Его амплитуда также уменьшилась в

два раза.

11. Дискретное преобразование Фурье ДПФ. Вывод формулы ДПФ.

Дискретное преобразование Фурье ДПФ

Преобразование Фурье в общем случае:

Для периодических функций:

Заметим, что периодическая функция в формуле обозначена как sg(t) – сигнал, хотя, строго говоря, она им не является, так как сигнал – изменение физ. величины, несущее информацию, период. функция информации не несет.

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

Справа – спектр периодической функции (зависимость от частоты). Заметим, что спектр период. функции становится линейчатым.

Размер дискрета (выделен стрелками на графике) равен = 1/пер, если шкала , и равен ω = 2π, если шкала ω.

Преимущественно используем шкалу , но в некоторых случаях полезнее шкала ω (к примеру, при переходе от к ω в прямом (сигнал в спектр) и обратном (спектр в сигнал) преобразованиях Фурье формулы стали симметричны (c отличались коэффициентом)).

Отрицательная частота является лишь математической абстракцией для удобства работы со спектром.

Теперь продискретизируем сигнал:

Левая картинка:

Т – период дискретизации

N – количество периодов дискретизации в периоде сигнала Правая картинка:

пер = 1/ = – частота дискретизации

F – дискрет по частоте

N – количество дискретов в

Как понять, сколько дискретов (N) в fs (промежуток от 0 до fs или от fs до 2fs)? Поделим частоту дискретизации (fs) на размер одного дискрета (F), получим:

Вывод формулы ДПФ

Теперь строим Дискретное Преобразование Фурье:

Здесь:

, n – некоторое количество периодов T

, k – некоторое количество дискретов F

– эф дискретное(дискретное преобразование Фурье)

При дискретном представлении интеграл превращается в сумму

Преобразование границ интегрирования: пер = , поэтому складываем по N и домножаем на T

Преобразуем выражение в степени

(2 = ):

Тогда получим:

Конечная формула:

Пример поворота вектора:

Заметим, что угол поворота в операторе < 0, поэтому в ДПФ крутим по часовой.

Здесь:

ДПФ – дискретное преобразование Фурье, вертим в отрицательную сторону (по часовой)

ОДПФ – обратное дискретное преобразование Фурье, вертим в положительную сторону (против часовой) (так как в формуле ОДПФ степень имеет противоположный знак относительно степени в ДПФ)

12. Пример вычисления ДПФ периодических функций. Формула ОДПФ. ДПФ, ОДПФ комплексных функций.

Пример вычисления ДПФ периодических ф-ций

Входные данные: некоторый сигнал, определенный частотой f1 и начальной фазой 0

1. Продискретизировать сигнал.

Чтобы продискретизировать сигнал, нужна частота дискретизации fs(она обратно пропорциональна периоду дискретизации T как fs = 1/T). Выбираем частоту дискретизации (Здесь f1T и fT - одно и то же!)

**Принцип Найквиста-Котельникова: частота дискретизации должна быть по крайней мере вдвое больше, чем самая большая частота, содержащаяся в сигнале. (В данном примере взяли частоту дискр. в 4 раза больше!) fT-показывает во сколько раз исходная частота больше частоты дискретизации.

Какую константу выбрать для периода? Нужно, чтобы fT при умножении на эту константу было целым числом, потому что тогда в периоде, определенном константой N поместится целое кол-во волн.

В сигнал вместо f1t подставляем f1Tn, где n=0...N-1, далее подставляем остальные значения:

2. Выполнить предсказание спектра.

Надо предсказать спектр, амплитуду(A), фазу(Ф), действительную часть(Re) и мнимую часть(Im) спектра.

На каких отсчетах появится спектр? Первый отсчет появится, если f1 поделить на F. Первый отсчет появится на 1. Другой отсчет будет симметричен ему на -1:

Чтобы избавиться от -1, пользуясь тем, что ф-ция периодическая с периодом 4, к периоду прибавляем -1, получаем 3, и эта -1 переедет на место 3:

Теперь смотрим диапазон не от -k до +kF, а от 0 до +kF и период у нас 4( от 0 до 3, 0’ уже чужой). График симметричен относительно половины периода, половина периода-fн-частота Найквиста.

Какая величина линейки(вертикальная палочка на графике)? Допустим, что U=1, тогда амплитуда(A):

Фаза будет на тех же самых местах(1 и 3) и она задана как -π/4, поэтому на 1

сразу выставляем -π/4, и пользуясь, тем, что фазовый спектр(Ф)

антисимметричен относительно частоты Найквиста.

Соседние файлы в предмете Цифровая обработка сигналов