matlab-(含教程)基于MATLAB的HDB3编译码仿真,输出AMI码,AMI插V码,HDB3码以及解码

clc; clear; close all; %使每次产生的随机数是相同的。 rand('state',3) %a表示元素个数 len=20; %随机产生a-1个数据 rand_number=unidrnd(2,1,len-1); %产生01随机数列 rand_number2=rand_number-1; rand_number2=[1,rand_number2]; data1=rand_number2; disp('HDB3编码前的原始数据');data1 %首先计算AMI码 flag=1; for n=1:len if rand_number2(n)==1 ; flag = -flag; rand_number2(n)= flag; end end AMI = rand_number2; data2 = AMI; disp('AMI码');data2 m=-1; for k=1:len-3 if abs(rand_number2(k))==1 m=rand_number2(k); end %检测是否连续4个码元都是0 if rand_number2(k)==0 if rand_number2(k+1)==0 if rand_number2(k+2)==0 if rand_number2(k+3)==0 %用2表示V rand_number2(k+3)=2*m; end end end end end AMIV = rand_number2; data20 = AMIV; disp('AMI码插V后的值');data20 rand_number2 ; flag=1; for s=1:len if abs(rand_number2(s))==2 for d=s+1:len %V之间的非零码只有+1，-1遇到一个一使l变一次号，表示偶数或奇数 if abs(rand_number2(d))==1 flag=-flag; end if abs(rand_number2(d))==2 if flag~=1 flag=1; else %检测到下一个V时，若为偶数个，则插入B %3代表B rand_number2(d-3)=3*(-(rand_number2(d-4))/(abs(rand_number2(d-4)))); rand_number2(d:len)=-rand_number2(d:len); end end end end end HDB = rand_number2; data3 = HDB; disp('HDB3码');data3 %解码 for k=1:len if abs(rand_number2(k))==2 rand_number2(k-3)=0;rand_number2(k)=0; end if abs(rand_number2(k))==1 rand_number2(k)=1; end end %解码后的恢复序列 data4 = rand_number2; disp('HDB3解码后的值');data4 figure; subplot(511); stem(data1);title('HDB3编码前的原始数据'); subplot(512); stem(data2);title('AMI码'); subplot(513); stem(data20);title('AMI码插V后的值'); subplot(514); stem(data3);title('HDB3编码前的原始数据'); subplot(515); stem(data4);title('HDB3解码后的值');