cnv_encd1.m

基于matlab的卷积码编码译码程序

function output=cnv_encd(g,k0,input)

%k0=2;

%g=[0 0 1 0 1 0 0 1;0 0 0 0 0 0 0 1;1 0 0 0 0 0 0 1];

%input=[1 0 0 1 1 1 0 0 1 1 0 0 0 0 1 1 1];

%check to see if extra zero padding is necessary

if rem(length(input),k0)>0

    input=[input,zeros(size(1:k0-rem(length(input),k0)))];

end    

  %check the size of matrix g

n=length(input)/k0;

if rem(size(g,2),k0)>0

    error('error,g is not of the right size.')

end   

% determine l and n0

l=size(g,2)/k0;

n0=size(g,1);

%add extra zeros

u=[zeros(size(1:(l-1)*k0)),input,zeros(size(1:(l-1)*k0))];

% gerenate uu,a matrix whose columns are the contents of conv. encoder at

% various clock cycles.

u1=u(l*k0:-1:1);

for i=1:n+l-2

u1=[u1,u((i+l)*k0:-1:i*k0+1)];

end

uu=reshape(u1,l*k0,n+l-1);

% determine the output

output=reshape(rem(g*uu,2),1,n0*(l+n-1));