%ssst06
%Ammended code to run simulation for the sssto6 paper
%Build 3 signals- Add noise - Correlate - Display results.
%This code calls: compcor.m
%rjb 10/05/05
%r_short_new
%Edwin Spencer
%9/17/05

%Set up variables that will be used throughout the program
clear
clc
tic
N = 2*1024;
%N = 512;
%N = 256;
t = (0:(N-1))/N;

%Create, normalize, and plot the signals to be manipulated
dop = MakeSignal('Doppler',N);
dop = dop - mean(dop);
dop = dop./max(dop);

blo = MakeSignal('Blocks',N);
blo = blo - mean(blo);
blo = blo./max(blo);

rie = MakeSignal('Riemann',N);
rie = rie - mean(rie);
rie = rie./max(rie);

%Compute and store
snr = -10:15;
for q = 1:length(snr)
    M = 10;
    for m = 1:M
        %ndop =awgn(dop,snr(q));  ndop = ndop - mean(ndop);
        %[t,w,p,r] = compcor(ndop,dop,0);
        nblo =awgn(blo,snr(q));  nblo = nblo - mean(nblo);
        [t,w,p,r] = compcor(nblo,blo,0);
        %nrie =awgn(rie,snr(q));  nrie = nrie - mean(nrie);
        %[t,w,p,r] = compcor(nrie,rie,0);
        tt(m,q)=t;
        %ff(m,q)=f;
        ww(m,q)=w;
        pp(m,q)=p;
        rr(m,q)=r;
        %ss(m,q)=s;
    end
end

%save and plot
TT = mean(tt);  WW = mean(ww); %FF = mean(ff);
PP= mean(pp); RR = mean(rr); %SS=mean(ss);
CORR_OUT = [TT',WW',PP',RR',snr']
%save corr_out CORR_OUT, snr
%figure(2), plot(snr,CORR_OUT(:,1:4)), Title('NCC vs SNR')
figure(2), plot(snr,TT','+',snr,WW','*',snr,PP','o',snr,RR','x')
hold on, plot(snr,TT',snr,WW',snr,PP',snr,RR'), 
title('Peak NCC vs. SNR for Signal 1'),xlabel('SNR  {dB}'), 
axis([-10,15,0,1]),legend('TDC ','WDC ','PHAT ','ROTH ',2),grid
hold off
toc