1. 
   Цель работы
   

Разработать алгоритм и программу, реализующую спектральное разделение тестовой фонограммы на две области: НЧ (F_от и гармоники F_от) и ВЧ (формантную).

2. 
   Описание этапов выполнения работы
   

Для выполнения лабораторной работы используем фонограмму под названием Music.aac. Частота дискретизации исходного звукового сигнала Исходный звуковой сигнал во временной области представлен на рисунке 2.1.Частота основного тона таким образом частота среза для ФНЧ будет равняться 440 Гц. А так как форманта – максимум энергии в спектре звука, то, в соответствии со спектром исходного звукового сигнала, представленного на рисунке 2.2, частота среза для ФВЧ будет равняться 740 Гц. Рисунок 2.1 – Исходный звуковой сигнал во временной областиС помощью алгоритма, описанного в функции FFT разложим спектр исходного звукового сигнала. Спектр представлен на рисунке 2.2.Рисунок 2.2 – Спектр исходного звукового сигналаДалее, с помощью фильтра Баттерворта, спроектируем высокочастотный фильтр в среде Matlab. АЧХ и ФЧХ фильтра представлены на рисунке 2.3.Рисунок 2.3 – АЧХ и ФЧХ ВЧ фильтра Баттерворта Спектр звукового сигнала на выходе фильтра представлен на рисунке 2.4. Масштаб графика взят в соответствии со спектром исходного звукового сигнала, представленного на рисунке 2.1.

---

Рисунок 2.4 – Спектр звукового сигнала на выходе ВЧ фильтра БаттервортаТаким же образом, с помощью фильтра Баттерворта, спроектируем НЧ фильтр в среде Matlab. АЧХ и ФЧХ фильтра представлены на рисунке 2.5.Рисунок 2.5 – АЧХ и ФЧХ НЧ фильтра Баттерворта Спектр звукового сигнала на выходе фильтра представлен на рисунке 2.6.Рисунок 2.6 – Спектр звукового сигнала на выходе НЧ фильтра БаттервортаЛистинг программы представлен ниже:clc, close all, clear allfilename='Music.aac';[Audio, Fs]=audioread(filename);Time=(1:length(Audio))*1/Fs; %Зная частоту дискретизации и количество точек, рассчитываем вектор времениplot(Time, Audio'), grid on %Строим временной график сигналаtitle(['Исходный звуковой сигнал (Fs=48000 кГц)'])xlabel('Время, с')ylabel ('Амплитуда, В')[f1,Shift1]=FFT(Fs,Audio);xlim([0 Fs/5])title(['Спектр исходного звукового сигнала'])hold on%%% ВЧ фильтрацияfigurefreq_sreza=740;[b,a]=butter(5, freq_sreza/(Fs/2),'high');freqz(b,a);ylim([-50 0])Audio1=filter(b,a,Audio);[f2,Shift2]=FFT(Fs,Audio1);xlim([0 Fs/2]);ylim([0 1600]);title(['Спектр звукового сигнала на выходе ВЧ фильтра'])%%% НЧ фильтрацияfigurefreq_sreza=440;[b,a]=butter(5, freq_sreza/(Fs/2));freqz(b,a);ylim([-50 0])Audio2=filter(b,a,Audio);[f2,Shift2]=FFT(Fs,Audio2);xlim([0 Fs/2]);ylim([0 1600]);title(['Спектр звукового сигнала на выходе НЧ фильтра'])sound(Audio2,Fs);ВыводыАЧХ фильтра Баттерворта максимально гладкая на частотах полосы пропускания и снижается практически до нуля на частотах полосы подавления. При отображении частотного отклика фильтра Баттерворта на логарифмической АФЧХ, амплитуда снижается к минус бесконечности на частотах полосы подавления.В сравнении с фильтрами Чебышёва I и II типов или эллиптическим фильтром, фильтр Баттерворта имеет более пологий спад характеристики и поэтому должен иметь больший порядок (что более трудно в реализации) для того, чтобы обеспечить нужные характеристики на частотах полосы подавления. Однако фильтр Баттерворта имеет более линейную фазо-частотную характеристику на частотах полосы пропускания.

---

Смотрите также файлы

• 1Общий раздел.docx

• Занятие по гпп тема Участие прокурора в гражданском процессе Цели и задачи участия прокуратуры в гражданском процессе. Реформа роли прокурора в гражданском процессе причины и последствия (ретроспективный анализ законодательства!).docx

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

• Значение идеологического многообразия и политического плюрализма для отечественного государства.pdf

• Краткое содержание учебнопроизводственной деятельности Подпись Руководителя практики.docx