数据通信实验指导书.

1、数据通信实验指导书信息技术学院目录实验一数字通信的基带码型3实验二 2ASK数字调制与解调7实验三 2FSK数字调制与解调12实验四 2PSK数字调制与解调18实验一数字通信的基带码型一、实验目的：1.了解几种常见的数字基本信号2.掌握常用数字基带传输码型的编码规则3.掌握用MATLAB仿真技术实现目前用于数字通信的基带码型二、实验仪器：安装Matlab软件的PC机一台三、实验原理：1、单极性不归零码NRZ（Non Return Zero）脉冲宽度等于码元宽度T特点：（1）有直流，零频附近的低频分量一般信道难传输。（2）收端判决门限与信号功率有关，不方便。（3）要求传输线一端接地。（4）不能用

2、滤波法直接提取位定时信号。2、双极性非归零码（BNRZ） ，有正负电平特点：不能用滤波直接提取位定时信号。3、单极性归零码（RZ） 特点：（1）可用滤波法提取位同步信号 （2）NRZ的缺点都存在4、双极性归零码（BRZ）特点：（1）整流后可用滤波提取位同步信号 （2）NRZ的缺点都不存在四、实验步骤：1、单极性不归零码NRZ（Non Return Zero）(1) 建立文件，MATLAB实现程序如下： function y=snrz(x) %本函数实现将输入的一段二进制代码编为相应的单极性不归零码输出%输入x为二进制码，输出y为编好的码t0=200;t=0:1/t0:length(x); %给

3、出相应的时间序列for i=1:length(x) %计算码元的值 if x(i)=1 %如果输入信息为1 for j=1:t0 %该码元对应的点值取1 y(i-1)*t0+j)=1; end else for j=1:t0 %如果输入信息为0，码元对应的点值取0 y(i-1)*t0+j)=0; end end endy=y,x(i);plot(t,y); %采用title命令来实现标记出各码元对应的二元信息title('1 0 1 1 0 0 1 0'); grid on;axis(0,i, -0.1,1.1);(2) 在命令窗口中键入x的二进制代码和函数名，就可以得到所对应

4、的单极性不归零码输出，如输入以下指令，将出现图1-2所示结果。x=1 0 1 1 0 0 1 0;snrz(x)图1-2 单极性不归零码2、双极性非归零码（BNRZ）双极性非归零码的实现同单极性基本一样，只需将snrz.m中的判断得到0信息后的语句y(i-1)*t0+j)=0;中的0改为-1即可，将axis(0,i, -0.1,1.1);中的-0.1改为-1.1即可，即axis(0,i, -1.1,1.1);所以就不再给出MATLAB函数文件了，波形图如图1-3所示。 图1-3 双极性不归零码3、单极性归零码（RZ）function y=srz(x)%本函数实现将输入的一段二进制代码编为相应的

5、单极性归零码输出 %输入x为二进制码，输出y为编好的码t0=200;t=0:1/t0:length(x); %给出相应的时间序列for i=1:length(x) %计算码元的值 if x(i)=1 %如果输入信息为1 for j=1:t0/2 y(2*i-2)*t0/2+j)=1; %定义前半段时间值为1 y(2*i-1)*t0/2+j)=0; %定义后半段时间值为0 endelse for j=1:t0 %如果输入信息为0 y(i-1)*t0+j)=0; %定义所有时间值为0 end endendy=y,x(i);plot(t,y);title('1 0 1 1 0 0 1 0&#

6、39;);grid on;axis(0,i, -0.1,1.1);同理，在命令窗口中键入x的二进制代码和函数名，就可以得到所对应的单极性归零码输出，如输入以下指令，将出现图7-4所示结果。 x=1 0 1 1 0 0 1 0; srz(x)图1-4 单极性归零码4、双极性归零码（BRZ）双极性归零码的MATLAB实现同单极性也基本一样，只需将srz.m中的判断得到0信息后的语句for j=1:t0 y(i-1)*t0+j)=0;改为for j=1:t0/2 y(2*i-2)*t0/2+j)=-1; y(2*i-1)*t0/2+j)=0;即可，所以也就不再给出MATLAB函数文件了，其波形图如图

7、1-5所示。 图1-5 双极性归零码五、实验报告要求：用如下码型重复步骤，并做出相应的波形记录：a 单极性不归零码 b 单极性归零码c 双极性不归零码 d双极性归零码实验二 2ASK数字调制与解调一、实验目的：1用MATLAB仿真技术实现数字调制与解调、基带数字调制与解调。2掌握用键控法产生2ASK信号的方法3了解2ASK信号的频谱与数字基带信号频谱之间的关系。二、实验仪器：安装Matlab软件的PC机一台三、实验原理：本次设计采用的流程图如图2-1所示。图 2-1 2ASK调制解调框图1. ASK调制原理2ASK二进制振幅调制就是用二进制数字基带信号控制正弦载波的幅度，使载波振幅随着二进制数

8、字基带信号而变化，而其频率和初始相位保持不变。信息比特是通过载波的幅度来传递的。其信号表达式为：，S(t)为单极性数字基带信号。由于调制信号只有0或1两个电平，相乘的结果相当于将载频或者关断，或者接通，它的实际意义是当调制的数字信号“1”时，传输载波；当调制的数字信号为“0”时，不传输载波。2ASK信号的时间波形e2ASK(t)随二进制基带信号S(t)通断变化。所以又被称为通断键控信号。典型波形如图2-2所示。图 2-2 典型2ASK波形e2ASK(t)为已调信号，它的幅度受s(t)控制，也就是说它的幅度上携带有s(t)的信息。2ASK信号的产生方法通常有两种：模拟调制法（相乘器法）和键控法。

9、模拟调制法就是用基带信号与载波相乘，进而把基带信号调制到载波上进行传输。键控法由s(t)来控制电路的开关进而进行调制。两种方法的调制如图2-3和图2-4所示。图 2-3 模拟调制法（相乘器法）图 2-4 键控法2. ASK解调原理2ASK有两种基本解调方法：相干解调法（同步检测法）和非相干解调法（包络检波法）。相干解调需要将载频位置的已调信号频谱重新搬回原始基带位置，因此用相乘器与载波相乘来实现。相乘后的信号只要滤除高频部分就可以了。为确保无失真还原信号，必须在接收端提供一个与调制载波严格同步的本地载波，这是整个解调过程能否顺利完好进行的关键。本次设计采用相干解调法。两种解调原理图如图2-5和

10、图2-6所示。图 2-5 相干解调法（同步检测法）图 2-6 非相干解调法（包络检波法）四、实验步骤：通过编写M文件程序，产生随机信号，按图2-1顺序对每一模块编程后。程序中注有需注意语句及解释。运行程序，实现2ASK的调制与解调过程。本次设计采用模拟调制法（相乘器法）和相干解调法。2ASK调制解调程序及注释clear allclose alli=10;%10个码元j=5000;t=linspace(0,5,j);%0-5之间产生5000个点行矢量，即分成5000份fc=10;%载波频率fm=i/5;%码元速率%产生基带信号x=(rand(1,i)%rand函数产生在0-1之间随机数，共1-1

11、0个figure(2)plot(x)a=round(x);%随机序列，round取最接近小数的整数figure(3)stem(a)%火柴梗状图st=t;for n=1:10 if a(n)<1; for m=j/i*(n-1)+1:j/i*n st(m)=0; end else for m=j/i*(n-1)+1:j/i*n st(m)=1; end endendfigure(1);subplot(421);plot(t,st);axis(0,5,-1,2);title('基带信号st');%载波s1=cos(2*pi*fc*t);subplot(422);plot(s1

12、);title('载波信号s1');%调制e_2ask=st.*s1;subplot(423);plot(t,e_2ask);title('已调信号');noise =rand(1,j);e_2ask=e_2ask+noise;%加入噪声subplot(424);plot(t,e_2ask);title('加入噪声的信号');%相干解调at=e_2ask.*cos(2*pi*fc*t);at=at-mean(at);%因为是单极性波形，还有直流分量，应去掉subplot(425);plot(t,at);title('与载波相乘后信号

13、9;);f,af = T2F(t,at);%通过低通滤波器t,at = lpf(f,af,2*fm);subplot(426);plot(t,at);title('相干解调后波形');%抽样判决for m=0:i-1; if at(1,m*500+250)+0.5<0.5; for j=m*500+1:(m+1)*500; at(1,j)=0; end else for j=m*500+1:(m+1)*500; at(1,j)=1; end endendsubplot(427);plot(t,at);axis(0,5,-1,2);title('抽样判决后波形

14、9;)五、实验报告要求：1根据实验仿真结果，画出相应波形图2.对仿真结果加以分析说明实验三 2FSK数字调制与解调一、实验目的：1.用MATLAB仿真技术实现数字调制与解调、基带数字调制与解调。2.掌握用键控法产生2FSK信号的方法；3.了解2FSK信号的频谱与数字基带信号频谱之间的关系。二、实验仪器：安装Matlab软件的PC机一台三、实验原理：本次设计所用流程图如图 2-9所示。图 2-9 2FSK调制解调框图1. FSK调制原理一个FSK信号可以看成是两个不同载波的2ASK信号的叠加。其解调和解调方法和FSK差不多。2FSK信号的频谱可以看成是f1和f2的两个2ASK频谱的组合。频移键控

15、是利用载波的频率来传递数字信号，在2FSK中，载波的频率随着二进制基带信号在f1和f2两个频率点间变化，频移键控是利用载波的频移变化来传递数字信息的。在2FSK中，载波的频率随基带信号在f1和f2两个频率点间变化。故其表达式为：典型波形如图 2-10所示。图 2-10 2FSK典型波形图2FSK的调制方式有两种，即模拟调频法和键控法。本次设计采用键控法。键控法中可以用二进制“1”来对应于载频f1，而“0”用来对应于另一频率f2，而这个可以用受矩形脉冲序列控制的开关电路对两个不同的独立的频率源f1、f2进行选择通。键控法原理图如图2-11示图 2-11 2FSK键控法原理图2. FSK解调原理

16、2FSK的解调方式有两种：相干解调方式和非相干解调方式。非相干解调是经过调制后的2FSK数字信号通过两个频率不同的带通滤波器f1、f2滤出不需要的信号，然后再将这两种经过滤波的信号分别通过包络检波器检波，最后将两种信号同时输入到抽样判决器同时外加抽样脉冲，最后解调出来的信号就是调制前的输入信号。其原理图如图2-12所示。图 2-12 2FSK非相干解调原理图相干解调是根据已调信号由两个载波f1、f2调制而成，则先用两个分别对f1、f2带通的滤波器对已调信号进行滤波，然后再分别将滤波后的信号与相应的载波f1、f2相乘进行相干解调，再分别低通滤波、用抽样信号进行抽样判决器即可。原理图如图3-14所

17、示。图 2-13 2FSK相干解调原理图四、实验步骤：通过编写M文件程序，产生随机信号，按流程图2.2.1顺序对每一模块编程后。程序中注有需注意语句及解释。运行程序，实现2FSK的调制与解调过程。本次设计中采用键控法调制法和相干解调法。2FSK调制解调程序及注释clear allclose alli=10;%基带信号码元数j=5000;a=round(rand(1,i);%产生随机序列t=linspace(0,5,j);f1=10;%载波1频率f2=5;%载波2频率fm=i/5;%基带信号频率%产生基带信号st1=t;for n=1:10 if a(n)<1; for m=j/i*(n-

18、1)+1:j/i*n st1(m)=0; end else for m=j/i*(n-1)+1:j/i*n st1(m)=1; end endendst2=t;%基带信号求反for n=1:j; if st1(n)>=1; st2(n)=0; else st2(n)=1; endend;figure(1);subplot(411);plot(t,st1);title('基带信号st1');axis(0,5,-1,2);subplot(412);plot(t,st2);title('基带信号反码st2');axis(0,5,-1,2);%载波信号s1=cos

19、(2*pi*f1*t)s2=cos(2*pi*f2*t)subplot(413),plot(s1);title('载波信号s1');subplot(414),plot(s2);title('载波信号s2');%调制F1=st1.*s1;%加入载波1F2=st2.*s2;%加入载波2figure(2);subplot(411);plot(t,F1);title('F1=s1*st1');subplot(412);plot(t,F2);title('F2=s2*st2');e_fsk=F1+F2;subplot(413);plot(t

20、,e_fsk);title('2FSK信号')%键控法产生的信号在相邻码元之间相位不一定连续nosie=rand(1,j);fsk=e_fsk+nosie;subplot(414);plot(t,fsk);title('加噪声后信号')%相干解调st1=fsk.*s1;%与载波1相乘f,sf1 = T2F(t,st1);%通过低通滤波器t,st1 = lpf(f,sf1,2*fm);figure(3);subplot(311);plot(t,st1);title('与s1相乘后波形');st2=fsk.*s2;%与载波2相乘f,sf2 = T2F

21、(t,st2);%通过低通滤波器t,st2 = lpf(f,sf2,2*fm);subplot(312);plot(t,st2);title('与s2相乘后波形');%抽样判决for m=0:i-1; if st1(1,m*500+250)<st2(1,m*500+250); for j=m*500+1:(m+1)*500; at(1,j)=0; end else for j=m*500+1:(m+1)*500; at(1,j)=1; end endend;subplot(313);plot(t,at);axis(0,5,-1,2);title('抽样判决后波形&

22、#39;)五、实验报告要求：1根据实验仿真结果，画出相应波形图2.对仿真结果加以分析说明实验四 2PSK数字调制与解调一、实验目的：1.用MATLAB仿真技术实现数字调制与解调、基带数字调制与解调。2.掌握用键控法产生2PSK信号的方法；3.了解2PSK信号的频谱与数字基带信号频谱之间的关系。二、实验仪器：安装Matlab软件的PC机一台三、实验原理：本次设计所用流程图如图2-17所示。图 2-17 2PSK调制解调框图1. PSK调制原理在二进制数字调制中，当正弦载波的相位随二进制数字基带信号离散变化时，则产生二进制移相键控(2PSK)信号。2PSK信号调制有两种方法，即模拟调制法和键控法。

23、通常用已调信号载波的 0°和 180°分别表示二进制数字基带信号的 1 和 0，模拟调制法用两个反相的载波信号进行调制。2PSK以载波的相位变化作为参考基准的，当基带信号为0时相位相对于初始相位为0°，当基带信号为1时相对于初始相位为180°。键控法，是用载波的相位来携带二进制信息的调制方式。通常用0°和180°来分别代表0和1。其时域表达式为：其中，2PSK的调制中an必须为双极性码。本次设计中采用模拟调制法。两种方法原理图分别如图2-18和图2-19所示。图 2-18 模拟调制法原理图图 2-19 键控法原理图2. PSK解调原理

24、由于2PSK的幅度是恒定的，必须进行相干解调。经过带通滤波的信号在相乘器中与本地载波相乘，然后用低通滤波器滤除高频分量，在进行抽样判决。判决器是按极性来判决的。即正抽样值判为1，负抽样值判为0。2PSK信号的相干解调原理图如图2-20所示，各点的波形如图2-21所示。由于2PSK信号的载波回复过程中存在着180°的相位模糊，即恢复的本地载波与所需相干载波可能相同，也可能相反，这种相位关系的不确定性将会造成解调出的数字基带信号与发送的基带信号正好相反，即“1”变成“0”吗“0”变成“1”，判决器输出数字信号全部出错。这种现象称为2PSK方式的“倒”现象或“反相工作”。但在本次仿真中是直

25、接给其同频同相的载波信号，所以不存在此问题。图 2-20 2PSK的相干解调原理图图 2-21 相干解调中各点波形图四、实验步骤：通过编写M文件程序（见附录），产生随机信号，按流程图2-17所示顺序对每一模块编程后。程序中注有需注意语句及解释。运行程序，实现2PSK的调制与解调过程。2PSK调制解调程序及注释clear allclose alli=10;j=5000;fc=4;%载波频率fm=i/5;%码元速率B=2*fm;t=linspace(0,5,j);%产生基带信号a=round(rand(1,i);%随机序列,基带信号figure(3);stem(a);st1=t;for n=1:1

26、0 if a(n)<1; for m=j/i*(n-1)+1:j/i*n st1(m)=0; end else for m=j/i*(n-1)+1:j/i*n st1(m)=1; end endendfigure(1);subplot(411);plot(t,st1);title('基带信号st1');axis(0,5,-1,2);%&&&&&&%基带信号求反%由于PSK中的是双极性信号，因此对上面所求单极性信号取反来与之一起构成双极性码st2=t; for k=1:j; if st1(k)>=1; st2(k)=0;

27、 else st2(k)=1; endend;subplot(412);plot(t,st2);title('基带信号反码st2');axis(0,5,-1,2);st3=st1-st2;subplot(413);plot(t,st3);title('双极性基带信号st3');axis(0,5,-2,2);%载波信号s1=sin(2*pi*fc*t);subplot(414);plot(s1);title('载波信号s1');%调制e_psk=st3.*s1;figure(2);subplot(511);plot(t,e_psk);title(&

28、#39;e_2psk');noise=rand(1,j);psk=e_psk+noise;%加入噪声subplot(512);plot(t,psk);title('加噪后波形');%相干解调psk=psk.*s1;%与载波相乘subplot(513);plot(t,psk);title('与载波s1相乘后波形');f,af = T2F(t,psk);%通过低通滤波器t,psk = lpf(f,af,B);subplot(514);plot(t,psk);title('低通滤波后波形');%抽样判决for m=0:i-1; if psk(1

29、,m*500+250)<0; for j=m*500+1:(m+1)*500; psk(1,j)=0; end else for j=m*500+1:(m+1)*500; psk(1,j)=1; end endendsubplot(515);plot(t,psk);axis(0,5,-1,2);title('抽样判决后波形')五、实验报告要求：1根据实验仿真结果，画出相应波形图2.对仿真结果加以分析说明附录实验二、三、四用到的傅立叶T2F函数%利用FFT计算信号的频谱并与信号的真实频谱的抽样比较。%脚本文件T2F.m定义了函数T2F，计算信号的傅立叶变换。function f,sf= T2F(t,st)%This is a function using the FFT function to calculate a signal's Fourier%T