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));