mpsk.m

mpsk matlab

%MPSK系统仿真
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
        %产生信源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('仿真结果','理论计算结果');