Hi Wessel Lubberhuizen,
Thanks for your response!
In fact, I've carried out some investigations about filter banks.
Now I can say that dct/idct filter bank differs from fft. Sorry for asking about such functionality :)
Actually in case of real-valued filter-banks one should use two different prototype filters and two parrallel phase-shifted dct.
It is told everywhere that real-valued oversampled filter banks has more complex analysis/synthesys than complex-valued.
By the way, I've found the most interesting filter bank implementation in the "Multirate Signal Processing" by Crochiere and Rabiner. It is called Weighted Overlap-Add (WOLA) and it allows to use non-integer oversampling rates.
Hello,
I don't know the reason, but those files don't work at MATLAB 7 (R14).
It looks like .mat files miss variables like N (number of iterations), B+A (plant model), and many others.
Thus the code itself is unusable. It's possible that the book may help if carefully reading it...
Hi, I'd like to thank you for this job, but I have a question about real-valued signals too.
Is it possible to replace fft/ifft pair with dct/idct?
Due to theory the only difference: fft shifts prototype only to the 0 .. pi/2 range, while dct shifts prototypes to 0..pi/2 range, and its copy to -pi/2 .. 0.
In fact I don't quite understand why you multiply x2 by exp(-1..) before the fft. Again theoretically fft contains the same exp(-1) which shifts prototype filter.