main.m

交织码的突发信道仿真研究

M=input('输入初始信号长度');
K=4;
N=7;
a=round(rand(1,M));%产生随机序列
b=hanming(a,N,K);%经过汉明编码器


%%%%%%%%%%%%%%%%%%%%%%%%%%%%研究交织矩阵规模对误码率的影响%%%%%%%%%%%%%%%%%%%%%%%%%%%
for i=1:6
    m=[5,5,10,20,25,25];n=[2,5,10,10,10,20];
    for i1=0:(length(b)/(m(i)*n(i))-1)%经过交织器，m为交织矩阵行数，n为交织矩阵列数
        c1([1:m(i)*n(i)]+i1*m(i)*n(i))=jiaozhi(b([1:m(i)*n(i)]+i1*m(i)*n(i)),m(i),n(i));
    end;

    d1=tufazaosheng(100,10,length(b));
    e1=xor(c1,d1);%经过突发性噪声信道

    for i11=0:(length(b)/(m(i)*n(i))-1)
        f1([1:m(i)*n(i)]+i11*m(i)*n(i))=fanjiaozhi(e1([1:m(i)*n(i)]+i11*m(i)*n(i)),m(i),n(i));%经过反交织器
    end;

    g1=fanhanming(f1,N,K);%经过汉明译码器
    h1(i)=length(find(a~=g1))/M;%纠错性能的研究
end;
figure(1)
plot([5*2,5*5,10*10,20*10,25*10,25*20],h1);
xlabel('交织器容量 5*2    5*5    10*10    20*10    25*10    25*20');
ylabel('误码率');


%%%%%%%%%%%%%%%%%%%%%%%研究突发噪声的突发间隔对误码率的影响%%%%%%%%%%%%%%%%%%%%%%%%%%%
    for i2=0:(length(b)/100-1)
        c2([1:100]+i2*100)=jiaozhi(b([1:100]+i2*100),10,10);%经过交织器
    end;

for ii=1:20
    d2=tufazaosheng(ii*10,5,length(b));%10--200的突发间隔
    e2=xor(c2,d2);%经过突发性噪声信道

    for i22=0:(length(b)/100-1)
        f2([1:100]+i22*100)=fanjiaozhi(e2([1:100]+i22*100),10,10);%经过反交织器
    end;

    g2=fanhanming(f2,N,K);
    h2(ii)=length(find(a~=g2))/M;%纠错性能的研究
end;
figure(2)
plot([1:20]*10,h2);
xlabel('突发间隔');
ylabel('误码率');

%%%%%%%%%%%%%%%%%%%%%%%研究突发噪声的突发长度对误码率的影响%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    for i3=0:(length(b)/100-1)
        c3([1:100]+i3*100)=jiaozhi(b([1:100]+i3*100),10,10);%经过交织器
    end;

for iii=1:30
    d3=tufazaosheng(100,iii,length(b));
    e3=xor(c3,d3);%经过突发性噪声信道

    for i33=0:(length(b)/100-1)
        f3([1:100]+i33*100)=fanjiaozhi(e3([1:100]+i33*100),10,10);%经过反交织器
    end;

    g3=fanhanming(f3,N,K);
    h3(iii)=length(find(a~=g3))/M;%纠错性能的研究
end;
figure(3)
plot([1:30],h3);
xlabel('突发长度');
ylabel('误码率');

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