mathpack 0.3.0 → 0.3.1

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA1:
-  metadata.gz: e0f0d47957524843b49cef7c52ac64230e6f7d68
-  data.tar.gz: cb662a79d6977fb3f2d891f2b2eb2881b2c00cf9
+  metadata.gz: 990d138b20708752c754e294a22627435a3a5ce9
+  data.tar.gz: dbb63f6b4cea510030e05c819cd2bac9a5c7376a
 SHA512:
-  metadata.gz: 1d134a1000e9bf0068a36923484f373cf546492e8a6795539ee8616f06d09ec6e31886e096929f6c463a23a8d48dbcbf973b035309d3ccac92dceec6973f88e0
-  data.tar.gz: 0a8fc14ec5126df7b9afe59ca2141b52539cea8a710875c6ca6bb10e5aae7fc4238cca77a0f6b359be6bc6d0b215f85d220439f0f5b99475592616d8ed264e30
+  metadata.gz: 36318534efa67765b74962c5d6b102b7e0e4c4b2e2bbda35a795c443be5d6e7b09f6ec9be97f014ae7e8db88a6831ef710afb42714164ead92322b35f3c9433a
+  data.tar.gz: 38db3fd63908159911ef2ba99d8b6a0b899ae26010482104db5aa384f42f61f1176dcec6c59672fe664fb342f4f6817a5540a3d3289950cfaf4d382212194338

data/.gitignore

@@ -21,3 +21,4 @@ tmp
 *.a
 mkmf.log
 .rubocop.yml
+*.idea

data/README.md

@@ -19,7 +19,7 @@ Or install it yourself as:
     $ gem install mathpack
 ## Information
 Gem `mathpack` allows to count statistical functions through `Statistics` class, solve nonlinear equations through `Equation` module
-`solve` method, solve systems of linear equations through `SLE` module `solve` method.
+`solve` method, solve systems of linear equations through `SLE` module `solve` method, approximate functions by polynom through `approximate_by_polynom` method of `Approximation` module.
 ## Statistics
 `Statistics` class have following methods
 - **number** - returns a number of elements in series
@@ -35,6 +35,7 @@ Gem `mathpack` allows to count statistical functions through `Statistics` class,
 - **empirical_pdf** - returns *empirical probability density function* value in some point
 - **print_empirical_cdf_to_csv** - allows to print empirical_cdf line chart values to `.csv` file with name *filename*
 - **print_empirical_cdf_to_csv** - allows to print empirical_cdf line chart values to `.csv` file with name *filename*
+- **trend** - returns trend polynom
 
 ### Usage
 ```ruby
@@ -52,13 +53,14 @@ stat.empirical_cdf(5.5) #=> 0.75
 stat.empirical_pdf(3) #=> 0.07639393483317147
 stat.print_empirical_cdf_to_csv('cdf.csv') #=> nil
 stat.print_empirical_pdf_to_csv('pdf.csv') #=> nil
+stat.trend(polynom_power: 1) #=> 1.7999999999999996*x - 0.9999999999999987
 ```
 
-##Equation
-`Equation` module has only one method
+## Equation
+`Equation` module has only one method:
 - **solve** - method, which allows to solve *nonlinear equations*. Neccessary params are `eps` representing calculations accuraccy and `start` representing point to start root search
 
-###Usage
+### Usage
 Now you have no problems solving **nonlinear equations**. If you want, for example, to solve equation ![equation](http://latex.codecogs.com/gif.latex?x%5E%7B2%7D%20%3D%20%5Csin%28%7Bx+1%7D%29)
 
 You need to complete the following steps:
@@ -80,10 +82,10 @@ Mathpack::Equation.solve(start: 0.01, eps: 0.00001){|x| Math.exp(x-2) - Math.sin
 ```
 
 ## SLE
-`SLE` module has only one method
+`SLE` module has only one method:
 - **solve** - method, which allows to solve *system of linear equations*. Neccessary params are `matrix` representing system matrix and `f`. Params can be **Array** or **Matrix** class. If the system can't be solved, method raise exception
 
-## Usage
+### Usage
 Let's solve some system of linear equations. It can be written as
 
 ![equation](http://latex.codecogs.com/gif.latex?%5Csmall%20AX%20%3D%20B)
@@ -104,6 +106,36 @@ b = Matrix.row_vector [15, 30, 15]
 Mathpack::SLE.solve(matrix: a, f: b) #=> Matrix[[-1.0, 2.0, 4.0]]
 ```
 
+## Approximation
+May be sometimes you need to approximate some function, which is analytic or represented by table of values, by polynom. You can
+manage with this problem using `Approximation` module. This module uses least squares approximation method.
+`Approximation` module has the following methods:
+- **approximate_by_polynom** - main method of the module, that allows to approximate function. It returns array of coefficients of polynom.
+- **generate_nodes** - method, that allows to generate nodes for approximation with some step.
+- **print_polynom** - method, that returns a string representing polynom with given coefficients.
+
+### Usage
+```ruby
+# Function to print polynom having it's coefficients
+Mathpack::Approximation::print_polynom([1, -2, 1]) #=> x^2 - 2*x + 1
+
+# Function to generate nodes for approximation. Only choose start value, end value and step
+x = Mathpack::Approximation::generate_nodes(start: 0, end: 10, step: 0.25)
+#=> [0, 0.25, ..., 10]
+
+# Function that approximate given function by polynom with power polynom_power in approximation nodes x and returns coefficients of approximating polynom
+result = Mathpack::Approximation::approximate_by_polynom(x: x, polynom_power: 5){ |x| Math.sin(x) }
+#=> [0.0008009744982571882, -0.030619986086758588, 0.3805927651948083, -1.8481035413353888, 2.985465287759307, -0.3873066069630569]
+
+# May be you want to print this polynom, let's call print_polynom function
+Mathpack::Approximation::print_polynom(result)
+#=> 0.0008009744982571882*x^5 - 0.030619986086758588*x^4 + 0.3805927651948083*x^3 - 1.8481035413353888*x^2 + 2.985465287759307*x - 0.3873066069630569
+
+# If you have a table of values, where x - array of arguement values and f - array of function values, call approximate by polynom function with parameter f instead of block
+result = Mathpack::Approximation::approximate_by_polynom(x: [1, 2, 3], f: [1, 4, 9], polynom_power: 2) #=> [1, 0, 0]
+Mathpack::Approximation::print_polynom(result) #=> x^2
+```
+
 ## Contributing
 
 1. Fork it ( https://github.com/[my-github-username]/mathpack/fork )

data/lib/mathpack/statistics.rb

@@ -82,6 +82,12 @@ module Mathpack
       end
     end
 
+    def trend(params = {})
+      numbers = Array.new(number){ |i| i + 1 }
+      polynom = Mathpack::Approximation::approximate_by_polynom(x: numbers, f: @series, polynom_power: params[:polynom_power])
+      Mathpack::Approximation.print_polynom(polynom)
+    end
+
     private
 
     def heaviside(x)

data/lib/mathpack/version.rb

@@ -1,3 +1,3 @@
 module Mathpack
-  VERSION = '0.3.0'
+  VERSION = '0.3.1'
 end

data/mathpack.gemspec

@@ -19,5 +19,5 @@ Gem::Specification.new do |spec|
   spec.require_paths = ["lib"]
 
   spec.add_development_dependency 'bundler', '~> 1.6'
-  spec.add_development_dependency 'rake'
+  spec.add_development_dependency 'rake', '~>0'
 end

data/spec/approximation_spec.rb

@@ -9,30 +9,30 @@ describe 'Approximation' do
     let(:values) { Mathpack::Approximation.fill_f(nodes){ |val| val**2 } }
 
     it 'should generate nodes' do
-      expect(nodes.length).to be_eql(41)
-      expect(nodes.first).to be_eql(x_start)
-      expect(nodes.last).to be_eql(x_end)
+      expect(nodes.length).to eql(41)
+      expect(nodes.first).to eql(x_start)
+      expect(nodes.last).to eql(x_end)
     end
 
     it 'should fill function table' do
-      expect(values.length).to be_eql(nodes.length)
-      expect(values.last).to be_eql(nodes.last**2)
+      expect(values.length).to eql(nodes.length)
+      expect(values.last).to eql(nodes.last**2)
     end
 
     it 'should print polynoms' do
-      expect(Mathpack::Approximation.print_polynom([5])).to be_eql('5')
-      expect(Mathpack::Approximation.print_polynom([2, 3])).to be_eql('2*x + 3')
-      expect(Mathpack::Approximation.print_polynom([1, -1, 1])).to be_eql('x^2 - x + 1')
-      expect(Mathpack::Approximation.print_polynom([1, 3, 5])).to be_eql('x^2 + 3*x + 5')
-      expect(Mathpack::Approximation.print_polynom([1, -2, 3, -4])).to be_eql('x^3 - 2*x^2 + 3*x - 4')
-      expect(Mathpack::Approximation.print_polynom([1, 0, 1])).to be_eql('x^2 + 1')
+      expect(Mathpack::Approximation.print_polynom([5])).to eql('5')
+      expect(Mathpack::Approximation.print_polynom([2, 3])).to eql('2*x + 3')
+      expect(Mathpack::Approximation.print_polynom([1, -1, 1])).to eql('x^2 - x + 1')
+      expect(Mathpack::Approximation.print_polynom([1, 3, 5])).to eql('x^2 + 3*x + 5')
+      expect(Mathpack::Approximation.print_polynom([1, -2, 3, -4])).to eql('x^3 - 2*x^2 + 3*x - 4')
+      expect(Mathpack::Approximation.print_polynom([1, 0, 1])).to eql('x^2 + 1')
     end
 
     it 'should approximate by polynom' do
       coefficients = Mathpack::Approximation.approximate_by_polynom(x: nodes, f: values, polynom_power: 2)
-      expect(coefficients).to be_eql([1.0, 0.0, 0.0])
+      expect(coefficients).to eql([1.0, 0.0, 0.0])
       coefficients = Mathpack::Approximation.approximate_by_polynom(x: nodes, f: values, polynom_power: 0)
-      expect(coefficients).to be_eql([33.75])
+      expect(coefficients).to eql([33.75])
     end
   end
 end

data/spec/statistics_spec.rb

@@ -59,5 +59,9 @@ describe 'Statistics' do
     it 'should calculate empirical pdf values in points' do
       expect(stat.empirical_pdf(stat.mean)).to eql(0.1882412842233359)
     end
+
+    it 'should find trend' do
+      expect(stat.trend(polynom_power: 2)).to eql('0.0075757575757576055*x^2 + 0.2681818181818179*x + 2.233333333333334')
+    end
   end
 end

metadata

@@ -1,14 +1,14 @@
 --- !ruby/object:Gem::Specification
 name: mathpack
 version: !ruby/object:Gem::Version
-  version: 0.3.0
+  version: 0.3.1
 platform: ruby
 authors:
 - maxmilan
 autorequire: 
 bindir: bin
 cert_chain: []
-date: 2015-01-25 00:00:00.000000000 Z
+date: 2015-02-04 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: bundler
@@ -28,14 +28,14 @@ dependencies:
   name: rake
   requirement: !ruby/object:Gem::Requirement
     requirements:
-    - - ">="
+    - - "~>"
       - !ruby/object:Gem::Version
         version: '0'
   type: :development
   prerelease: false
   version_requirements: !ruby/object:Gem::Requirement
     requirements:
-    - - ">="
+    - - "~>"
       - !ruby/object:Gem::Version
         version: '0'
 description: Includes collection of mathematical methods