Matlab通信仿真系列——差分PSK(DPSK)仿真

微信公众号上线，搜索公众号小灰灰的FPGA,关注可获取相关源码，定期更新有关FPGA的项目以及开源项目源码，包括但不限于各类检测芯片驱动、低速接口驱动、高速接口驱动、数据信号处理、图像处理以及AXI总线等

本节目录

一、差分PSK信号的产生和调制
二、8-DPSK信号在awgn信道的性能仿真
三、Matlab源码

本节内容
一、差分PSK信号的产生和调制
在实际中，载波相位是从接收信号通过某些非线性运算提取的，会引入相位模糊。
载波相位估计中相位模糊，可以通过相连信号传输间的相位差来进行信息编码。
在二进制PSK中，信息比特1通过载波相位对前一载波相位180°相移发送，信息比特0则是通过与前一载波相位相同发送。
在四相PSK中，相继区间之间的相对相移0°,90°,180°和270°，分别对应信息比特00,01,11和10。
M>4的多相情况，采用该种编码处理产生的PSK信号称为差分编码PSK信号。
差分编码的相位调制和解调过程，检测器的接收信号被解调成M个可能发送的相位中的一个。在检测器之后是一个相位比较器，比较相邻信号间隔上已调信号的相位，以便提取信息。
具体框图如下：

二相DPSK的差错概率为

理论上在信噪比Eb/N0较高的时候，二相DPSK比二相PSK损失小3dB。
二、8-DPSK信号在awgn信道的性能仿真

matlab代码示例：
仿真8-DPSK载波调制信号在awgn信道下的误码率和误比特性能，
并与理论值比较。

三、Matlab源码

clear all;
clc;
close all;
%采样点数
numsamp=10;
%发送符合数
num_symbol=100000;
%信噪比
EsN0=0:20;
%信噪比转化为线性值
snrl=10.^(EsN0/10);
%8-DPSK
M=8;
%Gray编码规则
graycode=[0 1 2 3 6 7 4 5];
%消息数据
msg=randi([0,7],1,num_symbol);
%Gray映射
msg1=graycode(msg+1);
%8-DPSK调制
msgmod=dpskmod(msg1,M);
%求每个符号的平均功率
spow=norm(msgmod).^2/num_symbol;
for indx=1:length(EsN0)
    %根据符号功率求噪声功率
    sigma=sqrt(spow/(2*snrl(indx)));
    %加入高斯白噪声
    rx=msgmod+sigma*(randn(1,length(msgmod))+j*randn(1,length(msgmod)));   
    %8-DPSK解调
    y1=dpskdemod(rx,M);
    %Gray逆映射
    decmsg=graycode(y1+1);
    %求误比特率和误符号率
    [err,ber(indx)]=biterr(msg(2:end),decmsg(2:end),log2(M));
    [err,ser(indx)]=symerr(msg(2:end),decmsg(2:end));
end
%理论误符号率
ser1=2*qfunc(sqrt(snrl)*sin(pi/M));
%理论误比特率
ber1=1/log2(M)*ser1;


semilogy(EsN0,ber,'-ko',EsN0,ser,'-k*',EsN0,ber1,'-r',EsN0,ser1,'-b');
title('8-PSK载波信号在awgn信道下的性能');
xlabel('Es/N0');
ylabel('误比特率和误符号率');
legend('误比特率','误符号率','理论误比特率','理论误符号率');