gp_plot.m

王小平《遗传算法——理论、应用与软件实现》随书光盘

% For generating residual error plots for GP journal paper.
%

% Need to edit ver_ga
ver_ga
prob_ga

% Which model ?
mod_num=2;

num_points=length(y_v);

figure(1)
hold off
clg

subplot(211)

plot(y_v,'w')
hold on
plot(Y_P_v(:,I_tol(I_suc(mod_num))),'w:')

%legend('w','Actual','w:','Predicted')

xlabel('Data Points')
ylabel('Composition')


subplot(212)

plot(y_v-Y_P_v(:,I_tol(I_suc(mod_num))),'w')
hold on
mn=mean(y_v-Y_P_v(:,I_tol(I_suc(mod_num))));
sd=std(y_v-Y_P_v(:,I_tol(I_suc(mod_num))));
plot([0 num_points],[0 0],'w--')
plot([0 num_points],[mn+3*sd mn+3*sd],'w--')
plot([0 num_points],[mn-3*sd mn-3*sd],'w--')

xlabel('Data Points')
ylabel('Residual Error')