function intercov(x,y,L,s);

% intercov(x,y,L,s);
% calcul de l'intercovariance Cxy
% x et y signaux
% L nombre de valeurs représentées, donc entre -L et L
% s paramètre optionnel. 'norm' si normalisée

N = length(x);

v = xcov(x,y,'biased');
v = flipud(v);
if nargin==4,
    sx = sum((x-mean(x)).*(x-mean(x)))/N;
    sy = sum((y-mean(y)).*(y-mean(y)))/N;
    v = v/sqrt(sx*sy);
end

subplot(311)
plot(x)
title('signal 1')
subplot(312)
plot(y)
title('signal 2')
subplot(313)
stem(-L:L,v(N-L:N+L),'fill');
set(gca,'Ylim',[-1 1])
if nargin==4,
    title('intercovariance normalisee')
else
    title('intercovariance')
end
subplot(111)

[bid,c] = max(abs(v));
disp(cat(2,'indice maximum intercovariance: ',num2str(c-N)))
disp(cat(2,'maximum intercovariance: ',num2str(v(c),'%6.2f')))