source: trunk/src/sigmoide2.m @ 327

function [result]=sygmoide()

%+
%
% ==========
% sygmoide.m
% ==========
%
% .. function:: sygmoide()
%
% TODO
% =====
%
% name of function vs name of file
%
% corrige la  trace les deux branches de la fonction sigmoïde avec elimination des points suplémentaires
% (refaire l'initialisation des variables)
% automatisation des parametres d'ajustement
%
% gestion  erreur load
%
% EVOLUTIONS
% ==========
% 
% $Id$
%
% $URL$
%
% - fplod 20110705T134725Z aedon.locean-ipsl.upmc.fr (Darwin)
%
%   * correction header (chapter level)
%
% - soclod ?
%
%   * charger les donnees datelai.csv sous format d'une ligne qui
%     contient les jours de l'annee
%   * definition manuelle des parametres d'ajustement pour les deux branches 
%     de la fonction (mc,nc,md,et nd)
%   * definition des indices correspondant à l'amplitude total et la
%     faible valeur du LAI  
%   * description deux branches de la fonction sigmoïde, se basant
%     sur l'evolution du LAI pour boucler sur les numeros de dates 
%
%-

result=[];

%clear all
%close all
% date=load('datelai.csv');
% ifigure=1;
% figure(ifigure);
global PROJECT_ID
% build the datelai filename to be read
fullfilename=[PROJECT_ID  'datelai.csv'];
disp(['iii : opening for reading ', fullfilename]);
% opening file
date=load(fullfilename);
clear fullfilename

nbdate=size(date,2);
% parametres d'ajustement correspondant à la croissance 
mc=30.11;
nc=0.13;
% les valeurs du LAI correspondant a l'amplitude totale et à la
% faible valeur 
a1c=1.34;
a2c=4.10;
% parametres d'ajustement de la branche correspondant à la decroissance
md=37.97;
nd=0.13;
%les valeurs du LAI correspondant a l'amplitude totale et à la
% faible valeur 
a1d=1.18;
a2d=4.25;
%j1=153;
%j2=281;
%j3=233;
%j4=361;
nbdate=45;

% initialisation des varaiables 
ac=zeros(nbdate);
ad=zeros(nbdate);
%for i=1:nbdate;
%indc=find((date == 153) & (date == 281));
%indd=find((date == 233) & (date == 361));
% description de la branche de la croissance de la fonction
% sigmoïde 
for i=20:35;
  ac(i)=((a1c+a2c*(1/(1+exp(mc-nc*date(i))))))/5;
end
% description de la branche de la decroissance de la fonction
% sigmoïde
  for j=30:nbdate;
     ad(j)=(a1d+a2d*(1/(1+exp(-(md-nd*date(j))))))/5;
end

plot(date,ac,'rx');
hold on
plot(date,ad,'b*');
set(gca,'xtick',0:50:361)
xlim([0 361]);
xlabel('jours(2006)');
ylabel('LAI');
title('fonction sigmoïde');
%legend('lai','verdure ajustee','senescence ajustée');

 result=0;
end