QPSK数字调制系统在AWGN信道下的BER仿真

四相移相键控（QPSK）

QPSK将每两个比特流映射到四个矩形星座点上，本文的映射规则如下

双比特码	符号
00	1+1j
10	-1+1j
11	-1-1j
01	1-1j

假设信源0和1比特等概出现，则QPSK系统的误比特率在AWGN信道下的理论误比特率同BPSK，即

$$P_b=\frac{1}{2}erfc(\sqrt{\frac{E_b}{N_0}})$$

$SNR$与$\frac{E_b}{N_0}$的换算关系

$$\frac{S}{N}=\frac{E_b}{N_0}\ast \frac{R_c{R}_m}{(1+\alpha )d}$$

其中，$\frac{S}{N}$是信号噪声功率比，$\frac{E_b}{N_0}$是比特信噪比，$R_c$是信道编码速率，$R_m$是调制率，$\alpha$是滚降因子，$d$是扩频倍数。

两边取$10lo{g}_{10}$转化为$dB$形式，即

$$SNR=(\frac{E_b}{N_0}{)}_{dB}+10lo{g}_{10}(R_c)+10lo{g}_{10}(R_m)-10lo{g}_{10}(1+\alpha )-10lo{g}_{10}(d)$$

具体推导过程见文献1

Matlab仿真与理论误比特率对比

close all;clc;clear
qpsk_code=[1+1j  1-1j  -1+1j   -1-1j]./sqrt(2);%归一化使每个符号的功率为1
Num=100;
SNR_dB = -10:2:16;
SNR = 10.^(SNR_dB/10);
NoiPower = 1./SNR;
sigma = sqrt(NoiPower/2);
for snr_N = 1:length(SNR)
    ber = zeros(1,Num);
    for tmp = 1:Num
        qpsk_bit = randi([0,1],5e4,1);
        qpsk_bit2 = reshape(qpsk_bit,2,[]);
        data_signal = qpsk_code([2 1]*qpsk_bit2 + 1);%比特映射到符号
        noise = sigma(snr_N).*((randn(1,length(data_signal)))+ 1i*randn(1,length(data_signal)));;
        R_sig=data_signal+noise;
        bit_dec=zeros(length(R_sig),1);
        for loop_d = 1:length(R_sig)
            [~,position] = min(abs(R_sig(loop_d)-qpsk_code()));%MAP判决
            bit_dec(2*loop_d-1) = floor((position-1)/2);
            bit_dec(2*loop_d) = mod((position-1),2);
        end
        err_bit = sum(bit_dec~=qpsk_bit);
        ber(tmp) = err_bit/length(qpsk_bit);
    end
    ber_Mu(snr_N,:) = mean(ber,2).';
end
semilogy(SNR_dB,ber_Mu,'ro-')
hold on;grid on
BER = 1/2.*erfc(sqrt(SNR/2));
semilogy(SNR_dB,BER,'bo-')
legend('QPSK仿真','理论误码率')
xlabel('E_b/N_0(dB)')
ylabel('BER')

仿真结果

参考文献

张少侃,吕聪敏,甘浩.数字通信系统中Eb/N0与SNR转换方法的研究[J].现代计算机,2019(12):33-36.