test.m

GSM信道编码设计

%
cx=[];
ux=[];
ccx=[];
sys_info='';
tmp=0;

for i=1:184,
   tmp=rand;
   if(tmp<0.5)
      cx(i)=0;
   else
      cx(i)=1;
   end
end

ux=fire_sys_coding(cx);

tmp=fix(rand*184);
for i=1:12,
   ux(tmp+i)=rem(ux(tmp+i)+fix(rand*2),2);
end

[ccx,sys_info]=fire_sys_encoder(ux);

if length(sys_info)>0,
   sys_info
   break
end

if ccx==cx,
   'Correct!'
else
   'Error Exist!'
end

break

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%division circuit

%ccx(x)=uu(x)/g(x)
%%g(x)=x^40+x^26+x^23+x^17+x^3+1
uu=ux;

ccx=[];
gx_out=[];
gx_reg=zeros(1,40);
gx_gg=zeros(1,41);
gx_gg(41)=1;
gx_gg(27)=1;
gx_gg(24)=1;
gx_gg(18)=1;
gx_gg(4)=1;
gx_gg(1)=1;
for i=1:length(uu)
   feedback=gx_gg(41)*gx_reg(40);
   gx_out(i)=gx_reg(40);
   for j=40:-1:2
      gx_reg(j)=rem(gx_reg(j-1)+feedback*gx_gg(j),2);
   end
   gx_reg(1)=rem(uu(i)+feedback*gx_gg(1),2);
end

ccx=gx_out(41:length(uu));
gx_reg

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%