水轮发电机组盘车数据处理及软件开发
合集下载
相关主题
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
(0)
பைடு நூலகம்
bi bi(0) i
(0)
（8）
这样，确定 bi 的问题转化为确定修正值 i 的问题。为确定 bi 在 bi 附近对f作泰勒级数 展开，并略去 i 的高阶项得：
f ( xk , b1 , b2 ,..., bm ) f k 0
(0)
f k 0 f f 1 k 0 2 .... k 0 m b1 b2 bm
按水导处的净摆度计算有：
J ca D L1 L2
（17）
式中：L1为法兰摆度测量处至水导摆度测量处的距离。 但轴线往往不是一条直线，按式（16）和式（17）计算出的结果是不一样的。根据轴线 曲线程度的不同，单独按式（16）或式（17）进行计算，有时得不到好的效果。 当轴线如图2实线所示时，若以法兰处的摆度值进行处理，水导处的摆度会加大（图2 中虚线）；若以水导处的摆度值进行处理，则法兰处的摆度没有明显改善（图2中点线）。
图 1 摆度曲线与测点关系
3 机组轴线调整
当机组摆度过大时， 可通过刮削推力头底面或绝缘垫或法兰面的方法进行调整。 当机组 轴线是一条直线时， 可按法兰处的净摆度和水导处的净摆度计算刮削推力头底面或绝缘垫的 最大厚度。 按法兰处的净摆度计算有：
J ba D L1
（16）
式中：L1为推力头底面至法兰摆度测量处的距离，D为推力头或绝缘垫外径。
Ff 14.9036 sin( 142.95) 8.125 Fs 17.6918 sin( 150.59) 9.25
(15)
最大摆度方位角分别为232.95和240.59。摆度值和最大摆度方位角的计算相比式（4） 和式（5）更为精准。 计算出摆度曲线和测点值如图1。 由图1可以判断测点数据是否符合正弦规律。 若测点数 据严重偏离正弦规律，说明盘车数据不能真实反映轴线状态，需重新进行盘车。
参考文献
[1] [3] [5] 华北水利水电学院等，水力机组安装与检修，电力工业出版社，1981. 罗清萍 . 电动盘车装置的原理与应用实例剖析［J］. 黑龙江水利科技，2013，41（4）：13-16 . 刘润根，杨虹，李辉，等. 水电站机组智能盘车系统设计与应用[J]. 中国水利水电科学研究院学 报. 2015(01): 74-80. [6］ 李则泉， 陈苗 . 锦屏二级水电站7号机组盘车与轴线调整新工艺 ［J ］ . 水电站机电技术， 2014 （4） . [7］ 王胜超 . 自动盘车装置在周宁水电站的应用［J］. 防爆电机，2013（2）. [8］ 孟安波，徐海波 . 基于遗传算法的水电机组盘车数据处理方法研究［J］. 水力发电，2011，37 （10）：61-62 . [9］ 方建新，谭啸 . 浅析状态监测系统的故障分析与盘车应用［J］ . 水电站机电技术，2009，32（3）： 96-98 . [2] 王胜超. 简析水轮发电机组盘车[J]. 电工文摘. 2012(04): 54-58. [4］ 王胜超 . 自动盘车装置在周宁水电站的应用［J］. 防爆电机，2013，48（2）：46-47 .
表 1 发电机盘车记录
测量部位 上导 法兰 水导 1 1 -12 -10 2 1 -24 -30 3 1 -19 -24 4 0 -11 -13 5 -1 0 -2 6 -2 8 9 7 -1 -1 1 8 0 -7 -6
上导处的全摆度计算如下：
a 51 (1) 1 2 a 62 (2) 1 3 a 73 (1) 1 2 a 84 (0) 0 0
J mc X c2 Yc2
（21）
（22） （23）
mc COS 1 ( X c / J mc )
刮垫最大厚度的计算公式改为：
4 盘车处理软件开发与应用
J mc D L1 L2
（24）
根据以上思路与方法，采用C语言开发了盘车助手软件，选择相应的机组形式，进入相 应的操作界面。悬式机组的操作界面如图4所示。只需输入相应的摆度测量值，以及相关的 几何尺寸，点击计算便可给出净摆度计算结果，摆度曲线与测点示意；刮垫计算结果及轴线 状态示意。为方便使用，分别给出了法兰、水导与综合计算结果。伞式机组盘车处理界面如 图5，与悬式机组相比，只是多了一组输入集电环可上导的输入数据。
记：
aij
n
f k 0 f k 0 k 1 bi b j
n
（12）
aiy
有：
f k 0 ( yk f k 0 ) k 1 bi
（13）
a111 a12 2 ... a1m m a1 y a211 a22 2 ... a2 m m a2 y ......... am11 am 2 2 ... amm m amy
0 引言
对于立式水轮发电机组， 如果镜板摩擦面与机组轴线绝对垂直， 且组成轴线的各部件又 没有曲折和偏心，那么当机组回转时，机组的轴线将和理论回转中心相重合。但是，由于制 造与安装形成的误差， 实际的镜板摩擦面与机组轴线不会绝对垂直， 轴线本身也不会是一条 理想的直线。因而在机组回转时，机组中心线就要偏离理论中心线。那么不可避免地会使各 部件产生偏摆， 这样由导瓦确定的其理论运行空间变小， 在运转过程中轴系会产生径向平移， 导致导轴承负荷的加大(除正常水推力和不均衡磁拉力外的力)，当这种情况达到一定程度 时，导瓦的损伤就不可避免。因此，机组轴线的测量、分析和调整是机组检修中最为关键的 工序工序之一。 盘车时间的长短将直接关系到整个安装检修的直线工期控制， 盘车质量的好 [1] 坏对机组的运行状况也有着重要的影响 。 2] 机组轴线调整通常采用传统盘车人工等角测量[ 。 传统的等角盘车人工读数方法存在不 可克服的失真弊病，如：测数不准，无法一步调整到位，造成反复、无效劳动多；劳动强度 ［3-5］ 大，工作效率低。近年来，自动化盘车系统 大量应用于机组检修。但这些盘车系统仅解 决了自动转动问题，并不能提高盘车的速度与工艺，且已有自动化盘车系统造价高、维护复 ［6-7］ 杂，无法得到很好应用 。 无论是人工等角测量方法， 还是自动盘车方法， 轴线测量数据的分析与处理是其中重要 ［8］ 的一环。为了提高盘车的准确与适用性，文献 提出了一种基于改进遗传算法提取正弦信 ［ 9］ 号参数的求解，文献 建立了任意角盘车摆度公式。本文拟基于测量的相关数据，结合非 线性函数拟合、最小二乘算法，并采用高斯—牛顿法进行求解，求取各测量点的摆度向量， 并开发相应软件，自动计算机组轴线调整参数，给出合理的调整方案，加快机组轴线调整进 程，提高盘车质量。
J ba b 62 a 62 17.5 2
(3)
由上述计算值可知：法兰处的最大倾斜在测点号“6”处（方位角为225度），其值为： (4)
水导处的最大倾斜同样在测点号“6”处（方位角为225度），其值为：
J ca
c 62 a 62 21 2
(5)
2 盘车数据拟合
(0)
（14）
当数据点（xk,yk）和初始近似值 bi 给定后，由式（12）和式（13）可以算出各系数 的值，再由式（14）则可以计算出 i ，然后由式（8）可算出 bi 。当 i 的绝对值较大时， 重得上述计算过程。当 i 的绝对值小于某一事先确定的数值时，计算完成。 由表2的数据，由上述方法可得：
1 盘车数据处理
采用等角8点法，通过盘车已得到推力轴承处导轴承的摆度值、发电机主轴连接法兰处 的摆度值，以及水轮机导轴承处的摆度值如表1. 对于连续盘车方式，可以通过相位换算得 出如表1所示的等效的等角8点的数据。
_______________________ 作者简介：谢建国（1974-），男，四川双流，大学，高级工程师，主要从事水电站生产技术管理工作。E_mail: 70015680@qq.com
法兰处的全摆度计算如下：
(1)
b 51 0 (12) 12 b 62 8 (24) 32 b 73 (1) (19) 18 b84 (7) (11) 4
水导处的全摆度计算如下：
(2)
c 51 (2) (10) 8 c 62 9 (30) 39 c 73 1 (24) 25 c84 (6) (13) 7
以法兰处为例，假定测量的净摆度值符合正弦变化规律，则有
F J m sin( 0 ) J 0
(6)
式（6）为一非线性方程，可按最小二乘意义对一般的非线性函数进行曲线拟合。对于 给定的n 对数据点（xi，yi）,i=1,2,…n,要求确定函数y=f(x,B)中的非线性参数B，使得式 （7）为最小。
n
Q [ yk f ( xk , B)]2
k 1
（7）
对于由式（6）所表述的过程，变量x为单变量 x 。参数 B (b1 , b2 ,..., bm ) B为三个， 即0 、 J 0 和 J m 。 对于这种多参数的非线性问题，这里采用高斯—牛顿法，通过逐次“线性化”的间接方 法处理。其求解方法如下： 先给bi一个近似初值，记为 bi ,并记初值与真值之差为 i ，有：
水轮发电机组盘车数据处理及软件开发
谢建国1 方仲超1 程远楚2 王豪2 郑小贺2
（1.宝珠寺水力发电厂，四川广元 6280032；2.武汉大学，湖北武汉 430072）
摘 要：盘车是立式水轮发电机组安装检修过程中一项十分重要的工作，良好的盘车质量是机组健康 运行的重要保证。传统的人工等角测量方法，劳动强度大，工作效率低下。自动盘车方法不能从根本 上提高盘车速度及工艺。本文采用非线性函数拟合方法进行盘车数据处理，并提出了一种综合考虑法 兰处摆度和水导处摆度的刮垫厚度和最大摆度方位计算方法，并开发了相应的盘车数据处理工具软 件。实际应用表明，该方法应用方便，加快了盘车速度，提高了盘车质量。 关键词：盘车；数据拟合；软件开发 中图分类号：TK730.8 文献标识码：A
图 2 轴线内折
图 3 轴线外折
当轴线如图3实线所示时，若以法兰处的摆度值进行处理，水导处的摆度会改善，但仍 然较大 （图3中虚线）；若以水导处的摆度值进行处理，则法兰处的摆度同样改善不明显 （图 3中的点线）。 显然，应综合考虑水导处的摆度与法兰处的摆度进行处理。为此，本文提出如下算法： 已知法兰处摆度值换算到水导处的向量为：
理论上讲， 盘车时在法兰处和水导处所测得的净摆度值应符合正弦变化规律， 并可在正 弦曲线中找到最大摆度值及其对应的方位。但是，由于加工误差、测量误差及在盘车的受力 不均匀等诸多因素的影响，使得测量到的数据不完全满足正弦曲线规律。因此，需提取测量 [8] 数据中的有效成分，文献 采用遗传算法实现了连续采样数据的正弦拟合，突破了传统等角 8点盘车的局限。 本文结合等角8点盘车和连续盘车法， 利用非线性函数拟合和最小二乘算法， 给出一种盘车数据的拟合方法。 对于表1所述的数据，可以计算出法兰处和水导处的净摆度值如表2。 表 2 净摆度计算数据 测量部位 法兰 水导 1 -13 -11 2 -25 -31 3 -20 -25 4 -11 -13 5 1 -1 6 10 11 7 0 2 8 -7 -6
Ffc J mfc mfc X f jY f
水导处的摆度向量为：
.
.
（18）
Fs J ms ms X s jYs
设综合摆度向量为：
（19）
Fc J mc mc X c jYc
则有：
.
(20)
Xc (X f Xs) / 2 Yc (Y f Ys ) / 2
图 4 悬式机组盘车处理
图 5 伞式机组盘车处理
5 结论
本文采用非线性函数拟合方法进行盘车数据处理， 并提出了一种综合考虑法兰处摆度和 水导处摆度的刮垫厚度和最大摆度方位计算方法。 根据所提出方法， 开发了相应的盘车数据 处理工具软件，并在水电机组的安装与检修中得到具体应用，提高了轴线处理的精准性，加 快了盘车时度，取得了良好的效果。
（9）
当 bi 给定时，其值为直接算出。根据最小二乘法原理，由式（9）可推出：
Q [ yk f ( xk , b1 , b2 ,..., bm )]
k 1
n
（10）
n n n Q Q f f f f f 2[1 k 0 k 0 ... m k 0 k 0 k 0 ( yk f k 0 )] （11） bi i k 1 b k 1 bm bi k 1 bi 1 bi
பைடு நூலகம்
bi bi(0) i
(0)
（8）
这样，确定 bi 的问题转化为确定修正值 i 的问题。为确定 bi 在 bi 附近对f作泰勒级数 展开，并略去 i 的高阶项得：
f ( xk , b1 , b2 ,..., bm ) f k 0
(0)
f k 0 f f 1 k 0 2 .... k 0 m b1 b2 bm
按水导处的净摆度计算有：
J ca D L1 L2
（17）
式中：L1为法兰摆度测量处至水导摆度测量处的距离。 但轴线往往不是一条直线，按式（16）和式（17）计算出的结果是不一样的。根据轴线 曲线程度的不同，单独按式（16）或式（17）进行计算，有时得不到好的效果。 当轴线如图2实线所示时，若以法兰处的摆度值进行处理，水导处的摆度会加大（图2 中虚线）；若以水导处的摆度值进行处理，则法兰处的摆度没有明显改善（图2中点线）。
图 1 摆度曲线与测点关系
3 机组轴线调整
当机组摆度过大时， 可通过刮削推力头底面或绝缘垫或法兰面的方法进行调整。 当机组 轴线是一条直线时， 可按法兰处的净摆度和水导处的净摆度计算刮削推力头底面或绝缘垫的 最大厚度。 按法兰处的净摆度计算有：
J ba D L1
（16）
式中：L1为推力头底面至法兰摆度测量处的距离，D为推力头或绝缘垫外径。
Ff 14.9036 sin( 142.95) 8.125 Fs 17.6918 sin( 150.59) 9.25
(15)
最大摆度方位角分别为232.95和240.59。摆度值和最大摆度方位角的计算相比式（4） 和式（5）更为精准。 计算出摆度曲线和测点值如图1。 由图1可以判断测点数据是否符合正弦规律。 若测点数 据严重偏离正弦规律，说明盘车数据不能真实反映轴线状态，需重新进行盘车。
参考文献
[1] [3] [5] 华北水利水电学院等，水力机组安装与检修，电力工业出版社，1981. 罗清萍 . 电动盘车装置的原理与应用实例剖析［J］. 黑龙江水利科技，2013，41（4）：13-16 . 刘润根，杨虹，李辉，等. 水电站机组智能盘车系统设计与应用[J]. 中国水利水电科学研究院学 报. 2015(01): 74-80. [6］ 李则泉， 陈苗 . 锦屏二级水电站7号机组盘车与轴线调整新工艺 ［J ］ . 水电站机电技术， 2014 （4） . [7］ 王胜超 . 自动盘车装置在周宁水电站的应用［J］. 防爆电机，2013（2）. [8］ 孟安波，徐海波 . 基于遗传算法的水电机组盘车数据处理方法研究［J］. 水力发电，2011，37 （10）：61-62 . [9］ 方建新，谭啸 . 浅析状态监测系统的故障分析与盘车应用［J］ . 水电站机电技术，2009，32（3）： 96-98 . [2] 王胜超. 简析水轮发电机组盘车[J]. 电工文摘. 2012(04): 54-58. [4］ 王胜超 . 自动盘车装置在周宁水电站的应用［J］. 防爆电机，2013，48（2）：46-47 .
表 1 发电机盘车记录
测量部位 上导 法兰 水导 1 1 -12 -10 2 1 -24 -30 3 1 -19 -24 4 0 -11 -13 5 -1 0 -2 6 -2 8 9 7 -1 -1 1 8 0 -7 -6
上导处的全摆度计算如下：
a 51 (1) 1 2 a 62 (2) 1 3 a 73 (1) 1 2 a 84 (0) 0 0
J mc X c2 Yc2
（21）
（22） （23）
mc COS 1 ( X c / J mc )
刮垫最大厚度的计算公式改为：
4 盘车处理软件开发与应用
J mc D L1 L2
（24）
根据以上思路与方法，采用C语言开发了盘车助手软件，选择相应的机组形式，进入相 应的操作界面。悬式机组的操作界面如图4所示。只需输入相应的摆度测量值，以及相关的 几何尺寸，点击计算便可给出净摆度计算结果，摆度曲线与测点示意；刮垫计算结果及轴线 状态示意。为方便使用，分别给出了法兰、水导与综合计算结果。伞式机组盘车处理界面如 图5，与悬式机组相比，只是多了一组输入集电环可上导的输入数据。
记：
aij
n
f k 0 f k 0 k 1 bi b j
n
（12）
aiy
有：
f k 0 ( yk f k 0 ) k 1 bi
（13）
a111 a12 2 ... a1m m a1 y a211 a22 2 ... a2 m m a2 y ......... am11 am 2 2 ... amm m amy
0 引言
对于立式水轮发电机组， 如果镜板摩擦面与机组轴线绝对垂直， 且组成轴线的各部件又 没有曲折和偏心，那么当机组回转时，机组的轴线将和理论回转中心相重合。但是，由于制 造与安装形成的误差， 实际的镜板摩擦面与机组轴线不会绝对垂直， 轴线本身也不会是一条 理想的直线。因而在机组回转时，机组中心线就要偏离理论中心线。那么不可避免地会使各 部件产生偏摆， 这样由导瓦确定的其理论运行空间变小， 在运转过程中轴系会产生径向平移， 导致导轴承负荷的加大(除正常水推力和不均衡磁拉力外的力)，当这种情况达到一定程度 时，导瓦的损伤就不可避免。因此，机组轴线的测量、分析和调整是机组检修中最为关键的 工序工序之一。 盘车时间的长短将直接关系到整个安装检修的直线工期控制， 盘车质量的好 [1] 坏对机组的运行状况也有着重要的影响 。 2] 机组轴线调整通常采用传统盘车人工等角测量[ 。 传统的等角盘车人工读数方法存在不 可克服的失真弊病，如：测数不准，无法一步调整到位，造成反复、无效劳动多；劳动强度 ［3-5］ 大，工作效率低。近年来，自动化盘车系统 大量应用于机组检修。但这些盘车系统仅解 决了自动转动问题，并不能提高盘车的速度与工艺，且已有自动化盘车系统造价高、维护复 ［6-7］ 杂，无法得到很好应用 。 无论是人工等角测量方法， 还是自动盘车方法， 轴线测量数据的分析与处理是其中重要 ［8］ 的一环。为了提高盘车的准确与适用性，文献 提出了一种基于改进遗传算法提取正弦信 ［ 9］ 号参数的求解，文献 建立了任意角盘车摆度公式。本文拟基于测量的相关数据，结合非 线性函数拟合、最小二乘算法，并采用高斯—牛顿法进行求解，求取各测量点的摆度向量， 并开发相应软件，自动计算机组轴线调整参数，给出合理的调整方案，加快机组轴线调整进 程，提高盘车质量。
J ba b 62 a 62 17.5 2
(3)
由上述计算值可知：法兰处的最大倾斜在测点号“6”处（方位角为225度），其值为： (4)
水导处的最大倾斜同样在测点号“6”处（方位角为225度），其值为：
J ca
c 62 a 62 21 2
(5)
2 盘车数据拟合
(0)
（14）
当数据点（xk,yk）和初始近似值 bi 给定后，由式（12）和式（13）可以算出各系数 的值，再由式（14）则可以计算出 i ，然后由式（8）可算出 bi 。当 i 的绝对值较大时， 重得上述计算过程。当 i 的绝对值小于某一事先确定的数值时，计算完成。 由表2的数据，由上述方法可得：
1 盘车数据处理
采用等角8点法，通过盘车已得到推力轴承处导轴承的摆度值、发电机主轴连接法兰处 的摆度值，以及水轮机导轴承处的摆度值如表1. 对于连续盘车方式，可以通过相位换算得 出如表1所示的等效的等角8点的数据。
_______________________ 作者简介：谢建国（1974-），男，四川双流，大学，高级工程师，主要从事水电站生产技术管理工作。E_mail: 70015680@qq.com
法兰处的全摆度计算如下：
(1)
b 51 0 (12) 12 b 62 8 (24) 32 b 73 (1) (19) 18 b84 (7) (11) 4
水导处的全摆度计算如下：
(2)
c 51 (2) (10) 8 c 62 9 (30) 39 c 73 1 (24) 25 c84 (6) (13) 7
以法兰处为例，假定测量的净摆度值符合正弦变化规律，则有
F J m sin( 0 ) J 0
(6)
式（6）为一非线性方程，可按最小二乘意义对一般的非线性函数进行曲线拟合。对于 给定的n 对数据点（xi，yi）,i=1,2,…n,要求确定函数y=f(x,B)中的非线性参数B，使得式 （7）为最小。
n
Q [ yk f ( xk , B)]2
k 1
（7）
对于由式（6）所表述的过程，变量x为单变量 x 。参数 B (b1 , b2 ,..., bm ) B为三个， 即0 、 J 0 和 J m 。 对于这种多参数的非线性问题，这里采用高斯—牛顿法，通过逐次“线性化”的间接方 法处理。其求解方法如下： 先给bi一个近似初值，记为 bi ,并记初值与真值之差为 i ，有：
水轮发电机组盘车数据处理及软件开发
谢建国1 方仲超1 程远楚2 王豪2 郑小贺2
（1.宝珠寺水力发电厂，四川广元 6280032；2.武汉大学，湖北武汉 430072）
摘 要：盘车是立式水轮发电机组安装检修过程中一项十分重要的工作，良好的盘车质量是机组健康 运行的重要保证。传统的人工等角测量方法，劳动强度大，工作效率低下。自动盘车方法不能从根本 上提高盘车速度及工艺。本文采用非线性函数拟合方法进行盘车数据处理，并提出了一种综合考虑法 兰处摆度和水导处摆度的刮垫厚度和最大摆度方位计算方法，并开发了相应的盘车数据处理工具软 件。实际应用表明，该方法应用方便，加快了盘车速度，提高了盘车质量。 关键词：盘车；数据拟合；软件开发 中图分类号：TK730.8 文献标识码：A
图 2 轴线内折
图 3 轴线外折
当轴线如图3实线所示时，若以法兰处的摆度值进行处理，水导处的摆度会改善，但仍 然较大 （图3中虚线）；若以水导处的摆度值进行处理，则法兰处的摆度同样改善不明显 （图 3中的点线）。 显然，应综合考虑水导处的摆度与法兰处的摆度进行处理。为此，本文提出如下算法： 已知法兰处摆度值换算到水导处的向量为：
理论上讲， 盘车时在法兰处和水导处所测得的净摆度值应符合正弦变化规律， 并可在正 弦曲线中找到最大摆度值及其对应的方位。但是，由于加工误差、测量误差及在盘车的受力 不均匀等诸多因素的影响，使得测量到的数据不完全满足正弦曲线规律。因此，需提取测量 [8] 数据中的有效成分，文献 采用遗传算法实现了连续采样数据的正弦拟合，突破了传统等角 8点盘车的局限。 本文结合等角8点盘车和连续盘车法， 利用非线性函数拟合和最小二乘算法， 给出一种盘车数据的拟合方法。 对于表1所述的数据，可以计算出法兰处和水导处的净摆度值如表2。 表 2 净摆度计算数据 测量部位 法兰 水导 1 -13 -11 2 -25 -31 3 -20 -25 4 -11 -13 5 1 -1 6 10 11 7 0 2 8 -7 -6
Ffc J mfc mfc X f jY f
水导处的摆度向量为：
.
.
（18）
Fs J ms ms X s jYs
设综合摆度向量为：
（19）
Fc J mc mc X c jYc
则有：
.
(20)
Xc (X f Xs) / 2 Yc (Y f Ys ) / 2
图 4 悬式机组盘车处理
图 5 伞式机组盘车处理
5 结论
本文采用非线性函数拟合方法进行盘车数据处理， 并提出了一种综合考虑法兰处摆度和 水导处摆度的刮垫厚度和最大摆度方位计算方法。 根据所提出方法， 开发了相应的盘车数据 处理工具软件，并在水电机组的安装与检修中得到具体应用，提高了轴线处理的精准性，加 快了盘车时度，取得了良好的效果。
（9）
当 bi 给定时，其值为直接算出。根据最小二乘法原理，由式（9）可推出：
Q [ yk f ( xk , b1 , b2 ,..., bm )]
k 1
n
（10）
n n n Q Q f f f f f 2[1 k 0 k 0 ... m k 0 k 0 k 0 ( yk f k 0 )] （11） bi i k 1 b k 1 bm bi k 1 bi 1 bi