BPSK和DPSK Matlab仿真

1、1. BPSK调制1.系统原理图1 BPSK调制系统原理框图BPSK调制系统的原理框图如图1所示，其中脉冲成形的作用是抑制旁瓣，减少邻道干扰，通常选用升余弦滤波器；加性高斯白噪声模拟信道特性，这是一种简单的模拟；带通滤波器BPF可以滤除有效信号频带以外的噪声，提高信噪比；在实际通信系统中相干载波需要使用锁相环从接收到的已调信号中恢复，这一过程增加了系统的复杂度，同时恢复的载波可能与调制时的载波存在180度的相位偏差，即180度相位反转问题，这使得BPSK系统在实际中无法使用；低通滤波器LPF用于滤除高频分量，提高信噪比；抽样判决所需的同步时钟需要从接收到的信号中恢复，即码元同步，判决门限跟码元

2、的统计特性有关，但一般情况下都为0。2.参数要求码元速率2400波特，载波频率4800Hz，奈奎斯特频率19200Hz。3.仿真参数采样频率：为码元速率的16倍，即奈奎斯特频率为码元速率的8倍，在计算误码率的时，采样频率为码元速率的8倍，即奈奎斯特频率为码元速率的4倍。脉冲成形滤波器参数：脉冲成形滤波器选用FIR型升余弦滤波器，滚降系数为0.5，阶数为21，调用MATLAB中的rcosfir函数来设计滤波器。低通滤波器：低通滤波器选用FIR型低通滤波器，参数为：f2=0 0.125 0.4 1; w2= 1 0.95 0.1 0;b2 = fir2(30,f2,w2);调用fir2函数来生成。

3、 两种调制/解调方式的各种参数都相同，后面不再赘述。4.脉冲成形滤波器的冲击响应和频率特性 图2 脉冲成形滤波器的冲击响应和频率特性5.低通滤波器的冲击响应和频率特性 图3 低通滤波器的冲击响应和频率特性 6.仿真过程各点波形：说明：在画图时，进行了归一化处理。噪声的方差为1，随机序列的幅度为1，调制加噪后信噪比为-3dB，按照理想情况计算，解调后信噪比为0dB。 图2 随机序列（16个）图3 脉冲成形后的波形脉冲成形的作用是抑制旁瓣，减少邻道干扰，在实际通信过程中必不可少，但是为了与理论误码率进行对比，在后面计算误码率时，没有进行脉冲成形。图4 调制后波形图5 加噪后波形 噪声的方差即平均功

4、率为1，信噪比为0dB，从图中可以看出，噪声比较大。图6 低通滤波后波形图7 误码情况（0dB） 从图中可以看出在0dB的信噪比情况下有明显的误码出现。2.2DPSK调制1.系统原理图8 2DPSK系统原理框图 与BSPK调制相比，2DPSK在调制前要进行差分编码，通过前后相邻的两个码元的相对相位来表征信息，这样做的好处是，在接收端不需要提取相干载波，也就不存在180度相位反转问题，所以是一种实用的系统。在接收端，对接收到的已调信号延迟一个码元间隔Ts，然后与接收信号相乘，再进行低通滤波、抽样判决，这种非相干解调方式实现简单，但误码率要高于BPSK调制，这可以通过后面的误码率统计看过。2. 仿

5、真过程各点波形：说明：在画图时，进行了归一化处理。噪声的方差为1，随机序列的幅度为1，调制加噪后信噪比为-3dB，按照理想情况计算，解调后信噪比为0dB。图9 随机序列（15个）图10 脉冲成形后的波形图11 调制后的波形图12 加噪后的波形 图13 移位相乘后波形 非相干解调的一个重要过程就是对接收到的信号进行一个码元周期的移位，将移位后的信号与接收到的信号相乘。再进行低通滤波，然后抽样判决。图 14 低通滤波后的波形图15 误码情况（0dB）3.误码率在统计误码率时，考虑到程序运行时间，码元个数为：65536，在高信噪比下统计的误码率可信度较低。判决前的信噪比按下式计算：SNR=A2/2

6、，其中A为已调信号幅度，2 为噪声方差。BPSK的理论误码率为：0.5*erfc(r)，其中r=SNR=A2/2为判决时的信噪比。 2DPSK的理论误码率为：0.5*exp(-r)。BPSK在-10dB,-8dB,-6dB,-4dB,-2dB,0dB,2dB,4dB,6dB,8dB,10dB的信噪比下的误码率依次为：0.3346252441406250.2932128906250000.2441711425781250.1925964355468750.1409759521484380.08575439453125000.04336547851562500.01483154296875000.

7、003280639648437500.0004882812500000004.57763671875000e-05 2DPSK在-10dB,-8dB,-6dB,-4dB,-2dB,0dB,2dB,4dB,6dB,8dB,10dB的信噪比下的误码率依次为： 0.4639047837033650.4448004882887010.4107270923933780.3701380941481650.3099107347219040.2290836957351030.1444266422522320.06889448386358440.02119478141451130.002670328831921

8、870.000183108262760357根据这两组数据以及对应的理论数据，绘制的误码率曲线如下图所示： 图16 误码率从图中可以看出实际仿真数据接近理论数据，并且在相同的信噪比情况下，BPSK调制的误码率小于2DPSK调制。附录：源程序1.BPSK.m% Bpsk调制，相干解调% 参数：% 载波频率4800Hz，码元速率2400Hz% 采样率是载波频率的4倍，即11920Hz% %Nyquist_Freq=19200;%奈奎斯特速率CW_Freq=4800;%载波频率Sample_Ratio=0.125;%载波频率/采样频率m=16;%码元数目N=2*m/Sample_Ratio;%采样点

9、数n=1:1:N;N0=randn(1,N); %噪声,均值为0，方差为1h0=zeros(1,N);% 21阶整形滤波器b1=rcosfir(0.5,10,2/Sample_Ratio,1);% 30阶低通滤波器 h2f2=0 0.125 0.4 1; w2= 1 0.95 0.1 0;b2 = fir2(30,f2,w2);% 随机序列a=rand(1,m);for i=1:m if(a(1,i)>0.5) a(1,i)=1; else a(1,i)=-1; end;end;% 对随机序列进行采样spread_a=zeros(1,N);k=1;j=0;for i=1:N if(j=N

10、/m) j=0; k=k+1; end j=j+1; spread_a(1,i)=a(1,k);end%脉冲成形fft_spread_a=fft(spread_a,N);fft_b1=fft(b1,N);wave_shape_a=ifft(fft_spread_a.*abs(fft_b1);%功率归1化energy_sum=0;for i=1:N energy_sum=energy_sum+wave_shape_a(1,i)2;endaverage_power=energy_sum/N;wave_shape_a=wave_shape_a/sqrt(average_power);%功率归1化%随

11、机序列%figure(1);subplot(2,1,1);plot(spread_a/max(spread_a),'k');% 随机序列时域波形title(' 随机序列时域波形');axis(1 N -1.5 1.5);subplot(2,1,2);plot(0:N-1)*2*pi/N,abs(fft_spread_a)/max(abs(fft_spread_a),'k');%随机序列幅频特性title(' 随机序列频域波形');%整形以后的信号%fft_wave_shape_a=fft(wave_shape_a,N);figur

12、e(2);subplot(2,1,1);plot(wave_shape_a/max(wave_shape_a),'k');% 整形以后的信号时域波形title('脉冲成形后时域波形');subplot(2,1,2);plot(0:N-1)*2*pi/N,abs(fft_wave_shape_a)/max(abs(fft_wave_shape_a),'k');%整形以后的信号幅频特性title('脉冲成形后频域波形');% 调制以后生成BPSK信号bpsk_m=zeros(1,N);j=1;k=1;for i=1:N if(j=(N

13、/m+1) j=1; k=k+1; end; bpsk_m(1,i)=wave_shape_a(1,i)*cos(2*pi*Sample_Ratio*i);%调制 j=j+1;end;%调制以后的信号%fft_bpsk_m=fft(bpsk_m,N);figure(3);subplot(2,1,1);plot(bpsk_m/max(bpsk_m),'k');% 调制以后的信号时域波形title('调制后时域波形');subplot(2,1,2);plot(0:N-1)*2*pi/N,abs(fft_bpsk_m)/max(abs(fft_bpsk_m),'

14、;k');%调制以后的信号幅频特性title('调制后频域波形');% 信号加噪声，模拟信道bpsk_m=bpsk_m+N0;%加噪以后的信号%fft_bpsk_m=fft(bpsk_m,N);figure(4);subplot(2,1,1);plot(bpsk_m/max(bpsk_m),'k');% 加噪以后的信号时域波形title('加噪后时域波形');subplot(2,1,2);plot(0:N-1)*2*pi/N,abs(fft_bpsk_m)/max(abs(fft_bpsk_m),'k');%加噪以后的信号

15、幅频特性title('加噪后频域波形');% 相干解调，用正交本振与信号相乘bpsk_dm=bpsk_m.*cos(2*pi*Sample_Ratio*n);%滤波fft_dm=fft(bpsk_dm,N);fft_b2=fft(b2,N);%N点的FFT变换fft_y1=fft_dm.*(abs(fft_b2);y1=ifft(fft_y1); %低通滤波后得到的信号%滤波以后的信号%figure(5);subplot(2,1,1);plot(y1/max(y1),'k');% 滤波以后的信号时域波形title('滤波后时域波形');subpl

16、ot(2,1,2);plot(0:N-1)*2*pi/N,abs(fft_y1)/max(abs(fft_y1),'k');%滤波以后的信号幅频特性title('滤波后频域波形');%采样和判决an=zeros(1,m); for i=1:m an(1,i)=y1(1,N/m*(i-1)+N/(2*m);end;for i=1:m if(an(1,i)>0) % 判决，得到解调结果 result(1,i)=1; else result(1,i)=-1; end;end;%误码情况%figure(6);plot(result,'ko');ho

17、ld on;plot(a,'kx');axis(1 m -1.5 1.5);legend('判决后','原码');title('误码情况');%3. DPSK.m% dpsk调制，非相干解调% 参数：% 载波频率4800Hz，码元速率2400Hz% 采样率是载波频率的8倍 %Nyquist_Freq=19200;%奈奎斯特速率CW_Freq=4800;%载波频率Sample_Ratio=0.125;%载波频率/采样频率m=15;%码元数目N=2*(m+1)/Sample_Ratio;%采样点数，因为差分编码后会多出一个码元，所以是

18、m+1n=1:1:N;N0=randn(1,N); %噪声,均值为0，方差为1% 21阶整形滤波器b1=rcosfir(0.5,10,2/Sample_Ratio,1);% 30阶低通滤波器 h2f2=0 0.125 0.4 1; w2= 1 0.95 0.1 0;b2 = fir2(30,f2,w2);% 随机序列a=rand(1,m);for i=1:m if(a(1,i)>0.5) a(1,i)=1; else a(1,i)=-1; end;end;%差分编码d_a=zeros(1,m+1);d_a(1,1)=-1;for i=1:m if(a(1,i)=d_a(1,i) d_a(

19、1,i+1)=-1; else d_a(1,i+1)=1; endend% 对原序列进行采样spread_a=zeros(1,N-2/Sample_Ratio);k=1;j=0;for i=1:N-2/Sample_Ratio if(j=2/Sample_Ratio) j=0; k=k+1; end j=j+1; spread_a(1,i)=a(1,k);end% 对差分编码后的序列进行采样spread_d_a=zeros(1,N);k=1;j=0;for i=1:N if(j=2/Sample_Ratio) j=0; k=k+1; end j=j+1; spread_d_a(1,i)=d_a

20、(1,k);end%脉冲整形fft_spread_d_a=fft(spread_d_a,N);fft_b1=fft(b1,N);wave_shape_d_a=ifft(fft_spread_d_a.*abs(fft_b1);%功率归1化energy_sum=0;for i=1:N energy_sum=energy_sum+wave_shape_d_a(1,i)2;endaverage_power=energy_sum/N;wave_shape_d_a=wave_shape_d_a/sqrt(average_power);%功率归1化%随机序列和差分编码后的波形%figure(1);subpl

21、ot(2,1,1);plot(spread_a/max(spread_a),'k');% 随机序列时域波形title(' 随机序列时域波形');axis(1 N -1.5 1.5);subplot(2,1,2);plot(spread_d_a/max(spread_d_a),'k');%随机序列幅频特性title(' 随机序列差分编码后波形');axis(1 N -1.5 1.5);%整形以后的信号%fft_wave_shape_d_a=fft(wave_shape_d_a,N);figure(2);subplot(2,1,1);

22、plot(wave_shape_d_a/max(wave_shape_d_a),'k');% 整形以后的信号时域波形title('脉冲成形后时域波形');subplot(2,1,2);plot(0:N-1)*2*pi/N,abs(fft_wave_shape_d_a)/max(abs(fft_wave_shape_d_a),'k');%整形以后的信号幅频特性title('脉冲成形后频域波形');% 调制以后生成dpsk信号dpsk_m=zeros(1,N);j=0;k=1;for i=1:N if(j=2/Sample_Ratio

23、) j=1; k=k+1; end dpsk_m(1,i)=wave_shape_d_a(1,i)*cos(2*pi*Sample_Ratio*i);%调制 j=j+1;end%调制以后的信号%fft_dpsk_m=fft(dpsk_m,N);figure(3);subplot(2,1,1);plot(dpsk_m/max(dpsk_m),'k');% 调制以后的信号时域波形title('调制后时域波形');subplot(2,1,2);plot(0:N-1)*2*pi/N,abs(fft_dpsk_m)/max(abs(fft_dpsk_m),'k&#

24、39;);%调制以后的信号幅频特性title('调制后频域波形');% 信号加噪声，模拟信道dpsk_m=dpsk_m+N0;%加噪以后的信号%fft_dpsk_m=fft(dpsk_m,N);figure(4);subplot(2,1,1);plot(dpsk_m/max(dpsk_m),'k');% 加噪以后的信号时域波形title('加噪后时域波形');subplot(2,1,2);plot(0:N-1)*2*pi/N,abs(fft_dpsk_m)/max(abs(fft_dpsk_m),'k');%加噪以后的信号幅频特性title('加噪后频域波形');% 非相干解调dpsk_m_shift=zeros(1,N);for i=1:N if(i<=2/Sample_Ratio) dpsk_m_shift(1,i)=0; else dpsk_m_shift(1,i)=dpsk_m(1,i-2/Sample_Ratio); end enddpsk_mul=dpsk_m_shift.*dpsk_m;fft_dpsk_mul=fft(dpsk_mul,N);%滤波fft_b2=fft(b2,N);%N点的FFT变换f