一、
M=16;
k=log2(M);
x=randint(30000,1);%产生二进制随机数
y=modulate(modem.qammod('M',16,'InputType','Bit'),x);%调制
EbNo=-5:1:10;%信噪比
s_b2d=bi2de(reshape(x,k,length(x)/k).','left-msb');%二进制变为十进制
for n=1:length(EbNo)
snr(n)=EbNo(n)+10*log10(k);%Ratio of symbol energy to noise power spectral density
ynoisy=awgn(y,snr(n),'measured');%加入高斯白噪声
z=demodulate(modem.qamdemod('M',16,'OutputType','Bit'),ynoisy);%解调
r_b2d=bi2de(reshape(z,k,length(z)/k).','left-msb');%二进制变为十进制
[sym(n),sym_rate(n)]=symerr(s_b2d,r_b2d);%计算仿真误码率,不是误比特率。
theo_sym_rate(n)=(3/2*erfc(sqrt(k*0.1*(10.^(EbNo(n)/10)))))*(1-3/8*erfc(sqrt(k*0.1*(10.^(EbNo(n)/10)))));
%《移动通信》P55的式子,已化简。
end
disp(sym);
disp(sym_rat