clear

clc

a=1024;                           %a表示元素个数，可以修改码元为任意个数如1024

r=unidrnd(2,1,a-1);              %随机产生离散均匀二值分布(1，2)，共a-1个

rn=r-1;                          %每个元素减1变为(0，1)二值随机数列

rn=[1,rn];                        %跟书上一致，从1开始随机产生离散均匀二值分布(1，2)，共a个

ori=rn;                           %将原序列保存起来，便于后面解码后比较

ori

l=1;

for n=1:a

if rn(n)==1 ;

l=-l;                            %每检测到一个1变一次符号

rn(n)=l;

end

end

AMI=rn;                             %变为AMI码：每过一个1就使原来的1变一次号

AMI

m=-1;                           %符号初始化

for k=1:a-3

if abs(rn(k))==1

m=rn(k);                 %m记录可能存在的连续4个0前面非零码元的符号

end

if rn(k)==0

if rn(k+1)==0

if rn(k+2)==0

if rn(k+3)==0  %检测是否连续4个码元都是0

rn(k+3)=2*m; %用2表示书上的V

end

end

end

end

end

rn ;                                 %插入V(2)后   checked

l=1;                                 %l=1表示两个V之间有偶数个非零码，0是偶数，所以初始为1，l=-1表示两个V之间有奇数个非零码

for s=1:a

if abs(rn(s))==2

for d=s+1:a

if abs(rn(d))==1           %V之间的非零码只有+1，-1遇到一个一使l变一次号，表示偶数或奇数

l=-l;

end

if abs(rn(d))==2

if l~=1

l=1;                       %不需要插入B时若检测到V也应把计数器清零

%最初版本的升级处1

else                           %检测到下一个V时，若为偶数个，则插入B

rn(d-3)=3*(-(rn(d-4))/(abs(rn(d-4))));               %3代表B,插入不带符号的B

rn(d:a)=-rn(d:a);

end                           %V后的符号再交替

end

end

end

end

rn                                      %到此处完成了插入不带符号的B

%for s=1:a

% if abs(rn(s))==3                                %找到B

%rn(s)=rn(s)*(-(rn(s-1))/(abs(rn(s-1))));           %B(3)与前一个位置的带符号的归一值相乘再取反，实现符号B(3)的极性与前一非符号的相反

%B后面第三个就是接下来的V,从它开始非零码正负号交替变化

%end

%end

HDB=rn;                                        %给B(3)添加了符号，并且实现了V后的符号再交替

HDB ;

%以上便实现了HDB3的编码

%下面进行解码

for k=1:a

if abs(rn(k))==2

rn(k-3)=0;rn(k)=0;                       %每个V(2)都变成0，V前面第三个有可能是B(3)有可能是0也恢复为0

end

if abs(rn(k))==1

rn(k)=1;

end

end

rn;                                            %解码后的恢复序列

rn-ori                                         %解码与原码比较全为0则解码正确