基于MATLAB的二二进制数字系统的调制(包括2ask-2fsk-2psk-2dpsk)

1、基于matlab实现三种二进制数字调制ask、fsk、psk课 程现代通信技术任课老师 学生姓名 学 号 专 业 电子信息科学与技术 联系方式 二零一七年十二月目录一、实验要求4二、实验仪器4三、实验原理4实验一 二进制幅度键控（2ask）5一、2ask信号的调制5（1）2ask定义5（2）2ask调制原理5（3）2ask信号的解调6二、matlab编程实现仿真7实验二 二进制频移键控（2fsk）10一、2fsk调制原理101、2fsk定义102、典型波形图103、2fsk调制114、2fsk解调原理112fsk解调原理11二、 matlab编程实现仿真12（1）生成序列13（2）时域波形13

2、（3）三者频谱图14实验三 二进制相移键控（2psk）15一、psk调制解调原理15（1）2psk定义15（2）2psk调制15（3）2psk解调16（4）2psk解调时域波形16二、matlab编程实现仿真16实验四 用simulink实现2ask/2fsk/2psk的仿真20一、2ask的仿真20二、2ask仿真实验结果及分析211、 2ask调制前和调制后的信号对比212、解调过程223、2ask调制前、解调后和零阶保持后的编码数据22三、2fsk的simulimk仿真23四、2fsk的仿真分析23（1）基带数字信号和载波23（2）调制与解调24五、2psk的simulink仿真25六、

3、2psk仿真结果分析25（1）基带信号及载波25（2）调制与解调26一、实验要求使用matlab软件仿真实现三种二进制数字调制ask、fsk、psk，满足下列要求：1.使用随机序列产生10位二进制基带信号an，显示基带信号的时域波形及频谱。2.载波频率设为基带信号速率的2倍，产生三种调制信号（fsk另一个频率设为4倍基带信号速率）并显示时域波形和频谱。3.分析比较基带信号和三种已调信号的带宽。二、实验仪器1. pc机一台（系统：windows 7 x64）；2. matlab2016a仿真软件。三、实验原理二进制数字调制系统的原理及实现数字通信系统, 按调制方式可以分为基带传输和带通传输。数字

4、基带信号的功率一般处于从零开始到某一频率（如06m）低频段，因而在很多实际的通信（如无线信道）中就不能直接进行传输，需要借助载波调制进行频谱搬移，将数字基带信号变换成适合信道传输的数字频带信号进行传输，这种传输方式，称为数字信号的频带传输或调制传输、载波传输。所谓调制，是用基带信号对载波波形的某参量进行控制，使该参量随基带信号的规律变化从而携带消息。对数字信号进行调制可以便于信号的传输；实现信道复用；改变信号占据的带宽；改善系统的性能。数字基带通信系统中四种基本的调制方式分别称为振幅键控（ask，amplitude-shift keying）、移频键控（ fsk，frequency-shift

5、 keying）、移相键控(psk，phase-shift keying )和差分移相键（dpsk，different phase-shift keying）。本次课程设计对前三种调制方式进行了仿真。实验一 二进制幅度键控（2ask）一、2ask信号的调制（1）2ask定义幅移键控是利用载波的幅度变化来传递数字信息，而其频率和初始相位保持不变。在2ask中，载波的幅度只有两种变化状态，分别对应二进制信息“0”或“1”。二进制振幅键控的表达式为：st=atcos0+, 0tt.式中，w0=2f0为载波的角频率；a(t)是随基带调制信号变化的时变振幅，即at=a1 ,当发送“1”时a2,当发送“0

6、”时（2）2ask调制原理典型波形如图所示：2ask信号的产生方法通常有两种：相乘法和键控法，相应的调制器如图：相乘法：由2ask信号的表达式可知2ask信号是调制信号和载波信号的乘积，因此与一般的模拟幅度调制方法类似，可用相乘法产生2ask信号，如图2（a）所示。键控法：由2ask信号的波形图可以看出2ask信号是用载波信号的有无来表示的，因此可用开关电路来控制载波的通断来产生2ask信号，如图2（b）所示，这里的开关电路受控制。（3）2ask信号的解调在接收端，2ask有两种基本的解调方法：非相干解调（包络检波法）相干解调（同步检测法）非相干解调方式（a）相干解调方式（b）相干解调法相干解

7、调也叫同步解调，就是利用相干波和接收到的2ask信号相乘分离出包含原始信号的低频信号，再进行抽样判决恢复数字序列。相干波必须是与发送端同频同相的正弦信号：z（t）=y(t)cos()=m(t)cos2()=m(t)1+cos()=m(t)+m(t)cos().式中1/2m(t)是基带信号，1/2m(t)cos(2)是频率为2的高频信号，利用低通滤波器可检测出基带信号，再经过抽样判决，即可恢复出原始数字信号序列an,2ask信号带宽为码元速率的2倍，即：b2ask=2rb.式中rb为信息速率。相干解调的原理图如下：非相干解调2ask信号与模拟调制中的am信号类似。所以,对2ask信号也能够采用非

8、相干解调(包络检波法)。在本次仿真中，由于2ask相干解调的误码率小于非相干解调，所以选取相干解调方式。二、matlab编程实现仿真采用matlab编程实现2ask调制，用随机函数生成10位二进制基带信号，同时生成3hz的载波，用相乘得出2ask信号，使用傅里叶变换分析基带信号、载波信号和2ask信号的频谱，得出结论。程序如下：n=10;x=;%调制信号生成a=(randint(1,10,2); for i=1:nif a(i)=1x(i*80-79:i*80)=ones(1,80);elsex(i*80-79:i*80)=zeros(1,80);endendfigure(1);subplot

9、(3,1,1)plot(x);title(调制信号);axis(0 800 -2 2);%载波信号t=0.0125:0.0125:n;carry=cos(2*pi*3*t);subplot(3,1,2)plot(carry);title(载波);axis(0 800 -2 2);%调制生成2ask信号st=x.*carry;subplot(3,1,3)plot(st);title(2ask信号);axis(0 800 -1.5 1.5);%频谱分析m=length(st);t=t(end);df=1/t;f=(-m/2:m/2-1)*df;xf=fftshift(abs(fft(x);figu

10、re(2);subplot(3,1,1);plot(f,xf);title(调制信号频谱);cf=abs(fft(carry);subplot(3,1,2);plot(f,cf);title(载波信号频谱);st1=fftshift(abs(fft(st);sf=fftshift(st1);subplot(3,1,3);plot(f,sf);title(2ask信号频谱);有如下实验结果：（1）生成序列为：0 1 1 1 1 0 0 1 0 0（2）调制信号、载波和已调信号时域波形如下：（3）三者频谱图如下：波形分析：由基带信号和已调信号对比可知，一个码元间隔tb内，有两个正弦波形，即载波频率

11、是码元速率的两倍；随机信号时域波形上升下降沿会稍有陡峭，说明产生基带信号an的方法还需再优化；由频谱分析可知，基带信号位于低频附近，而ask让基带频谱搬移到载波频率处。若基带信号的带宽为b，则ask信号的带宽应为；实验二 二进制频移键控（2fsk）一、2fsk调制原理1、2fsk定义载波的幅度不变，频率随an在两个频率中变化，用两个不同的频率携带传递二进制数字信息，当发送“1 ”时对应于某个载波频率c1发送 “0”时对应于另 个载频 c0 我们把这种调制方式称之为二进制频移键控（2fsk）。当二进制数字信息序列是具有逻辑电平1和0的单极性不归零码时，二进制频移键控信号可以看成是两个不同载频的二

12、进制幅度键控信号之和，所以，二进制频移键控信号的时域表达式为：一个fsk信号可以看成是两个不同载波的2ask信号的叠加。其解调和解调方法和ask差不多。2fsk信号的频谱可以看成是f1和f2的两个2ask频谱的组合。频移键控是利用载波的频率来传递数字信号，在2fsk中，载波的频率随着二进制基带信号在f1和f2两个频率点间变化，频移键控是利用载波的频移变化来传递数字信息的。在2fsk中，载波的频率随基带信号在f1和f2两个频率点间变化。故其表达式为：2、典型波形图图 2fsk典型波形图3、2fsk调制2fsk的调制方式有两种，即模拟调频法和键控法。本次设计采用键控法。键控法中可以用二进制“1”来

13、对应于载频f1，而“0”用来对应于另一频率f2，而这个可以用受矩形脉冲序列控制的开关电路对两个不同的独立的频率源f1、f2进行选择通。二者原理图如图所示：2fsk调制原理图4、2fsk解调原理2fsk信号的常用解调方法是相干解调（包络检波）和相干解调解调。原理是将2fsk信号分解为上下两路2ask信号分别解调，然后进行判决。这里的抽样判决是直接比较两路信号抽样值的大小，可以不专门设置门限判决规则。但应与调制规则相呼应，调制时若规定“1”符号对应载波频率f1，则接收时上支路的样值较大，应判为“1”；反之则判为“0”，即择大判决。图 2fsk解调原理2、 matlab编程实现仿真用randint函

14、数生成10个随机数，形成基带数字信号an。由于一个fsk信号可以看成是两个不同载波的2ask信号的叠加，故可利用：生成已调信号，再通过傅里叶变换得到频谱，分析数据。n=10;x=;%调制信号生成a=(randint(1,10,2); for i=1:nif a(i)=1x(i*80-79:i*80)=ones(1,80);elsex(i*80-79:i*80)=zeros(1,80);endendfigure(1)subplot(4,1,1)plot(x);title(调制信号);axis(0 800 -2 2);%载波信号t=0.0125:0.0125:n;ca1=sin(2*pi*2*t)

15、;ca2=sin(2*pi*4*t);subplot(4,1,2)plot(ca1);title(载波1);axis(0 800 -2 2);subplot(4,1,3)plot(ca2);title(载波2);axis(0 800 -2 2);%调制生成2ask信号x1=x;fsk=x.*ca1+x1.*ca2;subplot(4,1,4)plot(fsk);title(2fsk信号);axis(0 800 -1.5 1.5);%频谱分析m=length(fsk);t=t(end);df=1/t;f=(-m/2:m/2-1)*df;xf=fftshift(abs(fft(x);figure(

16、2);subplot(4,1,1);plot(f,xf);title(调制信号频谱);cf1=abs(fft(ca1);subplot(4,1,2);plot(f,cf1);title(载波1信号频谱);cf2=abs(fft(ca2);subplot(4,1,3);plot(f,cf2);title(载波2信号频谱);st1=fftshift(abs(fft(fsk);sf=fftshift(st1);subplot(4,1,4);plot(f,sf);title(2fsk信号频谱);有如下实验结果：（1）生成序列为： 0 1 0 1 0 1 1 0 1 0（2）调制信号、载波1信号、载波2

17、信号和2fsk信号时域波形如下：（3）三者频谱图如下：波形分析：由基带信号和已调信号对比可知：时，一个码元间隔tb内有两个正弦波形，即载波频率是码元速率的两倍，为；时，tb内有四个正弦波形，即载波频率是码元速率的四倍，为；可知，2fsk满足两载波频率为：由频谱分析可知，基带信号位于低频附近，而fsk让基带频谱搬移到载波频率处。若基带信号的带宽为b，则2fsk信号的带宽应为；实验三二进制相移键控（2psk）一、psk调制解调原理（1）2psk定义当am(t)=1，载波的幅度不变，瞬时相位偏移m(t)随数字信息序列m(t)而变化，则称为数字调相或相移键控，其时域表达式可一般地表示为：（2）2psk

18、调制在二进制数字调制中，当正弦载波的相位随二进制数字基带信号离散变化时，则产生二进制移相键控(2psk)信号。2psk信号调制有两种方法，即模拟调制法和键控法。通常用已调信号载波的 0和 180分别表示二进制数字基带信号的 1 和 0，模拟调制法用两个反相的载波信号进行调制。2psk以载波的相位变化作为参考基准的，当基带信号为0时相位相对于初始相位为0，当基带信号为1时相对于初始相位为180。键控法，是用未调载波的相位来携带二进制信息的调制方式。通常用0和180来分别代表0和1。其时域表达式为：其中，2psk的调制中an必须为双极性码。两种方法原理图分别如图2-18和图2-19所示。图 模拟调

19、制法原理图图 键控法原理图（3）2psk解调2psk信号的解调通常不能采用非相干解调方法，只能采用相干解调法，解调器原理框图如下：2psk信号相干解调各点波形如图所示。图中假设想干载波的基准相位与2psk信号的调制载波的基准相位一致。但是，由于在2psk信号的载波恢复过程中存在180的相位模糊，也就是当恢复的相干载波产生180倒相时，解调出的数字基带信号将与发送的数字基带信号正好是相反，解调器输出数字基带信号全部出错。这种现象通常称为“倒”现象。由于在2psk信号的载波恢复过程中存在着180的相位模糊，所以2psk信号的相干解调存在随机的“倒”现象，从而使得2psk方式在实际中很少采用。（4）

20、2psk解调时域波形图 2psk信号相干解调时域波形二、matlab编程实现仿真采用matlab编程实现2psk调制，用随机函数生成10位二进制基带信号，同时生成载波1和载波2（和载波1有的相差），用相乘之和得出2psk信号，使用傅里叶变换分析基带信号、载波信号和2ask信号的频谱，得出结论。n=10;x=;%调制信号生成a=(randint(1,10,2); for i=1:nif a(i)=1x(i*80-79:i*80)=ones(1,80);elsex(i*80-79:i*80)=zeros(1,80);endendfigure(1)subplot(4,1,1)plot(x);titl

21、e(调制信号);axis(0 800 -2 2);%载波信号t=0.0125:0.0125:n;ca1=sin(2*pi*2*t);ca2=sin(2*pi*2*t+pi);subplot(4,1,2)plot(ca1);title(载波1);axis(0 800 -2 2);subplot(4,1,3)plot(ca2);title(载波2);axis(0 800 -2 2);%调制生成2ask信号x1=x;psk=x.*ca1+x1.*ca2;subplot(4,1,4)plot(psk);title(2psk信号);axis(0 800 -1.5 1.5);%频谱分析m=length(p

22、sk);t=t(end);df=1/t;f=(-m/2:m/2-1)*df;xf=fftshift(abs(fft(x);figure(2);subplot(4,1,1);plot(f,xf);title(调制信号频谱);cf1=abs(fft(ca1);subplot(4,1,2);plot(f,cf1);title(载波信号0相位频谱);cf2=abs(fft(ca2);subplot(4,1,3);plot(f,cf2);title(载波信号1相位频谱);st1=fftshift(abs(fft(psk);sf=fftshift(st1);subplot(4,1,4);plot(f,sf);title(2psk信号频谱);有如下实验结果：（1）生成序列为：1 1 1 0 0 1 1 1 0 1（2）调制信号、载波1信号、载波2信号和2fsk信号时域波形如下：相位调制方式：，时，有；，时，有；（3）三者频谱如下：波形