3)type=mahalanobis:表示用马氏距离进行判别分析.
§9.4 判别分析示例
例1 某地大气样品污染分类表如下:
在MATLAB中写一个名为opt_linkage_1的M文件: 气体 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 氯 0.056 0.04 0.05 0.045 0.038 0.03 0.034 0.03 0.084 0.085 0.064 0.054 0.048 0.045 0.069 硫化氢 0.084 0.055 0.074 0.05 0.13 0.11 0.095 0.09 0.066 0.076 0.072 0.065 0.089 0.092 0.087 二氧化硫 0.031 0.1 0.041 0.11 0.079 0.07 0.058 0.068 0.029 0.019 0.02 0.022 0.062 0.072 0.027 碳4 0.038 0.11 0.048 0.1 0.17 0.16 0.16 0.18 0.32 0.3 0.25 0.28 0.26 0.2 0.05 环氧氯丙烷 0.0081 0.022 0.0071 0.025 0.058 0.05 0.2 0.22 0.012 0.01 0.028 0.021 0.038 0.035 0.089 环己烷 0.022 0.0073 0.02 0.0063 0.043 0.046 0.029 0.039 0.041 0.04 0.038 0.04 0.036 0.032 0.021 污染分类 1 1 1 1 2 2 1 1 2 2 2 2 2 2 1 在此地某大型化工厂的厂区及临近地区挑选4个有代表性的大气样本取样点,获取数据如下: 气体 样品1 样品2 样品3 样品4 氯 0.052 0.041 0.03 0.074 硫化氢 0.084 0.055 0.112 0.083 二氧化硫 0.021 0.11 0.072 0.105 碳4 0.037 0.11 0.16 0.19 环氧氯丙烷 0.0071 0.021 0.0056 0.02 环己烷 0.022 0.0073 0.021 1 污染分类 求它们的污染分类.
在MATLAB软件包中写一个名为opt_class_1的M文件: training=[0.056 0.084 0.031 0.038 0.0081 0.022
0.04 0.055 0.1 0.11 0.022 0.0073 0.05 0.074 0.041 0.048 0.0071 0.02 0.045 0.05 0.11 0.1 0.025 0.0063 0.038 0.13 0.079 0.17 0.058 0.043 0.03 0.11 0.07 0.16 0.05 0.046 0.034 0.095 0.058 0.16 0.2 0.029 0.03 0.09 0.068 0.18 0.22 0.039
0.084 0.066 0.029 0.32 0.012 0.041 0.085 0.076 0.019 0.3 0.01 0.04 0.064 0.072 0.02 0.25 0.028 0.038 0.054 0.065 0.022 0.28 0.021 0.04 0.048 0.089 0.062 0.26 0.038 0.036 0.045 0.092 0.072 0.2 0.035 0.032 0.069 0.087 0.027 0.05 0.089 0.021];
group=[1;1;1;1;2;2;1;1;2;2;2;2;2;2;1]; sample=[0.052 0.084 0.021 0.037 0.0071 0.022
0.041 0.055 0.11 0.11 0.021 0.0073 0.03 0.112 0.072 0.16 0.0056 0.021 0.074 0.083 0.105 0.19 0.02 1];
class=classify(sample,training,group) 存盘后按F5键执行,得到结果如下:
class = 1 1 2 2
即样品1、样品2为1类污染,样品3、样品4为2类污染.
第八章 聚类分析和判别分析修正版



