puncture.m

国外大学的同志作出来的好程序

function outbits=puncture(inbits, punct_vec);
%
%Puncture the parity bits via punct_vec.
%
%
punct_blen=length(punct_vec);

%zero-padding
tval=mod(length(inbits), punct_blen);
if tval,
    inbits=[inbits; zeros(punct_blen-tval, 1)];
end

%reshape inbits into matrix
punct_blks=length(inbits)/punct_blen;
inbits_mat=reshape(inbits, punct_blen, punct_blks);

%get the punctured matrix.
punct_index=find(punct_vec~=0);
outbits_mat=inbits_mat(punct_index, :);

outbits=outbits_mat(:);