t_lib 0.1.0 → 0.1.0.1

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA1:
-  metadata.gz: 71a976dc879dcb427fc7547d10efa0b3b9894924
-  data.tar.gz: 13889f6b40136de79fce29832389def051ba3772
+  metadata.gz: 9821928a9b695eb9a02d9508c512dc82baa21119
+  data.tar.gz: 815e887eeba7749e2a80aafb3944bb20898e8d4b
 SHA512:
-  metadata.gz: a6c04a9b669e1bd9c58d058bf82269822355626016ff91201a7224ce58562a6294000d8f0404b0ffb68bb9ced23933726a89c3be08049c5327e89b380eea5954
-  data.tar.gz: cae4d1a937666d5a5f6eb6bfeb7118b94171f2e343d4d796699aa05a011800c88b5ca71566c8a0e47bf7356dbad1afbf949edb2d8dedb190a79019bbfcb51047
+  metadata.gz: 7a26f0ffcc3d7da5b84934772261a1919332dd6cbd56673b1693dfdb24317f4a015b6246eb96dfd0406c3350905d57ac9daf9936d2da71de6c289f0a876fa437
+  data.tar.gz: e28be8515a0263556be90f0c719e90d07c9d928619f60c2bcad2e7022f56b086ea73360500c8d3c183d7a81373f8afd21486eb2c7ba924b1d00eb5da98c82e9d

data/lib/t_lib/math.rb

@@ -0,0 +1,119 @@
+#!/usr/bin/ruby
+# -*- encoding: utf-8 -*-
+
+module TLib
+  module Math
+
+    #
+    # === gauusian distribution 
+    # @todo
+    #
+    def gauusian(x, mu, sigma)
+      f1 = 1.0/(Math.sqrt(2.0*Math::PI)*Math.sqrt(sigma))
+      f2 = Math.exp(-(((x-mu)**2)/((2.0*sigma))))
+      return f1 * f2
+    end
+
+    #
+    # === gauusian distribution 
+    # @todo
+    #
+    def gauusian_array(x, mu, sigma)
+      if x[0].size <= 1
+        x = x[0]
+        mu = mu[0]
+        sigma = sigma[0][0]
+        return gauusian(x, mu, sigma)
+      else
+        return gauusian_over_2dim(x, mu, sigma)
+      end
+    end
+
+
+    #
+    # === gauusian distribution over 2 dim version
+    #
+    def gauusian_over_2dim(x, mu, conv)
+      x = Matrix[x]
+      mu = Matrix[mu]
+      conv = Matrix[*conv]
+      f1 = 1.0/(((2.0 * Math::PI)**(@dim/2.0)) * ( conv.det**(0.5) ))
+      f2 = Math.exp((-1.0/2.0)*((x-mu) * conv.inverse * (x-mu).transpose)[0, 0])
+
+      return (f1 * f2)
+    end
+
+    #
+    # ===ボックス―ミューラー法をよる正規分布乱数発生
+    # @param mu flout 平均
+    # @param sigma flout 標準偏差
+    # @return ボックスミューラー法に従う正規分布に従う乱数を生成
+    #
+    def normal_rand(mu = 0,sigma = 1.0)
+      a, b = rand(), rand() ;
+      return (Math.sqrt(-2*Math.log(rand()))*Math.sin(2*Math::PI*rand()) * sigma) + mu
+    end
+
+    #
+    # === 指数分布に従った乱数を返す
+    #
+    def exp(ave=0.5)
+      x = rand() ;
+      y = rand() ;
+      while ((1.0/ave)*Math.exp(-(x/ave)) < y )
+        x = rand() ;
+        y = rand() ;
+      end
+      return x ;
+    end
+
+
+
+    #
+    # ===ポアソン分布に従う乱数を発生する
+    #
+    def poisson_rand(mu=0.0)
+      lambda = Math.exp(-mu)
+      k = 0
+      p = 1.0
+      while p >= lambda
+        p *= rand()
+        k += 1
+      end
+      return k - 1
+    end
+    def scale(x)
+      sum_each_vec = []
+      ave_list = []
+      std_list = []
+      x.each{|vec| 
+        vec.each_with_index{|data, i|
+          sum_each_vec[i] = (sum_each_vec[i] == nil) ? data : sum_each_vec[i]+data
+        }
+      }
+      x[0].size.times{|i|
+        ave_list.push(sum_each_vec[i]/x.size)
+      }
+
+      sum_each_vec = []
+      x.each{|vec| 
+        vec.each_with_index{|data, i|
+          sum_each_vec[i] = (sum_each_vec[i] == nil) ? (ave_list[i]-data)**2 : (sum_each_vec[i]+(ave_list[i]-data)**2)
+        }
+      }
+      x[0].size.times{|i|
+        std_list.push(Math.sqrt(sum_each_vec[i]/x.size))
+      }
+
+      scaled_x = []
+      x.each_with_index{|vec, i| 
+        scaled_x[i] ||= []
+        vec.each_with_index{|data, j|
+          scaled_x[i][j] ||= (data-ave_list[j])/std_list[j]
+        }
+      }
+      return scaled_x
+    end
+
+  end
+end

data/lib/t_lib/version.rb

@@ -1,3 +1,3 @@
 module TLib
-  VERSION = "0.1.0"
+  VERSION = "0.1.0.1"
 end

data/lib/t_lib.rb

@@ -2,9 +2,9 @@ require "t_lib/version"
 
 module TLib
   # Your code goes here...
-  self.write_json(file_name, data) {
+  def self.write_json(file_name, data) 
     File.open(file_name) do |file|
       JSON.dump(data, file)
     end
-  }
+  end
 end

metadata

@@ -1,14 +1,14 @@
 --- !ruby/object:Gem::Specification
 name: t_lib
 version: !ruby/object:Gem::Version
-  version: 0.1.0
+  version: 0.1.0.1
 platform: ruby
 authors:
 - tcom
 autorequire: 
 bindir: exe
 cert_chain: []
-date: 2017-06-24 00:00:00.000000000 Z
+date: 2017-06-27 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: bundler
@@ -70,6 +70,7 @@ files:
 - bin/console
 - bin/setup
 - lib/t_lib.rb
+- lib/t_lib/math.rb
 - lib/t_lib/version.rb
 - t_lib.gemspec
 homepage: https://github.com/Tcom242242/t_lib