ricebtest.m

来自「MATLBA 信道仿真,完成信道卷积码循环码的编译码」· M 代码 · 共 25 行

echo off
clear all
close all
snr_in_DB=-4:1:6;
SNR_in_DB=-4:0.1:6;

for i=1:length(snr_in_DB)
   [py(i),pb3(i),pb5(i)]=RicePeB(snr_in_DB(i));%py为突发长度为3没有经过信道编码的误码率，pb3为突发长度为3经过信道编码的误码率，pb5为突发长度为5经过信道编码的误码率，
end;
for i=1:length(SNR_in_DB),
   SNR=10^(SNR_in_DB(i)/10);
   p(i)=Qfunct(sqrt(2*SNR));%计算理想误码率p
end; 
%画出误比特率图
semilogy(snr_in_DB,py,'b*');
hold on
semilogy(snr_in_DB,pb3,'r+');
semilogy(snr_in_DB,pb5,'g+');
semilogy(SNR_in_DB,p);
py
pb3
pb5