c9_mcbpsk.asv

motecalo仿真实现源代码

snrdB_min=-3;
snrdB_max=8;

snrdB=snrdB_min:snrdB_max;
Nsymbols=input('enter number of symbols:');
snr=10.^(snrdB/10);
h=waitbar(0,'SNR Iteration');
len_snr=length(snrdB);
for j=1:len_snr
    waitbar(j/len_snr)
    sigma=sqrt(1/(2*snr(j)));
    error_count=0;
    for k=1:Nsymbols
        d=round(rand(1));
        x_d=2*d-1;
        n_d=sigma*randn(1);
        y_d=x_d+n_d;
        if y_d>0
            d_est=1;
        else
            d_est=0;
        end
        if (d_est~=d)
            error_count=error_count+1;
        end
    end
    errors(j)=error_count;
    en