XIRR 实现中的 bigdecimal/newton 错误

Question

我想使用 Ruby 中的 bigdecimal/newton 模块来实现 XIRR。我编写了一个脚本来尝试一下 这个例子。

当我运行代码（Mac OS X 10.6 上的 Ruby 1.9.2）时，出现以下错误：

/opt/local/lib/ruby1.9/1.9.1/bigdecimal/ludcmp.rb:84:in `div': wrong number of arguments(2 for 1) (ArgumentError)
   from /opt/local/lib/ruby1.9/1.9.1/bigdecimal/ludcmp.rb:84:in `block in lusolve'
   from /opt/local/lib/ruby1.9/1.9.1/bigdecimal/ludcmp.rb:78:in `downto'
   from /opt/local/lib/ruby1.9/1.9.1/bigdecimal/ludcmp.rb:78:in `lusolve'
   from /opt/local/lib/ruby1.9/1.9.1/bigdecimal/newton.rb:60:in `nlsolve'
   from gistfile1.rb:52:in `<main>'

这是因为我的程序中出现错误，还是由于 BigDecimal 库？

#!/opt/local/bin/ruby1.9
require 'bigdecimal'
require 'bigdecimal/newton'
require 'time'
include Newton

transactions = [
  {amount:BigDecimal::new("-1000"), date:Time.parse("2011-01-01")},
  {amount:BigDecimal::new("100"), date:Time.parse("2011-02-11")},
  {amount:BigDecimal::new("100"), date:Time.parse("2011-03-11")},
  {amount:BigDecimal::new("100"), date:Time.parse("2011-04-11")},
  {amount:BigDecimal::new("100"), date:Time.parse("2011-05-11")},
  {amount:BigDecimal::new("100"), date:Time.parse("2011-06-11")},
  {amount:BigDecimal::new("100"), date:Time.parse("2011-07-11")},
  {amount:BigDecimal::new("100"), date:Time.parse("2011-08-11")},
  {amount:BigDecimal::new("100"), date:Time.parse("2011-09-11")},
  {amount:BigDecimal::new("100"), date:Time.parse("2011-10-11")},
  {amount:BigDecimal::new("100"), date:Time.parse("2011-11-11")},
  {amount:BigDecimal::new("100"), date:Time.parse("2011-12-11")},
  {amount:BigDecimal::new("100"), date:Time.parse("2012-01-11")},
]

class Function
  values ={ eps: "1.0e-16", one: "1.0", two: "2.0", ten: "10.0", zero: "0.0" }

  values.each do |key, value|
    define_method key do
      BigDecimal::new value
    end
  end

  def initialize(transactions)
    @transactions = transactions
  end

  def values(x)
    start = @transactions[0][:date]
    value = []

    value << @transactions.reduce(0) do |sum, t|
      pwr = (t[:date] - start) / 365
      r = t[:amount] / (1.0 + x[0]) ** pwr
      sum + r
    end

    value
  end
end

f = Function.new(transactions)
x = [f.zero]
n = nlsolve(f,x)
puts x

任何调试帮助将不胜感激。

谢谢！

Answer 1

问题是您在计算中将 BigDecimals 与“正常”数字混合在一起。
尝试将 values 替换为：

values = { 
  eps: BigDecimal.new("1.0e-16"), 
  one: BigDecimal.new("1.0"), 
  two: BigDecimal.new("2.0"), 
  ten: BigDecimal.new("10.0"), 
  zero: BigDecimal.new("0.0") 
}

并将 def value(x) 中的块替换为

value << @transactions.reduce(0) do |sum, t|
  pwr = (t[:date] - start) / BigDecimal.new("365.0")
  r = t[:amount] / (BigDecimal.new("1.0") + x[0]) ** pwr
  sum + r
end