function output=cnv_encd(G,k0,input)
% cnv_encd(G,k0,input),k0是每一時鐘周期輸入編碼器的bit數，
% G是決定輸入序列的生成矩陣，它有n0行L*k0列n0是輸出bit數，
% 參數n0和L由生成矩陣G導出，L是約束長度。L之所以叫約束長度
% 是因為編碼器在每一時刻里輸出序列不但與當前輸入序列有關，
% 而且還與編碼器的狀態有關，這個狀態是由編碼器的前(L-1)k0。
% 個輸入決定的,通常卷積碼表示為(n0,k0,m)，m=(L-1)*k0是編碼
% 器中的編碼存貯個數，也就是分為L-1段，每段k0個
% 有些人將m=L*k0定義為約束長度，有的人定義為m=(L-1)*k0

% 查看是否需要補0，輸入input必須是k0的整數部 
if rem(length(input),k0)>0
    input=[input,zeros(size(1:k0-rem(length(input),k0)))];
end
n=length(input)/k0;

% 檢查生成矩陣G的維數是否和k0一致 
if rem(size(G,2),k0)>0
    error('Error,G is not of the right size.')
end

% 得到約束長度L和輸出比特數n0
 
L=size(G,2)/k0;
n0=size(G,1);
% 在信息前后加0，使存貯器歸0，加0個數為(L-1)*k0個
u=[zeros(size(1:(L-1)*k0)),input,zeros(size(1:(L-1)*k0))];

% 得到uu矩陣,它的各列是編碼器各個存貯器在各時鐘周期的內容 
u1=u(L*k0:-1:1);
%將加0后的輸入序列按每組L*k0個分組，分組是按一比特增加
%從1到L*k0比特為第一組，從2到L*k0+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);
% 得到輸出，輸出由生成矩陣G*uu得到
output=reshape(rem(G*uu,2),1,n0*(L+n-1));