MATLAB实现HDB3编码

function [code,varargout]=HDB3(signal,varargin)
%-------------------------使用须知-------------------------------
%功能   ,先得到AMI码,即把消息码中的"1"交替变成"-1"和"+1"将"0"仍保持为"0"
%        然后检查AMI码中连0的情况,
%            若AMI码中连0的个数小于4,此时的AMI码就是HDB3码;
%            若AMI码中连0的个数大于4,则将每4个连0小段的第4个0变换成与前一个非0符号(+1或-1)同极性的符号,
%            用V表示(+1 === +V,-1 === -V);
%
%使用   ,example:
%       s=[0 1 0 1 1 0 0 0 1 0 0 0 0 1 1 1 1 1 1];
%       code=HDB3(s);                        %全部按默认设置,返回二进制序列码
%       code=HDB3(s,2);                      %初始连"0"个数为2,返回二进制序列码
%       code=HDB3(s,2,-1);                   %初始连"0"个数为2,信号码序中的第一个"1"设为-1,返回二进制序列码和最后一个"1"的状态
%       [code,pre]=HDB3(s,-1);               %信号码序中的第一个"1"设为-1,返回二进制序列码和最后连"0"个数
%       [code,pre,status]=HDB3(s,2);         %初始连"0"个数为2,返回二进制序列码,最后连"0"个数
%                                            %和最后一个"1"的状态
%       [code,pre,status]=HDB3(s,2,-1);      %初始连"0"个数为2,信号码序中的第一个"1"设为-1,返回二进制序列码,
%                                            %最后连"0"个数和最后一个"1"的状态
%输入参数
%signal ,表示输入信号
%pre    ,初始连"0"个数,默认值0
%        最终存储信号码序中的最后连"0"个数,并可返回供其它程序使用
%status ,表示信号码序中的第一个"1"的状态,默认值1
%        最终存储信号码序中的最后一个"1"的状态,并可返回供其它程序使用
%输出参数
%code   ,返回变换后的二进制序列串
%pre
%status
%局部变量
%sf     ,16倍输入信号码速