AHP层次分析法Python代码:让AHP帮你选一个最符合你的礼物

1,241 阅读5分钟

本文已参与「新人创作礼」活动,一起开启掘金创作之路。

前言

送什么礼物才能让女人满意,男人苦不堪言。像我这种有选择恐惧症的,每当节日来临一堆东西摆在面前都不知道挑啥(bushi,确实我还没有女友,落泪)。干脆就整个AHP帮咱挑选一个适合送给对象或者自己(也要爱自己喔~)的礼物,AHP原理其实都不需要理解,咱们会实际运用就贴切了,想要理解更多的欢迎看看我的博客有详解喔:层次分析法(AHP)原理以及应用。废话不多说了咱们开始吧!

一、建立层次结构模型

1.确定目标层

第一步我们需要确定我们的目标是什么,是选择个高端大气上档次的礼物呢,还是选择个低调奢华有内涵的礼物。 确定了我们需要挑选啥礼物的目标,会影响后续挑选礼物品质的打分。所以我们需要确定要挑选哪种品质的礼物。像我这种纯纯码农来说一般喜欢那种不贵且看上去很有内涵也上档次的礼物(是不是要求太高了)。当然大家可以自己决定想要什么,这里我先根据本人的目标来,大家看会一遍后自己就能做出来了:

image.png

 2.确定准则层

第二步我们需要思考一些礼物的属性,也就是我们一般衡量礼物的价值要考虑到的东西。比如贵不贵啊、颜值咋样、精准小巧、实用性高啊等等。这些将属性影响到我们最终选择什么礼物,当然就是建立在各个因素对比之上。首先我们确定要选择礼物的考虑品质:

image.png 当然大家也可以自行选择不同的考虑因素。

3.确定方案层

第三步也就是供我们选择具体送的哪种类型的礼物,如巧克力、鲜花、口红、化妆品、游戏(、手表包包等等。

image.png 这几个礼物应该是目前最热门的吧,现在鲜花都不能当面送了有点可惜,当然大家可以根据自己想送的礼物来衡量。

这样一来我们就建立了层次模型了:

image.png

二、构造判断(成对比较)矩阵

1.构建对比矩阵 咱们提出了考虑因素和礼物当然要进行对比了,这里我们不是把所有因素加起来一起比较,而是两两进行比较:

不把所有因素放在一起比较,而是两两相互比较。 对此时采用相对尺度,以尽可能减少性质不同的诸因素相互比较的困难,以提高准确度。 这里我们需要用到两两因素对比之间衡量二者重要性的标度:

image.png

现在我们要根据考虑因素来衡量一些各个因素之间到底哪个对于我们最终选择的目标最值得考虑:

image.png 这里需要自己主观对比打分。

2.获取权向量

接下来我们需要通过列向量归一化和行向量归一化获取权向量,过程其实很简单并不复杂:

(1).列向量归一化

我们依据矩阵来看:矩阵为ij列,这第一行第一列就是a_{11},那么列向量归一化运算就是把第一行第一个元素进行:a_{11}=\frac{a_{11}}{a_{11}+a_{21}+a_{31}+...a_{i1}}. 根据此运算我们把上述矩阵进行列向量归一化:

import numpy as np
a=np.array([1,4,2,0.5,1/3])
b=np.array([1/4,1,0.2,0.5,0.2])
c=np.array([1/2,5,1,1,0.5])
d=np.array([2,2,1,1,0.2])
e=np.array([3,6,2,5,1])
a/a.sum()
b/b.sum()
c/c.sum()
d/d.sum()
e/e.sum()

image.png

(2).行和归一化

行和归一化为每行的数相加除以每行的个数,这里进行降维转为i行1列的矩阵,计算公式为:

image.png

np.sum(ep,axis=1)/5
 
array([0.16109798, 0.45525528, 0.15045596, 0.17535918, 0.05783161])

这便得到了我们的权向量。

3.计算特征值

我们假设特征值为\lambda,则线性代数特征公式为A\omega =\lambda \omega,其中\omega为权向量。

image.png

b=np.array([1/4,1,0.2,0.5,0.2])
c=np.array([1/2,5,1,1,0.5])
d=np.array([2,2,1,1,0.2])
e=np.array([3,6,2,5,1])
em=np.array([a,b,c,d,e])
a=a/a.sum()
b=b/b.sum()
c=c/c.sum()
d=d/d.sum()
e=e/e.sum()
ep=np.array([a,b,c,d,e]).T
ex=np.sum(ep,axis=1)/5
lamda=np.matmul(ex,em)/ex
lamda.sum()*1/5

得到最大特征根:

三、一致性检验

现在我们需要巩固计算出来的成果,是否具有科学性和可靠性。这时候我们需要进行一致性检验,验证我们的模型是否具有说服力。

一致性检验原理: 检验我们构造的判断矩阵和一致矩阵是否有太大的差别。 1.第一步计算CI 我们需要计算衡量一致性的指标:

我们需要计算衡量一致性的指标CICI=\frac{\lambda _{max}-n}{n-1}

  • CI=0,有完全的一致性;
  • CI接近于0,有满意的一致性;
  • CI越大,不一致越严重 为了衡量CI的大小,引入随机一致性指标RI

2.第二步计RI

RI的值根据n的大小来决定:

image.png

 3.第三步计算CR

一般,当一致性比率CR<0.1时,认为A的不一致程度在容许范围之内,有满意的一致性,通过一致性检验。可用其归一化特征向量作为权向量,否则要重新构造成对比较矩阵A,对a_{ij}加以调整。

image.png 这里我们的判断矩阵计算结果为:

image.png CR=CI/RI=0.0996<0.1,通过一次性检验。

四、填充权重矩阵

接下来我们要根据三种选择方案:口红、巧克力和手工DIY礼品来进行逐个影响因素对比。计算方法是和计算准则层方法一样的,这里便不再逐步进行计算:

image.png

image.png

image.png

image.png

image.png

image.png 接下来大家可自行对比打分,这里不逐个演示了,直接进行填充权重矩阵:

image.png 最终计算各个方案的得分: 口红:

image.png 巧克力:0.2164 DIY:0.2848


点关注,防走丢,如有纰漏之处,请留言指教,非常感谢

以上就是本期全部内容。我是fanstuck ,有问题大家随时留言讨论 ,我们下期见。