数学建模———层次分析法及其matlab语法,函数和代码实现

  1. 层次分析法思想登场

建模比赛中最基础的模型之一,其主要用于解决评价类问题(例如:选择那种方案最好,哪位运动员或者员工表现的更优秀。)

评价类问题字眼:

  1. 评价的目标是什么?
  2. 达到这个目标有那几种方案?
  3. 评价准则或指标是什么?

确定权重的方方法——分而治之

(两两比较推算出权重)

比较的标度:

 

例子引用(本文引用的例子及博主跟的课中举的例子)

旅游景点的选择

 正互反矩阵/判断矩阵:
根据标度所填的矩阵有着一定的规律矩阵中坐标[I,j]与坐标[j,i]的两个元素相乘为1且坐标i==j的元素为1.

 一个不成问题的问题:

在填入矩阵的时候我们填入的数值可能会出现矛盾,此时需要做出修改

 

一致矩阵:

特点:各行/各列成倍数关系且一定是正互反矩阵

 

一致性矩阵的定义:

 

一致性矩阵的检验:

 

 

计算CI中λmax 是最大特征值(后面将特征值求权重时会解释),n是方形矩阵的行/列(因为为方形矩阵所以行=列)

RI是课中给出的

计算权重:

(由于计算权重方法的原理文字或口头叙述过于枯燥,晦涩,我们直接告诉方法,后面在代码实现的时候我们会详细讲解。)

 

 

 

  1. 算数平均法求权重

 

 

 

2.几何平均法求权重

  1. 3.特征值法求权重

 

 

                                        

  1. MATLAB基本语法函数以及层次分析法的实现

(博主是计算机专业所以很多语法函数会与c/cpp中的语法做比较方便计算机专业的理解)

  1. Matlab中的每个语句不需要必须加 ‘ ;’ ,若加上分号表示在执行中此语句执行但不显示出来,若不加分号则执行后语句结果会在命令行窗口(c中的控制台)显示

 

(两个语句都执行了却只显示了a=20)

2.注释:matlab中的注释符号是%(类似c中的//)

 

3.matlab中的命令行窗户(c中的控制台)和工作区(c中监视区的监视变量)每次执行后都会有保存,不会清空,所以使用clear(清楚工作区)clc(清除命令行窗口)

clear;clc  平时都会这样两个连用清楚所有:起到了初始化的作用,防止之前的结果对新的脚本产生干扰

4输入函数disp()

()中为要输出的内容,若为字符串需要用“”/’’引起来(类似c/cpp中的printf,cout)

注意:disp后面加不加;都会输出

5.向量(c中数组)一般赋值矩阵使用,为了方便理解我们成为矩阵。

矩阵中间列的符号是“ ,”(英文逗号)“ ”(空格)

间隔行的符号是“;”(分号)

Eg:

 

6.合1并字符串

a.strcat函数   strcat(‘str1’,’str2’)

b. [‘str1’,’str2’](看成一个数组)

7.num2str()将数字转换成字符串(因为直接把数字当字符串输出会根据ASCII码转换成对应的字符如:97对应‘a’)

如果使用disp()函数输出的时候在disp括号中直接把数字转换成字符串且与已有字符串相连接

8.输入函数 intput()(类似与c中的scanf)

9.求和函数sum(),若括号中为矩阵,默认只输入矩阵是按列求和,最后得到一个一列n行的向量。Sum(A,1)后面加1是按列求和sum(A,2)后面加2是按行求和

10.提取矩阵中指定位置元素

a.取指定行和列A(行,列)

b.取某一行A(行,:)加冒号

c.取某一列A(:,列)

d取某写全部元素A([2,5],:)取第二行和第五行

补充matlab中的等差数列 a1 : d : max(max是数列能取到的最大范围)

如果不加d则默认d为1。

e.取全部元素A( :)

11.size()函数求矩阵的行和列

Size(A,1)返回行size(A,2)返回列

12.repmat()函数复制矩阵函数

B=remat(A,2,1)将A矩阵赋值两行A,一列A成为一个新的矩阵B

13.矩阵之间的运算

* /是矩阵的* /。.*  ./是矩阵中对应的每个元素相互*/。

14.开方 ^    a^n意味a的n次方(同样应用于矩阵,若相对矩阵中每个元素开方应用  .^)

15.eig(A)求矩阵的特征值和特征向量的函数

[V,D]=eig(A)  V表示特征向量   D表示特征值(若只有一个接收则只返回特征值)

16.find()函数  返回向量或矩阵中不为0的元素的下标

若后面加数字n则表示返回前n个不为0的元素。

若x为矩阵:

[r,c]=find(x) 

r:位置不为0的元素的行坐标

c:位置不为0的元素的列坐标

17.矩阵中判断大小关系:

>  <  ==  ~=(对应c中的> ,< ,== ,!=)

返回一个logical数组(其中元素为1代表符合条件,为0不符合)

18.判断语句

 

19.for循环(一定要以end结尾)

For i=1:10(等同于c中的for(int i=1;i<=10;i++))

~~

~~

End

20.相乘函数prod()后面加1为按列相乘,加2为按行相乘

21.ones与zeros函数

 层次分析法matlab代码的实现

算数平均法求权重

%填入方形矩阵
A=[1,1,4,1/3,3;
1,1,4,1/3,3;
1/4,1/4,1,1/3,1/2;
3,3,3,1,3;
1/3,1/3,2,1/3,1]
%求出矩阵每列的和 
SUM_A=sum(A)
%求出矩阵的行数
n=size(A,1)
%将求出来的每列的和扩大n行倍,方便归一化(每个元素除以其所在列的和)
SUM_cA=repmat(SUM_A,n,1)
%第一步归一化处理
A_guiyi=A./SUM_cA
%第二部将归一化元素各列相加(按行求和得到一个列向量)
SUM_guiyiA=sum(A_guiyi,2)
LAST_A=SUM_guiyiA./n

几何平均法求权重

%填入方形矩阵
A=[1,1,4,1/3,3;
1,1,4,1/3,3;
1/4,1/4,1,1/3,1/2;
3,3,3,1,3;
1/3,1/3,2,1/3,1]
%第一步将A元素各行相乘得到一个列向量
Prod_A=prod(A,2)
%第二部将新向量中每个元素开n次方(成1/n次方)
Prod_nA=Prod_A.^(1/n)
%将开n次方后的向量归一化处理得到权重(归一化就是将向量中的每个元素除向量和)
disp('几何平均法求权重结果:')
disp(Prod_nA/sum(Prod_nA))

特征值法求权重

%填入方形矩阵
A=[1,1,4,1/3,3;
1,1,4,1/3,3;
1/4,1/4,1,1/3,1/2;
3,3,3,1,3;
1/3,1/3,2,1/3,1]
%第一步求出A的最大特征值,以及其对应的特征向量(V向量,Val值)
[V,val]=eig(A)
%求最大特征值
MAX_val=max(max(val))%第一个max求每列的最大特征值第二个求真个矩阵最大特征值
%求最大特征值对应的特征向量
% 1.找到最大特征值的位置,并保存它的位置
[r,c]=find(val==MAX_val,1)
%2.(***)求出列所在那一行的向量就是最大特征值向量
V(:,c)
%第二步对求出的特征向量进行归一化处理(让最大特征向量的每个元素除其和)得到我们的权重
disp('特征值法求权重的结果为')
disp(V(:,c)./sum(V(:,c)))

博主主要跟着清风数学建模的课程学习,其中里面的一些图片都来源于上课视频的截图。

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.hqwc.cn/news/22941.html

如若内容造成侵权/违法违规/事实不符,请联系编程知识网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

【已解决】JavaScript为控件绑定事件死活不生效

本博文源于最近操作的JavaScript一个诡异的事情。我在用js一个动态生成的按钮&#xff0c;然后想要为其绑定一个alert的事件&#xff0c;结果发现控制台给我报未定义的错误&#xff0c;我慌的不要要的&#xff0c;经过长时间胡乱尝试终于解决了。 问题再现 dataHtml "&…

uniapp 打包安卓apk (原生App)云打包

uniapp 打包安卓apk (原生App)云打包 hbuilder中操作 项目的一些配置appid DCloud appid 用途/作用/使用说明&#xff1a; https://ask.dcloud.net.cn/article/35907 右键我们项目目录-》发行-》原生APP-云打包 说明&#xff1a; 1. 打包安卓&#xff0c;只选择安卓打包项&…

【Redis】—— Redis的RDB持久化机制

&#x1f4a7; 【 R e d i s 】—— R e d i s 的 R D B 持久化机制 \color{#FF1493}{【Redis】 —— Redis的RDB持久化机制} 【Redis】——Redis的RDB持久化机制&#x1f4a7; &#x1f337; 仰望天空&#xff0c;妳我亦是行人.✨ &#x1f984; 个人主页——微风撞…

Jupyter notebook添加与删除kernel

目录 1 添加虚拟环境的kernel 2 删除jupyter notebook已有的kernal 3 切换内核与查看当前内核 4 添加C语言的kernel 5 添加python2的kernel 6 添加java语言的kernel 6.1 sudo apt install default-jre 6.2 下载并安装 ijava 6.3 sudo apt install openjdk-11…

RabbitMQ ---- Hello World

RabbitMQ ---- Hello World 1. 依赖2. 消息生产者3. 信息消费者 本节使用 Java 编写两个程序。发送单个消息的生产者和接收消息并打印出来的消费者。 1. 依赖 <!--指定 jdk 编译版本--><build><plugins><plugin><groupId>org.apache.maven.plu…

前端开发中的微服务架构设计

前端服务化和小程序容器技术为前端应用带来了更好的组织结构、可维护性和可扩展性。这些技术的应用将促进前端开发的创新和发展&#xff0c;使团队能够更好地应对复杂的前端需求和业务挑战。通过将前端视为一个服务化的架构&#xff0c;我们能够构建出更强大、可靠且可持续的前…

独立看门狗 IWDG

独立看门狗介绍 Q&#xff1a;什么是看门狗&#xff1f; A&#xff1a;可以理解为对于一只修勾的定时投喂&#xff0c;如果不给它吃东西就会狂叫&#xff0c;因此可以通过观察修勾的状态来判断喂它的人有没有正常工作。 在由单片机构成的微型计算机系统中&#xff0c;由于单…

Pixi + Tone 实现简单midi音频可视化

依赖库 Pixi.js 是一个前端图形渲染库&#xff0c;使用精灵技术绘制高性能的图形。Tone.js是一个前端音频框架&#xff0c;对web audio api进行了封装&#xff0c;可以快速创建音频样本、音频效果、进行音频分析和音频播放。tonejs/midi是tonejs的一个插件&#xff0c;可以讲m…

解决小程序 scroll-view 里面的image有间距、小程序里面的图片之间有空隙的问题。

1&#xff09;小程序 image跟view标签上下会有间隙&#xff0c;解决方法如下&#xff1a; 在image那里设置vertical-align:top/bottom/text-top/text-bottom 原因&#xff1a;图片文字等inline元素默许是跟父级元素的baseline对齐&#xff0c;而baseline又和父级底边有必定间距…

http1.0、http1.1 http 2.0

HTTP/1.0是无状态、无连接的应用层协议。 无连接 无连接&#xff1a;每次请求都要建立连接&#xff0c;需要使用 keep-alive 参数建立长连接、HTTP1.1默认长连接keep-alive   无法复用连接&#xff0c;每次发送请求都要进行TCP连接&#xff0c;TCP的连接释放都比较费事&…

【测试开发】测试用例的设计方法

目录 一. 测试用例的基本要素 二. 测试用例的设计方法 1. 测试用例设计的万能公式 水杯测试用例 2. 基于需求的设计方法 邮箱注册测试用例 3. 等价类方法 有效等价类和无效等价类 等价类思想设计测试用例步骤 4. 边界值方法 边界值思想设计测试用例步骤 5. 判定表方法…

120、仿真-51单片机温湿度光照强度C02 LCD1602 报警设计(Proteus仿真+程序+元器件清单等)

方案选择 单片机的选择 方案一&#xff1a;STM32系列单片机控制&#xff0c;该型号单片机为LQFP44封装&#xff0c;内部资源足够用于本次设计。STM32F103系列芯片最高工作频率可达72MHZ&#xff0c;在存储器的01等等待周期仿真时可达到1.25Mip/MHZ(Dhrystone2.1)。内部128k字节…