cauchyrandomgenerator.java

来自「JGAP（发音"jay-gap"）是一款用Java编写的遗传算法包。提供了基本的」· Java 代码 · 共 104 行

/*
 * This file is part of JGAP.
 *
 * JGAP offers a dual license model containing the LGPL as well as the MPL.
 *
 * For licencing information please see the file license.txt included with JGAP
 * or have a look at the top of class org.jgap.Chromosome which representatively
 * includes the JGAP license policy applicable for any file delivered with JGAP.
 */
package org.jgap.impl;

import java.util.*;

import org.jgap.*;

/**
 * Cauchy probability density function (cumulative distribution function)
 *
 * @author Klaus Meffert
 * @since 1.1
 */
public class CauchyRandomGenerator
    implements RandomGenerator {
  /** String containing the CVS revision. Read out via reflection!*/
  private final static String CVS_REVISION = "$Revision: 1.12 $";

  private double m_scale;

  private double m_location;

  private Random m_rn;

  /**
   * @author Klaus Meffert
   * @since 1.1
   */
  public CauchyRandomGenerator() {
    this(0.0d, 1.0d);
  }

  /**
   * @param a_location double
   * @param a_scale double
   *
   * @author Klaus Meffert
   * @since 1.1
   */
  public CauchyRandomGenerator(final double a_location, final double a_scale) {
    m_location = a_location;
    m_scale = a_scale;
    m_rn = new Random();
  }

  public int nextInt() {
    return Math.min(Integer.MAX_VALUE - 1,
                    (int) Math.round(nextCauchy() * Integer.MAX_VALUE));
  }

  public int nextInt(final int a_ceiling) {
    return Math.min(a_ceiling - 1,
                    (int) Math.round(nextCauchy() * a_ceiling));
  }

  public long nextLong() {
    return Math.min(Long.MAX_VALUE - 1,
                    Math.round(nextCauchy() * Long.MAX_VALUE));
  }

  public double nextDouble() {
    return nextCauchy();
  }

  public float nextFloat() {
    return Math.min(Float.MAX_VALUE - 1,
                    (float) (nextCauchy() * Float.MAX_VALUE));
  }

  public boolean nextBoolean() {
    return nextCauchy() >= 0.5d;
  }

  /**
   * Calculate Cumulative Cauchy distribution function.
   * @return the probability that a stochastic variable x is less than X
   *
   * @author Klaus Meffert
   * @since 1.1
   */
  public double nextCauchy() {
    return 0.5
        + Math.atan( (m_rn.nextDouble() - m_location) / m_scale) / Math.PI;
  }

  /**
   * @return the cauchy standard deviation
   *
   * @author Klaus Meffert
   * @since 1.1
   */
  public double getCauchyStandardDeviation() {
    return m_scale;
  }
}