cnv_encd.m

来自「这是关于数字处理信号的MATLAB源代码」· M 代码 · 共 36 行

% 卷积码编码函数
function output=cnv_encd(G,k0,input)
%G  生成矩阵
%k0 输入码长
%input  输入信源序列
%output 输出卷积编码序列

%看填充的0是否必须
if rem(length(input),k0)>0
   input=[input,zeros(size(1:k0-rem(length(input),k0)))];
end
n=length(input)/k0;% 把输入比特按k0分组，n为所得的组数。

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

%li,n0大小
li=size(G,2)/k0;
n0=size(G,1);

%以保证编码器是从全0开始，并回到全0状态。
u=[zeros(size(1:(li-1)*k0)),input,zeros(size(1:(li-1)*k0))];

%uu lie i*T,i=1,2...
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
output=reshape(rem(G*uu,2),1,n0*(n+li-1));