Kernel functions

Radial Basis Function (RBF)

Given

 0 

import numpy as np
from sklearn.gaussian_process.kernels import RBF

np.random.seed(42)

xlim = (-4, 4)
X = np.expand_dims(np.linspace(*xlim, num=75), 1)
zero = np.array([0.](0.))

amplitude = 1.0
Σ1 = (amplitude**2) * RBF(length_scale=1)(X)
Σ2 = (amplitude**2) * RBF(length_scale=0.5)(X)
amplitude = 0.5
Σ3 = (amplitude**2) * RBF(length_scale=1)(X)