文档介绍：卷积码的编译码MATLAB程序————————————————————————————————作者:————————————————————————————————日期: %survivorstate是一个矩阵,它显T了通过网格的最优路径,这个矩阵通过一个单独的函数metric(x,y)给出。%其中G是一个矩阵,%这里,我们做了一个简单的(2,1,7)卷积码编码器。k=1;G=[1011011;1111001];%G1=133,G2=171%以下3种输入序列,可任选一种%%input=[0000000];%全0输入%input=[1111111];%全1输入input=[round(rand(1,7)*1)];%随机系列输入,也可用randint(1,7,[01])figure;plot(input,'*r')%figure1:画图:目标input,红色(red,r),形状为*s=input;g1=G(1,:);g2=G(2,:);c1=conv(s,g1);%作卷积%disp(c1);c2=conv(s,g2);%disp(c2);n=length(c1);%7位输入时n=13c=zeros(1,2*n);%生成全0矩阵,1*26%disp(c);fori=1:nc(2*i-1)=c1(i);c(2*i)=c2(i);%两个模2加法器分别输出卷积结果序列后,由旋转开关读取的结果(此时仅为卷积结果,非2进制0/1)endfori=1:2*nif(mod(c(i),2)==0)%mod(c(i),2)==0意思:c(i)除以2,余数为0c(i)=0;elsec(i)=1;endendoutput=c;channel_output=output;%输出矩阵%disp(channel_output);figure;plot(output,'*b')%画图:目标:卷积码编码输出,蓝色(blue,b)*%————————————————以上为编码部分,以下为维特比译码————————————————n=size(G,1);%取矩阵G的行数,故n=2。即得到输出端口,即2个模2加法器%检验G的维数ifrem(size(G,2),k)~=0%当矩阵G的列数不为k的整数倍时,rem为求余函数error('SizeofGandkdonotagree')%报错endifrem(size(channel_output,2),n)~=0%当输出矩阵的列数不是输出端口n的整数倍时。(注:size(channel_output,2)=26,2个模2加法器合成的输出)error('channleoutputnotoftherightsize')endL=size(G,2)/k;%得出移位数,即寄存器个数,此例程为7%由于L-1个寄存器的状态即可表示出输出状态,所以总的状态数number_of_states可由前L-1个寄存器的状态组合来确定number_of_states=2^((L-1)*k);%此例程中2^6,移位寄存器组的状态数为64个%产生状态转移矩阵,输出矩阵和输入矩阵forj=0:number_of_states-1%表示当前寄存器组的状态。因状态从0开始,所以循环为从0到number_of_states-