function vvec = veltest(n,Wd,N)

vvec = [];
p = doppler(8,32e-6,4e6,8,4096,0,0,0);
r = doppler(8,32e-6,4e6,8,4096,0,Wd,0);

for I=1:n,
   noise = randn(1,length(r));
   rn =N*noise+r;
   pn =N*noise+p;
   vvec = [vvec velcalc(pn,rn,4e6,32768)];
   x=I
end

subplot(1,1,1),plot(vvec)