BFSK/BPSK-BER 蒙特卡洛仿真程序

this is a homework of a course. hope it can help somebody~

%File: c9_MCBFSK.m
clear all
snrdB_min=0; snrdB_max=12;   %SNR（in dB）limits
snrdB = snrdB_min:1:snrdB_max;
Nsymbols = input（‘Enter number of symbols > ‘）;
snr = 10.^（snrdB/10）;       % convert from dB
h = waitbar（0‘SNR Iteration‘）;
len_snr = length（snrdB）;
for j=1:len_snr              %invrement SNR
    waitbar（j/len_snr）
    sigma = sqrt（1/（2*snr（j）））;   %noise standart deviation
    error_count=0;
    for k=1:Nsymbols       %simulation loop begins
        d = round（rand（1））;     %data
        if d ==0
            x_d = 1;
            x_q = 0;
        else
            x_d = 0;
            x_q = 1;
        end
        n_d = sigma*randn（1）;     %noise
        n_q = sigma*randn（1）;
        y_d = x_d+n_d;    %receiver input
        y_q = x_q+n_q;
        if y_d>y_q
            d_est = 0;
        else
            d_est = 1;
        end
        if （d_est ~= d）
            error_count = error_count+1;  %error counter
        end
    end %simulation loop ends
    errors（j） = error_count;   %store error count for plot
end
close（h）
ber_sim = errors/Nsymbols;     %BER estimate
ber_theor = erfc（sqrt（snr/2））/2;
semilogy（snrdBber_theorsnrdBber_sim‘o‘）
axis（[snrdB_min snrdB_max 0.0001 1]）
xlabel（‘SNR in dB‘）
ylabel（‘BER‘）
legend（‘Theoretical‘‘Simulation‘）
%End of script file

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1161  2008-10-17 15:43  BPSK\BPSK.m

     文件       1178  2008-10-17 14:22  BPSK\BPSK.asv

     文件       1331  2008-10-18 13:25  BFSK\bfsk.m

     文件       7959  2008-10-17 15:47  BFSK\nsymbols100000.fig

     文件       7949  2008-10-17 15:49  BFSK\nsymbols10000.fig

     文件       7915  2008-10-17 15:49  BFSK\nsymbols1000.fig

     文件       7888  2008-10-17 15:50  BFSK\nsymbols100.fig

     文件      18266  2008-10-17 16:11  BFSK\nsymbols100.jpg

     文件      18199  2008-10-17 16:11  BFSK\nsymbols1000.jpg

     文件      17983  2008-10-17 16:11  BFSK\nsymbols10000.jpg

     文件      18004  2008-10-17 16:11  BFSK\nsymbols100000.jpg

     文件       1331  2008-10-18 13:25  BFSK\bfsk.txt

     目录          0  2008-10-15 16:40  BPSK

     目录          0  2008-10-17 15:14  BFSK

----------- ---------  ---------- -----  ----

               109164                    14