1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57
| clc
close
clear
////////////////////////////////////CHOIX DU REPERTOIRE////////////////////////////////////////////
path=uigetdir("C:\Users\salle_laser\Desktop", "Choix du répertoire")
path = path + "/"
//////////////////////////////////////////////////////////////////////////////////////////////////
///////////////////////////////////CHOIX DU FICHIER///////////////////////////////////////////////
fichier=uigetfile('.txt',path,'Choix du fichier contenant les données',%t)
//////////////////////////////////////////////////////////////////////////////////////////////////
///////////////////////////////DETERMINATION DE I_max ET SON ANTECEDENT///////////////////////////
donnees=fscanfMat(fichier)
[m,n]=size(donnees)
[value_Intensity,k]=max(donnees(1:m,2))
value_Lambda=donnees(k,1)
////fprintfMat(path + "résultats.txt",[value_Lambda, value_Intensity],"%f")//Fonctionne bien
//////////////////////////////////////////////////////////////////////////////////////////////////
/////////////////////////////////RECHERCHE DES VALEURS DE Y ET Z/////////////////////////////////
n = fileinfo(fichier)
fd1=mopen(fichier,'rt')
//C = mget(n(1), 'c', fd1)
//mseek(0,fd1,'end');
datas = mfscanf(fd1,'%s %s %s %s')
y=datas(1)
y_not_string = part(y,3:10)
z=datas(2)
z_not_string = part(z,3:10)
//donnees=fscanfMat(fichier)
mclose(fd1)
/////////////////////////////////////////////////////////////////////////////////////////////////
///////////////////////////INSERTION DES VALEURS DE Y, Z, Lambda ET Intensity///////////////////
fd3=mopen(path + 'résultats_avec_valeurs_en_Y_et_Z.txt','at')
mfprintf(fd3,'%s\t%s\t%f\t%f\n',y,z,value_Lambda,value_Intensity)
mclose(fd3)
//donnees_resultats=fscanfMat(path + "résultats_avec_valeurs_en_Y_et_Z.txt")
fd2=mopen(path + 'résultats_sans_valeurs_en_Y_et_Z.txt','at')
mfprintf(fd2,'%f\t%f\n',value_Lambda,value_Intensity)
mclose(fd2)
//[m2,n2]=size(donnees_resultats)
///////////////////////////////////////////////////////////////////////////////////////////////
///////////////////////////////AFFICHAGE GRAPHIQUE/////////////////////////////////////////////
//plot(donnees_resultats(1:m2,1),donnees_resultats(1:m2,n2),'*')
Sgrayplot(y_not_string,z_not_string,941.0000000)
//xtitle("I = f (périodicité spatiale)","$\text{\lambda (nm)}$","$\text{Intensity (a.u.)}$")**
mclose("all")
/////////////////////////////////////////////////////////////////////////////////////////////// |
Partager