放大转发AF模式的matlab仿真

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 =============================
        % **************************