matlab进行HARQ协议仿真

matlab进行HARQ协议仿真，采用汉明码纠错编码，BPSK调制，分析BER和SNR的关系

%modify date: 08.12.10
% 模拟的假设，信道情况稳定，不会有丢包现象，不需要处理超时的情况。
% 不需要处理重复帧的情况


clear all
% N = 100000; 

%计算汉明码，m为数据位
m=3;
%总长计算，2^m-1 > k + m
n = 2^m-1;
% 发包相关计算
k = 4;%（nk）

length_block = k*20;
num_block = 1500; %分组数
N = length_block*num_block;

num_minerr = 10;  % 测试过程中，出错数达到最小错误停止
num_maxARQ = 1; %允许最大传输次数

snr_array = 0:11;%14;                % 0 1 2 3 4 5 6 ... 14  信噪比序列，递进测试
ber_array = zeros（size（snr_array））;  % 0 0 0 0 0 0 0 ... 0

seed = [487944 44984];
% 这个是想干嘛？感觉没用
rand（‘state‘ seed（1））;  %===>randsrc
randn（‘state‘ seed（2））; %===>awgn

%产生传输随机信息
alphabet = [0 1];
% 每行为一个block 共1500个block，每个80bit
msg_orig = randsrc（num_block length_block alphabet）;

%（74）汉明编码 
msg_encode = zeros（num_blocklength_block*n/k）;
for msg_row = 1:num_block
    % 每行分别汉明码 nk为参数 k
    msg_encode（msg_row:） = encode（msg_orig（msg_row:）nk‘hamming‘）;
end;

%BPSK调制
msg_tx = msg_encode*2-1;


%加噪ARQ过程
% 对不同的信噪比（递增的）进行测试
for i = 1:length（snr_array）
    n_trans = 1;
    % 一次完整的ARQ-SAW（停等）过程 送进去的数据是15000 * 80大小的数据 
    n_biterr = transmit（msg_origmsg_txnum_blocklength_blocknum_maxARQsnr_array（i）nk）;
    % 测试达到10次错误所经过的总传输次数（这里是指针对同一个信噪比）
    while（n_biterr < num_minerr）
        n_trans = n_trans + 1;
        n_biterr = n_biterr + transmit（msg_origmsg_txnum_blocklength_blocknum_maxARQsnr_array（i）nk）;
    end;
    % 统计这个信噪比的误码率  （总的错误比特 / 总比特数）
    ber_array（i） = n_biterr/（N * n_trans）;
end;

%作图
%0  ‘r*‘ ‘r-‘50
%1  ‘b+‘ ‘b:‘20
%2  ‘g.‘ ‘g-.‘10
%3  ‘kx‘ ‘k--‘40
hold on;%
semilogy（snr_arrayber_array‘kx‘snr_arrayber_array‘k--‘）; 
axis tight; 
axis（[0 15 0 1]）;                   %             
xlabel（‘SNR（dB）‘）;                                           
ylabel（‘BER‘）; 
legend（‘Simulation‘‘Analytical‘）; %
hold off; %

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        2060  2012-07-17 19:40  harq.m
     文件        2994  2012-07-17 19:40  transmit.m