Создать 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)