一个卷积码编码程序

newbabyborn

function  output = cnv_encd(g,k0,input)
%   cnv_encd(g,k0,input)
%   determines the output sequence of a binary convolutional encoder
%   g is the generator matric of the convolutional code
%    with n0 rows and l*k0 columns, Its rows are g1,g2….gn
%    k0 is the number of bus entering the encoder at each clock cycle.
%    input The binary input seq
g = [0 0 1 0 0 1 0 0; 0 0 0 0 0 0 0 1; 1 0 0 0 0 0 0 1; 0 1 0 0 1 1 0 1];
k0 = 2;
input=[10011100110000111];
Output=cnv_encd(g,k0,input);
% check to see if extra zero padding is necessary    看填充的0是否必须
if rem(length(input), k0) > 0
   input = [input, zeros(size(1:k0-rem(length(input),k0)))];
end
n = length(input)/k0;          % 把输入比特按k0分组，n为所得的组数。
% check the size of matrix g
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，以保证编码器是从全0开始，并回到全0状态。
u = [zeros(size(1L-1)*k0)), input, zeros(size(1L-1)*k0))];
% generate 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));
len_tal = n0*(L + n - 1);
% write the output to the encodetext
result = fopen(encodetext, 'w');
for i = 1:n0*(L+n -1)
    fwrite(result, output(i), 'bit1');
end
fclose(result);



这是我的一个卷积编码程序，可是无法运行，也没有错误提示，请高人指点，帮忙改下程序！小妹在此谢过了！o(∩_∩)o...

[ 本帖最后由 newbabyborn 于 2008-5-13 15:27 编辑 ]

superhero

good!good!

tw178103556

一定要顶

fishdf

wo ding

fishdf

xue xi

jianzhiqi

ding~~~~~~~~~

bashenan

好东西啊，谢谢了啊