曲面搜索法推求暴雨强度公式参数.docx
- 文档编号:26124650
- 上传时间:2023-06-17
- 格式:DOCX
- 页数:8
- 大小:51.70KB
曲面搜索法推求暴雨强度公式参数.docx
《曲面搜索法推求暴雨强度公式参数.docx》由会员分享,可在线阅读,更多相关《曲面搜索法推求暴雨强度公式参数.docx(8页珍藏版)》请在冰豆网上搜索。
曲面搜索法推求暴雨强度公式参数
2006年11月第22卷第6期沈阳建筑大学学报(自然科学版
JournalofShenyangJianzhuUniversity(NaturalScienceNov.2006Vol122,No16
文章编号:
1671-2021(200606-0959-03
曲面搜索法推求暴雨强度公式参数
傅金祥1,林齐1,王俊凯2
(1.沈阳建筑大学市政与环境工程学院,辽宁沈阳110168;2.大连理工大学数学系,辽宁大连116023
摘要:
目的为解决暴雨强度公式参数在传统推求方法中所存在的非线性函数拟合的问题.方法应用黄金分割法将参数b和n的值固定,即将非线性函数进行线性化,并在由这两个参数所确定的曲面上,通过计算机的特性来搜索参数A1和C的最优值.结果通过对线性函数的拟合,可以依次得到4个参数的值,近而推得暴雨强度公式.结论经过对多个城市的暴雨强度公式推求与比较,曲面搜索法能够快速、精确地得出计算结果,有效地解决了传统推求方法存在的诸多弊端,计算误差远远小于国家规范要求.
关键词:
暴雨强度公式;曲面搜索法;黄金分割;最小二乘法中图分类号:
X143文献标识码:
A
目前我国各城市应用的暴雨强度公式为超定非线性方程[1],参数的求解过程属于非线性函数的优化问题.传统推求公式的方法一般存在着难以直接优化,或计算较为烦琐、拟合误差较大等诸多问题.现将曲面搜索法应用到暴雨强度公式的推求中.此法可以有效地对方程进行拟合并解得最优值,在对传统推求方法改进的基础上,大大提高了公式推求的速度和精度.使用VisualBasic610计算机语言编制了曲面搜索法推求暴雨强度公式参数的程序,计算机程序的应用大大提高了参数推求的速度,并且可以避免人为的计算错误,其推求结果远远小于国家规范所要求的0105mm#min
-1
n.则暴雨强度公式变为线性函数[3]
i=f(A1,C,b,n=
fA1(b,n(1+fC(b,nlgP
(t+b
g(b,n
=
(1
式中:
降雨强度i、重现期P以及降雨历时t均为已知参数.A1为重现期1年的设计降雨的雨力,C为雨力变动参数,b为历时附加参数,n历时参数函数g(b,n即为参数b、n构成的二维曲面函数,使用计算机程序,根据拟合式
m
EF=
的绝对误差和5%的相对误差
[2]
k=1
6[g(b,n-
ik]
2
(2
.
1曲面搜索法基本原理
曲面搜索法是在由已知参数所确定的二维曲
面上,将曲面无限网格化,并对未知参数进行全局搜索,根据误差最小原则,得出参数的最优值.暴雨强度公式为四参数超定非线性方程,对未知参数进行拟合,可将参数A1和C分别看作是参数b、n的函数,并分别表示为fA1(b,n和fC(b,
在曲面上对参数进行全局搜索,便可以将参数b、
n的最优值解出.式中m为降雨历时的个数.在搜索过程中,还使用了黄金分割法,不断的大范围压缩搜索空间,以达到使未知参数快速接近最优值的目的.拟合得最优参数值后,便可以通过最小二乘法分别计算求得参数A1和C的值[4].
由于传统的曲面最小二乘法所使用的推导基础曲面,是由重现期P和降雨历时t所构成的二元空间曲面[5],所以这种方法只能推导出暴雨强
收稿日期:
2006-03-07
基金项目:
辽宁省自然基金项目(20072005
作者简介:
傅金祥(1955-,男,教授,博士后,主要从事水处理理论与技术研究.
960
沈阳建筑大学学报(自然科学版第22卷
度公式的总公式.如果所提供的重现期只有一个,则公式推导过程中就会出现计算错误.而在曲面搜索法中,寻优基础是由降雨历时修正参数b和降雨衰减指数n所构成的曲面,这种方法不但可以推导出总公式,还可以推导出分公式[6].
雨强度公式的4个参数值.曲面搜索法推求暴雨强度公式程序流程图如图1所示
.
2曲面搜索法推求步骤
为了提高推求速度和计算精度,笔者使用计算机VisualBasic610语言编制了曲面搜索法推求暴雨强度公式的程序.现就其推求步骤进行说明.
第一步:
使用计算机程序导入降雨强度值i、降雨历时t和重现期P,即已知的i-t-P表.
第二步:
用黄金分割法在区间[0,20]内搜索降雨历时修正参数b的值;在区间[0,2]内搜索降雨衰减指数n的值[7].
第三步:
使用公式推求出雨力变动参数C的值.
n
C=ilgPi2
6lgP6-6lgP6(t+b(t+b第四步:
使用公式推出重现期为一年的设计
6lgP6
n-
m6
降雨的雨力参数A1的值.
A1=
6
(t+b第五步:
计算公式的拟合误差EF
6
2
EF=k6[g(b,n-ik]=1
m
图1曲面搜索法推求暴雨强度公式程序流程图
3应用与分析
分别使用曲面搜索法和传统使用的最小二乘法推求4个城市的暴雨强度公式,并计算出公式
的绝对误差R绝和相对误差R相并在表1中(4个城市分别定义为A、B、C、D进行对比.其中,R绝12
/m,R相=R绝/i表,(t+b
式中的i表为所收集的降雨强度值[8].=
6(i表-
第六步:
比较误差的大小;如果符合要求,则返回参数的值;否则返回第二步,直到满足误差要求.
至此,便得到了使用曲面搜索法推求出的暴
表1曲面搜索法和最小二乘法计算误差比较
计算法
R绝
最小二乘法曲面搜索法
010*********
A市
R相41692195
R绝010*********
B市
R相31612169
R绝010*********
C市
R相21391175
R绝010*********
D市
R相21431149
注:
表1中绝对均方差R绝的单位为mm#min-1;相对均方差R相的单位为%.
从表1中可以看到,曲面搜索法相对于最小二乘法在推求暴雨强度公式时,无论是绝对误差还是相对误差都有明显地下降[9].4个城市的误
差最大下降了3816%,最小下降了2516%,平均下降26195%.
可见,曲面法在推求暴雨强度公式时,计算精
第22卷傅金祥等:
曲面搜索法推求暴雨强度公式参数
961
度得到了显著地提高.此外,在计算速度,即搜索速度方面也有很大地提高.
[2][3]
4结论
(1曲面搜索法将超定非线性式化为关于未
知参数b、n的线性函数,利用计算机快速的全局搜索能力,在二维空间上寻得最优值[10].所解得的值对目标值具有良好的拟合性,其计算速度相对于传统推求方法也有显著提高.
(2曲面搜索法可以作为一种新的方法来推求暴雨强度公式的参数,经过大量实际计算与应用,其计算速度与精度均高于其他方法,计算误差远远小于国家规范标准.
(3使用VisualBasic610语言编制的计算机程序,能够快速准确地计算出结果,可以节省大量烦琐的推求计算过程,并能避免人工推求出现的错误.可广泛的应用于暴雨强度公式参数的推求中.参考文献:
[1]陈芸.GBJ14-87室外排水设计规范[S].北京:
中
[4][5][6]
[7]
[8][9][10]
国建筑工业出版社,1998.
周玉文,赵洪宾.排水管网理论与计算[M].北京:
中国建筑工业出版社,2000.
LorenteJ.Relationbetweenmaximalrainfallratesfordifferenttimeintervalsinthecourseofastorm[J].AtmosRes,1991(7:
61-66.
金光淡.水文统计原理与方法[M].北京:
中国工业出版社,1994.
GuestPG.Numericalmethodsofcurvefitting[M].Camb:
CambridgeUniversityPress,1999.
刘钧哲,马兴冠,傅金祥,等.皮尔逊III型分布曲线的快速简便算法[J].沈阳建筑工程学院学报,2004,20(1:
60-62.
SrinivasM,PatnaikLM.Adaptiveprobabilitiesofcrossoverandmutationingeneticalgorithm[J].IEEETransactiononSystems,ManandCybernetics,1994,24(4:
28-32.
HollandJH.GeneticAlgorithm[J].ScientificAmer2ican,1992(4:
44-50.王小平.遗传算法理论、应用与软件实现[M].西安:
西安交通大学出版社,2002.
AmrnellV.Reviewofrainfalldataapplicationforde2signandanalysis[J].WaterScienseandTechnology,1994(16:
8-9.
CalculatingRainfallIntensityFormula
ThroughtheMethodofSearchinginCurvedFace
FUJinxiang,LINQi,WANGJunkai
1
1
2
(1.SchoolofMunicipalandEnvironmentalEngineering,ShenyangJianzhuUniversity,Shenyang110168,China;2.Depart2mentofMathematics,DalianUniversityofTechnology,Dalian116023,China
Abstract:
Toresolvetheproblemofnon2linearfunctionduringtheprocessoftraditionalcalculatingrainfallintensityformula.Wesuggestthemethodofsearchingincurvedface.Itfixesuptheparametersbandn,inanotherwords,itletthenon2linearfunctionbethelinearone.Onthefaceofparametersbandn,searchthebestanswersofparametersA1andC.Bycalculatingthelinearfunction,wecangetthefourparametersor2derlyandtherainfallintensityformula.Bycalculatingandcomparingseveralcitiescrainfallintensityformu2la,themethodofsearchingincurvedfacecouldgottheanswersfleetlyandaccurately,anditcouldresolvethedisadvantageoftraditionalcalculatingmethods.Theerrorisfarlessthanthecriterionrequest.
KeyWords:
rainfallintensityformula;methodofsearchingincurvedface;goldsegmentation;leasttwicemultiplication
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 曲面 搜索 推求 暴雨 强度 公式 参数