卷积编码和viterbi译码的matlab程序

function vit R=[1 1 0 1 1 0 0 1 1 0 0 1 0 1 0 1 0 0 0 0 0 0 1 0 1 0 0 0 1 1 0 0 1 0 1 1 0 1 1 0 0 1 0 0 1 1 1 1 0 1 0 0 1 1 1 0 0 1 1 1 0 0 0 0 1 0 1 0 0 1 0 1 0 1 0 0 0 0 0 0 0 1 0 0 0 0 0 1 0 1 1 1 1 1 1 1 0 1 0 0 0 1 0 1 0 1 1 0 1 0 1 1 ]; N=length(R); statedst=[0,15*(ones(1,63))]; %存储汉明码距初始化 statepath=zeros(64,N/2); %存储估计路径初始化 code=zeros(1,2); minstate=zeros(1,N/2); STBM=[0 0; %buf1 %存储32个蝶型的最初始分支度量值 1 0; 1 1; 0 1; 1 1; 0 1; %buf6 0 0; 1 0; 0 0; 1 0; 1 1; 0 1; %buf12 1 1; 0 1; 0 0; 1 0; 0 1; 1 1; %buf18 1 0; 0 0; 1 0; 0 0; 0 1; 1 1; %buf24 0 1; 1 1; 1 0; 0 0; 1 0; 0 0; %buf30 0 1; 1 1]; for k=1:N/2 code=R(2*k-1:2*k); %读取接受到的两位码子 dst_ram=statedst; %读取前一时刻码距 path_ram=statepath; %读取前一时刻路径估计 for i=1:32 j=32+i; %i，j，p，q分别是每一个蝶型的状态 p=2*i-1; q=2*i; BM10=STBM(i,; %BM10，BM20，BM11，BM21分别是每一个蝶型的分支度量值 BM20=~BM10; BM11=~BM10; BM21=BM10; dst1=dst_ram(i); %dst1，dst2是每一个蝶型两个输入的码距 dst2=dst_ram(j); [dst3,path3,dst4,path4]=BUF(code,dst1,BM10,BM11,dst2,BM20,BM21); %调用蝶型算法 statedst(p)=dst3; % dst3,path3,dst4,path4是每一个蝶型两个输出的码距和路径估计 statepath(p,k)=path3; statedst(q)=dst4; statepath(q,k)=path4; end end statepath %traceback dec=zeros(1,56); t=6; currsta_bin=zeros(1,6); inista_bin=zeros(1,6); nextsta_bin=zeros(1,6); inista_dec=0; currsta_dec=0; nextsta_dec=0; for i=42:1:56 inista_dec=minstate(i); inista_bin=deci2bin(inista_dec-1,t); inivalue=statepath(inista_dec,i); for j=42:-1:1 if j==42 currsta_dec=inista_dec; currsta_bin=inista_bin; currvalue=inivalue; else currsta_dec=nextsta_dec; currsta_bin=deci2bin(currsta_dec-1,t); currvalue=statepath(currsta_dec,i+j-42); end nextsta_bin=[currsta_bin(2:6),currvalue]; nextsta_dec=bin2deci(nextsta_bin)+1; end dec(i-41)=currsta_bin(1); end function c=comp(a,b) M=length(a); c=0; for i=1:M if a(i)~=b(i) c=c+1; end end ％％％％％％％％％％％％％％％％％％％％％％％％％％％％％％％％％％ 函数BUF是蝶型算法。 function [dst3,path3,dst4,path4]=BUF(code,dst1,BM10,BM11,dst2,BM20,BM21) d1=dst1+comp(code,BM10); d2=dst2+comp(code,BM20); tmp1=0; tmp2=1; if d1<d2 dst3=d1; path3=tmp1; else dst3=d2; path3=tmp2; end d1=dst1+comp(code,BM11); d2=dst2+comp(code,BM21); if d1<d2 dst4=d1; path4=tmp1; else dst4=d2; path4=tmp2; end ％％％％％％％％％％％％％％％％％％％％％％％％％％％％％％％％％％ function y=deci2bin(x,t) % 十进制数x转化为二进制数，二进制数至少表示为t位 y=zeros(1,t); i=1; while x>=0&i<=t y(i)=rem(x,2); x=(x-y(i))/2; i=i+1; end y=y(t:-1:1); ％％％％％％％％％％％％％％％％％％％％％％％％％％％％％％％％％％％％％ function y=bin2deci(x) % 将二进制数转化为十进制数 t=length(x); y=(t-1:-1:0); y=2.^y; y=x*y'; 问题基本解决，原因有两个： （1），在于蝶型算法的地址控制，输入的和输出节点的关系为 p=1:32 q=32+p; i=2*p-1; j=2*p; i，j分别为输入，p，q为输出。 （2），分支度量矩阵要重新生成一个： STBM=[0 0 1 1 1 1 0 0 0 1 1 0 1 0 0 1 0 0 1 1 1 1 0 0 0 1 1 0 1 0 0 1 1 1 0 0 0 0 1 1 1 0 0 1 0 1 1 0 1 1 0 0 0 0 1 1 1 0 0 1 0 1 1 0 1 1 0 0 0 0 1 1 1 0 0 1 0 1 1 0 1 1 0 0 0 0 1 1 1 0 0 1 0 1 1 0 0 0 1 1 1 1 0 0 0 1 1 0 1 0 0 1 0 0 1 1 1 1 0 0 0 1 1 0 1 0 0 1 1 0 0 1 0 1 1 0 1 1 0 0 0 0 1 1 1 0 0 1 0 1 1 0 1 1 0 0 0 0 1 1 0 1 1 0 1 0 0 1 0 0 1 1 1 1 0 0 0 1 1 0 1 0 0 1 0 0 1 1 1 1 0 0 0 1 1 0 1 0 0 1 0 0 1 1 1 1 0 0 0 1 1 0 1 0 0 1 0 0 1 1 1 1 0 0 1 0 0 1 0 1 1 0 1 1 0 0 0 0 1 1 1 0 0 1 0 1 1 0 1 1 0 0 0 0 1 1] 其他算法实现仍然按照以上程序。