s3_2_q1.m

来自「Duda《模式分类》第二版第1、3、5章部分课后习题和上机题的解答和程序代码」· M 代码 · 共 24 行

% 画出sita=1时，p(x|sita)关于x的函数图像
x1=0:0.1:10;
sita1=1;
px=sita1*exp(-sita1*x1);
plot(x1,px,'r','LineWidth',2);hold on;
% 画出x=2时，p(x|sita)关于sita的函数图像
sita2=0:0.1:5;
x2=2;
ps=sita2.*exp(-sita2.*x2);
plot(sita2,ps,'b','LineWidth',2);
% 生成n个服从p(x|1)分布的样本点
plist=[];
xlist=[];
n=100;
for i=1:n
    p=rand(1);
    plist(end+1)=p;
    x=-log(p/sita1)/sita1;
    xlist(end+1)=x;
end
% 画出这n个样本点与p(x|1)的关系，求出关于sita的最大似然估计值theta
plot(xlist,plist,'*');
theta=1/(mean(xlist));
legend({['p(x|sita) sita=',num2str(sita1)];['p(x|sita) x=',num2str(x2)];['样本点 n=',num2str(n),',sita~=',num2str(theta)]});