区间值直觉模糊麦克劳林聚合算子的近红外反射波段划分方法
任维佳1,2,3, 杜向军1,3, 杜玉琴4,*, 孙荣禄1,3, 李学良1,3
1. 天津工业大学机械工程学院, 天津 300387
2. 天津科技大学机械工程学院, 天津 300222
3. 天津市现代机电装备技术重点实验室, 天津 300387
4. 中国社会科学院大学经济学院, 北京 102401
*通讯作者 e-mail: duyugin@ucass.edu.cn

作者简介: 任维佳, 1993年生, 天津工业大学机械工程学院博士研究生 e-mail: rwj_2020@163.com;
杜向军, 1980年生, 天津工业大学机械工程学院博士研究生 e-mail: dxjun216@163.com;任维佳, 杜向军: 并列第一作者

摘要

采用不同波段光谱对不同类异性纤维进行识别, 可有效剔除异性纤维, 提高检出率。 针对波段传统划分方法存在的评价指标单一性问题, 研究不同异性纤维多种属性间相互影响, 结合多属性群决策(MAGDM)方法的优势, 提出利用多属性群决策方法实现棉花中异性纤维最佳检测波段的选取, 根据不同异性纤维多种属性指标的关联性, 综合考虑多种属性指标, 确定波间可分度、 相关性和ABS指标作为异性纤维属性的评价指标。 首先, 针对多属性群决策方法中评价准则不准确问题, 构建评价准则函数线性方程组, 使其增广矩阵的秩等于未知量个数, 保证方程组具有唯一解, 从而提高决策结果的准确性。 接着, 使用幂等(PA)算子消除不合理评价信息数据值对决策结果的负面影响, 结合麦克劳林对称平均(MSM)算子综合考量多个数据间的相互关系, 推导加权区间直觉模糊幂等麦克劳林对称平均(WIVIFPMSM)聚合算子, 拓展TOPSIS方法确定异性纤维权重信息, 聚合不同异性纤维多种属性的评价信息, 根据所建立的评价准则确定决策结果, 进而构建基于区间值直觉模糊集(IVIFS)的多属性决策方法, 实现多种异性纤维多种属性指标的最佳波段选取。 将区间值直觉模糊幂等麦克劳林聚合算子与类间可分性波段选择(ISBC)方法和自适应选择(ABS)方法进行了对比, 分析了不同波段划分方法对结果的影响, 总结了现有研究中存在的问题和不足。 为提高多属性群决策方法的决策精度, 分析了参数 k对决策结果的影响, 证明区间值直觉模糊幂等麦克劳林聚合方法稳定性更好, 为复杂环境下异性纤维波段划分研究提供了新思路。 最后, 通过实验验证近红外波段 W3: 780~1 100 nm为最佳检测波段。 该方法对波段选取领域的理论扩展和多属性群决策方法的应用具有一定的指导意义。

关键词: 异性纤维; 反射率光谱; WIVIFPMSM聚合算子; 多属性决策; 波段选择
中图分类号:TP391.41 文献标志码:A
Near Infrared Reflection Band Division Method Based on Interval-Valued Intuitionistic Fuzzy Maclaurin Symmetric Mean Aggregation Operator
REN Wei-jia1,2,3, DU Xiang-jun1,3, DU Yu-qin4,*, SUN Rong-lu1,3, LI Xue-liang1,3
1. School of Mechanical Engineering, Tiangong University, Tianjin 300387, China
2. College of Mechanical Engineering, Tianjin University of Science and Technology, Tianjin 300222, China
3. Tianjin Key Laboratory of Advanced Mechatronics Equipment Technology, Tiangong University, Tianjin 300387, China
4. School of Economics, University of Chinese Academy of Social Sciences, Beijing 102401, China
*Corresponding author

REN Wei-jia and DU Xiang-jun: joint first authors

Abstract

Using different wavelength spectra to identify different types of foreign fibers can effectively eliminate foreign fibers and increase the detection rate. Given the problem of a single evaluation index existing in traditional band division methods, the interaction of attribute indexes of different fibers is studied, combined with the advantages of the multi-attribute group decision-making (MAGDM) method, this paper proposes to use the MAGDM method to realize the selection of the optimal detection band for foreign fibers in cotton. According to the relation of attribute indexes of different fibers, the inter-class separability, correlation and ABS index are determined as attribute evaluation indexes. Firstly, to solve the problem of inaccurate evaluation criteria in the MAGDM method, a system of evaluation criteria function linear equations is constructed so that the rank of the augmented matrix is equal to the number of unknowns, ensuring that the equation system has a unique solution, thereby improving the accuracy of the decision result. Next, the power mean (PA) operator is used to eliminate the adverse effects of unreasonable evaluation information values on decision results, combined with the Maclaurin symmetric mean (MSM) operator to comprehensively consider the relationship between input arguments, deriving the weighted interval-valued intuitionistic fuzzy power Maclaurin symmetric mean (WIVIFPMSM) aggregation operator.Then, the TOPSIS method is used to determine the weight information of foreign fibers, the evaluation information of various attributes of different foreign fibers is aggregated, and the decision results are chosenaccording to the established evaluation criteria. Thus, a MAGDM method based on interval-valued intuitionistic fuzzy sets (IVIFSs) is constructed to realize the optimal band selection of various attributes of foreign fibers. Moreover, the WIVIFPMSM aggregation operator is compared with the inter-class separability band selection (ISBC) method and adaptive band selection (ABS) method, the influence of different band division methods on the results are analyzed, and the existing problems and deficiencies in existing research are summarised. To improve the decision accuracy of the MAGDM method, parameter ks influence on decision results is analysed, and it is proved that the IVIFPMSM aggregation method has better stability, which provides a new idea for the study of band division of foreign fibers in complex environments. Finally, it is verified through experiments that the near-infrared band W3: 780~1 100 nm is the optimal detection band. In addition, this paper has a specificguiding significance for the theoretical extension of band selection and the application of MAGDM methods.

Keyword: Foreign fiber; Reflectance spectrum; WIVIFPMSM operator; Multi-attribute decision making; Band selection
引言

“ 十四五” 规划的棉花行业战略方向研究建设报告中要求, 棉花质量检测要进一步加强。 棉花中的异性纤维(简称异纤)不仅会造成纱线的断裂, 更影响织物的质量, 造成纺织企业巨大经济损失[1]。 异性纤维是指棉花中混入的非棉纤维和非本色棉纤维等有害物质。 目前主要采用图像处理的方法[2, 3]实现棉花中异纤的检测, 利用特定光源照射棉花, 通过将CCD采集到的图像进行数据处理, 找出与棉纤维的差异识别异纤。 异性纤维种类繁多, 主要有羽毛、 塑料薄膜、 化学纤维、 丙纶丝。 采用不同波段的光源对异纤检测识别精度和分类效果具有显著差别。

波段选择是通过研究波段间的相对距离、 可分度和相关性, 划分出可分性大、 冗余小的最佳波段; 常用方法为类间可分性和信息量划分方法。 类间可分性波段选择是通过分析波段间相对距离研究其可分度。 Bao等[4]分析单因素方差和Jeffries-Matusita(J-M)距离光谱一致性检验, 并利用J-M距离选取最佳波段; Charifi等[5]利用Bhattacharyya距离进行可分性检验, 根据可分度选择波段组合; 杨文柱等[6]利用波段间的类间可分性, 认定780~1 800 nm的近红外波段为最佳检测波段。 郏东耀[7]采用基于区域信息相关度权值小波分析算法将多个波段的信息进行融合, 提出了一种利用多波段光谱信息融合成像检测的新方法。 信息量波段选择方法分析波段相关性, 去除波间冗余量, 搜寻信息量最大和相关性最小的波段。 Jiang等[8]采用最小冗余最大相关性(mRMR)算法选择最优波长, 提高棉花中异纤的分类效果; Wang等[9]提出了一种基于自适应子空间划分策略(ASPS)波段选择方法; 杨赛等[10]发现平均光谱和熵联合的分类精度高于单一特征模型, 为波段融合方法提供了理论基础; 马梦茹等[11]利用OIF和ABS方法进行最佳波段选取, ABS比OIF方法提升精度9.21%。

类间可分性和信息量划分方法均未考虑综合划分指标, 当一种指标达到最优时其他指标无法达到最佳效果, 影响了最终结果。 本工作引入多属性决策法, 利用MSM算子[12]综合考量多个输入参数交互关系, 使用PA算子[13]削弱极端输入参数的优势, 提出WIVIFPMSM聚合算子[14, 15, 16], 结合TOPSIS确定各类异性纤维权重向量[17], 构建最佳波段的多属性决策框架, 实现最佳波段的选取。 并将所提方法与类间可分性波段选择(ISBC)和自适应选择(ABS)算法的实验结果进行对比, 验证该方法的合理性和有效性。

1 实验部分
1.1 仪器与样本

实验设备为日立UH4150紫外-可见近红外分光光度计, 可以实现异纤的全光谱检测, 选用的光谱检测范围为240~2 600 nm, 采样间隔为2 nm, 配套软件是UV Solutions 4.2, 参比和样品放置处都要放置Al2O3白板, 异纤使用基地和夹片夹持。 另外, 处理器采用Intel(R)i7-7400CPU@3.00GH, RAM大小为16 GB。 实验样本来源于湖北谷城银纺公司等全国10个棉纺企业实际生产中剔除的异纤, 经初步检查有布块、 树叶、 白纸、 羊毛、 羽毛、 白色泡沫、 麻绳、 黑色塑料薄膜、 头发、 白色塑料薄膜、 丙纶丝、 化学纤维等多种异纤, 而羽毛、 塑料薄膜、 化学纤维、 丙纶丝含量较多, 作为主要研究对象[17]

1.2 数据采集

波段选择方法一般用一组测量数据测验, 缺乏全面性。 针对12类异纤进行光谱数据采集, 每类异纤采集10组数据信息, 数据采集过程中按颜色、 厚度、 疏密度对样本进行调整[8]。 研究光谱信息的变化趋势, 分析数据的可靠性、 稳定性, 构造一组合理的综合异纤光谱数据, 以确保波段选择的合理性和科学性。

1.3 波段划分方法

1.3.1 基于类间可分性波段选择算法

为实现波段划分, 研究各类异纤到棉花之间的距离, 采用欧式平均距离进行数据计算。 基于波段间距离, 结合各类异纤的可分度, 实现波段划分。 平均欧式距离公式为

dis(V, V')=1ni=1n(vi-v'i)2(1)

式(1)中, viv'i表示目标光谱信息, dis(V, V')表示目标间距离。

根据异纤光谱信息的分布及波动情况, 定义光谱反射率的均值 fik和样本标准差 Sik

fik=1ni=1nRik(λ)(2)

Sik=1n-1i=1n[Rik(λ)-fik]2(3)

建立第i种与第j种异纤反射光谱可区分度, 模型为

Ji, jk=α[dis(X, Y)]β(Sik+Sjk)(4)

式(4)中, α 为距离系数, β 为方差系数, 通过调整α β 使可区分函数值 Ji, jk基本位于[0, 1], 可区分函数随着dis(X, Y)增加而增大。 Ji, jk越大说明可分度越大, 有利于异纤检测。

1.3.2 自适应波段选择算法

传统自适应波段选择算法中仅考量了相邻波段间的相关性, 未对其他波段间的相关性做出说明; 本工作对自适应波段选择算法进行改进, 引入子波段平均相关性完善自适应指标, 进而建立ABS指数Ii

Ii=Si(Ri-1, i, Ri, i+1, R̅i)/3(5)

Rij=i=1mj=1n[wi-w̅i][wj-w̅j](i=1mj=1n[wi-w̅i]2)i=1mj=1n[wj-w̅j]2(6)

w̅i=i=mnwi(n-m)+1(7)

其中, Ri, j为波段i与波段j的相关系数, R̅i为波段i与非相邻波段相关系数的均值, w̅i为波段i的均值。

1.3.3 区间值直觉模糊幂等麦考劳林对称平均聚合算子

为解决多种异纤的最佳波段选取问题, 综合考虑多类属性评价指标, 引入多属性群决策方法。 由于异纤光谱特性中的吸收率和反射率之和小于1, 而IVIFS的适用范围完全满足, 同时考虑到输入参量间关系, 利用PA算子消若极端输入参数产生的负面影响, 使用MSM算子综合考量数据间的交互关系, 并将PA和MSM算子融入到区间值直觉模糊框架中, 提出基于WIVIFPMSM算子的异纤波段综合因素确定方法, 实现最佳波段选取。

令(α 1, α 2, …, α n)为一个非负实数集, α i=([ uiL, uiU], [ viU, viL]), α i> 0, i=1, 2, …, n, WIVIFPMSM聚合算子的映射定义为

WIVIFPMSM(k)(α1, α2, , αn)=1i1i2inj=1knηijαijCnk1k=1-1i1i2ikn1-j=1k1-(1-uijL)nηij1Cnk1k, 1-1i1i2ikn1-j=1k1-(1-uijU)nηij1Cnk1k1-1-1i1i2ikn1-j=1k(1-(vijL)nηij)1Cnk1k, 1-1-1i1i2ikn1-j=1k(1-(vijU)nηij)1Cnk1k(8)

式(8)中, k是一个非负整数, i1, i2, , ik1, 2, , n的任意k阶组合, ηij=wi(1-T(αij))m=1nwm(1-T(αm)), 0wm1, m=1nwm=1, T(αi)=j=1jinSup(αi, αj), Sup(αi, αj)=1-D(αi, αj), Cnk是二项式系数。

分析两个区间值直觉模糊数(IVIFNs), 当隶属度区间和非隶属度区间完全一致, 两者才能够相等。 仅靠得分函数S(α i)和精度函数H(α i)不能保证评价准则结果的准确性, 无法满足评价函数个数等于未知参数个数的要求, 方程组的解不唯一。 考虑区间范围对结果的影响, 基于区间Pythagorean模糊信息进行评价函数诱导[17], 其中, 隶属度范围越大说明非隶属度范围越小、 不确定性越大, 引入IVIFNs的隶属度不确定函数Mu(α i), 结合犹豫度π (ai)=( πiL, πiU), 提出犹豫度不确定函数Hu(α i), 具体公式为

S(αi)=12(μiL+μiU-viL-viU)(9)

H(αi)=12(μiL+μiU+viL+viU)(10)

Mu(αi)=uiU-uiL-(viU-viL)(11)

Hu(αi)=uiU-uiL+viU-viL(12)

对于任意一组评价函数S(α i)、 H(α i)、 Mu(α i)和Hu(α i), 构造非齐次线性方程组, 通过研究方程组解的情况, 判定评价准则的严谨性, 为构建精准的评价准则提供前提, 有利于决策结果的准确性。 增广矩阵如式(13)所示, 增广矩阵的秩Rank(A)=4, 等于变量 μiLμiUviLviU的个数, 因此评价函数方程组有唯一解, 可以区分任意两个IVIFNs。

A=11-1-12S0202Mu-2S00222H-2S00042H+Hu-Mu-2S(13)

根据IVIFNs的四个指标函数S(α i)、 H(α i)、 Mu(α i)和Hu(α i), 定义任意两个IVIFNs的评价准则: 若S(α 1)> S(α 2), 则α 1> α 2; 若S(α 1)< S(α 2), 则α 1< α 2; 若S(α 1)=S(α 2), 使用H(α i)函数比较: 若H(α 1)> H(α 2), 则α 1> α 2; 若H(α 1)< H(α 2), 则α 1< α 2; 若H(α 1)=H(α 2), 使用Mu(α i)函数比较: 若Mu(α 1)> Mu(α 2), 则α 1< α 2; 若Mu(α 1)< Mu(α 2), 则α 1> α 2; 若Mu(α 1)=Mu(α 2), 使用Hu(α i)函数比较: 若Hu(α 1)> Hu(α 2), 则α 1< α 2; 若Hu(α 1)< Hu(α 2), 则α 1> α 2; 若Hu(α 1)=Hu(α 2), 则α 1=α 2

1.3.4 拓展优劣解距离法

优劣解距离法(TOPSIS)[16]可实现多输入参数的综合评价, 为了确定出异纤光谱决策信息的相关权重分布, 将TOPSIS方法进行扩展, 与区间值直觉模糊集结合, 诱导相应决策数据下的理想矩阵和贴近度。 令X={X1, X2, …, Xm}为方案集, A={A1, A2, …, An}为属性集, D={D1, D2, …, Ds}为决策集合。 任意XiAj的评价信息为 aijk=([ uijLk, uijUk], [ vijUk, vijLk]), IVIFS决策信息简记为Dm× n

根据决策矩阵 Dm×n获得正理想矩阵D* 和3个负理想矩阵D-D+和(D* )c, 矩阵为

Dij* =(Dij* )m×n=([μijL* , μijU* ], [vijL* , vijU* ])m×n(14)

D-=(Dij-)m×n=([μijL-, μijU-], [νijL-, νijU-])m×n(15)

D+=(Dij+)m×n=([μijL+, μijU+], [νijL+, νijU+])m×n(16)

(D* )c=((Dij* )c)m×n=([νijL* , νijU* ], [μijL* , μijU* ])m×n(17)

式中

μijL* =1-k=1t(1-(μijL(k))2)1/t1/2, μijU* =1-k=1t(1-(μijU(k))2)1/t1/2, vijL* =k=1t(μijL(k))1/t, vijU* =k=1t(μU(k)ij)1/t, μijL-=mink=1, 2, , t{μijL(k)}, μijU-=mink=1, 2, , t{μijU(k)}, vijL-=maxk=1, 2, , t{vijL(k)}, vijU-=maxk=1, 2, , t{vijU(k)}, μijL+=maxk=1, 2, , t{μijL(k)}, μijU+=maxk=1, 2, , t{μijU(k)}, vijL+=mink=1, 2, , t{vijL(k)}

根据正理想矩阵和负理想矩阵, 计算其贴近度。

RCk=NprojD* (D(k))NprojD* (D(k))+NprojD-(D(k))+NprojD+(D(k))+Nproj(D* )C(D(k))(18)

决策矩阵离正理想矩阵越近, 离负理想矩阵就越远, 相对贴近度就越大, 通过值的归一化计算权重矩阵。

wk=RCkk=1tRCk(19)

2 结果与讨论
2.1 异性纤维波段划分

对采集的异性纤维波段数据进行统计分析, 使用数据的均值作为有效数据, 以保证后续实验的可靠性, 异纤光谱具体信息如图1所示。 丙纶丝、 化学纤维和白色塑料薄膜的反射率明显高于其他异性纤维, 并且在波段380~1 100 nm有较好的辨识率, 但是整体差异性较小。 波段240~715 nm棉花与布块、 树叶的反射率差在5%之内, 可识别性小不利于检测。 黑色塑料薄膜和头发的反射率收敛较慢, 在波段240~1 700 nm持续增长。 另外, 可以看到棉花与主要异性纤维的可辨识性明显, 有利于检测识别。 而棉花与羊毛在波段300~700 nm可分性和棉花与布块在波段300~715 nm和波段2 200~2 450 nm的可分性较差。 根据光谱学划分方法, 将波段划分为: W1: 240~380 nm、 W2: 380~780 nm、 W3: 780~1 100 nm、 W4: 1 100~1 700 nm、 W5: 1 700~2 100 nm、 W6: 2 100~2 600 nm。

图1 棉花与异性纤维反射率图像Fig. 1 Reflectance spectra of cotton and foreign fibers

2.2 基于类间可分的异性纤维波段划分方法

为提高异纤检测效果, 对异纤与棉花间的可区分性进行分析, 如图2所示, 羽毛、 塑料薄膜、 丙纶丝以及化学纤维均在波段W3取得最大可分度, 其中塑料薄膜的可分度为1, 说明塑料薄膜在波段W3达到完全可分的效果。 另外, 主要异纤在波段W3的可分效果明显高于波段W1W2W4W5, 有利于棉花中异纤检测。

图2 棉花与主要异性纤维可分度Fig.2 Separabilities of cotton and main foreign fibers

2.3 基于信息量的异性纤维波段划分方法

自适应波段选择算法通过考察波段间的相关性, 对信息量做出保证, 同时该方法不需要预先确定波段个数, 提供多种波段划分选择。 异纤经过ABS算法处理获得各波段ABS指标图像, 如图3所示。 由于波段间的信息量与波段数据的标准差有关, 在检测初期, 反射率发生剧烈变化, 造成初始阶段W1的波间信息较大。 当反射率趋于稳定, 信息量在波段W3达到最大, 另外, 塑料薄膜在波段W3的信息量大于在波段W1处的信息量。 因此, 在保证数据稳定的前提下, 波段W3各类异性纤维的信息量达到最大化。

图3 主要异性纤维在波段W1W6信息量Fig.3 The main foreign fibers information in W1W6 wave bands

2.4 基于多属性决策方法在异性纤维波段划分中的应用

针对棉花中异纤检测最佳波段选取问题, 采用WIVIFPMSM算子聚合多种异纤的波段决策信息, 拓展TOPSIS方法确定异纤权重向量, 利用评价准则确定最佳波段信息。 使用羽毛、 塑料薄膜、 丙纶丝、 化学纤维四类异纤, 构建F={F1, F2, F3, F4}, 对6个子波段W={W1, W2, W3, W4, W5, W6}, 在3种属性A={A1, A2, A3}下进行评价, 其中A1表示类间可分度, A2表示相关性, A3表示信息量, 异纤权重向量为w'={w'1, w'2, w'3, w'4, w'5, w'6}, 属性权重为w={w1, w2, w3}, 考虑到属性的相对独立性, 设定w={0.333, 0.333, 0.333}。 利用评价信息获得决策矩阵, 使用评价值作为直觉模糊集的隶属度, 非隶属度取直觉模糊环境下的最大值, 由于数据采集及处理过程中存在误差, 将直觉模糊集转变为区间值直觉模糊集, 如表1表4所示。

表1 羽毛决策矩阵 Fm×n1 Table 1 The decision matrix Fm×n1 of the feather
表2 塑料薄膜决策矩阵 Fm×n2 Table 2 The decision matrix Fm×n2 of the plastic film
表3 丙纶丝决策矩阵 Fm×n3 Table 3 The decision matrix Fm×n3 of the PP Yarn
表4 化学纤维决策矩阵 Fm×n4 Table 4 The decision matrix Fm×n4 of the chemical fiber

考虑到异纤之间的相互作用, 避免极端数据的负面影响, 利用TOPSIS方法确定各类异性纤维的权重向量。 使用决策矩阵, 计算决策矩阵F(k)在正理想矩阵F* 和负理想矩阵F+F-、 (F* )c上的投影, 其结果如表5所示。

表5 Fm×nk在理想决策矩阵上的投影 Table 5 The projection of Fm×nk in ideal decision matrix

根据式(19)计算各类异纤的贴近度RC1=0.281, RC2=0.319, RC3=0.331, RC4=0.328, 将贴近度进行归一化, 获得权重w'1=0.223, w'2=0.253, w'3=0.263, w'4=0.261, 所以异纤权重向量为w'1=(0.223, 0.253, 0.263, 0.261)。

计算支持矩阵Sup( αijd, αikl), 为方便将 αijdαikl的支持矩阵Sup( αijd, αikl)记为Sdl

S12=S21=0.9990.8830.7930.9820.8160.9450.3280.9580.7930.8920.9550.9600.7590.9900.9050.3560.9950.968S13=S31=0.9980.8500.9740.5300.8300.8890.3460.9640.9920.5200.9990.9960.3490.8650.9100.3560.9990.941S14=S41=0.3751.0000.9990.3780.6960.9080.5810.9640.9680.3340.9630.9750.5970.9440.9940.9990.9690.857S23=S32=0.9970.9670.7670.5120.9860.9440.9820.9950.8010.4110.9560.9650.5900.8550.9951.0000.9960.909S24=S42=1.0000.8530.6500.5790.7610.9510.9940.9950.7680.4730.9910.9920.6200.6780.9970.9810.9950.967S34=S43=0.9970.8190.8830.9330.7750.8950.9890.9990.9680.9390.9650.9720.9710.8310.9980.9810.9990.942

根据支持函数Sup( αijd, αikl)求解加权支持矩阵T( αijd), 计算幂权重矩阵 λijd, 记作λ d

λ1=0.2140.2150.2190.2140.2150.2140.1770.2120.2170.2160.2140.2140.1990.2160.2090.1790.2140.215λ2=0.2530.2540.2430.2510.2530.2540.2660.2530.2420.2450.2520.2520.2570.2550.2540.2660.2520.252λ3=0.2610.2600.2680.2590.2630.2580.2750.2620.2660.2610.2610.2610.2660.2050.2630.2740.2610.259λ4=0.2600.2580.2570.2630.2560.2610.2730.2600.2620.2650.2600.2600.2670.2500.2610.2720.2600.261

利用IVIFPMSM算子将决策矩阵D1, D2, D3, D4聚合成一个IVIFS决策矩阵D'=(α ij)6× 3, 综合决策矩阵D'表6

表6 综合决策矩阵D' Table 6 The integrated decision matrix D'

计算D'm× n(α ij)中α ij, α ik(j, k=1, 2, 3; jk)支持函数Sup'(α ij, α ik), 简记为S'jk

S'12=S'21=0.6930.9860.1400.6550.7290.175S'13=S'31=0.2240.9610.5860.8810.6100.349S'23=S'32=0.5310.9750.5540.7740.8810.826

根据支持函数Sup'(α ij, α ik)求解加权支持矩阵T'(α ij), 计算幂权重矩阵η ij

ηij=0.3290.3550.3160.3330.3340.3330.3220.3200.3580.3330.3250.3420.3230.3430.3340.3010.3420.357

利用IVIFPMSM算子将决策矩阵D'中的α ij聚合为α i, 聚合矩阵如表7所示。

表7 聚合矩阵 Table 7 The aggregation matrix

计算S(α i)、 H(α i)、 Mu(α i)和Hu(α i), 按评价准则进行排序。 得分函数计算结果为: S(W1)=-0.273, S(W2)=-0.401, S(W3)=-0.074, S(W4)=-0.404, S(W5)=-0.142, S(W6)=-0.466, 根据第一评价准则进行排序, 排序结果为W3> W5> W1> W2> W4> W6, 即最佳波段为近红外波段W3

2.5 敏感度分析

参数k的取值会对最终的决策结果产生影响, 研究不同k值下WIVIFPMSM算子决策结果的敏感性, 如表8所示。 尽管k被赋值为不同的数值, 但最佳的波段均为W3, 不同参数k下决策结果的一致性说明了所提方法具有一定的稳定性。 在实际决策中, 可选择合适的k值, 即根据决策风险的程度选择参数k的自适应值。 对参数k进行敏感性分析, 若优先考虑风险, 则WIVIFPMSM算子的参数k将尽可能大; 如不考虑风险, 则尽量减小参数k

表8 不同k值对决策结果的影响 Table 8 The influences of different k values on decision results
2.6 方法比较

分析基于WIVIFPMSM算子异纤最佳波段选择结果与ISBC和ABS方法的差异, 如表9所示。 ISBC方法表示异纤与棉花的类间可分程度, 当使用该方法时, 4种主要异纤的最佳波段均为W3, 最差波段均为W1; ABS方法可以显示波段所包含的信息量, 波段间的标准差对最终的结果存在影响, 在使用ABS方法进行波段选择时, 只有塑料薄膜的最佳波段为W3, 羽毛、 丙纶丝和化学纤维的最佳波段为W1, 其主要原因是异纤光谱信息在采集初期的波段波动较大, 较大标准差造成了异纤在波段W1包含信息量较大, 但是在波段稳定后, 各类异纤最佳波段均为W3。 基于WIVIFPMSM算子将类间可分度、 相关性和信息量综合考虑, 对异纤的最佳波段进行多方面分析, 发现最佳波段为W3, 与上述ISBC和ABS方法结果一致。 进一步比较最佳决策结果同次等结果的差异突显所提多属性决策方法优势, 发现所提方法较ISBC方法下的主要异纤平均差异提升88.89%, 在不考虑波段W1下所提方法较ABS方法下的主要异纤平均差异提升67.23%。 通过将多属性决策方法与传统划分方法比较, 说明了所提方法的有效性, 且多属性决策方法较传统方法能实现波段信息综合评价, 聚合多种评价信息, 为后续多角度指标的精确波长确定提供一种可行性方案。

表9 多属性决策方法与传统划分方法比较 Table 9 Comparison the decision-making method with the traditional method
3 结论

针对波段传统划分方法指标单一缺陷, 综合考虑可分度、 相关性和信息量等属性评价指标, 提出一种异纤波段多属性决策方法, 并将决策结果与传统波段划分方法结果进行对比, 说明所提方法的合理性和有效性。 研究了传统波段划分方法并进行改进, 通过使用平均欧式距离进行距离度量优化ISBC法, 利于增加目标波段与非相邻波段的相关系数改进基于信息量的ABS法。 诱导WIVIFPMSM模糊聚合算子, 使用拓展TOPSIS法优化决策过程中的权重分布, 构建多属性决策机制, 并将其应用到异纤最佳波段的选取中, 以类间可分性、 相关性和信息量等性质为属性集, 以4种主要异纤为决策集, 实现最佳波段的选取。 通过决策结果的敏感度分析和方法比较体现决策机制优越性, 发现不同k值下的最佳决策结果一致, 说明所提多属性决策方法具有较好稳定性。 将基于WIVIFPMSM算子的决策结果与传统波段划分方法对比, 决策结果一致, 处理结果有效, 然而所提方法能够改善单一波段划分方法信息量的局限性, 从多角度、 多属性进行最佳光谱特性的选取, 为后续精确波长确定提供一种有效的决策方案。

参考文献
[1] REN Wei-jia, DU Yu-hong, ZUO Heng-li, et al(任维佳, 杜玉红, 左恒力, ). Journal of Textile Research(纺织学报), 2021, 42(12): 196. [本文引用:1]
[2] ZHANG Lin, WEI Ping, WU Jian-bo, et al(张林, 韦平, 伍剑波, ). Transactions of the Chinese Society of Agricultural Engineering(农业工程学报), 2016, 32(15): 289. [本文引用:1]
[3] SHI Hong-yu, GUAN Sheng-qi (师红宇, 管声启). Journal of Donghua University(Natural Science) [东华大学学报(自然科学版)], 2016, 42(3): 400. [本文引用:1]
[4] Bao S, Cao C, Chen W, et al. Geocarto International, 2021, 36(2): 226. [本文引用:1]
[5] Charifi R, Es-Sbai N, Zennayi Y, et al. Computers & Geosciences, 2021, 150: 104696. [本文引用:1]
[6] YANG Wen-zhu, LI Dao-liang, WEI Xin-hua, et al(杨文柱, 李道亮, 魏新华, ). Transactions of the Chinese Society of Agricultural Engineering(农业工程学报), 2009, 25(10): 186. [本文引用:1]
[7] JIA Dong-yao, DING Tian-huai (郏东耀, 丁天怀). Journal of Tsinghua University(Science and Technology) [清华大学学报(自然科学版)], 2005, 2: 193. [本文引用:1]
[8] Jiang Y, Li C. Computers and Electronics in Agriculture, 2015, 119: 191. [本文引用:2]
[9] Wang Q, Li Q, Li X. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2019, 12(12): 4940. [本文引用:1]
[10] YANG Sai, ZHU Qi-bing, HUANG Min(杨赛, 朱启兵, 黄敏). Spectroscopy and Spectral Analysis(光谱学与光谱分析), 2017, 37(3): 990. [本文引用:1]
[11] MA Meng-ru, ZHANG Yong-bin, WANG Yi-dan (马梦茹, 张永彬, 王奕丹). Journal of North China University of Science and Technology(Natural Science Edition) [华北理工大学学报(自然科学版)], 2021, 43(3): 20. [本文引用:1]
[12] Ali Z, Mahmood T. Computational and Applied Mathematics, 2020, 39: 161. [本文引用:1]
[13] Yager R R. IEEE Transactions on Systems, 2001, 31(6): 724. [本文引用:1]
[14] Zindani D, Maity S R, Bhowmik S. Soft Computing, 2020, 24(18): 114091. [本文引用:1]
[15] Liu Z, Teng F, Liu P, et al. International Journal for Uncertainty Quantifications, 2018, 8(3): 211. [本文引用:1]
[16] Garg H, Ali Z, Mahmood T. Expert Systems, 2021, 38(1): e12609. [本文引用:2]
[17] Ren W, Du Y. Journal of Natural Fibers, 2023, 20(2): 2219906. [本文引用:3]