spcore 0.2.0 → 0.2.1

metadata

@@ -1,7 +1,7 @@
 --- !ruby/object:Gem::Specification
 name: spcore
 version: !ruby/object:Gem::Version
-  version: 0.2.0
+  version: 0.2.1
   prerelease: 
 platform: ruby
 authors:
@@ -9,7 +9,7 @@ authors:
 autorequire: 
 bindir: bin
 cert_chain: []
-date: 2013-05-21 00:00:00.000000000 Z
+date: 2013-07-03 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: hashmake
@@ -123,6 +123,38 @@ dependencies:
     - - ! '>='
       - !ruby/object:Gem::Version
         version: '0'
+- !ruby/object:Gem::Dependency
+  name: pry-nav
+  requirement: !ruby/object:Gem::Requirement
+    none: false
+    requirements:
+    - - ! '>='
+      - !ruby/object:Gem::Version
+        version: '0'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    none: false
+    requirements:
+    - - ! '>='
+      - !ruby/object:Gem::Version
+        version: '0'
+- !ruby/object:Gem::Dependency
+  name: wavefile
+  requirement: !ruby/object:Gem::Requirement
+    none: false
+    requirements:
+    - - ! '>='
+      - !ruby/object:Gem::Version
+        version: '0'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    none: false
+    requirements:
+    - - ! '>='
+      - !ruby/object:Gem::Version
+        version: '0'
 description: ! 'A library of signal processing methods and classes for resampling,
   frequency domain transform (FFT and DFT), windowing (Blackman, Hamming, etc.), FIR
   filtering (windowed sinc), interpolation (linear and polynomial), plotting data,
@@ -144,14 +176,13 @@ files:
 - README.rdoc
 - Rakefile
 - lib/spcore.rb
-- lib/spcore/analysis/correlation.rb
-- lib/spcore/analysis/envelope.rb
-- lib/spcore/analysis/extrema.rb
-- lib/spcore/analysis/signal.rb
+- lib/spcore/analysis/calculus.rb
+- lib/spcore/analysis/features.rb
+- lib/spcore/analysis/frequency_domain.rb
+- lib/spcore/analysis/statistics.rb
 - lib/spcore/core/circular_buffer.rb
 - lib/spcore/core/constants.rb
 - lib/spcore/core/delay_line.rb
-- lib/spcore/core/oscillator.rb
 - lib/spcore/filters/fir/dual_sinc_filter.rb
 - lib/spcore/filters/fir/fir.rb
 - lib/spcore/filters/fir/sinc_filter.rb
@@ -161,6 +192,9 @@ files:
 - lib/spcore/filters/iir/cookbook_highpass_filter.rb
 - lib/spcore/filters/iir/cookbook_lowpass_filter.rb
 - lib/spcore/filters/iir/cookbook_notch_filter.rb
+- lib/spcore/generation/comb_filter.rb
+- lib/spcore/generation/oscillator.rb
+- lib/spcore/generation/signal_generator.rb
 - lib/spcore/interpolation/interpolation.rb
 - lib/spcore/resampling/discrete_resampling.rb
 - lib/spcore/resampling/hybrid_resampling.rb
@@ -173,7 +207,7 @@ files:
 - lib/spcore/util/plotter.rb
 - lib/spcore/util/saturation.rb
 - lib/spcore/util/scale.rb
-- lib/spcore/util/signal_generator.rb
+- lib/spcore/util/signal.rb
 - lib/spcore/version.rb
 - lib/spcore/windows/bartlett_hann_window.rb
 - lib/spcore/windows/bartlett_window.rb
@@ -191,16 +225,20 @@ files:
 - lib/spcore/windows/triangular_window.rb
 - lib/spcore/windows/tukey_window.rb
 - spcore.gemspec
-- spec/analysis/correlation_spec.rb
-- spec/analysis/envelope_spec.rb
-- spec/analysis/extrema_spec.rb
-- spec/analysis/signal_spec.rb
+- spec/analysis/calculus_spec.rb
+- spec/analysis/features_spec.rb
+- spec/analysis/frequency_domain_spec.rb
+- spec/analysis/piano_C4.wav
+- spec/analysis/statistics_spec.rb
+- spec/analysis/trumpet_B4.wav
 - spec/core/circular_buffer_spec.rb
 - spec/core/delay_line_spec.rb
-- spec/core/oscillator_spec.rb
 - spec/filters/fir/dual_sinc_filter_spec.rb
 - spec/filters/fir/sinc_filter_spec.rb
 - spec/filters/iir/cookbook_filter_spec.rb
+- spec/generation/comb_filter_spec.rb
+- spec/generation/oscillator_spec.rb
+- spec/generation/signal_generator_spec.rb
 - spec/interpolation/interpolation_spec.rb
 - spec/resampling/discrete_resampling_spec.rb
 - spec/resampling/hybrid_resampling_spec.rb
@@ -213,7 +251,7 @@ files:
 - spec/util/gain_spec.rb
 - spec/util/limiters_spec.rb
 - spec/util/saturate_spec.rb
-- spec/util/signal_generator_spec.rb
+- spec/util/signal_spec.rb
 - spec/windows/window_spec.rb
 homepage: https://github.com/jamestunnell/spcore
 licenses:
@@ -230,7 +268,7 @@ required_ruby_version: !ruby/object:Gem::Requirement
       version: '0'
       segments:
       - 0
-      hash: 688071795
+      hash: 2316548233979027039
 required_rubygems_version: !ruby/object:Gem::Requirement
   none: false
   requirements:
@@ -239,7 +277,7 @@ required_rubygems_version: !ruby/object:Gem::Requirement
       version: '0'
       segments:
       - 0
-      hash: 688071795
+      hash: 2316548233979027039
 requirements: []
 rubyforge_project: 
 rubygems_version: 1.8.23
@@ -247,16 +285,20 @@ signing_key:
 specification_version: 3
 summary: A library of signal processing methods and classes.
 test_files:
-- spec/analysis/correlation_spec.rb
-- spec/analysis/envelope_spec.rb
-- spec/analysis/extrema_spec.rb
-- spec/analysis/signal_spec.rb
+- spec/analysis/calculus_spec.rb
+- spec/analysis/features_spec.rb
+- spec/analysis/frequency_domain_spec.rb
+- spec/analysis/piano_C4.wav
+- spec/analysis/statistics_spec.rb
+- spec/analysis/trumpet_B4.wav
 - spec/core/circular_buffer_spec.rb
 - spec/core/delay_line_spec.rb
-- spec/core/oscillator_spec.rb
 - spec/filters/fir/dual_sinc_filter_spec.rb
 - spec/filters/fir/sinc_filter_spec.rb
 - spec/filters/iir/cookbook_filter_spec.rb
+- spec/generation/comb_filter_spec.rb
+- spec/generation/oscillator_spec.rb
+- spec/generation/signal_generator_spec.rb
 - spec/interpolation/interpolation_spec.rb
 - spec/resampling/discrete_resampling_spec.rb
 - spec/resampling/hybrid_resampling_spec.rb
@@ -269,6 +311,6 @@ test_files:
 - spec/util/gain_spec.rb
 - spec/util/limiters_spec.rb
 - spec/util/saturate_spec.rb
-- spec/util/signal_generator_spec.rb
+- spec/util/signal_spec.rb
 - spec/windows/window_spec.rb
 has_rdoc: 

data/lib/spcore/analysis/envelope.rb

@@ -1,76 +0,0 @@
-module SPCore
-# Determines the envelope of given samples. Unlike the EnvelopeDetector, this
-# operates on a signal (time-series data) after it is recieved, not
-# sample-by-sample. As a result, it provides the envelope as an entire signal.
-#
-# @author James Tunnell
-class Envelope
-  
-  attr_reader :data
-  
-  def initialize samples
-    # combine absolute values of positive maxima and negative minima
-    extrema = Extrema.new(samples)
-    starting_outline = {}
-    extrema.minima.each do |idx,val|
-      if val <= 0.0
-        starting_outline[idx] =val.abs
-      end
-    end
-    
-    extrema.maxima.each do |idx,val|
-      if val >= 0.0
-        starting_outline[idx] = val.abs
-      end
-    end
-    
-    # add in first and last samples so the envelope follows entire signal
-    starting_outline[0] = samples[0].abs
-    starting_outline[samples.count - 1] = samples[samples.count - 1].abs
-    
-    # the extrema we have now are probably not spaced evenly. Upsampling at
-    # this point would lead to a time-distorted signal. So the next step is to
-    # interpolate between all the extrema to make a crude but properly sampled
-    # envelope.
-    
-    proper_outline = Array.new(samples.count, 0)
-    indices = starting_outline.keys.sort
-    
-    for i in 1...indices.count
-      l_idx = indices[i-1]
-      r_idx = indices[i]
-      
-      l_val = starting_outline[l_idx]
-      r_val = starting_outline[r_idx]
-      
-      proper_outline[l_idx] = l_val
-      proper_outline[r_idx] = r_val
-      
-      idx_span = r_idx - l_idx
-      
-      for j in (l_idx + 1)...(r_idx)
-        x = (j - l_idx).to_f / idx_span
-        y = Interpolation.linear l_val, r_val, x
-        proper_outline[j] = y
-      end
-    end
-    
-    # Now downsample by dropping samples, back to the number of starting_outline we had
-    # with just the extrema, but this time with samples properly spaced so as
-    # to avoid time distortion after upsampling.
-    
-    downsample_factor = (samples.count / starting_outline.count).to_i
-    downsampled_outline = []
-
-    (0...proper_outline.count).step(downsample_factor) do |n|
-      downsampled_outline.push proper_outline[n]
-    end
-    
-    # finally, use polynomial interpolation to upsample to the original sample rate.
-    
-    upsample_factor = samples.count / downsampled_outline.count.to_f
-    @data = PolynomialResampling.upsample(downsampled_outline, upsample_factor)
-  end
-
-end
-end

data/lib/spcore/analysis/extrema.rb

@@ -1,55 +0,0 @@
-module SPCore
-# Finds extrema (minima and maxima).
-class Extrema
-  
-  attr_reader :minima, :maxima, :extrema
-  
-  def initialize samples
-    @minima = {}
-    @maxima = {}
-    
-    global_min_idx = 0
-    global_min_val = samples[0]
-    global_max_idx = 0
-    global_max_val = samples[0]
-    
-    diffs = []
-    for i in (1...samples.count)
-      diffs.push(samples[i] - samples[i-1])
-      
-      if samples[i] < global_min_val
-        global_min_idx = i
-        global_min_val = samples[i]
-      end
-      
-      if samples[i] > global_max_val
-        global_max_idx = i
-        global_max_val = samples[i]
-      end
-    end
-    @minima[global_min_idx] = global_min_val
-    @maxima[global_max_idx] = global_max_val
-    
-    is_positive = diffs.first > 0.0 # starting off with positive difference?
-    
-    # at zero crossings there is a local maxima/minima    
-    for i in (1...diffs.count)
-      if is_positive
-        # at positive-to-negative transition there is a local maxima
-        if diffs[i] <= 0.0
-          @maxima[i] = samples[i]
-          is_positive = false
-        end
-      else
-        # at negative-to-positive transition there is a local minima
-        if diffs[i] > 0.0
-          @minima[i] = samples[i]
-          is_positive = true
-        end
-      end
-    end
-    
-    @extrema = @minima.merge(@maxima)    
-  end
-end 
-end

data/spec/analysis/correlation_spec.rb

@@ -1,28 +0,0 @@
-require File.expand_path(File.dirname(__FILE__) + '/../spec_helper')
-
-describe SPCore::Correlation do
-  context 'image => triangular window' do
-    before :all do
-      @size = size = 48
-      @triangle = TriangularWindow.new(size * 2).data
-    end
-    
-    context 'feature => rising ramp (half size of triangular window)' do
-      it 'should have maximum correlation at beginning' do
-        rising_ramp = Array.new(@size){|i| i / @size.to_f }
-        correlation = Correlation.new(@triangle, rising_ramp)
-        correlation.data.first.should eq(correlation.data.max)
-      end
-    end
-    
-    context 'feature => falling ramp (half size of triangular window)' do
-      it 'should have maximum correlation at end' do
-        falling_ramp = Array.new(@size){|i| (@size - i) / @size.to_f }
-        correlation = Correlation.new(@triangle, falling_ramp)
-        correlation.data.last.should eq(correlation.data.max)
-        
-        #Plotter.plot_1d "correlate falling ramp" => correlation.data
-      end
-    end
-  end  
-end

data/spec/analysis/envelope_spec.rb

@@ -1,50 +0,0 @@
-require File.expand_path(File.dirname(__FILE__) + '/../spec_helper')
-
-describe SPCore::Envelope do
-  before :all do
-    sample_rate = 1000
-    sample_count = 512 * 2
-    generator = SignalGenerator.new(:size => sample_count, :sample_rate => sample_rate)
-    
-    @modulation_signal = generator.make_signal [4.0], :amplitude => 0.1
-    @modulation_signal.multiply! BlackmanWindow.new(sample_count).data
-
-    @base_signal = generator.make_signal [64.0]
-    @base_signal.multiply! @modulation_signal
-  end
-
-  it 'should produce an output that follows the amplitude of the input' do
-    envelope = @base_signal.envelope
-    check_envelope(envelope)
-  end
-  
-  def check_envelope envelope    
-    #signals = {
-    #  "signal" => @base_signal,
-    #  "modulation (abs)" => @modulation_signal.abs,
-    #  "envelope" => envelope,
-    #}
-    #
-    #Plotter.new(
-    #  :title => "signal and envelope",
-    #  :xlabel => "sample",
-    #  :ylabel => "values",
-    #).plot_signals(signals)
-    
-    #Plotter.new().plot_2d("envelop freq magnitudes" => envelope.freq_magnitudes)
-    
-    begin
-      ideal = @modulation_signal.energy
-      actual = envelope.energy
-      error = (ideal - actual).abs / ideal
-      error.should be_within(0.1).of(0.0)
-    end
-    
-    begin
-      ideal = @modulation_signal.rms
-      actual = envelope.rms
-      error = (ideal - actual).abs / ideal
-      error.should be_within(0.1).of(0.0)
-    end
-  end
-end

data/spec/analysis/extrema_spec.rb

@@ -1,42 +0,0 @@
-require File.expand_path(File.dirname(__FILE__) + '/../spec_helper')
-
-describe SPCore::Extrema do
-  context '#minima' do
-    it 'should return points where local and global minima occur' do
-      cases = {
-        [3.8, 3.0, 2.9, 2.95, 3.6, 3.4, 2.8, 2.3, 2.1, 2.0, 2.5] => { 2 => 2.9, 9 => 2.0 },
-        [3.2, 3.5, 2.9, 2.7, 2.8, 2.7, 2.5, 2.2, 2.4, 2.3, 2.0] => { 3 => 2.7, 7 => 2.2, 10 => 2.0 },
-      }
-        
-      cases.each do |samples, minima|
-        Extrema.new(samples).minima.should eq minima
-      end
-    end
-  end
-  
-  context '#maxima' do
-    it 'should return points where local and global maxima occur' do
-      cases = {
-        [3.8, 3.0, 2.9, 2.95, 3.6, 3.4, 2.8, 2.3, 2.1, 2.0, 2.5] => { 0 => 3.8, 4 => 3.6 },
-        [3.2, 3.5, 2.9, 2.7, 2.8, 2.7, 2.5, 2.2, 2.4, 2.3, 2.0] => { 1 => 3.5, 4 => 2.8, 8 => 2.4},
-      }
-        
-      cases.each do |samples, maxima|
-        Extrema.new(samples).maxima.should eq maxima
-      end
-    end
-  end
-  
-  context '#extrema' do
-    it 'should return points where local and global extrema occur' do
-      cases = {
-        [3.8, 3.0, 2.9, 2.95, 3.6, 3.4, 2.8, 2.3, 2.1, 2.0, 2.5] => { 0 => 3.8, 2 => 2.9, 4 => 3.6, 9 => 2.0},
-        [3.2, 3.5, 2.9, 2.7, 2.8, 2.7, 2.5, 2.2, 2.4, 2.3, 2.0] => { 1 => 3.5, 3 => 2.7, 4 => 2.8, 7 => 2.2, 8 => 2.4, 10 => 2.0},
-      }
-        
-      cases.each do |samples, extrema|
-        Extrema.new(samples).extrema.should eq extrema
-      end
-    end
-  end
-end