Java UnivariateInterpolator类代码示例

import org.apache.commons.math3.analysis.interpolation.UnivariateInterpolator; //导入依赖的package包/类
public IDecimalAggregate<E> cubicBSplineInterpolation(final InterpolationConfig config) {
    if (values.isEmpty()) {
        return DummyDecimalAggregate.getInstance();
    }

    if (values.size() < 4) {
        return parent;
    }

    final Pair<List<Double>, List<Double>> pair = fillInterpolationPoints(config, null);
    final List<Double> xval = pair.getFirst();
    final List<Double> yval = pair.getSecond();

    final ControlPath cp = new ControlPath();
    for (int i = 0; i < xval.size(); i++) {
        cp.addPoint(PointFactory.create(xval.get(i), yval.get(i)));
    }
    final GroupIterator gi = new GroupIterator("0:n-1", cp.numPoints());
    final CubicBSpline curve = new CubicBSpline(cp, gi);
    curve.setInterpolateEndpoints(true);
    calculateCurve(xval, yval, curve);

    final UnivariateInterpolator interpolator = new SplineInterpolator();
    return interpolate(config, xval, yval, interpolator);
}
 

import org.apache.commons.math3.analysis.interpolation.UnivariateInterpolator; //导入依赖的package包/类
public IDecimalAggregate<E> bezierCurveInterpolation(final InterpolationConfig config) {
    if (values.isEmpty()) {
        return DummyDecimalAggregate.getInstance();
    }

    final Pair<List<Double>, List<Double>> pair = fillInterpolationPoints(config, BEZIER_CURVE_MAX_SIZE);
    final List<Double> xval = pair.getFirst();
    final List<Double> yval = pair.getSecond();

    final ControlPath cp = new ControlPath();
    for (int i = 0; i < xval.size(); i++) {
        cp.addPoint(PointFactory.create(xval.get(i), yval.get(i)));
    }
    final GroupIterator gi = new GroupIterator("0:n-1", cp.numPoints());
    final BezierCurve curve = new BezierCurve(cp, gi);
    calculateCurve(xval, yval, curve);

    final UnivariateInterpolator interpolator = new SplineInterpolator();
    return interpolate(config, xval, yval, interpolator);
}
 

import org.apache.commons.math3.analysis.interpolation.UnivariateInterpolator; //导入依赖的package包/类
public IDecimalAggregate<E> bSplineInterpolation(final BSplineInterpolationConfig config) {
    if (values.isEmpty()) {
        return DummyDecimalAggregate.getInstance();
    }

    final Pair<List<Double>, List<Double>> pair = fillInterpolationPoints(config, null);
    final List<Double> xval = pair.getFirst();
    final List<Double> yval = pair.getSecond();

    final ControlPath cp = new ControlPath();
    for (int i = 0; i < xval.size(); i++) {
        cp.addPoint(PointFactory.create(xval.get(i), yval.get(i)));
    }
    final GroupIterator gi = new GroupIterator("0:n-1", cp.numPoints());
    final BSpline curve = new BSpline(cp, gi);
    curve.setDegree(config.getDegree());
    final int maxDegree = cp.numPoints() - 1;
    if (curve.getDegree() > maxDegree) {
        curve.setDegree(maxDegree);
    }
    calculateCurve(xval, yval, curve);

    final UnivariateInterpolator interpolator = new SplineInterpolator();
    return interpolate(config, xval, yval, interpolator);
}
 

import org.apache.commons.math3.analysis.interpolation.UnivariateInterpolator; //导入依赖的package包/类
public IDecimalAggregate<E> loessInterpolation(final LoessInterpolationConfig config) {
    if (values.isEmpty()) {
        return DummyDecimalAggregate.getInstance();
    }

    if (values.size() < 3) {
        return parent;
    }

    final Pair<List<Double>, List<Double>> pair = fillInterpolationPoints(config, null);
    final List<Double> xval = pair.getFirst();
    final List<Double> yval = pair.getSecond();
    double bandwidth = config.getSmoothness().getValue(PercentScale.RATE).doubleValue();
    if (bandwidth * values.size() < 2) {
        bandwidth = Decimal.TWO.divide(values.size()).doubleValue();
    }
    final UnivariateInterpolator interpolator = new LoessInterpolator(bandwidth,
            LoessInterpolator.DEFAULT_ROBUSTNESS_ITERS);
    return interpolate(config, xval, yval, interpolator);
}
 

import org.apache.commons.math3.analysis.interpolation.UnivariateInterpolator; //导入依赖的package包/类
private UnivariateInterpolator createInterpolator()
{	
switch(this.gcDepthInterpolation)
	{
	case loess: return new LoessInterpolator(0.5,4);
	case neville: return new NevilleInterpolator();
	case difference  : return new DividedDifferenceInterpolator();
	case linear: return new LinearInterpolator();
	case spline : return new SplineInterpolator();
	case identity : return new UnivariateInterpolator()
			{
			@Override
			public UnivariateFunction interpolate(double[] arg0, double[] arg1)
					throws MathIllegalArgumentException, DimensionMismatchException {
				return new Identity();
				}
			};
	default: throw new IllegalStateException("Not implemented");
	}
}
 

import org.apache.commons.math3.analysis.interpolation.UnivariateInterpolator; //导入依赖的package包/类
private IDecimalAggregate<E> interpolate(final InterpolationConfig config, final List<Double> xval,
        final List<Double> yval, final UnivariateInterpolator interpolator) {
    final UnivariateFunction interpolated = interpolator.interpolate(Doubles.toArray(xval), Doubles.toArray(yval));
    final List<E> interpolatedValues = new ArrayList<E>();
    interpolateAndMaybeReverseMultiplier(config, interpolated, interpolatedValues);
    Assertions.assertThat(interpolatedValues).hasSameSizeAs(values);
    return new DecimalAggregate<E>(interpolatedValues, converter);
}
 

import org.apache.commons.math3.analysis.interpolation.UnivariateInterpolator; //导入依赖的package包/类
/**
  * Create a {@link FromFileTimeseriesParameter} object using the specified
  * file.<br><br>
  * 
  * See also {@link FromFileTimeseriesParameter} and {@link Parameter}.
  * 
  * @param filename
  *        The name of the file from which to read timeseries data.
  * @param parameterName
  *        The name of the resulting {@link Parameter}. It is acceptable
  *        for this argument to be the empty string.
  * @param interpolator
  *        An instance of a {@link UnivariateInterpolator}. This interpolator
  *        is used to convert timeseries keypoints into a queryable function
  *        {@code f(T)} for all {@code T} in the domain of the interpolation
  *        keys. Possible examples, among others, include the {@link SplineInterpolator}
  *        and the {@link LinearInterpolator}.
  * 
  * @throws IOException
  *         This method raises {@link IOException} if, for any reason, the stated
  *         file cannot be processed. Reasons for such failure include: (a) the
  *         file does not exist; (b) the file is inaccessible to the read stream;
  *         (c) the file does not contain any data or contains fewer than two 
  *         records; (d) the file exists and can be accessed but does not have
  *         the required format. See {@link FromFileTimeseriesParameter}.
  */
@Inject
public FromFileTimeseriesParameter(
@Named("FROM_FILE_TIMESERIES_PARAMETER_FILENAME")
   final String filename,
@Named("FROM_FILE_TIMESERIES_PARAMETER_NAME")
   final String parameterName,
@Named("FROM_FILE_TIMESERIES_PARAMETER_INTERPOLATOR")
   final UnivariateInterpolator interpolator
   ) throws IOException {
   super(parameterName);
   Preconditions.checkNotNull(filename);
   final BufferedReader
      reader = new BufferedReader(new FileReader(filename));
   String
      line = null;
   final List<Double>
      x = new ArrayList<Double>(),
      y = new ArrayList<Double>();
   int
      numInterpolationPoints = 0;
   try {
      while((line = reader.readLine()) != null) {
         if(line.isEmpty()) continue;
         final String[]
            entries = line.split("[\\s:;,\t]+");
         if(entries.length != 2)
            throw new IOException(
               getClass().getSimpleName() + ": serialized timeseries are expected to "
             + "consist of exactly two data columns. The input data: " + entries + " does "
             + "not conform to this template. The data in file " + filename + " could not "
             + "be parsed."
             );
         x.add(Double.parseDouble(entries[0]));
         y.add(Double.parseDouble(entries[1]));
         ++numInterpolationPoints;
      }
   }
   catch(final NumberFormatException e) {
      throw new IOException(e.getMessage());
   }
   finally {
      reader.close();
   }
   
   if(numInterpolationPoints <= 1)
      throw new IllegalArgumentException(
         getClass().getSimpleName() + ": file " + filename + " exists and has valid format,"
       + " but contains fewer than 2 records. At least two interpolation points are required."
         );
   final UnivariateFunction function =
      interpolator.interpolate(ArrayUtil.toPrimitive(x), ArrayUtil.toPrimitive(y));
   this.interpolation = function;
   this.bounds = new SimpleDomainBounds(x.get(0), x.get(x.size()-1));
   this.outOfRangeEvaluation = new SimpleDomainBounds(y.get(0), y.get(y.size() - 1));
}
 

import org.apache.commons.math3.analysis.interpolation.UnivariateInterpolator; //导入依赖的package包/类
@Override
protected void addBindings() {
   bind(UnivariateInterpolator.class)
      .annotatedWith(Names.named(getScopeString()))
      .to(SplineInterpolator.class);
}
 

import org.apache.commons.math3.analysis.interpolation.UnivariateInterpolator; //导入依赖的package包/类
@Override
protected void addBindings() {
   bind(UnivariateInterpolator.class)
      .annotatedWith(Names.named(getScopeString()))
      .to(FloorInterpolator.class);
}
 

import org.apache.commons.math3.analysis.interpolation.UnivariateInterpolator; //导入依赖的package包/类
@Override
protected final void setRequiredBindings() {
   requireBinding(Key.get(UnivariateInterpolator.class, Names.named(getScopeString())));
   expose(Key.get(UnivariateInterpolator.class, Names.named(getScopeString())));
}
 

import org.apache.commons.math3.analysis.interpolation.UnivariateInterpolator; //导入依赖的package包/类
/**
 * Constructs a polynomial interpolator to wrap the given univariate interpolator.
 * 
 * @param univariantInterpolator univariant interpolator to wrap
 */
public PolynomialInterpolator( UnivariateInterpolator univariantInterpolator ) {
    _univariateInterpolator = univariantInterpolator;
}