CN107368928A - 一种古建筑沉降的组合预测方法及系统 - Google Patents
一种古建筑沉降的组合预测方法及系统 Download PDFInfo
- Publication number
- CN107368928A CN107368928A CN201710654232.2A CN201710654232A CN107368928A CN 107368928 A CN107368928 A CN 107368928A CN 201710654232 A CN201710654232 A CN 201710654232A CN 107368928 A CN107368928 A CN 107368928A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- sequence
- msup
- mtd
- 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
- 238000004062 sedimentation Methods 0.000 title claims abstract description 79
- 238000013277 forecasting method Methods 0.000 title claims abstract description 10
- 238000000034 method Methods 0.000 claims abstract description 30
- 230000001537 neural effect Effects 0.000 claims abstract description 10
- 230000001186 cumulative effect Effects 0.000 claims description 21
- 238000010606 normalization Methods 0.000 claims description 18
- 238000001914 filtration Methods 0.000 claims description 17
- 230000004044 response Effects 0.000 claims description 11
- 239000011888 foil Substances 0.000 claims description 10
- 238000004422 calculation algorithm Methods 0.000 claims description 7
- 238000011161 development Methods 0.000 claims description 7
- 230000001373 regressive effect Effects 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 3
- 230000008901 benefit Effects 0.000 abstract description 3
- 230000007935 neutral effect Effects 0.000 abstract description 2
- 238000005457 optimization Methods 0.000 abstract description 2
- 230000006870 function Effects 0.000 description 9
- 238000013528 artificial neural network Methods 0.000 description 7
- 238000004590 computer program Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 7
- 238000012545 processing Methods 0.000 description 6
- 239000003016 pheromone Substances 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 238000012986 modification Methods 0.000 description 4
- 230000004048 modification Effects 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 238000006073 displacement reaction Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 230000035772 mutation Effects 0.000 description 2
- 239000002689 soil Substances 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 238000003860 storage Methods 0.000 description 2
- 241000257303 Hymenoptera Species 0.000 description 1
- LTXREWYXXSTFRX-QGZVFWFLSA-N Linagliptin Chemical compound N=1C=2N(C)C(=O)N(CC=3N=C4C=CC=CC4=C(C)N=3)C(=O)C=2N(CC#CC)C=1N1CCC[C@@H](N)C1 LTXREWYXXSTFRX-QGZVFWFLSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000006378 damage Effects 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 239000003673 groundwater Substances 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000003062 neural network model Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 238000005096 rolling process Methods 0.000 description 1
- 238000012163 sequencing technique Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000001960 triggered effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/004—Artificial life, i.e. computing arrangements simulating life
- G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/045—Combinations of networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Tourism & Hospitality (AREA)
- Health & Medical Sciences (AREA)
- Economics (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- Life Sciences & Earth Sciences (AREA)
- Evolutionary Computation (AREA)
- Software Systems (AREA)
- Marketing (AREA)
- Mathematical Physics (AREA)
- Development Economics (AREA)
- Artificial Intelligence (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Computational Linguistics (AREA)
- Data Mining & Analysis (AREA)
- General Business, Economics & Management (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- Educational Administration (AREA)
- Primary Health Care (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- External Artificial Organs (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供了一种古建筑沉降的组合预测方法及系统,涉及古建筑沉降预测技术领域。利用Verhulst灰色模型通过采用等维递补滚动方法进行多次预测,可以利用较少的样本进行趋势预测,对沉降数据的趋势进行总体把握。然后,然后将初步预测结果输入RBF神经网络,充分发挥RBF局部寻优的优势。该模型既可以有效避免BP等神经网络的全局优化过拟合现象,提高预测的精度的同时,又可以规避Verhuslt灰色模型多步预测的低精度问题,真正实现古建筑沉降数据的高精度多步预测。
Description
技术领域
本发明涉及古建筑沉降预测技术领域,特别涉及一种古建筑沉降的组合预测方法及系统。
背景技术
当今,文物保护模式逐渐从“抢救性”保护过渡到“预防性”保护,利用物联网技术实现文物多角度全方位监测,为文物病害发展趋势预测提供了丰富的数据支持,为文物保护和修复方案的制定提供了重要依据,避免因决策不当造成文化遗产价值遗失。文化遗产包含古墓葬、古建筑等不可移动文物。古建筑沉降是由于建筑物荷载通过基础底面传递给地基,使天然土层原有的应力状态发生变化,当地基荷载不均匀时,就会造成不均匀沉降,从而给古建筑完整性和安全性带来隐患。如果沉降变形估计不足,轻则出现裂缝影响古建筑观瞻性和稳定性,重则会引发结构坍塌和严重人员伤亡事故。因此,提高古建筑沉降的预测精度至关重要。
造成古建筑沉降的因素包括地下水开采、人为破坏、自身结构等多种随机性、不确定性因素,因此,沉降数据具有非线性、复杂性及多模态性等特点,致使沉降预测存在诸多困难。目前计算与预测沉降形变的理论方法可以分为两类:一类是模型法,即基于土体本构关系的数值方法,从地面沉降的整个发展过程来分析考虑,建立沉降量与时间的关系模型;另一类是基于实测数据的智能分析方法。
模型法主要采用渗透固结理论、静态预测法和灰色系统法等。其中,渗透固结理论由于参数获得、理论假设等方面的问题,其预测误差往往较大;静态预测法是建立建筑物沉降与时间的拟合关系,常用的拟合函数有双曲线、指数函数等,由于拟合函数事先给定,使得该方法适用范围受到较大限制,且难以反映全过程的沉降与时间的关系;灰色预测模型适用于对非线性、不确定系统的时间序列数据进行预测,尤其对“贫信息”的小样本量作总体趋势预测。张庆伟、刘毅、王亚光等采用灰色模型对地面沉降数据进行了预测,从预测结果来看,灰色模型初期预测较好,但后期收敛较慢,因沉降的多变性和复杂性以及试验参数获取的困难性,不能达到预测精度的要求。尤其是当系统中出现了突变、故障等情况时,由于数据的突变,破坏了数据本身的平稳性,造成预测误差大幅上升。
在智能分析方法中,很多学者利用BP人工神经网络、小波分析等方法对沉降预测进行了研究。BP神经网络模型具有自学习、容错性强、计算简单、并行处理速度快等优点,在理论上可以任意逼近任何非线性映射。Yong、Carles等采用人工神经网络(ArtificialNeural Network,ANN)、高斯过程(Gaussian processes,GP)法进行沉降预测。OCAK等通过BP神经网络对隧道围岩变形了预测。由于BP神经网络是一种前馈神经网络,而沉降会有一定的滞后性,与前几年的沉降量存在一定关联性,使得BP神经网络在中短期预测方面存在很大局限性。郭健运用小波分析对沉降监测数据去噪处理后,构建了W-RBF预测模型,验证了运用RBF神经网络进行沉降预测的可行性。从现有模型对路基沉降的预测情况看,多数模型只适合中短期预测,而中长期的预测能力(预测长度和预测精度)较弱,同时在很多情况下缺少对模型可靠性的验证。
发明内容
本发明实施例提供了一种古建筑沉降的组合预测方法及系统,用以解决现有技术中存在的问题。
一种古建筑沉降的组合预测方法,所述方法包括:
对应变片采集的原始沉降序列进行五点三次平滑滤波,将所述原始沉降序列中的噪声滤除,得到相应的初始沉降序列;
对所述初始沉降序列采用Verhulst灰色模型进行预测,得到相应的初步预测沉降序列;
对所述初步预测沉降序列进行归一化处理,得到所述初步预测沉降序列归一化后的初步预测沉降序列;
采用基于蚁群聚类算法的RBF神经网络对所述归一化后的初步预测沉降序列进行预测,获得沉降预测序列;
将所述沉降预测序列中的元素进行反归一化操作,获得反归一化后的沉降预测序列。
优选地,步骤对应变片采集的原始沉降序列进行五点三次平滑滤波,将所述原始沉降序列中的噪声滤除,得到相应的初始沉降序列包括:
对所述应变片采集的k组原始沉降序列
Xi={xi1,xi2,...,xip,...,xim},i=1,2,...,k;p=1,2,...,m进行按照以下公式进行五点三次平滑滤波,得到所述初始沉降序列
Xi′={xi1′,xi2′,...,xip′,...,xim′},i=1,2,...,k;p=1,2,...,m;
其中m为每个所述原始沉降序列中数据的数量。
优选地,步骤对所述初始沉降序列采用Verhulst灰色模型进行预测,得到相应的初步预测沉降序列包括:
对于所述初始沉降序列X′i进行累加运算,获得累加序列
Xi (1)=(xi1 (1),xi2 (1),...,xip (1),...,xim (1)),i=1,2,...,k;p=1,2,...,m,其中
计算所述累加序列Xi (1)的紧邻均值生成数,获得紧邻均值生成序列
Zi (1)=(zi2 (1),zi3 (1),...,zip (1),...,zim (1)),i=1,2,...,k;p=1,2,...,m,其中
zip (1)=0.5xip (1)+0.5xi(p-1) (1);
建立所述累加序列Xi (1)和紧邻均值生成序列Zi (1)之间的灰色Verhulst微分方程:
xip (1)+azip (1)=b(zip (1))2,i=1,2,...,k;p=1,2,...,m
式中,a为发展系数,b为灰色系数,求解所述微分方程得到所述微分方程的时间响应序列:
其中,所述发展系数a和灰色系数b通过最小二乘法确定,即
式中,为灰色参数的估计向量,B为自变量矩阵,Y为估计预测值向量,其表达式分别为:
将所述时间响应序列进行累减运算,获得所述初始沉降序列Xi′的初步预测沉降序列其中
优选地,步骤对所述初步预测沉降序列进行归一化处理,得到所述初步预测沉降序列归一化后的初步预测沉降序列包括:
设为中的最大元素,为中的最小元素,则:
其中,为所述归一化后的初步预测沉降序列中的元素;
步骤将所述沉降预测序列中的元素进行反归一化操作,获得反归一化后的沉降预测序列包括:
设为所述沉降预测序列中的最大元素,为中的最小元素,则:
其中,为所述反归一化后的沉降预测序列中的元素,r为聚类半径。
优选地,在计算所述初步预测沉降序列时,每预测一个初步预测沉降值后,在所述初始沉降序列X′i中添加一个新的初始沉降值,并将所述初始沉降序列X′i中的首个元素去除,保持所述初始沉降序列X′i的序列长度不变,再次进行预测。
一种古建筑沉降的组合预测系统,所述系统包括:
预处理模块,用于对应变片采集的k组原始沉降序列进行五点三次平滑滤波,将所述原始沉降序列中的噪声滤除,得到相应的初始沉降序列;
灰色模型初步预测模块,用于对所述初始沉降序列采用Verhulst灰色模型进行预测,得到相应的初步预测沉降序列;
归一化模块,用于对所述初步预测沉降序列进行归一化处理,得到所述初步预测沉降序列归一化后的初步预测沉降序列;
沉降预测模块,用于采用基于蚁群聚类算法的RBF神经网络对所述归一化后的初步预测沉降序列进行预测,获得沉降预测序列;
反归一化模块,用于将所述沉降预测序列中的元素进行反归一化操作,获得反归一化后的沉降预测序列。
优选地,所述灰色模型初步预测模块包括:
累加子模块,用于对于所述初始沉降序列进行累加运算,获得累加序列;
紧邻均值计算子模块,计算所述累加序列的紧邻均值生成数,获得紧邻均值生成序列;
灰色微分方程建立子模块,用于建立所述累加序列和紧邻均值生成序列之间的灰色Verhulst微分方程,求解所述微分方程得到所述微分方程的时间响应序列;
序列估计子模块,用于对所述时间响应序列进行累减运算,获得所述初始沉降序列的初步预测沉降序列。
本发明的有益效果为:利用Verhulst灰色模型通过采用等维递补滚动方法进行多次预测,可以利用较少的样本进行趋势预测,对沉降数据的趋势进行总体把握。然后,然后将初步预测结果输入RBF神经网络,充分发挥RBF局部寻优的优势。该模型既可以有效避免BP等神经网络的全局优化过拟合现象,提高预测的精度的同时,又可以规避Verhuslt灰色模型多步预测的低精度问题,真正实现古建筑沉降数据的高精度多步预测。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的一种古建筑沉降的组合预测方法的步骤流程图;
图2为本发明实施例提供的一种古建筑沉降的组合预测系统的功能模块图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
在详细说明本发明的技术方案前,首先对一些定义进行解释:
定义1:设序列X=(x(1),x(2),...,x(k),x(k+1),...,x(n)),其中x(k)和x(k+1)为X的一对紧邻值。
定义2:X(0)为原始序列,对X(0)中的每个元素进行运算,若则称为X(0)的一次累加生成算子,记为1-AGO;
若x(0)(k)=x(0)(k)-x(0)(k-1),k=1,2,...,n,则称为X(0)的一次累减生成算子,记为1-IAGO。
定义3:设序列X=(x(1),x(2),...,x(n)),x*(k)=0.5x(k)+0.5x(k-1),则称x*(k)为X的紧邻均值生成数,由紧邻均值生成数构成的序列称为紧邻均值生成序列。
参照图1,本发明实施例提供了一种古建筑位移的组合预测方法,该方法包括:
步骤100,对应变片采集的k组原始沉降序列
Xi={xi1,xi2,...,xip,...,xim},i=1,2,...,k;p=1,2,...,m进行五点三次平滑滤波,将所述原始沉降序列中的噪声滤除,得到相应的初始沉降序列
X′i={x′i1,x′i2,...,x′ip,...,x′im},i=1,2,...,k;p=1,2,...,m。
具体地,Xi′中的各个元素通过以下公式计算得到:
步骤120,对所述初始沉降序列X′i采用Verhulst灰色模型进行预测,得到相应的初步预测沉降序列在计算所述初步预测沉降序列时,采用等维递补滚动预测方法,即每预测一个初步预测沉降值后,在所述初始沉降序列X′i中添加一个新的初始沉降值,并将所述初始沉降序列X′i中的首个元素去除,保持所述初始沉降序列X′i的序列长度不变,再次进行预测。
具体地,步骤120包括:
子步骤121,对于所述初始沉降序列X′i按照定义2进行累加运算,获得累加序列Xi (1)=(xi1 (1),xi2 (1),...,xip (1),...,xim (1)),i=1,2,...,k;p=1,2,...,m;
子步骤122,对于所述累加序列Xi (1)按照定义3运算,获得紧邻均值生成序列Zi (1)=(zi2 (1),zi3 (1),...,zip (1),...,zim (1)),i=1,2,...,k;p=1,2,...,m;
子步骤123,建立所述累加序列Xi (1)和紧邻均值生成序列Zi (1)之间的灰色Verhulst微分方程:
xip (1)+azip (1)=b(zip (1))2,i=1,2,...,k;p=1,2,...,m
式中,a为发展系数,b为灰色系数,求解所述微分方程得到其时间响应序列:
其中,所述发展系数a和灰色系数b通过最小二乘法确定,即
式中,为灰色参数的估计向量,B为自变量矩阵,Y为估计预测值向量,其表达式分别为:
子步骤124,将所述时间响应序列按照定义2进行累减运算,获得所述初始沉降序列X′i的初步预测沉降序列
步骤140,对所述初步预测沉降序列进行归一化处理,得到所述初步预测沉降序列归一化后的初步预测沉降序列
具体地,设为中的最大元素,为中的最小元素,则:
其中,为中的元素,通过上式即可将中的元素归一化到0.1~0.9之间。
步骤160,采用基于蚁群聚类算法的RBF神经网络对所述归一化后的初步预测沉降序列进行预测,获得沉降预测序列
具体地,步骤130包括:
子步骤161,将所述归一化后的初步预测沉降序列作为样本点集合,通过聚类算法将中的样本点聚为n个聚类,n个聚类的集合为C={cq|q=1,2,...,n},其中且任意两个聚类cq的交集为空。
子步骤162,初始化信息启发因子α、期望启发因子β、信息素挥发系数ρ和信息素强度Q、聚类半径r,以及蚂蚁个数M,最大进化代数N,令聚类中心cq的初始值为聚类中心cq对应的径向基函数的样本中心点初始值。
子步骤163,计算所述样本点集中每个样本点在t时刻到每个聚类cq的欧氏距离giq(t)以及启发函数
子步骤164,计算t时刻样本点到聚类中心cq的路径上的信息素:
子步骤165,将每只蚂蚁爬过的所述样本点集中的样本点构成的聚类结果构成一个解,即所述径向基函数的解。其中,蚂蚁随机选择一个所述样本点集中的样本点该样本点被聚类到聚类中心cq的概率为:
其中,
子步骤166,采用轮盘赌选择法确定所属的聚类中心。
子步骤167,更新聚类中心cq,通过以下公式确定新的聚类中心:
其中,为更新后的聚类中心,J为原聚类中心cq中的样本点数量。
子步骤168,当所有蚂蚁完成一次周游后,各路径上的信息素更新为τiq(t)=ρτiq(t)+Δτiq(t),Δτiq(t)为信息素的增加量,其表达式为:
其中lmb为最优聚类结果中各模式样本到其聚类中心的距离之和,表示
子步骤169,计算n个聚类中每个聚类的样本点到各自更新后的聚类中心的欧式距离之和即若F达到预设值,则聚类成功,否则重新聚类。
子步骤170,设RBF网络的径向基向量为H(x)=[h1,h2,...,hm]T,选取高斯函数作为径向基函数,并将子步骤169中计算得到的聚类中心作为径向基函数的中心,则径向基函数表达式为σq为扩展常数,cq为RBF网络的第q个结点的中心矢量。
子步骤171,根据cq以及以下公式计算所述RBF神经网络的输出:
式中,为沉降预测序列中的元素,ωiq为径向基函数的连接权值。根据梯度下降算法,其连接权值ωiq、径向基函数的中心点以及扩展常数σq(t)的优化表达式为:
式中,γ为学习率,取值范围为0<γ<2,yj为预测目标值,η2和η3表示梯度下降算法中的学习率,表示误差函数。
步骤180,将所述沉降预测序列中的元素按照以下公式进行反归一化操作,获得反归一化后的沉降预测序列
其中,为中的元素,为中的最大元素,为中的最小元素。
基于同一发明构思,本发明实施例还提供了一种古建筑位移的组合预测系统,如图2所示,由于该系统解决技术问题的原理和一种古建筑位移的组合预测方法相似,因此该系统的实施可参照方法的实施,重复之处不再赘述。
预处理模块200,用于对应变片采集的k组原始沉降序列进行五点三次平滑滤波,将所述原始沉降序列中的噪声滤除,得到相应的初始沉降序列。
灰色模型初步预测模块210,用于对所述初始沉降序列采用Verhulst灰色模型进行预测,得到相应的初步预测沉降序列。
归一化模块220,用于对所述初步预测沉降序列进行归一化处理,得到所述初步预测沉降序列归一化后的初步预测沉降序列。
沉降预测模块230,用于采用基于蚁群聚类算法的RBF神经网络对所述归一化后的初步预测沉降序列进行预测,获得沉降预测序列。
反归一化模块240,用于将所述沉降预测序列中的元素进行反归一化操作,获得反归一化后的沉降预测序列。
其中所述灰色模型初步预测模块210包括以下子模块:
累加子模块211,用于对于所述初始沉降序列进行累加运算,获得累加序列。
紧邻均值计算子模块212,计算所述累加序列的紧邻均值生成数,获得紧邻均值生成序列。
灰色微分方程建立子模块213,用于建立所述累加序列和紧邻均值生成序列之间的灰色Verhulst微分方程,求解所述微分方程得到所述微分方程的时间响应序列。
序列估计子模块214,用于对所述时间响应序列进行累减运算,获得所述初始沉降序列的初步预测沉降序列。
本领域内的技术人员应明白,本发明的实施例可提供为方法、系统、或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
尽管已描述了本发明的优选实施例,但本领域内的技术人员一旦得知了基本创造性概念,则可对这些实施例作出另外的变更和修改。所以,所附权利要求意欲解释为包括优选实施例以及落入本发明范围的所有变更和修改。
显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包含这些改动和变型在内。
Claims (7)
1.一种古建筑沉降的组合预测方法,其特征在于,所述方法包括:
对应变片采集的原始沉降序列进行五点三次平滑滤波,将所述原始沉降序列中的噪声滤除,得到相应的初始沉降序列;
对所述初始沉降序列采用Verhulst灰色模型进行预测,得到相应的初步预测沉降序列;
对所述初步预测沉降序列进行归一化处理,得到所述初步预测沉降序列归一化后的初步预测沉降序列;
采用基于蚁群聚类算法的RBF神经网络对所述归一化后的初步预测沉降序列进行预测,获得沉降预测序列;
将所述沉降预测序列中的元素进行反归一化操作,获得反归一化后的沉降预测序列。
2.如权利要求1所述的方法,其特征在于,步骤对应变片采集的原始沉降序列进行五点三次平滑滤波,将所述原始沉降序列中的噪声滤除,得到相应的初始沉降序列包括:
对所述应变片采集的k组原始沉降序列Xi={xi1,xi2,...,xip,...,xim},i=1,2,...,k;p=1,2,...,m进行按照以下公式进行五点三次平滑滤波,得到所述初始沉降序列Xi′={xi1′,xi2′,...,xip′,...,xim′},i=1,2,...,k;p=1,2,...,m;
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msup>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>&prime;</mo>
</msup>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>70</mn>
</mfrac>
<mo>&lsqb;</mo>
<mn>69</mn>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>+</mo>
<mn>4</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mn>4</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mn>6</mn>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mn>3</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mn>5</mn>
</mrow>
</msub>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>&prime;</mo>
</msup>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>35</mn>
</mfrac>
<mo>&lsqb;</mo>
<mn>2</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mn>5</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mn>27</mn>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>+</mo>
<mn>12</mn>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mn>3</mn>
</mrow>
</msub>
<mo>-</mo>
<mn>8</mn>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mn>4</mn>
</mrow>
</msub>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mi>p</mi>
</mrow>
</msub>
<mo>&prime;</mo>
</msup>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>35</mn>
</mfrac>
<mo>&lsqb;</mo>
<mo>-</mo>
<mn>3</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mrow>
<mo>(</mo>
<mi>p</mi>
<mo>-</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mrow>
<mo>(</mo>
<mi>p</mi>
<mo>+</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mn>12</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mrow>
<mo>(</mo>
<mi>p</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mrow>
<mo>(</mo>
<mi>p</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mn>17</mn>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mi>p</mi>
</mrow>
</msub>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>&prime;</mo>
</msup>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>35</mn>
</mfrac>
<mo>&lsqb;</mo>
<mn>2</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mrow>
<mo>(</mo>
<mi>p</mi>
<mo>-</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mn>12</mn>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>-</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>-</mo>
<mn>8</mn>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>-</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>-</mo>
<mn>12</mn>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
</mrow>
</msub>
<mo>&prime;</mo>
</msup>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>70</mn>
</mfrac>
<mo>&lsqb;</mo>
<mo>-</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mrow>
<mo>(</mo>
<mi>p</mi>
<mo>-</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>+</mo>
<mn>4</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>-</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mn>6</mn>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>-</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>+</mo>
<mn>69</mn>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
</mrow>
</msub>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<mi>p</mi>
<mo>=</mo>
<mn>3</mn>
<mo>,</mo>
<mn>4</mn>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<mi>m</mi>
<mo>-</mo>
<mn>2</mn>
</mrow>
其中m为每个所述原始沉降序列中数据的数量。
3.如权利要求2所述的方法,其特征在于,步骤对所述初始沉降序列采用Verhulst灰色模型进行预测,得到相应的初步预测沉降序列包括:
对于所述初始沉降序列Xi′进行累加运算,获得累加序列Xi (1)=(xi1 (1),xi2 (1),...,xip (1),...,xim (1)),i=1,2,...,k;p=1,2,...,m,其中
<mrow>
<msup>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mi>p</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msup>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>l</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>p</mi>
</munderover>
<msup>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mi>l</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msup>
<mo>,</mo>
<mi>l</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<mi>p</mi>
<mo>;</mo>
</mrow>
计算所述累加序列Xi (1)的紧邻均值生成数,获得紧邻均值生成序列Zi (1)=(zi2 (1),zi3 (1),...,zip (1),...,zim (1)),i=1,2,...,k;p=1,2,...,m,其中
zip (1)=0.5xip (1)+0.5xi(p-1) (1);
建立所述累加序列Xi (1)和紧邻均值生成序列Zi (1)之间的灰色Verhulst微分方程:
xip (1)+azip (1)=b(zip (1))2,i=1,2,...,k;p=1,2,...,m
式中,a为发展系数,b为灰色系数,求解所述微分方程得到所述微分方程的时间响应序列:
<mrow>
<msup>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>i</mi>
<mrow>
<mo>(</mo>
<mi>p</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msup>
<mo>=</mo>
<mfrac>
<mrow>
<mi>a</mi>
<msup>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>i</mi>
<mn>0</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msup>
</mrow>
<mrow>
<mi>b</mi>
<msup>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>i</mi>
<mn>0</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msup>
<mo>+</mo>
<mrow>
<mo>(</mo>
<mi>a</mi>
<mo>-</mo>
<mi>b</mi>
<msup>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>i</mi>
<mn>0</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msup>
<mo>)</mo>
</mrow>
<msup>
<mi>e</mi>
<mrow>
<mi>a</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mi>p</mi>
<mo>)</mo>
</mrow>
</mrow>
</msup>
</mrow>
</mfrac>
<mo>,</mo>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<mi>k</mi>
<mo>;</mo>
<mi>p</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<mi>m</mi>
</mrow>
其中,所述发展系数a和灰色系数b通过最小二乘法确定,即
<mrow>
<mover>
<mi>a</mi>
<mo>^</mo>
</mover>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>a</mi>
<mo>,</mo>
<mi>b</mi>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<msup>
<mi>B</mi>
<mi>T</mi>
</msup>
<mi>B</mi>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<msup>
<mi>B</mi>
<mi>T</mi>
</msup>
<mi>Y</mi>
</mrow>
式中,为灰色参数的估计向量,B为自变量矩阵,Y为估计预测值向量,其表达式分别为:
<mrow>
<mi>B</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msup>
<msub>
<mi>z</mi>
<mrow>
<mi>i</mi>
<mn>2</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msup>
</mrow>
</mtd>
<mtd>
<msup>
<mrow>
<mo>(</mo>
<msup>
<msub>
<mi>z</mi>
<mrow>
<mi>i</mi>
<mn>2</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msup>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msup>
<msub>
<mi>z</mi>
<mrow>
<mi>i</mi>
<mn>3</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msup>
</mrow>
</mtd>
<mtd>
<msup>
<mrow>
<mo>(</mo>
<msup>
<msub>
<mi>z</mi>
<mrow>
<mi>i</mi>
<mn>3</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msup>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msup>
<msub>
<mi>z</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msup>
</mrow>
</mtd>
<mtd>
<msup>
<mrow>
<mo>(</mo>
<msup>
<msub>
<mi>z</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msup>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<mi>Y</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msup>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>&prime;</mo>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mn>3</mn>
</mrow>
</msub>
<mo>&prime;</mo>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mi>m</mi>
</mrow>
</msub>
<mo>&prime;</mo>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>;</mo>
</mrow>
将所述时间响应序列进行累减运算,获得所述初始沉降序列Xi′的初步预测沉降序列其中
4.如权利要求3所述的方法,其特征在于,步骤对所述初步预测沉降序列进行归一化处理,得到所述初步预测沉降序列归一化后的初步预测沉降序列包括:
设为中的最大元素,为中的最小元素,则:
<mrow>
<msub>
<mover>
<mi>x</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>i</mi>
<mi>p</mi>
</mrow>
</msub>
<mo>=</mo>
<mn>0.1</mn>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>i</mi>
<mi>p</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>i</mi>
<mi>min</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>i</mi>
<mi>max</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>i</mi>
<mi>min</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>&times;</mo>
<mrow>
<mo>(</mo>
<mn>0.9</mn>
<mo>-</mo>
<mn>0.1</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<mi>k</mi>
<mo>;</mo>
<mi>p</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<mi>m</mi>
</mrow>
其中,为所述归一化后的初步预测沉降序列中的元素;
步骤将所述沉降预测序列中的元素进行反归一化操作,获得反归一化后的沉降预测序列包括:
设为所述沉降预测序列中的最大元素,为中的最小元素,则:
<mrow>
<msub>
<mover>
<mi>y</mi>
<mo>~</mo>
</mover>
<mi>j</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mi>j</mi>
</msub>
<mo>-</mo>
<mn>0.1</mn>
<mo>)</mo>
<mo>(</mo>
<msub>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>j</mi>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>j</mi>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>0.8</mn>
</mfrac>
<mo>,</mo>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<mi>r</mi>
</mrow>
其中,为所述反归一化后的沉降预测序列中的元素,r为聚类半径。
5.如权利要求1所述的方法,其特征在于,在计算所述初步预测沉降序列时,每预测一个初步预测沉降值后,在所述初始沉降序列Xi′中添加一个新的初始沉降值,并将所述初始沉降序列Xi′中的首个元素去除,保持所述初始沉降序列Xi′的序列长度不变,再次进行预测。
6.一种古建筑沉降的组合预测系统,其特征在于,所述系统包括:
预处理模块,用于对应变片采集的k组原始沉降序列进行五点三次平滑滤波,将所述原始沉降序列中的噪声滤除,得到相应的初始沉降序列;
灰色模型初步预测模块,用于对所述初始沉降序列采用Verhulst灰色模型进行预测,得到相应的初步预测沉降序列;
归一化模块,用于对所述初步预测沉降序列进行归一化处理,得到所述初步预测沉降序列归一化后的初步预测沉降序列;
沉降预测模块,用于采用基于蚁群聚类算法的RBF神经网络对所述归一化后的初步预测沉降序列进行预测,获得沉降预测序列;
反归一化模块,用于将所述沉降预测序列中的元素进行反归一化操作,获得反归一化后的沉降预测序列。
7.如权利要求6所述的系统,其特征在于,所述灰色模型初步预测模块包括:
累加子模块,用于对于所述初始沉降序列进行累加运算,获得累加序列;
紧邻均值计算子模块,计算所述累加序列的紧邻均值生成数,获得紧邻均值生成序列;
灰色微分方程建立子模块,用于建立所述累加序列和紧邻均值生成序列之间的灰色Verhulst微分方程,求解所述微分方程得到所述微分方程的时间响应序列;
序列估计子模块,用于对所述时间响应序列进行累减运算,获得所述初始沉降序列的初步预测沉降序列。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710654232.2A CN107368928B (zh) | 2017-08-03 | 2017-08-03 | 一种古建筑沉降的组合预测方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710654232.2A CN107368928B (zh) | 2017-08-03 | 2017-08-03 | 一种古建筑沉降的组合预测方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107368928A true CN107368928A (zh) | 2017-11-21 |
CN107368928B CN107368928B (zh) | 2021-05-04 |
Family
ID=60309143
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710654232.2A Expired - Fee Related CN107368928B (zh) | 2017-08-03 | 2017-08-03 | 一种古建筑沉降的组合预测方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107368928B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108304649A (zh) * | 2018-01-31 | 2018-07-20 | 安庆师范大学 | 一种高层建筑形变预测方法 |
CN109492793A (zh) * | 2018-09-29 | 2019-03-19 | 桂林电子科技大学 | 一种动态灰色费尔哈斯特神经网络滑坡形变预测方法 |
CN111473768A (zh) * | 2020-03-20 | 2020-07-31 | 淮阴工学院 | 一种建筑物安全智能检测系统 |
CN115140685A (zh) * | 2022-05-16 | 2022-10-04 | 安徽合力股份有限公司 | 一种叉车数据驱动稳定控制方法 |
CN115358475A (zh) * | 2022-08-29 | 2022-11-18 | 河南农业大学 | 基于支持向量机和灰色bp神经网络的灾害预测方法和系统 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102425148A (zh) * | 2011-09-02 | 2012-04-25 | 铁道第三勘察设计院集团有限公司 | 基于静力触探和bp神经网络的路基沉降快速预测方法 |
CN104102853A (zh) * | 2014-08-08 | 2014-10-15 | 武汉理工大学 | 一种利用灰色原理改进的边坡位移分形预测方法 |
CN105069763A (zh) * | 2015-09-17 | 2015-11-18 | 天津理工大学 | 基于五点三次平滑算法的荧光光谱数据噪声滤波方法 |
CN105809249A (zh) * | 2016-03-09 | 2016-07-27 | 浙江工业大学 | 一种基于双神经网络的pm2.5浓度检测与预测系统及方法 |
CN106295873A (zh) * | 2016-08-08 | 2017-01-04 | 西安科技大学 | 一种面向采煤机记忆截割的灰色马尔可夫链轨迹预测方法 |
CN106295869A (zh) * | 2016-08-04 | 2017-01-04 | 辽宁工程技术大学 | 一种基于改进无偏灰色模型的建筑物沉降预测方法 |
CN106407593A (zh) * | 2016-10-10 | 2017-02-15 | 哈尔滨理工大学 | 一种基于gep的高速铁路路基沉降预测系统 |
BR102015009815A2 (pt) * | 2015-04-30 | 2017-03-21 | Inst Fed De Educação Ciencia E Tecnologia Da Paraíba - Ifpb | método de estimativa e previsão de movimento em séries temporais de imagens |
CN106529185A (zh) * | 2016-11-24 | 2017-03-22 | 西安科技大学 | 一种古建筑位移的组合预测方法及系统 |
US20170177992A1 (en) * | 2014-04-24 | 2017-06-22 | Conocophillips Company | Growth functions for modeling oil production |
-
2017
- 2017-08-03 CN CN201710654232.2A patent/CN107368928B/zh not_active Expired - Fee Related
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102425148A (zh) * | 2011-09-02 | 2012-04-25 | 铁道第三勘察设计院集团有限公司 | 基于静力触探和bp神经网络的路基沉降快速预测方法 |
US20170177992A1 (en) * | 2014-04-24 | 2017-06-22 | Conocophillips Company | Growth functions for modeling oil production |
CN104102853A (zh) * | 2014-08-08 | 2014-10-15 | 武汉理工大学 | 一种利用灰色原理改进的边坡位移分形预测方法 |
BR102015009815A2 (pt) * | 2015-04-30 | 2017-03-21 | Inst Fed De Educação Ciencia E Tecnologia Da Paraíba - Ifpb | método de estimativa e previsão de movimento em séries temporais de imagens |
CN105069763A (zh) * | 2015-09-17 | 2015-11-18 | 天津理工大学 | 基于五点三次平滑算法的荧光光谱数据噪声滤波方法 |
CN105809249A (zh) * | 2016-03-09 | 2016-07-27 | 浙江工业大学 | 一种基于双神经网络的pm2.5浓度检测与预测系统及方法 |
CN106295869A (zh) * | 2016-08-04 | 2017-01-04 | 辽宁工程技术大学 | 一种基于改进无偏灰色模型的建筑物沉降预测方法 |
CN106295873A (zh) * | 2016-08-08 | 2017-01-04 | 西安科技大学 | 一种面向采煤机记忆截割的灰色马尔可夫链轨迹预测方法 |
CN106407593A (zh) * | 2016-10-10 | 2017-02-15 | 哈尔滨理工大学 | 一种基于gep的高速铁路路基沉降预测系统 |
CN106529185A (zh) * | 2016-11-24 | 2017-03-22 | 西安科技大学 | 一种古建筑位移的组合预测方法及系统 |
Non-Patent Citations (7)
Title |
---|
WANG JIN-SONG等: "Battlefield Situation Analysis and Prediction Based on Grey Neural Network", 《ELECTRONICS OPTICS & CONTROL》 * |
ZHANG XIAOHONG等: "Research on anti-collision algorithm of block ALOHA in static environment", 《COMPUTER ENGINEERING AND APPLICATIONS》 * |
张小红等: "文物裂隙趋势预测模型研究", 《文物保护与考古科学》 * |
李中凯: "《产品平台设计与产品族开发》", 31 July 2013 * |
杨根新、杨建文: "灰色Verhulst模型在高层建筑沉降预测中的应用", 《勘察科学技术》 * |
王美霞: "高速铁路软土路基沉降预测研究", 《中国优秀硕士学位论文全文数据库 工程科技II辑》 * |
高红: "支持向量机优化组合模型在基坑变形预测中的应用研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108304649A (zh) * | 2018-01-31 | 2018-07-20 | 安庆师范大学 | 一种高层建筑形变预测方法 |
CN108304649B (zh) * | 2018-01-31 | 2020-11-20 | 安庆师范大学 | 一种高层建筑形变预测方法 |
CN109492793A (zh) * | 2018-09-29 | 2019-03-19 | 桂林电子科技大学 | 一种动态灰色费尔哈斯特神经网络滑坡形变预测方法 |
CN111473768A (zh) * | 2020-03-20 | 2020-07-31 | 淮阴工学院 | 一种建筑物安全智能检测系统 |
CN111473768B (zh) * | 2020-03-20 | 2021-10-29 | 淮阴工学院 | 一种建筑物安全智能检测系统 |
CN115140685A (zh) * | 2022-05-16 | 2022-10-04 | 安徽合力股份有限公司 | 一种叉车数据驱动稳定控制方法 |
CN115140685B (zh) * | 2022-05-16 | 2024-02-13 | 安徽合力股份有限公司 | 一种叉车数据驱动稳定控制方法 |
CN115358475A (zh) * | 2022-08-29 | 2022-11-18 | 河南农业大学 | 基于支持向量机和灰色bp神经网络的灾害预测方法和系统 |
Also Published As
Publication number | Publication date |
---|---|
CN107368928B (zh) | 2021-05-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107368928A (zh) | 一种古建筑沉降的组合预测方法及系统 | |
Vazirizade et al. | Seismic reliability assessment of structures using artificial neural network | |
Xu et al. | Improving performance of retaining walls under dynamic conditions developing an optimized ANN based on ant colony optimization technique | |
CN114707227A (zh) | 一种基于数字孪生的大坝安全预警消警方法及系统 | |
Tesfamariam et al. | Risk-based seismic evaluation of reinforced concrete buildings | |
CN105260786B (zh) | 一种电力推进系统仿真可信度评估模型综合优化方法 | |
Jan et al. | Neural network forecast model in deep excavation | |
Kazemi et al. | Predicting seismic response of SMRFs founded on different soil types using machine learning techniques | |
CN106529820A (zh) | 一种运营指标的预测方法及系统 | |
Jalayer et al. | Effects of two alternative representations of ground‐motion uncertainty on probabilistic seismic demand assessment of structures | |
Mahmoudi et al. | Predicting the earthquake magnitude using the multilayer perceptron neural network with two hidden layers | |
CN106529185B (zh) | 一种古建筑位移的组合预测方法及系统 | |
Fathnejat et al. | A data-driven structural damage identification approach using deep convolutional-attention-recurrent neural architecture under temperature variations | |
Fei et al. | Research on tunnel engineering monitoring technology based on BPNN neural network and MARS machine learning regression algorithm | |
Nikose et al. | Dynamic wind response of tall buildings using artificial neural network | |
CN108416458A (zh) | 一种基于bp神经网络的隧道富水岩体综合超前地质预报方法 | |
CN106096115A (zh) | 一种基于自组织沙堆模型的人群疏散仿真方法 | |
CN111709072A (zh) | 一种地下地震动幅值参数预测方法 | |
Kaveh et al. | Feasibility of pso-anfis-pso and ga-anfis-ga models in prediction of peak ground acceleration | |
Karimi Ghaleh Jough et al. | Uncertainty analysis through development of seismic fragility curve for an SMRF structure using an adaptive neuro-fuzzy inference system based on fuzzy C-means algorithm | |
Noureldin et al. | Machine learning-based seismic assessment of framed structures with soil-structure interaction | |
Fayaz et al. | Assessment of ground motion amplitude scaling using interpretable Gaussian process regression: Application to steel moment frames | |
Aloisio et al. | Damage classification after the 2009 L’Aquila earthquake using multinomial logistic regression and neural networks | |
CN112666605B (zh) | 基于主成分分析和多目标遗传算法挑选地震动的方法 | |
KR102589525B1 (ko) | 지진 취약 정보 제공 장치 및 방법, 지진 취약 정보 제공을 위한 학습 모델의 학습 방법 |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210504 Termination date: 20210803 |
|
CF01 | Termination of patent right due to non-payment of annual fee |