%仿真時(shí)假設(shè)似然比檢測門限η為1

%信號幅度A不同的仿真
clear?all;?clc;

?M=1e5;??N=10;
?A=0.5:0.5:3;??%不同信號幅度A
?L=length（A）;
for?k=1:L;
????X0=zeros（MN）;
????X1=A（k）*ones（MN）;
????n=randn（MN）;
????Y0=X0+n;??%疊加了均值為0，方差為1的加性高斯噪聲
????Y1=X1+n;??%疊加了均值為0，方差為1的加性高斯噪聲
????y0=sum（Y02）/N;
????y1=sum（Y12）/N;
????
????y0（find（y0>A（k）/2））=A（k）;
????y0（find（y0????y1（find（y1>A（k）/2））=A（k）;
????y1（find（y1????
????
????err_num0=length（find（abs（y0-X0（:1））>1e-2））;
????err_num1=length（find（abs（y1-X1（:1））>1e-2））;
????
????BER10（k）=err_num0/M;
????ber10_theory（k）=qfunc（sqrt（N*A（k）^2）/2）;
????BER01（k）=err_num1/M;
????ber01_theory（k）=1-qfunc（-sqrt（N*A（k）^2）/2）;
????
????Pe（k）=（BER10（k）+BER01（k））/2;
????Pe_theory（k）=（ber10_theory（k）+ber01_theory（k））/2;
end

figure（1）
semilogy（ABER10‘-bo‘Aber10_theory‘-r*‘）;?grid?on;
xlabel（‘信號幅度A?‘）;
ylabel（‘誤碼率‘）;
legend（‘仿真的H0為真時(shí)的誤碼率‘‘理論的H0為真時(shí)的誤碼率‘）;
title（‘二元信號檢測中H0為真時(shí)判決H1成立的概率?‘）;

figure（2）
semilogy（APe‘-bo‘APe_theory‘-r*‘）;?grid?on;
xlabel（‘信號幅度A‘）;
ylabel（‘誤碼率‘）;
legend（‘仿真的誤碼率‘‘理論的誤碼率‘）;
title（‘二元數(shù)字通信系統(tǒng)的信號檢測性能?‘）;

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????文件????????1264??2014-11-25?18:40??symbol2_A.m
?????文件????????1219??2014-11-25?18:15??symbol2_N.m
?????文件????????3308??2014-11-25?19:24??symbol4_N.m
?????文件?????????424??2014-11-25?19:19??symbol4_picture.m
?????文件????????3227??2014-11-25?21:06??symbol4_sigma.m