ricegtest.m

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

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),pb(i)]=RicePeG(snr_in_DB(i));%py为没有经过信道编码的误码率，pb为经过信道编码的误码率
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,pb,'r+');
semilogy(SNR_in_DB,p);
py
pb