最大后验概率准则matlab仿真

%版权所有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