%------------------------------------------------------------------------
%  Plot time-varying transfer function
%------------------------------------------------------------------------
clear all 
tx_sig  = ones(1,1000);
fadeBW1 = 0.025; 
fadeBW2 = 0.025;
KdB1    = -100; 
KdB2    = -100; 
tau1     = 0.0; 
tau2     = 0.25; 
RelPow  = 0.5;
Fs      = 8E0;

%[rx_sig,chansamps1] = RayleighFade(tx_sig, fadeBW1,KdB1,tau, fadeBW2,KdB2, RelPwdB);
[rx_sig,chansamps1] = RayleighFade(tx_sig, fadeBW1,KdB1);
[rx_sig,chansamps2] = RayleighFade(tx_sig, fadeBW1,KdB1);

blksz = 24;
freq  = [-blksz/2 : blksz/2-1]/blksz * Fs; 

for kk = 16:26
    kkk = kk -15;
    A1 = chansamps1(kk*10);
    A2 = chansamps2(kk*10);
    FreqResp = A1*exp(j*2*pi*freq*tau1) + sqrt(RelPow)*A2*exp(j*2*pi*freq*tau2);
    spec = abs(FreqResp).^2;
    spec = spec/blksz;
    plot3(0.1*kk*ones(blksz,1), freq, spec); hold on
    XX(kkk,:) = (kk-15)*2*ones(1,blksz);
    YY(kkk,:) = freq*100;
    ZZ(kkk,:) = 4*spec;
    pause
end

xlabel('Time (s)')
ylabel('Frequency (Hz)')
zlabel('Spectral Response')
grid
view(-37.5,30)
hold off
pause
mesh(XX,YY,ZZ);