线性插值技术在水情工作中的应用

1、线性插值技术在水情工作中的应用摘要：本文以昭平台水库水位库容查询表格、昭平台水库降雨径 流查询表格的制作方法为例，阐述一下线性插值在水情工作中的应 用。关键词：水情工作；线性插值技术；电子表格abstract: taking the production methods of the zhaopingtai reservoir water storage query form and reservoir runoff query forms as an example, this paper explains to the application of the linear interpola

2、tion in the hydrological workkey words: water regimen; linear interpolation technology; electronic form中图分类号：tv211.2文献标识码：a文章编号:2095-2104(2012) exclc电子表格在我们的日常学习、工作中得到了广泛应用，它可 以输入输出、显示数据，利用其提供的丰富的函数，可以进行复杂 繁琐的数据计算，并能把运算后的结果显示为我们需要的图表，大 大提高了我们的工作效率。1昭平台水库水位库容查询表格的制作方法水位库容曲线是一条单一线，需要用一维插值技术来解决，什么 是一维插

3、值呢？打开excel,我们首先建立一个工作薄，取名“水位库容查算”， 在工作薄里新建两个工作表，一个取名“水位库容查算冬一个取 名"zw”。打开“zw”表，从已知的昭平台水库的水位库容关系线摘录节 点数据，输入到表格中，为了计算的方便，要按照两列输入，一列 水位数据，一列库容数据。在f2单元格(也可以是别的单元格) 中输入“z”,在g2单元格输入“w”,把f3单元格链接到“水位库 容查算表的a3单元格，也就是对输入水位数据的一个引用，在g3 单元格中输入公式"二if(f3二offset($a$4, last,), offset($b$4, last,), trend(off

4、s ct($b$4, match(f3, offset($a$4, 1, last), 2), offset($a$4, ma tch (f3, offset ($a$4, 1, , last), , 2), f3) 公式的意义就是:通 过判定给定数据在数据表中的位置，通过线性拟合，返回一个我们 需要的值，这样查算表格就算做好了。实际应用时，在“水位库容查算”表屮输入值，库容w值自 动显示，十分快捷方便。还有一个方法就是：我们新建一个工作表，在工作表的a列与b 列分别输入库容、水位节点数据，在c7单元格(也可以是别的单 元格)我们输入公式“二if($b$2二b8, 0, b7*($b$2-b8

5、)*($b$2-b9)/(b7-b8)/(b7-b9)+b 8* ($b$2-b7)* ($b$2-b9)/(b8-b7)/(b8b9)+b9*($b$2-b7)*($b$2- b8)/(b9-b7)/(b9-b8)n；其实公式的意义就是如果判定给定的数 据不在区间，则返回值“0”，如果在给定的区间则利用插值计算公 式进行计算,在b3单元格中输入“=max(c7:c3156) ”,就是找出计 算出数据列的最大值，因为不在给定区间返冋的都是“0”值，这 样就可以得到给定水位值的库容值。这种方法其实是利用拉格朗日 线性插值公式，由于是利用直线去模拟曲线，所以节点数据一点要 多。应用时，在b2单元格

6、中输入水位数据，在b3中显示库容值。比较以上两种方法，第二种有些繁琐，运行效率较低，并且很难 实现同时对多个值进行查询，所以只是作为一种方法来探讨，实际 工作中，还是选用笫种方法，更为简便快捷。2、昭平台水库降雨径流关系查询表格过程打开excel,我们首先建立一个工作薄，取名“径流查算系统”， 在工作薄里新建两个工作表，一个取名“径流查算”；一个取名“p pa r”。打开“ppa严表，从已知的昭平台水库的降雨径流关系线摘 录节点数据，输入到表格中。其实这是一个二维线性插值的问题，什么是二维线性插值呢？对 于一个二维数组r,我们假设对于任意一个浮点数i,r(a, i)到 r (a+1, i)之间

7、是线性变化的,r (i, b)到r (i, b+1)之间也是线性变化 的(a, b都是整数),那么对于浮点数的坐标(x, y)满足(ar(x, b) = ra+l*(x-s)+ 讥s*(l+3-x);r(x, b+1)=讥a+1 b+l*(x-a)+ra b+l*(l+a-x);好，现在已经知道r(x, b)和r(x,b+l) 了,而根据假设r(x, b)到r(x, b+1)也是线性变化的，所以：r (x, y)二 r (x, b+l)*(y-b)+r (x, b)*(l+b-y)这就是二维线性插值，在数学上，二维线性插值是有两个变量的 插值函数的线性插值扩展，其核心思想是在两个方向分别进行一

8、次 线性插值。在人工查询降雨径流关系线的时候，我们其实是在图上 判断给定条件降雨量p、前期影响雨量w在图上的位置，用二维线 性插值进行模拟计算，返回我们需要的径流深r。根据模拟曲线的需要，我们多摘录一些数据，摘录的时候，一定 要再三核对，确保数据正确，因为这牵涉到以后查算结果的准确性。 在“i3”里输入“p”,在“i4”里输入“pa”,在“i5”里输入“严, 节点数据录入完后，在中输入“=c7”，在j4中输入c8”(当然也可以是别的单元格)，在“j5”单元格里输入以下公式jtrend (if (j4=lasty, offset (a7, wy, wx-(j3=lastx), 2), choos

9、 e (1, 2, trend (offset (a7, wy, wx, 2), offset (a7, wy, 2), j4), tre nd (offset (a7, wy, wx+ (j3olastx), 2), offset (a7, wy, , 2), j4), offset (a7, wx-(j3=lastx), 2), j3) ”。这里用到几个函数,trend 函数是返回线性回归拟合线的纵坐标，offset函数的功能为以指定 的引用为参照系，通过给定偏移量得到新的引用；choose函数从参 数列表中选择并返回一个值;还有"lastx> lasty> wy&g

10、t; wx”儿个 自定义的名称。我们先判定给定的查询条件“p”与“pa”值在节 点数据表格中的位置，然后利用插值公式进行计算。“径流查算表” 中的r值是链接了 “ppar”表中的数据，这是为了美观和节点 数据安全的需要，当然，为了进一步提高安全性，我们还可以为工 作表加密码，这样，昭平台水库降雨径流查询系统就完成了。实际应用时，在“径流查算”表屮输入与“pa”值，径流深 值自动显示。本查询系统利用excel电子表格开发，技术成熟，适用范围广， 拓展性好。洪水预报方案中有许多关系线，比如水位一-库容关系 线、径流深一峰量关系线等，都可以用该技术去进行插值计算。完 全可以利用该技术开发出开发出一个功能更强大的洪水预报软件， 彻底抛开原来的方案图表，比如预报入库的水量，水位的涨幅、洪 峰流量等等，大大提高了工作效率。3.结束语利