meritfactor_comp.asv

calculation of the autocorrelation peak sidelobe levels and corsscorrelation. CE sequence, frank aut

% clear, close all;

p=1;
clear,close all;

sqL=10;
L=128;
M=3;
s=deng_seq(L);
% nfft=2*L+1;
% x=zeros(M,201,nfft);
for m=1:M
    [x,nfft]=autocorrelationwin(L,s(m,:),0);
    
    xnew            = x(1,:);
    xnew((end+1)/2) = 0;
    SL1             = sum(abs(xnew).^2);
    
   PSL(m)=10*log10(max(abs(xnew)) );

    F(m)=10*log10(1/SL1);
    R(m)=10*log10(1/abs(PSL(m)));

    % end of merit factor calculations
    
    % start of Doppler merit factor calculations
    
    fdT_t=0;
    j=1;
    
    Ffd=F(m);
    [r,c] = size(x);
    while Ffd>3
        fdT_t           =   fdT_t+0.01;
        j               =   j+1;
        % first find tha value of mainlobe
        xnew            =   x(j,:);
        ml              =   xnew((end+1)/2);
        xnew((end+1)/2) =   0;
        sl              =   sum(abs(xnew).^2);
        Ffd=10*log10(ml^2/sl);
    end
   
    fd_t(m)=fdT_t;
    Am(m,:,:)=x;
end


%Mean_F=10*log10( sum(10.^(F/10))/M );
%Mean_R=10*log10( sum(10.^(R/10))/M );
Mean_F          =   sum(F)/M;
Mean_fd_t=sum(fd_t)/M;


time=(-(nfft/2)):1:(nfft/2-1);
fdT=0.:.01:2.;
% box off;
figure;
Am=20*log10(Am);
% subplot(2,2,1);
for n=1:201
A1(n,:)=Am(1,n,:);
end

image(time,fdT,A1,'CDataMapping','scaled');
% pcolor ()
% shading interp;
xlabel('sample shift');
ylabel('Doppler shift fd*T');
zlabel('Autocorrelation function');
colorbar;
caxis([-40 0]);

% subplot(2,2,2);
% for n=1:201
% A2(n,:)=Am(2,n,:);
% end
% image(time,fdT,A2,'CDataMapping','scaled');
% 
% xlabel('sample shift');
% ylabel('Doppler shift fd*T');
% zlabel('Autocorrelation function');
% colorbar('YLim',[-40 0]);

% subplot(2,2,3);
% for n=1:201
% A3(n,:)=Am(3,n,:);
% end
% image(time,fdT,A3,'CDataMapping','scaled');
% 
% xlabel('sample shift');
% ylabel('Doppler shift fd*T');
% zlabel('Autocorrelation function');
% colorbar('YLim',[-40 0]);
% 
% subplot(2,2,4);
% for n=1:201
% A4(n,:)=Am(4,n,:);
% end
% image(time,fdT,A4,'CDataMapping','scaled');
% 
% xlabel('sample shift');
% ylabel('Doppler shift fd*T');
% zlabel('Autocorrelation function');
% colorbar('YLim',[-40 0]);