single_packet.m

来自「OFDM系统的仿真代码MATLAB环境下」· M 代码 · 共 47 行

function r1 = single_packet

global sim_consts;

for k = 1:length(sim_consts.SNR)
    
    % Generate tx signal, returns also information bits and raw bits
    [txsignal, tx_inf_bits, tx_raw_bits] = transmitter;

    % Four-path Rayleigh Channel
    delay = [0,5,8,11];
    [path_1, path_2, path_3, path_4] = multipath_channel;

    % Signals that have passed through the channel
    number_txsignal = length(txsignal);
    for i = 1:delay(2)-delay(1)
        rxsignal(i) = path_1(i).*txsignal(i-delay(1));
    end
    for i = delay(2)-delay(1)+1:delay(3)-delay(1)
        rxsignal(i) = path_1(i).*txsignal(i-delay(1))+path_2(i).*txsignal(i-delay(2));
    end
    for i = delay(3)-delay(1)+1:delay(4)-delay(1)
        rxsignal(i) = path_1(i).*txsignal(i-delay(1))+path_2(i).*txsignal(i-delay(2))+...
            path_3(i).*txsignal(i-delay(3));
    end
    for i = delay(4)-delay(1)+1:number_txsignal
        rxsignal(i) = path_1(i).*txsignal(i-delay(1))+path_2(i).*txsignal(i-delay(2))+...
            path_3(i).*txsignal(i-delay(3))+path_4(i).*txsignal(i-delay(4));
    end

    % Add AWGN to signal
    len = length(txsignal);
    rxsignal = awgn(txsignal,sim_consts.SNR(k),'measured');
 
    % Receiver, return data bits and undecoded bits
    [rx_inf_bits, rx_raw_bits] = receiver(rxsignal, tx_raw_bits);

    % Calculate bit errors
    traceback_len = 2;
    [n1(k),r1(k)] = biterr(rx_raw_bits,tx_raw_bits);
    [n2(k),r2(k)] = biterr(rx_inf_bits(traceback_len+1:end),tx_inf_bits(1:end-traceback_len));

    % disp(['The bit error rates are:   ',num2str([r1 r2])]);
end