CN116226604B 一种水文序列总体变异度诊断度量方法、装置及设备 （长安大学）

本发明公开了一种水文序列总体变异度诊为多个IMF序列和一个趋势项序列；判断趋势项项序列的趋势变异度；将多个IMF序列叠加重组断均值归零化处理后的IMFs序列是否存在方差变异点，若存在，则计算均值归零化处理后的变异度和方差变异度计算原始水文序列的综合2采用自分段线性拟合法对所述趋势项序列进行分段根据所述均值变异点两侧子序列均值对所述IMFs序列进行均值3)应用二分法在方差变异点处将Yt分为两段子序列，再按照上述步骤对每段子序列进根据所述趋势变异度、所述均值变异度和所述方差变异度计算所1)向所述原始水文序列x(t)中加入一段同时间尺度且符合正态分布的白噪声信号⑴得到多个IMF序列和一个趋势项序列。3rlη=4.19ln(n)_11.546)对5)中得到的两段子序列分别进行以上操作，检验子序列是否还存在均值变异点，4)I22将所述趋势变异度、所述均值变异度和所述方差变异度相加，得到所述原始水采用自分段线性拟合法对所述趋势项序列进行分段方差变异度计算模块，用于根据所述均值变异点两侧子序列均值3)应用二分法在方差变异点处将Yt分为两段子序列，再按照上述步骤对每段子序列进综合变异度计算模块，用于根据所述趋势变异度、所述均值变异58.一种设备，包括存储器、处理器以及存储在所6[0009]将多个所述IMF序列叠加重组得到IMFs序列，判断所述IMFs序列是否存在均值变[0012]进一步地，所述将所述原始水文序列分解为多个IMF序列和一个趋势项序列，包[0013]1)向所述原始水文序列x(t)中加入一段同时间尺度且符合正态分布的白噪声信782个趋势项序列；9[0069]图2为本发明一种水文序列总体变异度诊断度量方法中总体经验模态分解法示意[0070]图3为本发明实施例中无定河流域丁家沟水文站1960_2010年的年径流序列示意[0071]图4为本发明实施例中原始水文序列应用总体经验模态分解法的分解结果示意[0081]S4、根据所述均值变异点两侧子序列均值对所述IMFs序[0085]1)向所述原始水文序列x(t)中加入一段同时间尺度且符合正态分布的白噪声信[0091]2.2)计算所述新水文序列y(t)的上包络线Ymax(t)和所述新水文序列y(t)的下包络nRn把y(t)分解为多个分量序列：[0104]经验模态分解法(EMD)在复杂的非线性信号处理上具有显著的优势，可以将复杂[0109]1.2)按最小二乘法对所述趋势项序列为F(ti)进行分段回归线拟合，得到一组拟[0112]1.3)计算所述趋势项序列为F(ti)的回归拟合线与序列值之间的偏差平方和S，S变异点；更加具体地说，选择显著性水平α,比较χ2检验量与上侧分位数值的大小，若[0119]在上述实施方式的基础上，作为更加优选的实施方式，关于步骤S3中判断所述2个趋势项序列。[0185]本发明提供的一种水文序列总体变异度诊断度量装置用于实现前述的一种水文方式可见于前文中的一种水文序列总体变异度诊断度量方法的执行所述计算机存储介质存储的程序指令。处理器可以是中央处理单元(CentralProcessingUnit，CPU)，还可以是其他通用处理器、数字信号处理器(DigitalSignal流程图一个流程或多个流程和/或方框图一个方框或多个方框中[0191]这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一