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 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151
| %***********************************
%INIT
%***********************************
importfile signalbsb1.txt;
importfile signalbsb2.txt;
importfile signalbsb3.txt;
importfile signalbsb4.txt;
importfile signalbsba.txt;
importfile signalbsbb.txt;
norm_bord1=0;
norm_bord2=0;
norm_bord3=0;
norm_bord4=0;
norm_borda=0;
norm_bordb=0;
norm_bouc1=0;
norm_bouc2=0;
norm_bouc3=0;
norm_bouc4=0;
norm_bouca=0;
norm_boucb=0;
%***********************************
%***********************************
%***********************************
%SCRIPT
%***********************************
%***********************************
%MISE EN FORME DES DONNEES
%***********************************
image1 = reshape(signalbsb1,445,[]).';
image2 = reshape(signalbsb2,445,[]).';
image3 = reshape(signalbsb3,445,[]).';
image4 = reshape(signalbsb4,445,[]).';
imagea = reshape(signalbsba,445,[]).';
imageb = reshape(signalbsbb,445,[]).';
mat_bord1=image1(600:800,90:355);
mat_bord2=image2(600:800,90:355);
mat_bord3=image3(600:800,90:355);
mat_bord4=image4(600:800,90:355);
mat_borda=imagea(600:800,90:355);
mat_bordb=imageb(600:800,90:355);
mat_bouc1=image1(900:1100,90:355);
mat_bouc2=image2(900:1100,90:355);
mat_bouc3=image3(900:1100,90:355);
mat_bouc4=image4(900:1100,90:355);
mat_bouca=imagea(900:1100,90:355);
mat_boucb=imageb(900:1100,90:355);
for i=1:266
vect_bord1((i-1)*201+1:i*201,1)=mat_bord1(:,i);
vect_bord2((i-1)*201+1:i*201,1)=mat_bord2(:,i);
vect_bord3((i-1)*201+1:i*201,1)=mat_bord3(:,i);
vect_bord4((i-1)*201+1:i*201,1)=mat_bord4(:,i);
vect_borda((i-1)*201+1:i*201,1)=mat_borda(:,i);
vect_bordb((i-1)*201+1:i*201,1)=mat_bordb(:,i);
vect_bouc1((i-1)*201+1:i*201,1)=mat_bouc1(:,i);
vect_bouc2((i-1)*201+1:i*201,1)=mat_bouc2(:,i);
vect_bouc3((i-1)*201+1:i*201,1)=mat_bouc3(:,i);
vect_bouc4((i-1)*201+1:i*201,1)=mat_bouc4(:,i);
vect_bouca((i-1)*201+1:i*201,1)=mat_bouca(:,i);
vect_boucb((i-1)*201+1:i*201,1)=mat_boucb(:,i);
end
%***********************************
%***********************************
%***********************************
%NORME
%***********************************
norm_bord1=sqrt(sum(vect_bord1.*vect_bord1));
norm_bord2=sqrt(sum(vect_bord2.*vect_bord2));
norm_bord3=sqrt(sum(vect_bord3.*vect_bord3));
norm_bord4=sqrt(sum(vect_bord4.*vect_bord4));
norm_borda=sqrt(sum(vect_borda.*vect_borda));
norm_bordb=sqrt(sum(vect_bordb.*vect_bordb));
norm_bouc1=sqrt(sum(vect_bouc1.*vect_bouc1));
norm_bouc2=sqrt(sum(vect_bouc2.*vect_bouc2));
norm_bouc3=sqrt(sum(vect_bouc3.*vect_bouc3));
norm_bouc4=sqrt(sum(vect_bouc4.*vect_bouc4));
norm_bouca=sqrt(sum(vect_bouca.*vect_bouca));
norm_boucb=sqrt(sum(vect_boucb.*vect_boucb));
%***********************************
%***********************************
%***********************************
%PRODUIT SCALAIRE
%***********************************
ps_bord2=norm_bord2+(vect_bord1.*vect_bord2);
ps_bord3=norm_bord3+(vect_bord1.*vect_bord3);
ps_bord4=norm_bord4+(vect_bord1.*vect_bord4);
ps_borda=norm_borda+(vect_bord1.*vect_borda);
ps_bordb=norm_bordb+(vect_bord1.*vect_bordb);
ps_bouc2=norm_bouc2+(vect_bouc1.*vect_bouc2);
ps_bouc3=norm_bouc3+(vect_bouc1.*vect_bouc3);
ps_bouc4=norm_bouc4+(vect_bouc1.*vect_bouc4);
ps_bouca=norm_bouca+(vect_bouc1.*vect_bouca);
ps_boucb=norm_boucb+(vect_bouc1.*vect_boucb);
%***********************************
%***********************************
%***********************************
%COEFFICIENT DE CORRELATION
%***********************************
% la bouteille 1 est suppoée bonne et donc prise comme reference %
cor_bord2=ps_bord2/(norm_bord2*norm_bord1);
cor_bord3=ps_bord3/(norm_bord3*norm_bord1);
cor_bord4=ps_bord4/(norm_bord4*norm_bord1);
cor_borda=ps_borda/(norm_borda*norm_bord1);
cor_bordb=ps_bordb/(norm_bordb*norm_bord1);
cor_bouc2=ps_bouc2/(norm_bouc2*norm_bouc1);
cor_bouc3=ps_bouc3/(norm_bouc3*norm_bouc1);
cor_bouc4=ps_bouc4/(norm_bouc4*norm_bouc1);
cor_bouca=ps_bouca/(norm_bouca*norm_bouc1);
cor_boucb=ps_boucb/(norm_boucb*norm_bouc1);
%***********************************
%***********************************
%***********************************
%AFFICHAGE
%***********************************
colormap(gray);
image(mat_bord1*256);
% image(mat_bord2*256);
% image(mat_bord3*256);
% image(mat_bord4*256);
% image(mat_borda*256);
% image(mat_bordb*256);
%
% image(mat_bouc1*256);
% image(mat_bouc2*256);
% image(mat_bouc3*256);
% image(mat_bouc4*256);
% image(mat_bouca*256);
% image(mat_boucb*256);
%***********************************
%***********************************
%***********************************
%*********************************** |
Partager