Commit d62ddeed authored by Fnadi Mohamed's avatar Fnadi Mohamed
Browse files

Coverge rates between simulated and experiments data

parent ad4aa6cb
This source diff could not be displayed because it is too large. You can view the blob instead.
This source diff could not be displayed because it is too large. You can view the blob instead.
This source diff could not be displayed because it is too large. You can view the blob instead.
This source diff could not be displayed because it is too large. You can view the blob instead.
This source diff could not be displayed because it is too large. You can view the blob instead.
This source diff could not be displayed because it is too large. You can view the blob instead.
This source diff could not be displayed because it is too large. You can view the blob instead.
This source diff could not be displayed because it is too large. You can view the blob instead.
This source diff could not be displayed because it is too large. You can view the blob instead.
This source diff could not be displayed because it is too large. You can view the blob instead.
This source diff could not be displayed because it is too large. You can view the blob instead.
This source diff could not be displayed because it is too large. You can view the blob instead.
This source diff could not be displayed because it is too large. You can view the blob instead.
This source diff could not be displayed because it is too large. You can view the blob instead.
This source diff could not be displayed because it is too large. You can view the blob instead.
This source diff could not be displayed because it is too large. You can view the blob instead.
%%avec export j0 (les fichiers doivent être modifiés avec suppression [ et , )
%50
%A=load('Coordonees1_50.txt');
%tube=load('Taux_de_couverture_/IA_q1_50_mod.txt');
%tube=load('Taux_de_couverture_/LS_q1_50_mod.txt');
%100
%A=load('Coordonees1_100.txt');
%tube=load('Taux_de_couverture_/IA_q2_100_mod.txt');
%tube=load('Taux_de_couverture_/LS_q2_100_mod.txt');
%01
A=load('Coordonees_0_1.txt');
tube=load('Taux_de_couverture_/IA_q2_0_1_mod.txt');
%tube=load('Taux_de_couverture_/LS_q2_0_1_mod.txt');
%017
%A=load('Coordonees1_0_17.txt');
%tube=load('Taux_de_couverture_/IA_q2_0_17_mod.txt');
%tube=load('Taux_de_couverture_/LS_q2_0_17_mod.txt');
dec=A(1,1);
for i=1:size(A,1)
A(i,1)=A(i,1)-dec;
A(i,2)=A(i,2)*(3.14159/180.0);
A(i,3)=A(i,3)*(3.14159/180.0);
endfor
eps=0;%1e-6;
k=1;
inclu=0;
time_inclu=0;
for j=1:size(A,1)
k=1;
while (k<size(tube,1))
if ((tube(k,1) <= A(j,1)+eps) && (A(j,1) -eps <= tube(k,2)))
%minb=tube(k,7);
%maxb=tube(k,8);
minb = min([tube(k,7), tube(max(k-1,1),7), tube(max(k-2,1),7), tube(max(k-3,1),7),tube(min(k+1,size(tube,1)),7),...
tube(min(k+2,size(tube,1)),7), tube(min(k+3,size(tube,1)),7),tube(k,8), tube(max(k-2,1),8),tube(max(k-3,1),8), ...
tube(max(k-1,1),8), tube(min(k+1,size(tube,1)),8), tube(min(k+2,size(tube,1)),8),tube(min(k+3,size(tube,1)),8)]);
maxb = max([tube(k,7), tube(max(k-1,1),7), tube(max(k-2,1),7), tube(max(k-3,1),7),tube(min(k+1,size(tube,1)),7),...
tube(min(k+2,size(tube,1)),7), tube(min(k+3,size(tube,1)),7),tube(k,8), tube(max(k-2,1),8),tube(max(k-3,1),8), ...
tube(max(k-1,1),8), tube(min(k+1,size(tube,1)),8), tube(min(k+2,size(tube,1)),8),tube(min(k+3,size(tube,1)),8)]);
%if ((tube(k,1) <= A(j,3)+eps) && (A(j,3)-eps <= tube(k,2)))
%disp(minb)
%disp(maxb)
if ((minb <= A(j,3)) && (A(j,3) <= maxb))
inclu++;
%disp(A(j,1))
endif
time_inclu++;
break;
else
k++;
endif
endwhile
endfor
disp(time_inclu)
disp(inclu)
disp('couverture : ')
disp(inclu*100/time_inclu)
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment