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‘）;