%HDB3prog.m
%本函数实现将输入的一段二进制代码进行HDB3编译码
xn=[1 0 1 1 0 0 0 0 0 0 0 1 1 0 0 0 0 0 0 1 0]; %输入单极性码
yn=xn; %输出yn初始化
num=0; %计数器初始化
for k=1:length(xn)
if xn(k)==1
num=num+1; %“1”计数器
if num/2== fix(num/2) %奇数个1时输出-1，进行极性交替
yn(k)=1;
else
yn(k)=-1;
end
end
end
%HDB3编码
num=0; %连0计数器初始化
yh=yn; %输出初始化
sign=0; %极性标致初始化为0
V=zeros(1,length(yn)); %V脉冲位置记录变量
B=zeros(1,length(yn)); %B脉冲位置记录变量
for k=1:length(yn)
if yn(k)==0
num=num+1; %连“0”个数计数
if num== 4 %如果4连0
num=0; %计数器清零
yh(k)=1*yh(k-4);
%让0000的最后一个0改变为与前一个非零符号相同极性的符号
V(k)=yh(k); %V脉冲位置记录
if yh(k)==sign %如果当前V符号与前一V符号的极性相同
yh(k)=-1*yh(k); %则让当前V符号极性反转，以满足V符号间相互极性反转要求
yh(k-3)=yh(k); %添加B符号，与V符号同极性
B(k-3)=yh(k); %B脉冲位置记录
V(k)=yh(k); %V脉冲位置记录
yh(k+1:length(yn))=-1*yh(k+1:length(yn));%让后面的非零符号从V符号开始再交替变化
end
sign=yh(k); %记录前一个V符号的极性
end
else
num=0; %当前输入为1，则连0计数器清零
end
end %编码完成
re=[xn',yn',yh',V',B']; %编码结果输出
%HDB3译码
input=yh; %HDB3码输入
decode=input; %输入初始化
sign=0; %极性标志初始化
for k=1:length(yh)
if input(k)~=0
if sign==yh(k) %如果当前码与前一个非零符号的极性相同
decode(k-3:k)=[0 0 0 0]; %则该码判为V码并将*000V清零
end
sign=input(k); %极性标志
end
end
decode=abs(decode); %整流
error=sum([xn'-decode']); %解码的正确性检验，作图
subplot(3,1,1);stairs([0:length(xn)-1],xn);axis([0 length(xn) -2 2]);
subplot(3,1,2);stairs([0:length(xn)-1],yh);axis([0 length(xn) -2 2]);
subplot(3,1,3);stairs([0:length(xn)-1],decode);axis([0 length(xn) -2 2]);
