返回信息流M=16;
k=log2(M);
n=30000;
x=randint(n,1);
xsym=bi2de(reshape(x,k,length(x)/k).','left-msb');
y=qammod(xsym,M);
ytx=y;
EbNo=-5:0.5:10;
for i=1:length(EbNo)
snr=(i-1)*0.5-5+10*log10(k);
ynoisy=awgn(ytx,snr,'measured');
yrx=ynoisy;
zsym=qamdemod(yrx,M);
z=de2bi(zsym,'left-msb');
z=reshape(z.',prod(size(z)),1);
[number(i),Pe(i)]=biterr(x,z);
end
theoryBer=(1/k)*3/2*erfc(sqrt(k*0.1*(10.^(EbNo/10))));
semilogy(EbNo,Pe,'bs-','LineWidth',2);
hold on;
semilogy(EbNo,theoryBer,'ms-', 'LineWidth',2);
legend('理论', '仿真');
xlabel('信噪比',dB')
ylabel('误码率')
title('16QAM误码率')
这是一条镜像帖。来源:北邮人论坛 / matlab / #9980同步于 2013/5/18
该镜像源已超过 30 天没有更新,可能在源站已被删除。
Matlab机器人发帖
加性高斯信道换成瑞丽信道(有代码),瑞丽信道怎么用——求大
yubin1315
2013/5/18镜像同步1 回复
订阅后,新回复会通过你的通知中心匿名送达。
1 条回复