testtchleave2_8.m

GSM信道编码设计

echo off
clear all
close all
% 产生随机的01序列
c=round(rand(1,456));
%根据公式重排并分成8组
for i=1:8
   	for j=1:57
   	   a=mod(57*(i-1)+64*(j-1)+1,456);
          if a==0
           a=456;
          end
  	    d(i,j)=c(a);
   	end
end
%交织，将前四组加到前四个子组的奇数位，后四组加到后四个子组的偶数位	
   B=zeros(8,114);
	for i=1:4
   	B(i,1:2:114)=d(i,:);
   	B(i+4,2:2:114)=d(i+4,:);
	end
%反交织    
    newD=zeros(8,57);
    for i=1:4
        newD(i,:)=B(i,1:2:114);
        newD(i+4,:)=B(i+4,2:2:114);
    end
%去重排
for i=1:8
   	for j=1:57
   	   a=mod(57*(i-1)+64*(j-1)+1,456);
          if a==0
           a=456;
          end
  	    newC(a)=newD(i,j);
   	end
end
%检测得到序列是否和原序列相同
if(c==newC)
    display "good";
else 
    display "error";
end