Nsound  0.9.4
Public Member Functions | Protected Types | Protected Attributes | List of all members
Nsound::FilterSlinky Class Reference

Base class for IIR Filters, defines the interface. More...

#include <Nsound/FilterSlinky.h>

Inheritance diagram for Nsound::FilterSlinky:
Inheritance graph
[legend]

Public Member Functions

 FilterSlinky (const float64 &sample_rate, const float64 &delay_time, const float64 &frequency_window)
 
virtual ~FilterSlinky ()
 
AudioStream filter (const AudioStream &x)
 
AudioStream filter (const AudioStream &x, const Buffer &frequencies)
 
Buffer filter (const Buffer &x)
 
Buffer filter (const Buffer &x, const Buffer &frequencies)
 
float64 filter (const float64 &x)
 
float64 filter (const float64 &x, const float64 &frequency)
 
void reset ()
 
void setRealtime (bool flag)
 
AudioStream filter (const AudioStream &x, const float64 &frequency)
 
Buffer filter (const Buffer &x, const float64 &frequency)
 
Buffer getFrequencyAxis (const uint32 n_fft=8192)
 
Buffer getFrequencyResponse (const uint32 n_fft=8192)
 
Buffer getImpulseResponse (const uint32 n_samples=8192)
 
virtual uint32 getKernelSize () const
 
Buffer getPhaseResponse ()
 
float64 getSampleRate () const
 
void plot (boolean show_phase=false)
 

Protected Types

typedef std::vector
< FilterDelay * > 
DelayVector
 
typedef std::vector
< FilterBandPassIIR * > 
FilterVector
 

Protected Attributes

DelayVector delay_lines_
 
FilterVector filters_
 
float64 sample_rate_
 
float64 two_pi_over_sample_rate_
 
float64 sample_time_
 
uint32 kernel_size_
 
bool is_realtime_
 

Detailed Description

Base class for IIR Filters, defines the interface.

Definition at line 53 of file FilterSlinky.h.

Member Typedef Documentation

typedef std::vector< FilterDelay * > Nsound::FilterSlinky::DelayVector
protected

Definition at line 87 of file FilterSlinky.h.

typedef std::vector< FilterBandPassIIR * > Nsound::FilterSlinky::FilterVector
protected

Definition at line 88 of file FilterSlinky.h.

Constructor & Destructor Documentation

FilterSlinky::FilterSlinky ( const float64 sample_rate,
const float64 delay_time,
const float64 frequency_window 
)

Definition at line 45 of file FilterSlinky.cc.

References delay_lines_, filters_, and Nsound::Filter::sample_rate_.

49  :
50  Filter(sample_rate),
51  delay_lines_(),
52  filters_()
53 {
54  float64 f = 0.0;
55 
56  while( f < sample_rate_ / 2.0)
57  {
58  if(f == 0.0)
59  {
60  f = frequency_octave;
61  }
62  else
63  {
64  f += frequency_octave;
65  }
66  }
67 
68  // Create the highest frequency pass band filter.
69 
70  float64 f_high = f;
71  f -= frequency_octave;
72  float64 f_low = f;
73 
74  filters_.push_back(
77  6,
78  f_low,
79  f_high,
80  0.01));
81 
82  // Create pass band filters for the remaing frequency octaves.
83  float64 dt = delay_octave;
84  while(f >= frequency_octave)
85  {
86  f_high = f;
87  f -= frequency_octave;
88  f_low = f;
89 
90  if(f_high == frequency_octave)
91  {
92  f_low = 0.0;
93  }
94 
95  filters_.push_back(
98  6,
99  f_low,
100  f_high,
101  0.01));
102 
103  delay_lines_.push_back(
104  new FilterDelay(sample_rate_, dt));
105 
106  dt += delay_octave;
107  }
108 
109 //~ cout << "n_filters = " << filters_.size() << endl
110 //~ << "n_delays = " << delay_lines_.size() << endl;
111 }
FilterVector filters_
Definition: FilterSlinky.h:91
DelayVector delay_lines_
Definition: FilterSlinky.h:90
double float64
Definition: Nsound.h:146
Filter(const float64 &sample_rate)
Definition: Filter.cc:41
A class for filtering audio in the frequecy domain.
Definition: FilterDelay.h:47
float64 sample_rate_
Definition: Filter.h:113
FilterSlinky::~FilterSlinky ( )
virtual

Definition at line 115 of file FilterSlinky.cc.

References delay_lines_, filter(), and filters_.

116 {
117  FilterVector::iterator filter = filters_.begin();
118  FilterVector::iterator end = filters_.end();
119 
120  DelayVector::iterator delay = delay_lines_.begin();
121 
122  delete *filter;
123  ++filter;
124 
125  while(filter != end)
126  {
127  delete *filter;
128  delete *delay;
129 
130  ++filter;
131  ++delay;
132  }
133 }
FilterVector filters_
Definition: FilterSlinky.h:91
DelayVector delay_lines_
Definition: FilterSlinky.h:90
AudioStream filter(const AudioStream &x)

Member Function Documentation

AudioStream FilterSlinky::filter ( const AudioStream x)

Definition at line 138 of file FilterSlinky.cc.

References Nsound::AudioStream::getNChannels(), and Nsound::AudioStream::getSampleRate().

Referenced by filter(), reset(), and ~FilterSlinky().

139 {
140  uint32 n_channels = x.getNChannels();
141 
142  AudioStream y(x.getSampleRate(), n_channels);
143 
144  for(uint32 channel = 0; channel < n_channels; ++channel)
145  {
146  y[channel] = FilterSlinky::filter(x[channel]);
147  }
148 
149  return y;
150 }
unsigned int uint32
Definition: Nsound.h:153
float64 getSampleRate() const
Returns the sample rate of the stream.
Definition: AudioStream.h:217
uint32 getNChannels(void) const
Returns the number of audio channels in the stream.
Definition: AudioStream.h:212
AudioStream filter(const AudioStream &x)
AudioStream FilterSlinky::filter ( const AudioStream x,
const Buffer frequencies 
)

Definition at line 154 of file FilterSlinky.cc.

References filter().

155 {
156  return FilterSlinky::filter(x);
157 }
AudioStream filter(const AudioStream &x)
Buffer FilterSlinky::filter ( const Buffer x)

Definition at line 161 of file FilterSlinky.cc.

References filter(), and Nsound::Buffer::getLength().

162 {
163  Buffer y;
164 
165  uint32 n_samples = x.getLength();
166 
167  for(uint32 n = 0; n < n_samples; ++n)
168  {
169  y << FilterSlinky::filter(x[n]);
170  }
171 
172  return y;
173 }
unsigned int uint32
Definition: Nsound.h:153
uint32 getLength() const
Returns the number of samples in the Buffer.
Definition: Buffer.h:587
A Buffer for storing audio samples.
Definition: Buffer.h:60
AudioStream filter(const AudioStream &x)
Buffer FilterSlinky::filter ( const Buffer x,
const Buffer frequencies 
)

Definition at line 177 of file FilterSlinky.cc.

References filter().

178 {
179  return FilterSlinky::filter(x);
180 }
AudioStream filter(const AudioStream &x)
float64 FilterSlinky::filter ( const float64 x)
virtual

Implements Nsound::Filter.

Definition at line 184 of file FilterSlinky.cc.

References delay_lines_, filter(), and filters_.

185 {
186  float64 y = 0.0;
187 
188  FilterVector::iterator filter = filters_.begin();
189  FilterVector::iterator end = filters_.end();
190 
191  DelayVector::iterator delay = delay_lines_.begin();
192 
193  // There is always one more filter than delay lines, this is because the
194  // highest frequency have zero delay.
195 
196  y += (*filter)->filter(x);
197 
198  ++filter;
199 
200  while(filter != end)
201  {
202  y += (*delay)->filter((*filter)->filter(x));
203 
204  ++filter;
205  ++delay;
206  }
207 
208  return y;
209 }
FilterVector filters_
Definition: FilterSlinky.h:91
DelayVector delay_lines_
Definition: FilterSlinky.h:90
double float64
Definition: Nsound.h:146
AudioStream filter(const AudioStream &x)
float64 FilterSlinky::filter ( const float64 x,
const float64 frequency 
)
virtual

Implements Nsound::Filter.

Definition at line 213 of file FilterSlinky.cc.

References filter().

214 {
215  return FilterSlinky::filter(x);
216 }
AudioStream filter(const AudioStream &x)
void FilterSlinky::reset ( )
virtual

Implements Nsound::Filter.

Definition at line 220 of file FilterSlinky.cc.

References delay_lines_, filter(), and filters_.

221 {
222  FilterVector::iterator filter = filters_.begin();
223  FilterVector::iterator end = filters_.end();
224 
225  DelayVector::iterator delay = delay_lines_.begin();
226 
227  (*filter)->reset();
228  ++filter;
229 
230  while(filter != end)
231  {
232  (*filter)->reset();
233  (*delay)->reset();
234 
235  ++filter;
236  ++delay;
237  }
238 
239 }
FilterVector filters_
Definition: FilterSlinky.h:91
DelayVector delay_lines_
Definition: FilterSlinky.h:90
AudioStream filter(const AudioStream &x)
void Nsound::Filter::setRealtime ( bool  flag)
inlineinherited

Definition at line 57 of file Filter.h.

References Nsound::Filter::is_realtime_.

57 {is_realtime_ = flag;}
bool is_realtime_
Definition: Filter.h:118
AudioStream Filter::filter ( const AudioStream x,
const float64 frequency 
)
inherited

Definition at line 76 of file Filter.cc.

References Nsound::Filter::filter(), Nsound::AudioStream::getNChannels(), Nsound::AudioStream::getSampleRate(), Nsound::Filter::is_realtime_, M_THROW, and Nsound::Filter::reset().

77 {
78  if(!is_realtime_) reset();
79 
80  uint32 n_channels = x.getNChannels();
81 
82  if(is_realtime_ && n_channels > 1)
83  {
84  M_THROW("In real-time mode, a filter per audio channel must be used!");
85  }
86 
87  AudioStream y(x.getSampleRate(), n_channels);
88 
89  for(uint32 channel = 0; channel < n_channels; ++channel)
90  {
91  y[channel] = filter(x[channel], frequency);
92  }
93 
94  return y;
95 }
unsigned int uint32
Definition: Nsound.h:153
float64 getSampleRate() const
Returns the sample rate of the stream.
Definition: AudioStream.h:217
virtual void reset()=0
uint32 getNChannels(void) const
Returns the number of audio channels in the stream.
Definition: AudioStream.h:212
AudioStream filter(const AudioStream &x)
Definition: Filter.cc:53
#define M_THROW(message)
Definition: Macros.h:108
bool is_realtime_
Definition: Filter.h:118
Buffer Filter::filter ( const Buffer x,
const float64 frequency 
)
inherited

Definition at line 140 of file Filter.cc.

References Nsound::Buffer::begin(), Nsound::Buffer::end(), Nsound::Filter::filter(), Nsound::Buffer::getLength(), Nsound::Filter::is_realtime_, and Nsound::Filter::reset().

141 {
142  if(!is_realtime_) reset();
143 
144  Buffer::const_iterator itor = x.begin();
145  Buffer::const_iterator end = x.end();
146 
147  Buffer y(x.getLength());
148 
149  while(itor != end)
150  {
151  y << filter(*itor, frequency);
152 
153  ++itor;
154  }
155 
156  return y;
157 }
uint32 getLength() const
Returns the number of samples in the Buffer.
Definition: Buffer.h:587
virtual void reset()=0
iterator end()
Retruns the itreator at the end of the Buffer.
Definition: Buffer.h:348
FloatVector::const_iterator const_iterator
Definition: Buffer.h:70
iterator begin()
Retruns the itreator at the start of the Buffer.
Definition: Buffer.h:285
AudioStream filter(const AudioStream &x)
Definition: Filter.cc:53
bool is_realtime_
Definition: Filter.h:118
A Buffer for storing audio samples.
Definition: Buffer.h:60
Buffer Filter::getFrequencyAxis ( const uint32  n_fft = 8192)
inherited

Definition at line 185 of file Filter.cc.

References Nsound::FFTransform::roundUp2(), and Nsound::Filter::sample_rate_.

Referenced by main(), Nsound::Filter::plot(), and Nsound::FilterIIR::savePlot().

186 {
187  uint32 fft_chunk_size = FFTransform::roundUp2(
188  static_cast<int32>(n_fft));
189 
190  uint32 n_samples = fft_chunk_size / 2 + 1;
191 
192  float64 f_step = (1.0 / (static_cast<float64>(fft_chunk_size) / 2.0))
193  * (sample_rate_ / 2.0);
194 
195  Buffer f_axis;
196 
197  float64 f = 0.0;
198 
199  for(uint32 i = 0; i < n_samples; ++i)
200  {
201  f_axis << f;
202  f += f_step;
203  }
204 
205  return f_axis;
206 };
unsigned int uint32
Definition: Nsound.h:153
double float64
Definition: Nsound.h:146
static int32 roundUp2(int32 raw)
Returns nearest power of 2 >= raw.
Definition: FFTransform.cc:274
A Buffer for storing audio samples.
Definition: Buffer.h:60
float64 sample_rate_
Definition: Filter.h:113
Buffer Filter::getFrequencyResponse ( const uint32  n_fft = 8192)
inherited

Definition at line 210 of file Filter.cc.

References Nsound::FFTransform::fft(), Nsound::Filter::getImpulseResponse(), and Nsound::Filter::sample_rate_.

Referenced by Nsound::FilterBandPassIIR::FilterBandPassIIR(), FilterLeastSquaresFIR_UnitTest(), Nsound::FilterIIR::getRMS(), main(), Nsound::Filter::plot(), and Nsound::FilterIIR::savePlot().

211 {
213 
214 //~ fft.setWindow(HANNING);
215 
216  FFTChunkVector vec;
217 
218  vec = fft.fft(getImpulseResponse(), n_fft);
219 
220  return vec[0].getMagnitude();
221 }
Buffer getImpulseResponse(const uint32 n_samples=8192)
Definition: Filter.cc:225
A Class that performes the Fast Fouier Transfrom on a Buffer.
Definition: FFTransform.h:57
std::vector< FFTChunk > FFTChunkVector
Definition: FFTChunk.h:119
float64 sample_rate_
Definition: Filter.h:113
Buffer Filter::getImpulseResponse ( const uint32  n_samples = 8192)
inherited

Definition at line 225 of file Filter.cc.

References Nsound::Filter::filter(), Nsound::Filter::is_realtime_, and Nsound::Filter::reset().

Referenced by Nsound::Filter::getFrequencyResponse(), Nsound::FilterHighPassFIR::getImpulseResponse(), Nsound::FilterLowPassFIR::getImpulseResponse(), Nsound::FilterLeastSquaresFIR::getImpulseResponse(), Nsound::FilterIIR::getImpulseResponse(), and Nsound::Filter::getPhaseResponse().

226 {
227  if(!is_realtime_) reset();
228 
229  Buffer response(n_samples);
230 
231  response << filter(1.0);
232 
233  for(uint32 i = 1; i < n_samples; ++i)
234  {
235  response << filter(0.0);
236  }
237 
238  if(!is_realtime_) reset();
239 
240  return response;
241 }
unsigned int uint32
Definition: Nsound.h:153
virtual void reset()=0
AudioStream filter(const AudioStream &x)
Definition: Filter.cc:53
bool is_realtime_
Definition: Filter.h:118
A Buffer for storing audio samples.
Definition: Buffer.h:60
virtual uint32 Nsound::Filter::getKernelSize ( ) const
inlinevirtualinherited
Buffer Filter::getPhaseResponse ( )
inherited

Definition at line 245 of file Filter.cc.

References Nsound::FFTransform::fft(), Nsound::Filter::getImpulseResponse(), Nsound::Buffer::getLength(), Nsound::Filter::sample_rate_, and Nsound::Buffer::subbuffer().

Referenced by Nsound::Filter::plot().

246 {
247  uint32 n_samples = static_cast<uint32>(sample_rate_ * 2);
248 
249  FFTransform fft(n_samples);
250 
251  FFTChunkVector vec;
252 
253  vec = fft.fft(getImpulseResponse(), n_samples);
254 
255  Buffer phase = vec[0].getPhase();
256 
257  return phase.subbuffer(0, phase.getLength() / 2 + 1);
258 }
Buffer subbuffer(uint32 start_index, uint32 n_samples=0) const
Slice the Buffer.
Definition: Buffer.cc:2073
unsigned int uint32
Definition: Nsound.h:153
uint32 getLength() const
Returns the number of samples in the Buffer.
Definition: Buffer.h:587
Buffer getImpulseResponse(const uint32 n_samples=8192)
Definition: Filter.cc:225
A Class that performes the Fast Fouier Transfrom on a Buffer.
Definition: FFTransform.h:57
A Buffer for storing audio samples.
Definition: Buffer.h:60
std::vector< FFTChunk > FFTChunkVector
Definition: FFTChunk.h:119
float64 sample_rate_
Definition: Filter.h:113
float64 Nsound::Filter::getSampleRate ( ) const
inlineinherited

Definition at line 102 of file Filter.h.

References Nsound::Filter::sample_rate_.

102 { return sample_rate_; };
float64 sample_rate_
Definition: Filter.h:113
void Filter::plot ( boolean  show_phase = false)
inherited

Definition at line 262 of file Filter.cc.

References Nsound::Plotter::figure(), Nsound::Buffer::getdB(), Nsound::Filter::getFrequencyAxis(), Nsound::Filter::getFrequencyResponse(), Nsound::Buffer::getMax(), Nsound::Filter::getPhaseResponse(), Nsound::Plotter::plot(), Nsound::Plotter::subplot(), Nsound::Plotter::xlabel(), Nsound::Plotter::ylabel(), and Nsound::Plotter::ylim().

Referenced by main(), Nsound::FilterLowPassMoogVcf::plot(), Nsound::FilterPhaser::plot(), Nsound::FilterTone::plot(), Nsound::FilterHighPassIIR::plot(), Nsound::FilterCombLowPassFeedback::plot(), Nsound::FilterLowPassIIR::plot(), Nsound::FilterLowPassFIR::plot(), Nsound::FilterHighPassFIR::plot(), Nsound::FilterAllPass::plot(), Nsound::FilterBandPassVocoder::plot(), Nsound::FilterBandPassFIR::plot(), Nsound::FilterLeastSquaresFIR::plot(), Nsound::FilterFlanger::plot(), Nsound::FilterBandRejectFIR::plot(), Nsound::FilterBandPassIIR::plot(), Nsound::FilterBandRejectIIR::plot(), and Nsound::FilterParametricEqualizer::plot().

263 {
264  Buffer x = getFrequencyAxis();
266 
267  Plotter pylab;
268 
269  pylab.figure();
270 
271  uint32 n_rows = 1;
272  uint32 n_cols = 1;
273 
274  if(show_phase)
275  {
276  n_rows = 2;
277  }
278 
279  pylab.subplot(n_rows, n_cols, 1);
280 
281  // Frequency response
282  pylab.plot(x,fr, "blue");
283 
284  pylab.xlabel("Frequency (Hz)");
285  pylab.ylabel("Frequency Response (dB)");
286 
287  // Phase response
288  if(show_phase)
289  {
290  pylab.subplot(n_rows, n_cols, 2);
291 
292  Buffer pr = getPhaseResponse().getdB();
293 
294  pylab.plot(x,pr);
295 
296  pylab.xlabel("Frequency (Hz)");
297  pylab.ylabel("Phase Response (dB)");
298  }
299 
300  float64 ymax = fr.getMax();
301  float64 height = ymax - -60.0;
302 
303  pylab.ylim(-60.0, ymax + 0.05 * height);
304 }
unsigned int uint32
Definition: Nsound.h:153
void xlabel(const std::string &label, const std::string &kwargs="")
Add a label x axis.
Definition: Plotter.cc:1154
void plot(const Buffer &y, const std::string &fmt="", const std::string &kwargs="")
Plots the Buffer on the current figure.
Definition: Plotter.cc:765
void figure(const std::string &kwargs="") const
Creates a new figure window to plot in.
Definition: Plotter.cc:455
Buffer getPhaseResponse()
Definition: Filter.cc:245
double float64
Definition: Nsound.h:146
Axes subplot(const uint32 n_rows, const uint32 n_cols, const uint32 n, const std::string &kwargs="", Axes *sharex=NULL, Axes *sharey=NULL)
Creates a figure in a subplot, subplot(A, B, C, **kwargs)
Definition: Plotter.cc:1031
void ylim(const float64 &ymin, const float64 &ymax)
Definition: Plotter.cc:422
Buffer getFrequencyAxis(const uint32 n_fft=8192)
Definition: Filter.cc:185
void ylabel(const std::string &label, const std::string &kwargs="")
Add a label y axis.
Definition: Plotter.cc:1180
A Buffer for storing audio samples.
Definition: Buffer.h:60
Buffer getFrequencyResponse(const uint32 n_fft=8192)
Definition: Filter.cc:210
float64 getMax() const
Returns the maximum sample value in the Buffer.
Definition: Buffer.cc:951
Buffer getdB() const
Returns the Buffer in dB.
Definition: Buffer.h:487

Member Data Documentation

DelayVector Nsound::FilterSlinky::delay_lines_
protected

Definition at line 90 of file FilterSlinky.h.

Referenced by filter(), FilterSlinky(), reset(), and ~FilterSlinky().

FilterVector Nsound::FilterSlinky::filters_
protected

Definition at line 91 of file FilterSlinky.h.

Referenced by filter(), FilterSlinky(), reset(), and ~FilterSlinky().

float64 Nsound::Filter::sample_rate_
protectedinherited

Definition at line 113 of file Filter.h.

Referenced by Nsound::FilterPhaser::filter(), Nsound::FilterCombLowPassFeedback::filter(), Nsound::FilterDelay::filter(), Nsound::FilterAllPass::FilterAllPass(), Nsound::FilterCombLowPassFeedback::FilterCombLowPassFeedback(), Nsound::FilterDelay::FilterDelay(), Nsound::FilterFlanger::FilterFlanger(), Nsound::FilterPhaser::FilterPhaser(), FilterSlinky(), Nsound::Filter::getFrequencyAxis(), Nsound::Filter::getFrequencyResponse(), Nsound::Filter::getPhaseResponse(), Nsound::Filter::getSampleRate(), Nsound::FilterStageIIR::makeIIRKernelHelper(), Nsound::FilterHighPassFIR::makeKernel(), Nsound::FilterStageIIR::makeKernel(), Nsound::FilterBandPassVocoder::makeKernel(), Nsound::FilterLeastSquaresFIR::makeKernel(), Nsound::FilterParametricEqualizer::makeKernel(), Nsound::FilterPhaser::operator=(), Nsound::FilterLeastSquaresFIR::operator=(), Nsound::FilterFlanger::operator=(), Nsound::FilterIIR::operator=(), Nsound::FilterLowPassMoogVcf::plot(), Nsound::FilterPhaser::plot(), Nsound::FilterTone::plot(), Nsound::FilterHighPassIIR::plot(), Nsound::FilterLowPassIIR::plot(), Nsound::FilterCombLowPassFeedback::plot(), Nsound::FilterLowPassFIR::plot(), Nsound::FilterHighPassFIR::plot(), Nsound::FilterAllPass::plot(), Nsound::FilterBandPassFIR::plot(), Nsound::FilterLeastSquaresFIR::plot(), Nsound::FilterFlanger::plot(), Nsound::FilterBandRejectFIR::plot(), Nsound::FilterBandRejectIIR::plot(), Nsound::FilterBandPassIIR::plot(), Nsound::FilterParametricEqualizer::plot(), Nsound::FilterLowPassIIR::setCutoff(), and Nsound::FilterLeastSquaresFIR::setWindow().

float64 Nsound::Filter::two_pi_over_sample_rate_
protectedinherited
float64 Nsound::Filter::sample_time_
protectedinherited

Definition at line 115 of file Filter.h.

Referenced by Nsound::FilterLowPassMoogVcf::_make_filter().

uint32 Nsound::Filter::kernel_size_
protectedinherited
bool Nsound::Filter::is_realtime_
protectedinherited

The documentation for this class was generated from the following files: