dt.m

来自「MATLAB实现的直接传输通信的统计误码率仿真。」· M 代码 · 共 27 行

clear;
P0=0.5;
Dsd=1;
Qsd=sqrt(10^(Dsd/10));
N=100;
SNR_dB=0:35; 
SNR=10.^(SNR_dB/10);
sigma2=sqrt(1./(10.^(SNR./10)));
simutimes=10000;%simulation times
for k=1:length(sigma2)
    for m=1:simutimes
    X=rand(1,N)>P0;%resource
    S=1-2*X;%BPSKmapping
    noise=sigma2(k)*(randn(1,N)+j*randn(1,N));%noise
    H=Qsd*(randn+j*randn);;% Channel
    y=H.*S+noise;%receive signal
    Y=y./H;
    z=(Y<0);
    Pe(m)=sum(abs(X-z))/N;
end
PE(k)=sum(Pe)/simutimes;
end
semilogy(SNR_dB,PE,'-*');
grid on;
ylabel('BER','FontSize',12);
xlabel('SNR(dB)','FontSize',12);