function s = dchirp(T,W,p)
%DCHIRP      generate a sampled chirp 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)

%W=20/p;
%T=TW/W;
N=round(p*T*W);
n= [0:1:N-1];
 
s= exp(j*pi*(W/T)*(((n/(p*W))-.5*T).^2) );

%plot(n,s);

%f = fft(s);
%a=sqrt((real(f)).^2 + (imag(f)).^2)
%plot(n,a);