function????BER?=?AF_Simulation（max_snr）
%?*************************************************************************
%?Time:?21/9/2010
%?Programmer:?（modified?by）?Liu?Huifeng?
%?Function?file:?AF_Simulation.m
%?Description:?This?file?simulates?the?BER?performance?of?AF?cooperative
%??????????????transmission?according?to?SNR.
%?Input:
%???????max_snr:?the?maximum?SNR
%?Output:
%???????BER:?the?bit?error?rate
%?Usage:
%???????BER?=?AF_Simulation（max_snr）
%?=========================================================================

%?*****************************?定義信道增益?*******************************
Dsd=5;??????????????????????????????????%?用dB表示的信道增益
Dsr=4;
Drd=1;

Qsd=sqrt（10^（Dsd/10））;??????????????????%?由dB值換算得到信道增益值
Qsr=sqrt（10^（Dsr/10））;
Qrd=sqrt（10^（Drd/10））;
%?====================================?END?================================
N=10000;????????????????????????????????%?輸入符號數
max_snr=100;
BER=zeros（1max_snr）;???????????????????%?用于存儲后面計算出的BER的值
ber=zeros（1max_snr）;
for?SNR=0:max_snr
????N0=1/（10^（SNR/10））;?????????????????%?N0為方差，噪聲能量歸一化
????sig?=?randsrc（1N[0?1]）;???????????%?產生二進制輸入序列
????sig_mod=QpskMapping（sig）;???????????%?對輸入的二進制序列進行QPSK調制
????sig_det=[];?????????????????????????%?用于存諸對接收信號判決后的QPSK符號?
????
????M=length（sig_mod）;????
????for?i=1:M
????????%?************************?產生高斯白噪聲?**************************
????????Nsd=sqrt（N0/2）*（randn+j*randn）;?????????????%?均值為0，方差為N0
????????Nsr=sqrt（N0/2）*（randn+j*randn）;
????????Nrd=sqrt（N0/2）*（randn+j*randn）;
????????%?=============================?END?===============================
????????%?****************************?信道?********************************
????????Hsd=（Qsd/sqrt（2））*（randn+j*randn）;
????????Hsr=（Qsr/sqrt（2））*（randn+j*randn）;
????????Hrd=（Qrd/sqrt（2））*（randn+j*randn）;
????????%?===============================?END?=============================
????????%?**************************