Разложение пилообразного сигнала в ряд фурье. Представление периодических сигналов рядом фурье

ЛАБОРАТОРНАЯ РАБОТА № 1

РАЗЛОЖЕНИЯ СИГНАЛОВ В РЯД ФУРЬЕ

Цель задания

Ознакомиться с примерами разложения сигналов в ряд Фурье и практически реализовать разложение различного вида сигналов в системе MatLab.

Постановка задачи

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

Для каждого варианта и каждого вида сигнала заданы параметры:

для последовательности прямоугольных импульсов – амплитуда, период повторения и длительность импульсов ;

для меандра, пилообразного сигнала и последовательности треугольных импульсов – амплитуда и период повторения импульсов.

Для всех видов сигналов задано число ненулевых гармоник.

Cоставить программы в системе MatLab и построить графики.

Методические указания

Ряд Фурье

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

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

Синусно-косинусная форма

В этом варианте ряд Фурье имеет следующий вид:

Здесь
– круговая частота, соответствующая периоду повторения сигнала , равному . Входящие в формулу кратные ей частоты
называются гармониками, гармоники нумеруются в соответствии с индексом ; частота
называется –й гармоникой сигнала. Коэффициенты ряда и рассчитываются по формулам:

,

.

Константа рассчитывается по общей формуле для . Само же это слагаемое представляет собой среднее значение сигнала на периоде:

.
Если
является четной функцией , то все будут равны нулю и в формуле ряда Фурье будут присутствовать только косинусные слагаемые. Если является нечетной функцией , равны нулю будут, наоборот, косинусные коэффициенты и в формуле останутся лишь синусные слагаемые.

ПОСЛЕДОВАТЕЛЬНОСТЬ ПРЯМОУГОЛЬНЫХ ИМПУЛЬСОВ

Последовательность прямоугольных импульсов с амплитудой , длительностью и периодом повторения .

Рис. 1 Периодическая последовательность прямоугольных импульсов
Данный сигнал является четной функцией , поэтому для его представления удобнее использовать синусно-косинусную форму ряда Фурье – в ней будут присутствовать только косинусные слагаемые , равные

.

Отношение периода к длительности импульсов называют скважностью последовательности импульсов и обозначают буквой :
.

Представление последовательности прямоугольных импульсов в виде ряда Фурье:

.

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

МЕАНДР

Частным случаем предыдущего сигнала является меандр – последовательность прямоугольных импульсов со скважностью, равной двум, когда длительности импульсов и промежутков между ними становятся равными (рис.2).

Рис. 2 Меандр

При
, получим

Здесь m – произвольное целое число.

При разложении в ряд Фурье четные составляющие будут отсутствовать.

ПИЛООБРАЗНЫЙ СИГНАЛ

В пределах периода он описывается линейной функцией:

Рис. 3. Пилообразный сигнал
Данный сигнал является нечетной функцией, поэтому его ряд Фурье в синусно-косинусной форме будет содержать только синусные слагаемые:

.

Сам ряд Фурье для пилообразного сигнала выглядит следующим образом:

ПОСЛЕДОВАТЕЛЬНОСТЬ ТРЕУГОЛЬНЫХ ИМПУЛЬСОВ

Рис.4. Последовательность треугольных импульсов
Сигнал является четной функцией, поэтому будут присутствовать косинусные составляющие.

Вычислим коэффициенты ряда Фурье:

Сам ряд Фурье имеет следующий вид:

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

Код программы для меандра

N = 8; % число ненулевых гармоник

t = -1:0.01:1; % вектор моментов времени

A = 1; % амплитуда

harmonics = cos(2*pi*nh"*t/T);

Am = 2/pi./nh; % амплитуды гармоник

Am(2:2:end) = -Am(2:2:end); % чередование знаков

s1 = harmonics .* repmat(Am", 1, length(t));

% строки-частичные суммы гармоник

s2 = cumsum(s1);

for k=1:N, subplot(4, 2, k), plot(t, s2(k,:)), end

Р
езультат работы программы

Комментарии : repmat – создание блочной матрицы или многомерного блочного массива из одинаковых блоков. repmat(Am", 1, length(t)) – матрица состоит из 1 блока по вертикали и length(t) блоков по горизонтали, каждый блок является матрицей Am".

Cumsum – расчет частичных сумм элементов.

Subplot (Rows , Cols , N ) – команда для вывода нескольких графиков. Графическое окно разбивается на клетки в виде матрицы, имеющей Rows – строк, Cols – столбцов, и N – клетка становится текущей.

Варианты

№ варианта	Параметры для сигналов
	– амплитуда сигнала	– период повторения сигналов	– длительность сигнала	– число ненулевых гармоник
1	7	3	2	10
2	5	4	3	12
3	4	5	4	14
4	3	6	5	16
5	2	8	6	18
6	5	3	2	14
7	4	4	3	16
8	3	5	4	18
9	2	6	5	10
10	7	8	6	12
11	4	4	3	18
12	3	5	4	10
13	2	6	5	12
14	7	8	6	14
15	5	3	2	16
16	7	3	2	12
17	5	4	3	14
18	4	5	4	16
19	3	6	5	18
20	2	8	6	10
21	5	3	2	16
22	4	4	3	18
23	3	5	4	10
24	2	6	5	12
25	7	8	6	14
26	4	4	3	10
27	3	5	4	12
28	2	6	5	14
29	7	8	6	16
30	5	3	2	18

1.3 Сделать общие выводы.

Часть 2

Цель работы: углубление теоретических знаний, полученных в ходе изучения преобразования Фурье (Fourier Transform)

Необходимые теоретические сведения.

Изменяя период Т и длительность импульса как показано на рис. 7, можно изменять спектр сигнала. С увеличением периода гармоники сближаются, не изменяя форму огибающей.

Рис.7 – Изменение спектра

Смоделируем одиночный прямоугольный импульс, периодическую последовательность импульсов с периодом Т и 10Т .

t = 0:.0314:25; y= square(2*pi*t/10, pi*pi); z = rectpuls(2*pi*t1/10); subplot(4,2,1); plot(t,x) subplot(4,2,2); plot(t,y) subplot(4,2,3); plot(t1,z)

Проведем спектральный анализ полученных сигналов. Непериодические процессы - таковыми являются информационные сигналы , одиночные импульсы , хаотические колебания (шумы ) - обладают сплошным или непрерывным спектром. Интуитивно к такому выводу можно прийти, представляя одиночный импульс частью периодической последовательности, период которой неограниченно увеличивается. Действительно, при увеличении интервала между импульсами гармоники на спектральных диаграммах периодических последовательностей импульсов сближаются: чем реже следуют импульсы, тем меньше расстояние между соседними гармониками (оно равно 1/T ). Спектр одиночного импульса (предельный случай увеличения периода) становится непрерывным, и вводится он не рядами, а интегралами Фурье .

Преобразование Фурье (Fourier transform) является инструментом спектрально­го анализа непериодических сигналов.

В описанных ниже функциях реализован особый метод быстрого преобразования Фурье (БПФ) - Fast Fourier Transform (FFT ), позволяющий резко уменьшить число арифметических операций в ходе приведенных выше преобразований. Метод особенно эффективен, если число обрабатываемых элементов (отсчетов) составляет 2 n , где n - целое положительное число. В MatLab используются следующие функции:

fft(X ) - возвращает для вектора X дискретное преобразование Фурье, по возможности используя алгоритм быстрого преобразования Фурье. Если X - матрица, функция fft возвращает преобразование Фурье для каждого столбца матрицы;

fft(X.n) - возвращает n-точечное преобразование Фурье. Если длина вектора X меньше n, то недостающие элементы заполняются нулями. Если длина X больше п, то лишние элементы удаляются. Когда X - матрица, длина столбцов корректируется аналогично;

ft(X,[ Ldirn) и fft(X,n,dim) - применяют преобразование Фурье к одной из размерностей массива в зависимости от значения параметра dim .

Возможно одномерное обратное преобразование Фурье, реализуемое следующими функциями:

ifft(F) - возвращает результат дискретного обратного преобразования Фурье вектора F . Если F - матрица, то ifft возвращает обратное преобразование Фурье для каждого столбца этой матрицы;

ifft(F.n) - возвращает результат n-точечного дискретного обратного преобразования Фурье вектора F ;

ifft(F.,dim) иу = ifft(X,n,dim) - возвращают результат обратного дискретного преобразования Фурье массива F по строкам или по столбцам в зависимости от значения скаляра dim .

Для любого X результат последовательного выполнения прямого и обратного преобразований Фурье ifft(fft(x)) равен X с точностью до погрешности округления. Если X - массив действительных чисел, ifft(fft(x)) может иметь малые мнимые части.

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

Вызовем программу SPTool (Signal Processing Tool) . Импортируем смоделированные сигналы и рассчитаем спектр сигнала. С этой целью выделяем сигнал в списке сигналов и нажмите кнопку Create , расположенную под списком спектров. В окне Spectrum Viewer в поле Parameters нужно указать метод спектрального анализа. Указываем метод ДПФ (используется быстрое преобразование Фурье БПФ (FFT)). Указав метод, следует щёлкнуть мышью по кнопке Apply . Будет выведен график спектральной плотности мощности. Имеется возможность выводить спектры в линейном или в логарифмическом масштабе (меню Options ).

Непрерывным (сплошным) является спектр хаотических (шумовых ) колебаний . В этом случае спектральная характеристика, как функция частоты, также представляет собой хаотический (случайный ) процесс , статистические параметры которого определяются спецификой конкретного случайного временного процесса. Сформируем сигнал, содержащий регулярные составляющие с частотами 50 Гц и 120 Гц и случайную аддитивную компоненту с нулевым средним.

ЗАДАНИЕ 2

ЛАБОРАТОРНАЯ РАБОТА № 1

РАЗЛОЖЕНИЯ СИГНАЛОВ В РЯД ФУРЬЕ

Цель задания

Ознакомиться с примерами разложения сигналов в ряд Фурье и практически реализовать разложение различного вида сигналов в системе MatLab.

Постановка задачи

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

Для каждого варианта и каждого вида сигнала заданы параметры:

для последовательности прямоугольных импульсов – амплитуда, период повторения и длительность импульсов;

для меандра, пилообразного сигнала и последовательности треугольных импульсов – амплитуда и период повторения импульсов.

Для всех видов сигналов задано число ненулевых гармоник.

Cоставить программы в системеMatLabи построить графики.

Постановка задачи.

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

Результаты выполнения программ – графики промежуточных стадий суммирования.

Методические указания

Ряд Фурье

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

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

Синусно-косинусная форма

В этом варианте ряд Фурье имеет следующий вид:

Здесь
– круговая частота, соответствующая периоду повторения сигнала, равному. Входящие в формулу кратные ей частоты
называются гармониками, гармоники нумеруются в соответствии с индексом ; частота
называется –й гармоникой сигнала. Коэффициенты ряда ирассчитываются по формулам:

,

.

Константа рассчитывается по общей формуле для. Само же это слагаемое представляет собой среднее значение сигнала на периоде:

.

Если
является четной функцией, то всебудут равны нулю и в формуле ряда Фурье будут присутствовать только косинусные слагаемые. Если
является нечетной функцией, равны нулю будут, наоборот, косинусные коэффициентыи в формуле останутся лишь синусные слагаемые.

ПОСЛЕДОВАТЕЛЬНОСТЬ ПРЯМОУГОЛЬНЫХ ИМПУЛЬСОВ

Последовательность прямоугольных импульсов с амплитудой , длительностьюи периодом повторения.

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

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

.

Отношение периода к длительности импульсов называют скважностью последовательности импульсов и обозначают буквой :
.

Представление последовательности прямоугольных импульсов в виде ряда Фурье:

.

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

МЕАНДР

Частным случаем предыдущего сигнала является меандр – последовательность прямоугольных импульсов со скважностью, равной двум, когда длительности импульсов и промежутков между ними становятся равными (рис.2).

Рис. 2 Меандр

При
, получим

Здесь m – произвольное целое число.

При разложении в ряд Фурье четные составляющие будут отсутствовать.

ПИЛООБРАЗНЫЙ СИГНАЛ

В пределах периода он описывается линейной функцией:

Рис. 3. Пилообразный сигнал

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

.

Сам ряд Фурье для пилообразного сигнала выглядит следующим образом:

ПОСЛЕДОВАТЕЛЬНОСТЬ ТРЕУГОЛЬНЫХ ИМПУЛЬСОВ

Рис.4. Последовательность треугольных импульсов

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

Вычислим коэффициенты ряда Фурье:

Сам ряд Фурье имеет следующий вид:

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

Код программы для меандра

N= 8; % число ненулевых гармоник

t= -1:0.01:1; % вектор моментов времени

A= 1; % амплитуда

T= 1; % период

nh= (1:N)*2-1; % номера ненулевых гармоник

harmonics = cos(2*pi*nh"*t/T);

Am= 2/pi./nh; % амплитуды гармоник

Am(2:2:end) = -Am(2:2:end); % чередование знаков

s1 = harmonics .* repmat(Am", 1, length(t));

% строки-частичные суммы гармоник

for k=1:N, subplot(4, 2, k), plot(t, s2(k,:)), end

Р
езультат работы программы

Комментарии :repmat – создание блочной матрицы или многомерного блочного массива из одинаковых блоков.repmat(Am", 1,length(t)) – матрица состоит из 1 блока по вертикали иlength(t) блоков по горизонтали, каждый блок является матрицейAm".

Cumsum – расчет частичных сумм элементов.

Subplot (Rows , Cols , N ) – команда для вывода нескольких графиков. Графическое окно разбивается на клетки в виде матрицы, имеющейRows – строк,Cols – столбцов, иN – клетка становится текущей.

Варианты

№ варианта	Параметры для сигналов
	–амплитуда сигнала	–период повторения сигналов	–длительность сигнала	–число ненулевых гармоник

Общие описания

Французский математик Фурье (Ж. Б. Ж. Фурье 1768-1830) провоз гласил достаточно смелую для своего времени гипотезу. Согласно этой гипотезе не существует функции, которую нельзя было бы разложить в тригонометрический ряд. Однако, к сожалению, в то время такая идея не была воспринята всерьез. И это естественно. Сам Фурье не смог привести убедительных доказательств, а интуитивно поверить в гипотезу Фурье очень трудно. Особенно нелегко представить тот факт, что при сложении простых функций, подобных тригонометрическим, воспроизводятся функции, совершенно на них не похожие. Но если предположить, что гипотеза Фурье верна, то периодический сигнал любой формы можно разложить на синусоиды различных частот, или наоборот, посредством соответствующего сложения синусоид с разными частотами возможно синтезировать сигнал какой угодно формы. Следовательно, если эта теория верна, то ее роль в обработке сигналов может быть очень велика. В этой главе первым делом попы­таемся проиллюстрировать правильность гипотезы Фурье.

Рассмотрим функцию

f(t)= 2sin t – sin 2t

Простой тригонометрический ряд

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

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

Коэффициенты этого тригонометрического ряда обозначим как b k , где k - целые числа. Если внимательно посмотреть на последнее соотношение, то видно, что коэффициенты можно описать следующим выражением:

Тогда функцию f(t) можно представить следующим образом:

Коэффициенты b k - это амплитуды синусоид с угловой частотой к. Иначе говоря, они задают величину частотных составляющих.

Рассмотрев случай, когда верхний индекс к равен 10, т.е. М= 10. Увеличив значение М до 100, получим функцию f(t).

Эта функция, будучи тригонометрическим рядом, по форме приближается к пилообразному сигналу. И, похоже, гипотеза Фурье совершенно верна по отноше­нию к физическим сигналам, с которыми мы имеем дело. К тому же в этом примере форма сигнала не гладкая, а включает точки разрыва. И то, что функция воспроизводится даже в точках разрыва, выглядит многообещающим.

В физическом мире действительно много явлений, которые можно представить как суммы колебаний различных частот. Типичным примером этих явлений является свет. Он представляет собой сумму электромагнитных волн с длиной волны от 8000 до 4000 ангстрем (от красного цвета свечения до фиолетового). Вы, конечно, знаете, что если белый свет пропустить через призму, то появится спектр из семи чистых цветов. Это происходит потому, что коэффициент преломления стекла, из которого сделана призма, изменяется в зависимости от длины электромагнитной волны. Это как раз и является доказательством того, что белый свет - это сумма световых волн различной дли­ны. Итак, пропустив свет через призму и получив его спектр, мы можем проанализировать свойства света, исследуя цветовые комбинации. Подобно этому, посредством разложения принятого сигнала на различные частотные составляющие, мы можем узнать, как возник первоначальный сигнал, по какому пути он следовал или, наконец, какому внешнему влиянию он подвергался. Одним словом, мы можем получить информацию для выяснения происхождения сигнала.

Подобный метод анализа называется спектральным анализом или анализом Фурье.

Рассмотрим следующую систему ортонормированных функций:

Функцию f(t) можно разложить по этой системе функций на отрезке [-π, π] следующим образом:

Коэффициенты α k , β k , как было показано ранее, можно выразить через скалярные произведения:

В общем виде функцию f(t) можно представить следующим образом:

Коэффициенты α 0 , α k , β k называют коэффициентами Фурье, а подобное представление функции называется разложением в ряд Фурье. Иногда такое представление называют действительным разложением в ряд Фурье, а коэффициенты - действительными коэффициентами Фурье. Термин «действительный» вводится для того, чтобы отличить представленное разложение от разложения в ряд Фурье в комплексной форме.

Как уже было сказано ранее, произвольную функцию можно разложить по системе ортогональных функций, даже если функции из этой системы не представляются в виде тригонометрического ряда. Обычно под разложением в ряд Фурье подразумевается разложение в тригонометрический ряд. Если коэффициенты Фурье выразить через α 0 , α k , β k получим:

Поскольку при k = 0 coskt = 1, то константа а 0 /2 выражает общий вид коэффициента а k при k = 0.

В соотношении (5.1) колебание самого большого периода, представленное суммой cos t и sin t, называют колебанием основной частоты или первой гармоникой. Колебание с периодом, равным половине основного периода, называют второй гармоникой. Колебание с периодом, равным 1/3 основного периода, называют третьей гармоникой и т.д. Как видно из соотношения (5.1) a 0 является постоянной величиной, выражающей среднее значение функции f{t) . Если функция f(t) представляет собой электрический сигнал, то а 0 представляет его постоянную составляющую. Следовательно, все остальные коэффициенты Фурье выражают его переменные составляющие.

На Рис. 5.2 представлен сигнал и его разложение в ряд Фурье: на постоянную составляющую и гармоники различных частот. Во временной области, где переменной величиной является время, сигнал выражается функцией f(t), а в частотной области, где переменной величиной является частота, сигнал представляется коэффициен­тами Фурье (a k , b к).

Первая гармоника является периодической функцией с периодом 2 π.Прочие гармоники также имеют период, кратный 2 π. Исходя из этого, при формировании сигнала из составляющих ряда Фу­рье мы, естественно, получим периодическую функцию с периодом 2 π. А если это так, то разложение в ряд Фурье - это, собственно говоря, способ представления периодических функций.

Разложим в ряд Фурье сигнал часто встречающегося вида. Например, рассмотрим упомянутую ранее пилообразную кривую (Рис. 5.3). Сигнал такой формы на отрезке - π < t < π я выражается функцией f(t) = t , поэтому коэффициенты Фурье могут быть выражены следующим образом:

Пример 1.

Разложение в ряд Фурье сигнала пилообразной формы

f(t) = t,

Примеры разложения в ряд Фурье.

а) Последовательность прямоугольных импульсов .

Рис 2. Последовательность прямоугольных импульсов.

Данный сигнал является четной функцией и для его представления удобно использовать синусно-косинусную форму ряда Фурье:

. (17)

Длительность импульсов и период их следования входят в полученную формулу в виде отношения, ĸᴏᴛᴏᴩᴏᴇ принято называть скважностью последовательности импульсов :.

. (18)

Значение постоянного слагаемого ряда с учетом соответствует:

.

Представление последовательности прямоугольных импульсов в виде ряда Фурье имеет вид:

. (19)

График функции носит лепестковый характер.
Размещено на реф.рф
Горизонтальную ось градуируют в номерах гармоник и в частотах.

Рис 3. Представление последовательности прямоугольных импульсов

в виде ряда Фурье.

Ширина лепестков , измеренная в количестве гармоник, равна скважности (при , имеем , в случае если ). Отсюда следует важное свойство спектра последовательности прямоугольных импульсов – в нем отсутствуют гармоники с номерами, кратными скважности . Расстояние по частоте между сосœедними гармониками равно частоте следования импульсов . Ширина лепестков, измеренная в единицах частоты, равна , ᴛ.ᴇ. обратно пропорциональна длительности сигнала. Можно сделать вывод: чем короче импульс, тем шире спектр .

б) Пилообразный сигнал.

Рис 4. Пилообразный сигнал.

Пилообразный сигнал в пределах периода описывается линœейной функцией

, . (20)

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

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

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

в) Последовательность треугольных импульсов .

Ряд Фурье имеет вид:

Рис 5. Последовательность треугольных импульсов.

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

Лекция №3. Преобразование Фурье.

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

Примеры разложения в ряд Фурье. - понятие и виды. Классификация и особенности категории "Примеры разложения в ряд Фурье." 2017, 2018.