CN105718426B - 一种基于多元线性回归分析的疏浚产量数学模型建立方法 - Google Patents

一种基于多元线性回归分析的疏浚产量数学模型建立方法 Download PDF

Info

Publication number
CN105718426B
CN105718426B CN201610046597.2A CN201610046597A CN105718426B CN 105718426 B CN105718426 B CN 105718426B CN 201610046597 A CN201610046597 A CN 201610046597A CN 105718426 B CN105718426 B CN 105718426B
Authority
CN
China
Prior art keywords
decision
variable
dredging
yield
parameter
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.)
Expired - Fee Related
Application number
CN201610046597.2A
Other languages
English (en)
Other versions
CN105718426A (zh
Inventor
王祥冰
许焕敏
孟冬冬
吕品
李凯凯
孔德强
王亮
吴鹏程
吴渊凯
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Changzhou Campus of Hohai University
Original Assignee
Changzhou Campus of Hohai University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Changzhou Campus of Hohai University filed Critical Changzhou Campus of Hohai University
Priority to CN201610046597.2A priority Critical patent/CN105718426B/zh
Publication of CN105718426A publication Critical patent/CN105718426A/zh
Application granted granted Critical
Publication of CN105718426B publication Critical patent/CN105718426B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
    • 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

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)
  • Bioinformatics & Cheminformatics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Computing Systems (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Evolutionary Biology (AREA)
  • Operations Research (AREA)
  • Probability & Statistics with Applications (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Devices For Executing Special Programs (AREA)

Abstract

本发明公开了一种基于多元线性回归分析的疏浚产量数学模型建立方法，包括如下步骤：（1）收集影响疏浚作业施工工艺相关决策参量的数据资料，建立决策参量集；（2）构造决策参量矩阵，并对其进行规范化；（3）确立多元线性回归方程的一般形式；（4）建立数学模型系数的方程组，求解系数；（5）建立疏浚产量与参量的数学模型。本发明基于疏浚设备在施工作业中已获得的众多参量调控数据，建立疏浚产量与众多参量关系的数学模型，可为使疏浚设备达到高效作业而调节各个工艺参量提供最优决策参考，达到疏浚作业高产量的目的。

Description

一种基于多元线性回归分析的疏浚产量数学模型建立方法
技术领域
本发明涉及疏浚工程技术领域，特别是一种基于多元线性回归分析的疏浚产量数学模型建立方法。
背景技术
疏浚作为水下作业，工艺调控参量众多，其中众多参数变量之间常常存在一定的相关性。这势必增加了分析问题的复杂性。鉴于疏浚涉及多学科理论知识的综合性、疏浚作业过程的复杂性和实验条件的限制,疏浚工艺决策逻辑的基础理论研究远远落后于实际的需要；改变长期以来因为疏浚工艺知识与理论缺乏导致疏浚作业困难、效率不高的现状，指导疏浚作业参数调控，降低操作难度，为疏浚工艺决策智能化与自动化提供新的思路与理论依据。
回归分析的基本思想和方法以及“回归”名称的由来归功于英国统计学家F.高尔顿(F.Galton：1822—1911),该方法的根本特点是：回归分析是处理变量x与y之间的关系的一种统计方法和技术。
发明内容
本发明要解决的技术问题为：基于多元线性回归分析技术，研究一种疏浚产量数学模型建立方法，分析调控不同参变量之间的关系，以便对疏浚产量进行优化。
本发明采取的技术方案具体为：一种基于多元线性回归分析的疏浚产量数学模型建立方法，包括以下步骤：
步骤一，基于影响疏浚作业施工工艺相关决策变量的数据资料，建立决策变量集：
影响疏浚产量的决策参量为多个，所述决策变量集中的元素包括多个决策参量分别在不同决策方案中的取值，以及各决策方案对应的疏浚产量值；
步骤二，基于决策变量集构造决策变量矩阵，并对决策变量矩阵中的元素分别进行规范化；
步骤三，确立基于多元线性回归方程关于决策参量与疏浚产量之间关系的一般公式：
y＝β0+β1x1+β2x2+β3x3+...+βpxp (1)
式(1)中y为因变量，表示疏浚产量；x1，x2，...，xp为自变量，表示p个决策参量；β0，β1，…，βp是p+1个未知系数，β0为回归常数，β1，…，βp为回归系数；
步骤四，建立疏浚产量数学模型系数的方程组，将规范化后的决策变量矩阵中的元素代入方程组，求解式(1)中的未知系数；
步骤五，将求得的未知系数代入公式(1)，即得到疏浚产量的数学模型。
具体的，本发明步骤一中，定义决策变量集为U＝(u1,u2,…up)，其中，u1,u2,…up为决策变量，p为决策变量的个数；
设对p个决策变量进行n次观测，决策方案集为X＝(x1,x2,…xn)，其中xi＝|ai1ai2…aip|，表示第i次观测得到的决策方案，i＝1，2，……，n；aij表示第i次观测得到的第j个决策参量即自变量的原始数据，j＝1，2，……，p；
与决策方案集对应的，疏浚产量也即决策因变量Y＝(yi)，yi＝|bi|，i＝1，2，……，n；bi表示第i次观测得到的决策方案对应的疏浚产量值。
步骤二中，基于决策变量集建立的决策变量矩阵为：
对决策变量矩阵中的每个变量进行规范化处理的规范化公式为：
得到规范化后的决策变量矩阵：
可消除不同物理量纲对决策结果的影响，。
基于式(2)，根据微积分求极值原理，应满足以下方程组：
将规范化后的决策变量矩阵中的因变量和自变量对应的元素值对应代入方程组(3)，即可得到的取值，也即未知系数β0，β1，…，βp的取值。
有益效果
基于疏浚设备在施工作业中已获得的众多参量调控数据，建立疏浚产量与众多参量关系的数学模型，可为使疏浚设备达到高效作业而调节各个工艺参量提供最优决策参考，达到疏浚作业高产量的目的。为后来的各调控参量与产量的优化研究打下理论基础，对挖泥船的产量优化具有重要意义。
附图说明
图1所示为本发明方法流程示意图。
具体实施方式
以下结合附图和具体实施例进一步描述。
参考图1所示，本发明基于多元线性回归分析的疏浚产量数学模型建立方法，包括以下步骤：
步骤一，基于影响疏浚作业施工工艺相关决策变量的数据资料，建立决策变量集：
影响疏浚产量的决策参量为多个，所述决策变量集中的元素包括多个决策参量分别在不同决策方案中的取值，以及各决策方案对应的疏浚产量值；
步骤二，基于决策变量集构造决策变量矩阵，并对决策变量矩阵中的元素分别进行规范化；
步骤三，确立基于多元线性回归方程关于决策参量与疏浚产量之间关系的一般公式：
y＝β0+β1x1+β2x2+β3x3+...+βpxp (1)
式(1)中y为因变量，表示疏浚产量；x1，x2，...，xp为自变量，表示p个决策参量；β0，β1，…，βp是p+1个未知系数，β0为回归常数，β1，…，βp为回归系数；
步骤四，建立疏浚产量数学模型系数的方程组，将规范化后的决策变量矩阵中的元素代入方程组，求解式(1)中的未知系数；
步骤五，将求得的未知系数代入公式(1)，即得到疏浚产量的数学模型。
实施例
本实施例中，疏浚设备为绞吸式挖泥船，影响绞吸式挖泥船疏浚作业产量的主要工艺参量包括：绞刀转速，绞刀横移速度，绞刀切泥厚度，绞刀前进距离等，这几个工艺参量即为决策参量，p为8。
实施例中对上述8个决策参量进行了10次观测,得到决策变量集为如下表1：
表1
累计产量(y) 绞刀转速(x1) 台车行程(x2) 横移速度(x3) 管路平均浓度(x4) 管路流速(x5) 出口流速(x6) 绞刀切泥厚度(x7) 绞刀前进距离(x8)
682.53 27.89 2.76 13.46 57.45 4.7 12.02 17.86 2.74
698.53 27.81 2.76 13.28 57.8 4.7 12.02 17.71 2.74
709.97 27.87 2.76 12.77 55.86 4.65 11.9 17.73 2.74
720.82 27.86 2.75 0 52.56 4.58 11.74 17.81 0
732.86 27.25 3.17 1.64 49.03 4.55 11.64 17.89 0.42
745.49 27.47 4.15 10.75 45.93 4.5 11.53 17.69 1.41
754.27 27.47 4.16 11.78 42.5 4.43 11.33 17.79 1.41
761.45 27.66 4.14 12.01 38.59 4.38 11.21 17.86 1.39
771.92 27.34 4.14 12.9 34.9 4.35 11.13 17.73 1.39
785.72 26.37 4.13 12.06 31.84 4.35 11.14 17.93 1.39
步骤二中，对决策变量集进行规范化，具体如下：
根据决策方案集，为消除不同物理量纲对决策结果的影响，对每个决策变量aij，bi的原始数据进行如下规范化处理：
这样，就得到规范化后的决策变量集。如下表2：
表2
累计产量(y) 绞刀转速(x1) 台车行程(x2) 横移速度(x3) 管路平均浓度(x4) 管路流速(x5) 出口流速(x6) 绞刀切泥厚度(x7) 绞刀前进距离(x8)
0.02 0.75 0.00 0.94 0.04 0.13 0.13 0.40 0.00
0.03 0.57 0.00 0.97 0.08 0.28 0.27 0.76 0.00
0.04 0.38 0.00 0.96 0.12 0.45 0.46 0.53 0.00
0.05 0.61 0.00 0.98 0.14 0.59 0.60 0.29 0.00
0.05 0.74 0.15 0.12 0.16 0.75 0.74 0.04 0.20
0.06 0.19 0.28 0.96 0.18 0.83 0.83 0.64 0.40
0.06 0.64 0.28 0.91 0.20 0.89 0.88 0.49 0.40
0.07 0.86 0.28 0.95 0.24 0.95 0.95 0.00 0.40
步骤三中，确立基于多元线性回归方程的疏浚产量与决策参量关系模型的一般公式如下：
y＝β0+β1x1+β2x2+β3x3+...+βpxp (1)
其中，β0，β1，…，βp是p+1个未知参数，β0称为回归常数，β1，…，βp称为回归系数。y称为解释变量，即因变量，而x1，x2，...，xp是p个自变量。
步骤四中，建立数学模型系数的方程组，求解系数具体如下：
从(2)式，根据微积分求极值原理，应满足以下方
将标准化后多个决策方案中的自变量x和对应的因变量y值带入(3)式，解方程组可得多元线性回方程的系数。
解得系数如下：
β0＝0.03；β1＝-0.07；β2＝0.51；β3＝0.05；β4＝0.49；β5＝0.3；
β6＝-1.47；β7＝-0.02；β8＝-0.15
步骤五，建立关于疏浚产量与参量关系的数学模型：
将步骤四解得的系数取值带入式(1)即可得疏浚产量的数学模型为：
y＝0.03-0.07x1+0.51x2+0.05x3+0.49x4+1.3x5-1.47x6-0.02x7-0.15x8
其中自变量y对应产量，x1,x2,…xn对应各决策参量即调控参量。
以上所述仅是本发明的优选实施方式，应当指出：对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。

Claims (3)

1.一种基于多元线性回归分析的疏浚产量数学模型建立方法，其特征是，包括以下步骤：
步骤一，基于影响疏浚作业施工工艺相关决策变量的数据资料，建立决策变量集：
影响疏浚产量的决策参量为多个，所述决策变量集中的元素包括多个决策参量分别在不同决策方案中的取值，以及各决策方案对应的疏浚产量值；
步骤二，基于决策变量集构造决策变量矩阵，并对决策变量矩阵中的元素分别进行规范化；
步骤三，确立基于多元线性回归方程关于决策参量与疏浚产量之间关系的一般公式：
y＝β0+β1x1+β2x2+β3x3+...+βpxp (1)
式(1)中y为因变量，表示疏浚产量；x1，x2，...，xp为自变量，表示p个决策参量；β0，β1，…，βp是p+1个未知系数，β0为回归常数，β1，…，βp为回归系数；
步骤四，建立疏浚产量数学模型系数的方程组，将规范化后的决策变量矩阵中的元素代入方程组，求解式(1)中的未知系数；
步骤五，将求得的未知系数代入公式(1)，即得到疏浚产量的数学模型。
2.根据权利要求1所述的方法，其特征是，步骤一中，定义决策变量集为U＝(u1,u2,…up)，其中，u1,u2,…up为决策变量，p为决策变量的个数；
设对p个决策变量进行n次观测，决策方案集为X＝(x1,x2,…xn)，其中xi＝|ai1 ai2 …aip|，表示第i次观测得到的决策方案，i＝1，2，……，n；aij表示第i次观测得到的第j个决策参量即自变量的原始数据，j＝1，2，……，p；
与决策方案集对应的，疏浚产量也即决策因变量Y＝(yi)，yi＝|bi|，i＝1，2，……，n；bi表示第i次观测得到的决策方案对应的疏浚产量值。
3.根据权利要求1所述的方法，其特征是，步骤二中，基于决策变量集建立的决策变量矩阵为：
对决策变量矩阵中的每个变量进行规范化处理，得到规范化后的决策变量矩阵为：
矩阵R中，r10,r20,...,rn0为矩阵U'中决策因变量即疏浚产量值b1,b2,...,bn的规范化处理结果，rij(i＝1,2,...,n，j＝1,2,...,p)为矩阵U'中自变量即决策参量aij的规范化处理结果；
步骤四中，未知的回归常数β0和回归系数β1,β2,...,βp采用最小二乘法进行估计，即使得离差平方和达到最小，也即寻找未知数的取值，使其满足以下公式：
xi1,xi2,...,xip为疏浚产量即因变量yi＝|bi|时，对应的第i次观测得到的决策方案xi＝|ai1 ai2 … aip|中的决策参量即自变量ai1,ai2,...,aip；为离差平方和函数Q()中的自变量，满足公式(3)的自变量的取值即为β0,β1,...,βp的最小二乘法估计结果值；
基于式(2)，根据微积分求极值原理，应满足以下方程组：
将规范化后的决策变量矩阵中的因变量ri0和自变量ri1～rip对应的元素值对应代入方程组(3)，即可得到的取值，也即β0,β1,...,βp的取值。
CN201610046597.2A 2016-01-22 2016-01-22 一种基于多元线性回归分析的疏浚产量数学模型建立方法 Expired - Fee Related CN105718426B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610046597.2A CN105718426B (zh) 2016-01-22 2016-01-22 一种基于多元线性回归分析的疏浚产量数学模型建立方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610046597.2A CN105718426B (zh) 2016-01-22 2016-01-22 一种基于多元线性回归分析的疏浚产量数学模型建立方法

Publications (2)

Publication Number Publication Date
CN105718426A CN105718426A (zh) 2016-06-29
CN105718426B true CN105718426B (zh) 2018-07-27

Family

ID=56153899

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610046597.2A Expired - Fee Related CN105718426B (zh) 2016-01-22 2016-01-22 一种基于多元线性回归分析的疏浚产量数学模型建立方法

Country Status (1)

Country Link
CN (1) CN105718426B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108596489A (zh) * 2018-04-25 2018-09-28 河海大学常州校区 一种面向疏浚产量的工艺参量聚类分析方法
CN108573285A (zh) * 2018-04-25 2018-09-25 河海大学常州校区 一种面向疏浚能耗的工艺参量聚类分析方法
CN110309966B (zh) * 2019-06-27 2022-08-02 河海大学常州校区 基于偏最小二乘法的疏浚能耗预测方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7478019B2 (en) * 2005-01-26 2009-01-13 Kla-Tencor Corporation Multiple tool and structure analysis
CN102855221A (zh) * 2012-08-07 2013-01-02 中国科学院对地观测与数字地球科学中心 一种基于嵌套回归的有理函数模型参数自动优选及求解方法
CN104123451A (zh) * 2014-07-16 2014-10-29 河海大学常州校区 基于偏最小二乘回归的疏浚作业产量预测模型建立方法
CN104463359A (zh) * 2014-12-01 2015-03-25 河海大学常州校区 一种基于bp神经网络的疏浚作业产量预测模型分析方法
CN104573869A (zh) * 2015-01-16 2015-04-29 河海大学常州校区 基于bp神经网络与nsga-ⅱ实现疏浚作业的优化方法及系统
CN104794589A (zh) * 2015-05-11 2015-07-22 河海大学常州校区 一种基于层次分析的疏浚作业高效性分析方法
CN104850910A (zh) * 2015-05-25 2015-08-19 河海大学常州校区 一种基于owa的疏浚工艺参量多属性决策分析方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9355069B2 (en) * 2013-12-20 2016-05-31 Johnson Controls Technology Company Systems and methods for determining the uncertainty in parameters of an energy use model

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7478019B2 (en) * 2005-01-26 2009-01-13 Kla-Tencor Corporation Multiple tool and structure analysis
CN102855221A (zh) * 2012-08-07 2013-01-02 中国科学院对地观测与数字地球科学中心 一种基于嵌套回归的有理函数模型参数自动优选及求解方法
CN104123451A (zh) * 2014-07-16 2014-10-29 河海大学常州校区 基于偏最小二乘回归的疏浚作业产量预测模型建立方法
CN104463359A (zh) * 2014-12-01 2015-03-25 河海大学常州校区 一种基于bp神经网络的疏浚作业产量预测模型分析方法
CN104573869A (zh) * 2015-01-16 2015-04-29 河海大学常州校区 基于bp神经网络与nsga-ⅱ实现疏浚作业的优化方法及系统
CN104794589A (zh) * 2015-05-11 2015-07-22 河海大学常州校区 一种基于层次分析的疏浚作业高效性分析方法
CN104850910A (zh) * 2015-05-25 2015-08-19 河海大学常州校区 一种基于owa的疏浚工艺参量多属性决策分析方法

Also Published As

Publication number Publication date
CN105718426A (zh) 2016-06-29

Similar Documents

Publication Publication Date Title
CN105718426B (zh) 一种基于多元线性回归分析的疏浚产量数学模型建立方法
CN105447090B (zh) 一种自动化数据挖掘预处理方法
CN109262368A (zh) 一种刀具失效判定方法
CN103744293A (zh) 基于模糊神经网络的废水处理监控方法和系统
Zhang et al. Delay-dependent robust stability of neutral systems with mixed delays and nonlinear perturbations
CN108399298B (zh) 一种管材切割的套料算法
CN106292292A (zh) 基于案例推理的铁矿浮选药剂用量优化设定方法及系统
CN105574587A (zh) 一种塑料注塑过程的在线工况过程监控方法
CN109115992A (zh) 一种农田土壤健康评价系统
CN104050547A (zh) 一种油田开发规划方案非线性优选决策方法
CN104573869A (zh) 基于bp神经网络与nsga-ⅱ实现疏浚作业的优化方法及系统
CN107403196A (zh) 基于谱聚类分析的即时学习建模方法
CN102521690A (zh) 一种企业产品订单的优化方法
Hayashi et al. Generation of fuzzy decision trees by fuzzy ID3 with adjusting mechanism of and/or operators
CN105787113A (zh) 一种基于plm数据库面向dpipp工艺信息的挖掘算法
CN115879047A (zh) 一种基于梯度提升决策的煤矿井下钻进过程工况识别方法
CN104850910A (zh) 一种基于owa的疏浚工艺参量多属性决策分析方法
CN113023873B (zh) 用于污水处理的智能管理系统及方法
Baxter et al. Model-based advanced process control of coagulation
CN104881530A (zh) 一种基于优化工艺参数的滚齿干切加工方法
CN110134008A (zh) 一种基于经济性和统计性的mewma控制图的设计方法
CN103617726A (zh) 粗粒度交通状态转换为细粒度速度的方法
Simsek et al. ESTIMATION OF NUTRIENT CONCENTRATIONS IN RUNOFF FROM BEEF CATTLE FEEDLOT USING ADAPTIVE NEURO-FUZZY INFERENCE SYSTEMS.
CN112445256A (zh) 一种基于传感识别技术的植物种植管理系统
CN104573356A (zh) 一种基于稀疏T-S模糊的铁水Si含量建模方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20180727

Termination date: 20210122