夏 鹃 ( Ｒhododendron pulchrum) ，由 Ｒ． indicum 和 Ｒ．eriocarpum ( Hayata) Nakai 杂交或变异而来，属常绿灌木，为晚花杜鹃品种［1］。其种类繁多，花色艳丽，枝叶纤细，且树冠丰满，发枝力强，极耐修剪，特别适合于盆景的制作［2］。目前，在夏鹃的引种及盆景制作的品种选择上，多以人为喜好为主要判断依据，易造成盲目引种或造景效果不佳等问题。因此，亟需一套较客观、系统的定量分析方法，来比较各夏鹃品种的盆景应用价值，以提供指导［3］。
层次分析法( Analytic Hierarchy Process，简称 AHP) ，是运用多因素分级处理来确定权重的方法，是一种定性和定量相结合的系统化、层次化的分析方法［4］，适合于处理多因素问题中各评价指标权重因子的确定，具有良好的有效性、可靠性和实用性［5］。此外，K － Means 聚类法，是一种基于划分的经典聚类算法，该算法的基本目标是将 n 个点的数据集合划分为 k 类，使每一类都具有较高的相似性，是研究( 样品或指标) 分 类 问 题 的 一 种 统 计 分 析 方 法［6 － 7］。 在 应 用 上，K － Means 聚类法多与层次分析法、模糊数学法、灰色系统理论等评价体系相结合使用［8 － 10］。
近几年，层次分析法已被逐渐应用于观赏植物的引种、优良品种选择、景观评价及休闲农业等相关研究领域［11 － 14］，在紫薇［15］、芍药［16］、切花菊［17］、观赏甘薯［18］等品种中已有相关报道。AHP 在杜鹃的种质资源及观赏性等方面虽有相关研究［19 － 22］，但在盆景应用方面却鲜见相关报道。本研究采用AHP 及 K － Means 聚类分析法，以引种的 83 个夏鹃品种在盆景应用上的特点为评价重点，利用周年观测数据、市场调查并结合专家群体判断，对引种的夏鹃品种进行评价及分级，筛选出适合于盆景制作应用的夏鹃品种，为夏鹃品种的应用提供科学依据。
1 材料与方法
1． 1 试验区概况
夏鹃引种试验基地位于福建省东北部的周宁县卢下洋66． 67 hm2花卉区内，地处鹫峰山脉东麓，地势由西北向东南倾斜，平均海拔 800 m。地处 27°712″N，119°1948″E，属中亚热带季风山地气候，冬长夏短，云雾多，雨量充沛，日照少。年平均相对湿度为 83% ，年平均气温为 13 ～ 17℃ ，最冷的 1 月年平均气温为 3 ～ 7 ℃ ，极端最高气温为 35． 6 ℃ ，极端最低气温为 － 8． 9 ℃ ( 1983 年 12 月 31 日) ，年降水量在 1 800 ～2 200 mm，比较适合夏鹃的种植和生长。
1． 2 材料
以福建天蓝蓝生态有限公司引进的 83 个 3 年生的地栽夏鹃品种作为本研究的试验材料，杜鹃花品种及综合评价见表 1。
1． 3 方法
1． 3． 1 评价方法 层次分析法、K － Means 均值聚类分析法、最小欧式距离法。
1． 3． 2 调查方法 文献研究法、实地观测法、问卷调查法( 邮件咨询) 。
1． 3． 3 测量方法 参考徐忠等的方法略作修改。
以生长势中等，能够代表品种整体生长状况的单株作为观测对象，并挂牌标示。在相同生境条件下，每个夏鹃品种选择 3 ～ 5 株植株作为观测对象; 对于不同栽培生境，每一生境条件下选择 3 ～ 5 株作为观测对象。观测时候为 2014 年 6 月至 2016 年 6 月。4 月下旬至 6 月中旬 1 ～ 2 d 观测 1 次，其他时间 1 周观测 1 次。花色特征的观测以英国皇家园艺协会( ＲHS) 出版的比色卡( ＲHS color chart) 为标准，花朵盛开时以直尺测量 2 个垂直方向的长度，取其均值为花径。其中表型性状测量方法参照杜鹃 DUS 测试指南［23］的相关测量部位及测量方法。
由于夏鹃病虫害种类较多，同时受季节和气候等因素影响，难以在短期内对夏鹃的病虫害详细作出评价，本试验仅调查了夏鹃最主要的害虫( 军配虫) 情况，以夏鹃花期军配虫危害情况，作为夏鹃品种盆景应用病虫害抗性的评价指标。病虫害抗性的计算方法参考张冬菊的方法［24］，用品种病虫害指数 = 花期品种军配虫危害株数/品种总株数来衡量。目前，夏鹃盆景多以嫁接方式进行造型，因此，繁殖难易程度用品种嫁接成活率 = 品种嫁接成活数/品种总嫁接数 × 100% 来衡量。
1． 3． 4 调查问卷的设计 评价指标的确定。本试验邀请园艺、园林等相关领域的专家及从事杜鹃盆景行业的技术人员，填写专家意见征询表，以确定选取哪些指标作为评价指标。共发放调查问卷 30 份，收回问卷 30 份，回收率为 100% 。
评价指标重要度的确定。共包括 2 份问卷。问卷 1: 邀请园艺、园林、盆景造型等专业背景的教师及学生，采用 1 ～ 9标度法比较指标间的相对重要性，量化各项指标并构建判断矩阵。共发放调查问卷 70 份，收回问卷 70 份，回收率为100% 。
问卷 2: 邀请园艺及非园艺专业人士各 100 人，对 14项指标，参照评分标准及相关资料进行主观感受的赋分，采用1 ～ 5 分制进行打分。其中花径、花期长短、花期迟早、花丰富度、繁殖难易程度、病虫害抗性等指标，主要根据课题组观测的数据进行打分，其他指标主要根据所拍摄的 83 个夏鹃品种相关照片进行打分。共发放调查问卷 200 份，收回问卷 196份，回收率为 98% 。
1． 3． 5 统计分析方法 利用 Microsoft Excel 2003 软件［25 － 26］对 14 个评价指标的权重和各品种的综合评价分值进行计算，制作图表。利用 IBM SPSS Statistics 22． 0 统计软件［27］对 83个夏鹃品种进行聚类分析，划分等级。
2 结果与分析
2． 1 夏鹃品种盆景应用综合评价体系的构建
2． 1． 1 评价指标的确定 通过现场观测、资料调研和专家咨询等方法，最终筛选出花色、花型、花径、花期长短、花期迟早、花丰富度、残花状况、株型、叶形、叶色、繁殖难易程度、生长势、病虫害抗性、抗逆性等 14 个评价指标。
2． 1． 2 评分标准的确定 各评价指标具体评分标准，是在对引种夏鹃的生物学特性、观赏性以及生长适应性充分观察的基础上，结合调研结果及实际观测数据，最终拟订了 5 分制夏鹃品种盆景应用的综合评分标准( 表 2) 。
2． 1． 3 层次分析模型的建立 以筛选出的 14 个与夏鹃盆景应用密切相关的评价指标，并依据其相互关系建立递阶层次结构评价模型。模型分 4 层。目标层 A: 为对夏鹃品种资源的综合评价( 划分其等级) ; 约束层 B: 由花观赏性、整体观赏性、适应性 3 个部分组成，作为对 A 层的约束层; 标准层 C: 由花色、花型、花径等 14 个指标组成; 最底层 D: 由供试的 83 个夏鹃品种组成。由此构成了由总目标、主要性状、评价指标等组成的多层次评价体系，夏娟品种盆景评价指标及模型见图 1。
2． 1． 4 判断矩阵的构建及一致性检验 在充分利用相关专家、技术人员的经验和判断的基础上，对夏鹃品种定量和非定量因素进行统一度量，构建判断矩阵。根据已构建的层次模型，以发放问卷形式，邀请相关专家依据各因子间的相对重要性标度( 表 3) ，进行对立打分，后用求平均值四舍五入取整法建立 A － B，B1－ C，B2－ C，B3－ C 共 4 个矩阵，并计算出对应的判断矩阵的最大特征值( λmax) 和对应特征向量( W) ，结果见表 4 ～ 表 7。
根据 AHP 理论，当判断矩阵具有满意的一致性时，其特征根大于 n，且其余特征根近于 0［28］。所以，为保证结论的可靠性和合理性，需对判断矩阵进行一致性检验。根据定义:CＲ = CI / ＲI，其中一致性指标 CI = ( λmax－ n) / ( n － 1 ) ，平均随机一致性指标( ＲI) 系数取值见表 8。当 CＲ 值小于 0． 1，说明该矩阵具有满意一致性，否则需对判断矩阵进行调整［29］。2． 2 各评价指标层次单排序。
2． 2 各评价指标层次单排序
从表 4 可以看出，夏鹃品种应用综合评价体系中约束层三类因子权重值分别为花观赏性( B1) 0． 649 1，整体观赏性( B2) 0． 279 0，适应性( B3) 0． 071 9。其中花观赏性所占权重最大，为综合评价的重要指标，表明在夏鹃盆景适宜品种筛选过程中，首要目标是选择花观赏性佳的优异植株，其次要考虑到植株整体观赏性和适应性。由表 5 ～ 表 7 可以看出，花观赏性评价因子 B1－ C 层的各个具体评价指标的权重值大小排序为:花色( 0． 365 0) ＞ 花型( 0． 256 5) ＞ 花丰富度( 0． 169 4) ＞ 花期长短( 0． 091 4) ＞ 花径( 0． 055 8) ＞ 花期迟早( 0． 034 7) ＞ 残花状况( 0． 027 1) ，表明花色、花型和花丰富度是评价花观赏性的主要指标，花期长短、花径、花期迟早和残花状况相对而言。
只能作为辅助评价指标; 在整体观赏性评价因子 B2－ C 中，株型所占权重最大，叶色次之，叶形最小，权重值分别为:0. 661 2、0． 271 8、0． 067 0; 适应性评价因子 B3－ C 层的各个具体评价指标的权重值大小排序为: 生长势( 0． 531 8) ＞ 繁殖难易程 度 ( 0． 270 2) ＞ 病 虫 害 抗 性 ( 0． 122 1 ) ＞ 抗 逆 性( 0． 076 0) ，表明在盆景应用中对夏鹃的生长势强弱和繁殖难易程度要求较高，其次是病虫害抗性和抗逆性。评价结果与实践相比，具有较高的一致性，说明该模型及评价方法具有较佳的适合度、可信度。
2． 3 各评价指标层次总排序
同一层次所有因素对于最高层次的相对重要性权重值的排序为层次总排序，通过计算出各评价指标( C) 相对于所属性状( B) 的加权值后，再与该性状( B) 相对于总目标( A) 所得的权值进行加权综合，即可计算出各评价指标因素( C) 相对于总目标( A) 的权值，进而得到总排序［30］。根据层次总排序计算结果可知，在综合评价中，花色、株型、花型和花丰富度这4 个评价指标所占权重较大。其中花色所占的权重最大，为23． 70% ，株 型 占 18． 44% ，花 型 占 16． 65% ，花 丰 富 度 占10． 99% ，这几个评价指标所占权重总和为 69． 79% ( 表 9 ) 。
结果表明花色、株型、花型及花丰富度是夏鹃品种盆景应用价值综合评价的决定性因素; 叶色、花期长短、生长势、花径、花期迟早及繁殖难易程度等所占权重次之，分别为 7． 58% 、5. 93% 、3． 82% 、3． 62% 、2． 25% 、1． 94% ，总权重为 25． 16% ，为夏鹃品种盆景应用价值综合评价的重要因素; 而叶形、残花状况、病虫害抗性、抗逆性等所占权重最小，分别为 1． 87% 、1． 76% 、0． 88% 、0． 55% ，权重总和约为 5． 05% ，为评价夏鹃品种盆景应用价值综合评价的一般评价因素。通过建立AHP 评价体系，对影响夏鹃品种盆景应用的因素进行综合评价，可为夏鹃的盆景应用提供指导，避免因盲目开发而造成所制作盆景的景观效果差或植株生长不良等问题。
2． 4 夏鹃品种盆景应用综合评价得分及等级的确定
依据 14 项评价指标的评分标准，对 83 个供试夏鹃品种进行逐项打分，再根据各指标的综合权重值，计算出各品种的综合分值。进而应用 K － Means 聚类法对各品种的综合分值进行聚类分析，获得的聚类中心分别为 3． 824 9、3． 159 3、2． 304 8，再利用最小欧式距离法结合聚类中心值进行等级划分。最终将供试的夏鹃品种划分为 3 级，其中，Ⅰ级得分 ＞3． 51，为盆景应用价值最高的夏鹃品种，共 19 个品种; Ⅱ级得分 3． 51 ～ 2． 77，为盆景应用价值一般的夏鹃品种，共 33 个品种; Ⅲ级得分 ＜ 2． 77，为盆景应用价值较低的夏鹃品种，共 31个品种。
2． 5 综合评价结果的分析与验证
盆景应用价值较高、一般、较低的 3 个等级所占比率分别为 23% 、40% 、37% 。
Ⅰ级夏鹃品种包括了五宝绿珠、丽子、日之丸等综合评分值最高的品种，这些Ⅰ级夏鹃品种的综合性状优良，具有极佳的应用前景，可在今后盆景应用中作较多考虑。其特点是花型饱满，花色鲜艳有光泽，着花量大，同时株型紧凑，生长势强，较耐修剪，也具有较佳的病虫害抗性及抗逆性强。
Ⅱ级夏鹃品种为和、晃明、夏锦袍等综合评分值居中的品种，这些夏鹃品种的花观赏性与整体观赏性较Ⅰ级夏鹃品种差，生长势较弱，适应性一般。在综合评价中，表现为某些方面指标值较高，其他指标值较低，因此，这些Ⅱ级夏鹃品种可作为盆景应用的补充材料考虑。彩、光华、新日光等的综合评分值最低，为Ⅲ级夏鹃品种，这些品种一般存在较大缺陷，观赏性一般，生长缓慢，植株松散，叶片色泽变淡，不适于盆景应用，可作为种质资源保存。
从综合评价结果可以看出，本研究所选用的评价因子基本能体现夏鹃的品种特性，是夏鹃品种在盆景应用中的主要影响因素。与传统方法相比，本研究综合评价体系不仅衡量了植株的观赏价值，而且评估了各品种的抗性，能较全面地反映夏鹃在盆景应用中的综合应用价值。在评价体系中花观赏性和整体观赏性所占权重最大，这既符合观赏效果上对植株花色、花型、花丰富度和株型的要求，又符合盆景造型过程中对植株生长势强、嫁接后易于成活等的要求。
同时，由于盆景具有精细养护的特点，对适应性方面的要求会明显低于花观赏性和整体观赏性。在权重总排序中花色、株型、花型、花丰富度、叶色、花期长短、生长势等因子的权重值较大，体现了本次名优品种筛选的目的，符合预期的结果。因此，本模型具有一定的借鉴意义。
3 讨论
本研究参照国内外观赏植物资源评价指标，以夏鹃的夏季城市盆景应用为总目标，运用层次分析法构建夏鹃品种应用综合评价模型，运用层次分析法构建夏鹃品种盆景应用综合评价模型，可以使视觉美感、生态应用适应性、功能性等抽象性评价因素定量化，将定性与定量相结合，大大提高了评价的客观性和有效性。
在评价模型的建立过程中，通过多样本( 样本数 ＞ 30) 调查后取平均值的办法，使因子间的相对重要性标度的人为主观因素降到最小。此外，不同栽培条件及气候环境差异等因素对夏鹃品种综合观赏价值的表现有一定影响，如夏鹃的生长发育及观赏特性等性状在偏酸性土壤中能得到较充分的展现，但在偏碱性土壤中则相反，尤其是适应性较弱的夏鹃品种的观赏价值会随着栽培时间的延长而逐渐降低［21］。
因此，夏鹃品种的评价应在栽培条件、栽培环境、引种栽培时间等相对较一致的情况下进行，以提高评价的科学性。