digiproc 0.1.0

checksums.yaml

@@ -0,0 +1,7 @@
+---
+SHA256:
+  metadata.gz: a9beaf357290bb8d1ddfb103b8d39614ce93539ead7353e2053faf8a3d3c36c3
+  data.tar.gz: cf4a4a6c03e78ef312109654871c294c372f22bd9437813f61de092ae9453f89
+SHA512:
+  metadata.gz: 6891122f07044a1f0c20b424d5c16753f67ef4dbc71458ada887835bf39747d616bb6f1914a0c5d8fae155c1c74188daf601e88081c9922111ff11ec8b74f387
+  data.tar.gz: 5a39c3fdba680ce91dec8122f8d9fd48959bc9df550bdce413d2b8efabe67f6d0feab00919f360a1fb2247ad9e15c4399b0af86bf771da7d1c7c514d97e2a851

data/.gitignore

@@ -0,0 +1,12 @@
+/.bundle/
+/.yardoc
+/_yardoc/
+/coverage/
+/doc/
+/pkg/
+/spec/reports/
+/tmp/
+
+# rspec failure tracking
+.rspec_status
+.DS_Store

data/.rspec

@@ -0,0 +1,3 @@
+--format documentation
+--color
+--require spec_helper

data/.travis.yml

@@ -0,0 +1,7 @@
+---
+sudo: false
+language: ruby
+cache: bundler
+rvm:
+  - 2.5.1
+before_install: gem install bundler -v 1.16.6

data/CODE_OF_CONDUCT.md

@@ -0,0 +1,74 @@
+# Contributor Covenant Code of Conduct
+
+## Our Pledge
+
+In the interest of fostering an open and welcoming environment, we as
+contributors and maintainers pledge to making participation in our project and
+our community a harassment-free experience for everyone, regardless of age, body
+size, disability, ethnicity, gender identity and expression, level of experience,
+nationality, personal appearance, race, religion, or sexual identity and
+orientation.
+
+## Our Standards
+
+Examples of behavior that contributes to creating a positive environment
+include:
+
+* Using welcoming and inclusive language
+* Being respectful of differing viewpoints and experiences
+* Gracefully accepting constructive criticism
+* Focusing on what is best for the community
+* Showing empathy towards other community members
+
+Examples of unacceptable behavior by participants include:
+
+* The use of sexualized language or imagery and unwelcome sexual attention or
+advances
+* Trolling, insulting/derogatory comments, and personal or political attacks
+* Public or private harassment
+* Publishing others' private information, such as a physical or electronic
+  address, without explicit permission
+* Other conduct which could reasonably be considered inappropriate in a
+  professional setting
+
+## Our Responsibilities
+
+Project maintainers are responsible for clarifying the standards of acceptable
+behavior and are expected to take appropriate and fair corrective action in
+response to any instances of unacceptable behavior.
+
+Project maintainers have the right and responsibility to remove, edit, or
+reject comments, commits, code, wiki edits, issues, and other contributions
+that are not aligned to this Code of Conduct, or to ban temporarily or
+permanently any contributor for other behaviors that they deem inappropriate,
+threatening, offensive, or harmful.
+
+## Scope
+
+This Code of Conduct applies both within project spaces and in public spaces
+when an individual is representing the project or its community. Examples of
+representing a project or community include using an official project e-mail
+address, posting via an official social media account, or acting as an appointed
+representative at an online or offline event. Representation of a project may be
+further defined and clarified by project maintainers.
+
+## Enforcement
+
+Instances of abusive, harassing, or otherwise unacceptable behavior may be
+reported by contacting the project team at micah.shute@gmail.com. All
+complaints will be reviewed and investigated and will result in a response that
+is deemed necessary and appropriate to the circumstances. The project team is
+obligated to maintain confidentiality with regard to the reporter of an incident.
+Further details of specific enforcement policies may be posted separately.
+
+Project maintainers who do not follow or enforce the Code of Conduct in good
+faith may face temporary or permanent repercussions as determined by other
+members of the project's leadership.
+
+## Attribution
+
+This Code of Conduct is adapted from the [Contributor Covenant][homepage], version 1.4,
+available at [http://contributor-covenant.org/version/1/4][version]
+
+[homepage]: http://contributor-covenant.org
+[version]: http://contributor-covenant.org/version/1/4/

data/Gemfile

@@ -0,0 +1,8 @@
+source "https://rubygems.org"
+
+git_source(:github) {|repo_name| "https://github.com/#{repo_name}" }
+
+# Specify your gem's dependencies in digiproc.gemspec
+gemspec
+gem 'gruff'
+

data/Gemfile.lock

@@ -0,0 +1,48 @@
+PATH
+  remote: .
+  specs:
+    digiproc (0.1.0)
+      gruff (~> 0.7.0)
+
+GEM
+  remote: https://rubygems.org/
+  specs:
+    coderay (1.1.2)
+    diff-lcs (1.3)
+    gruff (0.7.0)
+      rmagick (~> 2.13, >= 2.13.4)
+    method_source (0.9.2)
+    pry (0.12.2)
+      coderay (~> 1.1.0)
+      method_source (~> 0.9.0)
+    rake (10.5.0)
+    rdoc (6.1.1)
+    rmagick (2.16.0)
+    rspec (3.8.0)
+      rspec-core (~> 3.8.0)
+      rspec-expectations (~> 3.8.0)
+      rspec-mocks (~> 3.8.0)
+    rspec-core (3.8.0)
+      rspec-support (~> 3.8.0)
+    rspec-expectations (3.8.3)
+      diff-lcs (>= 1.2.0, < 2.0)
+      rspec-support (~> 3.8.0)
+    rspec-mocks (3.8.0)
+      diff-lcs (>= 1.2.0, < 2.0)
+      rspec-support (~> 3.8.0)
+    rspec-support (3.8.0)
+
+PLATFORMS
+  ruby
+
+DEPENDENCIES
+  bundler
+  digiproc!
+  gruff
+  pry
+  rake (~> 10.0)
+  rdoc
+  rspec (~> 3.0)
+
+BUNDLED WITH
+   2.0.1

data/LICENSE.txt

@@ -0,0 +1,21 @@
+The MIT License (MIT)
+
+Copyright (c) 2018 micahshute
+
+Permission is hereby granted, free of charge, to any person obtaining a copy
+of this software and associated documentation files (the "Software"), to deal
+in the Software without restriction, including without limitation the rights
+to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
+copies of the Software, and to permit persons to whom the Software is
+furnished to do so, subject to the following conditions:
+
+The above copyright notice and this permission notice shall be included in
+all copies or substantial portions of the Software.
+
+THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN
+THE SOFTWARE.

data/README.md

@@ -0,0 +1,78 @@
+# Digiproc
+
+Use Digiproc to perform a number of digital signal processing calculations. Easily represent digital and analog signals, find their Discrete Fourier Transform in O(n*lgn) time, and create various systems from filtering, modulation, encoding, etc. Digiproc also enables you to easily solve a wide range of probability and linear algebra-based problems. 
+
+Easily visualize data in a graphing API built on top of gruff.
+
+See the `examples` folder as well as the `spec`s. 
+
+## Installation
+
+Add this line to your application's Gemfile:
+
+```ruby
+gem 'digiproc'
+```
+
+And then execute:
+
+    $ bundle
+
+Or install it yourself as:
+
+    $ gem install digiproc
+
+If you are having trouble downloading the rmagick gem on OSX, try the following:
+
+```
+    brew uninstall imagemagick
+    brew install imagemagick@6
+    export PATH="/usr/local/opt/imagemagick@6/bin:$PATH"
+    brew link --force imagemagick@6
+    gem install rmagick
+```
+
+## Usage
+
+Graphing capabillity built on top of [gruff](https://github.com/topfunky/gruff) gem, which uses the [rMagick](https://rubygems.org/gems/rmagick/versions/2.15.4) gem and ImageMagick. 
+
+__Look at specs and the `examples` folder for examples on how to use these classes in calculations__
+
+__Limitations to be aware of__:
+ - All signals must be causual -> negative values of n are not supported
+ - The FFT, as of now, uses a Radix 2 algorithm. So, the overall size must be a power of 2 - any dataset is zero-filled automatically to meet this necessity. This gives you less control over the exact size if needed and may cause slower runtimes (however the radix 2 algorithm is ϴ(nlgn)). This zero-fill _will_ increase the resolution of the FFT output, however.
+ - All filter data and windows must have an even number of datapoints. If not, this will be done automatically. This is to facilitate the ability to create any type of filter using the windowing method. (Odd numbers of values can preclude certain types of filters)
+ - All filters are FIR and are implemented via the windowing method. Although you can get pretty good results depending on your application, this is an older method compared to the optimal Parks-McClellan algorithm.
+ - IIR Filter design is not included
+ - The quantization process in Digiproc::AnalogSignal maps an analog signal to -1 to 1 before quantizing the result. A floating point number is accurate up to 7 decimal places, so this process will cause unwanted ACTUAL quantization errors (not simulated quantization errors) for an amplitude range greater than $2 x 10^16$ (ie 20 Quadrillion)
+
+This is still a work in progress. There are many usuable functions, and many that still require tuning. Functions tested in the `spec` file are very reliable.
+
+
+
+__Run Examples__
+
+Use ```rake examples:run[file]``` to run an example file to see how it works, ie: ```rake examples:run[binomial_distribution/dice]```
+
+
+## Development
+
+After checking out the repo, run `bin/setup` to install dependencies. Then, run `rake spec` to run the tests. You can also 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`, which will create a git tag for the version, push git commits and tags, and push the `.gem` file to [rubygems.org](https://rubygems.org).
+
+## Contributing
+
+Bug reports and pull requests are welcome on GitHub at https://github.com/micahshute/digiproc. 
+
+There is a long list of intended actions in the `TODO.md` file. Tests (rspec) covers many core capabilities, but more need to be written to cover the entire gem.
+
+This project is intended to be a safe, welcoming space for collaboration, and contributors are expected to adhere to the [Contributor Covenant](http://contributor-covenant.org) code of conduct. That being said, negative feedback is always welcome. 
+
+## License
+
+The gem is available as open source under the terms of the [MIT License](https://opensource.org/licenses/MIT).
+
+## Code of Conduct
+
+Everyone interacting in the Digiproc project’s codebases, issue trackers, chat rooms and mailing lists is expected to follow the [code of conduct](https://github.com/[USERNAME]/digiproc/blob/master/CODE_OF_CONDUCT.md).

data/Rakefile

@@ -0,0 +1,37 @@
+require "bundler/gem_tasks"
+require "rspec/core/rake_task"
+
+RSpec::Core::RakeTask.new(:spec)
+
+task :default => :spec
+
+
+task :environment do 
+    require_relative './config/environment'
+end
+
+desc "Test data structures in console environment"
+task :console => :environment do
+    require 'pry'
+    pry.start
+end
+
+desc "Loads data for console tests, runs console"
+task :console_tests => :environment do
+    require 'pry'
+    require_relative './console_tests.rb'
+ 
+end
+
+desc "Env playground, same as console tests but not for testing"
+task :play => :environment do 
+    require 'pry'
+    require_relative './playground.rb'
+end
+
+namespace :examples do 
+    desc "run an example file"
+    task :run, [:file] => :environment do |task, args|
+        require_relative "./examples/#{args[:file]}.rb"
+    end
+end

data/TODO.md

@@ -0,0 +1,50 @@
+# TODOS
+
+- [x] Change FFT `data` variable to hold fft vals instead of time domain vals. 
+    - [x] Make FFT class `convolvable`  
+
+- [x] Make new convolution strategy to replace Brute Force (BF) strategy
+    - [ ] This is what the FFT is for -> So, instead, clean up the FFT to allow ease of multiplying 2 together at the length necessary to ensure convolution signal can be fully regained (ie FFT sizes are length that convolution would be (ie data1.length + data2.length - 1))
+    - [ ] Figure out how to deal with ifft having very small imaginary values
+    - [ ] Figure out how to deal with zero padding from Radix2 strategy in outputs. 
+- [x] Add Gaussian distribution generator
+- [x] Add White noise generator
+- [ ] Add `analytic signal` function 
+    - [ ] In where? Used for bandpass signals
+- [x] Implement Analog signal, allow custom sampling rate, include nyquist freq (based on freq or bandwidth)
+    - [x] Implement quantization
+    - [x] Allow companding techniques in quantization process
+- [x] Look at using another graphical tool instead of gruff 
+    - [x] Sticking with gruff for now
+- [x] Ensure FFT class is fully disconnected from FFTStrategy (ie having to zero fill for Radix2)
+    - [x] Ensure custom strategies can be built and implemented
+- [x] Same as above for companding strategies and window strategies
+- [x] Add #process to digitalSignal and FFT which accepts a block and iterates through each item in `data`
+    - [ ] Consider making this a module `processible`, which `requires_data`
+- [x] Look into using a better quantization strategy (right now mapping max <-> min to -1 <-> +1 is limited by floating point accuraccy for large ranges)
+- [x] Move `#map_to` method from AnalogSignal and put it in Digiproc::Functions
+- [ ] Investigate changing Digiproc to Digiproc
+- [ ] Review / cleanup / investigate proper scope of methods in Digiproc::Functions
+    - [ ] If the `#map_to` methods are kept in Digiproc::Functions, remove them from `AnalogSignal`
+- [x] Add IFFT, make FFT class multiplyable, 
+- [ ] test to ensure FFT multiplication equals convolution in the time domain
+- [ ] Add Decorators to Digiproc::QuickPlot to allow vertical lines being added, and maybe color schemes too 
+- [x] Make an alternative name to `dot` which is confusing because it is not the dot product (it does not sum at the end)
+- [ ] Add power_spectral_density, etc to FourierTransformable module
+- [ ] Figure out why convolution in matlab runs as slow as here, but cross-correlation is much faster
+- [ ] Make transition_width and bandwidth, wo, etc consistant in the filter part of the project (ie right now some are normalized frequencies from 0 to 1 and some are normalized frequencies in radians)
+- [ ] Ensure both Arrays and Vectors are supported for all `data` properties
+- [ ] Digiproc::AnalogSignal => sample_rate: to sample_interval:
+- [ ] Fix Strategies for DifferentialEncoding (change psk to the one in the PSK class)
+    - [ ] Find out why the PSK Differnetial Decoder starts at time != 0 
+    - [ ] Put the differential decoder in a different class (ie a decoding strategy)
+- [ ] Make `DigitalFilter` convolvable, fourierTranformable 
+- [ ] Automatcially match FFT sizes
+- [ ] Fix sizing for ifft, etc
+- [ ] Fix output of ifft (ie small imaginary numbers)
+- [ ] Ensure no freq domain inconsitancies with radian vs cycle^-1 units (ie inputs in functions should all be radians or hz)
+- [ ] "INTRODUCTION TO DIGITAL COMMUNICATIONS, ZIEMER, PETERSON" -> pg 239 Figure 4-16 add reciever for M-ary PSK decoding
+- [ ] Put DigitalSignal and Windows (and any other necessary functions not currently in) inside of the Digiproc Module
+- [ ] Make a working VerticalLine (and other) plot decorators
+- [ ] Make graphs visually similar across all APIs, make easier custom visual changes
+

data/bin/console

@@ -0,0 +1,14 @@
+#!/usr/bin/env ruby
+
+require "bundler/setup"
+require "digiproc"
+
+# 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(__FILE__)

data/bin/setup

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

data/config/environment.rb

@@ -0,0 +1,118 @@
+
+require "bundler/setup"
+
+Bundler.require(:default)
+
+require 'matrix'
+require 'ostruct'
+#Namespace
+require 'digiproc'
+
+
+#Extensions
+require './lib/extensions/core_extensions'
+
+Array.include Digiproc::CoreExtensions::ArrayExtension::DotProduct
+Array.include Digiproc::CoreExtensions::ArrayExtension::Sum 
+Array.include Digiproc::CoreExtensions::ArrayExtension::Multiply
+Math.extend Digiproc::CoreExtensions::MathExtension::Decible
+Vector.include Digiproc::CoreExtensions::VectorExtension::Projection::InstanceMethods
+Vector.extend Digiproc::CoreExtensions::VectorExtension::Projection::ClassMethods
+# Float.include Digiproc::CoreExtensions::FloatExtension::OddPatch
+class Gruff::Base
+
+    def draw_label(x_offset, index)
+        return if @hide_line_markers
+        if !@labels[index].nil? && @labels_seen[index].nil?
+          y_offset = @graph_bottom + LABEL_MARGIN
+          y_offset += @label_stagger_height if ((index.respond_to? :odd?) && index.odd?)
+          label_text = labels[index].to_s
+          if label_text.size > @label_max_size
+            if @label_truncation_style == :trailing_dots
+              if @label_max_size > 3
+                # 4 because '...' takes up 3 chars
+                label_text = "#{label_text[0 .. (@label_max_size - 4)]}..."
+              end
+            else # @label_truncation_style is :absolute (default)
+              label_text = label_text[0 .. (@label_max_size - 1)]
+            end
+  
+          end
+  
+          if x_offset >= @graph_left && x_offset <= @graph_right
+            @d.fill = @font_color
+            @d.font = @font if @font
+            @d.stroke('transparent')
+            @d.font_weight = NormalWeight
+            @d.pointsize = scale_fontsize(@marker_font_size)
+            @d.gravity = NorthGravity
+            @d = @d.annotate_scaled(@base_image,
+                                    1.0, 1.0,
+                                    x_offset, y_offset,
+                                    label_text, @scale)
+          end
+          @labels_seen[index] = 1
+          debug { @d.line 0.0, y_offset, @raw_columns, y_offset }
+        end
+    end
+
+end
+
+#Strategies
+require './lib/strategies/strategies'
+require './lib/strategies/fft/inverse_fft_conjugate_strategy'
+require './lib/strategies/fft/brute_force_dft_strategy'
+require './lib/strategies/fft/radix2_strategy'
+require './lib/strategies/convolution/bf_conv'
+require './lib/strategies/window/window'
+require './lib/strategies/window/blackman_window'
+require './lib/strategies/window/hamming_window'
+require './lib/strategies/window/hanning_window'
+require './lib/strategies/window/rectangular_window'
+require './lib/strategies/companding/custom_companding_strategy'
+require './lib/strategies/gaussian/gaussian_generator'
+require './lib/strategies/orthogonalize/gram_schmidt'
+require './lib/strategies/code/gray_code'
+require './lib/strategies/code/differential_encoding_strategy'
+require './lib/strategies/code/xor_differential_encoding_strategy'
+require './lib/strategies/code/xor_differential_encoding_zero_angle_strategy'
+require './lib/strategies/modulation/phase_shift_keying_strategy'
+require './lib/strategies/linear_algebra/jacobi_strategy'
+require './lib/strategies/linear_algebra/gauss_seidel_strategy'
+require './lib/strategies/linear_algebra/sor_strategy'
+require './lib/strategies/linear_algebra/sor2_strategy'
+
+#Modules
+require './lib/concerns/plottable'
+require './lib/concerns/requires_data'
+require './lib/concerns/initializable'
+require './lib/concerns/convolvable'
+require './lib/concerns/fourier_transformable'
+require './lib/concerns/multipliable'
+require './lib/concerns/data_properties'
+require './lib/functions'
+require './lib/probability/probability'
+require './lib/concerns/os'
+
+
+
+#Classes
+require './lib/fft'
+require './lib/signals/digital_signal'
+require './lib/filters/digital_filter.rb'
+require './lib/filters/lowpass_filter.rb'
+require './lib/filters/highpass_filter.rb'
+require './lib/filters/bandpass_filter.rb'
+require './lib/filters/bandstop_filter.rb'
+require './lib/signals/analog_signal'
+require './lib/probability/gaussian_distribution'
+require './lib/probability/theoretical_gaussian_distribution'
+require './lib/probability/bit_generator'
+require './lib/probability/binomial_distribution'
+require './lib/quick_plot'
+require './lib/rbplot'
+
+#Factories
+require './lib/factories/factories'
+require './lib/factories/window_factory'
+require './lib/factories/filter_factory'

data/console_tests.rb

@@ -0,0 +1,44 @@
+# x1 = Digiproc::AnalogSignal.new(eqn: ->(t){ 8 * Math.cos(20 * Math::PI * t) }, sample_rate: 1.0 / 51.0, size: 6, quantization_bits: 4, quant_max: 8, quant_min: -8)
+# x2 = Digiproc::AnalogSignal.new(eqn: ->(t){  Math.cos(20 * Math::PI * t) }, sample_rate: 1.0 / 51.0, size: 6, quantization_bits: 4, quant_max: 8, quant_min: -8)
+
+# compression_eqn = ->(n){ n < 0 ? (-4 * ((-n) ** (1.0 / 3))) : (4 * n ** (1.0 / 3))}
+# expansion_eqn =->(n){ n < 0 ? ((-1.0 / 64) * ((-n) ** 3)) : ((1.0 / 64) * n ** 3)}
+# compander = CustomCompandingStrategy.new(compression_eqn, expansion_eqn)
+
+# x3 = Digiproc::AnalogSignal.new(eqn: ->(t){ 8 * Math.cos(20 * Math::PI * t) }, sample_rate: 1.0 / 51.0, size: 6, quantization_bits: 4, quant_max: 8, quant_min: -8, companding_strategy: compander)
+# x4 = Digiproc::AnalogSignal.new(eqn: ->(t){  Math.cos(20 * Math::PI * t) }, sample_rate: 1.0 / 51.0, size: 6, quantization_bits: 4, quant_max: 8, quant_min: -8, companding_strategy: compander)
+
+# puts "x(1): "
+# puts x1.digitize.to_s
+# puts "Raw samples: #{x1.raw_samples}"
+# puts "Normalized Quantization RMS Error: #{x1.normalized_quantization_rms_error.round(3)}"
+# puts "\n\nx(2):"
+# puts x2.digitize.to_s
+# puts "Raw samples: #{x2.raw_samples}"
+# puts "Normalized Quantization RMS Error: #{x2.normalized_quantization_rms_error.round(3)}"
+# puts "\n\nx(1) with companding: "
+# puts x3.digitize.to_s
+# puts "Raw samples: #{x3.raw_samples}"
+# puts "Normalized Quantization RMS Error: #{x3.normalized_quantization_rms_error.round(3)}"
+# puts "\n\nx(2) with companding:"
+# puts x4.digitize.to_s
+# puts "Raw samples: #{x4.raw_samples}"
+# puts "Normalized Quantization RMS Error: #{x4.normalized_quantization_rms_error.round(3)}"
+# binding.pry
+
+# signal = ->(t){Math.cos(2 * Math::PI * t)} #cos at 1 Hz
+# sample_times = (0..1000).map{ |int| int / 100.0} #Sample for 10 seconds every 0.01 seconds (100 times per period over 10 periods)
+# data = sample_times.map{ |time| signal.call(time) }
+# ft = []
+# for k in 0...data.length do
+#     tot = 0
+#     data.each_with_index do |x_n, n|
+#         tot += x_n * Math::E ** (Complex(0,-1) * 2.0 * Math::PI * k * n / data.length.to_f)
+#     end
+#     ft << tot
+# end
+
+# fft = Digiproc::FFT.new(data: data).calculate
+
+
+binding.pry

data/digiproc.gemspec

@@ -0,0 +1,49 @@
+
+lib = File.expand_path("../lib", __FILE__)
+$LOAD_PATH.unshift(lib) unless $LOAD_PATH.include?(lib)
+require "digiproc/version"
+
+
+Gem::Specification.new do |spec|
+  spec.name          = "digiproc"
+  spec.version       = Digiproc::VERSION
+  spec.authors       = ["micahshute"]
+  spec.email         = ["micah.shute@gmail.com"]
+
+  spec.summary       = %q{Perform basic Digital Signal Processing tasks, including convolution, fft, filtering.}
+  spec.description   = %q{Allows design of digital signals using the FFT, design of Digital Filters using the Windowing Method, creation of Digital Signals or Analog Signals sampled at a certain interval, convolution, cross-correlation, and visualization of the data. .}
+  spec.homepage      = "https://rubygems.org/gems/digiproc"
+  spec.license       = "MIT"
+
+  # Prevent pushing this gem to RubyGems.org. To allow pushes either set the 'allowed_push_host'
+  # to allow pushing to a single host or delete this section to allow pushing to any host.
+  if spec.respond_to?(:metadata)
+    spec.metadata["allowed_push_host"] = "https://rubygems.org"
+
+    spec.metadata["homepage_uri"] = spec.homepage
+    spec.metadata["source_code_uri"] = "https://github.com/micahshute/digiproc"
+    spec.metadata["changelog_uri"] = "https://github.com/micahshute/digiproc"
+  else
+    raise "RubyGems 2.0 or newer is required to protect against " \
+      "public gem pushes."
+  end
+
+  # Specify which files should be added to the gem when it is released.
+  # The `git ls-files -z` loads the files in the RubyGem that have been added into git.
+  spec.files         = Dir.chdir(File.expand_path('..', __FILE__)) do
+    `git ls-files -z`.split("\x0").reject { |f| f.match(%r{^(test|spec|features|png)/}) }
+  end
+  spec.bindir        = "exe"
+  spec.executables   = spec.files.grep(%r{^exe/}) { |f| File.basename(f) }
+  spec.require_paths = ["lib"]
+
+  spec.add_development_dependency "bundler"
+  spec.add_development_dependency "rake", "~> 10.0"
+  spec.add_development_dependency "rspec", "~> 3.0"
+  spec.add_development_dependency 'rdoc'
+  spec.add_development_dependency 'pry'
+  spec.add_development_dependency 'gruff', "~> 0.7.0" 
+  spec.add_runtime_dependency "gruff", "~> 0.7.0"
+  # spec.add_dependency "gruff"
+end
+

data/examples/analog_signals/analog_to_digital.rb

@@ -0,0 +1,31 @@
+# Sample an analog signal, turn into digital signal and calculate the FFT
+plt = Digiproc::QuickPlot
+analog_signal_eqn = ->(t){ (t >= 0 and t <= 0.1) ? (100 * t * Math.cos(Math::PI*2.5*t)) : 0 }
+
+analog_signal_eqn = ->(t){ Math.cos(2*Math::PI * 50 * t) + Math.cos(2 * Math::PI * 1250 * t) }
+
+analog_signal = Digiproc::AnalogSignal.new(eqn: analog_signal_eqn, sample_rate: 0.001, size: 1000)
+
+digital_signal = analog_signal.to_ds
+
+fft = digital_signal.fft
+
+fft.plot_magnitude(path: './examples/analog_signals/')
+
+rate = 0.01
+size = 1000
+fn = ->(x){ Math.sin(x) }
+as64 = Digiproc::AnalogSignal.new(eqn: fn, sample_rate: rate, size: size, quantization_bits: 64, quant_max: 1, quant_min: -1)
+as8 = Digiproc::AnalogSignal.new(eqn: fn, sample_rate: rate, size: size, quantization_bits: 8, quant_max: 1, quant_min: -1)
+as4 = Digiproc::AnalogSignal.new(eqn: fn, sample_rate: rate, size: size, quantization_bits: 4, quant_max: 1, quant_min: -1)
+t = Digiproc::Functions.linspace(0, rate * size, size)
+plt = Digiproc::Rbplot.line(t, as64.digitize, "64 bits")
+plt.add_line(t, as8.digitize, '8 bits')
+plt.add_line(t, as4.digitize, '4 bits')
+plt.path('./examples/analog_signals/')
+plt.title('Quantization Outputs')
+plt.xlabel('time (s)')
+plt.ylabel('magnitude')
+plt.xsteps(10)
+plt.theme(:dark)
+plt.show