%本程序以課本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