Can anybody help me with the following code?

Question

Darsana P M 2017-3-16

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/330239-can-anybody-help-me-with-the-following-code

评论： Darsana P M 2017-3-17

[z,fs]=audioread('bird1.mp3'); 
bfil=fft(z); %fft of input signal
wn=[4000 8000]/(fs/2);   %bandpass
[b,a]=butter(6,wn);
fvtool(b,a);

I was able to draw the graph for this particular range. Can i draw wn=[8000 12000]/(fs/2); n the same graph itself, instead of another figure???

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

GEEVARGHESE TITUS 2017-3-16

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/330239-can-anybody-help-me-with-the-following-code#answer_259042

在 MATLAB Online 中打开

I not sure of obtaining multiple plots in the same window using fvtool. Instead you can get what you required, using functions freqz(), plot() and hold on.

[z,fs]=audioread('bird1.mp3'); 
bfil=fft(z); %fft of input signal
wn=[4000 8000]/(fs/2);   %bandpass
[b,a]=butter(6,wn);
[h,w]=freqz(b,a);
plot(w/pi,20*log10(abs(h)),'r');
hold on;
wn=[8000 12000]/(fs/2);
[b,a]=butter(6,wn);
[h,w]=freqz(b,a);
plot(w/pi,20*log10(abs(h)),'b');

3 个评论
显示 1更早的评论隐藏 1更早的评论

GEEVARGHESE TITUS 2017-3-17

Sorry , i was not able to understand your question? Elaborate please. Are you thinking of bank of butterworth filers?

Darsana P M 2017-3-17

yaa exactlty. How can i write a code for a bank of butterworth filters? Such that i can divide the given audio signal into different bands ie 10-20Khz, 20-30Khz and so on.

请先登录，再进行评论。