资源描述：
《卷积码的编译码matlab程序》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。
1、%survivorstate是一个矩阵，它显T了通过网格的最优路径，这个矩阵通过一个单独的函数metric(x,y)给出。%其中G足一个矩阵，它的任一行决定了从移位寄存器到模2加法器的连接方式.为生成矩阵%这里，我们做了一个简单的(2,1,7)卷积码编码器。k=l；G=[1011011;1111001];%G1=133,G2=171%以下3种输入序列，可任选一种％%input=[0000000]；0A%input=[1111111]；1$fi5Ainput=[round(rand(1,7)*1)];%随机系歹(J输入，也川用randint(1,7,[0
2、1])figure；plot(inputz'*r')%figurel：iffij•图：□标input，红色(red,r)，形状s=input;gi=G(i#：)；g2=G(2z：)；cl=conv(s,gl);%作卷积%disp(cl);c2=conv(s,g2);%disp(c2)；n=length(cl);%7位辅5入时n=13c=zeros(1,2*n);%生成全0矩阵，1*26%disp(c);fori=l：nc(2*i-l)=cl(i)；c(2*i)=c2(i);%两个模2加法器分别输出卷积结果序列后,由旋转开关读取的结果(此时仅为卷积结果，
3、非2进制0/1)endfori=l: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,1*b1)%岡图：目标：卷和码编码输出，蓝色(blue,b)*%以上为编码部分，以下为维特比译码resize(G,l);%取知:阵G的行数，故n=2。即得到输出端U,即2个模2加法器%检验G的维数ifrem(size
4、(G,2),k)〜=0%3矩陈0的列数不为k的整数倍时，rem为求余函数error(1SizeofGandkdonotagree1)%报错endifrem(size(channel_output#2),n)~=0%输岀知(阵的列数不足输Hi端Un的幣数倍时。(注:size(channel_output,2)=26,2个模2加法器合成的输出)error(1channleoutputnotoftherightsize1)endL=size(G,2)/k;%得Ji移位数，即寄存器个数，此例程为7%由于L-1个寄存器的状态即川'农示出输出状态，所以总的状态数n
5、umber_of_statesnJ由前L-1个寄存器的状态纟M合來确定number_of_states=2A((L-1)*k)；%此例ft!中2、，移位寄存器组的状态数力64个%产生状态转移矩阵，输出妬阵和输入矩阵forj=0：number_of_states-l%衣示广!府寄存器组的状态。因状态从0升始，所以循环为从0至ljnumber_of_states-1fort=0:2"k-l%1＜位输入端的信兮组成的状态，总的状态数为2"k，所以循环为从0到2人k-1[next_state,memory_contents]=nxt_stat(j,t,L,k)
6、;%nxt_stat究成从当前的状态和输入的矢量得出下寄存器组的一个状态input(j+1,next_state+l)=t；%input数姐值是用于记求当前状态到卜‘一个状态所要的输入信号矢量%input数组的维数：一维坐标x=j+l指当前状态的值，二维坐标y=next_state+l指卜*—个状态的值%由于Matlab中数组的下标足从1开始的，而状态值是从0开始的，所以以上少标值为：状态值+1branch_output=rem(memory_contents*G',2)；%branch_output用pi己录在状态j卜输入1时的输出nextstate
7、(j+1,t+1)=next_state；%nextstate状态转移矢P:陈，id求了当前状态j下输入1吋的下一个状态output(j+1zt+1)=bin2deci(branch_output)；%output记55；了当前状态jl、®5入1时的输出(十进制)endendinput;state_metric=zeros(number_of_states,2)；%state_metric数糾.川十记55：详过在毎个状态时的汉明距离，大小为number_of_states,2%(:,1)为当前状态位置的汉明距离，为确定值；(：，2)为当前状态加输入得到
8、的下一个状态汉明距离，为临时值depth_of_trellis=length(channel_
