2024年大数据最新数学建模之层次分析法_数学建模的层次分析法,大数据开发程序员

img
img

网上学习资料一大堆,但如果学到的知识不成体系,遇到问题时只是浅尝辄止,不再深入研究,那么很难做到真正的技术提升。

需要这份系统化资料的朋友,可以戳这里获取

一个人可以走的很快,但一群人才能走的更远!不论你是正从事IT行业的老鸟或是对IT行业感兴趣的新人,都欢迎加入我们的的圈子(技术交流、学习资源、职场吐槽、大厂内推、面试辅导),让我们一起学习成长!

i

,

j

×

a

j

,

i

=

1

a_{i,j}×a_{j,i}=1

ai,j​×aj,i​=1 (我们称满足这一条件的矩阵为

正互反矩阵

{\color{Red}正互反矩阵 }

正互反矩阵)

❗️ 注意:

  • 在层次分析法中,我们构造的判断矩阵都是正互反矩阵

实际上,上面这个矩阵就是层次分析法中的

判断矩阵

{\color{Red}判断矩阵 }

判断矩阵

2️⃣ 判断矩阵C-P

在这里插入图片描述

❗️注意:

  • 判断矩阵中的元素只能是1至9和它们的倒数,有些文章中填入了其他的数(例如3/2、5/4),这是不可以的。

2.3 一致性检验

在填写判断矩阵的时候,可能会出现一个问题,如下:

在这里插入图片描述

若正互反矩阵满足

a

i

,

j

×

a

j

,

k

=

a

i

,

k

a_{i,j}×a_{j,k}=a_{i,k}

ai,j​×aj,k​=ai,k​,我们称其为

一致矩阵

{\color{Red}一致矩阵 }

一致矩阵

在这里插入图片描述

另外,如果一个正互反矩阵的各行(列)之间呈倍数关系,那么它一定是一致矩阵

❗️注意:

  • 在使用判断矩阵求权重之前,必须对其进行一致性检验
  • 只有非一致矩阵才需要进行一致性检验,如果我们的判断矩阵本身就是一个一致矩阵(明显的各行各列成倍数关系),那么没必要进行一致性检验

对判断矩阵进行

一致性检验

{\color{Red}一致性检验 }

一致性检验,其主要思想是,检验我们构造的判断矩阵和一致矩阵是否有太大的差别。

在这里插入图片描述

在这里插入图片描述

判断矩阵越不一致时,最大特征值与 n 相差就越大

一致性检验的步骤:

  • 计算一致性指标CI

C

I

=

λ

m

a

x

n

n

1

CI=\frac{\lambda_{max}-n }{n-1}

CI=n−1λmax​−n​

  • 查找对应的平均随机一致性指标RIRI我们只需要会查表即可,不用管怎么来的)
    在这里插入图片描述
  • 计算一致性比例CR

C

R

=

C

I

R

I

CR=\frac{CI}{RI}

CR=RICI​
+ 如果CR < 0.1, 则可认为判断矩阵的一致性可以接受,继续计算被比较元素对于该准则的相对权重;
+ 否则,需要对判断矩阵进行修正(往一致矩阵上调整,一致矩阵的各行或列成倍数关系)

2.4 计算权重

上述,我们已经得到了判断矩阵,那么我们如何根据判断矩阵计算权重呢?

三种方法计算权重,强烈建议在比赛时三种方法都使用,在论文中加入:

**以往的论文利用层次分析法解决实际问题时,都是采用其中某一种方法求权重,而不同的计算方法可能会导致结果有所偏差。为了保证结果的

稳健性

{\color{Red} 稳健性}

稳健性 ,本文采用了三种方法分别求出了权重后计算平均值,再根据得到的权重矩阵计算各方案的得分,并进行排序和综合分析,这样避免了采用单一方法所产生的偏差,得出的结论将更全面、更有效。**

2.4.1 算术平均法

在这里插入图片描述

使用数学符号表示如下(可以放到论文中):

在这里插入图片描述

2.4.2 几何平均法

在这里插入图片描述

2.4.3 特征值法

在这里插入图片描述

在这里插入图片描述

权重计算完成后,我们就可以将计算结果填入权重表中:

在这里插入图片描述

在这里插入图片描述

汇总结果得到权重矩阵:

在这里插入图片描述

我们可以得到使用特征值法求得的权重矩阵,根据此矩阵,我们可以计算出每个旅游景点的得分。

2.5 计算各方案的得分

在这里插入图片描述

类似的,我们可以得到北戴河得分为

0.245

0.245

0.245,桂林得分为

0.455

0.455

0.455,因此最佳的旅游景点是桂林。

❗️注意:这里用EXCEL计算可大大减轻工作量,提示F4锁定单元格

3 编写代码

本例中的判断矩阵:

[1, 1/2, 4, 3, 3;
 2, 1,   7, 5, 5;
 1/4, 1/7, 1, 1/2, 1/3;
 1/3, 1/5, 2, 1, 1;
 1/3, 1/5, 3, 1, 1]
也可以写成一行(最后一行后面的分号可要可不要):
[1, 1/2, 4, 3, 3;2, 1,   7, 5, 5;1/4, 1/7, 1, 1/2, 1/3;1/3, 1/5, 2, 1, 1;1/3, 1/5, 3, 1, 1;]


[1,2,5;
1/2,1,2;
1/5,1/2,1;]

[1,1/3,1/8;
3,1,1/3;
8,3,1;]

[1,1,3;
1,1,3;
1/3,1/3,1;]

[1,3,4;
1/3,1,1;
1/4,1,1;]

[1,1,1/4;
 1,1,1/4;
 4,4,1;]

代码:

%% 注意:在论文写作中,应该先对判断矩阵进行一致性检验,然后再计算权重,因为只有判断矩阵通过了一致性检验,其权重才是有意义的。
%% 在下面的代码中,我们先计算了权重,然后再进行了一致性检验,这是为了顺应计算过程,事实上在逻辑上是说不过去的。
%% 因此大家自己写论文中如果用到了层次分析法,一定要先对判断矩阵进行一致性检验。
%% 而且要说明的是,只有非一致矩阵的判断矩阵才需要进行一致性检验。
%% 如果你的判断矩阵本身就是一个一致矩阵,那么就没有必要进行一致性检验。


disp('请输入判断矩阵A')
A=input('A=');
[n,n] = size(A);
% % % % % % % % % % % % %方法1: 算术平均法求权重% % % % % % % % % % % % %
Sum_A = sum(A);
SUM_A = repmat(Sum_A,n,1);
Stand_A = A ./ SUM_A;

disp('算术平均法求权重的结果为:');
disp(sum(Stand_A,2)./n)
% % % % % % % % % % % % %方法2: 几何平均法求权重% % % % % % % % % % % % %
Prduct_A = prod(A,2);
Prduct_n_A = Prduct_A .^ (1/n);
disp('几何平均法求权重的结果为:');
disp(Prduct_n_A ./ sum(Prduct_n_A))
% % % % % % % % % % % % %方法3: 特征值法求权重% % % % % % % % % % % % %
[V,D] = eig(A);
Max_eig = max(max(D));
[r,c]=find(D == Max_eig , 1);
disp('特征值法求权重的结果为:');
disp( V(:,c) ./ sum(V(:,c)) )
% % % % % % % % % % % % %下面是计算一致性比例CR的环节% % % % % % % % % % % % %
CI = (Max_eig - n) / (n-1);
RI=[0 0.0001 0.52 0.89 1.12 1.26 1.36 1.41 1.46 1.49 1.52 1.54 1.56 1.58 1.59];  %注意哦,这里的RI最多支持 n = 15
% 这里n=2时,一定是一致矩阵,所以CI = 0,我们为了避免分母为0,将这里的第二个元素改为了很接近0的正数
CR=CI/RI(n);
disp('一致性指标CI=');disp(CI);
disp('一致性比例CR=');disp(CR);
if CR<0.10
    disp('因为CR<0.10,所以该判断矩阵A的一致性可以接受!');
else
    disp('注意:CR >= 0.10,因此该判断矩阵A需要进行修改!');
end

🌟 代码优化:

  • 如果我们输入的是一个二阶的判断矩阵,请观察结果有什么问题?怎么改进代码来修正这个问题。(提示:二阶判断矩阵一定是一致矩阵
    • 会造成CRNANn=2 时,一定是一致矩阵,所以CI = 0,我们为了避免分母为0,将 RI 第二个元素改为了很接近0的正数

4 模型局限性

  1. 评价的决策层不能太多,太多的话

img
img
img

既有适合小白学习的零基础资料,也有适合3年以上经验的小伙伴深入学习提升的进阶课程,涵盖了95%以上大数据知识点,真正体系化!

由于文件比较多,这里只是将部分目录截图出来,全套包含大厂面经、学习笔记、源码讲义、实战项目、大纲路线、讲解视频,并且后续会持续更新

需要这份系统化资料的朋友,可以戳这里获取

🌟 代码优化:

  • 如果我们输入的是一个二阶的判断矩阵,请观察结果有什么问题?怎么改进代码来修正这个问题。(提示:二阶判断矩阵一定是一致矩阵
    • 会造成CRNANn=2 时,一定是一致矩阵,所以CI = 0,我们为了避免分母为0,将 RI 第二个元素改为了很接近0的正数

4 模型局限性

  1. 评价的决策层不能太多,太多的话

[外链图片转存中…(img-nZumBp2a-1715270650800)]
[外链图片转存中…(img-bvOLB2jh-1715270650800)]
[外链图片转存中…(img-aKOB4ptc-1715270650801)]

既有适合小白学习的零基础资料,也有适合3年以上经验的小伙伴深入学习提升的进阶课程,涵盖了95%以上大数据知识点,真正体系化!

由于文件比较多,这里只是将部分目录截图出来,全套包含大厂面经、学习笔记、源码讲义、实战项目、大纲路线、讲解视频,并且后续会持续更新

需要这份系统化资料的朋友,可以戳这里获取

  • 3
    点赞
  • 10
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值