matlabHDB3编码译码数字信号调制解调

《matlabHDB3编码译码数字信号调制解调》由会员分享，可在线阅读，更多相关《matlabHDB3编码译码数字信号调制解调（13页珍藏版）》请在装配图网上搜索。

1、、HDB3码的编码和译码1、实验要求：掌握HDB3码的编码规则，利用MATLAB设计并实现HDB3码的编码和译码。2、原理简述：编码规则:1）先将消息代码变换成AMI码，若AMI码中连0的个数小于4,此时的AMI码就是HDB3码；2）若AMI码中连0的个数大于3,则将每4个连0小段的第4个0变换成与前一个非0符号（+1或-1）同极性的符号,用表示（+1+,-1-）；3）为了不破坏极性交替反转，当相邻符号之间有偶数个非0符号时，再将该小段的第1个0变换成+B或-B,符号的极性与前一非零符号的相反，并让后面的非零符号从符号开始再交替变化。例如:消息代码:100001000011000011AMI码

2、：+10000-10000+1-10000+1-1HDB3码：+1000+V-1000-V+1-1+B00+V-1+1简易编码方法：原理：HDB3码既要包含AMI的交替特性使输出无直流特性，又要不出现四个以上的连0，因此可以先满足后者。1）把0000换为取代节。规则：先将0000分离开来，第一个0000直接变为000V,然后数相邻两个0000之间1的个数，奇数则变为000V,偶数则变为B00V。2）更新符号。根据教材有：B总是与其前面的1或V符号相反，V总是与前面的1或B相符号相同,1总是与前面的V或B符号相反，就可以编符号了。例如：消息代码：1000010000110000110000000

3、011100001中间码：1000V1000V11B00V11B00VB00V111000V1HDB3码：+1000+V-1000-V+1-1+B00+V-1+1-B00-V+B00+V-1+1-1000-V+1解码规则：1）虽然编码很复杂，但解码规则很简单，若3连“0”前后非零脉冲同极性，则将最后一个非零元素译为零，如+1000+1就应该译成“10000”；若2连“0”前后非零脉冲极性相同，则两零前后都译为零，如-100-1，就应该译为0000.2）再将所有的-1变换成+1后，就可以得到原消息代码。3、程序运行结果运行程序可自定义二进制序列，第一个窗口显示所输入的二进制码元，第二个窗口为HD

4、B3编码后的结果，最后一个窗口为解码后的窗口。4、程序源代码%ch7example6prog1.m%AMI码的编码xn=101100000001100000010;%输入单极性码yn=xn;%输出yn初始化num=0;%计数器初始化fork=1:length(xn)ifxn(k)=1num=num+1;%1计数器ifnum/2=fix(num/2)%奇数个1时输出-1,进行极性交替yn(k)=1;elseyn(k)=-1;endendend%HDB3编码num=0;%连零计数器初始化yh=yn;%输出初始化sign=0;%极性标志初始化为0V=zeros(1,length(yn);%V脉冲位置

5、记录变量B=zeros(1,length(yn);%B脉冲位置记录变量fork=1:length(yn)ifyn(k)=0num=num+1;%连“0”个数计数ifnum=4%如果4连“0”num=0;%计数器清零yh(k)=1*yh(k-4);%让0000的最后一个0改变为与前一个非零符号相同极性的符号V(k)=yh(k);ifyh(k)=sign%V脉冲位置记录%如果当前V符号与前一个V符号的极性相同yh(k)=-1*yh(k);%则让当前V符号极性反转，以满足V符号间相互极性反转要求yh(k-3)=yh(k);%添加B符号，与V符号同极性B(k-3)=yh(k);%B脉冲位置记录V(k)

6、=yh(k);%V脉冲位置记录yh(k+1:length(yn)=-1*yh(k+1:length(yn);%并让后面的非零符号从V符号开始再交替变化sign=yh(k);endendelsenum=0;endendre=xn,yn,yh,V,B;%当前输入为“1”则连“0”计数器清零%编码完成%结果输出：xnAMIHDB3V&B符号%HDB3解码input=yh;decode=input;sign=0;fork=1:length(yh)%HDB3码输入%输出初始化%极性标志初始化ifinput(k)=0ifsign=yh(k)%如果当前码与前一个非零码的极性相同%记录前一个V符号的极性dec

7、ode(k-3:k)=0000;%则该码判为V码并将*00V清零endsign=input(k);%极性标志endenddecode=abs(decode);%整流error=sum(xn-decode);%解码的正确性检验,作图subplot(3,1,1);stairs(0:length(xn)-1,xn);axis(0length(xn)-22);subplot(3,1,2);stairs(0:length(xn)-1,yh);axis(0length(xn)-22);subplot(3,1,3);stairs(0:length(xn)-1,decode);axis(0length(xn)

8、-22);二、数字带通系统的调制解调1、实验要求：利用MATLAB程序设计语言实现数字带通系统的调制解调。2、原理简述：数字调制：用数字基带信号控制载波，把数字基带信号变换为数字带通信号(已调信号)的过程成为数字调制。数字调制与模拟调制的基本原理相同，但数字信号有离散取值的特点。因此数字调制技术有两种方法：利用模拟调制的方法实现数字调制，即把数字调制看成是模拟调制的特例把数字基带信号当成是模拟信号的特殊情况处理;利用数字信号的离散取值特点通过开关键控制载波，从而实现数字调制。这种方法称为键控法，比如对载波的振幅，频率和相位进行键控，便可得到振幅键控ASK,频移键控FSK和相移键控PSK三种基本

9、的数字调制方式。振幅键控ASK振幅键控利用载波的幅度变化来传递数字信息，而其频率和初始相位保持不变。在2ASK中，载波的幅度只有两种变化状态，分别为对应二进制信息“0”或“1”。常用的二进制振幅键控方式称为通断键控，2ASK信号的一般表达式为s(E)=口(f)沁他I其中2ASK有两种基本的解调方法：非相干解调和相干解调，解调方式如下：非相干解调方式：带通滤波器全波整流器低通滤波器抽样判决器输出想干解调方式：带通滤波器相乘器低通滤波器抽样判决器输出频移键控FSK频移键控是利用载波的频率变化来传递数字信息。在2FSK中，载波的频率随二进制基带信号在fl和f2两个频率点检变化。一个2FSK信号可以看

10、成是两个不同载频的2ASK信号的叠加。其产生方式有两种：采用模拟调频来实现也可采用键控法实现，即在二进制基带矩形序列的控制下通过开关电路对两个不同的独立频率源进行选通。但相邻码元之间的相位不一定连续。2FSK信号的常用解调方法也采用相干解调和非相干解调两种方式。其解调原理是将2FSK信号分解为上下两路2ASK信号分别进行解调，然后进行判决。相移键控PSK二进制相移键控2PSK相移键控是利用载波的相位变化来传递数字信息，而振幅和频率保持不变。在2PSK中，通常用初始相位0和n分别表示二进制“1”和“0”。与2ASK信号的产生方式相比，只是对s(t)的要求不同，在2ASK中是单极性的，而在2PSK

11、中是双极性的基带信号。2PSK信号的解调通常采用相干解调法。二进制差分相移键控2DPSK2PSK相干解调时，由于载波恢复中相位有0、n模糊性，导致解调过程出现“反向工作”现象，恢复出的数字信号“1”“0”倒置，使得2PSK难以实际应用。而2DPSK是利用前后相邻码元的载波相对相位变化传递数字信息，对于相同的基带数字信息序列，由于初始相位不同，2DPSK的相位可以不同。也就是说，2DPSK信号的相位并不直接代表基带信号，而前后码元相对相位的差才唯一决定信号符号。2DPSK信号的解调方法为相干解调(极性比较法)加码反变换法。其原理为：带通滤波器相乘器低通滤波器抽样判决器码反变换输出3、程序运行结果

12、2ASK调制后的信号以及2ASK解调后的信号D1D2FSK调制后的信号以及使用2FSK解调后的信号原诟诣号谓制牯的信号随机产生得二进制原始信号以及2DPSK调制之后得信号使用2DPSK解调之后的信号%调制序列的长度%载波频率%时域采样间隔%频率分辨率%设每个码元长度为1ms时间矢量：0：N%产生调制序列4、程序源代码2ASKclearall;clear;N=10;fc=20;datat=0.01;data=1/datat;df=0.01;t=0:datat:N;m=randint(1,N);fs=1/datat;c=cos(2*pi*fc*t);fori=1:length(m)if(m(i)=

13、1)forj=1:1/datat;ask_pre(j+(i-1)*(1/datat)=1;endelseif(m(i=0)ask_pre(j+(i-1)*(1/datat)=0;endendask_pre=ask_pre,m(i);ask=ask_pre.*c;ask_F,ask,df1=fftseq(ask,datat,df);C,c,df1=fftseq(c,datat,df);M,ask_pre,df1=fftseq(ask_pre,datat,df);f=0:df1:df1*(length(c)-1)-fs/2;subplot(2,1,1);plot(t,ask_pre(1:lengt

14、h(t);axis(0,N,-1.5,1.5);xlabel(时间ms);title(未调信号);ylabel(幅值：v);subplot(2,1,2);plot(t,ask(1:length(t);axis(0,N,-1.5,1.5);xlabel(时间ms);title(2ASK信号);ylabel(幅值：v);%解调%ask_n=(ask_pre+rand(1,length(ask_pre).*c;%(高斯白噪声)ask_NF,ask_n,df1=fftseq(ask_n,datat,df);f_strain=8;n_off=floor(f_strain/df1);f=0:df1:df1

15、*(length(ask)-1)-fs/2;n0=floor(fc-8)/df1):floor(fc+8)/df1);BP=zeros(1,length(f);BP(min(n0):max(n0)=ones(1,length(n0);%带通滤波器设计BP(length(BP)-max(n0):length(BP)-min(n0)=ones(1,length(n0);H=zeros(1,length(f);H(1:n_off)=2*ones(1,n_off);%矩形窗截取H(length(f)-n_off+1:length(f)=2*ones(1,n_off);%H(1:n_off)=2*ham

16、ming(n_off);%hamning窗截取%H(length(f)-n_off+1:length(f)=2*hamming(n_off);BP1=BP.*ask_F;BPN=BP.*ask_NF;ask_re1=ifft(BP1);ask_ren1=ifft(BPN);ask_re=2*ask_re1.*c;%相干解调法ask_ren=2*ask_ren1.*c;%相干解调(高斯白噪声)ask_reF,ask_re,df1=fftseq(ask_re,datat,df);ask_reNF,ask_ren,df1=fftseq(ask_ren,datat,df);D1=H.*ask_reF;

17、DN=H.*ask_reNF;d1=ifft(D1);dn=ifft(DN);forlength(m)if(d1(i-1)*(1/datat)+(1/(2*datat)=0.5)决规则:d10.5=1m_de(i)=1;d10forj=1:1/datatd11(j+(i-1)*(1/datat)=1;endelsem_de(i)=0;forj=1:1/datatd11(j+(i-1)*(1/datat)=0;endendendforlength(m)if(dn(i-1)*(1/datat)+(1/(2*datat)=0.5)决规则:d10.5=1m_den(i)=1;d10forj=1:1/d

18、atatdnn(j+(i-1)*(1/datat)=1;endelsem_den(i)=0;forj=1:1/datatdnn(j+(i-1)*(1/datat)=0;endi=1:%判决%判%为了画解调波形需要%为了画解调波形需要%否则i=1:%判决%判%为了画解调波形需要%为了画解调波形需要%否则endendfigure(2);subplot(4,1,1);plot(t,ask(1:length(t);axis(0,N,-1.5,1.5);xlabel(时间);title(解调器接收波形)；subplot(4,1,2);plot(t,ask_re1(1:length(t);axis(0,N

19、,-1.5,1.5);xlabel(时间);title(经过带通滤波器波形);subplot(4,1,3);plot(t,d1(1:length(t);axis(0,N,-0.5,2.5);xlabel(时间);title(经过低通滤波器波形);subplot(4,1,4);plot(t(1:length(d11),d11);axis(0,N,-1.5,1.5);xlabel(时间);title(经过判决器波形);2FSK%调制序列的长度%载波基频clear;N=20;fc=5;datat=0.01;fs=1/datat;df=0.01t=0:datat:N;m=randint(1,N);fo

20、ri=1:length(m)if(m(i)=1)forj=1:1/datat;fsk_f(j+(i-1)*(1/datat)=4;mm(j+(i-1)*(1/datat)=1;endelseif(m(i)=0)forj=1:1/datat;fsk_f(j+(i-1)*(1/datat)=1;mm(j+(i-1)*(1/datat)=0;endendendif(m(i)=1)nn=4;elsenn=1;endfsk_f=fsk_f,nn;fsk_c=fsk_f*fc;c=cos(2*pi*fsk_c.*t);C,c,df1=fftseq(c,datat,df);MM,mm,df1=fftseq(

21、mm,datat,df);f=0:df1:df1*(length(mm)-1)-fs/2;%时域采样间隔%频率分辨率%频率缩放系数%设每个码元长度为1ms时间矢量：0：N%产生调制序列%1码的载波基波倍频系数20hz%为了画码元时域波形需要%0码的载波基波倍频系数5hz%为了画码元时域波形需要%使fsk_f的长度与t相同%根据不同的倍频系数产生不同的载波频率%2fsk信号产生%2fsk信号的傅立叶变换%基带信号的傅立叶变换%频率作图矢量figure(1);subplot(2,1,1);plot(t,mm(1:length(t);axis(0,N,-1.5,1.5);xlabel(时间ms);t

22、itle(未调信号);ylabel(幅值：v);subplot(2,1,2);plot(t,c(1:length(t);axis(0,N,-1.5,1.5);xlabel(时间ms);title(2FSK信号);ylabel(幅值：v);解调%B0=fc-2,fc+2;%0码频带宽B1=fc*4-2,fc*4+2;%1码频带宽n0=floor(fc-2)/df1):floor(fc+2)/df1);n1=floor(4*fc-2)/df1):floor(4*fc+2)/df1);H0=zeros(1,length(C);H0(min(n0):max(n0)=ones(1,length(n0);

23、%滤波器设计(0码带通滤波器)H0(length(H0)-max(n0):length(H0)-min(n0)=ones(1,length(n0);H1=zeros(1,length(C);H1(min(n1):max(n1)=ones(1,length(n1);%滤波器设计(1码带通滤波器)H1(length(H1)-max(n1):length(H1)-min(n1)=ones(1,length(n1);n_off0=fc/df1;n_off1=4*fc/df1;L0=zeros(1,length(C);%(0码低通滤波器)L0(1:n_off0)=2*ones(1,n_off0);L0(

24、length(L0)-n_off0+1:length(L0)=2*ones(1,n_off0);L1=zeros(1,length(C);%(1码低通滤波器)L1(1:n_off1)=2*ones(1,n_off1);L1(length(L1)-n_off1+1:length(L1)=2*ones(1,n_off1);C0=C.*H0;C1=C.*H1;m0=ifft(C0);m1=ifft(C1);mm0=abs(m0);%包络检波mm1=abs(m1);MM0,mm0,df1=fftseq(mm0,datat,df);MM1,mm1,df1=fftseq(mm1,datat,df);MM0

25、=MM0.*L0;MM1=MM1.*L1;mm0=ifft(MM0);mm1=ifft(MM1);fori=1:length(m)%比较器判决if(abs(m0(i-1)*(1/datat)+(1/(2*datat)=abs(m1(i-1)*(1/datat)+(1/(datat*2)%判决规贝y：mom1=0m_de(i)=0;%否则=1forj=1：1/datatm_anti(j+(i-1)*(1/datat)=0;%为了画解调波形需要endelsem_de(i)=1;forj=1：1/datatm_anti(j+(i-1)*(1/datat)=1;%为了画解调波形需要endendendm

26、_anti=m_anti,m_de(i);figure(2)subplot(5,1,1);plot(t,m0(1:length(t);xlabel(时间);title(经过带通滤波器1后波形);subplot(5,1,2);plot(t,m1(1:length(t);xlabel(时间);title(经过带通滤波器2后波形);subplot(5,1,3);plot(t,mm0(1:length(t);xlabel(时间);axis(0,N,-0.1,2);title(经过低通滤波器1后波形);subplot(5,1,4);plot(t,mm1(1:length(t);xlabel(时间);ax

27、is(0,N,-0.1,2);title(经过低通滤波器2后波形);subplot(5,1,5);plot(t,m_anti(1:length(t);xlabel(时间);axis(0,N,-0.5,1.5);title(经过比较起判决后波形);2PSKclear;N=10;%调制序列的长度fc=10;%载波频率datat=0.01;%时域采样间隔fs=1/datat;%频率分辨率df=0.01%频率缩放系数t=0:datat:N;%设每个码元长度为1ms时间矢量：0：Nm=randint（1,N）;%产生调制序列%调制%fori=1:length(m)if(m(i)=1)forj=1:1/d

28、atat;psk_p(j+(i-1)*(1/datat)=pi;mm(j+(i-1)*(1/datat)=1;end%1码的对应的相位为pi（采用a方式）%为了画码元时域波形需要endelseif(m(i)=0)forj=1:1/datat;psk_p(j+(i-1)*(1/datat)=0;mm(j+(i-1)*(1/datat)=0;end%0码的对应的相位为0（采用a方式）%为了画码元时域波形需要endif(m(i)=1)nn=pi/2;elsenn=-(pi/2);endpsk_p=psk_p,nn;调制相位c=cos(2*pi*fc*t+psk_p);C,c,df1=fftseq(c

29、,datat,df);MM,mm,df1=fftseq(mm,datat,df);f=0:df1:df1*(length(mm)-1)-fs/2;%使p的长度与t相同,psk_p为根据不同码元产生不同的%2psk信号产生%2fsk信号的傅立叶变换%基带信号的傅立叶变换%频率作图矢量subplot(2,1,1);plot(t,mm(1:length(t);axis(0,N,-1.5,1.5);xlabel(时间ms);title(未调信号);ylabel(幅值：v);subplot(2,1,2);plot(t,c(1:length(t);axis(0,N,-1.5,1.5);xlabel(时间m

30、s);title(2PSK信号);ylabel(幅值：v);%解调carwav=2*cos(2*pi*fc.*t);CF,carwav,df1=fftseq(carwav,datat,df);n_off=(fc-1)/df;n0=floor(fc-2)/df1):floor(fc+2)/df1);BP=zeros(1,length(f);BP(min(n0):max(n0)=ones(1,length(n0);%带通滤波器设计BP(length(BP)-max(n0):length(BP)-min(n0)=ones(1,length(n0);H=zeros(1,length(f);H(1:n_

31、off)=2*ones(1,n_off);%低通滤波器设计(矩形窗截取)H(length(f)-n_off+1:length(f)=2*ones(1,n_off);CC=C.*BP;%接受时经过一个带通滤波器cc=ifft(CC);%反福利也变换c_re=cc.*carwav;%乘以载波进行相干解调cc1=c_re;c_reF,c_re,df1=fftseq(c_re,datat,df);c_reF=c_reF.*H;%经过低通滤波器c_ref=ifft(c_reF);forlength(m)if(c_ref(i-1)*(1/datat)+(1/(2*datat)0=0m_de(i)=1;1f

32、orj=1:1/datatc_de(j+(i-1)*(1/datat)=1;endelsem_de(i)=0;forj=1:1/datatc_de(j+(i-1)*(1/datat)=0;endendendc_de=c_de,c_re(i);figure(2);subplot(5,1,1);plot(t,c(1:length(t);axis(0,N,-2,2);i=1:%判决%判决%为了画解调波形需要%为了画解调波形需要%否则xlabel(时间ms);title(调制信号);ylabel(幅值：v);subplot(5,1,2);plot(t,carwav(1:length(t);axis(0

33、,N,-2,2);xlabel(时间ms);title(载波恢复);ylabel(幅值：v);subplot(5,1,3);plot(t,c_re(1:length(t);axis(0,N,-2,2);xlabel(时间ms);title(载波相乘);ylabel(幅值：v);subplot(5,1,4);plot(t,c_ref(1:length(t);axis(0,N,-2.5,2.5);xlabel(时间ms);title(低通滤波后信号);ylabel(幅值：v);subplot(5,1,5);plot(t,c_de(1:length(t);axis(0,N,-0.2,1.5);xlabel(时间ms);title(相干解调信号);ylabel(幅值：v);