数学建模灵敏度分析_数学建模中的灵敏度分析,到底在分析什么?

798f99d8c28409c2dc04e1e40dddfcd3.png

—2019/6/13—

概念

在数学建模中使用最优化方法时,我们常常会忽略对模型进行灵敏度分析,若缺少这一个步骤,会使得模型的可靠度受到质疑,那到底什么是灵敏度分析?

参考百度:灵敏度分析,是研究与分析一个系统(或模型)的状态或输出变化对系统参数或周围条件变化的敏感程度的方法。

说人话,就是改变模型(公式)的某个参数,引起这个模型输出的变化的程度

举例说明

一头牛重200斤,每天增重5斤,饲料每天花费45元。牛的市场价是每斤65元,但目前行情不好,每天下降1元,求出这头牛的最佳销售时间。

这是一个很简单的单变量优化问题。求这头牛的最佳销售时间,就是求卖出价格最高的时间,为了求出这个时间,我们假设:

销售时间

牛的重量

牛的价格

饲料花费

卖出牛的收益

净收益

于是有,

,其中,
,即得到
净收益模型如下:

…………………………………(1)

容易得出当

时,
取最大值13320元。

其实,大多数情况下,问题到这里就已划上句号。但现实生活中,稍微细心一点,我们会发现,在上述模型的参数中,牛的重量,当前市场价格、每天饲料钱3个参数都是很容易测量的,即确定性较大,但牛的生长率和市场价格的下降率2个参数则不是那么好确定,例子中虽然规定了市场价格下降率

,但实际中,每天的
都是不一样的。于是我们会有一个疑惑——
是不是对于所有的市场价格下降率,这个模型都适用?

为了消除疑惑,我们分别随便选取原下降率(1元/天)周围的几个数

来进行研究。

根据公式(1),分别得出

,得下图:

380c4887d60f6178a0c779be2dbdc835.png

从图中可以看出,牛最佳卖出时间

随着市场价格下降率
的变化而变化,这时候我们称:
牛的卖出最佳时间
对市场价格下降率
是敏感的!

但是究竟有多敏感,就要进行灵敏度分析敲黑板!!!划重点!!!)。

灵敏度分析,就是将市场价格下降率

作为未知参数,计算最佳卖出时间
关系
。这时候,将例子中的“每天下降1元”改成“每天下降
元”,净收益模型就变成:

求导,

的点为

………………………………………(2)

至此,我们已经得到了一种灵敏性关系——

关系,绘制
关系图如下:

26a01a101bc49db788c02d841e247586.png

公式(2)和上图的意义在于——作为牛贩子,我能确定一个

关系,确保我在不同
下,都能确定最佳卖牛时间
。当
时,最佳卖出时间就由公式(2)给出;当
时,
,这会导致
是一个始终下降的曲线,这意味着,最佳卖牛时间为
,也就是说,一旦市场价格下降率
,需要当天就把牛卖掉,这时候饲养牛已经没有任何经济价值了。

然而,在实际使用中,我们更多地是将这种灵敏性关系表示成相对改变量或百分比改变的形式。例如,

的10%的下降,导致了
的38%的增加。如果
的改变量为
,则
的相对改变量为
,百分比改变量为
,如果
改变了
,导致
的改变量,则相对改变量的比值为
的比值,令
,由导数定义,有

业内,称极限值

就是
灵敏性,记为

灵敏性有啥用?

比如,在这个例子中,有

在点

也就是说,在点

=
若市场价格下降率
增加 1%,则最佳卖牛时间
下降3.5%。

结语

本文仅分析了市场价格下降率

和最佳卖牛时间
的灵敏性关系,但可能还存在其他的关系,比如牛的重量增长率和最佳卖牛时间
的关系等,灵敏性分析的成功需要分析者具有较好的判断力,通常我们不需要对模型中的每个参数都进行灵敏性分析,只需要选择那些有
较大不确定性的参数进行灵敏性分析即可。

参考文献

数学建模方法与分析(原书第4版)/(美)Meerschaert, M. M. 著;刘来福,黄海洋,杨淳译.—北京:机械工业出版社, 2014.12:1-9.


—2019/6/18—更新

上述主要是对以下2个问题进行了说明:

  1. 一个模型为什么要进行灵敏度分析?
  2. 如何进行 灵敏度 分析 ?

这次更新目的是说明第3个问题:如何判断灵敏度分析的效果?

先摆出结果:

如果你改变了系统参数后,引起这个模型(公式)输出的变化的程度不大,则说明你的模型稳定性较强(即灵敏性较差),反之则反!

用上述例子说明,

3a25e9a7a98932236c3d1edf43dbc75d.png

改变了1%,
仅下降3.5%,这个改变很小,这就可以说明模型较为稳定;若
改变了1%,
下降了35%,那这个模型绝逼是不稳定的!

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

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

相关文章

虚函数表存放在哪

一个类只能有一个虚函数表。在编译时,一个类的虚函数表就确定了,所以虚函数表在 .rodata ( Linux g )

mysql数据对象

学习目标:了解掌握常见的几种数据库对象学会如何创建具体的数据对象mysql 常见的数据对象有哪些:DataBase/SchemaTableIndexView/Trigger/Function/Procedure多Database用途:业务的隔离资源的隔离表上的常用数据对象:索引约束视图&#xff0c…

赛迪研究院发布《2019量子计算发展白皮书》

来源:赛迪智库经国务院正式批准,由湖南省人民政府、工业和信息化部联合主办,中国电子信息产业发展研究院(简称“赛迪研究院”)、湖南省工业和信息化厅、长沙市人民政府承办的“2019世界计算机大会”在长沙召开。赛迪智…

安卓虚拟摄像头_iPhone 的「第四颗摄像头」位置,为什么给了激光雷达?

一个,两个,三个,四个…… 旗舰手机摄像头从 2016 年开始,以大约每年 1 个摄像头的数量在增长,有人调侃说麻将「八筒」将是所有智能手机摄像头的终极归宿。今年全新一代 iPhone 发布前,不少人推测 iPhone 12…

Mysql复制-Slave库设置复制延迟

mysql> stop slave; mysql> change master to master_delay10;#单位是秒 mysql> start slave; mysql> show slave status\G *************************** 1. row ***************************Slave_IO_State: Waiting for master to send event...SQL_Delay: 10SQL…

c++ vector用另一个vector初始化

vector<int> a(nums.begin(), nums.end()); 注意&#xff1a;vector<int> a(nums.begin()1, nums.end()-2);中a包含nums.begin()1&#xff0c;不包含nums.end()-2&#xff0c;也就是包含的区间是左闭右开区间 如果a不能在声明的时候初始化&#xff0c;可以先定义…

关于边缘计算和边云协同,看这一篇就够了~

来源&#xff1a;Deloitte物联网智库 编译导 读几年前&#xff0c;大多数人都期望将物联网部署至云端&#xff0c;这的确可以给个人用户带来便捷的使用体验&#xff0c;但构建企业级的物联网解决方案&#xff0c;仍然需要采用云计算和边缘计算的结合方案。与纯粹的云端解决方案…

删除当前文件夹下特定名称文件

以下命令删除当前文件夹下所有名称带log的文件 find . -name "*log*" |xargs rm -rf

ospf hello时间和dead_网络工程师_思科 | OSPF由简到难,配合命令学

交换机工作原理&#xff1a;基于源mac学习&#xff0c;基于目的mac转发。路由器的工作原理&#xff1a;收到一个数据包&#xff0c;拆掉帧头&#xff0c;拆开ip包头&#xff0c;提取目的ip地址&#xff0c;查找自己的路由表&#xff0c;有路由&#xff0c;转发&#xff0c;没有…

《自然》深度:中国AI能在2030年制霸世界吗?

来源&#xff1a;智东西,nature中国不仅拥有世界上最多的人口&#xff0c;而且即将成为最大的经济体&#xff0c;并且在人工智能&#xff08;AI&#xff09;方面也有着领先世界的野心。2017年&#xff0c;中国政府制定了《新一代人工智能发展规划》&#xff0c;对人工智能发展进…

object-c 代理反向传值

A.h A.m B.h B.m 转载于:https://www.cnblogs.com/ChouDanDan/p/5692578.html

linux高亮查找关键字

日志文件中通常需要高亮查找某个关键字&#xff0c;方法为 tail 日志文件 | grep 关键字 --colorauto

群晖pxe安装windows_使用win-server2016 进行 pxe 网络装机

首先,安装 windows server 2016此处省略若干字..........安装完成后,将网络设置成桥接模式,并为物理机网卡配置一个静态ip为server16 设置静态ip(要和物理机处于同一网段)c打开开始菜单-->服务器管理器(默认会弹出)2.选择添加角色3.一直点下一步,选择安装DHCP和Windows部署服…

CNCC技术论坛丨联邦学习冲刺人工智能“最后一公里”!

来源&#xff1a;中国计算机学会本论坛将于 CNCC2019 中国计算机大会第一天&#xff08;10月17日&#xff09;在苏州金鸡湖国际会议中心 A305 会议室举行&#xff0c;探讨号称人工智能“最后一公里”的联邦学习&#xff0c;届时业内专家将会碰撞出什么样的火花呢&#xff1f;敬…

javascript语法速查表

转载于:https://www.cnblogs.com/diantao/p/5693877.html

git删除分支

git branch -d 分支名

python和pytorch是什么关系_【PyTorch】Tensor和tensor的区别

本文列举的框架源码基于PyTorch1.0&#xff0c;交互语句在0.4.1上测试通过import torch在PyTorch中&#xff0c;Tensor和tensor都能用于生成新的张量&#xff1a;>>> atorch.Tensor([1,2])>>> atensor([1., 2.])>>> atorch.tensor([1,2])>>&g…

中国芯片将靠此超车!RISC-V架构神在哪全解构

来源&#xff1a; 智东西RISC-V&#xff08;发音为“risk-five”&#xff09;是一个基于精简指令集&#xff08;RISC&#xff09;原则的开源指令集架构&#xff08;ISA&#xff09;&#xff0c;简易解释为开源软件运动相对应的一种“开源硬件”。该项目2010年始于加州大学伯克利…

git撤销分支上代码的修改

新建了一条分支&#xff0c;在这条分之上进行开发&#xff0c;开发完之后突然不想要这条分支了&#xff0c;想要撤销这些修改&#xff0c;前提是这些文件都还没有git add&#xff0c;方法&#xff1a; 切换到这条分支上 git checkout -- filepathname //放弃单个文件修改 gi…

Petya and Staircases CF212div.2B

http://codeforces.com/problemset/problem/362/B 水题。 1 #include<cstdio>2 #include<iostream>3 #include<cstring>4 #include<algorithm>5 using namespace std;6 int a[5000];7 int main()8 {9 int n,m,i; 10 scanf("%d%d",&a…