softviterbiqaminrayleigh.m

来自「Hard Viterbi QPSK in AWGN, Rayleight so」· M 代码 · 共 15 行

SNR=[-5:1:10];
snr = 10.^(SNR/10); 
result_16qam = fopen('result_16qam_ViterbiSoft_Rayleigh.txt','r');
BER_16qam = fscanf(result_16qam,'%g');
fclose(result_16qam);

semilogy(SNR,BER_16qam,'c-s')
hold on
grid on
for k = 1: length(SNR)
theo_qam(k)=3/8*(1-sqrt(6*snr(k)/(15+6*snr(k))));
end
semilogy(SNR,theo_qam,'k->')
hold on
grid on