DFT Filter Bank Script

14 views (last 30 days)
Stefano
Stefano on 23 Jan 2016
Answered: Stefano on 25 Jan 2016
Hi guys, i've a problem with a script of maximally decimated DFT FBs. The first implementation as parallel filters works, but when i do the script with polyphase and fft/ifft i obtain an aliasing between channel! For generate prototype i used a simple method written by Mitra for pseudo QMF, and i've modulated with complex exponential in the 1st implementation, and used the fft in 2nd script.
In the 1st figure is shown the input signal hilbert(chirp), while in 2nd the output of the synthesis polyphase bank. I attach also the file for who wants try my filters, and i would be grateful if you help me. :(
The script for parallel bank is DFT_FB.m
The script for polyphase bank is DFT_FB_POLY.m, the others files are subfunction.
Thanks
<<
>>
  3 Comments
emehmetcik
emehmetcik on 24 Jan 2016
I get an error saying that the "spectrum_scope" function is missing.
Stefano
Stefano on 24 Jan 2016
True,sorry. The function spectrum scope just contains plot(20*log10(abs(fft(signal,length(signal)))))

Sign in to comment.

Answers (1)

Stefano
Stefano on 25 Jan 2016
This is function spectrum_scope:

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!