1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
| %% Premier et dernier intervalle
idx1N3 = x<0.3 & y>=0.9 | 0.7<x& 0.2<y&y<0.3 % Indexage logique
%% Deuxième intervalle
x2 = x(0.3<=x & x<=0.7); % Sélection des valeurs x telles que 0.3<= x <=0.7
y2 = y(0.3<=x & x<=0.7); % Sélection des valeurs y telles que 0.3<= x <=0.7
inf = linspace(1,0.3,length(x2)); % Pente du dessus
inf = interp1(linspace(0.3,0.7,length(x2)),inf,x2); % Ajustement par rapport aux valeurs de x2
sup = linspace(0.9,0.2,length(x2)); % Pente du dessous
sup = interp1(linspace(0.3,0.7,length(x2)),sup,x2); % Ajustement par rapport aux valeurs de x2
idx2 = y2<sup | y2>inf; % Indexage logique
plot(x,y,'+'),hold on
plot(x2(idx2),y2(idx2),'g.')
plot(x(idx1N3),y(idx1N3),'r.')
plot([0 0.3 0.3 0.7 0.7],[0.9 0.9 1 0.3 1],'k:')
plot([0.3 0.3 0.7 1],[0 0.9 0.2 0.2],'k:') |
Partager