matlab 水文断面,用MATLAB处理相关分析在水文计算中的应用

第 20 卷第 4 期 2014 年 4 月 水利科技与经济 Water Conservancy Science and Technology and Economy Vol. 20 No. 4 Apr. ，2014 ［收稿日期］ 2014 －03 －03 ［作者简介］ 崔琳琳( 1981 － ) ，女，黑龙江尚志人，工程师 . 用 MATLAB 处理相关分析在水文计算中的应用 崔琳琳 ( 哈尔滨市水务科学研究院，哈尔滨 150001) ［摘 要］ 生产实践中的许多实际问题需要对数据进行统计分析，建立适当的数学模型，过去通常使用 SAS、SPSS 软件进行分析。利用 MATLAB 软件处理相关分析来解决水文计算中的有关问题，如单位线的推求、相关分析及系列插补延长的计算等。实践证明，该方法计算准确、快捷，有一定的实用价值。 ［关键词］ MATLAB 软件; 单位线; 相关分析; 系列的插补延长 ［中图分类号］ TV122 ［文献标识码］ B ［文章编号］ 1006 －7175( 2014) 04 －0106 －03 0 引 言 水文计算是应用水文学的一个分支，是为在自然水体上建设各类工程如水库、桥涵、引水灌溉、堤防、运河等的工程规划、设计及管理运用提供所需水文依据的计算工作，是水利工程规划的一个前提条件［1］。水文系统一般包括动态物理规律和系统统计规律，在水文计算中对这些现象的研究通常采用统计方法，如水文频率计算、相关分析等［2］。由于 MATLAB 软件能够提供可靠的数值计算、图像及表格显示与处理高水平的图形界面功能［3］，将其应用于水文计算中，不仅可以方便地解决实际计算中的问题，而且具有运算速度快、精度高的优点，避免了传统的水文计算方法工作量大，速度慢，且精度低的不足［4］。 1 单位线的推求 单位过程线法是 L·K 谢尔曼 1932 年提出来的，至今仍是工程水文中最常用的方法。因为它是由汇流断面处的实测水文资料经过反复的验证而求来的，因此它又可以称为是经验单位线法。 1. 1 单位线的基本概念 流域上单位时段内均匀分布的单位地面净雨，汇流到流域出口断面处所形成的地面径流过程线称为单位线。单位线一般有 3 个基本要素: ①单位线洪峰流量; ② 洪峰滞时; ③单位线总历时。 在计算过程中，单位净雨深一般取为10 mm，单位时段可取 1、3、6 h等等，具体要根据流域的特征和精度的要求来确定。 1. 2 最小二乘法推求单位线 单位线通常可以根据实际的降雨量来推求，常用的 方法有分析法、试错法等。二乘法的基本原则是: 令实测的地面流量为 Q，由单位线推算的流量为 Q'。记误差 ε = Q － Q'，若取∑ε = ( Q － Q') = min，流域出口断面流量过程线的表达式为: Qi = ∑ m j =1 hj 10 qi － j +1 i =1，2，…，l j =1，2，…，m i － j +1 =1，2，…，{ n ( 1) 式中: Q 为流域出口断面各时刻流量值，m3 /s; hj 为各时段净雨量，mm; qi － j +1为单位线各个时间点的函数值; l 为流域出口断面流量过程线段数; m 为净雨历时段数; n 为单位线历时段数。该方程转变为数量相等的正规式，然后在最小二乘法求解思想下，解出 q 最优估计值。将式( 1) 写成矩阵的形式: h·q = Q ( 2) 式中: h = h1  h        m ，q = q1 q2   q             n ，Q = Q1 Q2  