awgnpeg.m

AWGN信道下高斯噪声干扰下GALILEO编码译码系统性能仿真

%高斯信道中误码率
function [p1,p2]=AwgnPeG(snr_in_dB)
E=1;
snr=10^(snr_in_dB/10);
sgma=sqrt((E/snr)/2);%噪声方差
N=200;%桢个数
input=round(rand(1,N*25));%信息序列

u1=2*input-1;
u1=sign(u1+sgma*randn(size(u1))+eps);
u1=(u1+1)/2;
p1=(length(find(input~=u1)))/length(input);%没有经过信道编译码的误码率

%信道编码
for n=0:(N-1)
	%外编码
	out_en=T_SCH_EX_EN(input([1:25]+n*25));
	%内编码
    in_en=SCH_IN_EN(out_en);
    %合成桢
    u([1:78]+n*78)=in_en;
end;

%经信道加高斯燥声
u=2*u-1;
u=sign(u+sgma*randn(size(u))+eps);
u=(u+1)/2;

%信道译码
for n=0:(N-1)
	%分离桢
    de=u([1:78]+n*78);
    %内编码解码
	in_de=SCH_IN_DE(de);
	%外编码解码
   [cc,output]=cyc_sch_decoder(in_de);
   %桢合成
   receive([1:25]+n*25)=output;
end;

p2=(length(find(input~=receive)))/length(input);%经过信道编译码的误码率