原创力文档- 1、本文档共5页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
.............
根据下表统计结构计算weibull函数三参数，并求解设计波高。
有义波高
下限
0
0.5
1
1.5
2
2.5
3
3.5
4
4.5
5
上限
0.49
0.99
1.49
1.99
2.49
2.99
3.49
3.99
4.49
4.99
5.49
数量
1639
3216
1949
1247
786
576
326
224
40
4
2
Total
10000
一、计算原理
海上定点波浪观测短期子样的统计特征，有义波高的长期累计子样为：
大量观测分析表明，子样代表的随机过程仍然是一个平稳的随机过程，可以寻求适当的概率密度函数来拟合观测结果，而三参数的 Weibull 函数就能用来恰当的拟合观测结果。
为 三参数的 Weibull 函数。其中 H0 为最小阈限水平，HC 为尺度因子，
相应的累计概率函数为：
对 Weibull 累计概率函数线性化，移项，取对数：
移动负号，再取对数：
做变量置换：
则有：
为截距 B 和斜率ξ的直线方程。
给定子样：
考虑到直线方程只能解两个未知数，对第三个未知数必须假定。通常，取 H0=0 (第一次近似)。可以用作图法和最小二乘法计算得到另外两个参数：
B和ξ则：
同时， 可以得到本次近似计算的拟合误差平方和：
迭代计算：在一系列 H0的假定下，重复上述计算，得到相应的拟合误差平方和集(子样)
由此得到其最小值σmin和对应的 H0，再重复一次计算，
最后，得到关于三参数的最优解：
将三参数最优解带入相应的公式可以求出设计波波高HD。
二、计算程序及程序说明
有义波高
下限
0
0.5
1
1.5
2
2.5
3
3.5
4
4.5
5
上限
0.49
0.99
1.49
1.99
2.49
2.99
3.49
3.99
4.49
4.99
5.49
数量
1639
3216
1949
1247
786
576
326
224
40
4
2
Total
10000
计算过程说明：
我们取有义波高HS的在每段范围内的平均值，如HS在（0—0.49）范围内我们取HS等于0.245 ，算出p(HS) ,依次求出各个有义波高下的概率，然后累计求和，求出各个波高下的累计概率P（HS）。
然后代入公式（先假设H0=0）
求出对应得不同有义波高下得Xi，Yi的，然后利用EXCEL中最小二乘法拟合公式求：
斜率ξ：=INDEX(LINEST(known_ys,known_xs),1)
Y 轴截距B：=INDEX(LINEST(known_ys,known_xs),2)
再利用公式 ，带入原来的X值， 求出新的Y′然后求出H0下得各点误差σi，进而求出相应的误差平方和Σσi2
改变H0 ，用上诉相同的方法求出不同H0下对应的误差平方和Σσi2，找出一个误差平方和Σσi2 最小下得H0，再代入上述公式求出关于三参数的最优解
[ H0 , Hc , ξ]σmin
代入下列公式
求出设计波波高HD= 9.28168m。EXCEL列表如下（以H 0=0 为例）。
H0=0时的线性化Weibull函数图
在不同H0下的误差平方和Σσi2
三、计算结果
由EXCEL表格我们可以求出当H0=0时 ，σmin最小，这时候所得出的三参数最优解是
H0=0 ,HC= 1.047735 ,ξ= 1.363397，代入下列公式
求出设计波波高HS= 9.28168m。
文档评论（0）