radio 0.0.1

data/README.md

@@ -0,0 +1 @@
+Amateur radio software

data/lib/radio.rb

@@ -0,0 +1,14 @@
+%w{
+  radio/psk31/*.rb
+}.each do |glob|
+  Dir.glob(File.expand_path(glob, File.dirname(__FILE__))).each do |filename|
+    require filename
+  end
+end
+
+class Radio
+
+  PI = Math::PI.freeze
+  PI2 = (8.0 * Math.atan(1.0)).freeze
+  
+end

data/lib/radio/psk31/bit_detect.rb

@@ -0,0 +1,46 @@
+class Radio
+  
+  class PSK31
+    
+    class BitDetect
+      
+      AVG_SAMPLES = 50.freeze
+      CHANGE_DELAY = 5
+      
+      def initialize
+        @averages = Array.new 21
+        reset
+      end
+      
+      def reset
+        @averages.fill 0.0
+        @phase = 0
+        @peak = 0
+        @next_peak = 0
+        @change_at = 0
+      end
+      
+      def call sample_x, sample_y
+        yield if @phase == @peak
+        @peak = @next_peak if @phase == @change_at
+        energy = sample_x**2 + sample_y**2
+        @averages[@phase] = (1.0-1.0/AVG_SAMPLES)*@averages[@phase] + (1.0/AVG_SAMPLES)*energy
+        @phase += 1
+        if @phase > 15
+          @phase = 0
+          max = -1e10
+          for i in 0...16
+            energy = @averages[i]
+            if energy > max
+              @next_peak = i
+              @change_at = (i + CHANGE_DELAY) & 0x0F
+              max = energy
+            end
+          end
+        end
+      end
+      
+    end
+    
+  end
+end

data/lib/radio/psk31/decoder.rb

@@ -0,0 +1,65 @@
+class Radio
+  class PSK31
+    
+    class Decoder
+      
+      attr_accessor :mode # :bpsk or :qpsk or :qpsklsb
+      
+      def initialize
+        @mode = :bpsk 
+        reset
+      end
+      
+      def reset
+        @prev_i = 0.0
+        @prev_q = 0.0
+        @this_i = 0.0
+        @this_q = 0.0
+        @code = 0
+        @prev_bit = false
+      end
+      
+      def call new_i, new_q
+        @prev_i = @this_i
+        @prev_q = @this_q
+        @this_i = new_i
+        @this_q = new_q
+        vect_y = @prev_i*@this_i + @prev_q*@this_q
+        if @mode == :bpsk
+          bit = vect_y >= 0.0 ? 1 : 0
+        else
+          vect_x = @prev_i*@this_q - @this_i*@prev_q
+          if vect_y == 0.0 # atan2 errors on (0,0)
+            angle = PI
+          elsif @mode == :qpsklsb
+            angle = PI + Math.atan2(vect_y, -vect_x)
+          else # :qpsk or :bpsk
+            angle = PI + Math.atan2(vect_y, vect_x)
+          end
+          bit = viterbi angle
+        end
+        if bit==0 and @prev_bit==0
+          if @code != 0
+            @code >>= 2
+            @code &= 0x07FF
+            ch = VARICODE_DECODE_TABLE[@code]
+            yield ch if ch
+            @code = 0
+          end
+        else
+          @code <<= 1
+          @code |= bit
+          @prev_bit = bit
+        end
+      end
+      
+      private
+      
+      def viterbi
+        raise 'todo'
+      end
+      
+    end
+
+  end
+end

data/lib/radio/psk31/filters.rb

@@ -0,0 +1,220 @@
+class Radio
+  
+  class PSK31
+    
+    # An instance of Filters handles everything from 48kHz to 500Hz.
+    # PSK31/PSK63/PSK125 emit at 500Hz/1000Hz/2000Hz respectively.
+    
+    class Filters
+      
+      attr_accessor :frequency
+      attr_accessor :clock # 8000.0 +/-
+      attr_accessor :phase_inc
+      attr_accessor :speed # 31 | 63 | 125
+      
+      # Format of the input data stream is specified in the same
+      # format as String#unpack.  Here are some examples:
+      #   'C' 8-bit unsigned mono
+      #   'xxv' Little-endian 16-bit right channel
+      # Not everything is supported outside this native Ruby implementation.
+      # The following are guaranteed to be in a C or Java implementation.
+      #   x - ignores a byte (in an unused channel)
+      #   C/c - Unsigned/signed bytes
+      #   n/v - Unsigned 16-bit words in big/little endian format
+      
+      # You must supply at least one of the dec16, dec8, or dec4 filters and
+      # it must be appropriate for the speed you desire.  The dec4 filter can
+      # be used for all speeds, dec8 only for PSK63, and dec16 only for PSK31.
+
+      def initialize sample_rate, format, frequency, filters
+        @do_dec6 = case sample_rate
+        when 8000 then false
+        when 48000 then true
+        else raise 'only 8000 and 48000 Hz sample rates are supported'
+        end
+        @format = format
+        @sample_size = [0].pack(format).size
+        case ("\x80"*16).unpack(format)[0]
+        when 128
+          @max = 128
+          @offset = -128
+        when -128
+          @max = 128
+          @offset = 0
+        when 32768 + 128
+          @max = 32768
+          @offset = -32768
+        when -32768 + 128
+          @max = 32768
+          @offset = 0
+        else
+          raise 'unable to interpret format'
+        end
+        @frequency = frequency
+        @clock = 8000.0
+        @phase = 0.0
+        @speed = 31
+        @pulse = 0
+        if filters[:dec6]
+          @dec6_coef = NArray.to_na filters[:dec6]*2
+          @dec6_data = NArray.float filters[:dec6].size
+          @dec6_pos = 0
+          @dec6_pulse = 6
+        elsif @do_dec6
+          raise 'no 48000 Hz filter found'
+        end
+        if filters[:dec16]
+          @dec16_coef = NArray.to_na filters[:dec16]*2
+          @dec16_sin = NArray.float filters[:dec16].size
+          @dec16_cos = NArray.float filters[:dec16].size
+          @dec16_pos = 0
+        else
+          @dec16_coef = nil
+        end
+        if filters[:dec8]
+          @dec8_coef = NArray.to_na filters[:dec8]*2
+          @dec8_sin = NArray.float filters[:dec8].size
+          @dec8_cos = NArray.float filters[:dec8].size
+          @dec8_pos = 0
+        else
+          @dec8_coef = nil
+        end
+        if filters[:dec4]
+          @dec4_coef = NArray.to_na filters[:dec4]*2
+          @dec4a_sin = NArray.float filters[:dec4].size
+          @dec4a_cos = NArray.float filters[:dec4].size
+          @dec4a_pos = 0
+          @dec4b_sin = NArray.float filters[:dec4].size
+          @dec4b_cos = NArray.float filters[:dec4].size
+          @dec4b_pos = 0
+        else
+          @dec4_coef = nil
+        end
+        @bit_coef = NArray.to_na filters[:bit]*2
+        @bit_sin = NArray.float filters[:bit].size
+        @bit_cos = NArray.float filters[:bit].size
+        @bit_pos = 0
+        recalc_phase_inc
+      end
+      
+      def reset
+        if @dec6_coef
+          @dec6_data.fill 0.0
+        end
+        if @dec16_coef
+          @dec16_sin_data.fill 0.0 
+          @dec16_cos_data.fill 0.0
+        end
+        if @dec8_coef
+          @dec8_sin_data.fill 0.0
+          @dec8_cos_data.fill 0.0
+        end
+        if @dec4_coef
+          @dec4a_sin.fill 0.0
+          @dec4a_cos.fill 0.0
+          @dec4b_sin.fill 0.0
+          @dec4b_cos.fill 0.0
+        end
+        @bit_sin.fill 0.0
+        @bit_cos.fill 0.0
+      end
+      
+      def recalc_phase_inc
+        @phase_inc = PI2 * @frequency / @clock
+      end
+      
+      def call sample_data
+        raise 'alignment error' unless sample_data.size % @sample_size == 0
+        sample_data.force_encoding('binary') # Ensure slice is fast like Ruby 1.9.3 byteslice
+        mod16_8 = @speed == 63 ? 8 : 16
+        pos = 0
+        while pos < sample_data.size
+          pos += @sample_size
+          sample = sample_data.slice(pos,@sample_size).unpack(@format)[0] || 0
+          sample = (sample + @offset).to_f / @max
+          if @do_dec6
+            @dec6_pos = @dec6_data.size if @dec6_pos == 0
+            @dec6_pos -= 1
+            @dec6_data[@dec6_pos] = sample
+            @dec6_pulse -= 1
+            next unless @dec6_pulse == 0
+            @dec6_pulse = 6
+            sample = @dec6_data.fir(@dec6_coef, @dec6_pos)
+          end
+          @phase += @phase_inc
+          @phase -= PI2 if @phase > PI2
+          if @dec16_coef and @speed == 31
+            @dec16_pos = @dec16_sin.size if @dec16_pos == 0
+            @dec16_pos -= 1
+            @dec16_sin[@dec16_pos] = sample * Math.sin(@phase)
+            @dec16_cos[@dec16_pos] = sample * Math.cos(@phase)
+            next unless ((@pulse = @pulse + 1 & 0xFF) % 16) == 0
+            ival = @dec16_sin.fir(@dec16_coef, @dec16_pos)
+            qval = @dec16_cos.fir(@dec16_coef, @dec16_pos)
+          elsif @dec8_coef and @speed == 63
+            @dec8_pos = @dec8_sin.size if @dec8_pos == 0
+            @dec8_pos -= 1
+            @dec8_sin[@dec8_pos] = sample * Math.sin(@phase)
+            @dec8_cos[@dec8_pos] = sample * Math.cos(@phase)
+            next unless ((@pulse = @pulse + 1 & 0xFF) % 8) == 0
+            ival = @dec8_sin.fir(@dec8_coef, @dec8_pos)
+            qval = @dec8_cos.fir(@dec8_coef, @dec8_pos)
+          elsif @dec4_coef
+            @dec4a_pos = @dec4a_sin.size if @dec4a_pos == 0
+            @dec4a_pos -= 1
+            @dec4a_sin[@dec4a_pos] = sample * Math.sin(@phase)
+            @dec4a_cos[@dec4a_pos] = sample * Math.cos(@phase)
+            next unless ((@pulse = @pulse + 1 & 0xFF) % 4) == 0
+            @dec4b_pos = @dec4b_sin.size if @dec4b_pos == 0
+            @dec4b_pos -= 1
+            ival = @dec4b_sin[@dec4b_pos] = @dec4a_sin.fir(@dec4_coef, @dec4a_pos)
+            qval = @dec4b_cos[@dec4b_pos] = @dec4a_cos.fir(@dec4_coef, @dec4a_pos)
+            next unless @speed == 125 or (@pulse % mod16_8 == 0)
+            unless @speed == 125
+              ival = @dec4b_sin.fir(@dec4_coef, @dec4b_pos)
+              qval = @dec4b_cos.fir(@dec4_coef, @dec4b_pos)
+            end
+          else
+            raise 'no suitable filter found'
+          end
+          @bit_pos = @bit_sin.size if @bit_pos == 0
+          @bit_pos -= 1
+          @bit_sin[@bit_pos] = ival
+          @bit_cos[@bit_pos] = qval
+          yield @bit_sin.fir(@bit_coef, @bit_pos), @bit_cos.fir(@bit_coef, @bit_pos)
+        end
+      end
+      
+      # NArray can easily double filter performance
+      begin
+        require 'narray'
+        class ::NArray
+          def fir filter, pos
+            (self * filter[size-pos..-1-pos]).sum
+          end
+        end
+      rescue LoadError => e
+        # Pure Ruby fake NArray
+        class NArray < Array
+          def self.float arg
+            new arg, 0.0
+          end
+          def self.to_na arg
+            new(arg).freeze
+          end
+          def fir filter, pos
+            acc = 0.0
+            index = size - pos
+            each do |val|
+              acc += val * filter[index]
+              index += 1
+            end
+            acc
+          end
+        end
+      end
+      
+    end
+  end
+end
+

data/lib/radio/psk31/fir_coef.rb

@@ -0,0 +1,287 @@
+class Radio
+  class PSK31
+
+    FIR_DEC6 = [
+       -0.001026086585,
+       -0.002131424398,
+       -0.003473651912,
+       -0.005415991401,
+       -0.007094896149,
+       -0.008685456592,
+       -0.009325598368,
+       -0.009096149588,
+       -0.007453841473,
+       -0.004787743788,
+       -0.001155271744,
+        0.002594565142,
+        0.005930466661,
+        0.007902971637,
+        0.008117998323,
+        0.006209532708,
+        0.002565394219,
+       -0.002211033166,
+       -0.006928668876,
+       -0.010456837484,
+       -0.011638240040,
+       -0.009915015895,
+       -0.005261804435,
+        0.001469200215,
+        0.008865519609,
+        0.015027698257,
+        0.018157576105,
+        0.016869453853,
+        0.010745489111,
+        0.000474997075,
+       -0.012015135186,
+       -0.023901641616,
+       -0.031878023716,
+       -0.032876349281,
+       -0.024689692215,
+       -0.006586586958,
+        0.020417464549,
+        0.053543134171,
+        0.088652482149,
+        0.120891532319,
+        0.145553464727,
+        0.158910442030,
+        0.158910442030,
+        0.145553464727,
+        0.120891532319,
+        0.088652482149,
+        0.053543134171,
+        0.020417464549,
+       -0.006586586958,
+       -0.024689692215,
+       -0.032876349281,
+       -0.031878023716,
+       -0.023901641616,
+       -0.012015135186,
+        0.000474997075,
+        0.010745489111,
+        0.016869453853,
+        0.018157576105,
+        0.015027698257,
+        0.008865519609,
+        0.001469200215,
+       -0.005261804435,
+       -0.009915015895,
+       -0.011638240040,
+       -0.010456837484,
+       -0.006928668876,
+       -0.002211033166,
+        0.002565394219,
+        0.006209532708,
+        0.008117998323,
+        0.007902971637,
+        0.005930466661,
+        0.002594565142,
+       -0.001155271744,
+       -0.004787743788,
+       -0.007453841473,
+       -0.009096149588,
+       -0.009325598368,
+       -0.008685456592,
+       -0.007094896149,
+       -0.005415991401,
+       -0.003473651912,
+       -0.002131424398,
+       -0.001026086585
+    ].freeze
+    
+    
+    # 64-tap raised-cosine
+    FIR_DEC16 = [
+      0.000000,
+      0.000038,
+      0.000150,
+      0.000336,
+      0.000595,
+      0.000922,
+      0.001317,
+      0.001773,
+      0.002288,
+      0.002856,
+      0.003472,
+      0.004130,
+      0.004823,
+      0.005545,
+      0.006288,
+      0.007047,
+      0.007812,
+      0.008578,
+      0.009337,
+      0.010080,
+      0.010802,
+      0.011495,
+      0.012153,
+      0.012769,
+      0.013337,
+      0.013852,
+      0.014308,
+      0.014703,
+      0.015030,
+      0.015289,
+      0.015475,
+      0.015587,
+      0.015625,
+      0.015587,
+      0.015475,
+      0.015289,
+      0.015030,
+      0.014703,
+      0.014308,
+      0.013852,
+      0.013337,
+      0.012769,
+      0.012153,
+      0.011495,
+      0.010802,
+      0.010080,
+      0.009337,
+      0.008578,
+      0.007813,
+      0.007047,
+      0.006288,
+      0.005545,
+      0.004823,
+      0.004130,
+      0.003472,
+      0.002856,
+      0.002288,
+      0.001773,
+      0.001317,
+      0.000922,
+      0.000595,
+      0.000336,
+      0.000150,
+      0.000038
+    ].freeze
+    
+    
+    #  Design method: Parks-McClellan method
+    #  Number of taps  =  35
+    #  Number of bands =  2 
+    #  Band     Lower       Upper       Value      Weight
+    #            edge        edge
+    #   1         0.0       .0125         1.0         1
+    #   2         .125      .5          .000001       10  
+    FIR_DEC4 = [
+     -0.00021203644,
+     -0.00070252426,
+      -0.0016680526,
+      -0.0031934799,
+      -0.0051899752,
+      -0.0072862086,
+      -0.0087714235,
+      -0.0086272102,
+      -0.0056735648,
+       0.0011784719,
+         0.01261353,
+        0.028615709,
+        0.048280707,
+        0.069812051,
+        0.090735013,
+         0.10830381,
+         0.12001897,
+         0.12413413,
+         0.12001897,
+         0.10830381,
+        0.090735013,
+        0.069812051,
+        0.048280707,
+        0.028615709,
+         0.01261353,
+       0.0011784719,
+      -0.0056735648,
+      -0.0086272102,
+      -0.0087714235,
+      -0.0072862086,
+      -0.0051899752,
+      -0.0031934799,
+      -0.0016680526,
+     -0.00070252426,
+     -0.00021203644
+    ].freeze
+
+    
+    # 16 Hz bw LP filter for data recovery
+    #  Filter type:   Multiband filter
+    #  Design method: Parks-McClellan method
+    #  Number of taps  =  65
+    #  Number of bands =  2 
+    #  Band   Lower       Upper       Value      Weight
+    #         edge        edge
+    #
+    #	1       0.0        .03125         1.         1.
+    #	2        .0625     .5          .0000       286.
+    FIR_BIT = [
+     4.3453566e-005,
+     -0.00049122414,
+     -0.00078771292,
+      -0.0013507826,
+      -0.0021287814,
+       -0.003133466,
+       -0.004366817,
+      -0.0058112187,
+      -0.0074249976,
+      -0.0091398882,
+       -0.010860157,
+       -0.012464086,
+       -0.013807772,
+       -0.014731191,
+       -0.015067057,
+       -0.014650894,
+       -0.013333425,
+        -0.01099166,
+      -0.0075431246,
+      -0.0029527849,
+       0.0027546292,
+       0.0094932775,
+        0.017113308,
+        0.025403511,
+        0.034099681,
+        0.042895839,
+        0.051458575,
+        0.059444853,
+        0.066521003,
+        0.072381617,
+        0.076767694,
+        0.079481619,
+        0.080420311,
+        0.079481619,
+        0.076767694,
+        0.072381617,
+        0.066521003,
+        0.059444853,
+        0.051458575,
+        0.042895839,
+        0.034099681,
+        0.025403511,
+        0.017113308,
+       0.0094932775,
+       0.0027546292,
+      -0.0029527849,
+      -0.0075431246,
+        -0.01099166,
+       -0.013333425,
+       -0.014650894,
+       -0.015067057,
+       -0.014731191,
+       -0.013807772,
+       -0.012464086,
+       -0.010860157,
+      -0.0091398882,
+      -0.0074249976,
+      -0.0058112187,
+       -0.004366817,
+       -0.003133466,
+      -0.0021287814,
+      -0.0013507826,
+     -0.00078771292,
+     -0.00049122414,
+     4.3453566e-005
+    ].freeze
+
+
+  end
+end

data/lib/radio/psk31/rx.rb

@@ -0,0 +1,52 @@
+class Radio
+  
+  class PSK31
+    
+    class Rx
+
+      def initialize sample_rate, format, frequency
+        @filter = Filters.new sample_rate, format, frequency, 
+          :dec6 => FIR_DEC6,
+          :dec16 => FIR_DEC16,
+          :dec4 => FIR_DEC4,
+          :bit => FIR_BIT
+        @bit_detect = BitDetect.new
+        @decoder = Decoder.new
+      end
+      
+      def call sample_data
+        decoded = ''
+        @filter.call sample_data do |i, q|
+          @bit_detect.call i, q do
+            @decoder.call i, q do |symbol|
+              decoded += symbol
+            end
+          end
+        end
+        decoded
+      end
+      
+      def frequency= frequency
+        if frequency != @filter.frequency
+          @filter.frequency = frequency
+          @filter.recalc_phase_inc
+          reset
+        end
+      end
+      
+      # To compensate for bad clock in A-D conversion
+      def adjust_sample_clock ppm
+        @filter.clock = 8000.0 * ppm / 1000000 + 8000
+        @filter.recalc_phase_inc
+      end
+      
+      def reset
+        @filter.reset
+        @bit_detect.reset
+        @decoder.reset
+      end
+      
+    end
+    
+  end
+end

data/lib/radio/psk31/varicode.rb

@@ -0,0 +1,275 @@
+class Radio
+  
+  class PSK31
+
+    VARICODE_TABLE = [ 
+      0xAAC0, #   0 1010101011
+      0xB6C0, #   1 1011011011
+      0xBB40, #   2 1011101101
+      0xDDC0, #   3 1101110111
+      0xBAC0, #   4 1011101011
+      0xD7C0, #   5 1101011111
+      0xBBC0, #   6 1011101111
+      0xBF40, #   7 1011111101
+      0xBFC0, #   8 1011111111
+      0xEF00, #   9 11101111
+      0xE800, #  10 11101
+      0xDBC0, #  11 1101101111
+      0xB740, #  12 1011011101
+      0xF800, #  13 11111
+      0xDD40, #  14 1101110101
+      0xEAC0, #  15 1110101011
+      0xBDC0, #  16 1011110111
+      0xBD40, #  17 1011110101
+      0xEB40, #  18 1110101101
+      0xEBC0, #  19 1110101111
+      0xD6C0, #  20 1101011011
+      0xDAC0, #  21 1101101011
+      0xDB40, #  22 1101101101
+      0xD5C0, #  23 1101010111
+      0xDEC0, #  24 1101111011
+      0xDF40, #  25 1101111101
+      0xEDC0, #  26 1110110111
+      0xD540, #  27 1101010101
+      0xD740, #  28 1101011101
+      0xEEC0, #  29 1110111011
+      0xBEC0, #  30 1011111011
+      0xDFC0, #  31 1101111111
+      0x8000, # ' ' 1
+      0xFF80, # '!' 111111111
+      0xAF80, # '"' 101011111
+      0xFA80, # '#' 111110101
+      0xED80, # '$' 111011011
+      0xB540, # '%' 1011010101
+      0xAEC0, # '&' 1010111011
+      0xBF80, # ''' 101111111
+      0xFB00, # '(' 11111011
+      0xF700, # ')' 11110111
+      0xB780, # '*' 101101111
+      0xEF80, # '+' 111011111
+      0xEA00, # ',' 1110101
+      0xD400, # '-' 110101
+      0xAE00, # '.' 1010111
+      0xD780, # '/' 110101111
+      0xB700, # '0' 10110111
+      0xBD00, # '1' 10111101
+      0xED00, # '2' 11101101
+      0xFF00, # '3' 11111111
+      0xBB80, # '4' 101110111
+      0xAD80, # '5' 101011011
+      0xB580, # '6' 101101011
+      0xD680, # '7' 110101101
+      0xD580, # '8' 110101011
+      0xDB80, # '9' 110110111
+      0xF500, # ':' 11110101
+      0xDE80, # ';' 110111101
+      0xF680, # '<' 111101101
+      0xAA00, # '=' 1010101
+      0xEB80, # '>' 111010111
+      0xABC0, # '?' 1010101111
+      0xAF40, # '@' 1010111101
+      0xFA00, # 'A' 1111101
+      0xEB00, # 'B' 11101011
+      0xAD00, # 'C' 10101101
+      0xB500, # 'D' 10110101
+      0xEE00, # 'E' 1110111
+      0xDB00, # 'F' 11011011
+      0xFD00, # 'G' 11111101
+      0xAA80, # 'H' 101010101
+      0xFE00, # 'I' 1111111
+      0xFE80, # 'J' 111111101
+      0xBE80, # 'K' 101111101
+      0xD700, # 'L' 11010111
+      0xBB00, # 'M' 10111011
+      0xDD00, # 'N' 11011101
+      0xAB00, # 'O' 10101011
+      0xD500, # 'P' 11010101
+      0xEE80, # 'Q' 111011101
+      0xAF00, # 'R' 10101111
+      0xDE00, # 'S' 1101111
+      0xDA00, # 'T' 1101101
+      0xAB80, # 'U' 101010111
+      0xDA80, # 'V' 110110101
+      0xAE80, # 'W' 101011101
+      0xBA80, # 'X' 101110101
+      0xBD80, # 'Y' 101111011
+      0xAB40, # 'Z' 1010101101
+      0xFB80, # '[' 111110111
+      0xF780, # '\' 111101111
+      0xFD80, # ']' 111111011
+      0xAFC0, # '^' 1010111111
+      0xB680, # '_' 101101101
+      0xB7C0, # '`' 1011011111
+      0xB000, # 'a' 1011
+      0xBE00, # 'b' 1011111
+      0xBC00, # 'c' 101111
+      0xB400, # 'd' 101101
+      0xC000, # 'e' 11
+      0xF400, # 'f' 111101
+      0xB600, # 'g' 1011011
+      0xAC00, # 'h' 101011
+      0xD000, # 'i' 1101
+      0xF580, # 'j' 111101011
+      0xBF00, # 'k' 10111111
+      0xD800, # 'l' 11011
+      0xEC00, # 'm' 111011
+      0xF000, # 'n' 1111
+      0xE000, # 'o' 111
+      0xFC00, # 'p' 111111
+      0xDF80, # 'q' 110111111
+      0xA800, # 'r' 10101
+      0xB800, # 's' 10111
+      0xA000, # 't' 101
+      0xDC00, # 'u' 110111
+      0xF600, # 'v' 1111011
+      0xD600, # 'w' 1101011
+      0xDF00, # 'x' 11011111
+      0xBA00, # 'y' 1011101
+      0xEA80, # 'z' 111010101
+      0xADC0, # '{' 1010110111
+      0xDD80, # '|' 110111011
+      0xAD40, # '}' 1010110101
+      0xB5C0, # '~' 1011010111
+      0xED40, # 127 1110110101
+      0xEF40, # 128 1110111101
+      0xEFC0, # 129 1110111111
+      0xF540, # 130 1111010101
+      0xF5C0, # 131 1111010111
+      0xF6C0, # 132 1111011011
+      0xF740, # 133 1111011101
+      0xF7C0, # 134 1111011111
+      0xFAC0, # 135 1111101011
+      0xFB40, # 136 1111101101
+      0xFBC0, # 137 1111101111
+      0xFD40, # 138 1111110101
+      0xFDC0, # 139 1111110111
+      0xFEC0, # 140 1111111011
+      0xFF40, # 141 1111111101
+      0xFFC0, # 142 1111111111
+      0xAAA0, # 143 10101010101
+      0xAAE0, # 144 10101010111
+      0xAB60, # 145 10101011011
+      0xABA0, # 146 10101011101
+      0xABE0, # 147 10101011111
+      0xAD60, # 148 10101101011
+      0xADA0, # 149 10101101101
+      0xADE0, # 150 10101101111
+      0xAEA0, # 151 10101110101
+      0xAEE0, # 152 10101110111
+      0xAF60, # 153 10101111011
+      0xAFA0, # 154 10101111101
+      0xAFE0, # 155 10101111111
+      0xB560, # 156 10110101011
+      0xB5A0, # 157 10110101101
+      0xB5E0, # 158 10110101111
+      0xB6A0, # 159 10110110101
+      0xB6E0, # 160 10110110111
+      0xB760, # 161 10110111011
+      0xB7A0, # 162 10110111101
+      0xB7E0, # 163 10110111111
+      0xBAA0, # 164 10111010101
+      0xBAE0, # 165 10111010111
+      0xBB60, # 166 10111011011
+      0xBBA0, # 167 10111011101
+      0xBBE0, # 168 10111011111
+      0xBD60, # 169 10111101011
+      0xBDA0, # 170 10111101101
+      0xBDE0, # 171 10111101111
+      0xBEA0, # 172 10111110101
+      0xBEE0, # 173 10111110111
+      0xBF60, # 174 10111111011
+      0xBFA0, # 175 10111111101
+      0xBFE0, # 176 10111111111
+      0xD560, # 177 11010101011
+      0xD5A0, # 178 11010101101
+      0xD5E0, # 179 11010101111
+      0xD6A0, # 180 11010110101
+      0xD6E0, # 181 11010110111
+      0xD760, # 182 11010111011
+      0xD7A0, # 183 11010111101
+      0xD7E0, # 184 11010111111
+      0xDAA0, # 185 11011010101
+      0xDAE0, # 186 11011010111
+      0xDB60, # 187 11011011011
+      0xDBA0, # 188 11011011101
+      0xDBE0, # 189 11011011111
+      0xDD60, # 190 11011101011
+      0xDDA0, # 191 11011101101
+      0xDDE0, # 192 11011101111
+      0xDEA0, # 193 11011110101
+      0xDEE0, # 194 11011110111
+      0xDF60, # 195 11011111011
+      0xDFA0, # 196 11011111101
+      0xDFE0, # 197 11011111111
+      0xEAA0, # 198 11101010101
+      0xEAE0, # 199 11101010111
+      0xEB60, # 200 11101011011
+      0xEBA0, # 201 11101011101
+      0xEBE0, # 202 11101011111
+      0xED60, # 203 11101101011
+      0xEDA0, # 204 11101101101
+      0xEDE0, # 205 11101101111
+      0xEEA0, # 206 11101110101
+      0xEEE0, # 207 11101110111
+      0xEF60, # 208 11101111011
+      0xEFA0, # 209 11101111101
+      0xEFE0, # 210 11101111111
+      0xF560, # 211 11110101011
+      0xF5A0, # 212 11110101101
+      0xF5E0, # 213 11110101111
+      0xF6A0, # 214 11110110101
+      0xF6E0, # 215 11110110111
+      0xF760, # 216 11110111011
+      0xF7A0, # 217 11110111101
+      0xF7E0, # 218 11110111111
+      0xFAA0, # 219 11111010101
+      0xFAE0, # 220 11111010111
+      0xFB60, # 221 11111011011
+      0xFBA0, # 222 11111011101
+      0xFBE0, # 223 11111011111
+      0xFD60, # 224 11111101011
+      0xFDA0, # 225 11111101101
+      0xFDE0, # 226 11111101111
+      0xFEA0, # 227 11111110101
+      0xFEE0, # 228 11111110111
+      0xFF60, # 229 11111111011
+      0xFFA0, # 230 11111111101
+      0xFFE0, # 231 11111111111
+      0xAAB0, # 232 101010101011
+      0xAAD0, # 233 101010101101
+      0xAAF0, # 234 101010101111
+      0xAB50, # 235 101010110101
+      0xAB70, # 236 101010110111
+      0xABB0, # 237 101010111011
+      0xABD0, # 238 101010111101
+      0xABF0, # 239 101010111111
+      0xAD50, # 240 101011010101
+      0xAD70, # 241 101011010111
+      0xADB0, # 242 101011011011
+      0xADD0, # 243 101011011101
+      0xADF0, # 244 101011011111
+      0xAEB0, # 245 101011101011
+      0xAED0, # 246 101011101101
+      0xAEF0, # 247 101011101111
+      0xAF50, # 248 101011110101
+      0xAF70, # 249 101011110111
+      0xAFB0, # 250 101011111011
+      0xAFD0, # 251 101011111101
+      0xAFF0, # 252 101011111111
+      0xB550, # 253 101101010101
+      0xB570, # 254 101101010111
+      0xB5B0  # 255 101101011011
+    ].freeze
+
+    VARICODE_DECODE_TABLE = []
+    for i in 0...256
+      wTemp = VARICODE_TABLE[i]
+      wTemp >>= 4
+      wTemp >>= 1 until (wTemp&1 == 1)
+      wTemp >>= 1
+      VARICODE_DECODE_TABLE[wTemp] = i.chr
+    end
+    VARICODE_DECODE_TABLE.freeze
+    
+  end
+end

metadata

@@ -0,0 +1,54 @@
+--- !ruby/object:Gem::Specification
+name: radio
+version: !ruby/object:Gem::Version
+  version: 0.0.1
+  prerelease: 
+platform: ruby
+authors:
+- David Turnbull
+autorequire: 
+bindir: bin
+cert_chain: []
+date: 2011-11-11 00:00:00.000000000Z
+dependencies: []
+description: 
+email:
+- dturnbull@gmail.com
+executables: []
+extensions: []
+extra_rdoc_files: []
+files:
+- lib/radio/psk31/bit_detect.rb
+- lib/radio/psk31/decoder.rb
+- lib/radio/psk31/filters.rb
+- lib/radio/psk31/fir_coef.rb
+- lib/radio/psk31/rx.rb
+- lib/radio/psk31/varicode.rb
+- lib/radio.rb
+- README.md
+homepage: https://github.com/dturnbull/radio
+licenses: []
+post_install_message: 
+rdoc_options: []
+require_paths:
+- lib
+required_ruby_version: !ruby/object:Gem::Requirement
+  none: false
+  requirements:
+  - - ! '>='
+    - !ruby/object:Gem::Version
+      version: '0'
+required_rubygems_version: !ruby/object:Gem::Requirement
+  none: false
+  requirements:
+  - - ! '>='
+    - !ruby/object:Gem::Version
+      version: '0'
+requirements: []
+rubyforge_project: 
+rubygems_version: 1.8.6
+signing_key: 
+specification_version: 3
+summary: Amateur radio software
+test_files: []
+has_rdoc: