mfcc 0.0.1

checksums.yaml

@@ -0,0 +1,7 @@
+---
+SHA1:
+  metadata.gz: 13ccb341c81931d49deda928b632b66f1c3530ea
+  data.tar.gz: 27c2b78a7c3a4d7927de7e50e4771fad21710bbf
+SHA512:
+  metadata.gz: 917228307517bd92ef1228b8fc52346e783f29e019f5d801436ef8b47c049092ad218b869e9c0102e33ecf593cb994c1ebf71ca1c4742ead1abbb054c72bdbc6
+  data.tar.gz: f462979ee7b29a2f3c434f33550ecd75f1480694243f95c2afdf26179e4ce1b1368a75b4e9cbc904a0c2d8f80a92caeec827919b6b8d1a6d5d6f809d4690905e

data/.gitignore

@@ -0,0 +1,9 @@
+/.bundle/
+/.yardoc
+/Gemfile.lock
+/_yardoc/
+/coverage/
+/doc/
+/pkg/
+/spec/reports/
+/tmp/

data/.rspec

@@ -0,0 +1,2 @@
+--format documentation
+--color

data/.rubocop.yml

@@ -0,0 +1,6 @@
+inherit_from:
+  - .rubocop_todo.yml
+
+AllCops:
+  Exclude:
+    - "*.gemspec"

data/.rubocop_todo.yml

@@ -0,0 +1,15 @@
+# This configuration was generated by `rubocop --auto-gen-config`
+# on 2015-05-25 08:46:53 -0700 using RuboCop version 0.31.0.
+# The point is for the user to remove these configuration records
+# one by one as the offenses are removed from the code base.
+# Note that changes in the inspected code, or installation of new
+# versions of RuboCop, may require this file to be generated again.
+
+# Offense count: 6
+# Configuration parameters: AllowURI, URISchemes.
+Metrics/LineLength:
+  Max: 156
+
+# Offense count: 2
+Style/Documentation:
+  Enabled: false

data/.travis.yml

@@ -0,0 +1,3 @@
+language: ruby
+rvm:
+  - 2.2.1

data/Gemfile

@@ -0,0 +1,4 @@
+source 'https://rubygems.org'
+
+# Specify your gem's dependencies in mfcc.gemspec
+gemspec

data/README.md

@@ -0,0 +1,39 @@
+# Mfcc
+
+Welcome to your new gem! In this directory, you'll find the files you need to be able to package up your Ruby library into a gem. Put your Ruby code in the file `lib/mfcc`. To experiment with that code, run `bin/console` for an interactive prompt.
+
+TODO: Delete this and the text above, and describe your gem
+
+## Installation
+
+Add this line to your application's Gemfile:
+
+```ruby
+gem 'mfcc'
+```
+
+And then execute:
+
+    $ bundle
+
+Or install it yourself as:
+
+    $ gem install mfcc
+
+## Usage
+
+TODO: Write usage instructions here
+
+## Development
+
+After checking out the repo, run `bin/setup` to install dependencies. Then, run `bin/console` for an interactive prompt that will allow you to experiment.
+
+To install this gem onto your local machine, run `bundle exec rake install`. To release a new version, update the version number in `version.rb`, and then run `bundle exec rake release` to create a git tag for the version, push git commits and tags, and push the `.gem` file to [rubygems.org](https://rubygems.org).
+
+## Contributing
+
+1. Fork it ( https://github.com/[my-github-username]/mfcc/fork )
+2. Create your feature branch (`git checkout -b my-new-feature`)
+3. Commit your changes (`git commit -am 'Add some feature'`)
+4. Push to the branch (`git push origin my-new-feature`)
+5. Create a new Pull Request

data/Rakefile

@@ -0,0 +1 @@
+require 'bundler/gem_tasks'

data/bin/console

@@ -0,0 +1,14 @@
+#!/usr/bin/env ruby
+
+require 'bundler/setup'
+require 'mfcc'
+
+# You can add fixtures and/or initialization code here to make experimenting
+# with your gem easier. You can also use a different console, if you like.
+
+# (If you use this, don't forget to add pry to your Gemfile!)
+# require "pry"
+# Pry.start
+
+require 'irb'
+IRB.start

data/bin/setup

@@ -0,0 +1,7 @@
+#!/bin/bash
+set -euo pipefail
+IFS=$'\n\t'
+
+bundle install
+
+# Do any other automated setup that you need to do here

data/lib/mfcc.rb

@@ -0,0 +1,61 @@
+require 'mfcc/version'
+
+# MFCCs are commonly derived as follows:
+#
+# - Take the Fourier transform of (a windowed excerpt of) a signal.
+# - Map the powers of the spectrum obtained above onto the mel scale, using triangular overlapping windows.
+# - Take the logs of the powers at each of the mel frequencies.
+# - Take the discrete cosine transform of the list of mel log powers, as if it were a signal.
+# - The MFCCs are the amplitudes of the resulting spectrum.
+# (http://en.wikipedia.org/wiki/Mel-frequency_cepstrum)
+module Mfcc
+  class Calculator
+    require 'mfcc/preemphasis'
+    require 'mfcc/frame'
+    require 'mfcc/hamming'
+    require 'mfcc/dft'
+    require 'mfcc/mel'
+    require 'mfcc/compressor'
+    require 'mfcc/dct'
+
+    include Mfcc::Preemphasis
+    include Mfcc::Frame
+    include Mfcc::Hamming
+    include Mfcc::Dft
+    include Mfcc::Mel
+    include Mfcc::Compressor
+    include Mfcc::Dct
+
+    attr_reader :data, :options, :frame_size, :frame_step, :alpha, :emphasis,
+                :mel_filters_length, :fft_size, :sample_rate, :dct_order
+
+    def initialize(data, options = {})
+      @data = data
+      @frame_size = options.fetch(:frame_size, 400)
+      @frame_step = options.fetch(:frame_step, 160)
+      @alpha = options.fetch(:alpha, 0.46)
+      @emphasis = options.fetch(:emphasis, 0.97)
+      @mel_filters_length = options.fetch(:mel_filters_length, 20)
+      @fft_size = options.fetch(:fft_size, 512)
+      @sample_rate = options.fetch(:sample_rate, 16_000)
+      @dct_order = options.fetch(:dct_order, 13)
+    end
+
+    def map
+      return to_enum(:map) { self.data.size || Float::Infinity } unless block_given?
+
+      data = preemphasis(self.data)
+
+      data = frame(data)
+
+      data = data.map do |frame|
+        frame = hamming(frame)
+        frame = dft(frame)
+        frame = Mfcc.magnitude(frame)
+        frame = compress(frame)
+        frame = dct(frame)
+        yield frame
+      end
+    end
+  end
+end

data/lib/mfcc/compressor.rb

@@ -0,0 +1,23 @@
+module Mfcc
+  module Compressor
+    def compress(data)
+      Mfcc.compress(data)
+    end
+  end
+
+  LOG10_ERROR_VALUE = -0.00000001
+
+  def self.compress(data)
+    return to_enum(:compress, data) { data.size } unless block_given?
+
+    data.each do |d|
+      v = if d == 0
+            LOG10_ERROR_VALUE
+          else
+            Math.log10(d)
+          end
+
+      yield v
+    end
+  end
+end

data/lib/mfcc/dct.rb

@@ -0,0 +1,33 @@
+module Mfcc
+  # Takes the discrete cosing transform.  Converts an n x m matrix to an n x
+  # order matrix.  ncol should be set to m.
+  module Dct
+    def dct(data)
+      Mfcc.dct(data, dct_order)
+    end
+  end
+
+  def self.dct(data, order = 13, orthogonalize = true)
+    length = data.size
+
+    scales = if orthogonalize
+               [Math.sqrt(1.0 / (4 * length)), Math.sqrt(1.0 / (2 * length))]
+             else
+               [1, 1]
+             end
+
+    dct_matrix(order, length).each_with_index.map do |row, index|
+      scale = index == 0 ? scales[0] : scales[1]
+      scale * row.zip(data).inject(0) { |memo, (a, b)| memo + (2 * a * b) }
+    end
+  end
+
+  def self.dct_matrix(n, m)
+    n.times.map do |i|
+      freq = (Math::PI * i) / m
+      m.times.map do |j|
+        Math.cos(freq * (j + 0.5))
+      end
+    end
+  end
+end

data/lib/mfcc/dft.rb

@@ -0,0 +1,52 @@
+require 'complex'
+
+module Mfcc
+  module Dft
+    def dft(data)
+      Mfcc.dft(data)
+    end
+  end
+
+  def self.dft(data)
+    return to_enum(:dft, data) { data.size } unless block_given?
+
+    n = data.size
+
+    data = plus_imaginary(data)
+
+    data.each_with_index do |_, k|
+      sumreal = 0
+      sumimag = 0
+
+      data.each_with_index do |d, t|
+        angle = 2 * Math::PI * t * k / n
+
+        sumreal += d.real * Math.cos(angle) + d.imaginary * Math.sin(angle)
+        sumimag += -1 * d.real * Math.sin(angle) + d.imaginary * Math.cos(angle)
+      end
+
+      yield Complex(sumreal, sumimag)
+    end
+  end
+
+  def self.plus_imaginary(data)
+    return to_enum(:plus_imaginary, data) { data.size } unless block_given?
+
+    data.each do |d|
+      yield case d
+            when Complex
+              d
+            else
+              Complex(d, 0)
+            end
+    end
+  end
+
+  def self.magnitude(data)
+    return to_enum(:magnitude, data) { data.size } unless block_given?
+
+    data.each do |d|
+      yield d.respond_to?(:magnitude) ? d.magnitude : d
+    end
+  end
+end

data/lib/mfcc/frame.rb

@@ -0,0 +1,33 @@
+module Mfcc
+  module Frame
+    def frame(data)
+      Mfcc.frame(data, frame_size, frame_step)
+    end
+  end
+
+  def self.frame(data, size, step)
+    return to_enum(:frame, data, size, step) unless block_given?
+
+    buffer = []
+
+    data.each do |d|
+      buffer.push(d)
+
+      if buffer.size == size
+        yield buffer
+
+        buffer = buffer.slice(step..buffer.length)
+      end
+    end
+
+    length = buffer.length
+
+    (length / step.to_f).ceil.times do
+      buffer.fill(0, buffer.size...size)
+
+      yield buffer
+
+      buffer = buffer.slice(step..buffer.length)
+    end
+  end
+end

data/lib/mfcc/hamming.rb

@@ -0,0 +1,22 @@
+module Mfcc
+  module Hamming
+    def hamming(data)
+      Mfcc.hamming(data, alpha)
+    end
+  end
+
+  def self.hamming(data, alpha = 0.46)
+    return to_enum(:hamming, data, alpha) { data.size } unless block_given?
+
+    beta = 1 - alpha
+    length = data.size
+
+    data.each_with_index do |d, i|
+      yield d * h(alpha, beta, length, i)
+    end
+  end
+
+  def self.h(alpha, beta, size, i)
+    alpha - beta * Math.cos(2 * Math::PI * i / (size - 1))
+  end
+end

data/lib/mfcc/mel.rb

@@ -0,0 +1,60 @@
+module Mfcc
+  module Mel
+    def mel(data)
+      Mfcc.mel(data, filter_banks)
+    end
+
+    def filter_banks
+      @filter_banks ||= Mfcc::FilterBanks.new(mel_filters_length, fft_size, sample_rate)
+    end
+  end
+
+  def self.mel(data, filter_banks)
+    filter_banks.process(data)
+  end
+
+  class FilterBanks
+    attr_reader :filters, :fft_size, :samplerate, :range
+    def initialize(filters = 20, fft_size = 512, samplerate = 16_000, lowfreq = 0, highfreq = nil)
+      @filters = filters
+      @fft_size = fft_size
+      @samplerate = samplerate
+      @range = lowfreq..(highfreq || (samplerate / 2))
+    end
+
+    def process(fft)
+      (0..filters).map { |n| fft.zip(bank[n]).reject { |(_, b)| b.nil? }.inject(0) { |memo, (a, b)| memo + a * b } }
+    end
+
+    def bank(i)
+      lower = fft_frequencies.map { |x| (x - mel_frequencies[i]) / (mel_frequencies[i + 1] - mel_frequencies[i]).to_f }
+      upper = fft_frequencies.map { |x| (mel_frequencies[i + 2] - x) / (mel_frequencies[i + 2] - mel_frequencies[i + 1]).to_f }
+
+      (0..(fft_size / 2)).map { |n| [0, [lower[n], upper[n]].min].max }
+    end
+
+    def r(i)
+      mel_frequencies[i]
+    end
+
+    def to_hz(mel)
+      700 * (10**(mel / 2595.0) - 1)
+    end
+
+    def to_mel(hz)
+      2595 * Math.log10(1 + hz / 700.0)
+    end
+
+    def mel_frequencies
+      @mel_frequencies ||= range.step(mel_step).to_a
+    end
+
+    def mel_step
+      ((range.end - range.begin) / filters)
+    end
+
+    def fft_frequencies
+      @fft_frequencies ||= [0] + (fft_size / 2).times.map { |i| (i + 1) * (samplerate / fft_size.to_f) }
+    end
+  end
+end

data/lib/mfcc/preemphasis.rb

@@ -0,0 +1,20 @@
+module Mfcc
+  module Preemphasis
+    def preemphasis(data)
+      Mfcc.preemphasis(data, emphasis)
+    end
+  end
+
+  def self.preemphasis(data, emph = 0.97)
+    return data if emph == 0
+
+    return to_enum(:preemphasis, data, emph) { |d, _| d.size || Float::INFINITY } unless block_given?
+
+    prev = 0
+
+    data.each do |y|
+      yield (y - emph * prev)
+      prev = y
+    end
+  end
+end

data/lib/mfcc/version.rb

@@ -0,0 +1,3 @@
+module Mfcc
+  VERSION = '0.0.1'
+end

data/mfcc.gemspec

@@ -0,0 +1,23 @@
+# coding: utf-8
+lib = File.expand_path('../lib', __FILE__)
+$LOAD_PATH.unshift(lib) unless $LOAD_PATH.include?(lib)
+require 'mfcc/version'
+
+Gem::Specification.new do |spec|
+  spec.name          = 'mfcc'
+  spec.version       = Mfcc::VERSION
+  spec.authors       = ['Chris Beer']
+  spec.email         = ['chris@cbeer.info']
+
+  spec.summary       = 'Mel-frequency cepstral coefficients for Ruby'
+  spec.homepage      = 'https://github.com/cbeer/mfcc'
+
+  spec.files         = `git ls-files -z`.split("\x0").reject { |f| f.match(%r{^(test|spec|features)/}) }
+  spec.bindir        = 'exe'
+  spec.executables   = spec.files.grep(%r{^exe/}) { |f| File.basename(f) }
+  spec.require_paths = ['lib']
+
+  spec.add_development_dependency 'bundler', '~> 1.9'
+  spec.add_development_dependency 'rake', '~> 10.0'
+  spec.add_development_dependency 'rspec', '~> 3.0'
+end

metadata

@@ -0,0 +1,105 @@
+--- !ruby/object:Gem::Specification
+name: mfcc
+version: !ruby/object:Gem::Version
+  version: 0.0.1
+platform: ruby
+authors:
+- Chris Beer
+autorequire: 
+bindir: exe
+cert_chain: []
+date: 2015-05-28 00:00:00.000000000 Z
+dependencies:
+- !ruby/object:Gem::Dependency
+  name: bundler
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '1.9'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '1.9'
+- !ruby/object:Gem::Dependency
+  name: rake
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '10.0'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '10.0'
+- !ruby/object:Gem::Dependency
+  name: rspec
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '3.0'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '3.0'
+description: 
+email:
+- chris@cbeer.info
+executables: []
+extensions: []
+extra_rdoc_files: []
+files:
+- ".gitignore"
+- ".rspec"
+- ".rubocop.yml"
+- ".rubocop_todo.yml"
+- ".travis.yml"
+- Gemfile
+- README.md
+- Rakefile
+- bin/console
+- bin/setup
+- lib/mfcc.rb
+- lib/mfcc/compressor.rb
+- lib/mfcc/dct.rb
+- lib/mfcc/dft.rb
+- lib/mfcc/frame.rb
+- lib/mfcc/hamming.rb
+- lib/mfcc/mel.rb
+- lib/mfcc/preemphasis.rb
+- lib/mfcc/version.rb
+- mfcc.gemspec
+homepage: https://github.com/cbeer/mfcc
+licenses: []
+metadata: {}
+post_install_message: 
+rdoc_options: []
+require_paths:
+- lib
+required_ruby_version: !ruby/object:Gem::Requirement
+  requirements:
+  - - ">="
+    - !ruby/object:Gem::Version
+      version: '0'
+required_rubygems_version: !ruby/object:Gem::Requirement
+  requirements:
+  - - ">="
+    - !ruby/object:Gem::Version
+      version: '0'
+requirements: []
+rubyforge_project: 
+rubygems_version: 2.4.5
+signing_key: 
+specification_version: 4
+summary: Mel-frequency cepstral coefficients for Ruby
+test_files: []