Question

关于matlab程序代码的问题

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));
这是段卷几编码的程序
以上这段代码在matlab7．0中运行出错
Error in ==> cnv_encd at 12
if rem(length(input),k0)>0

Answer 1

if rem(length(input),k0)>0
rem的参数问题，k0是你整个函数的形参，在编译时是未知的，所以报错，你先不要管它，保存后在cmd里面输入cnv_encd(G,k0,input)，就没有问题了。

>> cnv_encd(a,2,0)

ans =

0 0 0