%程序说明：对一个二进制序列进行hdb3编码运算
%function y=hdb3(x)
x=[0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 1 1 0 0 0 0 ]
n=length(x);%输入系列的长度
last_V=-1;%初始化V的符号
last_one=-1;%初始化1的符号
y=zeros(size(x)); %初始化输出序列
count=0; %连0计数器
for i=1:n
if x(i)==1 %遇1则极性反转
y(i)=-last_one;
last_one=y(i);
count=0;
else
count=count+1;%遇0则零的个数加1
if count==4 %遇到连续4个0
count=0;
y(i)=-last_V;%遇4连零则置为V，相邻的V极性反转
last_V=y(i);
if y(i)*last_one==-1 %V与前一个非0符号必须同极性（相邻V之间如果间隔偶数个1，则y(i)*last_one==-1，插入B符号，否者，不插入。）
y(i-3)=y(i);
end
last_one=y(i);
end
end
end
y
