survivorstate是一个矩阵,它显T了通过网格的最优路径,这个矩阵通过一个单独的函数metricxy给出。其中G是一个矩阵,它的任一行决定了从移位寄存器到模2加法器的连接方式为生成矩阵这里,我们做了一个简单的217卷积码编码器。k1G10110111111001G1133G2171以下3种输入序列,可任选一种i
put0000000全0输入i
put1111111全1输入i
putrou
dra
d171随机系列输入,也可用ra
di
t1701figureploti
putrfigure1:画图:目标i
put,红色(red,r),形状为si
putg1G1g2G2c1co
vsg1作卷积dispc1c2co
vsg2dispc2
le
gthc17位输入时
13czeros12
生成全0矩阵,126dispcfori1
c2i1c1ic2ic2i两个模2加法器分别输出卷积结果序列后,由旋转开关读取的结果(此时仅为卷积结果,非2进制01)e
dfori12
ifmodci20modci20意思:ci除以2,余数为0ci0elseci1
e
de
doutputccha
el_outputoutput输出矩阵dispcha
el_outputfigureplotoutputb画图:目标:卷积码编码输出,蓝色(blue,b)以上为编码部分,以下为维特比译码
sizeG1取矩阵G的行数,故
2。即得到输出端口,即2个模2加法器检验G的维数ifremsizeG2k0当矩阵G的列数不为k的整数倍时,rem为求余函数
errorSizeofGa
dkdo
otagree报错e
difremsizecha
el_output2
0当输出矩阵的列数不是输出端口
的整数倍时。(注:sizecha
el_output226,2个模2加法器合成的输出)
ferrorcha
leoutput
otoftherightsize
e
dLsizeG2k得出移位数,即寄存器个数,此例程为7由于L1个寄存器的状态即可表示出输出状态,所以总的状态数
umber_of_states可由前L1个寄存器的状态组合来确定
umber_of_states2L1k此例程中26,移位寄存器组的状态数为64个产生状态转移矩阵,输出矩阵和输入矩阵forj0
umber_of_states1表示当前寄存器组的状态。因状态从0开始,所以循环为从0到
umber_of_states1
fort02k1k位输入端的信号组成的状态,总的状态数为2k,所以循环为从0到2k1
ext_statememory_co
te
ts
xt_statjtLk
xt_stat完成从当前的状态和输入的矢量得出下寄存器组的一个状态
i
putj1
ext_state1ti
put数组值是用于记录当前状态到下一个状态所要的输入信号矢量
i
put数组的维数:一维坐标xj1指当前状态的值,二维坐标y
ext_state1指下一个状态的值
由于Matlab中数组的下标是从1开始的,而状态值是从0开始r