Скачиваний:
26
Добавлен:
01.05.2014
Размер:
291.84 Кб
Скачать

Ход выполнения работы

Примечание: для выполнения данной лабораторной работы (для формирования сигналов, проведения над ними математических операций, а также построения их графиков) будет использоваться математический пакетMathSoft ® MATLAB 7.0.1.

Наши ожидания

Наши ожидания перед началом какого-либо исследования – это почти всегда идеализация, ибо мы, естественно, предполагаем наиболее благоприятный для нас исход проводимых экспериментов. Этап построения подобных «идеализированных моделей» очень важен, ибо позволяет четко уяснить суть решаемой проблемы, а также установить желаемый наилучший результат ее решение (или же наилучший способ ее решения) для последующего сравнения с ним данных реальных экспериментов и оценки пригодности полученных практических результатов в сравнении с определенным нами эталоном.

Итак, что же мы хотели бы получить?

Было бы очень здорово, чтобы дискретное преобразование Фурье (далее - ДПФ) работало «максимально хорошо», т.е. не давало бы в процессе преобразования сигнала различных, скажем так, нежелательных явлений, ухудшающих конечный результат, а значит, усложняющих дальнейший анализ.

Если говорить более конкретно, в идеальном случае мы представляем ДПФ как параллельную систему фильтров следующего характера:

  • У каждого из этих фильтров ширина амплитудно-частотной характеристики (АЧХ) равна 1 Гц.

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

Дпф – практические исследования

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

Сначала зададим его с целой частотой () и выполним для этого сигнала ДПФ (Рис. 1):

Рис. 1.ДПФ для косинуса с целым значением частоты, равным 5.

Как видно из представленного выше рисунка, никаких посторонних элементов на графике спектра не наблюдается. Наша «идеализированная модель» в данном случае хорошо работает.

Продолжим экспериментировать. Построим теперь спектр той же функции, но на этот раз зададим для нее не целое значение частоты (Рис. 2):

Рис. 2.ДПФ для косинуса с дробным значением частоты, равным 5,5.

Вот тут-то и начинаются для нас «неприятности». Хотя по прежнему значение частоты 5 на спектре выделяется лучше всех остальных, на значение частоты 5,5 (судя по полученному спектру) отреагировали сразу несколько фильтров нашей «идеализированной модели». Идеального результата, увы, не получилось – ДПФ проявило свои нежелательные для нас дефекты.

Создадим теперь семейство гармонических сигналов (пусть это будут все те же функции косинуса). Для каждого сигнала из семейства выполним ДПФ. Наша цель здесь – построить зависимость амплитуды некоторого столбца в спектре от частоты подаваемого сигнала из созданного семейства гармонических сигналов. Для определенности будем подавать сигналы диапазона 10-20 Гц с шагом 0,1 Гц. А для дальнейшего исследования выберем 15-й и 16-й столбцы спектра. Графики указанных зависимостей приведены на Рис. 3(для 15-го столбца – график выделен сплошной линией, для 16-го столбца – пунктирной линией).

Рис. 3.Зависимость амплитуды 15-го и 16-го столбцов спектра

от частоты подаваемого сигнала.

Итак, что же можно сказать?

На представленном выше рисунке очень четко видны малые «лепестки» - это реакция 15-го фильтра из «идеализированной модели» на частоты из тех диапазонов, за которые он не отвечает и которые он (в идеальном случае) пропускать не должен.

Также наблюдается перекрытие АЧХ для 15-го и 16-го столбцов (соответственно, для 15-го и 16-го фильтров из «идеализированной модели»). Т.е. при подаче частоты 14,5 Гц на нее весьма «активно» реагирует не только 15-й фильтр, но и 16-й тоже (при этом в область действия последнего указанная частота не входит).

Описанные явления и есть дефекты ДПФ.

Соседние файлы в папке Лабораторная работа 31
  • #
    01.05.2014291.84 Кб26MOTS_Lab#3.doc
  • #
    01.05.20143.75 Кб20Чистовой вариант кода.m