%function [erro_dppm,erro_ppm,erro_ook]=mainMonteCarlo(N)
N=200;
M=50;
j=1;
[erro_dppm,erro_ppm,erro_ook]=mainappsimple(N);
erro_dppm1=erro_dppm;
erro_ppm1=erro_ppm;
erro_ook1=erro_ook;
erro_dppm=zeros(1,81);
erro_ppm=zeros(1,81);
erro_ook=zeros(1,81);

for j=1:M
[erro_dppm,erro_ppm,erro_ook]=mainappsimple(N);
erro_dppm1=erro_dppm1+erro_dppm;
erro_ppm1=erro_ppm1+erro_ppm;
erro_ook1=erro_ook1+erro_ook;
end
erro_dppm=erro_dppm1/M;
erro_ppm=erro_ppm1/M;
erro_ook=erro_ook1/M;
% create a figure
i=1;           
for SNR=0:0.2:16
    x(i)=SNR;
i=i+1;
end
    
figure('Name','DPPM、PPM和OOK的誤碼率比較','NumberTitle','off')
xlabel('SNR  信噪比 (dB)')
ylabel('Error 誤碼率 ')
title('DPPM,PPM和OOK判決信噪比--誤碼率')
axis([0 16 0 0.6])

hold all

% DPPM

plot(x,erro_dppm,':*')

% PPM

plot(x,erro_ppm,':d')

% OOK

plot(x,erro_ook,':v')

%
disp('計算中結束。')
%
% end