qpsk_berest.m

仿真具有三条固定信道的AWGN多径信道的QPSK系统的BER性能

% File:qpsk_berest.m
function [peideal,pesystem] = qpsk_berest(xx,yy,ebn0db,eb,tb,nbw)
% ebn0db:接收机输入端的信噪比的数组；
% tb:比特持续时间
% nbw：noise BW； 
% xx:标准(理想)输入  yy:扭曲的输出
[n1 n2] = size(xx); nx = n1*n2;
[n3 n4] = size(yy); ny = n3*n4;
[n5 n6] = size(ebn0db); neb = n5*n6;

% 为对比，设置理想(积分清除)接收机的 noise BW=rs/2
nbwideal = 1/(2*tb*2);
for m=1:neb
    peideal(m) = 0.0; pesystem(m) = 0.0;       % 初始化
    
    % find n0和噪声方差
    string1 = ['Eb/N0 = ',num2str(ebn0db(m))]; % num2str把向量表示成4个数字1个指数的形式
    disp(string1)                              % 跟踪执行
    ebn0(m) = 10^(ebn0db(m)/10);               
    n0 =eb/ebn0(m);
    sigma = sqrt(n0*nbw*2);                    % 方差
    sigma1 = sqrt(n0*nbwideal*2);              % 理想方差
    
    % 输入星座/信号乘比例因子，使其与接受滤波器的输入端的星座/信号有相同的平均功率 a=sqrt(2eb/(2tb))
    b = sqrt(2*eb/tb)/sqrt(sum(abs(xx).^2)/nx);
    for n = 1:nx
        theta = angle(xx(n));
        if(theta<0)
            theta = theta+2*pi;
        end
        
        % 将x,y旋转到第一象限，计算BER
        xxx(n) = b*xx(n)*exp(-i*(theta-(pi/4)));
        yyy(n) = yy(n)*exp(-i*(theta-(pi/4)));
        d1 = real(xxx(n)); d2 = imag(xxx(n));  % 参考
        d3 = real(yyy(n)); d4 = imag(yyy(n));  % 系统
        pe1 = q(d1/sigma1)+q(d2/sigma1);       % 参考
        pe2 = q(d3/sigma1)+q(d4/sigma1);       % 系统
        peideal(m) = peideal(m)+pe1;           % 参考SER
        pesystem(m) = pesystem(m)+pe2;          % 系统SER
    end
end
peideal = (1/2)*peideal./nx;                   % 转换成BER
pesystem = (1/2)*pesystem./nx;                 % 转换成BER
% END