Final comparison of techniques

To complete the comparison of the techniques, add the lines

[Wcg1, Ecg1]=pccg(W,50,1); %# 50 iters from W of cg versions 1, 2, 3 ...
[Wcg2, Ecg2]=pccg(W,50,2);
[Wcg3, Ecg3]=pccg(W,50,3);
immediately after the call to pclm(W,50) in your current version of ``sch.m'' and replace the current plotting call to semilogy(...) with
semilogy([1:length(Elm)],Elm-43.3371147782040,'-;lm;', ...
         [1:length(Epclm)],Epclm-43.3371147782040,'-;pclm;', ...
         [1:length(Ecg1)],Ecg1-43.3371147782040,'-;pccg-FR;', ...
         [1:length(Ecg2)],Ecg2-43.3371147782040,'-;pccg-PR;', ...
         [1:length(Ecg3)],Ecg3-43.3371147782040,'-;pccg-HS;');
%# Set W to best result for plotting
W=Wcg3;



Tomas Arias 2006-04-26