roctave 0.0.1

data/samples/firpm.rb

@@ -0,0 +1,78 @@
+#!/usr/bin/env ruby
+
+require_relative '../lib/roctave.rb'
+
+## Low pass type I
+n = 60 
+f = [(0 .. 0.2), (0.25 .. 0.5), (0.55 .. 1)]
+a = [0, 1, 0]
+b = Roctave.firpm(n, f, a, :even_symmetry, 16)
+exit
+
+# ## Low pass type I
+# n = 32 
+# f = [(0 .. 0.45), (0.55 .. 0.7), (0.75 .. 1)]
+# a = [1, 0, 0.5]
+# w = [1, 1, 20]
+# b = Roctave.firpm(n, f, a, w, :even_symmetry)
+
+# ## Low pass type II
+# n = 11 
+# f = [(0 .. 0.3), (0.7 .. 1)]
+# a = [1, 0]
+# w = [2, 1]
+# b = Roctave.firls(n, f, a, w, :even_symmetry)
+
+# ## Differentiator type III
+# n = 10 
+# f = [0.05, 0.75]
+# a = [f]
+# w = [Roctave.linspace(f.first, f.last, 16).collect{|v| 1/v}]
+# b = Roctave.firls(n, f, a, w, 16, :odd_symmetry)
+
+# ## Hilbert transformer type IV
+# n = 11 
+# f = [1.0/8, 7.0/8]
+# a = [1]
+# w = [1]
+# b = Roctave.firls(n, f, a, w, :odd_symmetry)
+
+
+
+r = Roctave.dft(b, 2048)
+h = r[0..1024].abs
+p = r[0..1024].arg
+p = Roctave.unwrap(p)
+w = Roctave.linspace(0, 1, 1025)
+
+# if (n & 1) == 0 then
+# 	Roctave.stem((-n/2..n/2).to_a, b, :filled, grid: true)
+# else
+# 	Roctave.stem((n+1).times.collect{|i| i - (n+1)/2.0 + 0.5}, b, :filled, grid: true)
+# end
+# Roctave.plot(w, p, grid: true, title: 'Phase')
+# Roctave.plot(w, h.collect{|v| 10*Math.log10(v)}, xlabel: 'Normalized frequency', ylabel: 'Magnitude response (dB)', grid: true, title: 'Magnitude')
+Roctave.plot([0, 0.3], [1, 1], '-b', [0.7, 1], [0, 0], '-b', w, h, 'r', title: 'Magnitude response', xlim: (0 .. 1), grid: true)
+
+
+
+
+
+exit
+n = 4
+f = [(0 .. 0.4), (0.6 .. 1)]
+a = [1, 0]
+b = Roctave.firls(n, f, a, 1)
+exit
+
+n = 32
+grid = 8
+f = [[0.2, 0.3], [0.4, 0.6], [0.7, 0.8]]
+a = [
+	0.5,
+	Roctave.linspace(0, Math::PI, 64).collect{|t| 1 - Math.sin(t)},
+	(0.5 .. 0)
+]
+
+b = Roctave.firls(n, f, a, grid)
+

data/samples/hilbert_transformer.rb

@@ -0,0 +1,20 @@
+#!/usr/bin/env ruby
+
+require_relative '../lib/roctave.rb'
+
+# H(e^(jw)) = -j when 0 <= w < pi, j when -pi <= w < 0
+# h(n) = 0 when n == 0, else (1 - (-1)^n) / (pi.n)
+
+n = 64
+
+b = Roctave.fir_hilbert(n, window: :nuttall)
+
+Roctave.stem((-n/2..n/2).to_a, b, :filled, grid: true, title: "FIR type III order #{n}")
+
+r = Roctave.dft(b, 2048)
+r = Roctave.fftshift(r)
+h = r.abs
+p = r.arg
+w = Roctave.linspace(-1, 1, 2048)
+Roctave.plot([-1, 1], [1, 1], '--;Ideal;', w, h, ";FIR type III order #{n};", grid: true, title: 'Magnitude')
+Roctave.plot(w, p, grid: true, title: 'Phase')

data/samples/hilbert_transformer_frequency_scaling.rb

@@ -0,0 +1,47 @@
+#!/usr/bin/env ruby
+
+require_relative '../lib/roctave.rb'
+
+
+
+n = 33
+f = [0, 1]
+m = [Complex(0, -1), Complex(0, -1)]
+
+b = Roctave.fir2(n, f, m, :odd_symetry, window: :blackman)
+
+Roctave.stem(b, :filled, grid: true)
+
+r = Roctave.dft(b, 2048)
+r = Roctave.fftshift(r)
+h = r.abs
+p = r.arg
+w = Roctave.linspace(-1, 1, 2048)
+Roctave.plot([-1, 1], [1, 1], '--;Ideal;', w, h, ';Actual;', grid: true, title: 'Magnitude')
+Roctave.plot(w, p, grid: true, title: 'Phase')
+
+=begin
+n = [2, order.to_i].max
+n += n & 1
+grid = 1024
+mgridl = Roctave.linspace(0, Math::PI, grid + 1).collect{|i| Complex(0.0, i)}
+mgridr = mgridl[1...-1].reverse.conj
+mgrid = mgridl + mgridr
+b = Roctave.idft(mgrid)
+mid = (n+1) / 2.0
+b = (b[-(mid.floor) .. -1] + b[0 ... mid.ceil]).collect{|e| e.real}
+window = Roctave.blackman(n+1)
+b = b.collect.with_index{|e, i| e*window[i]}
+
+indexes = (-n/2..n/2).to_a
+Roctave.stem(indexes, b, :filled, grid: true, title: "FIR type III order #{n}")
+
+r = Roctave.dft(b, 2048)
+h = r[0...1024].abs
+p = r[0...1024].arg
+w = Roctave.linspace(0, 1, 1024)
+Roctave.plot([0, Math::PI], '--;Ideal;', w, h, ";FIR type III order #{n};", grid: true, title: 'Magnitude')
+Roctave.plot(w, p, grid: true, title: 'Phase')
+=end
+
+

data/samples/plot.rb

@@ -0,0 +1,45 @@
+#!/usr/bin/env ruby
+
+require_relative '../lib/roctave.rb'
+
+# t = Roctave.linspace(0, 1)
+# s1 = t.collect{|e| Math.sin(2*Math::PI*2*e + 0*Math::PI/2)}
+# s2 = t.collect{|e| Math.sin(2*Math::PI*2*e + 1*Math::PI/2)}
+# s3 = t.collect{|e| Math.sin(2*Math::PI*2*e + 2*Math::PI/2)}
+# s4 = t.collect{|e| Math.sin(2*Math::PI*2*e + 3*Math::PI/2)}
+# Roctave.plot(t, s1, '-+;{/Symbol f} = 0;', t, s2, '--o;{/Symbol f} = {/Symbol p} / 2;', t, s3, ':*;{/Symbol f} = 2 {/Symbol p} / 2;', t, s4, '-..;{/Symbol f} = 3 {/Symbol p} / 2;', grid: true, title: 'y = sin(x + {/Symbol f})', xlabel: 'Time t^2 h^{titi} t_2 t_{titi} Z@^{en haut}_{en bas}  escape\\\_underscore', ylabel: 'Amplitude', terminal: 'pngcairo transparent truecolor rounded enhanced font "Sans,10" size 1280,720', output: 'titi.png')
+# #Roctave.plot(t, s1, '+-;{/Symbol f} = 0;', t, s2, '--o;{/Symbol f} = {/Symbol p} / 2;', t, s3, ':*;{/Symbol f} = 2 {/Symbol p} / 2;', t, s4, '-..;{/Symbol f} = 3 {/Symbol p} / 2;', grid: true, title: 'y = sin(x + {/Symbol f})', xlabel: 'Time t^2 h^{titi} t_2 t_{titi} Z@^{en haut}_{en bas}  escape\\\_underscore', ylabel: 'Amplitude', terminal: "wxt enhanced persist font \"Times New Roman,13\" size 1280,720")
+# Roctave.plot(t, s1, '+-;{/Symbol f} = 0;', t, s2, '--o;{/Symbol f} = {/Symbol p} / 2;', t, s3, ':*;{/Symbol f} = 2 {/Symbol p} / 2;', t, s4, '-..;{/Symbol f} = 3 {/Symbol p} / 2;', grid: true, title: 'y = sin(x + {/Symbol f})', xlabel: 'Time t^2 h^{titi} t_2 t_{titi} Z@^{en haut}_{en bas}  escape\\\_underscore', ylabel: 'Amplitude', terminal: "wxt enhanced persist font \"Sans,13\" size 1280,720")
+# Roctave.plot(
+# 	t, s1, 
+# 	:color,
+# 	'orange',
+# 	:dashtype,
+# 	'-_',
+# 	:pointtype,
+# 	:square,
+# 	:empty,
+# 	:pointsize,
+# 	2,
+# 	:linewidth,
+# 	2,
+# 	grid: true, title: 'y = sin(x + {/Symbol f})', xlabel: 'Time t^2 h^{titi} t_2 t_{titi} Z@^{en haut}_{en bas}  escape\\\_underscore', ylabel: 'Amplitude', terminal: "wxt enhanced persist size 1280,720"
+# )
+
+t = Roctave.linspace(0, 1, 1000)
+s1 = t.collect{|e| Math.sin(2*Math::PI*2*e + 0*Math::PI/2)}
+s2 = t.collect{|e| Math.sin(2*Math::PI*2*e + 1*Math::PI/2)}
+s3 = t.collect{|e| Math.sin(2*Math::PI*2*e + 2*Math::PI/2)}
+s4 = t.collect{|e| Math.sin(2*Math::PI*2*e + 3*Math::PI/2)}
+Roctave.plot(
+	t, s1, ';{/Symbol f} = 0;', 
+	t, s2, ';{/Symbol f} = {/Symbol p} / 2;', 
+	t, s3, ';{/Symbol f} = 2 {/Symbol p} / 2;', 
+	t, s4, ';{/Symbol f} = 3 {/Symbol p} / 2;', 
+	grid: true, 
+	title: 'y = sin(x + {/Symbol f})', 
+	xlabel: 'Time t^2 h^{titi} t_2 t_{titi} Z@^{en haut}_{en bas}  escape\\\_underscore', 
+	ylabel: 'Amplitude', 
+	terminal: "wxt enhanced persist size 854,480"
+)
+

data/samples/stem.rb

@@ -0,0 +1,8 @@
+#!/usr/bin/env ruby
+
+require_relative '../lib/roctave.rb'
+
+t = Roctave.linspace(0, 2*Math::PI, 50)
+y = t.collect{|x| Math.exp(0.3*x) * Math.sin(3*x)}
+Roctave.stem(t, y, :legend, 'Coucou', t, y.reverse, ';plop;', :filled, grid: true)
+

data/samples/type1.rb

@@ -0,0 +1,25 @@
+#!/usr/bin/env ruby
+
+require_relative '../lib/roctave.rb'
+
+n = 32
+fc = 0.5
+
+b1 = Roctave.fir_low_pass(n, fc)
+ns = (-n/2..n/2).to_a
+Roctave.stem(ns, b1, :filled, grid: true)
+
+h1 = Roctave.dft(b1, 2048).abs[0...1024]
+w = Roctave.linspace(0.0, 1.0, 1024)
+
+bb = b1[n/2+1..-1]
+h2 = (0...1024).collect do |i|
+	acc = 0.0
+	bb.each.with_index do |e, j|
+		acc += e*Math.cos((i*Math::PI/1024)*(j+1))
+	end
+	b1[n/2] + 2.0*acc
+end
+
+Roctave.plot(w, h1, ';real;', w, h2, ';test;', grid: true, ylim: [0, 1.4])
+

data/samples/type3.rb

@@ -0,0 +1,24 @@
+#!/usr/bin/env ruby
+
+require_relative '../lib/roctave.rb'
+
+n = 32
+
+b1 = Roctave.fir_differentiator(n, window: :blackman)
+ns = (-n/2..n/2).to_a
+Roctave.stem(ns, b1, :filled, grid: true)
+
+h1 = Roctave.dft(b1, 2048).abs[0...1024]
+w = Roctave.linspace(0.0, 1.0, 1024)
+
+bb = b1[n/2+1..-1]
+h2 = (0...1024).collect do |i|
+	acc = 0.0
+	bb.each.with_index do |e, j|
+		acc += -e*Math.sin((i*Math::PI/1024)*(j+1))
+	end
+	b1[n/2] + 2.0*acc
+end
+
+Roctave.plot(w, h1, ';real;', w, h2, ';test;', grid: true)
+

data/samples/windows.rb

@@ -0,0 +1,25 @@
+#!/usr/bin/env ruby
+
+require_relative '../lib/roctave.rb'
+require 'gnuplot'
+
+Gnuplot.open do |gp|
+	Gnuplot::Plot.new(gp) do |plot|
+		plot.terminal 'wxt'
+		plot.title "Roctave::Window"
+		plot.yrange "[0:1.1]"
+		plot.grid
+
+		n = 100
+		x = (0...n).to_a
+
+		[:hann, :hamming, :blackman, :nuttall, :blackman_nuttall, :blackman_harris, :gaussian].each do |name|
+			y = Roctave::Window.send(name, n)
+			plot.data << Gnuplot::DataSet.new([x, y]) { |ds|
+				ds.with = "lines"
+				ds.title = name.to_s.split('_').collect{|s| s.capitalize}.join(' ')
+			}
+		end	
+	end
+end
+

metadata

@@ -0,0 +1,123 @@
+--- !ruby/object:Gem::Specification
+name: roctave
+version: !ruby/object:Gem::Version
+  version: 0.0.1
+platform: ruby
+authors:
+- Théotime Bollengier
+autorequire: 
+bindir: bin
+cert_chain: []
+date: 2020-01-22 00:00:00.000000000 Z
+dependencies:
+- !ruby/object:Gem::Dependency
+  name: gnuplot
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '2.6'
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: 2.6.2
+  type: :runtime
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '2.6'
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: 2.6.2
+description: Roctave is an attempt to provide some Matlab / Gnu Octave functions in
+  a Ruby gem to generate and analyze digital filters.
+email: theotime.bollengier@gmail.com
+executables: []
+extensions:
+- ext/roctave/extconf.rb
+extra_rdoc_files: []
+files:
+- LICENSE
+- README.md
+- ext/roctave/cu8_file_reader.c
+- ext/roctave/cu8_file_reader.h
+- ext/roctave/extconf.rb
+- ext/roctave/fir_filter.c
+- ext/roctave/fir_filter.h
+- ext/roctave/freq_shifter.c
+- ext/roctave/freq_shifter.h
+- ext/roctave/iir_filter.c
+- ext/roctave/iir_filter.h
+- ext/roctave/roctave.c
+- ext/roctave/roctave.h
+- lib/roctave.rb
+- lib/roctave/bilinear.rb
+- lib/roctave/butter.rb
+- lib/roctave/cheby.rb
+- lib/roctave/cu8_file_reader.rb
+- lib/roctave/dft.rb
+- lib/roctave/filter.rb
+- lib/roctave/finite_difference_coefficients.rb
+- lib/roctave/fir.rb
+- lib/roctave/fir1.rb
+- lib/roctave/fir2.rb
+- lib/roctave/fir_design.rb
+- lib/roctave/firls.rb
+- lib/roctave/firpm.rb
+- lib/roctave/freq_shifter.rb
+- lib/roctave/freqz.rb
+- lib/roctave/iir.rb
+- lib/roctave/interp1.rb
+- lib/roctave/plot.rb
+- lib/roctave/poly.rb
+- lib/roctave/roots.rb
+- lib/roctave/sftrans.rb
+- lib/roctave/version.rb
+- lib/roctave/window.rb
+- roctave.gemspec
+- samples/butter.rb
+- samples/cheby.rb
+- samples/dft.rb
+- samples/differentiator.rb
+- samples/differentiator_frequency_scaling.rb
+- samples/fft.rb
+- samples/finite_difference_coefficient.rb
+- samples/fir1.rb
+- samples/fir2.rb
+- samples/fir2_windows.rb
+- samples/fir2bank.rb
+- samples/fir_low_pass.rb
+- samples/firls.rb
+- samples/firpm.rb
+- samples/hilbert_transformer.rb
+- samples/hilbert_transformer_frequency_scaling.rb
+- samples/plot.rb
+- samples/stem.rb
+- samples/type1.rb
+- samples/type3.rb
+- samples/windows.rb
+homepage: https://gitlab.com/theotime_bollengier/roctave
+licenses:
+- GPL-3.0+
+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: []
+rubygems_version: 3.0.6
+signing_key: 
+specification_version: 4
+summary: Some Matlab / Gnu Octave functions implemented in Ruby.
+test_files: []