matlab卷积码程序.doc

matlab卷积码程序
matlab卷积码程序
matlab卷积码程序
1、卷积码编码
fｕncｔｉon ［ｏutpuｔ］＝cnv_encd（inpｕｔ）
%output=ｃnv_ｅｎcd(ｇ，ｋ0,ｉｎput) 卷积码编码函数
％g 生成矩阵
％k０ 输入码长
％input 　输入信源序列
％ｏutｐut 输出卷积编码序列
g=[1　1 1；1 ０　1];编码矩阵
k０＝1；
iｎpｕt=［１ 1 ０ 1]；
if　ｒeｍ（length（ｉnｐｕt），k0)〉0
inpｕt=［inpｕｔ,zｅros(siｚe（1:ｋ０—rｅm(length（iｎpuｔ），ｋ0)））］;
eｎd
n＝length(inpｕt)/k0;
ｉf reｍ(sizｅ(ｇ，２），k0）＞0
erroｒ（’Ｅrror,g is noｔ of the right sｉze。'）
eｎｄ
li=size（g,2)/k0；
n0=size（g,1）;
ｕ=［zeros（size（１：(ｌｉ—1)*k0)），iｎｐuｔ,zerｏs（sｉze（１：(ｌｉ－1)*ｋ0)）］；
u1=u（li＊k０:－1:1）;
for i=1：ｎ+ｌｉ—2
　 　u1=［u1,u（(ｉ+ｌi）*k０:-1:i＊k0+1）］；
eｎd
uu＝reｓhape（ｕ1，li*ｋ０，n+lｉ－1);
oｕｔput=reshape（rem(g*ｕu,2）,１,n0*（n＋li—1））；
２、Vｉtｅrbi译码程序
1）
fuｎctiｏｎ y=bin2decｉ（x）
matlab卷积码程序
matlab卷积码程序
matlab卷积码程序
l=leｎgtｈ（ｘ)；
ｙ=(l－1：-１:0）；
y=2.^y；
y=ｘ＊y'；
2）
fuｎctioｎ y=ｄeｃｉ2bin（x，l）
y=zeros（1，l);
i=1;
ｗhile x>=0 & i<=l
y(i）＝rem（x，２)；
x＝（x－y（i））/２；
　 i=ｉ+1；
eｎd
y=y（l:－1:1）；
3)
fｕnction disｔance=ｍｅtrｉｃ(x,ｙ）
iｆ ｘ=＝y
　dｉstanｃｅ=０;
elsｅ
　 dｉstaｎce=1；
eｎd
4）
funｃtioｎ [nｅxt_ｓtaｔe,meｍory_ｃoｎtｅｎｔs]=ｎxｔ＿stat（curreｎｔ_state,inｐut,L，k)
ｂiｎary_statｅ=decｉ2ｂin（ｃurｒｅnt＿sｔａte，k＊（Ｌ－１)）;
biｎarｙ＿ｉｎput=deci２bin(iｎpｕｔ,k)；
next＿state_ｂinarｙ=[bｉnary_iｎput，bｉnａrｙ_stａte(１：(L—２）*k）］；
next_ｓｔａｔe=ｂｉn2deci(nexｔ＿stａte＿binarｙ）;
mｅmｏry＿conｔｅnts＝[ｂiｎａrｙ_input,binary_statｅ］；
５)
fｕｎｃtｉon ［decoder_oｕtｐut,ｓｕｒvｉvｏｒ_state,cｕｍulateｄ_metric]=vit