一种信号处理方法和装置及多普勒超声系统.pdf

《一种信号处理方法和装置及多普勒超声系统.pdf》由会员分享，可在线阅读，更多相关《一种信号处理方法和装置及多普勒超声系统.pdf（18页完整版）》请在专利查询网上搜索。

1、(10)授权公告号 CN 101991437 B (45)授权公告日 2013.09.04 CN 101991437 B *CN101991437B* (21)申请号 200910189727.8 (22)申请日 2009.08.24 A61B 8/06(2006.01) (73)专利权人 深圳迈瑞生物医疗电子股份有限 公司 地址 518057 广东省深圳市南山区高新技术 产业园区科技南十二路迈瑞大厦 (72)发明人 皮兴俊 (74)专利代理机构 广州华进联合专利商标代理 有限公司 44224 代理人 何平 (54) 发明名称 一种信号处理方法和装置及多普勒超声系统 (57) 摘要 本发明公开。

2、了一种信号处理方法和装置及多 普勒超声系统， 其将 A/D 采样器件的原始数据采 样率提高整数倍， 使其工作在过采样状态， 获得过 采样信号， 根据超声系统基线的自动调整频移量， 对所述过采样信号进行频移计算， 获得频移后的 信号 ； 进行滤波处理 ； 从滤波结果中抽取部分信 号， 恢复原始数据采样率对应的频谱成分， 获得降 采样结果 ； 对降采样结果进行反向频移调节。本 发明可消除 AD 器件引起的频谱黑线， 且去除过采 样引起的噪声干扰。 (51)Int.Cl. 审查员 桂叶晨 权利要求书 2 页 说明书 6 页 附图 9 页 (19)中华人民共和国国家知识产权局 (12)发明专利 权利要。

3、求书2页 说明书6页 附图9页 (10)授权公告号 CN 101991437 B CN 101991437 B *CN101991437B* 1/2 页 2 1. 一种应用于超声系统的信号处理方法， 其特征在于， 包括 ： 将 A/D 采样器件的原始数据采样率提高整数倍， 使其工作在过采样状态， 获得过采样 信号， 根据超声系统基线的自动调整频移量， 对所述过采样信号进行频移计算， 获得频移后 的信号 ； 进行低通滤波处理 ； 从滤波结果中抽取部分信号， 恢复原始数据采样率对应的频谱成分， 获得降采样结 果 ； 根据所述超声系统基线的自动调整频移量对降采样结果进行反向频移调节。 2. 根据权利。

4、要求 1 所述的方法， 其特征在于， 所述频移计算采用如下方式进行 ： x_shiftn xnexp(-j2a/k) 其中， x_shiftn表示频移后的数据序列， xn表示频移前的数据序列， a 为基线频移量， k 为原始数据采样率被提高的倍数， 且大于 1 的正整数。 3. 根据权利要求 1 所述的方法， 其特征在于， 所述低通滤波处理中， 截止角频率为 ： / k， 其中， k 为原始数据采样率被提高的倍数， 且大于 1 的正整数。 4. 根据权利要求 1 所述的方法， 其特征在于， 所述抽取过程是从滤波后的数据序列中 的每 k 个抽取一个， k 为原始数据采样率被提高的倍数， 且大于 。

5、1 的正整数。 5. 根据权利要求 1 所述的方法， 其特征在于， 所述根据超声系统基线的自动调整频移 量反向频移调节采用如下方式进行 ： x_minusn x_fsnexp(j2a) 其中， x_minusn表示反向频移后的数据序列， x_fsn表示降采样结果， a 为基线频移量。 6. 根据权利要求 1 所述的方法， 其特征在于， 所述方法还包括 ： 对反向频移后的数据进 行显示基线调节。 7. 一种应用于超声系统的信号处理装置， 其特征在于， 包括 ： 过采样模块， 用于以 A/D 采样器件的原始数据采样率的整数倍来进行数据采样， 获得 过采样信号 ； 频移模块， 用于根据超声系统基线的。

6、自动调整频移量， 对所述过采样信号进行频移计 算， 获得频移后的信号 ； 滤波模块， 用于对频移后的信号进行低通滤波处理 ； 降采样模块， 用于从所述滤波模块的输出结果中抽取部分信号， 恢复原始数据采样率 对应的频谱成分， 获得降采样结果 ； 反向频移调节模块， 用于根据所述超声系统基线的自动调整频移量对降采样结果进行 反向频移调节。 8. 根据权利要求 7 所述的装置， 其特征在于， 所述低通滤波处理中， 截止角频率为 ： / k， 其中， k 为原始数据采样率被提高的倍数， 且大于 1 的正整数。 9. 根据权利要求 7 所述的装置， 其特征在于， 所述降采样模块包括 ： 用于对所述滤波模。

7、块的输出结果进行循环计数的单元， 循环计数的最高值为 k， 其中， k 为原始数据采样率被提高的倍数， 且大于 1 的正整数 ； 用于缓存 k 个数据结果的单元 ； 权 利 要 求 书 CN 101991437 B 2 2/2 页 3 用于当单次计数到 k 时， 从缓存数据中任意抽取一个数据的单元。 10. 一种多普勒超声系统， 包括 ： 显示器、 AD 采样器件以及含有基线调节功能的数据处 理器， 其特征在于， 所述数据处理器包括 ： 过采样调节模块， 用于将所述 A/D 采样器件的原始数据采样率提高整数倍 ； 频移模块， 用于根据超声系统基线的自动调整频移量， 对从所述 A/D 采样器件获。

8、得的 过采样信号进行频移计算， 获得频移后的信号 ； 滤波模块， 用于对频移后的信号进行低通滤波处理 ； 降采样模块， 用于从所述滤波模块的输出结果中抽取部分信号， 恢复原始数据采样率 对应的频谱成分， 获得降采样结果 ； 反向频移调节模块， 用于根据所述超声系统基线的自动调整频移量对降采样结果进行 反向频移调节 ； 基线调节模块， 用于对所述反向频移调节模块的输出结果进行基线调节后输出至显示 器。 11. 根据权利要求 10 所述的系统， 其特征在于， 所述低通滤波处理中， 截止角频率为 ： /k， 其中， k 为原始数据采样率被提高的倍数， 且大于 1 的正整数。 12. 根据权利要求 1。

9、0 所述的系统， 其特征在于， 所述降采样模块包括 ： 用于对所述滤波模块的输出结果进行循环计数的单元， 循环计数的最高值为 k， 其中， k 为原始数据采样率被提高的倍数， 且大于 1 的正整数 ； 用于缓存 k 个数据结果的单元 ； 用于当单次计数到 k 时， 从缓存数据中任意抽取一个数据的单元。 权 利 要 求 书 CN 101991437 B 3 1/6 页 4 一种信号处理方法和装置及多普勒超声系统 技术领域 0001 本发明涉及超声系统信号处理技术领域， 具体涉及的是一种应用于多普勒超声系 统的信号处理方法和装置。 背景技术 0002 通常， 采用连续波多普勒技术测量血流速度的超声。

10、诊断仪中， 在应用连续波多普 勒技术检查心脏等器官的高速血流时， A/D 采样器件得到所需的采样率数据， 因 A/D 器件本 身的频谱特性， 当调节基线时， 频谱图上会出现很明显的水平黑线。 0003 如图 1 所示， 给出了传统的连续波超声多普勒频谱分析系统的基本信号处理流 程。根据焦点位置以及发射频率计算不同通道的发射相位以及不同通道的接收相位， 不同 的通道获取发射相位后发射方波， 同时送出所有通道的接收相位给解调器， 通道的接收信 号经过解调器后得到 IQ 分量， 不同通道 IQ 分量直接合并， 经过一系列的模拟处理后， AD 器 件将模拟信号变为数字信号。然后经过壁滤波处理。壁滤波是。

11、一个高通滤波器， 滤除因为 静止或者慢速的组织运动引起的强回波。经过壁滤波处理后， 将 IQ 两路分量送入谱估计模 块。谱估计模块一般采用快速傅立叶变换 (FFT) 估算功率谱。快速傅立叶变换的点数可以 是 128 点或者 256 点。由于估算出来的功率谱动态范围太大， 每次估算出来的功率谱需要 经过压缩处理， 压缩到灰度显示范围。 最后在屏幕上显示的多普勒频谱图是某段时间， 不同 频率偏移的强度。 0004 A/D 器件本身包含模拟滤波的功能， 根据输出的采样率， 决定模拟域的截止角频 率。 图2所示为该滤波器的幅频特性， 假设fs96K， 该模拟滤波期的通带截止频率为48K， 通带以及阻带。

12、均为等纹波， 过度带非常狭窄。图 3 给出了该滤波器的过渡带的频率衰减特 性。 0005 假设无信号输入， 频谱图输出全为噪声。从图 4 中可以看出在 fs/2 的频率处很 快衰减到 6dB， 正常情况下 fs/2 显示在屏幕的上下边， 频谱显示区域以外为黑色， 因此虽 然 fs/2 处的频谱衰减为黑色， 但是和频谱图显示外沿合在一起， 人眼并不敏感。大多数 连续波多普勒超声设备包含调节基线的功能时， 调节基线等同于将 FFT 的结果循环移位。 循环移位的结果使得零频基线不在频谱图的正中间位置， fs/2 频谱也不在频谱图的最上 沿， 最下沿显示， 而是在谱图的中间显示， 并且合在一起。如图 。

13、5 所示， 假设循环移位了 K 点 后的频谱图， 从图中可以看出， fs/2 的频谱值突然下降了 6dB， 人眼一下就感觉到该处的 灰度值比其它地方的灰度值低， 为一黑点， 分属不同竖直谱线的同一位置的黑点形成一条 水平黑线。 0006 由上可知， 当连续波多普勒采用带内等纹波， 过度带窄的模拟滤波器的 AD 器件 时， 在调节基线后， 若不采用特殊的处理技术， 必然存在水平频谱黑线。 0007 现有技术中， 通常按照实际采样率的 k 倍来提供 AD 采样时钟， 其中 k 为大于 1 的 正整数， 然后在后处理中简单求和平均后降 k 倍采样。不妨假设 k 2， 当采样率为 2fs时， 6dB 。

14、衰减处为 fs， 如图 6 所示， 如果不采取进一步的处理措施， 当调节基线时， fs的频率 说 明 书 CN 101991437 B 4 2/6 页 5 点移到谱线的中间， 同样会出现水平黑线。于是， 对输入的 2fs数据简单两点求和平均， 然 后间隔抽取一点， 得到采样率为fs的新的信号流。 求和平均的转移函数为 ： 0.5+0.5Z-， 其幅 频特性曲线如图 7 所示， 通过求和平均滤波后的频谱如图 8 所示。将平均之后的数据样点 每两个抽一个出来， 抽取之后采样率变为 fs， 从零频到 0.5fs， 幅频值逐渐减弱到 3dB， 频 谱显示如图 9 所示。当移动基线时， 0.5fs的频率。

15、点移到谱线的中间， 如图 10 所示， 和图 5 不同的是 ： 0.5fs处只有 3dB 的衰减， 而且是从零频的 0dB 逐渐过度到 3dB 衰减。因为 是逐渐过度， 人眼不敏感。消除了突然过度到 6dB 衰减引起的黑线。但是对于图 11 所示的 过采样信号频谱 ( 采样率为 2fs) 来说， 图中 fs/2， fs 段的频谱在经过基线频移后要求能 正确的显示出来。图 14 是简单求平均的情况， 从图中可以看出， 直接求和降采样时， fs/2， fs段的频谱失真严重 ； 从fs到2fs段属于负向频谱噪声， 是要求被过滤掉的。 但是， 简单求 和平均不能很好过滤掉该段的噪声干扰， 如图 12 。

16、所示， 在抽取时， 未完全抑制的负向噪声 会混叠到低频信号中。因此， 此种方法仍然不能很好的解决频谱黑线的问题， 从图 13 到 14 可以看出， 频谱突然衰减而引起的黑线仍然比较突出， 同时还引起了负向频谱噪声， 所以需 要进一步地改进现有技术， 以解决现有技术存在的问题。 发明内容 0008 本发明的目的在于提供一种信号处理方法和装置及多普勒超声系统， 其以过采样 后和调基线相对应的频移， 以及抽取滤波， 反向频移相结合的方式， 消除 AD 器件引起的频 谱黑线， 且去除过采样引起的噪声干扰。 0009 为了实现上述目的， 本发明采用如下技术方案 ： 0010 本发明提供了一种应用于超声系。

17、统的信号处理方法， 包括 ： 0011 将 A/D 采样器件的原始数据采样率提高整数倍， 使其工作在过采样状态， 获得过 采样信号， 0012 根据超声系统基线的自动调整频移量， 对所述过采样信号进行频移计算， 获得频 移后的信号 ； 0013 进行滤波处理 ； 0014 从滤波结果中抽取部分信号， 恢复原始数据采样率对应的频谱成分， 获得降采样 结果 ； 0015 对降采样结果进行反向频移调节。 0016 同时， 本发明还提供了一种应用于超声系统的信号处理装置， 包括 ： 过采样模块， 用于以 A/D 采样器件的原始数据采样率的整数倍来进行数据采样， 获得过采样信号 ； 频移 模块， 用于根。

18、据超声系统基线的自动调整频移量， 对所述过采样信号进行频移计算， 获得频 移后的信号 ； 滤波模块， 用于对频移后的信号进行滤波处理 ； 降采样模块， 用于从所述滤波 模块的输出结果中抽取部分信号， 恢复原始数据采样率对应的频谱成分， 获得降采样结果 ； 反向频移调节模块， 用于对降采样结果进行反向频移调节。 0017 基于上述装置， 本发明还提供了一种多普勒超声系统， 包括 ： 显示器、 AD 采样器件 以及含有基线调节功能的数据处理器， 所述数据处理器包括 ： 过采样调节模块， 用于将所述 A/D 采样器件的原始数据采样率提高整数倍 ； 频移模块， 用于根据超声系统基线的自动调 整频移量，。

19、 对从所述 A/D 采样器件获得的过采样信号进行频移计算， 获得频移后的信号 ； 滤 说 明 书 CN 101991437 B 5 3/6 页 6 波模块， 用于对频移后的信号进行滤波处理 ； 降采样模块， 用于从所述滤波模块的输出结果 中抽取部分信号， 恢复原始数据采样率对应的频谱成分， 获得降采样结果 ； 反向频移调节模 块， 用于对降采样结果进行反向频移调节 ； 基线调节模块， 用于对所述反向频移调节模块的 输出结果进行基线调节后输出至显示器。 0018 从上可见， 针对连续波多普勒采用 AD 器件的缺陷以及上述解决方法的缺点， 本发 明采用一种改进的过采样以及频移， 滤波， 反频移的处。

20、理方法， 来消除该水平黑线， 而且实 现方法简单， 同时也抑制了负向噪声。本发明特别适用于多普勒超声系统。 附图说明 0019 图 1 为连续波多普勒信号处理流图 ； 0020 图 2 为 AD 器件幅频特性图 ； 0021 图 3 为 AD 器件过度带详细幅频特性图 ； 0022 图 4 为基线位于频谱图正中的幅频特性图 ； 0023 图 5 为基线不在频谱正中的幅频特性图 ； 0024 图 6 为采样率提高两倍后的频谱图 ； 0025 图 7 为两点求和平均滤波器的幅频特性曲线 ； 0026 图 8 为求平均后的频谱特性 ； 0027 图 9 为降采样后的频谱特性 ； 0028 图 10 。

21、为移动基线后的频谱图 ； 0029 图 11 为假设信号频谱 ； 0030 图 12 为简单求和示意 ； 0031 图 13 为频移后的信号频谱 ； 0032 图 14 为低通滤波后的频谱 ； 0033 图 15 为抽取后的信号频谱 ； 0034 图 16 为反向频移以及调解显示基线后的频谱 ； 0035 图 17 为本发明流程图 ； 0036 图 18 为本发明应用于超声系统的信号处理装置 ； 具体实施方式 0037 以下将结合附图， 对本发明方法及其系统的具体实现方式进行详细说明。针对现 有技术方法的缺点， 本发明提供了一种信号处理方法， 如图 17 所示， 其包括 ： 0038 步骤 1。

22、00、 将 A/D 采样器件的原始数据采样率提高整数倍， 使其工作在过采样状 态， 获得过采样信号， 比如以系统要求的最高采样率的整数倍来完成对输入模拟信号的采 样 ； 0039 步骤 200、 根据超声系统基线的自动调整频移量， 对所述过采样信号进行频移计 算， 获得频移后的信号 ； 0040 步骤 300、 进行滤波处理 ； 0041 步骤 400、 从滤波结果中抽取部分信号， 恢复原始数据采样率对应的频谱成分， 获 得降采样结果 ； 说 明 书 CN 101991437 B 6 4/6 页 7 0042 步骤 500、 对降采样结果进行反向频移调节。 0043 以下以具体实施例详细说明上。

23、述方法的实现过程。 0044 首先， 假设原始数据采样率为 fs， 那么将原始数据采样率提高为 k 倍 fs， k 为大于 1 的正整数， 进行数据采样后， 获得信号序列为 xn。假设 k 2， 则信号频谱如图 11 所示， 此 时有用信号频谱位于 0 到 fs之间， fs到 2fs之间的频谱成份为负频率噪声。 0045 然后， 对信号序列 xn进行频移计算。对于降采样后的数据采样率 fs来说， 基线频 移量为 a， a 为一归一化因子， -0.5 a 0.5， 当 a 为负时， 表示负向频移， 当 a 为正时， 表 示正向频移， 最大频移量为 ： 0.5fs； 降采样之前相对于 kfs的归一。

24、化频移值为 ： a/k， 频移后 的数据序列为 ： x_shiftn 0046 则频移后的信号为 0047 x_shiftn xnexp(-j2a/k) n 0， 1， 2， ， ， ， ， ， 0048 频移后的信号频谱如图 13 所示。 0049 一般超声系统在显示前进行基线的自动调整， 同时计算基线的频移量， 将不同时 间的显示频移量大小进行统计， 生成旋转指数表， 然后就可以直接按照该旋转指数表的内 容， 对过采样后的信号执行时序的调制， 即频域的频移。 0050 其次， 进行滤波处理， 其目的是为了防止抽取混叠， 另外还可以滤出如图 13 所示 的在 fs/2 到 3fs/2 之间的。

25、频率噪声。因此， 此过程的最佳方式是采用低通滤波处理， 过滤 有用信号的带外干扰。假设将实际的数据采样率由 kfs降到 fs， 则低通滤波的截止角频率 为 ： /k， 滤波后的数据序列为 ： x_filtern， 滤波后的信号频谱如图 14 所示， 带外噪声被抑 制了。图 14 中向下箭头指示为因低通滤波狭窄的过渡带引起的黑线位置。 0051 再次， 从滤波后的数据序列 x_shiftn中抽取部分信号， 比如， 对数据序列 x_ filtern中的每 k 个抽取一个， 恢复原始数据采样率对应的频谱成分， 获得抽取后的数据序 列为 ： x_fsn， 信号频谱如图 15 所示， 图中向下箭头为抽取。

26、后的黑线位置。 0052 最后， 进行反向频移调节 ： 对数据序列 x_fsn执行反向频移， 归一化频移值为 ： a， 反向频移后的数据序列为 ： x_minusn， 则 0053 x_minusn x_fsnexp(j2a) n 0， 1， 2， ， ， ， ， ， 0054 反向频移以及调节显示基线后的信号频谱如图 16 所示。如图 17 所示， 为了适应 显示的要求， 通常还要对反向频移后的数据进行步骤 600 的基线调节或频移的操作。图 16 为最终显示在屏幕上的图谱， 从图中可以看出， 不论调节基线还是不调节基线， 频谱黑线始 终位于谱图显示区域的最上最下边。 0055 从上述过程可。

27、以看出， 本发明由过采样、 频移、 低通滤波、 抽取、 对应的反向频移 后， 使得显示基线始终位于谱图显示区域的最上最下边， 如图 16 所示， 从而保证了在进行 数据显示时消除频谱黑线的作用， 同时本发明还能避免频谱失真， 以及消除频谱噪声， 参见 图 12 和图 14 比较结果， 中间 fs/2 到 3fs/2 之间的频率噪声被消除了。 0056 基于上述方法， 通过软件编程可以对超声系统内的数据处理控制芯片或者中央处 理器进行内部功能的改进和重新设计， 从而获得如图 18 所示的应用于超声系统的信号处 理装置， 其包括 ： 0057 过采样模块， 用于以 A/D 采样器件的原始数据采样率。

28、的整数倍来进行数据采样， 获得过采样信号 ； 说 明 书 CN 101991437 B 7 5/6 页 8 0058 频移模块， 用于根据超声系统基线的自动调整频移量， 对所述过采样信号进行频 移计算， 获得频移后的信号 ； 0059 滤波模块， 用于对频移后的信号进行滤波处理 ； 0060 降采样模块， 用于从所述滤波模块的输出结果中抽取部分信号， 恢复原始数据采 样率对应的频谱成分， 获得降采样结果 ； 0061 反向频移调节模块， 用于对降采样结果进行反向频移调节。 0062 其中， 所述滤波模块采用低通滤波， 截止角频率为 ： /k， 其中， k 为原始数据采样 率被提高的倍数， 且大。

29、于 1 的正整数。 0063 其中， 所述降采样模块包括 ： 用于对所述滤波模块的输出结果进行循环计数的单 元， 循环计数的最高值为k， 其中， k为原始数据采样率被提高的倍数， 且大于1的正整数 ； 用 于缓存k个数据结果的单元 ； 用于当单次计数到k时， 从缓存数据中任意抽取一个数据的单 元。 0064 将上述装置的设计构思应用于传统的多普勒超声系统上时， 也可以通过软件编程 对多普勒超声系统内部含有基线调节功能的数据处理器进行性能的改进， 同时还能带来其 内部构造的改进， 从而获得如下所示的多普勒超声系统结构， 其包括 ： 显示器、 AD 采样器件 以及含有基线调节功能的数据处理器， 同。

30、时其数据处理器包括 ： 0065 过采样调节模块， 用于将所述 A/D 采样器件的原始数据采样率提高整数倍 ； 0066 频移模块， 用于根据超声系统基线的自动调整频移量， 对从所述 A/D 采样器件获 得的过采样信号进行频移计算， 获得频移后的信号 ； 0067 滤波模块， 用于对频移后的信号进行滤波处理 ； 0068 降采样模块， 用于从所述滤波模块的输出结果中抽取部分信号， 恢复原始数据采 样率对应的频谱成分， 获得降采样结果 ； 0069 反向频移调节模块， 用于对降采样结果进行反向频移调节 ； 0070 基线调节模块， 用于对所述反向频移调节模块的输出结果进行基线调节后输出至 显示器。

31、。 0071 上述各个功能模块所采用的结构和实现的方法可以参见上述相关说明， 其与上述 方法的对应步骤和上述应用于超声系统的信号处理装置的对应功能模块结构和实现方法 相同。 0072 综上所述， 本发明过采样模块按照系统要求的最高采样率的整数倍来完成对输入 模拟信号的采样， 频移模块按照显示频移的大小产生旋转指数表， 对采样后的信号执行时 序的调制， 即频域的频移 ； 为了防止抽取混叠， 滤波模块过滤有用信号的带外干扰 ； 完成抗 混叠凌波器后， 降低信号的采样率， 完成整数倍的抽取， 此时信号速率降为系统设计的采样 率 ； 根据频移值的大小， 反向频移处理， 恢复出信号原始的频谱成分 ； 在。

32、显示时， 按照频移 的大小， 将频谱分析后的数据移到合适的位置显示出来。 0073 本发明以过采样的方式来得到输入数据， 根据降采样后基线频移量的大小， 计算 降采样前的频移量， 然后对输入数据进行频移计算， 抽取之前滤波， 过滤掉带外干扰， 抽取 之后反向频移， 以恢复原先信号的成分， 该方法不仅抑制带外频谱成份混叠， 同时保证不论 是否调节基线， 频谱黑线始终位于显示屏的上方以及下方， 达到消除因调节基线产生的频 谱黑线的目的。 说 明 书 CN 101991437 B 8 6/6 页 9 0074 应当理解的是， 对本领域普通技术人员来说， 可以根据上述说明加以改进或变换， 而所有这些改。

33、进和变换都本应属于本发明所附权利要求的保护范围。 说 明 书 CN 101991437 B 9 1/9 页 10 图 1 图 2 图 3 说 明 书 附 图 CN 101991437 B 10 2/9 页 11 图 4 图 5 说 明 书 附 图 CN 101991437 B 11 3/9 页 12 图 6 图 7 说 明 书 附 图 CN 101991437 B 12 4/9 页 13 图 8 说 明 书 附 图 CN 101991437 B 13 5/9 页 14 图 9 图 10 说 明 书 附 图 CN 101991437 B 14 6/9 页 15 图 11 图 12 图 13 说 明 书 附 图 CN 101991437 B 15 7/9 页 16 图 14 图 15 图 16 说 明 书 附 图 CN 101991437 B 16 8/9 页 17 图 17 说 明 书 附 图 CN 101991437 B 17 9/9 页 18 图 18 说 明 书 附 图 CN 101991437 B 18 。