Changeset 16 for trunk/private

Timestamp:

08/18/09 14:38:44 (15 years ago)

Author:

pinsard

Message:

script_correc : save figures

File:

• trunk/private/script_correc.m (modified) (8 diffs)

trunk/private/script_correc.m

@@ -12 +12 @@
 %
 % ++
+%
+% ``script_correc`` plot ans save several figures :
 % 
+%  - script_correc_thetas.eps : 
+%
+%  - script_correc_thetav.eps : 
+%
+%  - script_correc_dphi.eps : 
+%
+%  - script_correc_rho_670_.eps : 
+%
+%  - script_correc_rho_765_.eps : 
+%
+%  - script_correc_rho_865_.eps : 
+%
+%  - script_correc_tau.eps : 
+%
+%  - script_correc_diagramme_dispersion.eps :
+%
+%  - script_correc_histogramme_erreur.eps :
+%
 % EXAMPLES
 % ========
@@ -39 +59 @@
 % improve description
 %
-% avoid overriding existing output file
+% avoid overriding existing output files
 %
 % make it launchable from everywhere (path and IO directories)
@@ -45 +65 @@
 % Running octave, ``more off`` is compulsary but I don't not yet why
 %
-% save output (sdout and figures) to make comparisons possible
+% save output (sdout and stderr) to make comparisons possible
+%
+% If printer is set to png, one might see following warning using octave :
+% ..
+%
+% gdImageStringFT: Could not find/open font while printing string ... with font Helvetica
+%
+% And no labels are visible on figure outputfiles.
+%
+% For now, find a portable solution (setting GDFONTPATH) is too tough, so 
+% printer is set to eps.
 %
 % EVOLUTIONS
@@ -51 +81 @@
 %
 % $Id$
-% 
+%
+% - fplod 2009-08-18T09:55:12Z aedon.locean-ipsl.upmc.fr (Darwin)
+%
+%   * save figures to make comparisons possible
+%   * replace error on suptitle using octave by a warning
+%
 % - fplod 2009-08-17T12:39:51Z aedon.locean-ipsl.upmc.fr (Darwin)
-% 
+%
 %   * add header
 %   * add octave intersect_octave
@@ -106 +141 @@
 
 %% QUESTION 3 : rien a faire
- 
+
 %% QUESTION 4
 moy=mean(base_app);
@@ -148 +183 @@
  suptitle(stitle{j});
 else
- error('eee : octave can not write supertitle');
-end
+ warning('www : octave can not write supertitle');
+end
+
+% save figure(j)
+printer='eps';
+print_printer=['-d', printer];
+% figure name is based on stitle{j}.
+% transformation is lower and replace ( and ) by _
+if (run_octave == 0)
+ figname=regexprep(lower(stitle{j}), '(', '_');
+else
+ figname=regexprep(lower(stitle{j}), '\(', '_');
+end
+if (run_octave == 0)
+ figname=regexprep(figname, ')', '_');
+else
+ figname=regexprep(figname, '\)', '_');
+end
+fullfilename=['./', mfilename, '_', figname , '.', printer];
+clear figname;
+clear printer;
+print(print_printer,fullfilename);
+clear print_printer;
+clear fullfilename;
 
 end
@@ -173 +230 @@
     [nets(n,i).W1 nets(n,i).W2]=MLPfit(base_val_n(:,1:end-1),base_val_n(:,end),base_app_n(:,1:end-1), ...
                                        base_app_n(:,end),nets(n,i).W1,nets(n,i).W2,'tah','lin',[1500,1e-6],[100,1]);
-        
-    
+
+
 %QUESTION 6
-    
+
     %Restitution du reseau
   tau_nn=MLPval(base_val_n(:,1:end-1),nets(n,i).W1,nets(n,i).W2);
-  
+
   %Performance du reseau
   delta=decenred(tau_nn,moy(end),et(end)) - decenred(base_val_n(:,end),moy(end),et(end));
   perf_opti=sqrt((1./length(delta))*sum(delta.^2));
-  
+
   fprintf(1,'\nRMS finale : %3.2e\n',perf_opti);
 
@@ -193 +250 @@
   xlabel('TAU(vrai)')
   ylabel('TAU(estime)');
-    
+
+  % save figure
+  printer='eps';
+  print_printer=['-d', printer];
+  fullfilename=['./', mfilename, '_diagramme_dispersion', '.', printer];
+  clear printer;
+  print(print_printer,fullfilename);
+  clear print_printer;
+  clear fullfilename;
+
   figure
   hist(delta,50);
   xlabel('TAU(estime)-TAU(vrai)');
   ylabel('Nombre');
+
+  % save figure
+  printer='eps';
+  print_printer=['-d', printer];
+  fullfilename=['./', mfilename, '_histogramme_erreur', '.', printer];
+  clear printer;
+  print(print_printer,fullfilename);
+  clear print_printer;
+  clear fullfilename;