下面以求矩适线法、三点适线法及权函数法为例详细说明目估适线法的具体过程。利用本章第4节的公式计算和Cv值，并假定Cs值，作为三个统计参数的初始值，通过适线来确定统计参数值的方法，就称为求矩适线法。按公式推算所求洪水频率的流量：Q2%=×1839=4119Q1%=×1839=4561Q0.33%=×1839=5241表3-6流量频率计算表续表表3-7频率曲线选配计算表三点适线法。 利用［例3-4］的年最大流量资料，用三点适线法确定其统计参数、Cv和Cs。......
2023-11-20
优化目标最优化算法之目估适线法
【摘要】：图X3.11频率计算适线法框图统计参数的初适值，可用矩法公式计算、Cv，并假定Cs/Cv比值，也可采用三点法、概率权重矩法和双权函数法等。图X3.11频率计算适线法框图调整参数。试用适线法推求该站年降水量的频率曲线，并确定相应于频率为10%、50%、90%的年降水量。图X3.13某站年降水量频率曲线类似地，计算F列的模比系数以及G列ki-1和H列2数据，并对各列求和，其中F列的模比系数的总和应等于n，G列的总和应等于0，据此可进行校核。
目估适线法的要点是：拟定理论频率曲线的线型，以样本经验点为依据，调试理论频率曲线的参数，用目估的方法使理论频率曲线与经验点配合良好。具体步骤如图X3.11所示。进一步说明以下几点。
（1）频率曲线线型一般选用皮尔逊Ⅲ型。
（2）统计参数的初适值，可用矩法公式计算
、Cv，并假定Cs/Cv比值，也可采用三点法、概率权重矩法和双权函数法等。三点法将在工作任务2中介绍，其他方法可参阅有关书籍。
图X3.11 频率计算适线法框图
（2）统计参数的初适值，可用矩法公式计算
、Cv，并假定Cs/Cv比值，也可采用三点法、概率权重矩法和双权函数法等。三点法将在工作任务2中介绍，其他方法可参阅有关书籍。
图X3.11 频率计算适线法框图
（3）调整参数。矩法计算
的抽样误差较小，一般可不作修改，主要调整Cv及Cs。
由以上可见，适线法层次清楚，图像明显，方法灵活，易于操作，它是一种能较好地满足水文频率计算要求的估计方法，在水文计算中广泛采用。此法的实质是以经验分布为基础，去估计总体的分布及统计参数。
【例X3.4】 某站年降水量资料见表X3.3。试用适线法推求该站年降水量的频率曲线，并确定相应于频率为10%、50%、90%的年降水量。
表X3.3 某站年降水量资料 单位：mm
具体步骤如下：
（1）计算经验频率并点绘经验点。利用Excel软件进行计算，建立名为“经验频率与统计参数计算.xls”的工作薄，在其工作表“Sheet1”中建立计算表。
1）输入样本系列并降序排列。在A4单元格中输入“1956”，然后移动鼠标选取A4单元格右下角的填充柄，再按下“Ctrl”键及鼠标左键向下拖动至“1984”，即得A列数据，然后在B列输入逐年的年降水量。在C4单元格输入序号“1”，移动鼠标选取C4单元格右下角的填充柄，再按下“Ctrl”键及鼠标左键向下拖动至序号“29”，即得C列数据，然后复制B列的年降水量数据，并粘贴到D列，在D列数据被选中的基础上，单击工具栏中的“降序排序”按钮，并在“排序警告”对话框中选择“以当前选定区域排序”，则得到递减排列的样本序列，如图X3.12中所示的D列数据。
2）计算经验频率。利用“数学与三角函数”中的“ROUND”函数进行计算，并按指定的位数对数值进行四舍五入。例如，保留1位小数，具体方法是，在E4单元格输入“＝ROUND（C4*100/（29＋1），1）”，然后按“Enter”键，即得到以百分数表示的经验频率“3.3”，接着选中E4单元格，鼠标指针变成黑十字时，按住鼠标左键，向下拖动填充E5～E32单元格，即得到各序号对应的经验频率，即E列数据。
图X3.12 经验频率与统计参数计算
3）根据图X3.12中D列数据与E列数据的对应数值，在频率格纸上点绘经验频率点，如图X3.13所示。
（3）调整参数。矩法计算
的抽样误差较小，一般可不作修改，主要调整Cv及Cs。
由以上可见，适线法层次清楚，图像明显，方法灵活，易于操作，它是一种能较好地满足水文频率计算要求的估计方法，在水文计算中广泛采用。此法的实质是以经验分布为基础，去估计总体的分布及统计参数。
【例X3.4】 某站年降水量资料见表X3.3。试用适线法推求该站年降水量的频率曲线，并确定相应于频率为10%、50%、90%的年降水量。
表X3.3 某站年降水量资料 单位：mm
具体步骤如下：
（1）计算经验频率并点绘经验点。利用Excel软件进行计算，建立名为“经验频率与统计参数计算.xls”的工作薄，在其工作表“Sheet1”中建立计算表。
1）输入样本系列并降序排列。在A4单元格中输入“1956”，然后移动鼠标选取A4单元格右下角的填充柄，再按下“Ctrl”键及鼠标左键向下拖动至“1984”，即得A列数据，然后在B列输入逐年的年降水量。在C4单元格输入序号“1”，移动鼠标选取C4单元格右下角的填充柄，再按下“Ctrl”键及鼠标左键向下拖动至序号“29”，即得C列数据，然后复制B列的年降水量数据，并粘贴到D列，在D列数据被选中的基础上，单击工具栏中的“降序排序”按钮，并在“排序警告”对话框中选择“以当前选定区域排序”，则得到递减排列的样本序列，如图X3.12中所示的D列数据。
2）计算经验频率。利用“数学与三角函数”中的“ROUND”函数进行计算，并按指定的位数对数值进行四舍五入。例如，保留1位小数，具体方法是，在E4单元格输入“＝ROUND（C4*100/（29＋1），1）”，然后按“Enter”键，即得到以百分数表示的经验频率“3.3”，接着选中E4单元格，鼠标指针变成黑十字时，按住鼠标左键，向下拖动填充E5～E32单元格，即得到各序号对应的经验频率，即E列数据。
图X3.12 经验频率与统计参数计算
3）根据图X3.12中D列数据与E列数据的对应数值，在频率格纸上点绘经验频率点，如图X3.13所示。
（2）计算统计参数初适值。利用算术平均值函数“AVERAGE”计算样本均值
。选择要输入平均值函数的单元格D34，然后输入“＝AVERAGE（D4：D32）”，再按下“Enter”键，即得到样本均值532.0mm。也可由第（4）列合计值
＝15427.1，除以样本容量29得到。
（2）计算统计参数初适值。利用算术平均值函数“AVERAGE”计算样本均值
。选择要输入平均值函数的单元格D34，然后输入“＝AVERAGE（D4：D32）”，再按下“Enter”键，即得到样本均值532.0mm。也可由第（4）列合计值
＝15427.1，除以样本容量29得到。
类似地，计算F列的模比系数
以及G列ki-1和H列（ki-1）2数据，并对各列求和，其中F列的模比系数
的总和应等于n，G列的总和
应等于0，据此可进行校核。图X3.12中F列、G列的总和，由于舍入影响，尾数略有误差是允许的。
图X3.13 某站年降水量频率曲线
类似地，计算F列的模比系数
以及G列ki-1和H列（ki-1）2数据，并对各列求和，其中F列的模比系数
的总和应等于n，G列的总和
应等于0，据此可进行校核。图X3.12中F列、G列的总和，由于舍入影响，尾数略有误差是允许的。
图X3.13 某站年降水量频率曲线
由H列的数据和
（ki-1）2＝2.3369，计算年降水量的变差系数：
还可通过计算样本均方差，然后除以样本均值，求得变差系数。方法是选择要输入样本均方差函数的单元格D35，然后输入“＝STDEV（D4：D32）”，再按下“Enter”键，即得到样本均方差153.7mm。在样本均方差和样本均值的基础上，依据式（X3.15），在单元格F35中，输入“＝D35/D34”，再按下“Enter”键，即得到变差系数0.29。
一般由式（X3.15）计算的Cv值偏小，取Cv＝0.32，并假定Cs＝2Cv＝0.64作为初适值。
（3）选配皮尔逊Ⅲ型频率曲线。
1）由参数初适值，查模比系数kp表，计算不同频率p相应的xp值，列入表X3.4中第（2）、（3）栏。在图X3.13中，由表X3.4第（1）、（3）两栏的对应数值点绘理论频率曲线，发现理论频率曲线在上部和下部均偏于经验点的下方，而中间部分略偏于经验点的上方，说明Cs偏小。
2）调整参数，重新配线。增大Cs，取Cs＝3Cv＝0.96，再次计算不同频率p相应的xp值列入表X3.4中第（5）栏。由（1）、（5）两栏的对应数值点绘理论频率曲线，如图X3.13所示，该线与经验点配合较好，将其作为该站采用的年降水量的频率曲线。
（4）推求频率10%、50%、90%的年降水量。由表X3.4第（1）、（5）栏或由采用的年降水量频率曲线，可确定频率10%、50%、90%对应的年降水量分别为x10%＝760.8mm、x50%＝505.4mm和x90%＝340.5mm。
表X3.4 理论频率曲线选配计算表
由H列的数据和
（ki-1）2＝2.3369，计算年降水量的变差系数：
还可通过计算样本均方差，然后除以样本均值，求得变差系数。方法是选择要输入样本均方差函数的单元格D35，然后输入“＝STDEV（D4：D32）”，再按下“Enter”键，即得到样本均方差153.7mm。在样本均方差和样本均值的基础上，依据式（X3.15），在单元格F35中，输入“＝D35/D34”，再按下“Enter”键，即得到变差系数0.29。
一般由式（X3.15）计算的Cv值偏小，取Cv＝0.32，并假定Cs＝2Cv＝0.64作为初适值。
（3）选配皮尔逊Ⅲ型频率曲线。
1）由参数初适值，查模比系数kp表，计算不同频率p相应的xp值，列入表X3.4中第（2）、（3）栏。在图X3.13中，由表X3.4第（1）、（3）两栏的对应数值点绘理论频率曲线，发现理论频率曲线在上部和下部均偏于经验点的下方，而中间部分略偏于经验点的上方，说明Cs偏小。
2）调整参数，重新配线。增大Cs，取Cs＝3Cv＝0.96，再次计算不同频率p相应的xp值列入表X3.4中第（5）栏。由（1）、（5）两栏的对应数值点绘理论频率曲线，如图X3.13所示，该线与经验点配合较好，将其作为该站采用的年降水量的频率曲线。
（4）推求频率10%、50%、90%的年降水量。由表X3.4第（1）、（5）栏或由采用的年降水量频率曲线，可确定频率10%、50%、90%对应的年降水量分别为x10%＝760.8mm、x50%＝505.4mm和x90%＝340.5mm。
表X3.4 理论频率曲线选配计算表
上述用矩法公式计算系列的均值和变差系数，亦可由袖珍计算器的统计功能键进行计算：进入统计计算“STAT”状态后，将样本数据依次输入“M＋”键后，就可通过
键得到样本均值和均方差的数值，然后由
即可算得变差系数Cv值，读者可按说明书或在教师指导下操作。
目前，适线法可通过开发的Excel水文频率计算软件来实现。图X3.14为利用河北省水文水资源勘测局工程技术人员[1]开发的Excel水文频率计算软件，对［例X3.4］资料进行频率计算的适线结果。
图X3.14 利用Excel水文频率计算软件的适线结果
使用水文频率计算软件进行适线，具有方便、规范、减小计算工作量等显著优点。
上述用矩法公式计算系列的均值和变差系数，亦可由袖珍计算器的统计功能键进行计算：进入统计计算“STAT”状态后，将样本数据依次输入“M＋”键后，就可通过
键得到样本均值和均方差的数值，然后由
即可算得变差系数Cv值，读者可按说明书或在教师指导下操作。
目前，适线法可通过开发的Excel水文频率计算软件来实现。图X3.14为利用河北省水文水资源勘测局工程技术人员[1]开发的Excel水文频率计算软件，对［例X3.4］资料进行频率计算的适线结果。
图X3.14 利用Excel水文频率计算软件的适线结果
使用水文频率计算软件进行适线，具有方便、规范、减小计算工作量等显著优点。
有关工程水文及水利计算的文章
- 详细阅读
-
多目标蚁群算法模型优化详细阅读
算法根据决策变量X和式～式的约束条件依次确定任务续传调度标识符和相应的执行时间，从而得到决策变量X的目标函数。图9-27可用时间窗口更新交叉；嵌入；包含；无关确定所有任务的续传调度标识和执行时间后，计算目标函数F=｛f 1，f 2，f 3｝，即可评价此调度方案优劣。......
2023-07-02
-
多目标蚁群算法的设计与优化详细阅读
为消除各类约束导致的大量任务间冲突，获得问题的Pareto解集，本书采用多目标蚁群算法优化任务的调度顺序。由于本章算法为多目标蚁群算法，因此本章算法采取与第3章的优化算法不同的启发策略。本章算法多目标蚁群算法采用自适应策略，在算法起始阶段，启发选择比例参数q 0［式（3-6）］取较大值利于加快收敛速度；在算法搜索后期，选择较小的q 0值可增加种群多样性。......
2023-07-02
-
适线法简介-桥涵水文学详细阅读
优化适线法是在一定的适线准则下，求解与经验点据拟合最优的频率曲线的统计参数的方法。随着计算机的推广普及，带有一定准则的计算机优化适线法也常为许多设计单位使用，优化适线法按不同的适线准则分为三种：离差平方和最小准则、离差绝对值和最小准则、相对离差平方和最小准则。本节只简要介绍离差平方和最小准则，使读者对该法有一个概略的了解，但不作为教学上的要求。离差平方和准则的适线法又称最小二乘估计法。......
2023-11-20
-
混沌粒子群优化算法优化方案探究详细阅读
混沌粒子群优化算法的基本思想是采用混沌序列初始化粒子的位置和速度,先对当前粒子群体中的最优粒子进行混沌寻优,然后把混沌寻优的结果随机替换粒子群体中的一个粒子。,z1N),根据式(4-6)得到N个向量z1,z2,…如果粒子适应度优于全局极值gbest,则将gbest设置为新位置。则混沌粒子群优化算法的流程为:初始化,设置最大允许迭代次数或适应误差限,以及CPSO算法的参数惯性权值和学习因子。,PgD)进行混沌优化:将Pgi(i=1,2,…计算其适应值,得到性能最好的可行解pb。......
2023-06-23
-
量子遗传优化算法的优化方法详细阅读
量子遗传算法是一种将遗传算法和量子计算相结合的概率优化方法，两者相互作用。量子遗传算法是一种将量子比特的概率幅用于染色体编码，用量子门的调整操作来实现染色体更新，以完成进化搜索的方法。量子遗传算法的流程如下：初始化种群Q，随机生成n个用量子比特编码的染色体。......
2023-06-29
-
模拟退火算法优化详细阅读
模拟退火算法是由N Metropolis等学者于1953年最早提出的。模拟退火算法从某一给定的初始温度开始，随着迭代的进行，温度参数不断下降，结合算法概率性的突跳能力在问题的可行解空间中随机寻找目标函数的解，并以Metropolis准则判定是否接受新解，如此迭代进行下去，逐步寻找问题的全局最优解。Metropolis接受准则是模拟退火算法的重要核心思想，Metropolis接受准则也是依据固体物质退火过程的特点而提出的。图11-2模拟退火算法流程图③算法运算产生新解。......
2023-06-28
-
中文分词算法的优化策略详细阅读
一般来说，中文分词在具体的算法实现上分为三种：基于字符串匹配的分词方法、基于理解的分词方法和基于统计的分词方法。在中文搜索引擎中，目前基本上是这三种算法混合使用。2）基于统计的分词方法基于统计的分词方法也叫最大概率分词方法。作为中文分词基础的词库，新词补充和老词删除就是非常重要的工作。比如“测试”在“每台计算机在出厂前都要经过严格的测试”这句话中是典型的动词，而在“软件测试领域”中是一个名词。......
2023-07-02
相关推荐