matlab HARQ仿真

对之前的一个代码进行了改进，可以随意更改重传次数，画出错误率随SNR变化曲线，采用BPSK调制，汉明码

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


clear all
% N = 100000; 

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

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

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

snr_array = 0:10;%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];
% ffy 确保每次产生的随机数都是确定的
rand（‘state‘ seed（1））;  %===>randsrc
randn（‘state‘ seed（2））; %===>awgn

%产生传输随机信息
alphabet = [0 1];
% ffy 每行为一个block 共1500个block，每个80bit
% ffy 产生一个0或1的随机矩阵
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;
end;

%作图
%0  ‘r*‘ ‘r-‘50
%1  ‘b+‘ ‘b:‘20
%2  ‘g.‘ ‘g-.‘10
%3  ‘kx‘ ‘k--‘40
semilogy（snr_arrayber_array）;
grid on;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1933  2014-03-18 09:43  CODE\harq.m

     文件       3120  2014-03-18 09:36  CODE\transmit.m

     目录          0  2014-03-18 22:10  CODE

----------- ---------  ---------- -----  ----

                 5053                    3