qpskwumalv.m

这是基于MATLAB的数字调相技术编程,包括对调制解调过程的编程以及其有效性和可靠性分析

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)+10000;
    end
%  drawnow
%  semilogy(EsN0dB,error./(s_data+eps));hold on;
end

Pe=error./s_data; 
Ps=erfc(sqrt(EsN0)*sin(pi/M)); %理论计算的误码率结果
semilogy(EsN0dB,Pe,'b*-');hold on;
semilogy(EsN0dB,Ps,'rd-');
xlabel('Es/N0(dB)');
ylabel('误码率');
legend('仿真结果','理论计算结果');