ArbitraryOctaveFilt

 Filters a signal with any arbitrary spectrum smoothed with any fractional octave band average
 
 Syntax:    Y = ARBITRARYOCTAVEFILT(X, SPECT, FREQS, N, FS, OCTBANDWIDTH)
 
 Inputs: 
     x - Input signal to filter as a vector
     SPECT - The spectrum to shape the input signal to
     FREQS - The frequencies of each SPECT element
     N - The length of the filter to usee
     fs - Description
     octBandwidth - Description
 
 Outputs: 
     y - Description

 Example: 
   fs = 16000;
   T = 10;
   N = 1000;
   f = linspace(0,fs/2,N);
   s = 1./f;
   x = wgn(T*fs,1,0);
   y = ArbitraryOctaveFilt(x,s,f,N,fs,1/3);
   pwelch([x y]);
 
 See also: fir2, filter