function output=cnv_encd(g,k0,input)
%output=cnv_encd(g,k0,input) 卷積碼編碼函數
%g  生成矩陣
%k0 輸入碼長
%input  輸入信源序列
%output 輸出卷積編碼序列

if rem(length(input),k0)>0
   input=[input,zeros(size(1:k0-rem(length(input),k0)))];
end
n=length(input)/k0;

if rem(size(g,2),k0)>0
   error('Error,g is not of the right size.')
end

li=size(g,2)/k0;
n0=size(g,1);

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

u1=u(li*k0:-1:1);
for i=1:n+li-2
   u1=[u1,u((i+li)*k0:-1:i*k0+1)];
end
uu=reshape(u1,li*k0,n+li-1);

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