2PSK和2DPSK系统仿真

1、一、 题目名称2PSK和2DPSK系统仿真二、题目意义 运用MATLAB编程和simulink框图实现2PSK和2DPSK调制解调过程，并且输出其调制及解调过程中的波形，讨论其调制和解调效果。三、设计原理1、2PSK系统（1）调制原理数字调相：如果两个频率相同的载波同时开始振荡，这两个频率同时达到正最大值，同时达到零值，同时达到负最大值，它们应处于"同相"状态；如果其中一个开始得迟了一点，就可能不相同了。如果一个达到正最大值时，另一个达到负最大值，则称为"反相"。一般把信号振荡一次（一周）作为360度。如果一个波比另一个波相差半个周期，我们说两个波的相位

2、差180度，也就是反相。当传输数字信号时，"1"码控制发0度相位，"0"码控制发180度相位。载波的初始相位就有了移动，也就带上了信息。相移键控是利用载波的相位变化来传递数字信息，而振幅和频率保持不变。在2PSK中，通常用初始相位0和分别表示二进制“1”和“0”。因此，2PSK信号的时域表达式为(t)=Acost+) 其中，表示第n个符号的绝对相位：= 因此，上式可以改写为 图2 2PSK信号波形 （2）解调原理2PSK信号的解调方法是相干解调法。由于PSK信号本身就是利用相位传递信息的，所以在接收端必须利用信号的相位信息来解调信号。下图2-3中给出了一

3、种2PSK信号相干接收设备的原理框图。图中经过带通滤波的信号在相乘器中与本地载波相乘，然后用低通滤波器滤除高频分量，在进行抽样判决。判决器是按极性来判决的。即正抽样值判为1，负抽样值判为0.2PSK信号相干解调各点时间波形如图 3 所示. 当恢复的相干载波产生180°倒相时,解调出的数字基带信号将与发送的数字基带信号正好是相反,解调器输出数字基带信号全部出错.图 32PSK信号相干解调各点时间波形 这种现象通常称为"倒"现象.由于在2PSK信号的载波恢复过程中存在着180°的相位模糊,所以2PSK信号的相干解调存在随机的"倒"现象,从

4、而使得2PSK方式在实际中很少采用. 2、2DPSK系统2DPSK方式即是利用前后相邻码元的相对相位值去表示数字信息的一种方式。现假设用表示本码元初相与前一码元初相之差，并规定：0表示0码，表示1码。则数字信息序列与2DPSK信号的码元相位关系可举例表示如2PSK信号是用载波的不同相位直接去表示相应的数字信号而得出的，在接收端只能采用相干解调，它的时域波形图如图2.1所示。图1.1 2DPSK信号在这种绝对移相方式中，发送端是采用某一个相位作为基准，所以在系统接收端也必须采用相同的基准相位。如果基准相位发生变化，则在接收端回复的信号将与发送的数字信息完全相反。所以在实际过程中一般不采用绝对移相

5、方式，而采用相对移相方式。定义 DF为本码元初相与前一码元初相之差，假设：DF=0数字信息“0”；DF=p数字信息“1”。则数字信息序列与2DPSK信号的码元相位关系可举例表示如下：数字信息： 1 0 1 1 0 1 1 1 0 1DPSK信号相位：0 p p 0 p p 0 p 0 0 p或：p 0 0 p 0 0 p 0 p p 0（1）调制原理一般来说，2DPSK信号有两种调试方法，即模拟调制法和键控法。2DPSK信号的的模拟调制法框图如图1.2.1所示，其中码变换的过程为将输入的单极性不归零码转换为双极性不归零码。码变换相乘载波s(t)eo(t)图1.2.1 模拟调制法2DPSK信号的

6、的键控调制法框图如图1.2.2所示，其中码变换的过程为将输入的基带信号差分，即变为它的相对码。选相开关作用为当输入为数字信息“0” 时接相位0，当输入数字信息为“1”时接pi。图1.2.2 键控法调制原理图（2）2DPSK信号的解调原理2DPSK信号最常用的解调方法有两种，一种是极性比较和码变换法，另一种是差分相干解调法。2DPSK信号解调的极性比较法：它的原理是2DPSK信号先经过带通滤波器，去除调制信号频带以外的在信道中混入的噪声，再与本地载波相乘，去掉调制信号中的载波成分，再经过低通滤波器去除高频成分，得到包含基带信号的低频信号，将其送入抽样判决器中进行抽样判决的到基带信号的差分码，再经

7、过逆差分器，就得到了基带信号。它的原理框图如图1.3.1所示。延迟T相乘器低通滤波器抽样判决器2DPSK带通滤波器图 1.3.1 极性比较解调原理图2DPSK信号解调的差分相干解调法：差分相干解调的原理是2DPSK信号先经过带通滤波器，去除调制信号频带以外的在信道中混入的噪声，此后该信号分为两路，一路延时一个码元的时间后与另一路的信号相乘，再经过低通滤波器去除高频成分，得到包含基带信号的低频信号，将其送入抽样判决器中进行抽样判决，抽样判决器的输出即为原基带信号。它的原理框图如图1.3.2所示。带通滤波器相乘器低通滤波器抽样判决器逆码变换本地载波2DPSK图 1.3.2 差分相干解调原理图四、源

8、程序1、2PSKclear all; close all;clc;max=10g=zeros(1,max);g=randint(1,max);%长度为max的随机二进制序列cp=;mod1=;f=2*2*pi;t=0:2*pi/199:2*pi;for n=1:length(g); if g(n)=0; A=zeros(1,200);%每个值200个点 else g(n)=1; A=ones(1,200); end cp=cp A; %s(t),码元宽度200 c=cos(f*t);%载波信号 mod1=mod1 c;%与s(t)等长的载波信号,变为矩阵形式endfigure(1);subpl

9、ot(4,2,1);plot(cp);grid on;axis(0 200*length(g) -2 2);title('二进制信号序列');cm=;mod=;for n=1:length(g); if g(n)=0; B=ones(1,200);%每个值200个点 c=cos(f*t); %载波信号 else g(n)=1; B=ones(1,200); c=cos(f*t+pi); %载波信号 end cm=cm B; %s(t),码元宽度200 mod=mod c; %与s(t)等长的载波信号endtiaoz=cm.*mod;%e(t)调制figure(1);subplo

10、t(4,2,2);plot(tiaoz);grid on;axis(0 200*length(g) -2 2);title('2PSK调制信号');figure(2);subplot(4,2,1);plot(abs(fft(cp);axis(0 200*length(g) 0 400);title('原始信号频谱');figure(2);subplot(4,2,2);plot(abs(fft(tiaoz);axis(0 200*length(g) 0 400);title('2PSK信号频谱');%带有高斯白噪声的信道tz=awgn(tiaoz,

11、10);%信号tiaoz中加入白噪声，信噪比为10figure(1);subplot(4,2,3);plot(tz);grid onaxis(0 200*length(g) -2 2);title('通过高斯白噪声信道后的信号');figure(2);subplot(4,2,3);plot(abs(fft(tz);axis(0 200*length(g) 0 400);title('加入白噪声的2PSK信号频谱');jiet=2*mod1.*tz;%同步解调figure(1);subplot(4,2,4);plot(jiet);grid onaxis(0 200

12、*length(g) -2 2);title('相乘后信号波形')figure(2);subplot(4,2,4);plot(abs(fft(jiet);axis(0 200*length(g) 0 400);title('相乘后信号频谱');%低通滤波器fp=500;fs=700;rp=3;rs=20;fn=11025;ws=fs/(fn/2); wp=fp/(fn/2);%计算归一化角频率n,wn=buttord(wp,ws,rp,rs);%计算阶数和截止频率b,a=butter(n,wn);%计算H(z）figure(4);freqz(b,a,1000,1

13、1025);subplot(2,1,1);axis(0 4000 -100 3 )title('LPF幅频相频图');jt=filter(b,a,jiet);figure(1);subplot(4,2,5);plot(jt);grid onaxis(0 200*length(g) -2 2);title('经低通滤波器后信号波形')figure(2);subplot(4,2,5);plot(abs(fft(jt);axis(0 200*length(g) 0 400);title('经低通滤波器后信号频谱');%抽样判决for m=1:200*l

14、ength(g); if jt(m)<0; jt(m)=1; else jt(m)>=0; jt(m)=0; endendfigure(1);subplot(4,2,6);plot(jt);grid onaxis(0 200*length(g) -2 2);title('经抽样判决后信号s(t)波形')figure(2);subplot(4,2,6);plot(abs(fft(jt);axis(0 200*length(g) 0 400);title('经抽样判决后信号频谱');2、2DPSK%- 2DPSK 调制与解调%-%>>>

15、>>>>>>>>>>>>>>>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

16、 );w = 300;A = 1;%-%>>>>>>>>>>>>>>>>>>Initial_The_Signal>>>>>>>>>>>>>>>%-Sign_Set = 0,1,1,0,1,0,0,1Lenth_Of_Sign = length ( Sign_Set );st = zeros ( 1, Num_Unit * Lenth_Of_Sign );sign_orign = zeros ( 1,

17、 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) = 1s

18、ign_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 :

19、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

20、);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%-%>>>>>>>>>>>>>>&g

21、t;>>>相乘>>>>>>>>>>>>>>>>>>>>>>>>>>>>>%-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); %映射为数字的d