paddlec 0.0.1

data/samples/stereo_chunk.rb

@@ -0,0 +1,161 @@
+#!/usr/bin/env ruby
+
+require 'roctave'
+require '../lib/paddlec'
+
+def show_bar(f, t = 80.0)
+	p = (f.pos * t / f.size).round
+	STDOUT.write "\r[#{'='*p}>#{'-'*(t - p)}]  "
+	STDOUT.flush
+end
+
+
+show_filters = false
+
+fs = 2.048e6
+
+# b = File.binread(File.expand_path '~/sdrrec/pouet.raw').unpack('C*').each_slice(2).collect{|a| Complex((a.first - 128) / 128.0, (a.last - 128) / 128.0)}.to_complex_buffer
+# puts b.length
+# File.binwrite('pouet_cf32le.raw', b.pack)
+# exit
+
+
+unless File.file? 'coefficients.f32le' then
+	STDOUT.write 'Generating filter coeficients... '; STDOUT.flush
+	now = Time.now
+
+	b_bbdecimator = Roctave.fir1(254, 95e3*2/fs)
+	Roctave.freqz(b_bbdecimator, :magnitude, nb_points: 2048, fs: fs) if show_filters
+
+	b_differentiator = Roctave.fir_differentiator(32)
+	Roctave.freqz(b_differentiator, :magnitude, nb_points: 2048, fs: fs/8) if show_filters
+
+	b_adecimator = Roctave.fir_low_pass(512, 15e3*2/256e3)
+	Roctave.freqz(b_adecimator, :magnitude, nb_points: 2048, fs: fs/8) if show_filters
+
+	tau = 50e-6
+	ts = 8*8/fs
+	b = [1 - Math.exp(-ts/tau)]
+	a = [1, -Math.exp(-ts/tau)]
+	h, w = Roctave.freqz(b, a, nb_points: 1024)
+	h.collect!{|v| v.abs}
+	w.collect!{|v| v/w.last}
+	n = 2048
+	b_deemphasis = Roctave.fir2(n, w, h, :odd_symmetry)
+	Roctave.freqz(b_deemphasis, :magnitude, nb_points: 2048, fs: fs/8/8) if show_filters
+
+	b_bppilote = Roctave.fir_band_pass(256, 2*19e3/(fs/8), 2*1e3/(fs/8))
+	Roctave.freqz(b_bppilote, :magnitude, nb_points: 2048, fs: fs/8) if show_filters
+
+	b_hilbert = Roctave.fir_hilbert(48, window: :blackman)
+	Roctave.freqz(b_hilbert, :magnitude, nb_points: 2048, fs: fs/8) if show_filters
+
+	b_bplmr = Roctave.fir_band_pass(256, 2*38e3/(fs/8), 2*30e3/(fs/8), window: :blackman)
+	Roctave.freqz(b_bplmr, :magnitude, nb_points: 2048, fs: fs/8) if show_filters
+
+	puts (Time.now - now).round(3)
+
+	File.open('coefficients.f32le', 'wb') do |ofile|
+		ofile.write b_bbdecimator.pack('f*')
+		ofile.write b_differentiator.pack('f*')
+		ofile.write b_adecimator.pack('f*')
+		ofile.write b_deemphasis.pack('f*')
+		ofile.write b_bppilote.pack('f*')
+		ofile.write b_hilbert.pack('f*')
+		ofile.write b_bplmr.pack('f*')
+	end
+else
+	ifile = File.open('coefficients.f32le', 'rb')
+	b = ifile.read.unpack('f*')
+	ifile.close
+	b_bbdecimator = b.shift(255)
+	b_differentiator = b.shift(33)
+	b_adecimator = b.shift(513)
+	b_deemphasis = b.shift(2049)
+	b_bppilote = b.shift(257)
+	b_hilbert = b.shift(49)
+	b_bplmr = b.shift(257)
+end
+
+
+
+bbdecimator = PaddleC::FirDecimator.new b_bbdecimator, 8
+differentiator = PaddleC::FirFilter.new b_differentiator
+adecimator_lpr = PaddleC::FirDecimator.new b_adecimator, 8
+adecimator_lmr = PaddleC::FirDecimator.new b_adecimator, 8
+deemphasis_lpr = PaddleC::FirFilter.new b_deemphasis
+deemphasis_lmr = PaddleC::FirFilter.new b_deemphasis
+bppilote       = PaddleC::FirFilter.new b_bppilote
+bplmr          = PaddleC::FirFilter.new b_bplmr
+hilbert        = PaddleC::FirTransformer.new b_hilbert
+delay_p = PaddleC::Delay.new 152
+delay_h = PaddleC::Delay.new 24
+
+time_chunk = 100e-3
+len = (32e3 * time_chunk).round
+
+decoding = 0.0
+waiting = 0.0
+fname = '~/sdrrec/102.4MHz_2.048Msps_cu8.raw'
+fname = '~/sdrrec/102.4MHz_2.048Msps_cu8_1.raw'
+fname = '~/sdrrec/102.4MHz_2.048Msps_cu8_2.raw'
+fname = '~/sdrrec/102.4MHz_2.048Msps_cu8_3.raw'
+fname = '~/sdrrec/pouet.raw'
+fname = '~/sdrrec/95.0MHz_2.048Msps_cu8.raw'
+PaddleC::PulseAudio::Simple::Sink.open(stream_name: fname, sample_spec: 'float32 2ch 32000Hz') do |sink|
+	File.open(File.expand_path(fname), 'rb') do |f|
+		loop do
+			show_bar(f)
+
+			now = Time.now
+
+			if2MHz_complex = PaddleC::ComplexBuffer.unpack(f.read(8*8*len*2), :u8)
+
+			baseband256kHz_complex = bbdecimator.decimate(if2MHz_complex)
+
+			baseband_diff_complex, baseband_delayed_complex = differentiator.filter(baseband256kHz_complex, delayed: true)
+
+			demodulated = (baseband_diff_complex.imag*baseband_delayed_complex.real - baseband_diff_complex.real*baseband_delayed_complex.imag) / baseband_delayed_complex.abs2
+
+			pilote19 = bppilote.filter demodulated
+			pilote19c = hilbert.transform pilote19
+			pilote_phase = pilote19c.phase * 2
+			pilote38 = pilote_phase.sin
+
+			lmrh = bplmr.filter demodulated
+			lmrhd = delay_h.delay lmrh
+			lmrm = lmrhd * pilote38
+			lmremph = adecimator_lmr.decimate(lmrm)
+			lmr = deemphasis_lmr.filter(lmremph) * 0.25
+
+			lpremphdelayed = delay_p.delay demodulated
+			lpremph = adecimator_lpr.decimate(lpremphdelayed)
+			lpr = deemphasis_lpr.filter(lpremph) * 0.25
+
+			left = (lpr + lmr) / 2
+			right = (lpr - lmr) / 2
+
+			stereo = PaddleC::ComplexBuffer.new(left, right)
+
+			decoding += Time.now - now
+
+			now = Time.now
+
+			sink << stereo
+
+			waiting += Time.now - now
+
+			break if f.eof?
+		end
+		show_bar(f)
+		puts
+	end
+end
+
+
+puts "decoding took #{decoding.round(3)}s"
+puts "waiting took #{waiting.round(3)}s"
+
+
+
+

metadata

@@ -0,0 +1,99 @@
+--- !ruby/object:Gem::Specification
+name: paddlec
+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: roctave
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '0.0'
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: 0.0.1
+  type: :runtime
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '0.0'
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: 0.0.1
+description: PaddleC is a Ruby C extension attempting to provide objects and methods
+  to rapidly set up real-time signal processing test benches in Ruby.
+email: theotime.bollengier@gmail.com
+executables: []
+extensions:
+- ext/paddlec/extconf.rb
+extra_rdoc_files: []
+files:
+- LICENSE
+- README.md
+- ext/libpaddlec/arithmetic.c
+- ext/libpaddlec/comparison.c
+- ext/libpaddlec/complex.c
+- ext/libpaddlec/delay.c
+- ext/libpaddlec/fir_filter.c
+- ext/libpaddlec/fir_filter_avx.c
+- ext/libpaddlec/fir_filter_neon.c
+- ext/libpaddlec/fir_filter_sse.c
+- ext/libpaddlec/libpaddlec.c
+- ext/libpaddlec/libpaddlec.h
+- ext/libpaddlec/math.c
+- ext/libpaddlec/no_fast_math.c
+- ext/libpaddlec/rounding.c
+- ext/paddlec/complex_buffer.c
+- ext/paddlec/complex_buffer.h
+- ext/paddlec/delay.c
+- ext/paddlec/delay.h
+- ext/paddlec/extconf.rb
+- ext/paddlec/fir_filter.c
+- ext/paddlec/fir_filter.h
+- ext/paddlec/float_buffer.c
+- ext/paddlec/float_buffer.h
+- ext/paddlec/paddlec.c
+- ext/paddlec/paddlec.h
+- ext/paddlec/pulseaudio.c
+- ext/paddlec/pulseaudio.h
+- lib/paddlec.rb
+- lib/paddlec/version.rb
+- paddlec.gemspec
+- samples/fmdemod.rb
+- samples/fmdemod_chunk.rb
+- samples/fmdemod_chunk_buffer.rb
+- samples/stereo_chunk.rb
+homepage: https://gitlab.com/theotime_bollengier/paddlec
+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: A Ruby C extension for real-time signal processing.
+test_files: []