% build dynamics model (part 1) %modele C Y=load('ListeDesVolumes.txt'); y_training=zeros(58,1); C=xlsread('points courbes.xls') % on normalise les patients J=1; I=1; ymax_training=Y; ymax_training=zeros(18,1); while (I<=18); % on cherche le max du patient J=1; while (J<=25); if Y(J,I) > ymax_training(I,1) ymax_training(I,1)=Y(J,I); end J=J+1; end % on normalise le patient J=1; while (J<=25); Y(J,I) = Y(J,I) / ymax_training(I,1); J=J+1; end I=I+1; end % on récupère les frames qui nous intéressent J=1; I=1; i=1; while (I<=18); J=C(I,2); while (J<=25); y_training(i,1)=Y(J,I); i=i+1; J=J+1; end I=I+1; end Y_training = y_training; %creation de X_training k=1; X_training=[]; while k<=18 X=[0:1/(25-C(k,2)):1]; X_training = [X_training, X]; k=k+1; end X_training=X_training'; %X_training=repmat(X_training,18,1); save('training.mat','X_training','Y_training')