2DPSK调制与解调系统的仿真
摘要设计了差分编码移相键控(2DPSK)调制解调系统的工作流程图,并利用Matlab 软件对该系统的动态进行了模拟仿真。利用仿真的结果,从基带信号的波形图可以衡量数字信号的传输质量;由系统的输入和输出波形图可以看出,仿真实验良好。2DPSK调制解调系统的仿真设计,为以后进一步研究基于Matlab的通信实验仿真系统奠定了坚实的基础。
关键词调制解调; 差分移相编码; 仿真设计
1、 2DPSK基本原理
1.1 2DPSK信号原理
2DPSK方式即是利用前后相邻码元的相对相位值去表示数字信息的一种方式。现假设用Φ表示本码元初相与前一码元初相之差，并规定：Φ＝0表示0码，Φ＝π表示1码。则数字信息序列与2DPSK信号的码元相位关系可举例表示如2PSK信号是用载波的不同相位直接去表示相应的数字信号而得出的，在接收端只能采用相干解调，它的时域波形图如图2.1所示。
图1.1 2DPSK信号
在这种绝对移相方式中，发送端是采用某一个相位作为基准，所以在系统接收端也必须采用相同的基准相位。如果基准相位发生变化，则在接收端回复的信号将与发送的数字信息完全相反。所以在实际过程中一般不采用绝对移相方式，而采用相对移相方式。
定义∆Φ为本码元初相与前一码元初相之差，假设：
∆Φ=0→数字信息“0”；
∆Φ=π→数字信息“1”。
则数字信息序列与2DPSK 信号的码元相位关系可举例表示如下： 数字信息： 1 0 1 1 0 1 1 1 0 1 DPSK 信号相位：0 π π 0 π π 0 π 0 0 π
或：π 0 0 π 0 0 π 0 π π 0
1.2 2DPSK 信号的调制原理
一般来说，2DPSK 信号有两种调试方法，即模拟调制法和键控法。2DPSK 信号的的模拟调制法框图如图1.2.1所示，其中码变换的过程为将输入的单极性不归零码转换为双极性不归零码。
图1.2.1 模拟调制法
2DPSK 信号的的键控调制法框图如图1.2.2所示，其中码变换的过程为将输入的基带信号差分，即变为它的相对码。选相开关作用为当输入为数字信息“0” 时接相位0，当输入数字信息为“1”时接pi 。
图1.2.2 键控法调制原理图
1.3 2DPSK 信号的解调原理
2DPSK 信号最常用的解调方法有两种，一种是极性比较和码变换法，另一种是差分相干解调法。 1.3.1 2DPSK 信号解调的极性比较法
它的原理是2DPSK 信号先经过带通滤波器，去除调制信号频带以外的在信道中混入的噪声，再与本地载波相乘，去掉调制信号中的载波成分，再经过低通滤波器去除高频成分，得到包含基带信号的低频信号，将其送入抽样判决器中进行抽样判决的到基带信号的差分码，再经过逆差分器，就
码变换
相乘
载波
s(t)
e o (t)
得到了基带信号。它的原理框图如图1.3.1所示。
图 1.3.1 极性比较解调原理图
1.3.2 2DPSK信号解调的差分相干解调法
差分相干解调的原理是2DPSK信号先经过带通滤波器，去除调制信号频带以外的在信道中混入的噪声，此后该信号分为两路，一路延时一个码元的时间后与另一路的信号相乘，再经过低通滤波器去除高频成分，得到包含基带信号的低频信号，将其送入抽样判决器中进行抽样判决，抽样判决器的输出即为原基带信号。它的原理框图如图1.3.2所示。
图 1.3.2 差分相干解调原理图
2、建立模型
2.1 差分和逆差分变换模型
差分变换模型的功能是将输入的基带信号变为它的差分码。逆码变换器原理图如下：
逆码变换器
（a）原理方框图
2.2 带通滤波器和低通滤波器的模型
带通滤波器模型的作用是只允许通过（fl，fh）范围内的频率分量、但将其他范围的频率分量衰减到极低水平。低通滤波器模型的作用是只允许通过（0，fh）范围内的频率分量，并且将其他范围
的频率分量衰减到极低水平。在Matlab中带通滤波器和低通滤波器的模型可以用编写程序来模拟。
2.3 抽样判决器模型
抽样判决器的功能是根据位同步信号和设置的判决电平来还原基带信号。在Matlab中抽样判决器可以用simulink中的模块来模拟。它的模型框图如图所示，它的内部结构图如图2.3所示。
图 3.3 抽样判决器
2.4 系统结构图
图2.4.1 系统结构图
2.4.2 2ＤＰＳＫ调制与解调总原理框图
图2.4.2 2ＤＰＳＫ调制与解调总原理框图
3、仿真
3.1 仿真程序
%- 2DPSK 调制与解调
%--------------------------------------------------- %>>>>>>>>>>>>>>>>>>Initial_Part>>>>>>>>>>>>>>>>>>>>> %---------------------------------------------------
function y=dpsk2()
fs = 30000;
Time_Hold_On = 0.1;
Num_Unit = fs * Time_Hold_On;
High_Level = ones ( 1, Num_Unit );
Low_Level = zeros ( 1, Num_Unit );
w = 300;
A = 1;
%--------------------------------------------------- %>>>>>>>>>>>>>>>>>>Initial_The_Signal>>>>>>>>>>>>>>> %---------------------------------------------------
Sign_Set = [0,1,1,0,1,0,0,1]
Lenth_Of_Sign = length ( Sign_Set );
st = zeros ( 1, Num_Unit * Lenth_Of_Sign );
sign_orign = zeros ( 1, Num_Unit * Lenth_Of_Sign );
sign_result = zeros ( 1, Num_Unit * Lenth_Of_Sign );
t = 0 : 1/fs : Time_Hold_On * Lenth_Of_Sign - 1/fs;
%--------------------------------------------------- %>>>>>>>>>>>Generate_The_Original_Signal>>>>>>>>>>>> %---------------------------------------------------
for I = 1 : Lenth_Of_Sign
