AWGN信道传输中QPSK调制系统的检测性能.docVIP

AWGN信道传输中QPSK调制系统的检测性能 Matlab程序源码如下： clear all; %对随机序列产生的方法初始化 s = RandStream(mt19937ar,Seed, 5489); RandStream.setDefaultStream(s); s = RandStream(mcg16807, Seed,0); RandStream.setDefaultStream(s); N = 10^5; %码元数 Eb_N0 = -4:20; %Eb/N0比特信噪比 s_data = zeros(1,N); ErrorCount=zeros(1,length(Eb_N0)); %预先分配内存，用于记录错误的码元数 for j = 1:length(Eb_N0) p = (2*(rand(1,N)0.5)-1) + 1i*(2*(rand(1,N)0.5)-1); s = (1/sqrt(2))*p; %归一化 n = 10^(-Eb_N0(j)/40)*1/sqrt(2)*(randn(1,N) + 1i*randn(1,N)); %加性高斯白噪声 r = s + n; %解调 r_re = real(r); %实部 r_im = imag(r); %虚部 s_data(r_re 0 r_im 0) = -1 + -1*1i; s_data(r_re = 0 r_im 0) = 1 + 1*1i; s_data(r_re 0 r_im = 0) = -1 + 1*1i; s_data(r_re = 0 r_im 0) = 1 - 1*1i; ErrorCount(j) = size(find(p- s_data),2); %错误的码元数 end sim_QPSK = ErrorCount/N; %仿真QPSK误码率 theory_QPSK = erfc(sqrt(0.5*(10.^(Eb_N0/20)))) - (1/4)*(erfc(sqrt(0.5*(10.^(Eb_N0/20))))).^2; %理论QPSK误码率 close all figure semilogy(Eb_N0,theory_QPSK,b.-); hold on semilogy(Eb_N0,sim_QPSK,mx-); axis([-4 20 10^-3 1]) grid on legend(理论QPSK误码率, 仿真QPSK误码率); xlabel(Eb/No(dB)) ylabel(BER) title(QPSK误码率的仿真与理论曲线比较) 运行结果： 分析比较： 有图可见QPSK仿真误码率曲线和理论误码率曲线几乎重合在一起，说明实验方法是正确可行的。 Box-Muller方法产生噪声和Matlab中AWGN方法产生噪声的差异： Box-muller方法的基本原理是先生成瑞利分布，然后乘以一个与之独立的(0，2pi)上均匀分布的余弦。无论是wgn还是awgn函数，实质都是由randn函数产生的噪声，即wgn函数中调用了randn函数，而awgn函数中调用了wgn函数。