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
| function res=inters(M,T,dx,dy)
res=zeros(size(T));
switch dy
case 1
res=T;
case 2
for n=1:length(T)-dy
v=intersect(T(n,:),T(n+1,:));
res(n,1:length(v))=v;
end
case 3
for n=1:length(T)-dy
v=intersect(T(n,:),intersect(T(n+1,:),T(n+2,:)));
res(n,1:length(v))=v;
end
case 4
for n=1:length(T)-dy
v=intersect(T(n,:),intersect(T(n+1,:),intersect(T(n+2,:),T(n+3,:))));
res(n,1:length(v))=v;
end
case 5
for n=1:length(T)-dy
v=intersect(T(n,:),intersect(T(n+1,:),intersect(T(n+2,:),intersect(T(n+3,:),T(n+4,:)))));
res(n,1:length(v))=v;
end
end
end |
Partager