《数据分析》实验报告三.docx

《数据分析》实验报告三.docx《数据分析》实验报告三
《数据分析》实验报告三
《数据分析》实验报告三
课程名称
数据剖析方法
课程编号
实验地址
系统建模与仿真切验室 SL110
实验时间
校外指导教师
无
校内指导教师
实验名称
鉴别准则。假如各整体间差异很小，做鉴别剖析的意义不大。
三、实验步骤
输入数据，鉴别两类整体的协方差阵能否相等，用马氏距离判断鉴别污染类型，计算
回代误判率与交错误判率，贝叶斯鉴别污染分类。
四、实验过程原始记录 (数据、图表、计算等 )
1、输入矩阵，计算协方差矩阵能否相等
A=[
]
>> x=[
]
G1=A([1:4 7:8 15],:);
G2=A([5:6 9:4 15],:);
n1=size(G1,1);
n2=size(G2,1);
n=n1+n2;
k=2;
p=6;
《数据分析》实验报告三
《数据分析》实验报告三
《数据分析》实验报告三
f=p*(p+1)*(k-1)/ 2;
d=(2*p^2+3*p-1)*(1 /(n1-1)+1/(n2-1)-1/(n-k))/(6*(p+1)*(k-1));
p1=n1/n;p2=n2 /n;
m1=mean(G1);m2=mean(G2);
s1=cov(G1);s2=cov(G2);
s=((n1-1)*s1+(n2-1)*s2)/(n-k);
M=(n-k)*log(det(s))-((n1-1)*log(det(s1))+(n2-1)*log(det(s2)));
T=(1-d)*M
=
+
C=chi2inv,f) C =
if T<chi2inv,f)
disp(' 两组数据协方差相等 ');
else
disp(' 两组数据协方差不相等 ');
end
两组数据协方差相等
2、结构鉴别函数，得出鉴别结果
>> for i=1:4
w(1)=m1*inv(s)*x(i,:)'-1 /2*m1*inv(s)*m1'+log(p1);
w(2)=m2*inv(s)*x(i,:)'-1 /2*m1*inv(s)*m2'+log(p2);
for j=1:2
if w(j)==max(w)
disp([' 待判样品属于第 ',num2str(j),' 类污染 ']);
end
end
end
待判样品属于第 2 类污染
《数据分析》实验报告三
《数据分析》实验报告三
《数据分析》实验报告三
待判样品属于第 2 类污染
待判样品属于第 2 类污染
待判样品属于第 2 类污染
3、计算回代误判率
n11=0;n22=0;
for i=1:n1
w1(i,1)=m1*inv(s)*G1(i,:)'-1 /2*m1*inv(s)*m1