第1页共13页编号：时间：2021年x月x日书山有路勤为径，学海无涯苦作舟页码：第1页共13页《海洋工程环境学》作业1班级:学号:姓名:1.计算条件给定有义波高统计值如下表所示：有义波高下限00.511.522.533.544.55上限0.490.991.491.992.492.993.493.994.494.995.49数量17653090184912978365853172204352Total10009平均周期.2.计算要求根据所给有义波高统计值和平均周期，用三参数Weibull函数预报此海域短期海况有义波高的长期分布规律，并计算百年一遇的最大波。最后作业应包含三部分:(1)计算原理说明;(2)计算程序说明;(3)计算结果:绘制线性化Weibull函数及其离散值图线.绘制曲线,其中为最小阈限高度试算值,为拟合误差的均方根平均值(标准差).第2页共13页第1页共13页编号：时间：2021年x月x日书山有路勤为径，学海无涯苦作舟页码：第2页共13页给出最佳拟合结果:和拟合(标准差).一、计算原理海上定点波浪观测短期子样的统计特征，有义波高的长期累计子样为：大量观测分析表明，子样代表的随机过程仍然是一个平稳的随机过程，可以寻求适当的概率密度函数来拟合观测结果。为三参数的Weibull函数。其中H0为最小阈限水平，HC为尺度因子，ζ为形状因子对Weibull累计概率函数线性化，移项，取对数：第3页共13页第2页共13页编号：时间：2021年x月x日书山有路勤为径，学海无涯苦作舟页码：第3页共13页移动负号，再取对数：做变量置换：则有：¾为截距B和斜率的直线方程给定子样：第4页共13页第3页共13页编号：时间：2021年x月x日书山有路勤为径，学海无涯苦作舟页码：第4页共13页考虑到直线方程只能解两个未知数，对第三个未知数必须假定。通常，取H0=0(第一次近似)。用作图法得到另外两个参数：和则：同时，可以得到本次近似计算的拟合误差平方和：迭代计算：在一系列H0的假定下，重复上述计算，得到相应的拟合误差平方和集(子样)由此得到其最小值和对应的H0，再重复一次计算，最后，得到关于三参数的最优解：二、计算程序本次计算使用了excel的强大的计算功能，并运用了其自带的画图工具。第5页共13页第4页共13页编号：时间：2021年x月x日书山有路勤为径，学海无涯苦作舟页码：第5页共13页三、计算结果Hmax数量p(Hs)P(Hs)-ln(1-P(Hs))X=lg(Hs-H0)Y=lg(-ln(1-P(Hs)))Y=kX+b误差平方0.4917650.17630.17630.1940-0.3098-0.7122-0.71950.00010.9930900.30870.4850.6636-0.0044-0.1781-0.23920.00371.4918490.18470.66971.10780.17320.04450.04010.00001.9912970.12960.79931.60590.29890.20570.23770.00102.498360.08350.88282.14400.39620.33120.39080.00352.995850.05840.94132.83460.47570.45250.51580.00403.493170.03170.97293.60920.54280.55740.62140.00413.992200.0220.99495.27950.60100.72260.71290.00014.49430.00430.99927.13190.65220.85320.79350.0036第6页共13页第5页共13页编号：时间：2021年x月x日书山有路勤为径，学海无涯苦作舟页码：第6页共13页4.9950.00050.99978.11270.69810.90920.86560.00195.4920.00020.99999.21030.73960.96430.93080.0011总量10009Ho=0k=1.5727b=-0.2323误差平方和0.0232Hmax数量p(Hs)P(Hs)-ln(1-P(Hs))X=lg(Hs-H0)Y=lg(-ln(1-P(Hs)))Y=kX+b误差平方0.4917650.17630.17630.1940-0.3565-0.7122-0.73240.00040.9930900.30870.4850.6636-0.0269-0.1781-0.23250.00301.4918490.18470.66971.10780.15840.04450.04850.00001.9912970.12960.79931.60590.28780.20570.24480.00152.498360.08350.88282.14400.38740.33120.39580.00422.995850.05840.94132.83460.46830.45250.51860.00443.493170.03170.97293.60920.53660.55740.62200.00423.992200.0220.99495.27950.59550.72260.71140.00014.49430.00430.99927.13190.64740.85320.79010.0040第7页共13页第6页共13页编号：时间：2021年x月x日书山有路勤为径，学海无涯苦作舟页码：第7页共13页4.9950.00050.99978.11270.69370.90920.86040.00245.4920.00020.99999.21030.73560.96430.92390.0016总量10009Ho=0.05k=1.5166b=-0.1917误差平方和0.0257Hmax数量p(Hs)P(Hs)-ln(1-P(Hs))X=lg(Hs-H0)Y=lg(-ln(1-P(Hs...
