Function Octbank in Matlab

I found function for one-third-octave filter:

http://www.mathworks.com/matlabcentral/fileexchange/69-octave/content/octave/oct3bank.m

Can someone help me and tell how it will look similar function for octave filter?

Answer 1

[Wayne King]

Do you have fdesign.octave in your installation?

You can design both octave and 1/3 octave filters with fdesign.octave

[Alexander]

Yes, I have. But I have a work function for the 1/3-octave filters (see link above), and I need function for octave filters for on the basis of this function.

[Wayne King]

but fdesign.octave can handle both 1/3 and octave very easily.

[Alexander]

It's easy for you, but not for me :(. Could you remake function Oct3bank for me using fdesign.octave?

Answer 2

[Wayne King]

Alexander, I don't see how using fdesign.octave with an input argument 1 (for octave band) vs. 3 (for 1/3 octave band) is harder than asking somebody to rewrite an entire function.

This is a 1/3 octave band design.

   d1 = fdesign.octave(3,'Class 1','N,F0',6,10000,48000);
   Hd1 = design(d1);
   fvtool(Hd1);

This is an octave band design.

   d2 = fdesign.octave(1,'Class 1','N,F0',6,1000,48000);
   Hd2 = design(d2);
   fvtool(Hd2);

[Alexander]

Wayne, I approximately understood how does work your function for octave filters, but i have one question: there are 17 octave filters in IEC 1260, but if i call following Matlab code

fm = 1000; % center freq

f = fdesign.octave(BandsPerOctave,'Class 1','N,F0',3,fm,65536);

fm = validfrequencies(f);

I get only 10 frequencies in fm variable. If you please can you tell me why it accuses?

[Alexander]

BandsPerOctave is equal 1

[David Santos]

Hi,

I'm using fdesign as Wayne suggest but the problem is that this function only reach to 20KHz and I need to be able to analyze also frequencies close to 50KHz. Anyone knows how to increase the bandwidth of fdesign or any other way of solving the problem?

All the best