Newenda
66Messages postés
1 décembre 2008Date d'inscription
9 avril 2009 à 11:04
ok j'en ai enlevé la moitié
-----------------------------------------------
--------------DEBUT PROGRAMME---------
-----------------------------------------------
for i = an_deb:an_fin;
for j = mois_deb:mois_fin;
for k = jour_deb:jour_fin;
for l = heure_deb:heure_fin;
try
dossierCourant = fullfile(num2str(i,'%02d'),num2str(j,'%02d'),num2str(k,'%02d'),num2str(l,'%02d'));
fichier = dir(fullfile(Directory_source,dossierCourant,'\*.pro'));
%%%%%%%%%%%%%% Boucle: Lecture de tous les fichiers du dossier %%%%%%%%%%%
for p=1:size(fichier,1);
fid=fopen(fullfile(Directory_source,dossierCourant,fichier(p).name),'r');
disp(fid)
%%%%%%%%%%%%%% Lecture des Données voulues %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
fseek(fid,312,'bof');
duree_avant = fread(fid,1,'single');
[...]
%%%%%%%%%%%%%% Lecture des voies 1 et voies 2 %%%%%%%%%%%%%%%%%%%%%%%%%%%
fseek(fid,premiere,'bof');
voie1 = fread(fid,nombre_echant,'int16');
size(voie1);
fseek(fid,premiere+decalage,'bof');
voie2 = fread(fid,nombre_echant,'int16');
size(voie2);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Corrélation entre voie 1 et voie 2 %%%%%%%%%
Correl = corr2(voie1, voie2);
%%%%%%%%%%%%%%%%%%%%CONDITION%%%%%%%%%%%%%%%%%%%%
if ( Correl > -1 & Correl < 0.0001 );
compteur_corr_bruit_neg = compteur_corr_bruit_neg + 1
fid = fopen('fichiers negatifs.txt','at');
fprintf(fid,'%s\n',fichier(p).name);
fclose(fid);
elseif( Correl > 0.000001 & Correl < 1 );
compteur_corr_95_1_neg = compteur_corr_95_1_neg + 1
fid = fopen('fichiers positifs.txt','at');
fprintf(fid,'%s\n',fichier(p).name);
fclose(fid)
end
if ( (abs(Correl) > 0.5) & (abs(Correl) < 0.6) );
fid = fopen('fichiers_correl_05_06.txt','at');
fprintf(fid,'%s\n',fichier(p).name);
fclose(fid);
compteur_correl_05_a_06 = compteur_correl_05_a_06 +1 ;
elseif ( (abs(Correl) > 0.9) & (abs(Correl) < 0.95) );
fid = fopen('fichiers_correl_09_95.txt','at');
fprintf(fid,'%s\n',fichier(p).name);
fclose(fid);
compteur_correl_09_a_095 = compteur_correl_09_a_095 +1 ;
elseif ( (abs(Correl) > 0.95) & (abs(Correl) < 1) );
if ( (abs(Correl) > 0.95) & (abs(Correl) < 0.97) );
compteur_correl_95_a_97 = compteur_correl_95_a_97 +1 ;
fid = fopen('fichiers_correl_95_97.txt','at');
fprintf(fid,'%s\n',fichier(p).name);
fclose(fid);
elseif ( (abs(Correl) > 0.97) & (abs(Correl) < 1) );
compteur_correl_97_a_1 = compteur_correl_97_a_1 +1 ;
end
end
fclose(fid) ;
end; %end de for
end; % de try
end
end
end
end
-----------------------------------------------
--------------FIN PROGRAMME--------------
-----------------------------------------------
Merci