psk16.m

仿真了QPSK

clear;
close all;
Fd=1;
Fs=3*Fd;
M=4;
SNR_db=0:12;
SNR1_db=0:12;
pm=zeros(13,1);
pm_theory=zeros(13,1);
for n=1:length(SNR_db)
    Eb_N0=10^(SNR_db(n)/10);
    sgma=sqrt(1/(8*Eb_N0));
    x=randint(10000,1,M);
    y=dmodce(x,Fd,Fs,'psk',M);
    ynoise=y+sqrt(Fs/Fd)*sgma*(randn(length(y),1)+j*randn(length(y),1));
    z=ddemodce(ynoise,Fd,Fs,'psk',M);
    [numbers,pm(n)]=symerr(x,z);
end;
for m=1:length(SNR1_db)
    Eb_N0=10^(SNR1_db(m)/10);
    pm_theory(m)=2*Qfunct(sqrt(2*log2(M)*Eb_N0)*sin(pi/M));
    pe_theory(m)=pm_theory(m)/log2(M)
end
semilogy(SNR_db,pm,'*',SNR1_db,pm_theory);
xlabel('信噪比(dB)');
ylabel('误符号率');
legend('仿真得到的误符号率','理论符号误码曲线');