差分编码移相键控(2DPSK)系统设计课程设计
- 格式:doc
- 大小:961.46 KB
- 文档页数:12
1.课程设计目的(1)加深对差分编码移相键控(2DPSK)调制解调系统基本理论知识的理解。
(2)掌握对MATLAB软件的应用。
2.课程设计要求(1)掌握课程设计的相关知识、概念清晰。
(2)原理图设计合理、能够正确运行。
3.相关知识3.1差分编码移相键控(2DPSK)信号原理差分编码移相键控(2DPSK)方式即是利用前后相邻码元的相对相位值去表示数字信息的一种方式。
现假设用Φ表示本码元初相与前一码元初相之差,并规定:Φ=0表示0码,Φ=π表示1码。
则数字信息序列与差分编码移相键控(2DPSK)信号的码元相位关系可举例表示如差分编码移相键控(2DPSK)信号是用载波的不同相位直接去表示相应的数字信号而得出的,在接收端只能采用相干解调,它的时域波形图如图3.1.1所示。
图3.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 π π 03.2差分编码移相键控(2DPSK)信号的调制原理一般来说,差分编码移相键控(2DPSK)信号有两种调试方法,即模拟调制法和键控法。
差分编码移相键控(2DPSK)信号的的模拟调制法框图如图3.1.2所示,其中码变换的过程为将输入的单极性不归零码转换为双极性不归零码。
图3.2 模拟调制法差分编码移相键控(2DPSK)信号的的键控调制法框图如图3.1.3所示,其中码变换的过程为将输入的基带信号差分,即变为它的相对码。
选相开关作用为当输入为数字信息“0” 时接相位0,当输入数字信息为“1”时接pi 。
码变换相乘载波s(t)e o (t)图3.3 键控法调制原理图3.3差分编码移相键控(2DPSK)信号的解调原理差分编码移相键控(2DPSK)信号最常用的解调方法有两种,一种是极性比较和码变换法,另一种是差分相干解调法。
3.4差分编码移相键控(2DPSK)信号解调的极性比较法它的原理是差分编码移相键控(2DPSK)信号先经过带通滤波器,去除调制信号频带以外的在信道中混入的噪声,再与本地载波相乘,去掉调制信号中的载波成分,再经过低通滤波器去除高频成分,得到包含基带信号的低频信号,将其送入抽样判决器中进行抽样判决的到基带信号的差分码,再经过逆差分器,就得到了基带信号。
它的原理框图如图3.1.4所示。
图3.4 极性比较解调原理图3.5差分编码移相键控(2DPSK)信号解调的差分相干解调法差分相干解调的原理是差分编码移相键控(2DPSK)信号先经过带通滤波器,去除调制信号频带以外的在信道中混入的噪声,此后该信号分为两路,一路延时一个码元的时间后与另一路的信号相乘,再经过低通滤波器去除高频成分,得到包含基带信号的低频信号,将其送入抽样判决器中进行抽样判决,抽样判决器的相乘器低通滤波器抽样判决器2DPSK带通滤波器延迟T输出即为原基带信号。
它的原理框图如图3.1.5所示。
图3.5 差分相干解调原理图4 建立模型4.1差分和逆差分变换模型差分变换模型的功能是将输入的基带信号变为它的差分码。
逆码变换器原理图3.2.1所示:4.2 带通滤波器和低通滤波器的模型带通滤波器模型的作用是只允许通过(fl ,fh )范围内的频率分量、但将其他范围的频率分量衰减到极低水平。
低通滤波器模型的作用是只允许通过(0,fh )范围内的频率分量,并且将其他范围的频率分量衰减到极低水平。
在Matlab 中带通滤波器和低通滤波器的模型可以用编写程序来模拟。
4.3 抽样判决器模型抽样判决器的功能是根据位同步信号和设置的判决电平来还原基带信号。
在Matlab 中抽样判决器可以用simulink 中的模块来模拟。
它的模型框图如图所示,它的内部结构图如图3.2.2所示。
带通 滤波器相乘器低通 滤波器 抽样 判决器逆码变换本地载波2DPSKcab微分整流脉冲展宽逆码变换器 图4.1 原理方框图图4.2 抽样判决器4.4 系统结构图图4.3系统结构图4.5 差分编码移相键控(2DPSK)调制与解调总原理框图图4.4 2DPSK调制与解调总原理框图5仿真5.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 = [1,1,0 ,0,1,0,1,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_Signif Sign_Set(I) == 1sign_orign( (I-1)*Num_Unit + 1 : I*Num_Unit) = High_Level;elsesign_orign( (I-1)*Num_Unit + 1 : I*Num_Unit) = Low_Level;endend%---------------------------------------------------%>>>>>>>>>>>>>>>>>>Modulation_Part>>>>>>>>>>>>>>>>>>%---------------------------------------------------for I = 1 : Lenth_Of_Signif Sign_Set(I) == 1st( (I-1)*Num_Unit + 1 : I*Num_Unit) = A * cos ( 2 * pi * w * t( (I-1)*Num_Unit + 1 : I*Num_Unit ) + ( pi / 2 ) );elsest( (I-1)*Num_Unit + 1 : I*Num_Unit) = A * cos ( 2 * pi * w * t( (I-1)*Num_Unit + 1 : I*Num_Unit ) );endendfiguresubplot ( 2, 1, 1 )plot(t, sign_orign);axis( [ 0 , Time_Hold_On *( Lenth_Of_Sign + 1), - (A / 2), A + (A / 2) ] );title ( '原始信号' );gridsubplot ( 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 );figuresubplot(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_Signif dt((2*I-1)*Num_Unit/2) < 0.25sign_result( (I-1)*Num_Unit + 1 : I*Num_Unit) = High_Level;elsesign_result( (I-1)*Num_Unit + 1 : I*Num_Unit) = Low_Level;endendfigureplot ( t, sign_result );axis( [ 0 , Time_Hold_On *( Lenth_Of_Sign + 1), - 3*(A / 2), 3*(A / 2) ] ); title ( '逆码变换后的波形' );grid5.2差分编码移相键控(2DPSK)模拟调制和差分相干解调法仿真图图4.4.2 2DPSK模拟调制和差分相干解调法仿真图6. 调试过程及结论6.1差分相干法差分编码移相键控(2DPSK)信号经相关模块调试后的波形图如下:调制过后加入高斯白噪声,连接到带通滤波器,去除调制信号以外的在信道中混入的噪声,再连接到相乘器。