高斯信道下的PSK调制

clear all
clc
N=1000;
M = 16;
tx = randint（1NM）;  %随机产生的发送信号tx
dpskSig = dpskmod（txM）;  %调制产生QPSK信号qpskSig
SNR = 0:4:20; 
BER = zeros（size（SNR））;
for i= 1:length（SNR）
    rxSig=awgn（dpskSigSNR（i））;
    rx = dpskdemod（rxSigM）; %解调后信号rx
    SER（i）=sum（tx~=rx）/N; %计算误码率BER
    BER（i） = SER（i）/log2（M）;
    capacity（i）=log2（1+SNR（i））;
end



scatterplot（dpskSig）;
xlabel（‘Re‘）;