bch.m

对数字信号进行分析

function [bch1]=bch(rand1,tb)
runum4=length(rand1);
bch1=zeros(1,runum4);
N4=752;
N5=fix(runum4/N4);
N3=752+261+10;
m4=zeros(1,N3-10);
for i=0:1:N5-1
m4(1,262:N3-10)=rand1(1,i*N4+1:i*N4+752);
g1=[1,0,0,0,0,0,0,1,0,0];
g2=zeros(N3-10,10);
g2(1,:)=g1;
for k=2:1:N3-10
    p1=g1(1,10);
    for j=10:-1:2
        g1(1,j)=g1(1,j-1);
    end
    g1(1,1)=p1;
    g2(k,:)=g1;
end

g3=eye(N3-10);
g4=zeros(N3-10,N3);
g4(:,1:N3-10)=g3;
g4(:,N3-9:N3)=g2;
bch2=m4*g4;

    for j=1:1:N3
      bch2(1,j)=rem(bch2(1,j),2);%进行模2运算
    end 
N6=N3-261;
bch1(1,i*N6+1:i*N6+N6)=bch2(:,262:N3);%bch编码
end
runum5=length(bch1);
t3=tb:tb:tb*runum5;
% figure
% stairs(t3,bch1)
% axis([0 tb*runum5 -0.5 1.5])