AMI及HDB3码波形及功率谱密度实验.docx
《AMI及HDB3码波形及功率谱密度实验.docx》由会员分享，可在线阅读，更多相关《AMI及HDB3码波形及功率谱密度实验.docx（15页珍藏版）》请在冰点文库上搜索。
AMI及HDB3码波形及功率谱密度实验
数字基带信号的波形与功率谱密度实验
一、实验目的
1、掌握数字基带码型有关概念及设计原则；
2、了解单极性码、双极性码、归零码和不归零码的波形特点；
3、掌握AMI和HDB3码的编码规则；
4、掌握各种基带码功率谱特性。
二、实验预习要求
1、复习?
数字通信原理?
第七章7.1节和7.2节——数字基带信号的码型与功率谱、AMI与HDB3码波形与功率谱密度；
2、学习MATLAB软件的使用；
3、认真阅读本实验容，熟悉实验步骤。
三、实验原理
通信的根本任务是远距离传递信息，因而如何准确地传输数字信息是数字通信的一个重要组成局部。
在数字传输系统中，其传输对象通常是二进制数字信息，它可能来自计算机、电传打字机或其它数字设备的各种数字代码，也可能来自数字终端的脉冲编码信号。
设计数字传输系统的根本考虑是选择一组有限的离散的波形来表示数字信息。
这些离散波形可以是未经调制的不同电平信号，也可以是调制后的信号形式。
由于未经调制的电脉冲信号所占据的频率带宽通常从直流和低频开场，因此称为数字基带信号。
而*些有线信道中，特别是传输距离不太远的情况下，数字基带信号可以直接传送，我们称之为数字信号的基带传输。
数字基带信号是数字信息的电脉冲表示，不同形式的数字基带信号〔又称码型〕具有不同的频谱构造，合理地设计数字基带信号以使数字信息变换为适合于给定信道传输特性的频谱构造，是基带传输首先要考虑的问题。
通常又把数字信息的电脉冲表示过程称为码型变换，在有线信道中传输的数字基带信号又称为线路传输码型。
事实上，在数字设备部用导线连接起来的各器件之间就是用一些最简单的数字基带信号来传送定时和信息的。
这些最简单的数字基带信号的频谱中含有丰富的低频分量乃到直流分量。
由于传输距离很近，高频分量衰减也不大。
但是数字设备之间长距离有线传输时，高频分量衰减随着距离的增加而增大，同时信道中往往还存在隔直流电容或耦合变压器，因而传输频带的高频和低频局部均受限。
此时必须考虑码型选择问题。
归纳起来，在设计数字基带信型时应考虑到以下原则：
1）在选用的码型的频谱中应该没有直流分量，低频分量也应尽量少。
这是因为终端机输出电路或再生中继器都是经过变压器与电缆相连接的，而变压器是不能通过直流分量和低频分量的。
2）传输型的频谱中高频分量要尽量少。
这是因为电缆号线之间的串话在高频局部更为严重，当码型频谱中高频分量较大时，就限制了信码的传输距离或传输质量。
3）码型应便于再生定时电路从码流中恢复位定时。
假设信号连“0〞较长，则等效于一段时间没有收脉冲，恢复位定时就困难，所以应该使变换后的码型中连“0〞较少。
4）设备简单，码型变换容易实现。
5）选用的码型应使误码率较低。
双极性基带信号波形的误码率比单极性信号低。
1、编码规则
〔1〕NRZ码
NRZ码的全称是单极性不归零码，在这种二元码中用高电平和低电平〔这里为零电平〕分别表示二进制信息“1〞和“0〞，在整个码元期间电平保持不变。
例如：
〔2〕RZ码
RZ码的全称是单极性归零码，与NRZ码不同的是，发送“1〞时在整个码元期间高电平只持续一段时间，在码元的其余时间则返回到零电平。
例如：
〔3〕AMI码
我们用“0〞和“1〞代表传号和空号。
AMI码的编码规则是“0〞码不变，“f0/2处（f0为码速率）。
信息代码：
10011000111……
AMI码：
+100-1+1000-1+1-1……
由于AMI码的传号交替反转，故由于它决定的基带信号将出现正负脉冲交替，而0电位保持不变的规律。
这种基带信号无直流成分，且只有很小的低频成分，因而它特别适宜在不允许这些成分通过的信道中传输。
除了上述特点以外，AMI码还有编译码电路简单以及便于观察误码情况等优点，它是以种根本的线路码，在高密度信息流得数据传输中，得到广泛采用。
但是，AMI码有一个重要缺点，即当它用来获取定时信息时，由于它可能出现长的连0串，因而会造成提取定时信号的困难。
〔4〕HDB3码
HDB3码是对AMI码的一种改良码，它的全称是三阶高密度双极性码。
其编码规则如下：
用B脉冲来保证任意两个相连取代节的V脉冲间“1〞的个数为奇数。
当相邻V脉冲间“1〞码数为奇数时，则用“000V〞取代，为偶数个时就用“B00V〞取代。
在V脉冲后面的“1〞码和B码都依V脉冲的极性而正负交替改变。
为了讨论方便，我们不管“0〞码，而把相邻的信码“1〞和取代节中的B码用B1B2......Bn表示，Bn后面为V，选取“000V〞或“B00V〞来满足Bn的n为奇数。
当信码中的“1〞码依次出现的序列为VB1B2B3...BnVB1时，HDB3码为＋－＋－...――＋或为－＋－＋...＋＋―。
由此看出，V脉冲是可以辩认的，这是因为Bn和其后出现的V有一样的极性，破坏了相邻码交替变号原则，我们称V脉冲为破坏点，必要时加取代节BOOV，保证n永远为奇数，使相邻两个V码的极性作交替变化。
由此可见，在HDB3码中。
相邻两个V码之间或是其余的“1〞码之间都符合交替变号原则，而取代码在整修码流中不符合交替变号原则。
经过这样的变换，既消除了直流成分，又防止了长连“0〞时位定时不易恢复的情况，同时也提供了取代信息。
图5-2给出了HDB3码的频谱，此码符合前述的对频谱的要求。
例如：
代码：
100001000011000011
AMI码：
-10000+10000-1+10000-1+1
HDB3码：
-1000-V+1000+V-1+1-B00-V+1-1
HDB3码的特点是明显的，它除了保持AMI码的优点外，还增加了使连0串减少到至多3个的优点，而不管信息源的统计特性如何。
这对于定时信号的恢复是十分有利的。
HDB3码是CCITT推荐使用的码型之一。
〔5〕CMI码
CMI码是传号反转码的简称，其编码规则为：
“1〞码交替用“11〞和“00〞表示；“0〞码用“01〞表示。
例如：
代码：
1101001
CMI码：
11000111010100
这种码型有较多的电平跃变，因此含有丰富的定时信息。
该码已被CCITT推荐为PCM〔脉冲编码调制〕四次群的接口码型。
在光缆传输系统中有时也用作线路传输码型。
〔6〕BPH码
BPH码的全称是数字双相码〔DigitalBiphase〕，又称Manchester码，即曼彻斯特码。
它是对每个二进制码分别利用两个具有2个不同相位的二进制新码去取代的码，编码规则之一是：
0→10〔零相位的一个周期的方波〕
1→01〔π相位的一个周期的方波〕
例如：
代码：
1100101
双相码：
01011010011001
双相码的特点是只使用两个电平，这种码既能提供足够的定时分量，又无直流漂移，编码过程简单，但这种码的带宽要宽些。
2、功率谱计算
数字基带信号一般是随机信号，因此不能用求确定信号频谱函数的方法来分析它的频谱特性。
随机信号的频谱特性必须用功率谱密度来描述。
对任意一种给定的数字基带信号，计算功率谱密度不是一件容易的事，往往需要繁复的数学计算。
假设数字基带信号以*种标准波形g（t）在码元周期Ts传送出去，则数字基带信号可用随机序列
〔2-1〕
表示。
这里，an是基带信号在nTs≤t≤（n+1）Ts时间间隔幅度值，由编码规律及输入信码决定；Ts为码元周期；g（t）为标准脉冲波形。
由于一般情况下，数字基带信号不是广义平稳随机过程，因此不能直接引用确定信号时自相关函数和功率谱密度之间存在的傅氏变换关系，即维纳-幸钦关系式。
但假设这种周期性平稳随机过程是各态历经性的，则可以导出它的平均功率谱密度计算式，于是：
〔2-2〕
其中，G（f）是脉冲波形g（t）的傅氏变换，
，〔n为任何值〕；
除了上式所定义的连续谱外，还在频率为
处存在如下所示的离散线谱：
〔2-3〕
这里，
为狄拉克
函数。
AMI码与HDB3码的频谱示意图如图2-1与2-2。
图2-1AMI码的频谱示意图
图2-2HDB3码的频谱示意图
四、实验仪器
WindowsNT/2000/*P/Windows7/VISTA；
MATLABV6.0以上。
五、实验容
1、利用MATLAB软件，编写.M文件，随机产生一组单极归零〔RZ〕码与单极非归零〔NRZ〕码，并计算出功率谱，如图2-1。
图2-1单极性RZ和NRZ的波形及功率谱
2、利用MATLAB软件，编写.M文件，随机产生一组双极性归零码与双极性非归零码，并计算出功率谱,如图2-2。
图2-2双极性RZ和NRZ码的波形与功率谱
3、根据已产生的RZ码，通过AMI和HDB3编码规则，利用.M文件产生AMI码与HDB3码，如图2-3。
图2-3RZ、AMI与HDB3双极性归零码波形
4、根据已产生的NRZ码，通过AMI和HDB3编码规则，利用.M文件产生AMI码与HDB3码，如图2-4。
图2-4NRZ、AMI与HDB3双极性非归零码波形
5、分析上述两种情况下产生的AMI码与HDB3的功率谱，如图2-5和2-6。
图2-5RZ码、AMI和HDB3双极性归零码功率谱
图2-6NRZ码、AMI和HDB3双极性非归零码功率谱
六、参考程序代码
%%%%%参数初始化%%%%%%%%%%%%%
Ts=1;%码元周期
N_sample=8;%每个码元的抽样点数
dt=Ts/N_sample;%抽样时间间隔
N=30;%码元数
t=0:
dt:
（N*N_sample-1）*dt;
T=t（end）;
Again=100;%迭代次数
%%%数据初始化%%%%%%%%%%%%%%%
rz_code=zeros（1,length（t））;%单极性RZ码初始化
nrz_code=zeros（1,length（t））;%单极性NRZ码初始化
drz_code=zeros（1,length（t））;%双极性RZ码初始化
dnrz_code=zeros（1,length（t））;%双极性NRZ码初始化
ami_code=zeros（1,length（t））;%双极性AMI归零码初始化
hdb3_code=zeros（1,length（t））;%双极性HDB3归零码初始化
ami_code_n=zeros（1,length（t））;%双极性AMI不归零码初始化
hdb3_code_n=zeros（1,length（t））;%双极性HDB3不归零码初始化
datab=zeros（1,2*N）;%
datad=zeros（1,2*N）;
ami_data=zeros（1,2*N）;%AMI码数据初始化
hdb3_data=zeros（1,2*N）;%HDB3码数据初始化
Sum_RZ=zeros（1,N*N_sample）;%RZ码功率值初始化
Sum_dRZ=zeros（1,N*N_sample）;%双极性RZ码功率值初始化
Sum_AMI=zeros（1,N*N_sample）;%AMI归零码功率值初始化
Sum_HDB3=zeros（1,N*N_sample）;%HDB3归零码功率值初始化
Sum_NRZ=zeros（1,N*N_sample）;%NRZ码功率值初始化
Sum_dNRZ=zeros（1,N*N_sample）;%双极性NRZ码功率值初始化
Sum_AMI_N=zeros（1,N*N_sample）;%AMI不归零码功率值初始化
Sum_HDB3_N=zeros（1,N*N_sample）;%HDB3不归零码功率值初始化
%%%%主程序单元%%%%%%%%%%%%%%%%%%%%%
r=1:
2:
2*N;
forii=1:
Again
dataa=round（rand（1,N））;
datab（r）=dataa;
datac=sign（randn（1,N））;
datad（r）=datac;
ami_data_n=ami（dataa）;
hdb3_data_n=hdb3（dataa）;
ami_data（r）=ami_data_n;
hdb3_data（r）=hdb3_data_n;
forj=1:
N
nrz_code（（j-1）*N_sample+1:
j*N_sample）=dataa（j）;
dnrz_code（（j-1）*N_sample+1:
j*N_sample）=datac（j）;
ami_code_n（（j-1）*N_sample+1:
j*N_sample）=ami_data_n（j）;
hdb3_code_n（（j-1）*N_sample+1:
j*N_sample）=hdb3_data_n（j）;
end
forj=1:
2*N
rz_code（（j-1）*N_sample/2+1:
j*N_sample/2）=datab（j）;
drz_code（（j-1）*N_sample/2+1:
j*N_sample/2）=datad（j）;
ami_code（（j-1）*N_sample/2+1:
j*N_sample/2）=ami_data（j）;
hdb3_code（（j-1）*N_sample/2+1:
j*N_sample/2）=hdb3_data（j）;
end
[f,RZ]=T2F（t,rz_code（1:
length（t）））;
[f,dRZ]=T2F（t,drz_code（1:
length（t）））;
[f,AMI]=T2F（t,ami_code（1:
length（t）））;
[f,HDB3]=T2F（t,hdb3_code（1:
length（t）））;
[f,NRZ]=T2F（t,nrz_code（1:
length（t）））;
[f,dNRZ]=T2F（t,dnrz_code（1:
length（t）））;
[f,AMI_N]=T2F（t,ami_code_n（1:
length（t）））;
[f,HDB3_N]=T2F（t,hdb3_code_n（1:
length（t）））;
Sum_RZ=Sum_RZ+10*log10（abs（RZ）.^2/T）;
Sum_dRZ=Sum_dRZ+10*log10（abs（dRZ）.^2/T）;
Sum_AMI=Sum_AMI+10*log10（abs（AMI）.^2/T）;
Sum_HDB3=Sum_HDB3+10*log10（abs（HDB3）.^2/T）;
Sum_NRZ=Sum_NRZ+10*log10（abs（NRZ）.^2/T）;
Sum_dNRZ=Sum_dNRZ+10*log10（abs（dNRZ）.^2/T）;
Sum_AMI_N=Sum_AMI_N+10*log10（abs（AMI_N）.^2/T）;
Sum_HDB3_N=Sum_HDB3_N+10*log10（abs（HDB3_N）.^2/T）;
end
PRZ=Sum_RZ/Again;
PdRZ=Sum_dRZ/Again;
PAMI=Sum_AMI/Again;
PHDB3=Sum_HDB3/Again;
PNRZ=Sum_NRZ/Again;
PdNRZ=Sum_dNRZ/Again;
PAMI_N=Sum_AMI_N/Again;
PHDB3_N=Sum_HDB3_N/Again;
figure
（1）;
subplot（2,2,1）;
plot（t,rz_code）;
*label（'时间'）;
ylabel（'幅值'）;
title（'单极性RZ码波形'）;
a*is（[030-22]）;
subplot（2,2,2）;
plot（t,nrz_code）;
*label（'时间'）;
ylabel（'幅值'）;
title（'单极性NRZ码波形'）;
a*is（[030-22]）;
subplot（2,2,3）;
plot（f,PRZ）;
*label（'频率'）;
ylabel（'幅值（dB）'）;
title（'单极性RZ码功率谱'）;
a*is（[-55-4020]）;
subplot（2,2,4）;
plot（f,PNRZ）;
*label（'频率'）;
ylabel（'幅值（dB）'）;
title（'单极性NRZ码功率谱'）;
a*is（[-55-4020]）;
figure
（2）;
subplot（2,2,1）;
plot（t,drz_code）;
*label（'时间'）;
ylabel（'幅值'）;
title（'双极性RZ码波形'）;
a*is（[030-22]）;
subplot（2,2,2）;
plot（t,dnrz_code）;
*label（'时间'）;
ylabel（'幅值'）;
title（'双极性NRZ码波形'）;
a*is（[030-22]）;
subplot（2,2,3）;
plot（f,PdRZ）;
*label（'频率'）;
ylabel（'幅值（dB）'）;
title（'双极性RZ码功率谱'）;
a*is（[-55-4020]）;
subplot（2,2,4）;
plot（f,PdNRZ）;
*label（'频率'）;
ylabel（'幅值（dB）'）;
title（'双极性NRZ码功率谱'）;
a*is（[-55-4020]）;
figure（3）;
subplot（3,1,1）;
plot（t,rz_code）;
*label（'时间'）;
ylabel（'幅值'）;
title（'RZ码波形'）;
a*is（[030-22]）;
subplot（3,1,2）;
plot（t,ami_code）;
*label（'时间'）;
ylabel（'幅值'）;
title（'AMI双极性归零码波形'）;
a*is（[030-22]）;
subplot（3,1,3）;
plot（t,hdb3_code）;
*label（'时间'）;
ylabel（'幅值'）;
title（'HDB3双极性归零码波形'）;
a*is（[030-22]）;
figure（4）;
subplot（3,1,1）;
plot（t,nrz_code）;
*label（'时间'）;
ylabel（'幅值'）;
title（'NRZ码波形'）;
a*is（[030-22]）;
subplot（3,1,2）;
plot（t,ami_code_n）;
*label（'时间'）;
ylabel（'幅值'）;
title（'AMI双极性不归零码波形'）;
a*is（[030-22]）;
subplot（3,1,3）;
plot（t,hdb3_code_n）;
*label（'时间'）;
ylabel（'幅值'）;
title（'HDB3双极性不归零码波形'）;
a*is（[030-22]）;
figure（5）;
subplot（3,1,1）;
plot（f,PRZ）;
*label（'频率'）;
ylabel（'幅值（dB）'）;
title（'RZ码功率谱'）;
a*is（[-55-6020]）;
subplot（3,1,2）;
plot（f,PAMI）;
*label（'频率'）;
ylabel（'幅值（dB）'）;
title（'AMI双极性归零码功率谱'）;
a*is（[-55-6020]）;
subplot（3,1,3）;
plot（f,PHDB3）;
*label（'频率'）;
ylabel（'幅值（dB）'）;
title（'HDB3双极性归零码功率谱'）;
a*is（[-55-6020]）;
figure（6）;
subplot（3,1,1）;
plot（f,PNRZ）;
*label（'频率'）;
ylabel（'幅值（dB）'）;
title（'NRZ码功率谱'）;
a*is（[-55-6020]）;
subplot（3,1,2）;
plot（f,PAMI_N）;
*label（'频率'）;
ylabel（'幅值（dB）'）;
title（'AMI双极性不归零码功率谱'）;
a*is（[-55-6020]）;
subplot（3,1,3）;
plot（f,PHDB3_N）;
*label（'频率'）;
ylabel（'幅值（dB）'）;
title（'HDB3双极性不归零码功率谱'）;
a*is（[-55-6020]）;
七、实验报告要求
1、表达数字信号基带传输码的特点；
2、掌握各种码型的产生规则及波形特点；
3、分析AMI与HDB3码的编码规则与功率谱特点；
4、对改良实验容有什么建议.
升级会员 