CN110110287A - 一种概率盒灵敏度分析方法 - Google Patents
一种概率盒灵敏度分析方法 Download PDFInfo
- Publication number
- CN110110287A CN110110287A CN201910306550.9A CN201910306550A CN110110287A CN 110110287 A CN110110287 A CN 110110287A CN 201910306550 A CN201910306550 A CN 201910306550A CN 110110287 A CN110110287 A CN 110110287A
- Authority
- CN
- China
- Prior art keywords
- probability box
- dss
- uncertain
- variables
- unc
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 17
- 238000010206 sensitivity analysis Methods 0.000 title claims abstract description 17
- 230000004927 fusion Effects 0.000 claims abstract description 21
- TVZRAEYQIKYCPH-UHFFFAOYSA-N 3-(trimethylsilyl)propane-1-sulfonic acid Chemical compound C[Si](C)(C)CCCS(O)(=O)=O TVZRAEYQIKYCPH-UHFFFAOYSA-N 0.000 claims abstract description 15
- 230000035945 sensitivity Effects 0.000 claims abstract description 11
- 239000008186 active pharmaceutical agent Substances 0.000 claims abstract description 7
- 239000000571 coke Substances 0.000 claims abstract description 6
- 238000004364 calculation method Methods 0.000 claims abstract description 5
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 239000006185 dispersion Substances 0.000 claims description 3
- 238000004458 analytical method Methods 0.000 abstract description 5
- 238000012360 testing method Methods 0.000 description 6
- 238000004422 calculation algorithm Methods 0.000 description 5
- 239000011159 matrix material Substances 0.000 description 4
- 238000009795 derivation Methods 0.000 description 3
- 239000012530 fluid Substances 0.000 description 2
- 230000005484 gravity Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 230000035939 shock Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Physics (AREA)
- Mathematical Analysis (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computational Mathematics (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Biology (AREA)
- Operations Research (AREA)
- Probability & Statistics with Applications (AREA)
- Complex Calculations (AREA)
Abstract
本发明涉及一种概率盒灵敏度分析方法，属于工业设备分析与应用技术领域。一种概率盒灵敏度分析方法，收集不确定变量的数据，建立K(K＝1，2，3…m)个不确定变量的概率盒模型；通过DS证据规则进行融合，然后对融合的焦元进行权重分配，最后进行焦元区间缩小，获得融合后新的DSS；将获得的新的DSS重新使用卷积计算获得融合后新的概率盒，定义融合后新的概率盒的累计不确定性统计量为T；定义去除变量A后得到的融合后概率盒累计不确定性统计量为R则变量A对融合后概率盒的灵敏度可以表达如下：100(1‑unc(R)/unc(T))％其中unc(R),unc(T)分别表示T和R的特征统计量。本发明使不确定变量的灵敏度分析结果更显著。
Description
技术领域
本发明涉及一种概率盒灵敏度分析方法，属于工业设备分析与应用技术领域。
背景技术
灵敏度分析有许多经典方法，如非参数法、方差分析法以及矩独立分析法。其中对可调参数不确定性描述和表达，尤其对随机不确定性表达方面，均值和方差较为常用，相当于用不确定可调参数数据的一个特征值描述不确定信息。在可调参数确定的情况下，直接求导法被用来计算灵敏度；在可调参数不确定且具有精确概率分布的情况下，常用方差分析法常被用来计算灵敏度。
在实际工况中，设备的可调参数通常是不确定的，用特征值描述不确定可调参数数据，会丢失丰富的统计信息。要获得不确定可调参数的精确概率分布，需要足量的试验样本。然而，由于工况条件和测试成本的限制，无法得到足量的试验样本。而利用有限的数据，容易获得其变化区间通常。但区间只能被定义不确定可调参数的范围，不能定义不确定可调参数的精确概率分布，这使得试验获取的统计信息不能被充分利用。概率盒被视为概率和区间的结合体，它实现了传统概率模型和区间的互补，能将诸多主观和客观的不确定因素被考虑进去，弥补了传统特征提取方法丢弃丰富概率统计信息的缺陷。然而，直接对概率盒用pinching算法进行灵敏度分析，存在指标不确定性区间的重叠交错现象。无法明显地区分各参数变量对目标变量的灵敏度。
概率盒理论的核心DS结构体(Dempster Shafer Structure，简称DSS)是由证据理论的作者提出来的，DSS和DS证据理论在主要观点上被保持了高度的一致性与互补性。为此，利用证据理论对随机集两两交叉组合，并进行卷积和计算，得到新的组合概率盒，用于计算单个可调参数概率盒对组合概率盒的灵敏度。因此，在概率盒的基础上，对不确定可调参数进行灵敏度分析具有一定的背景和意义。
发明内容
本发明针对现有技术存在的问题及不足，本发明提供一种概率盒灵敏度分析方法。本发明利用对不确定变量的概率盒建模来获取少量试验样本的概率分布，利用概率盒的外包络功能获取其变化区间。在不确定可调参数的概率盒模型基础上，利用证据理论对概率盒交叉组合，得到新的组合概率盒。然后通过pinching算法来分析不确定变量的灵敏度。本发明通过以下技术方案实现。
一种概率盒灵敏度分析方法，包括以下步骤：
步骤1、收集不确定变量的数据，建立K(K＝1，2，3…m)个不确定变量的概率盒模型，等离散率K＝m个概率盒模型变换转化后的DSS如下：
其中xi 、分别为DSS1中焦元的下上边界；y i、分别为DSS2中焦元的下上边界；z i、分别为DSSm中焦元的下上边界；mi为其对应的mass值，i为1，2，3…n；
步骤2、将DSS1中第i个焦元、DSS2中第j个焦元、...、DSSm中第k个焦元通过 DS证据规则进行融合，然后对融合的焦元进行权重分配，最后进行焦元区间缩小，获得融合后新的DSS；
步骤3、将步骤2获得的新的DSS重新使用卷积计算获得融合后新的概率盒，定义融合后新的概率盒的累计不确定性统计量为T；定义去除变量A后得到的融合后概率盒累计不确定性统计量为R则变量A对融合后概率盒的灵敏度可以表达如下：
100(1-unc(R)/unc(T))％
其中unc(R),unc(T)分别表示T和R的特征统计量。
上述概率盒灵敏度分析方法，以两个不确定变量的数据进行概率盒灵敏度分析方法举例，具体步骤包括：
步骤1、收集A和B两个不确定变量的数据，建立K(K＝1，2)个不确定变量的概率盒模型，等离散率K＝2个概率盒模型变换转化后的DSS如下：
其中xi 、分别为DSS1中焦元的下上边界；yi 、分别为DSS2中焦元的下上边界；mi为其对应的mass值，i为1，2，3…n；
步骤2、
步骤2.1、将DSSn1的第i个焦元与DSSn2的第j个焦元通过DS证据规则
进行融合，其中K为归一化因子，A、B分别对应DSS1和DSS2；m(A)为A焦元的mess 值，m(B)为B焦元的mess值，m(C)为A焦元与B焦元区间重叠处的mess值；对于任意的和m(C)≠0，其中
对融合的两焦元进行加权，定义DSS之间的加权矩阵表示为：
w＝[w1,w2...wn] (3)
其中：
w1+w2+…+wn＝1 (4)
wi为对应DSSi的加权权重，其中i为1，2，3…n。
步骤2.1、融合后得到新的焦元为对融合的两个概率盒进行权重分配，定义加权平均叠加因子矩阵为：
其中：
加权平均叠加因子矩阵作用示意如图1所示，从图1中可以看出相融合的两焦元矛盾越大，融合后得到的新的焦元越小。
步骤2.3、为了使数据叠加越大的焦元在融合时得到的焦元区间缩小，对焦元做如下改进：
z k＝z k·(1-c) (8)
将改进后所有的融合焦元组合获得融合后新的DSS；
步骤3、将步骤2.3得到融合后新的DSS重新使用卷积计算获得组合后新的概率盒。定义融合后概率盒的累计不确定性统计量为T；定义去除变量A(或B)后得到的融合后概率盒累计不确定性统计量为R；则变量A(或B)对融合后概率盒的灵敏度可以表达如下
100(1-unc(R)/unc(T))％ (10)
其中unc(R),unc(T)分别表示T和R的特征统计量。
本发明的有益效果是：本发明利用概率盒的外包络功能将诸多主观和客观的不确定因素都被概率盒建模考虑进去，弥补特征提取丢弃原始数据中的丰富概率统计信息的缺陷，让系统全面的不确定性能被真实反映。使不确定变量的灵敏度分析结果更显著。
附图说明
图1为本发明加权平均叠加因子矩阵作用示意图；
图2为本发明实施例1中超重力因子指标概率盒图；
图3为本发明实施例1中撞击角度指标概率盒图；
图4为本发明实施例1中喷嘴直径指标概率盒图；
图5为本发明实施例1中液体流量指标概率盒图；
图6为本发明实施例1中体积比指标概率盒图；
图7为本发明实施例1中撞击间距指标概率盒图；
图8为本发明实施例1中组合的概率盒图。
具体实施方式
下面结合附图和具体实施方式，对本发明作进一步说明。
实施例1
本实施例以超重力萃取器的参数调节数据为实验验证对象，分别获得了超重力因子、撞击角度、喷嘴直径、流量、体积比和撞击间距等六个可调参数的实测数据如表1 所示。
表1实验可调参数变量实测数据
上述概率盒灵敏度分析方法，包括以下步骤：
步骤1、收集六个可调参数不确定变量的数据，对每个可调参数进行等信度离散得到可调参数的区间变量，获得概率分布类型的不确定区间，建立6个不确定变量的概率盒模型，得到的超重力因子指标概率盒图如图2所示，撞击角度指标概率盒图如图3 所示，喷嘴直径指标概率盒图如图4所示，液体流量指标概率盒图如图5所示，体积比指标概率盒图如图6所示，撞击间距指标概率盒图如图7所示。
步骤2、通过DS证据规则将各指标的概率盒进行融合，得到融合后的概率盒如图 8所示，新的概率盒通过融合后新的DSS重新使用卷积计算获得，通过比较不确定参数常量化对融合概率盒带来的影响得到不确定参数的灵敏度结果。用求导法得到的灵敏度值和pinching算法得到的灵敏度值进行比较如表2所示。
表2灵敏度分析结果
从表2中可以看出采用现有技术中求导法得到的各指标之间灵敏度的差异不如本发明pinching算法得到的指标之间的灵敏度差异明显，可得本发明pinching算法得到的结果更为显著。
以上结合附图对本发明的具体实施方式作了详细说明，但是本发明并不限于上述实施方式，在本领域普通技术人员所具备的知识范围内，还可以在不脱离本发明宗旨的前提下作出各种变化。
Claims (1)
1.一种概率盒灵敏度分析方法，其特征在于：包括以下步骤：
步骤1、收集不确定变量的数据，建立K(K＝1，2，3…m)个不确定变量的概率盒模型，等离散率K＝m个概率盒模型变换转化后的DSS如下：
...
其中xi 、分别为DSS1中焦元的下上边界；yi 、分别为DSS2中焦元的下上边界；zi 、分别为DSSm中焦元的下上边界；mi为其对应的mass值，i为1，2，3…n；
步骤2、将DSS1中第i个焦元、DSS2中第j个焦元、...、DSSm中第k个焦元通过DS证据规则进行融合，然后对融合的焦元进行权重分配，最后进行焦元区间缩小，获得融合后新的DSS；
步骤3、将步骤2获得的新的DSS重新使用卷积计算获得融合后新的概率盒，定义融合后新的概率盒的累计不确定性统计量为T；定义去除变量A后得到的融合后概率盒累计不确定性统计量为R则变量A对融合后概率盒的灵敏度可以表达如下：
100(1-unc(R)/unc(T))％
其中unc(R),unc(T)分别表示T和R的特征统计量。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910306550.9A CN110110287A (zh) | 2019-04-17 | 2019-04-17 | 一种概率盒灵敏度分析方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910306550.9A CN110110287A (zh) | 2019-04-17 | 2019-04-17 | 一种概率盒灵敏度分析方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN110110287A true CN110110287A (zh) | 2019-08-09 |
Family
ID=67485613
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910306550.9A Pending CN110110287A (zh) | 2019-04-17 | 2019-04-17 | 一种概率盒灵敏度分析方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110110287A (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113051851A (zh) * | 2021-04-19 | 2021-06-29 | 电子科技大学 | 一种混合不确定性下的灵敏度分析方法 |
-
2019
- 2019-04-17 CN CN201910306550.9A patent/CN110110287A/zh active Pending
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113051851A (zh) * | 2021-04-19 | 2021-06-29 | 电子科技大学 | 一种混合不确定性下的灵敏度分析方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106549813A (zh) | 一种网络性能的评估方法及系统 | |
| WO2020248228A1 (zh) | 一种Hadoop平台计算节点负载预测方法 | |
| CN108446711A (zh) | 一种基于迁移学习的软件缺陷预测方法 | |
| CN106781486B (zh) | 基于浮动车数据的交通状态评价方法 | |
| CN107590565A (zh) | 一种构建建筑能耗预测模型的方法及装置 | |
| CN107038167A (zh) | 基于模型评估的大数据挖掘分析系统及其分析方法 | |
| CN105630743A (zh) | 一种光谱波数的选择方法 | |
| CN106407627A (zh) | 一种风速概率分布建模的方法和系统 | |
| CN109087277A (zh) | 一种基于图像特征和集成神经网络的空气细颗粒物pm2.5测量方法 | |
| CN108399255A (zh) | 一种分类数据挖掘模型的输入数据处理方法及装置 | |
| CN113435630B (zh) | 一种产流模式自适应的流域水文预报方法及系统 | |
| CN112926893A (zh) | 基于模糊综合评判和层次分析的水平井调剖效果评价方法 | |
| CN111832955B (zh) | 一种基于可靠度与多元统计的接触网状态评价方法 | |
| CN110110287A (zh) | 一种概率盒灵敏度分析方法 | |
| CN115526258A (zh) | 基于Spearman相关系数特征提取的电力系统暂稳评估方法 | |
| CN114580940A (zh) | 基于灰色关联度分析法的注浆效果模糊综合评价方法 | |
| CN107463532A (zh) | 一种电力统计数据的质量分析方法 | |
| CN108563201B (zh) | 一种dmaic驱动的小批量零件加工工艺过程质量改进方法 | |
| CN108711100A (zh) | 一种基于神经网络的p2p平台运营风险评估的系统 | |
| CN110210154B (zh) | 利用大坝测点数据判断测点表征大坝性态相似度的方法 | |
| CN105069214B (zh) | 一种基于非线性相关分析的工艺可靠性评估方法 | |
| CN108665090B (zh) | 基于主成分分析与Verhulst模型的城市电网饱和负荷预测方法 | |
| CN110909466A (zh) | 提高差孔隙分选储层平均毛管压力曲线计算精度的方法 | |
| CN110825583A (zh) | 一种针对云数据中心多指标融合的能效定性评估技术 | |
| CN108646688A (zh) | 一种基于回归学习的工艺参数优化分析方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190809 |
|
| RJ01 | Rejection of invention patent application after publication |