clear all;
close all;

b_plot_volt = 0;
b_plot_psd  = 0;

dname = '../capt';
fs = 4e6;
fc = 1.823e9;
Nfft = 1024;

f = fftfreq(Nfft,fs,1);

D = dir(fullfile(dname,'*.dat'));

for n = 1:length(D)

    fname = D(n).name;
    fid = fopen(fullfile(dname,fname));

    x = getchunk(fid,fs,0.4,0.9);  
    x = x/(2^15-1);
    
    if b_plot_volt == 1
        plot(real(x(1:100))); hold all; plot(imag(x(1:100))); legend('Real','Imag')
        pause; close;
    end

    X = avepsd(x,Nfft);
    X = fftshift(X);
    [v k] = max(X);

    if b_plot_psd == 1
        ppsd( avepsd(x,1024), fs)
        pause; close;
    end

    par(n) = mean(real(x).^2)/50;
    pai(n) = mean(imag(x).^2)/50;
    pave(n) = mean(abs(x).^2)/50;

    disp('---------------------------------------');
    disp(fname);
    disp(sprintf('Max real(x)  %g',max(abs(real(x)))));
    disp(sprintf('Max imag(x)  %g',max(abs(imag(x)))));
    disp(sprintf('RMS real(x)  %g',sqrt(mean(real(x).^2))));
    disp(sprintf('RMS imag(x)  %g',sqrt(mean(imag(x).^2))));
    disp(sprintf('Pave real(x) %g',par(n)));
    disp(sprintf('Pave imag(x) %g',pai(n)));
    disp(sprintf('Pave         %g',pave(n)));
    disp(sprintf('Peak Freq    %g',f(k)+fc));
    disp('---------------------------------------');

end