《数据分析》实验报告三

1、课程名称数据分析方法课程编号实验地点系统建模与仿真实验室SL110实验时间校外指导教师无校内指导教师实验名称实验3距离判别与贝叶斯判别分析评阅人签字成绩实验数据与内容我国山区某大型化工厂，在厂区及邻近地区挑选有代表性的15个大气取样点，每日4次冋时抽取大气样品，测定其中含有的6种气体的浓度，前后共4天，每个取样点每种气体实测16次，计算每个取样点每种气体的平均浓度，数据见表4-8。气体数据对应的污染地区分类见表4-8中最后一列。现有两个取自该地区的4个气体样本，气体指标见表4-8中后4行，试解决以下问题：1.判别两类总体的协方差矩阵是否相等,然后用马氏距离差别这 4个未知气体样本的污染类别，并

2、计算回代误判率与交叉误判率；若两类总体服从正态分布,第一类与第二类的先验概率分别为7/15、8/15 ,利用贝叶斯判别样本的污染分类。2先验概率为多少时，距离判别与贝时斯判别相同?调整先验概率对判别结果的影响是什么?3对第一类与第1类的先验概率分别为7/15、8/15，计算误判概率。裘札8大气样品數据舉吒祥硫代氧二啊化篠碳4环値氯丙烷开E堆I0. 0560. 0嗣0.0310.038QL 00B10.02212a 04Q0-0550+1000. UD0.02200.007 3130. 0500+0740. 040.0486 0071Q.020140. 0450* 0S00.1100. 1000

3、. 025 00.006 3150.0380. 1300.0790.1700. 058 00"432§0. 030Q.】Q0.070d 1600. 050 00.04 6270. 0340. 0950.0580. 1600. 2000.0沖1H化穗4JFWftS 烷环己烷0, 0300. 090仇1800. 2200.039190,0fl40.0.0290.3200.0126 04】2106 0850. 0760. 0190. 3000. 0100, 0402110.0640,0720.0200+2500. 0280. 038212CL 0540. 0450. 0220.

4、2800. 0210. 040211CLO4&0+0. 062a 2600. 0380. 036214a 0450.0920, 0720. 2000, 0350. 03221SD. 0690” 0B70. 027a 0500. 0&9Q. Q211样品】C>.0520-0840. 0210. 0370 007 10, 022持定样晶2<. 0410.0550. 1100.1100. 021 00, 007 3持定样駅3h 0300. n?0. 07201600 0560. 021待定样品4Q. 0740, 0830. 105o. m0. 020L 000待定一、实验

5、目的1 熟练掌握MATLA软件进行距离判别与贝叶斯判别的方法与步骤。2掌握判别分析的回代误判率与交叉误判率的编程。3掌握贝叶斯判别的误判率的计算。二、实验原理1) 在MATLAB，进行数据的判别分析命令为classify,其调用格式为：class= classify(sample,training,group' type ')将sample数据的每一行指定到训练集 training 的一个类中。Sample和training 必须具有相同的列数。group向量包含从1到组数的正整数，它指明训练营集中的每一 行属于哪一类。group和training 必须具有 相同的行 数。ty

6、pe '是可 选项，选'linear '表示总体为多元正态总体，选'quadratic '与'mahalanobis '。该函数 返回class，它是一个与sample具有相同行数的向量。Class的每一个元素指定sample 中对应元素的分类。通过计算sample和training 中每一行的马氏距离，classify 函数 决定sample中的每一行属于哪一个分类。2) 贝叶斯判别方法步骤第1步，验证两个总体服从二元正态分布；第 2步，检验两个总体的协方差矩阵相 等；估计两个总体的先验概率 p1、p2;利用MATLA软件计算。3) 回

7、代误判率设G1, G2为两个总体，x1,x2和y1,y2是分别来自G1, G2的训练样本，以全体 训练样本作为m+n个新样品，逐个代入已建立的判别准则中判别其归属，这个过程称为 回判。回判结果中若属于 G1的样品被误判为属于G2的个数为N1个，属于G2的样品被 误判为属于G1的个数为N2个，则误判估计为：PA= (N1+N2 /(m+n)误判率的回代估计易于计算。但是，卩八是由建立判别函数的数据反过来用作评估准 则的数据而得到的。所以有偏，往往比真实误判率小。当训练样本容量较大时，卩八可以作为真实误判率的一种估计。4) 交叉误判率估计是每次剔除一个样品，利用m+n-1个训练样本建立判别准则，再

8、 利用建立的准则对删除的样本进行判别。对每个样品做如上分析，以其误判的比例作为 误判率，步骤；从总体G1的训练样本开始，剔除其中一个样品，剩余的m-1个样品与G2中的全部样品建立判别函数；用建立的判别函数对剔除的样品进行判别；重复以上步骤，直到G1中的全部样本依次被删除又进行判别，其误判的样品个数记为 N1*；对G2的样品重复以上步骤，直到 G2中的全部样本依次被删除又进行判别，其误判 的样品个数记为N2*。于是交叉误判率估计为：pA*=(N1*+N2*)/(m+n)5) 贝叶斯判别的有效性可以通过平均误判率来确定。判别准则的误判率在一定程度上依赖于所考虑的各总体间的差异程度。各总体间差异越大

9、，就越有可能建立有效的判 别准则。如果各总体间差异很小，做判别分析的意义不大。三、实验步骤输入数据，判别两类总体的协方差阵是否相等，用马氏距离判断判别污染类别，计算 回代误判率与交叉误判率，贝叶斯判别污染分类。四、实验过程原始记录(数据、图表、计算等)1、输入矩阵，计算协方差矩阵是否相等>> A=0.05600.08400.03100.03800.00810.02200.0400 0.05500.1000 0.1100 0.0220 0.00730.06900.08700.02700.05000.08900.02100.05200.08400.02100.03700.00710.0

10、220>> x=0.052 0.084 0.021 0.037 0.0071 0.0220.04100.05500.11000.11000.02100.00730.03000.11200.07200.16000.05600.02100.07400.08300.10500.19000.02001.0000>> 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

11、=6;>> f=p*(p+1)*(k-1)/2;>> d=(2*pA2+3*p-1)*(1/( n1-1)+1/( n2-1)-1/( n-k)/(6*(p+1)*(k-1);>> p1= n1/n;p2=n2/n;>> m仁 mea n(G1);m2=mea n( G2);>> s仁 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(s

12、2);>> T=(1-d)*MT =-44.8237 + 0.9288i>> C=ch i2in v(0.95,f)C =32.6706>> if T<chi2i nv(0.95,f)disp('两组数据协方差相等');elsedisp('两组数据协方差不相等');end两组数据协方差相等2、构造判别函数，得出判别结果>> for i=1:4w(1)=m1*i nv(s)*x(i,:)'-1/2*m1* in v(s)*m1'+log(p1);w(2)=m2*i nv(s)*x(i,:)

13、9;-1/2*m1* in v(s)*m2'+log(p2);for j=1:2if w(j)=max(w)disp('待判样品属于第,num2str(j),'类污染');endend end 待判样品属于第2类污染待判样品属于第2类污染待判样品属于第2类污染待判样品属于第2类污染3、计算回代误判率>> n11=0;n 22=0;>> for i=1: n1w1(i,1)=m1*i nv(s)*G1(i,:)'-1/2*m1* in v(s)*m1'+log(p1); w1(i,2)=m2*i nv(s)*G1(i,:)&

14、#39;-1/2*m2* in v(s)*m2'+log(p2);for j=1:2if w1(i,j)=max(w1(i,:)&j=1n11=n11+1;endendend>> for i=1: n2w2(i,1)=m1*i nv(s)*G2(i,:)'-1/2*m1* in v(s)*m1'+log(p1);w2(i,2)=m2*i nv(s)*G2(i,:)'-1/2*m2* in v(s)*m2'+log(p2);for j=1:2if w2(i,j)=max(w2(i,:)&j=2n22=n22+1;endendend>> poo=( n11+ n22)/( n1+n2)poo =0.10004、计算交叉误判率>> N11=0;N22=0;>> for k=1: n1A=G1(1:k-1,k+1: n1,:);N1=le ngth(A(:,1);M1= mea n(A,1);s1仁cov(A);S1=(N1-1)*s11+( n2-1)*s2)/(N1+ n2-k);po1= N1/ (n-1);po2=n2/(n-1);for i=1: n1w1(i,1)=M1*i nv(S1)*G1(i,:)'-1/2*M1*i nv(S1)*