HelpSource/Classes/SpecPcile.schelp

   1 class:: SpecPcile
   2 summary:: Find a percentile of FFT magnitude spectrum
   3 categories:: UGens>FFT
   4 related:: Classes/SpecCentroid, Classes/SpecFlatness
   5
   6 description::
   7 Given an link::Classes/FFT:: chain this calculates the cumulative distribution of the frequency spectrum, and outputs the frequency value which corresponds to the desired percentile.
   8
   9 For example, to find the frequency at which 90% of the spectral energy lies below that frequency, you want the 90-percentile, which means the value of emphasis::fraction:: should be 0.9. The 90-percentile or 95-percentile is often used as a measure of strong::spectral roll-off::.
  10
  11 The optional third argument strong::interpolate:: specifies whether interpolation should be used to try and make the percentile frequency estimate more accurate, at the cost of a little higher CPU usage. Set it to 1 to enable this.
  12
  13 classmethods::
  14 method:: kr
  15
  16 argument:: chain
  17 an link::Classes/FFT:: chain.
  18 argument:: fraction
  19 argument:: interpolate
  20
  21 examples::
  22
  23 code::
  24 s = Server.internal.boot;
  25 b = Buffer.alloc(s,2048,1);
  26
  27 // Simple demo with filtering white noise, and trying to infer the cutoff freq.
  28 // Move the mouse.
  29 (
  30 {
  31 var in, chain, realcutoff, estcutoff;
  32 realcutoff = MouseX.kr(0.00001,22050);
  33 in = LPF.ar(WhiteNoise.ar, realcutoff);
  34 chain = FFT(b, in);
  35 estcutoff = Lag.kr(SpecPcile.kr(chain, 0.9), 1);
  36 realcutoff.poll(Impulse.kr(1), "real cutoff");
  37 estcutoff.poll(Impulse.kr(1), "estimated cutoff");
  38 Out.ar(0, in);
  39 Out.kr(0, estcutoff * 22050.0.reciprocal);
  40 }.scope;
  41 )
  42
  43 // Audio input - try different vowel/long-consonant sounds and see what comes out.
  44 // Specifically, change from "ssss" through to "aaaa" through to "wwww".
  45 (
  46 {
  47 var in, chain, perc;
  48 in = SoundIn.ar([0,1]).mean;
  49 chain = FFT(b, in);
  50 //Out.ar(0, in * 0.1);
  51 perc = SpecPcile.kr(chain, 0.5);
  52 Out.ar(1, LPF.ar(WhiteNoise.ar, perc)); //NB Outputting to right channel - handy on PowerBooks
  53 Out.kr(0, perc * 22050.0.reciprocal);
  54 }.scope;
  55 )
  56 ::