Создайте gammatoneFilterBank объект. Функции freqz чтобы получить комплексную частотную характеристику, H, банка фильтров и вектор частот, f, при котором вычисляется ответ. Постройте график частотной характеристики группы фильтров.

gammaFiltBank = gammatoneFilterBank;
[H,f] = freqz(gammaFiltBank);

plot(f,abs(H))
xlabel('Frequency (Hz)')

Чтобы получить частотную характеристику подмножества фильтров в группе фильтров, задайте второй аргумент как вектор-строка индексов между одним и количеством фильтров в группе фильтров. Получите частотную характеристику 10-го фильтра в наборе фильтров и постройте график частотной характеристики величины.

[H,f] = freqz(gammaFiltBank,10);

plot(f,abs(H))
xlabel('Frequency (Hz)')

Чтобы задать количество точек БПФ, используемых для вычисления частотной характеристики, используйте NFFT Пара "имя-значение". Задайте, что частотная характеристика вычисляется с помощью БПФ с 128 точками. Постройте график частотной характеристики величины.

[H,f] = freqz(gammaFiltBank,'NFFT',128);

plot(f,abs(H))
xlabel('Frequency (Hz)')

Чтобы визуализировать только частотную характеристику величины, вызовите freqz без выходных аргументов. Постройте график величины частотной характеристики в дБ фильтров 20, 21 и 22 с использованием ДПФ с 1024 точками.

freqz(gammaFiltBank,[20,21,22],'NFFT',1024)

Создайте octaveFilterBank объект. Функции freqz чтобы получить комплексную частотную характеристику, H, банка фильтров и вектор частот, f, при котором вычисляется ответ. Постройте график величины частотной характеристики в дБ.

octFiltBank = octaveFilterBank;
[H,f] = freqz(octFiltBank);

plot(f,20*log10(abs(H)))
xlabel('Frequency (Hz)')
ylabel('Magnitude (dB)')
set(gca,'XScale','log')
axis([10 octFiltBank.SampleRate/2 -100 2])

Чтобы получить частотную характеристику подмножества фильтров в группе фильтров, задайте второй аргумент как вектор-строка индексов между одним и количеством фильтров в группе фильтров. Получите частотную характеристику 5-го фильтра в наборе фильтров и постройте график частотной характеристики величины в дБ.

[H,f] = freqz(octFiltBank,5);

plot(f,20*log10(abs(H)))
xlabel('Frequency (Hz)')
ylabel('Magnitude (dB)')
set(gca,'XScale','log')
axis([10 octFiltBank.SampleRate/2 -100 2])

Чтобы задать количество точек БПФ, используемых для вычисления частотной характеристики, используйте NFFT Пара "имя-значение". Задайте, что частотная характеристика вычисляется с помощью БПФ с 8192 точками. Постройте график величины частотной характеристики в дБ.

[H,f] = freqz(octFiltBank,'NFFT',8192);

plot(f,20*log10(abs(H)))
xlabel('Frequency (Hz)')
ylabel('Magnitude (dB)')
set(gca,'XScale','log')
axis([10 octFiltBank.SampleRate/2 -100 2])

Чтобы визуализировать только частотную характеристику величины, вызовите freqz без выходных аргументов. Постройте график частотной характеристики величины в дБ фильтров 4, 5 и 6 с помощью ДПФ с 1024 точками.

freqz(octFiltBank,[4,5,6],'NFFT',1024)