qpsk_tj.txt

分析QPSK调制解调系统理论和存在噪声情况下的误码率的程序

QPSK调制的误码率信噪比matlab小程序关键词：matlab  EbN0  误码率  汉字输出
算一算4816PSK的在10e-3  10e-4  10e-5下的EbN0

%%%%---QPSK-----
% snr=0:0.01:10;
% x=10.^(snr/10);
% Q=0.5*erfc(sqrt(x));   %(sin(pi/M))^2
% Ps=2*Q.*(1-0.25*Q);                 %/log2(M)
% Pe=Ps/2;
%semilogy(snr,Pe),grid
  

M=input('please input M=4,8,16:M=');
if(M==4|M==8)
    snr=0:0.01:M+6;
else if(M==16)
        snr=0:0.01:M+2;
    end
end
%------------------------------------------------------------
x=10.^(snr/10);
er=erfc(sqrt(log2(M)*x*(sin(pi/M))^2));   %(sin(pi/M))^2
Ps=er;               %/log2(M)
Pe=Ps/log2(M);
semilogy(snr,Pe),grid
%------------------------------------------------------------
Pe1=10^(-3);Pe2=10^(-4);Pe3=10^(-5);
x_1=find(Pe>Pe1);x1=max(x_1)*0.01;
xx1=sprintf('in case M=%1.0f,the snr of Pe=10e-3 is snr=%3.2f',M,x1);disp(xx1);
x_2=find(Pe>Pe2);x2=max(x_2)*0.01;
xx2=sprintf('in case M=%1.0f,the snr of Pe=10e-4 is snr=%3.2f',M,x2);disp(xx2);
x_3=find(Pe>Pe3);x3=max(x_3)*0.01;
xx3=sprintf('in case M=%1.0f,the snr of Pe=10e-5 is snr=%3.2f',M,x3);disp(xx3);