s3_9_q47.m

Duda《模式分类》第二版第1、3、5章部分课后习题和上机题的解答和程序代码

% 画出数据点的位置
clear all;close all;clc;
Dg=[1,1;3,3];
plot(Dg(:,1),Dg(:,2),'o','MarkerFaceColor','r');hold on;
Dby=-1:0.1:4;Dbx=2*ones(1,length(Dby));
plot(Dbx,Dby,'r');
% 分别画出x1，x2的概率分布密度函数
sita0=[2,4];
x1=0:0.01:5;
x2=0:0.01:sita0(2);
l1=length(x1);
l2=length(x2);
p1=exp(-x1/sita0(1))/sita0(1);
p2=(1/sita0(2))*ones(1,l2);
zl=zeros(1,l1-l2);
x2=x1;
p2=[p2,zl];
plot(x1,p1,'m','LineWidth',2);
plot(x2,p2,'b','LineWidth',2);
% xlabel('0<=θ<=1');ylabel('p(D|θ)')
sita1=[0.5,1.134];
x1=0:0.01:5;
x2=0:0.01:sita1(2);
l1=length(x1);
l2=length(x2);
p1=exp(-x1/sita1(1))/sita1(1);
p2=(1/sita1(2))*ones(1,l2);
zl=zeros(1,l1-l2);
x2=x1;
p2=[p2,zl];
plot(x1,p1,'-.m','LineWidth',2);
plot(x2,p2,'-.b','LineWidth',2);