net.finmath.montecarlo.interestrate.models.covariance.ExponentialDecayLocalVolatilityModel Maven / Gradle / Ivy

Show more of this group Show more artifacts with this name
Show all versions of finmath-lib Show documentation
finmath lib is a Mathematical Finance Library in Java. It provides algorithms and methodologies related to mathematical finance.
The newest version!
/*
 * (c) Copyright Christian P. Fries, Germany. Contact: [email protected].
 *
 * Created on 26.05.2013
 */
package net.finmath.montecarlo.interestrate.models.covariance;

import java.util.Map;

import net.finmath.exception.CalculationException;
import net.finmath.montecarlo.RandomVariableFactory;
import net.finmath.stochastic.RandomVariable;
import net.finmath.stochastic.Scalar;

/**
 * Exponential decay model build on top of a given covariance model.
 *
 * The model constructed for the i-th factor loading is
 * (L_i(t) + d) F_i(t)
 * where d is the displacement and L_i is
 * the realization of the i-th component of the stochastic process and
 * F_i is the factor loading from the given covariance model.
 *
 * The parameter of this model is a joint parameter vector, consisting
 * of the parameter vector of the given base covariance model and
 * appending the displacement parameter at the end.
 *
 * If this model is not calibrateable, its parameter vector is that of the
 * covariance model, i.e., only the displacement parameter will be not
 * part of the calibration.
 *
 * @author Christian Fries
 * @version 1.0
 */
public class ExponentialDecayLocalVolatilityModel extends AbstractLIBORCovarianceModelParametric {

	private static final long serialVersionUID = 4522227972747028512L;
	private final RandomVariableFactory randomVariableFactory;
	private final AbstractLIBORCovarianceModelParametric covarianceModel;
	private final RandomVariable decay;

	private boolean isCalibrateable = false;

	/**
	 * Exponential decay model build on top of a standard covariance model.
	 *
	 * The model constructed for the i-th factor loading is
	 * exp(- a t) F_i(t)
	 * where a is the decay parameter and
	 * F_i is the factor loading from the given covariance model.
	 *
	 * The parameter of this model is a joint parameter vector, consisting
	 * of the parameter vector of the given base covariance model and
	 * appending the decay parameter at the end.
	 *
	 * If this model is not calibrateable, its parameter vector is that of the
	 * underlying covariance model, i.e., only the decay parameter will be not
	 * part of the calibration.
	 *
	 * @param randomVariableFactory A random variable factory (used when cloning with modifed parameters).
	 * @param covarianceModel The given covariance model specifying the factor loadings F.
	 * @param decay The decay a.
	 * @param isCalibrateable If true, the parameter a is a free parameter. Note that the covariance model may have its own parameter calibration settings.
	 */
	public ExponentialDecayLocalVolatilityModel(final RandomVariableFactory randomVariableFactory, final AbstractLIBORCovarianceModelParametric covarianceModel, final RandomVariable decay, final boolean isCalibrateable) {
		super(covarianceModel.getTimeDiscretization(), covarianceModel.getLiborPeriodDiscretization(), covarianceModel.getNumberOfFactors());
		this.randomVariableFactory = randomVariableFactory;
		this.covarianceModel	= covarianceModel;
		this.decay		= decay;
		this.isCalibrateable	= isCalibrateable;
	}

	/**
	 * Exponential decay model build on top of a standard covariance model.
	 *
	 * The model constructed for the i-th factor loading is
	 * exp(- a t) F_i(t)
	 * where a is the decay parameter and
	 * F_i is the factor loading from the given covariance model.
	 *
	 * The parameter of this model is a joint parameter vector, consisting
	 * of the parameter vector of the given base covariance model and
	 * appending the decay parameter at the end.
	 *
	 * If this model is not calibrateable, its parameter vector is that of the
	 * underlying covariance model, i.e., only the decay parameter will be not
	 * part of the calibration.
	 *
	 * @param randomVariableFactory A random variable factory (used for the given parameter and when cloning with modifed parameters).
	 * @param covarianceModel The given covariance model specifying the factor loadings F.
	 * @param decay The displacement a.
	 * @param isCalibrateable If true, the parameter a is a free parameter. Note that the covariance model may have its own parameter calibration settings.
	 */
	public ExponentialDecayLocalVolatilityModel(final RandomVariableFactory randomVariableFactory, final AbstractLIBORCovarianceModelParametric covarianceModel, final double decay, final boolean isCalibrateable) {
		super(covarianceModel.getTimeDiscretization(), covarianceModel.getLiborPeriodDiscretization(), covarianceModel.getNumberOfFactors());
		this.randomVariableFactory = randomVariableFactory;
		this.covarianceModel = covarianceModel;
		this.decay = randomVariableFactory != null ? randomVariableFactory.createRandomVariable(decay) : new Scalar(decay);
		this.isCalibrateable	= isCalibrateable;
	}

	/**
	 * Exponential decay model build on top of a standard covariance model.
	 *
	 * The model constructed for the i-th factor loading is
	 * exp(- a t) F_i(t)
	 * where a is the decay parameter and
	 * F_i is the factor loading from the given covariance model.
	 *
	 * The parameter of this model is a joint parameter vector, consisting
	 * of the parameter vector of the given base covariance model and
	 * appending the decay parameter at the end.
	 *
	 * If this model is not calibrateable, its parameter vector is that of the
	 * underlying covariance model, i.e., only the decay parameter will be not
	 * part of the calibration.
	 *
	 * @param covarianceModel The given covariance model specifying the factor loadings F.
	 * @param decay The displacement a.
	 * @param isCalibrateable If true, the parameter a is a free parameter. Note that the covariance model may have its own parameter calibration settings.
	 */
	public ExponentialDecayLocalVolatilityModel(final AbstractLIBORCovarianceModelParametric covarianceModel, final double decay, final boolean isCalibrateable) {
		this(null, covarianceModel, decay, isCalibrateable);
	}

	@Override
	public Object clone() {
		return new ExponentialDecayLocalVolatilityModel(randomVariableFactory, (AbstractLIBORCovarianceModelParametric) covarianceModel.clone(), decay, isCalibrateable);
	}

	/**
	 * Returns the base covariance model, i.e., the model providing the factor loading F
	 * such that this model's i-th factor loading is
	 * exp(- a t) F_i(t)
	 * where a is the decay parameter and
	 * F_i is the factor loading from the given covariance model.
	 *
	 * @return The base covariance model.
	 */
	public AbstractLIBORCovarianceModelParametric getBaseCovarianceModel() {
		return covarianceModel;
	}

	@Override
	public RandomVariable[] getParameter() {
		if(!isCalibrateable) {
			return covarianceModel.getParameter();
		}

		final RandomVariable[] covarianceParameters = covarianceModel.getParameter();
		if(covarianceParameters == null) {
			return new RandomVariable[] { decay };
		}

		// Append decay to the end of covarianceParameters
		final RandomVariable[] jointParameters = new RandomVariable[covarianceParameters.length+1];
		System.arraycopy(covarianceParameters, 0, jointParameters, 0, covarianceParameters.length);
		jointParameters[covarianceParameters.length] = decay;

		return jointParameters;
	}

	@Override
	public double[] getParameterAsDouble() {
		final RandomVariable[] parameters = getParameter();
		final double[] parametersAsDouble = new double[parameters.length];
		for(int i=0; i dataModified)
			throws CalculationException {
		RandomVariable newDecay = decay;
		boolean isCalibrateable = this.isCalibrateable;
		AbstractLIBORCovarianceModelParametric covarianceModel = this.covarianceModel;

		RandomVariableFactory newRandomVariableFactory = randomVariableFactory;

		if(dataModified != null) {
			if(dataModified.containsKey("randomVariableFactory")) {
				newRandomVariableFactory = (RandomVariableFactory)dataModified.get("randomVariableFactory");
				newDecay = newRandomVariableFactory.createRandomVariable(newDecay.doubleValue());
			}
			if (!dataModified.containsKey("covarianceModel")) {
				covarianceModel = covarianceModel.getCloneWithModifiedData(dataModified);
			}

			// Explicitly passed covarianceModel has priority
			covarianceModel = (AbstractLIBORCovarianceModelParametric)dataModified.getOrDefault("covarianceModel", covarianceModel);
			isCalibrateable = (boolean)dataModified.getOrDefault("isCalibrateable", isCalibrateable);

			if (dataModified.getOrDefault("decay", newDecay) instanceof RandomVariable) {
				newDecay = (RandomVariable) dataModified.getOrDefault("decay", newDecay);
			} else if (newRandomVariableFactory == null) {
				newDecay = new Scalar((double) dataModified.get("decay"));
			} else {
				newDecay = newRandomVariableFactory.createRandomVariable((double) dataModified.get("decay"));
			}
		}

		final AbstractLIBORCovarianceModelParametric newModel = new ExponentialDecayLocalVolatilityModel(newRandomVariableFactory, covarianceModel, newDecay, isCalibrateable);
		return newModel;
	}
}