%MPSK系統(tǒng)仿真
clear all;
close all;

M=4;%QPSK

EsN0dB=3:0.5:10;
EsN0=10.^(EsN0dB/10);
Es=1;
N0=10.^(-EsN0dB/10);
sigma=sqrt(N0/2);

error=zeros(1,length(EsN0dB));
s_data=zeros(1,length(EsN0dB));

for k=1:length(EsN0dB)
    error(k)=0;
    s_data(k)=0;
    while error(k)<1000
        %產(chǎn)生信源1,2,3,4均勻分布
        d=ceil(rand(1,10000)*M);
        %調制成QPSK信號(復基帶形式)
        s=sqrt(Es)*exp(j*2*pi/M*(d-1));
        %加入信道噪聲(復噪聲)
        r=s+sigma(k)*(randn(1,length(d))+j*randn(1,length(d)));
        %判決
        for m=1:M;%計算距離
        rd(m,:)=abs(r-sqrt(Es)*exp(j*2*pi/M*(m-1)));
        end
        for m=1:length(s) %判決距離最近的點
            dd(m)=find(rd(:,m)==min(rd(:,m)));
            if dd(m)~=d(m)
                error(k)=error(k)+1;
            end
        end
        s_data(k)=s_data(k)+1000;
    end
end
Pe=error./s_data;
%理論計算的誤碼率結果
Ps=erfc(sqrt(EsN0)*sin(pi/M));
semilogy(EsN0dB,Pe,'b*-');
semilogy(EsN0dB,Ps,'rd-');
xlabel('Es/N0(dB)');
ylabel('誤碼率');
legend('仿真結果','理論計算結果');