nti_tr_uwb.m

来自「这是一个较为完整的tr方式uwb接收机的M文件」· M 代码 · 共 69 行

load tx;
ber=[];
repetition_time = 100;

%--------------------------------------------------------------------------
%--------------------------------------------------------------------------
%--------------------------------------------------------------------------
for i = 1:5
    i
    snr = 5*i - 10;    
  
    %----------------------------------------------------------------------
    %basic settings
    sample_time = 1e-10;
    Ns = 6; %Ts/Tf
    nTp = 10;
    nTd = 50;
    nTf = 120;
    nTs = Ns*nTf;
    data_length = 10000;
   
    %----------------------------------------------------------------------
    %----------------------------------------------------------------------
    for ii = 1:11
        ii
        nTi = 3*ii + 7;
        
        %------------------------------------------------------------------
        berr = [];
        for iii = 1:repetition_time
            rx = awgn(tx,snr);
            rx = reshape(rx,nTs,data_length); 
        
            %--------------------------------------------------------------
            %small trick
            reference_index_select = [];
            for j=1:Ns
                reference_index_select = [reference_index_select,(i-1)*nTf+1:(i-1)*nTf+nTi];
            end

            data_index_select = [];
            for jj=1:Ns
                data_index_select = [data_index_select,(i-1)*nTf+nTp+nTd+1:(i-1)*nTf+nTp+nTd+nTi];
            end

            %--------------------------------------------------------------
            %seperate reference and data pulse
            rx_reference = rx(reference_index_select,:);
            rx_data = rx(data_index_select,:);

            %--------------------------------------------------------------
            %demodulation
            de = rx_reference.*rx_data;
            de = sum(de,1);
            de = de>0;
            source = source>0;
            berr(iii) = (data_length - sum(de==source))/(data_length);
        end
        %------------------------------------------------------------------
        
        ber(i,ii) = sum(berr)/repetition_time;
    end
    %----------------------------------------------------------------------
    %----------------------------------------------------------------------
    
end
%--------------------------------------------------------------------------
%--------------------------------------------------------------------------
%--------------------------------------------------------------------------