/*
* Encog(tm) Core v3.4 - Java Version
* http://www.heatonresearch.com/encog/
* https://github.com/encog/encog-java-core
* Copyright 2008-2016 Heaton Research, Inc.
*
* Licensed under the Apache License, Version 2.0 (the "License");
* you may not use this file except in compliance with the License.
* You may obtain a copy of the License at
*
* http://www.apache.org/licenses/LICENSE-2.0
*
* Unless required by applicable law or agreed to in writing, software
* distributed under the License is distributed on an "AS IS" BASIS,
* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
* See the License for the specific language governing permissions and
* limitations under the License.
*
* For more information on Heaton Research copyrights, licenses
* and trademarks visit:
* http://www.heatonresearch.com/copyright
*/
package org.encog.mathutil.rbf;
/**
* Multi-dimensional Mexican Hat, or Ricker wavelet, function.
*
* It is usually only referred to as the "Mexican hat" in the Americas, due to
* cultural association with the "sombrero". In technical nomenclature this
* function is known as the Ricker wavelet, where it is frequently employed to
* model seismic data.
*
* http://en.wikipedia.org/wiki/Mexican_Hat_Function
*
*/
public class MexicanHatFunction extends BasicRBF {
/**
* Serial id.
*/
private static final long serialVersionUID = 1L;
/**
* Construct a single-dimension Mexican hat function with the specified
* peak, centers and widths.
*
* @param peak
* The peak for all dimensions.
* @param center
* The centers for each dimension.
* @param width
* The widths for each dimension.
*/
public MexicanHatFunction(final double center, final double peak,
final double width) {
setCenters(new double[1]);
getCenters()[0] = center;
setPeak(peak);
setWidth(width);
}
/**
* Construct a multi-dimension Mexican hat function with the specified peak,
* centers and widths.
*
* @param peak
* The peak for all dimensions.
* @param center
* The centers for each dimension.
* @param width
* The widths for each dimension.
*/
public MexicanHatFunction(final double peak, final double[] center,
final double width) {
setCenters(center);
setPeak(peak);
setWidth(width);
}
/**
* Create centered at zero, width 0, and peak 0.
* @param dimensions Dimensions.
*/
public MexicanHatFunction(final int dimensions) {
setCenters(new double[dimensions]);
setPeak(1.0);
setWidth(1.0);
}
/**
* Calculate the result from the function.
*
* @param x
* The parameters for the function, one for each dimension.
* @return The result of the function.
*/
public final double calculate(final double[] x) {
final double[] center = getCenters();
final double width = getWidth();
// calculate the "norm", but don't take square root
// don't square because we are just going to square it
double norm = 0;
for (int i = 0; i < center.length; i++) {
norm += Math.pow(x[i] - center[i], 2)/ (2.0 * width * width);
}
// calculate the value
return getPeak() * (1 - norm) * Math.exp(-norm / 2);
}
}