Слайд 2
![Немного истории произвольные периодические функции - суммы простейших гармонических функций](/_ipx/f_webp&q_80&fit_contain&s_1440x1080/imagesDir/jpg/23609/slide-1.jpg)
Немного истории
произвольные периодические функции - суммы простейших гармонических функций – синусов
и косинусов кратных частот.
Эти суммы получили название рядов Фурье,
Французский инженер Жан Батист Фурье обосновал метод вычисления коэффициентов тригонометрического ряда, которым можно отображать с абсолютной точностью любую периодическую функцию, определенную на интервале одного периода T = b-a, и удовлетворяющую условиям Дирихле (ограниченная, кусочно-непрерывная, с конечным числом разрывов 1-го рода).
Слайд 3
![Задание функции](/_ipx/f_webp&q_80&fit_contain&s_1440x1080/imagesDir/jpg/23609/slide-2.jpg)
Слайд 4
![w 1 = 2p /T - частота повторения (или частота](/_ipx/f_webp&q_80&fit_contain&s_1440x1080/imagesDir/jpg/23609/slide-3.jpg)
w 1 = 2p /T - частота повторения
(или частота первой
гармоники);
k - номер гармоники.
Этот ряд содержит бесконечное число косинусных или синусных составляющих - гармоник, причем амплитуды этих составляющих ak и bk являются коэффициентами Фурье,
Слайд 5
![Построение графика](/_ipx/f_webp&q_80&fit_contain&s_1440x1080/imagesDir/jpg/23609/slide-4.jpg)
Слайд 6
![Формулы для коэффициентов](/_ipx/f_webp&q_80&fit_contain&s_1440x1080/imagesDir/jpg/23609/slide-5.jpg)
Формулы для коэффициентов
Слайд 7
![Вывод коэффициенты ряда](/_ipx/f_webp&q_80&fit_contain&s_1440x1080/imagesDir/jpg/23609/slide-6.jpg)
Слайд 8
![Термин "spectrum" ("спектр") впервые применил И. Ньютон в 1571 году](/_ipx/f_webp&q_80&fit_contain&s_1440x1080/imagesDir/jpg/23609/slide-7.jpg)
Термин "spectrum" ("спектр") впервые применил И. Ньютон в 1571 году при
описании разложения солнечного света, пропущенного через стеклянную призму, на многоцветную полосу. Он же дал и первую математическую трактовку периодичности волновых движений.
Слайд 9
![Вывод гармоник и функции](/_ipx/f_webp&q_80&fit_contain&s_1440x1080/imagesDir/jpg/23609/slide-8.jpg)
Слайд 10
![Гармонический синтез по 3 гармоникам Сравнение исходной и синтезированной функций](/_ipx/f_webp&q_80&fit_contain&s_1440x1080/imagesDir/jpg/23609/slide-9.jpg)
Гармонический синтез по 3 гармоникам
Сравнение исходной и синтезированной функций
Слайд 11
![Гармонический синтез по 10 гармоникам](/_ipx/f_webp&q_80&fit_contain&s_1440x1080/imagesDir/jpg/23609/slide-10.jpg)
Гармонический синтез по 10 гармоникам
Слайд 12
![Спектральный анализ Спектр амплитуд и спектр фаз](/_ipx/f_webp&q_80&fit_contain&s_1440x1080/imagesDir/jpg/23609/slide-11.jpg)
Спектральный анализ
Спектр амплитуд и спектр фаз
Слайд 13
![Спектр временной зависимости (функции) f(t) называется совокупность ее гармонических составляющих,](/_ipx/f_webp&q_80&fit_contain&s_1440x1080/imagesDir/jpg/23609/slide-12.jpg)
Спектр временной зависимости (функции) f(t) называется совокупность ее гармонических составляющих, образующих
ряд Фурье.
Спектр можно характеризовать некоторой зависимостью
Аk (спектр амплитуд) и
j k (спектр фаз) от частоты w k = kw 1.
Слайд 14
![Термин "spectrum" ("спектр") впервые применил И. Ньютон в 1571 году](/_ipx/f_webp&q_80&fit_contain&s_1440x1080/imagesDir/jpg/23609/slide-13.jpg)
Термин "spectrum" ("спектр") впервые применил И. Ньютон в 1571 году при
описании разложения солнечного света, пропущенного через стеклянную призму, на многоцветную полосу.
Слайд 15
![Спектральный синтез по 3 гармоникам](/_ipx/f_webp&q_80&fit_contain&s_1440x1080/imagesDir/jpg/23609/slide-14.jpg)
Спектральный синтез
по 3 гармоникам
Слайд 16
![](/_ipx/f_webp&q_80&fit_contain&s_1440x1080/imagesDir/jpg/23609/slide-15.jpg)
Слайд 17
![Спектральный анализ с использованием БПФ В Mathcad есть встроенные средства](/_ipx/f_webp&q_80&fit_contain&s_1440x1080/imagesDir/jpg/23609/slide-16.jpg)
Спектральный анализ с использованием БПФ
В Mathcad есть встроенные средства быстрого преобразования
Фурье (БПФ), которые существенно упрощают процедуру приближенного спектрального анализа
Слайд 18
![Встроенные в Mathcad средства быстрого преобразования Фурье (БПФ) fft(v) -](/_ipx/f_webp&q_80&fit_contain&s_1440x1080/imagesDir/jpg/23609/slide-17.jpg)
Встроенные в Mathcad средства быстрого преобразования Фурье (БПФ)
fft(v) - возвращает прямое
БПФ 2m-мерного вещественнозначного вектора v,
где v - вектор, элементы которого хранят отсчеты функции f(t).
Результатом будет вектор А размерности 1 + 2m - 1 с комплексными элементами - отсчетами в частотной области.
Фактически действительная и мнимая части вектора есть коэффициенты Фурье ak и bk,
Слайд 19
![ifft(v) - возвращает обратное БПФ для вектора v с комплексными](/_ipx/f_webp&q_80&fit_contain&s_1440x1080/imagesDir/jpg/23609/slide-18.jpg)
ifft(v) - возвращает обратное БПФ для вектора v с комплексными элементами.
Вектор v имеет 1 + 2m - 1 элементов.
Результатом будет вектор А размерности 2m с действительными элементами.
Слайд 20
![](/_ipx/f_webp&q_80&fit_contain&s_1440x1080/imagesDir/jpg/23609/slide-19.jpg)
Слайд 21
![](/_ipx/f_webp&q_80&fit_contain&s_1440x1080/imagesDir/jpg/23609/slide-20.jpg)
Слайд 22
![Обратное БПФ](/_ipx/f_webp&q_80&fit_contain&s_1440x1080/imagesDir/jpg/23609/slide-21.jpg)
Слайд 23
![Фильтрация аналоговых сигналов](/_ipx/f_webp&q_80&fit_contain&s_1440x1080/imagesDir/jpg/23609/slide-22.jpg)
Фильтрация аналоговых сигналов
Слайд 24
![](/_ipx/f_webp&q_80&fit_contain&s_1440x1080/imagesDir/jpg/23609/slide-23.jpg)
Слайд 25
![Фильтрация - выделение полезного сигнала из его смеси с мешающим](/_ipx/f_webp&q_80&fit_contain&s_1440x1080/imagesDir/jpg/23609/slide-24.jpg)
Фильтрация - выделение полезного сигнала из его смеси с мешающим сигналом
- шумом.
Наиболее распространенный тип фильтрации - частотная фильтрация.
Если известна область частот, занимаемых полезным сигналом, достаточно выделить эту область и подавить те области, которые заняты шумом
Слайд 26
![График полезного сигнала с шумом](/_ipx/f_webp&q_80&fit_contain&s_1440x1080/imagesDir/jpg/23609/slide-25.jpg)
График полезного сигнала с шумом
Слайд 27
![График сигнала после фильтрации](/_ipx/f_webp&q_80&fit_contain&s_1440x1080/imagesDir/jpg/23609/slide-26.jpg)
График сигнала после фильтрации
Слайд 28
![Результат фильтрации Сравнение временных зависимостей исходного и выходного сигналов, показывает,](/_ipx/f_webp&q_80&fit_contain&s_1440x1080/imagesDir/jpg/23609/slide-27.jpg)
Результат фильтрации
Сравнение временных зависимостей исходного и выходного сигналов, показывает, что выходной
сигнал почти полностью повторяет входной
и в значительной мере избавлен от высокочастотных шумовых помех, маскирующих полезный сигнал
Слайд 29
![Задание 1. Вычислить первые шесть пар коэффициентов разложения в ряд](/_ipx/f_webp&q_80&fit_contain&s_1440x1080/imagesDir/jpg/23609/slide-28.jpg)
Задание 1. Вычислить первые шесть пар коэффициентов разложения в ряд Фурье
функции f(t) на отрезке [0, 2p ].
Построить графики 1, 2 и 3 гармоник.
Выполнить гармонический синтез функции f(t) по 1, 2 и 3 гармоникам. Результаты синтеза отобразить графически.
Слайд 30
![Задание 2. Выполнить классический спектральный анализ и синтез функции f(t).](/_ipx/f_webp&q_80&fit_contain&s_1440x1080/imagesDir/jpg/23609/slide-29.jpg)
Задание 2. Выполнить классический спектральный анализ и синтез функции f(t). Отобразить
графически спектры амплитуд и фаз, результат спектрального синтеза функции f(t).
Задание 3. Выполнить численный спектральный анализ и синтез функции f(t). Для этого необходимо задать исходную функцию f(t) дискретно в 32 отсчетах. Отобразить графически спектры амплитуд и фаз, результат спектрального синтеза функции f(t).
Слайд 31
![Задание 4. Выполнить спектральный анализ и синтез функции f(t) с](/_ipx/f_webp&q_80&fit_contain&s_1440x1080/imagesDir/jpg/23609/slide-30.jpg)
Задание 4. Выполнить спектральный анализ и синтез функции f(t) с помощью
БПФ. Для этого необходимо:
задать исходную функцию f(t) дискретно в 128 отсчетах;
выполнить прямое БПФ с помощью функции fft и отобразить графически найденные спектры амплитуд и фаз первых шести гармоник;
выполнить обратное БПФ с помощью функции ifft и отобразить графически результат спектрального синтеза функции f(t).
Слайд 32
![Задание 5. Выполнить фильтрацию функции f(t) с помощью БПФ: синтезировать](/_ipx/f_webp&q_80&fit_contain&s_1440x1080/imagesDir/jpg/23609/slide-31.jpg)
Задание 5. Выполнить фильтрацию функции f(t) с помощью БПФ:
синтезировать функцию
f(t) в виде полезного сигнала, представленного 128 отсчетами вектора v;
к полезному сигналу v присоединить шум с помощью функции rnd (rnd(2) - 1) и сформировать вектор из 128 отсчетов зашумленного сигнала s;
преобразовать сигнал с шумом s из временной области в частотную, используя прямое БПФ (функция fft). В результате получится сигнал f из 64 частотных составляющих;
выполнить фильтрующее преобразование с помощью функции Хевисайда (параметр фильтрации a = 2);
с помощью функции ifft выполнить обратное БПФ и получить вектор выходного сигнала h;
построить графики полезного сигнала v и сигнала, полученного фильтрацией зашумленного сигнала s.