function a = match(T,W,p, rs)
%MATCH   finds a sampled chirp signal in a longer time signal
% usage s = dchirp ( TW, p )
%    s : sampes of a digital "chirp" signal
%      exp(j(W/T)pi*t^2)  -T/2 <= t <= T/2
%    TW: time-bandwidth product
%    p: sample at p times the Nyquist rate (W)

subplot(2,1,1),plot(abs(rs));
N=round(p*T*W);
n= [0:1:N-1];
D=zeros(1,length(rs)); 
for i = -32:32,
  s= doppler(1,T,W,p,T*W*p,0,(i*pi)/64,0);
  sflip=fliplr(s);

  %g= [zeros(1,N) doppler(8,T,W,p,4*p*W*T,0,0,1)];
  sflip = [zeros(1,length(rs)-length(sflip)) sflip];
  
  SF = fft(sflip);
  RS = fft(rs);
  D= SF.*RS + D;
end

d = ifft(D);
%d= d(N:N2);
%dmovavg=movavg(100,d);
a=d;
  
  
%a=dmovavg;

%for k=1:N2-N+1,
%   c = norm(rs(k:k+N-1));
%   if c==0
%      c=.000000000001;
%      end
%   d(k)= d(k)/c;
%end

%mham=m*hamming(length(m));
%m2=abs(d);

% Testing Bullshit Solution
%a(1:N) = zeros(1,N);
%for k=1:N2-3*N+1
%   a(k+N)=mean(m2(k:k+2*N-1));
%end

%a= [a(N+1:length(a)) .9*ones(1,2047)]; 



subplot(2,1,2),plot(abs(a));