2DPSK信号调制解调的matlab设计和仿真

2DPSK调制与解调系统的仿真摘要设计了差分编码移相键控(2DPSK)调制解调系统的工作流程图,并利用Matlab软件对该系统的动态进行了模拟仿真。利用仿真的结果,从基带信号的波形图可以衡量数字信号的传输质量;由系统的输入和输出波形图可以看出,仿真实验良好。2DPSK调制解调系统的仿真设计,为以后进一步研究基于Matlab的通信实验仿真系统奠定了坚实的基础。关键词调制解调;差分移相编码;仿真设计1、2DPSK基本原理1.12DPSK信号原理2DPSK方式即是利用前后相邻码元的相对相位值去表示数字信息的一种方式。现假设用Φ表示本码元初相与前一码元初相之差，并规定：Φ＝0表示0码，Φ＝π表示1码。则数字信息序列与2DPSK信号的码元相位关系可举例表示如2PSK信号是用载波的不同相位直接去表示相应的数字信号而得出的，在接收端只能采用相干解调，它的时域波形图如图2.1所示。图1.12DPSK信号在这种绝对移相方式中，发送端是采用某一个相位作为基准，所以在系统接收端也必须采用相同的基准相位。如果基准相位发生变化，则在接收端回复的信号将与发送的数字信息完全相反。所以在实际过程中一般不采用绝对移相方式，而采用相对移相方式。定义为本码元初相与前一码元初相之差，假设：→数字信息“0”；→数字信息“1”。1则数字信息序列与2DPSK信号的码元相位关系可举例表示如下：数字信息：1011011101DPSK信号相位：0或：1.22DPSK信号的调制原理一般来说，2DPSK信号有两种调试方法，即模拟调制法和键控法。2DPSK信号的的模拟调制法框图如图1.2.1所示，其中码变换的过程为将输入的单极性不归零码转换为双极性不归零码。图1.2.1模拟调制法2DPSK信号的的键控调制法框图如图1.2.2所示，其中码变换的过程为将输入的基带信号差分，即变为它的相对码。选相开关作用为当输入为数字信息“0”时接相位0，当输入数字信息为“1”时接pi。图1.2.2键控法调制原理图1.32DPSK信号的解调原理2DPSK信号最常用的解调方法有两种，一种是极性比较和码变换法，另一种是差分相干解调法。1.3.12DPSK信号解调的极性比较法它的原理是2DPSK信号先经过带通滤波器，去除调制信号频带以外的在信道中混入的噪声，再与本地载波相乘，去掉调制信号中的载波成分，再经过低通滤波器去除高频成分，得到包含基带信号的低频信号，将其送入抽样判决器中进行抽样判决的到基带信号的差分码，再经过逆差分器，就码变换相乘载波s(t)eo(t)2得到了基带信号。它的原理框图如图1.3.1所示。图1.3.1极性比较解调原理图1.3.22DPSK信号解调的差分相干解调法差分相干解调的原理是2DPSK信号先经过带通滤波器，去除调制信号频带以外的在信道中混入的噪声，此后该信号分为两路，一路延时一个码元的时间后与另一路的信号相乘，再经过低通滤波器去除高频成分，得到包含基带信号的低频信号，将其送入抽样判决器中进行抽样判决，抽样判决器的输出即为原基带信号。它的原理框图如图1.3.2所示。图1.3.2差分相干解调原理图2、建立模型2.1差分和逆差分变换模型差分变换模型的功能是将输入的基带信号变为它的差分码。逆码变换器原理图如下：2.2带通滤波器和低通滤波器的模型带通滤波器模型的作用是只允许通过（fl，fh）范围内的频率分量、但将其他范围的频率分量衰减到极低水平。低通滤波器模型的作用是只允许通过（0，fh）范围内的频率分量，并且将其他范围cab微分整流脉冲展宽逆码变换器（a）原理方框图相乘器低通滤波器抽样判决器2DPSK带通滤波器带通滤波器相乘器低通滤波器抽样判决器逆码变换本地载波2DPSK延迟T3的频率分量衰减到极低水平。在Matlab中带通滤波器和低通滤波器的模型可以用编写程序来模拟。2.3抽样判决器模型抽样判决器的功能是根据位同步信号和设置的判决电平来还原基带信号。在Matlab中抽样判决器可以用simulink中的模块来模拟。它的模型框图如图所示，它的内部结构图如图2.3所示。图3.3抽样判决器2.4系统结构图图2.4.1系统结构图42.4.22ＤＰＳＫ调制与解调总原理框图图2.4.22ＤＰＳＫ调制与解调总原理框图3、仿真3.1仿真程序%-2DPSK调制与解调%---------------------------------------------------%Initial_Part%---------------------------------------------------functiony=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%---------------------------------------------------forI=1:Lenth_Of_Sign5ifSign_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%---------------------------------------------------forI=1:Lenth_Of_SignifSign_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('相乘后的波形');grid6%---------------------------------------------------%低通滤波部分%---------------------------------------------------[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%---------------------------------------------------%抽样判决&逆码变换部分%---------------------------------------------------forI=1:Lenth_Of_Signifdt((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('逆码变换后的波形');grid3.2.12DPSK模拟调制和差分相干解调法仿真图78图3.4.22DPSK模拟调制和差分相干解调法仿真图4、调试过程及结论4.1差分相干法2DPSK信号经相关模块调试后的波形图如下：9调制过后加入高斯白噪声，连接到带通滤波器，去除调制信号以外的在信道中混入的噪声，再连接到相乘器。此相乘器是一路延时一个码元时间后与另一路信号相乘。作用是去除调制信号中的载波成分。信号经过低通滤波器后，去除高频成分，得到包含基带信号的低频信号。经过抽样判决，便还原成原始信号。105、心得体会本次课程设计在刚开始的过程中无从下手，手忙脚乱，时间又紧，最终决定用软件仿真来实现2DPSK调制解调的设计。通过这次课程设计我们能够比较系统的了解理论知识，掌握了2DPSK调制解调的工作原理及2DPSK调制解调系统的工作过程，学会了使用仿真软件Matlab，并学会通过应用软件仿真来实现某些通信系统的设计，对以后的学习和工作都起到了一定的作用，加强了动手能力和学业技能。通过这次课程设计还让我们知道了，我们平时所学的知识如果不加以实践的话等于纸上谈兵。课程设计主要是我们理论知识的延伸，它的目的主要是要在设计中发现问题，并且自己要能找到解决问题的方案，形成一种独立的意识。我们