matlab设计DPSK信号调制与解调.docx
《matlab设计DPSK信号调制与解调.docx》由会员分享，可在线阅读，更多相关《matlab设计DPSK信号调制与解调.docx（7页珍藏版）》请在冰豆网上搜索。
matlab设计DPSK信号调制与解调
2DPSK调制与解调系统的仿真
1、2DPSK基本原理
2DPSK信号原理
2DPSK方式即是利用前后相邻码元的相对相位值去表示数字信息的一种方式。
现假设用Φ表示本码元初相与前一码元初相之差，并规定：
Φ＝0表示0码，Φ＝π表示1码。
则数字信息序列与2DPSK信号的码元相位关系可举例表示如2PSK信号是用载波的不同相位直接去表示相应的数字信号而得出的，在接收端只能采用相干解调，它的时域波形图如图所示。
图2DPSK信号
在这种绝对移相方式中，发送端是采用某一个相位作为基准，所以在系统接收端也必须采用相同的基准相位。
如果基准相位发生变化，则在接收端回复的信号将与发送的数字信息完全相反。
所以在实际过程中一般不采用绝对移相方式，而采用相对移相方式。
定义为本码元初相与前一码元初相之差，假设：
数字信息“0”；→数字信息“1”。
→则数字信息序列与2DPSK信号的码元相位关系可举例表示如下：
数字信息：
1011011101
DPSK信号相位：
或：
0?
2DPSK信号的调制原理
一般来说，2DPSK信号有两种调试方法，即模拟调制法和键控法。
2DPSK信号的的模拟调制法框图如图1.2.1所示，其中码变换的过程为将输入的单极性不归零码转换为双极性不归零码。
．
s（t）
（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）原理方框图
带通滤波器和低通滤波器的模型）范围内的频率分量、但将其fh带通滤波器模型的作用是只允许通过（fl，）fh，0（低通滤波器模型的作用是只允许通过他范围的频率分量衰减到极低水平。
．
范围内的频率分量，并且将其他范围的频率分量衰减到极低水平。
在Matlab中带通滤波器和低通滤波器的模型可以用编写程序来模拟。
抽样判决器模型
抽样判决器的功能是根据位同步信号和设置的判决电平来还原基带信号。
在Matlab中抽样判决器可以用simulink中的模块来模拟。
它的模型框图如图所示，它的内部结构图如图所示。
图抽样判决器
系统结构图
图2.4.1系统结构图
2ＤＰＳＫ调制与解调总原理框图
图2ＤＰＳＫ调制与解调总原理框图
3、仿真
仿真程序
%-2DPSK调制与解调
%---------------------------------------------------
%>>>>>>>>>>>>>>>>>>Initial_Part>>>>>>>>>>>>>>>>>>>>>
%---------------------------------------------------
functiony=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>>>>>>>>>>>>
%---------------------------------------------------
forI=1:
Lenth_Of_Sign
ifSign_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>>>>>>>>>>>>>>>>>>
%---------------------------------------------------
forI=1:
Lenth_Of_Sign
ifSign_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）;
）;
]2）/3*（A2）,/3*（A-1）,+Lenth_Of_Sign*（Time_Hold_On,0[axis（
title（'低通滤波后的波形'）;
grid
%---------------------------------------------------
%>>>>>>>>>>>>>抽样判决&逆码变换部分>>>>>>>>>>>>>>>
%---------------------------------------------------
forI=1:
Lenth_Of_Sign
ifdt（（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））;
升级会员 