matlab设计DPSK信号调制与解调

1、 2DPSK基本原理 2DPSK信号原理

2DPSK方式即是利用前后相邻码元的相对相位值去表示数字信息的一种方式。现假设用Φ表示本码元初相与前一码元初相之差，并规定：Φ＝0表示0码，Φ＝π表示1码。则数字信息序列与2DPSK信号的码元相位关系可举例表示如2PSK信号是用载波的不同相位直接去表示相应的数字信号而得出的，在接收端只能采用相干解调，它的时域波形图如图所示。

图 2DPSK信号

在这种绝对移相方式中，发送端是采用某一个相位作为基准，所以在系统接收端也必须采用相同的基准相位。如果基准相位发生变化，则在接收端回复的信号将与发送的数字信息完全相反。所以在实际过程中一般不采用绝对移相方式，而采用相对移相方式。

定义 为本码元初相与前一码元初相之差，假设：

→数字信息“0”； →数字信息“1”。

则数字信息序列与2DPSK信号的码元相位关系可举例表示如下： 数字信息： 1 0 1 1 0 1 1 1 0 1 DPSK信号相位：

或：0 

2DPSK信号的调制原理

一般来说，2DPSK信号有两种调试方法，即模拟调制法和键控法。2DPSK信号的的模拟调制法框图如图1.2.1所示，其中码变换的过程为将输入的单极性不归零码转换为双极性不归零码。

s(t) 码变换 相乘 eo(t) 载波 图1.2.1 模拟调制法

2DPSK信号的的键控调制法框图如图1.2.2所示，其中码变换的过程为将输入的基带信号差分，即变为它的相对码。选相开关作用为当输入为数字信息“0” 时接相位0，当输入数字信息为“1”时接pi。

图1.2.2 键控法调制原理图

2DPSK信号的解调原理

2DPSK信号最常用的解调方法有两种，一种是极性比较和码变换法，另一种是差分相干解调法。

2DPSK信号解调的极性比较法

它的原理是2DPSK信号先经过带通滤波器，去除调制信号频带以外的在信道中混入的噪声，再与本地载波相乘，去掉调制信号中的载波成分，再经过低通滤波器去除高频成分，得到包含基带信号的低频信号，将其送入抽样判决器中进行抽样判决的到基带信号的差分码，再经过逆差分器，就得到了基带信号。它的原理框图如图1.3.1所示。

2DPSK 带通滤波器 延迟T 相乘器 低通滤波器 抽样判决器 图 1.3.1 极性比较解调原理图

2DPSK信号解调的差分相干解调法

差分相干解调的原理是2DPSK信号先经过带通滤波器，去除调制信号频带以外的在信道中混入的噪声，此后该信号分为两路，一路延时一个码元的时间后与另一路的信号相乘，再经过低通滤波器去除高频成分，得到包含基带信号的低频信号，将其送入抽样判决器中进行抽样判决，抽样判决器的输出即为原基带信号。它的原理框图如图1.3.2所示。

2DPSK 带通 相乘低通 滤波抽样 逆码变本地载器 图 1.3.2 差分相干解调原理图

2、建立模型

差分和逆差分变换模型

差分变换模型的功能是将输入的基带信号变为它的差分码。逆码变换器原理图如下：

a b c 微分整脉冲展逆码变换器 （a）原理方框图

带通滤波器和低通滤波器的模型

带通滤波器模型的作用是只允许通过（fl，fh）范围内的频率分量、但将其他范围的频率分量衰减到极低水平。低通滤波器模型的作用是只允许通过（0，fh）

范围内的频率分量，并且将其他范围的频率分量衰减到极低水平。在Matlab中带通滤波器和低通滤波器的模型可以用编写程序来模拟。 抽样判决器模型

抽样判决器的功能是根据位同步信号和设置的判决电平来还原基带信号。在Matlab中抽样判决器可以用simulink中的模块来模拟。它的模型框图如图所示，它的内部结构图如图所示。

图 抽样判决器

系统结构图

图 2.4.1 系统结构图 2ＤＰＳＫ调制与解调总原理框图

图 2ＤＰＳＫ调制与解调总原理框图

3、仿真 仿真程序

%- 2DPSK 调制与解调

%--------------------------------------------------- %>>>>>>>>>>>>>>>>>>Initial_Part>>>>>>>>>>>>>>>>>>>>> %--------------------------------------------------- function y=dpsk2() fs = 30000; Time_Hold_On = ;

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 if Sign_Set(I) == 1

sign_orign( (I-1)*Num_Unit + 1 : I*Num_Unit) = High_Level;

else

sign_orign( (I-1)*Num_Unit + 1 : I*Num_Unit) = Low_Level;

end end

%--------------------------------------------------- %>>>>>>>>>>>>>>>>>>Modulation_Part>>>>>>>>>>>>>>>>>> %--------------------------------------------------- for I = 1 : Lenth_Of_Sign if Sign_Set(I) == 1

st( (I-1)*Num_Unit + 1 : I*Num_Unit) = A * cos ( 2 * pi * w *

t( (I-1)*Num_Unit + 1 : I*Num_Unit ) + ( pi / 2 ) ); else

st( (I-1)*Num_Unit + 1 : I*Num_Unit) = A * cos ( 2 * pi * w *

t( (I-1)*Num_Unit + 1 : I*Num_Unit ) ); end end figure

subplot ( 2, 1, 1 ) plot(t, sign_orign);

axis( [ 0 , Time_Hold_On *( Lenth_Of_Sign + 1), - (A / 2), A + (A / 2) ] ); title ( '原始信号' ); grid

subplot ( 2, 1, 2 ); plot ( t, st );

axis( [ 0 , Time_Hold_On *( Lenth_Of_Sign + 1), - 3*(A / 2), 3*(A / 2) ] ); title ( '调制后的信号' );

grid

%--------------------------------------------------- %>>>>>>>>>>>>>>>>>>相乘>>>>>>>>>>>>>>>>>>>>>>>>>>>>> %--------------------------------------------------- dt = st .* cos ( 2 * pi * w * t ); figure

subplot(2,1,1) plot ( t, dt );

axis( [ 0 , Time_Hold_On *( Lenth_Of_Sign + 1), - 3*(A / 2), 3*(A / 2) ] ); title ( '相乘后的波形' ); grid

%--------------------------------------------------- %>>>>>>>>>>>>>>>>>>>低通滤波部分>>>>>>>>>>>>>>>>>>>> %---------------------------------------------------

[N,Wn] = buttord( 2*pi*50, 2*pi*150,3,25,'s'); %临界频率采用角频率表示

[b,a]=butter(N,Wn,'s');

[bz,az]=impinvar(b,a,fs); %映射为数字的 dt = filter(bz,az,dt); subplot(2,1,2) plot ( t, dt );

axis( [ 0 , Time_Hold_On *( Lenth_Of_Sign + 1), - 3*(A / 2), 3*(A / 2) ] );

title ( '低通滤波后的波形' ); grid

%--------------------------------------------------- %>>>>>>>>>>>>>抽样判决 & 逆码变换部分>>>>>>>>>>>>>>> %--------------------------------------------------- for I = 1 : Lenth_Of_Sign if dt((2*I-1)*Num_Unit/2) <

sign_result( (I-1)*Num_Unit + 1 : I*Num_Unit) = High_Level;

else

sign_result( (I-1)*Num_Unit + 1 : I*Num_Unit) = Low_Level;

end end figure

plot ( t, sign_result );

axis( [ 0 , Time_Hold_On *( Lenth_Of_Sign + 1), - 3*(A / 2), 3*(A / 2) ] ); title ( '逆码变换后的波形' ); grid

2DPSK模拟调制和差分相干解调法仿真图

图 2DPSK模拟调制和差分相干解调法仿真图

4、调试过程及结论 差分相干法

2DPSK信号经相关模块调试后的波形图如下：

调制过后加入高斯白噪声，连接到带通滤波器，去除调制信号以外的在信道中混入的噪声，再连接到相乘器。此相乘器是一路延时一个码元时间后与另一路信号相乘。作用是去除调制信号中的载波成分。

信号经过低通滤波器后，去除高频成分，得到包含基带信号的低频信号。经过抽样判决，便还原成原始信号。 附录：

如果调用Matlab中自带函数，也可完成：

symbols = [1:6];. p = [ ];

[dict,avglen] = huffmandict(symbols,p);. msg = randsrc(1,10000,[symbols; p]); code1 = huffmanenco(msg,dict); M =6;

code2 = dpskmod(code1,M);

ynoisy = awgn(code2,200,'measured'); decode2= dpskdemod(ynoisy,M);

[number1,ratio1]=symerr(code1(1:end),decode2(1:end)); msg1 = huffmandeco(decode2,dict);

[number,ratio]=biterr(msg(1:end),msg1(1:end));