psktest.m

来自「qpsk的simulink仿真 可看ber性能等参数」· M 代码 · 共 25 行

close all
clear all
no_symbols=10000; M=16;m=log2(M);N=2

EbN0_dB=0:14; EbN0=10.^(EbN0_dB./10); SER=zeros(1,length(EbN0_dB));

for kk=1:length(EbN0);
    s=randint(1,no_symbols,M);
    x=pskmod(s,M);
    n=(1/sqrt(2*log2(M)*EbN0(kk)))*(randn(1,no_symbols)+j*randn(1,no_symbols));
    y=x+n;
    z=pskdemod(y,M);
    SER(kk) = 1-sum(s==z,2)/size(s,2);
end

figure(1); semilogy(EbN0_dB,SER/m,'b*-'); hold

PB1= (N/m)*qfunc(sqrt(2*m*EbN0)*sin(pi/M));
semilogy(EbN0_dB,PB1,'ro-')

grid; title('BER Plot for PSK'); legend('Simulation', 'Theory'); 
xlabel('EbN0 - dB'); ylabel ('BER')