%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?%?????????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