%版權所有17021211285孟慶統
%本程序設定：H0:x=n;
%???????????H1:x=2+n;
%假設H1時候信源輸出是2；噪聲為均值為0方差為1的高斯噪聲；采樣次數為1;
%最大后驗概率準則是在c10-c00=c01-c11時的貝葉斯準則，所以本程序修改為c11=2;c00=1;c10=4;c01=5;N=1;A=2;

clc
clear?all
close?all?

c11=2;c00=1;c10=4;c01=5;N=1;A=2;%設定參數??
k=100;
p1=0.3;%先驗概率P（H1）設為0.3
p0=1-p1;
th=zeros（k1）;
pf_p1=zeros（k1）;
pm_p1=zeros（k1）;
pd_p1=zeros（k1）;
C=zeros（k1）;
gama=zeros（k1）;
delta2=1;
delta=sqrt（delta2）;
d=sqrt（N*（A^2）/delta2）;
l=-4:0.01:6;
p_l_h0=sqrt（N/（2*pi*delta2））*exp（-N*l.*l/（2*delta2））;???????????%p（l|H0）
p_l_h1=sqrt（N/（2*pi*delta2））*exp（-N*（l-A）.*（l-A）/（2*delta2））;???%p（l|H1）

%%?求平均代價C和各判決概率隨判決門限
for?i=1:100
????th（i1）=i/10;
????gama（i1）=delta2/（N*A）*log（th（i1））+A/2;
????p0=1-p1;
????pf_p1（i1）=qfunc（log（th（i1））/d+d/2）;
????pd_p1（i1）=qfunc（log（th（i1））/d-d/2）;
????pm_p1（i1）=1-pd_p1（i1）;
????p00（i1）=1-pf_p1（i1）;
????C（i1）=p0*c10+p1*c11+p1*（c01-c11）*pm_p1（i1）-p0*（c10-c00）*（1-pf_p1（i1））;
end

%%?最佳判決門限下的C和各判決概率
????thg=（p0*（c10-c