奈曼皮尔逊准则matlab仿真

%本程序以课本87页例题为例
%本程序设定：H0:x=n;
%           H1:x=1+n;
%假设H1时候信源输出是1；噪声为均值为0方差为1的高斯噪声；
clc
clear all
close all 
x=linspace（-34100）;
p_x_H0=sqrt（1/（2*pi））*exp（-x.*x/2）;
p_x_H1=sqrt（1/（2*pi））*exp（-（x-1）.*（x-1）/2）;
%化简得到检验统计量为l（x）=x，检测门限gama为log（th）+1/2th是似然比检测门限；
p_H1_H0=0.1;
syms P_H1_H0 l gama
P_H1_H0=int（1/sqrt（2*pi）*exp（-l.*l/2）lgama1000）;
gama=vpa（solve（P_H1_H0-0.1））
P_H1_H1=int（1/sqrt（2*pi）*exp（-（l-1）.*（l-1）/2）lgama1000）

   ssp1=1/（sqrt（2*pi））*exp（-（gama-1）.*（gama-1）/2）; %最佳检测门限的P（l|H1）值
   s=0:0.001:ssp1;
   ttp0=1/（sqrt（2*pi））*exp（-gama.*gama/2）;
   t=0:0.004:ttp0;
   y=0:0.001:3;
   figure
   plot（xp_x_H0‘b-‘xp_x_H1‘r-‘ga