数学建模——TOPSIS法

news2024/12/26 22:04:27

TOPSIS法（Technique for Order Preference by Similarity to ldeal Solution）

可翻译为逼近理想解排序法，国内常简称为优劣解距离法

TOPSIS法是一种常用的综合评价方法，其能充分利用原始数据的信息，其结果能精确的反应各评价方案之间的差距。

（上期我们刚讲述了层次分析法，但是对于层次分析法也有一定的限制）

层次分析法的局限性：

那TOPSIS法如何评分呢？

我们的成绩有时是分数越高对于我们越好，但当统计其它事物时，如情商（我们通过与人正常次数来做判断，争吵次数越多情商越低）这样就会出现了我们评价指标有时越大越好有时越小越好对此我们作出定义：

但每次的指标步同我们的统计的方法也会随之改变所以我们定义了将所有的指标都转换为极大型指标方便我们统计：

也就是我们TOPSOSIS法的第一步: 指标正向化

这里将数据正向化分为极小型正向化，中间型正向化，区间型正向化

（由于内容过多防止读者理解混乱我们先介绍大体的TOPSIS的步骤，最后再去讲正向化细节）

第二步：对正向化后的矩阵进行标准化

为了消去不同指标量纲的影响，需要对已经正向化的矩阵进行标准化处理。

标准化处理的解释就是：每一个元素/[（其所在列的元素的平方和）再开根号]

设X（n行，m列）为我们需要处理的数据矩阵。(n为评价目标的个数，m为评价指标的个数)

（接下来直接上matlab代码，为了大家理解的更清晰，将一行一行代码拆成小部分来解析）

  1. 求其所在列的平方和 sum(X.*X).
  2. 开根号sum(X.*X).^0.5

[a.b完成后我们得到了一行m列的向量。为了让每个元素都标准化（/每列平方和的根号）我们需要将的得到的这个向量复制成和X大小一样的矩阵，方便计算。]

c.复制矩阵Y=repmat(sum(X.*X).^0.5,n,1)将得到的结果复制n行使其成为与X相同形状。

d.标准化矩阵Z=X./Y;

 

第三步计算得分并归一化

构造评分方式，方便结果更加符合实际情况

文字描绘：

计算得分的公式（也是我们在写代码时实现的方法）

Z即是我们第二部求出来的标准化方阵

为了方便我们将z与最大值距离定义为D_P,最小值距离定义为D-N

D_P= sum([(Z-repmat(max(Z),n,1)) .^2],2) .^0.5;

Z与max的d：[（每个元素-所在列的max）的平方得到的矩阵按行求和]再开根号

D_P= sum([(Z-repmat(min(Z),n,1)] .^2],2) .^0.5;

Z与min的d：[（每个元素-所在列的min）的平方得到的矩阵按行求和]再开根号

得到D_P&&D_N后我们可以计算对评价对象未归一化的得分：

S=D_N/(D_P+D_N);

最后归一化S得到最后的评分：（S是一个n行1列的向量，对n个评价对象的评分）

Stand_S=S/sum(S)；

为了便于我们的观察我们最后通过sort函数将评分排序

(由于是分越高越好所以我们逆序排序（由大到小）所以矩阵后面要加上‘descend‘)

[sotred_s,index]=sort(stand_S,’descend)

(当sort函数有两个接受值的时候，第一个接收值是排序后的数列，第二个接收值是排序好的数列在原数列的下标，这样方便我们直接找到原数据中那个是第一名)

                                                                                

各种评价类型正向化的方法：

最常见的四种指标：

  1. 极小型指标=>正向化

矩阵中的最大值-每个元素形成的矩阵

  1. 中间型指标=>正向化

根据题目或者我们给定一个中间值：best

定义M：|每个元素-best|的绝对值，其中的最大值

1-（|每个元素-best的绝对值|/M）后得到的矩阵

  1. 区间型指标=>正向化

题目给出或我们天给定最佳区间的区间下限a,和区间上限b

定义M=max{（a-最小元素），（最大元素-b）}

每个元素：

  1. 如果 x<a , xi = 1 - (a-x) / M;
  2. 如果 a<=x<=b , xi = 1;
  3. 如果 x>b , xi = 1 -（x-b）/ M;

后的得到的矩阵。

                                                                            

列题以及代码实现和讲解：


题目分析：本题让我们根据给出的四个之变让我们对A-T这20个对象评分。

但由于指标2，3，4都不是极大型指标所以我们需要正向化。

代码实现：

  1. 我们需要将数据导入到matlab中：我们只需要将Excle中的B2：E21导入到matlab中

操作步骤：我们在matlab的工作区右击点新建形成有一个新文件我们点开后会出现一个跟exlce一样的表格我们点复制即可（别忘了重命名）

  1. 导入我们的数据，求其行和列
  2. 根据我们数据指标的类型判断是否需要正向化以及那些列需要正向化

这里我们正向化需要自定义函数。

Matlab中自定义函数不能和主函数在一个脚本中，需要新建一个脚本专门来写自定义函数，且自定义函数的脚本必须和主函数保存在一个文件夹中。

Matlab自定义函注意事项：和c一样主函数中写接受变量=函数名（参数）

但在自定义函数的脚本中所有函数都已function开头相当我我们c中的函数类型int void char等matlab中只有function，且所有函数都要以end结尾（相当于我们c中的break但是end不返回变量）

格式：function [返回变量] =函数名 （参数）

edn

3.对选定的列进项正向化处理

4.标准化处理

5.计算得分并归一化

代码参考：

主函数：

clear;clc
load data_water_quality.mat%导入的数据如果我们要引用的话记得前面添加load，但前提是你已经保存在和此代码同一文件夹中


%第一步判断指标是否需要正向化
[n,m] = size(X);%我们先获取我么导入数据的行行和列（n，m）
%对于有些指标不是极大型指标所以我们需要对其正向化
disp(['共有' num2str(n) '个评价对象' num2str(m) '个评价指标'])
judge = input(['这' num2str(m) '个指标是否需要正向化处理，需要请输入1，不需要请输出0：']);
%我们按评价指标的列排序做序号输入要正向化的指标对应的序号
Position = input('请输需要正向化处理的指标的所在的列，列如第1，2，3三列需要处理，那么你需要输入[1,2,3]: ');%[2,3,4]
%输入指标后我们position保存了要修改的指标序号
disp('请输入需要处理的这些列的指标类型（1：极小型，2：中间型，3：区间型）')
%我们type函数保存了对应指标的指标类型（如positon（1）是中间值那么type（i）就是3）
Type=input('列如：第二列是极小型，第二列是区间型，第三列是中间型，就输入 [1,3,2]:');%[2,1,3]
%我们自定义函数传参然后正向化
for i=1 : size(Position,2)
    %X(:,Position(i))就是对我们X矩阵的列进行传参和接收
    X(:,Position(i))=Positivization(X(:,Position(i)),Type(i),Position(i));
end
disp('正向化后的矩阵X=')
disp(X)
%end
%第二步:对正向化后的矩阵进行标准化
Z = X./repmat(sum(X.*X).^0.5,n,1)
disp('标准化矩阵Z = ')
disp(Z)
%第三步计算最大值的距离和最小值的距离，并算出得分
D_P = sum([(Z-repmat(max(Z),n,1)).^2],2).^0.5;
D_N = sum([(Z-repmat(min(Z),n,1)).^2],2).^0.5;
S = D_N ./ (D_P + D_N);
disp('最后得分为：')
stand_S = S/sum(S);
[sorted_S,index] = sort(stand_S,'descend')

正向化函数： 

function [posit_x]=Positivization(x,type,i)
if type==1%及小型
    disp(['第' num2str(i) '列是极小型，正在正向化'])
    posit_x = Min2Max(x);
    disp(['第' num2str(i) '列极小型正向化处理完成'])
    disp('~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~')
elseif type == 2%中间型
    disp(['第' num2str(i) '列是中间型'])
    best = input('请输入最佳的那个值：')
    posit_x = Mid2Max(x,best);
    disp(['第' num2str(i) '列中间型正向化处理完成'])
    disp('~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~')
elseif type == 3%区间型
    disp(['第' num2str(i) '列是中间型，正在正向化'])
    a = input('请输入区间的下界:');
    b = input('请输入区间的上界:');
    posit_x = Inter2Max(x,a,b);
    disp(['第' num2str(i) '列中间型正向化处理完成'])
    disp('~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~')
else
    disp('没有这种类型的指标，请检查Type向量中是否除了1，2，3之外的其他值')
end

极小型正向化函数：

%极小型正向化：公式 max-每个元素
function [posit_x] = Min2Max(x)
posit_x=max(x)-x;
end

中间型正向化函数：

%中间型正向化
%公式：先求M=max{|最大值-xi|的绝对值}，xi=1-|最大值-x|/M
function [posit_x]  = Mid2Max(x,best)
M = max(abs(x-best));
posit_x = 1 - abs(x-best) / M;
end

区间型正向化函数：

%区间型正向化
%公式：
%先求M=max{（a-最小值），（最大值-b）}
%xi:if xi<a,xi=1-(a-x)/M
  %if xi>b,xi=1-(b-x)/M
  %if a<=xi<=b,xi=1
function [posit_x] = Inter2Max(x,a,b)
r_x = size(x,1);
M = max([a-min(x),max(x)-b]);
posit_x =zeros(r_x,1);
for i=1 : r_x
    if x(i) < a
        posit_x(i) = 1-(a-x(i))/M;
    elseif x(i)>b
        posit_x(i) = 1-(x(i)-b)/M;
    else
        posit_x(i) = 1;
    end
  end
end

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

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

如若内容造成侵权/违法违规/事实不符，请联系多彩编程网进行投诉反馈，一经查实，立即删除！

相关文章

一款基于JAVA开发的Echarts后台生成框架

一款基于JAVA开发的Echarts后台生成框架

目录 前言 一、Echarts后台生成框架 1、简介 2、开源地址 3、支持类型 二、图表生成实战 1、项目使用 2、后台折线图生成 3、Echarts改造 4、最终效果 总结 前言 之前的博文主要分享了一些关于Echarts的基本开发知识&#xff0c;在之前的博客中主要介绍的是在前端页面进…
阅读更多...
魏副业而战：她又办了2套图书证

魏副业而战：她又办了2套图书证

我是魏哥&#xff0c;与其躺平&#xff0c;不如魏副业而战&#xff01; 社群成员董姐又办了2套图书证&#xff0c;加上之前的1套&#xff0c;她已经有3套图书证了。 3套图书证意味着可以开9个闲鱼图书店铺。是不是感觉很NB。 董姐说&#xff0c;2套图书证等了快一个多月&…
阅读更多...
Flink基本原理剖析讲解

Flink基本原理剖析讲解

1.Flink简介 Flink是一个批处理和流处理结合的统一计算框架&#xff0c;其核心是一个提供了数据分发以及并行化计算的流数据处理引擎。它的最大亮点是流处理&#xff0c;是业界最顶级的开源流处理引擎。 Flink最适合的应用场景是低时延的数据处理&#xff08;Data Processing…
阅读更多...
面试题更新之-伪元素和伪类

面试题更新之-伪元素和伪类

文章目录 伪元素和伪类是什么&#xff1f;伪元素&#xff08;Pseudo-elements&#xff09;:伪类&#xff08;Pseudo-classes&#xff09;: css伪元素和伪类的区别使用css伪元素的好处使用css伪类的好处 伪元素和伪类是什么&#xff1f; 在CSS中&#xff0c;伪元素&#xff08;…
阅读更多...
星云零售信贷基于 Apache Doris 的 OLAP 演进之路

星云零售信贷基于 Apache Doris 的 OLAP 演进之路

本文导读&#xff1a; 腾梭科技是国内领先的零售金融数字化及安全服务提供商&#xff0c;是腾讯投资且在金融领域的战略合作伙伴&#xff0c;并与腾讯联合研发了“星云智慧信贷解决方案。在其信贷业务转型过程中&#xff0c;随着系统规模不断扩大&#xff0c;早期架构无法再满…
阅读更多...
centos7根分区、文件系统扩容

centos7根分区、文件系统扩容

1、 输入lsblk&#xff0c;查看到新硬盘sde&#xff0c;根目录现71G. 2、 创建分区fidisk /dev/sde 3、 刷新分区 partprobe /dev/sde&#xff0c;并创建物理卷 pvcreate /dev/sde1 4、 查看卷组名 vgdisplay 5、 将物理卷扩展到卷组 vgextend centos /dev/sde1 6、 查看逻辑巻…
阅读更多...
linux系统安装步骤教程详解

linux系统安装步骤教程详解

linux是现在常用的一类操作系统&#xff0c;我们可以通过U盘、光驱、或者虚拟机进行安装。 具体步骤如下&#xff1a; 1.去https://www.centos.org/download/下载centos系统&#xff0c;其中CentOS-7.0-x86_64-DVD-1503-01.iso 是标准安装版&#xff0c;推荐下载这个; 2.使用U盘…
阅读更多...
OCR技术的昨天今天和明天！2023年最全OCR技术指南！

OCR技术的昨天今天和明天！2023年最全OCR技术指南！

OCR是一项科技革新&#xff0c;通过自动化大幅减少人工录入的过程&#xff0c;帮助用户从图像或扫描文档中提取文字&#xff0c;并将这些文字转换为计算机可读格式。这一功能在许多需要进一步处理数据的场景中&#xff0c;如身份验证、费用管理、自动报销、业务办理等都显得尤为…
阅读更多...
MySQL常见笔试题

MySQL常见笔试题

前言 数据库的考察在面试时可是十分常见的&#xff0c;MySQL作为一种常用的关系型数据库管理系统&#xff0c;对于它的介绍在面试时可是必不可少的&#xff0c;下面就是一些常见笔试题的模拟&#xff0c;希望可以帮助到你&#x1f642;&#x1f642; 所用到的表如下&#xff…
阅读更多...
国产MCU-CW32F030开发学习-BH1750模块

国产MCU-CW32F030开发学习-BH1750模块

国产MCU-CW32F030开发学习-BH1750模块 硬件平台 CW32_48F大学计划板CW32_IOT_EVA物联网开发评估套件BH1750数字型光照强度传感器 BH1750 BH1750是一款数字型光照强度传感器&#xff0c;能够获取周围环境的光照强度。其测量范围在0~65535 lx。lx勒克斯&#xff0c;是光照强…
阅读更多...
界面设计利器！10款Sketch插件合集，助你事半功倍！

界面设计利器！10款Sketch插件合集，助你事半功倍！

在界面设计领域&#xff0c;Sketch以其高效、小巧的优势获得了不少设计团队的喜爱&#xff0c;帮助全球设计师创造了许多不可思议的作品。在使用Sketch的过程中&#xff0c;辅助使用一些Sketch插件&#xff0c;可以让我们更加高效地完成设计任务。本篇文章&#xff0c;我们将揭…
阅读更多...
torch分布式通信基础

torch分布式通信基础

torch分布式通信基础 1. 点到点通信2. 集群通信 官网文档&#xff1a;WRITING DISTRIBUTED APPLICATIONS WITH PYTORCH 1. 点到点通信 # 同步&#xff0c;peer-2-peer数据传递 import os import torch import torch.distributed as dist import torch.multiprocessing as mpdef…
阅读更多...
热点探测技术架构设计与实践

热点探测技术架构设计与实践

1. 概述 说到热点问题&#xff0c;首先我们先理解一下什么是热点&#xff1f; 热点通常意义来说&#xff0c;是指在一段时间内&#xff0c;被广泛关注的物品或事件&#xff0c;例如微博热搜&#xff0c;热卖商品&#xff0c;热点新闻&#xff0c;明星直播等等&#xff0c;所以…
阅读更多...
领域知识图谱的医生推荐系统：利用BERT+CRF+BiLSTM的医疗实体识别，建立医学知识图谱，建立知识问答系统

领域知识图谱的医生推荐系统：利用BERT+CRF+BiLSTM的医疗实体识别，建立医学知识图谱，建立知识问答系统

项目设计集合&#xff08;人工智能方向&#xff09;&#xff1a;助力新人快速实战掌握技能、自主完成项目设计升级&#xff0c;提升自身的硬实力&#xff08;不仅限NLP、知识图谱、计算机视觉等领域&#xff09;&#xff1a;汇总有意义的项目设计集合&#xff0c;助力新人快速实…
阅读更多...
Xcode报错--访问keychain，出现弹窗处理方案

Xcode报错--访问keychain，出现弹窗处理方案

情景 访问keychain弹出弹窗&#xff0c;不想人工点击&#xff0c;比如自动化测试中使用keychain中的证书的情况 原因 Mac的保护机制 处理 1、人工&#xff1a;输入Password&#xff0c;点击Allow或者Always Allow 2、命令行处理 security unlock-keychain -p "<…
阅读更多...
Spring @RequestMapping 工作原理

Spring @RequestMapping 工作原理

Spring RequestMapping 工作原理 配置基础启动类及Controller类 SpringBootApplication public class DemoServiceApplication {public static void main(String[] args) {SpringApplication.run(DemoServiceApplication.class, args);} }RestController public class HelloC…
阅读更多...
列表定义状态比较不错的UI写法

列表定义状态比较不错的UI写法

<el-table-columnprop"status"label"状态"align"left":formatter"formatTd" ><template slot-scope"scope"><span class"grayStatus" v-if"scope.row.status 1">未开始</span>…
阅读更多...
Linux安装最新版的gcc13.1.0编译器，支持c++20、23

Linux安装最新版的gcc13.1.0编译器，支持c++20、23

Linux安装最新版的gcc13.1.0编译器&#xff0c;支持c20、23 最近在写c20的代码&#xff0c;所以需要升级支持c20及23的编译器&#xff0c;貌似gcc11就已经支持了c20了&#xff0c;但是我这里选择了最新的13.1版本。本文全程实操&#xff0c;上机验证通过。 查看gcc版本 gcc -v…
阅读更多...
MySql 高级-0711

MySql 高级-0711

3. 查询截取分析 分析 分析&#xff0c;至少跑一天&#xff0c;看看生产的慢 SQL 情况开启慢查询日志&#xff0c;设置阙值&#xff0c;比如超过5秒钟的就是慢SQL&#xff0c;并将它抓取出来。explain慢SQL分析Show Profile运维经理 or DBA 进行SQL 数据库服务器的参数调优 总…
阅读更多...
QT-QRegExp和QRegularExpression

QT-QRegExp和QRegularExpression

1.QRegExp qt5.0版本之前正则表示示类是QRegExp,通过它能够筛选出我们想要的数据,它的构造函数如下所示: QRegExp::QRegExp(const QString &pattern, Qt::CaseSensitivity cs Qt::CaseSensitive, QRegExp::PatternSyntax syntax); 其中QRegExp::PatternSyntax syntax用…
阅读更多...
推荐文章
最新文章

Copyright @ 2022~2023