insert_zero.asv

来自「convolution coding program using viterab」· ASV 代码 · 共 141 行

function [data_insert_zero] = insert_zero (data_soft_de, rate, hard_soft,
fl)
% check the sizes
if hard_soft == 1
p1 = length(data_soft_de);
if rate == 1
data_insert_zero = data_soft_de;
elseif rate == 2
if rem(p1, 3) == 0
p2 = p1/3;
for i = 1: p2
data_insert_zero(1 + (i-1)*4) = data_soft_de(1 + (i-1)*3);
data_insert_zero(2 + (i-1)*4) = data_soft_de(2 + (i-1)*3);
data_insert_zero(3 + (i-1)*4) = data_soft_de(3 + (i-1)*3);
data_insert_zero(4 + (i-1)*4) = 3;
end
else
p3 = rem(p1, 3);
p2 = (p1 - p3)/3;
for i = 1: p2
data_insert_zero(1 + (i-1)*4) = data_soft_de(1 + (i-1)*3);
data_insert_zero(2 + (i-1)*4) = data_soft_de(2 + (i-1)*3);
data_insert_zero(3 + (i-1)*4) = data_soft_de(3 + (i-1)*3);
data_insert_zero(4 + (i-1)*4) = 3;
end
for i= 1: p3
data_insert_zero(4*p2 + i) = data_soft_de(3*p2 + i);
end
end
elseif rate == 3
if rem(p1, 4)~= 0
p3 = rem(p1, 4);
p2 = (p1 - p3)/4;
for i = 1: p2
data_insert_zero(1 + (i-1)*6) = data_soft_de(1 + (i-1)*4);
data_insert_zero(2 + (i-1)*6) = data_soft_de(2 + (i-1)*4);
data_insert_zero(3 + (i-1)*6) = 3;
data_insert_zero(4 + (i-1)*6) = data_soft_de(3 + (i-1)*4);
data_insert_zero(5 + (i-1)*6) = 3;
data_insert_zero(6 + (i-1)*6) = data_soft_de(4 + (i-1)*4);
end
for i= 1: p3
data_insert_zero(6*p2 + i) = data_soft_de(4*p2 + i);
end
elseif rem(p1, 4) == 0 && fl == 0
p2 = p1/4;
for i = 1: p2
data_insert_zero(1 + (i-1)*6) = data_soft_de(1 + (i-1)*4);
data_insert_zero(2 + (i-1)*6) = data_soft_de(2 + (i-1)*4);
data_insert_zero(3 + (i-1)*6) = 3;
data_insert_zero(4 + (i-1)*6) = data_soft_de(3 + (i-1)*4) ;
data_insert_zero(5 + (i-1)*6) = 3;
data_insert_zero(6 + (i-1)*6) = data_soft_de(4 + (i-1)*4);
end
elseif rem(p1, 4) == 0 && fl == 1
p2 = p1/4;
for i = 1: (p2-1)
data_insert_zero(1 + (i-1)*6) = data_soft_de(1 + (i-1)*4);
data_insert_zero(2 + (i-1)*6) = data_soft_de(2 + (i-1)*4);
data_insert_zero(3 + (i-1)*6) = 3;
data_insert_zero(4 + (i-1)*6) = data_soft_de(3 + (i-1)*4) ;
data_insert_zero(5 + (i-1)*6) = 3;
data_insert_zero(6 + (i-1)*6) = data_soft_de(4 + (i-1)*4);
end
data_insert_zero((p2-1)*6+1) = data_soft_de((p2-1)*4+1);
data_insert_zero((p2-1)*6+2) = data_soft_de((p2-1)*4+2);
data_insert_zero((p2-1)*6+3) = data_soft_de((p2-1)*4+3);
data_insert_zero((p2-1)*6+4) = data_soft_de((p2-1)*4+4);
end
else error('The setup of rate is out of range!!!')
end
elseif hard_soft == 0
p1 = length(data_soft_de);
if rate == 1
data_insert_zero = data_soft_de;
elseif rate == 2
if rem(p1, 3) == 0
p2 = p1/3;
for i = 1: p2
data_insert_zero(1 + (i-1)*4) = data_soft_de(1 + (i-1)*3);
data_insert_zero(2 + (i-1)*4) = data_soft_de(2 + (i-1)*3);
data_insert_zero(3 + (i-1)*4) = data_soft_de(3 + (i-1)*3);
41
data_insert_zero(4 + (i-1)*4) = 0;
end
else
p3 = rem(p1, 3);
p2 = (p1 - p3)/3;
for i = 1: p2
data_insert_zero(1 + (i-1)*4) = data_soft_de(1 + (i-1)*3);
data_insert_zero(2 + (i-1)*4) = data_soft_de(2 + (i-1)*3);
data_insert_zero(3 + (i-1)*4) = data_soft_de(3 + (i-1)*3);
data_insert_zero(4 + (i-1)*4) = 0;
end
for i= 1: p3
data_insert_zero(4*p2 + i) = data_soft_de(3*p2 + i);
end
end
elseif rate == 3
if rem(p1, 4)~= 0
p3 = rem(p1, 4);
p2 = (p1 - p3)/4;
for i = 1: p2
data_insert_zero(1 + (i-1)*6) = data_soft_de(1 + (i-1)*4);
data_insert_zero(2 + (i-1)*6) = data_soft_de(2 + (i-1)*4);
data_insert_zero(3 + (i-1)*6) = 0;
data_insert_zero(4 + (i-1)*6) = data_soft_de(3 + (i-1)*4);
data_insert_zero(5 + (i-1)*6) = 0;
data_insert_zero(6 + (i-1)*6) = data_soft_de(4 + (i-1)*4);
end
for i= 1: p3
data_insert_zero(6*p2 + i) = data_soft_de(4*p2 + i);
end
elseif rem(p1, 4) == 0 & fl == 0
p2 = p1/4;
for i = 1: p2
data_insert_zero(1 + (i-1)*6) = data_soft_de(1 + (i-1)*4);
data_insert_zero(2 + (i-1)*6) = data_soft_de(2 + (i-1)*4);
data_insert_zero(3 + (i-1)*6) = 0;
data_insert_zero(4 + (i-1)*6) = data_soft_de(3 + (i-1)*4) ;
data_insert_zero(5 + (i-1)*6) = 0;
data_insert_zero(6 + (i-1)*6) = data_soft_de(4 + (i-1)*4);
end
elseif rem(p1, 4) == 0 && fl ~= 0
p2 = p1/4;
for i = 1: p2-2
data_insert_zero(1 + (i-1)*6) = data_soft_de(1 + (i-1)*4);
data_insert_zero(2 + (i-1)*6) = data_soft_de(2 + (i-1)*4);
data_insert_zero(3 + (i-1)*6) = 0;
data_insert_zero(4 + (i-1)*6) = data_soft_de(3 + (i-1)*4) ;
data_insert_zero(5 + (i-1)*6) = 0;
data_insert_zero(6 + (i-1)*6) = data_soft_de(4 + (i-1)*4);
end
data_insert_zero((p2-2)*6+1) = data_soft_de((p2-1)*4+1);
data_insert_zero((p2-2)*6+2) = data_soft_de((p2-1)*4+2);
data_insert_zero((p2-2)*6+3) = data_soft_de((p2-1)*4+3);
data_insert_zero((p2-2)*6+4) = data_soft_de((p2-1)*4+4);
end
else error('The setup of rate is out of range!!!')
end
end