CN105718426A - 一种基于多元线性回归分析的疏浚产量数学模型建立方法 - Google Patents
一种基于多元线性回归分析的疏浚产量数学模型建立方法 Download PDFInfo
- Publication number
- CN105718426A CN105718426A CN201610046597.2A CN201610046597A CN105718426A CN 105718426 A CN105718426 A CN 105718426A CN 201610046597 A CN201610046597 A CN 201610046597A CN 105718426 A CN105718426 A CN 105718426A
- Authority
- CN
- China
- Prior art keywords
- beta
- decision
- dredging
- variable
- yield
- 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.)
- Granted
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)
- Devices For Executing Special Programs (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (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次观测得到的决策方案对应的疏浚产量值。
步骤二中,基于决策变量集建立的决策变量矩阵为:
对决策变量矩阵中的每个变量进行规范化处理的规范化公式为:
得到规范化后的决策变量矩阵:
可消除不同物理量纲对决策结果的影响,。
步骤四中,未知系数β0,β1,…,βp采用最小二乘法进行估计,即使离差平方和 达到最小,也即寻找的取值,使其满足以下公式:
基于式(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个自变量。
步骤四中,建立数学模型系数的方程组,求解系数具体如下:
多元线性回归方程未知参数β0,β1,…,β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 (4)
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=|ai1ai2…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所述的方法,其特征是,步骤二中,基于决策变量集建立的决策变量矩阵为:
对决策变量矩阵中的每个变量进行规范化处理的规范化公式为:
得到规范化后的决策变量矩阵:
4.根据权利要求1所述的方法,其特征是,步骤四中,未知系数β0,β1,…,βp采用最小二乘法进行估计,即使离差平方和 达到最小,也即寻找的取值,使其满足以下公式:
基于式(2),根据微积分求极值原理,应满足以下方程组:
将规范化后的决策变量矩阵中的因变量和自变量对应的元素值对应代入方程组(3),即可得到的取值,也即未知系数β0,β1,…,βp的取值。
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 true CN105718426A (zh) | 2016-06-29 |
CN105718426B 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) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108573285A (zh) * | 2018-04-25 | 2018-09-25 | 河海大学常州校区 | 一种面向疏浚能耗的工艺参量聚类分析方法 |
CN108596489A (zh) * | 2018-04-25 | 2018-09-28 | 河海大学常州校区 | 一种面向疏浚产量的工艺参量聚类分析方法 |
CN110309966A (zh) * | 2019-06-27 | 2019-10-08 | 河海大学常州校区 | 基于偏最小二乘法的疏浚能耗预测方法 |
Citations (8)
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-ⅱ实现疏浚作业的优化方法及系统 |
US20150178633A1 (en) * | 2013-12-20 | 2015-06-25 | Johnson Controls Technology Company | Systems and methods for determining the uncertainty in parameters of an energy use model |
CN104794589A (zh) * | 2015-05-11 | 2015-07-22 | 河海大学常州校区 | 一种基于层次分析的疏浚作业高效性分析方法 |
CN104850910A (zh) * | 2015-05-25 | 2015-08-19 | 河海大学常州校区 | 一种基于owa的疏浚工艺参量多属性决策分析方法 |
-
2016
- 2016-01-22 CN CN201610046597.2A patent/CN105718426B/zh not_active Expired - Fee Related
Patent Citations (8)
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 | 中国科学院对地观测与数字地球科学中心 | 一种基于嵌套回归的有理函数模型参数自动优选及求解方法 |
US20150178633A1 (en) * | 2013-12-20 | 2015-06-25 | Johnson Controls Technology Company | Systems and methods for determining the uncertainty in parameters of an energy use model |
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的疏浚工艺参量多属性决策分析方法 |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108573285A (zh) * | 2018-04-25 | 2018-09-25 | 河海大学常州校区 | 一种面向疏浚能耗的工艺参量聚类分析方法 |
CN108596489A (zh) * | 2018-04-25 | 2018-09-28 | 河海大学常州校区 | 一种面向疏浚产量的工艺参量聚类分析方法 |
CN110309966A (zh) * | 2019-06-27 | 2019-10-08 | 河海大学常州校区 | 基于偏最小二乘法的疏浚能耗预测方法 |
CN110309966B (zh) * | 2019-06-27 | 2022-08-02 | 河海大学常州校区 | 基于偏最小二乘法的疏浚能耗预测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105718426B (zh) | 2018-07-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105718426A (zh) | 一种基于多元线性回归分析的疏浚产量数学模型建立方法 | |
Mat Deris et al. | Hybrid GR-SVM for prediction of surface roughness in abrasive water jet machining | |
Afshar et al. | Application of cellular automata to sewer network optimization problems | |
CN108388689B (zh) | 一种基于量纲分析的盾构刀盘扭矩计算方法 | |
CN106292292A (zh) | 基于案例推理的铁矿浮选药剂用量优化设定方法及系统 | |
CN104050547A (zh) | 一种油田开发规划方案非线性优选决策方法 | |
CN105956235A (zh) | 基于svr-pso的超声加工特种刀具优化设计方法 | |
CN113130014B (zh) | 一种基于多分支神经网络的稀土萃取模拟方法及系统 | |
CN106934074A (zh) | 全局最优涡扇发动机进气道减噪设计方法 | |
CN104077432A (zh) | 一种基于多维关联函数的调整工艺选择分析方法 | |
Charankumar et al. | An empirical study on modified second order response surface and Taguchi designs for optimizing the process parameters | |
CN105389442A (zh) | 一种耦合遗传算法、神经网络与数值模拟的反向设计方法 | |
CN105787113A (zh) | 一种基于plm数据库面向dpipp工艺信息的挖掘算法 | |
Xiao et al. | Ant colony system algorithm for the optimization of beer fermentation control | |
CN105512383A (zh) | 基于bp神经网络的疏浚工艺调控参量筛选方法 | |
CN104850910A (zh) | 一种基于owa的疏浚工艺参量多属性决策分析方法 | |
CN109989747B (zh) | 一种油田注采井组开发非均匀程度定量表征方法 | |
CN105894204A (zh) | 一种基于prolog的疏浚决策分析方法 | |
Ghasemi et al. | Classification of seismic vulnerability based on machine learning techniques for RC frames | |
CN104794589A (zh) | 一种基于层次分析的疏浚作业高效性分析方法 | |
Kariwala et al. | Efficient branch and bound methods for pairing selection | |
Dashkovskiy et al. | Application of the LISS Lyapunov-Krasovskii small-gain theorem to autonomously controlled production networks with time-delays | |
Ahlawat et al. | Fuzzy expert system for greenhouse rose yield prediction | |
CN104050370B (zh) | 基于二分法进行动力定位能力分析的方法 | |
Fan et al. | Typical machining process mining of servo valve parts based on self-adaptive affinity propagation clustering |
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 |
Granted publication date: 20180727 Termination date: 20210122 |
|
CF01 | Termination of patent right due to non-payment of annual fee |