fig2.m

Radial Basis Function

% Figure 2 in the rbf_fs_2 manual page.
clear

% Hermite training data.
[x, y] = get_data('hermite');

% Test data.
test.name = 'hermite';
test.p = 1000;
test.ord = 1;
test.std = 0;
[xt, yt] = get_data(test);

% Configure the method.
conf.scales = [1 0.5 0.2 0.1];
conf.type = 'cauchy';
conf.bias = 1;
conf.msc = 'BIC';

% Run the method.
[c, r, w, info] = rbf_fs_2(x, y, conf);

% Now do the test set predictions.
Ht = rbf_dm(xt, c, r, info.dmc);
ft = Ht * w;

% Get figure.
fig = get_fig('Figure 2');

% Plot.
hold off
plot(xt, yt, 'k--')
hold on
plot(xt, ft, 'r-', 'LineWidth', 2)

% Configure plot.
set(gca, 'FontSize', 16)
set(gca, 'Position', [0.1 0.15 0.85 0.8])
set(gca, 'XLim', [-4 4])
set(gca, 'YLim', [0 3])
set(gca, 'XTick', -4:2:4)
set(gca, 'YTick', 0:3)
xlabel('x', 'FontSize', 16)
ylabel('y', 'FontSize', 16, 'Rotation', 0)
legend('target', 'prediction')

% Save postscript.
print -depsc fig2