Changeset 4

Timestamp:

08/12/09 16:55:14 (15 years ago)

Author:

pinsard

Message:

dealing with CR/LF and encoding

Location:

trunk

Files:

• private/make_image_sw.m (modified) (1 diff)
• toolbox/MLPfit.m (modified) (1 diff)
• toolbox/MLPval.m (modified) (1 diff)
• toolbox/cenred.m (modified) (1 diff)
• toolbox/decenred.m (modified) (1 diff)

trunk/private/make_image_sw.m

@@ -13 +13 @@
 Tau_interp=interp2(Lat0,Lon0,Tau,Lat,Lon,'nearest');
 
-%mise au format des autres donnï¿œes
+%mise au format des autres données
 load([indir 'Valid_Pixels_251.dat']);
 Tau_1d=nan*ones(size(Valid_Pixels_251));

trunk/toolbox/MLPfit.m

@@ -460 +460 @@
         %%plotplot(Xa,MLPval(Xa,w1,w2,F1,F2),'b-',X,Yb,'r-')
         %plot(Xa,MLPval(Xa,w1,w2,F1,F2),'b-')
-        %S=sprintf('en rouge fonction utilisŽe pour simuler les donnŽes en bleu courbe apprise par le MLP');
+        %S=sprintf('en rouge fonction utilisée pour simuler les données en bleu courbe apprise par le MLP');
         %title(S);
       

trunk/toolbox/MLPval.m

@@ -1 +1 @@
 function [Y]=MLPval(x,W1,W2,F1,F2)
-
 %MLPval compute output values from a trained MLP
-
 %
-
 %      [Y]=MLPval(X,W1,W2,F1,F2)
-
 %
-
 % X  is the input data
-
 % W1 the weight matrix input -> hidden layer
-
 % W2 weight matrix hidden layer -> ouput
-
 % F1 (default 'tah') activation function of the hidden layer units
-
 % F2 (default 'lin') activation function of the output layer units
-
 %
-
 % Y the output ouf the network at point X
-
 %       
-
 %       21/04/97 S. Canu
-
-
-
 if nargin < 3;
-
    help MLPval
-
    error(sprintf('\n *** MLPval error: invalid call***\n\n\t[Y]=MLPval(X,W1,W2,F1,F2);\n\n'));
-
 end;
-
-
-
 % Check that matrix (X) and matrix (W1) have compatible dimensions
-
-
-
 [n,d] = size(x);
-
 [dp1,m] = size(W1);
-
 [mp1,q] = size(W2);
-
 if (d+1)~=dp1,
-
     error('The number of colum in x plus one must equal the number of rows in W1');
-
 end
-
 if (m+1)~=mp1,
-
     error('The number of colum in W1 plus one must equal the number of rows in W2');
-
 end
-
-
-
 % Complete unspecified name of activation function
-
-
-
 if nargin < 5;
-
    F2='lin';
-
    if nargin < 4;
-
       F1 ='tah';
-
    end;
-
 end;
-
-
-
-
-
 if F1=='tah' & F2=='lin'
-
   a1 = [x ones(n,1)]*W1;
-
   x1 = tanh(a1);
-
   a2 = [x1 ones(n,1)]*W2;
-
   Y =  a2;
-
-
-
 elseif F1=='sig' & F2=='lin'
-
   a1 = [x ones(n,1)]*W1;
-
   x1 = tanh(a1);
-
   a2 = [x1 ones(n,1)]*W2;
-
   Y =  a2;
-
-
-
 elseif F1=='tah' & F2=='tah'
-
   a1 = [x ones(n,1)]*W1;
-
   x1 = phi(a1);
-
   a2 = [x1 ones(n,1)]*W2;
-
   Y =  tanh(a2);
-
-
-
 elseif F1=='sig' & F2=='tah'
-
   a1 = [x ones(n,1)]*W1;
-
   x1 = phi(a1);
-
   a2 = [x1 ones(n,1)]*W2;
-
   Y =  tanh(a2);
-
-
-
 elseif F1=='tah' & F2=='sig'
-
   a1 = [x ones(n,1)]*W1;
-
   x1 = tanh(a1);
-
   a2 = [x1 ones(n,1)]*W2;
-
   Y =  phi(a2);
-
-
-
 elseif F1=='sig' & F2=='sig'
-
   a1 = [x ones(n,1)]*W1;
-
   x1 = phi(a1);
-
   a2 = [x1 ones(n,1)]*W2;
-
   Y =  phi(a2);
   
   
 elseif F1=='tah' & F2=='exp'
-
   a1 = [x ones(n,1)]*W1;
-
   x1 = tanh(a1);
-
   a2 = [x1 ones(n,1)]*W2;
-
   Y =  exp(a2);
-
-
-
 elseif F1=='sig' & F2=='exp'
-
   a1 = [x ones(n,1)]*W1;
-
   x1 = phi(a1);
-
   a2 = [x1 ones(n,1)]*W2;
-
   Y =  exp(a2);
-
-
 else
-
  error('The name of the activation function is incorect')
-
 end
-
-
-
-
-

trunk/toolbox/cenred.m

@@ -1 +1 @@
 function [X,moya,ecaa] = cenred(Z,moya,eca)
-
 % Fonction qui centree par moya et eca 
 %avec la formule X=(2/3)*(Z-moya)/eca
 %les invivudus etant en ligne
-
-
-
 UN = ones(size(Z));
-
 Moy =moya ;
 Std =eca;
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-
-
 Moy = diag(Moy);
-
 Moy = UN * Moy;
-
-
-
 taille=size(Z,1);
-
 Ecart_type = diag(sqrt((taille-1)/taille)*Std);
-
 Ecart_type = UN * Ecart_type;
-
  
  
-
-
-
-
-% Calcul de la matrice X des donnï¿œes centrï¿œes rï¿œduites
-
+% Calcul de la matrice X des données centrées réduites
 X = (2/3).*(Z - Moy)./Ecart_type ;
-
-
-

trunk/toolbox/decenred.m

@@ -1 +1 @@
 function X = decenred(Z,moya,ecaa)
-
 % Fonction pour denormaliser la matrice Z
-
-
 % en lignes
 %load moya;load ecaa
-
 UN = ones(size(Z));
 Moy = moya;
 Moy = diag(Moy);
 Moy = UN * Moy;
-
 taille=size(Z,1);
-
 Ecart_type = diag(sqrt((taille-1)/taille)*ecaa);
 Ecart_type = UN * Ecart_type;
-
  
-% Calcul de la matrice X des donnï¿œes centrï¿œes rï¿œduites
-
+% Calcul de la matrice X des données centrées réduites
 X = (3./2.)* Z .*Ecart_type + Moy;
-
-
-