Research on Gesture Recognition Based on Pressure-based Mechanomyogram and Improved Neural Fuzzy Inference System
-
摘要: 手势识别是人机交互领域的重要研究内容, 为截肢患者控制智能假肢手提供基础. 当前主流方法之一是利用表面肌电图(Electromyogram, EMG)识别手部运动意图, 但肌电信号存在信号弱和易受噪声、汗液、疲劳影响等缺点. 同时肌电图在识别准确率方面, 尤其是截肢患者手势识别方面仍然具有较大的提升空间. 针对这些问题, 设计了基于气压肌动图(Pressure-based mechanomyogram, pMMG)的穿戴式信号采集装置, 为手势识别提供了优质的信号源. 结合深度神经网络中全连接层结构、典型抽样和标准正则化技术, 提出了一种改进多类神经模糊推理系统(Improved multicalss neural fuzzy inference system, IMNFIS), 与传统自适应神经模糊推理系统(Adaptive neural fuzzy inference system, ANFIS)相比, 泛化能力得到显著提升. 招募了7名健康受试者和1名截肢受试者, 并用8种算法开展离线实验. 所提方法在残疾人手势识别实验中取得了97.25%的最高平均准确率, 在健康人手势识别实验中取得了98.18%的最高平均准确率. 与近年公开报道的多种手势识别研究相比, 所提方法的综合性能更优.Abstract: Gesture recognition, which provides a foundation for amputees to control smart prosthetic hands, is an important research content in the field of human-robot interaction. One of the current mainstream methods is using surface electromyogram (EMG) to identify the intention of hand motion, while EMG signals are weak and show some shortages of being interfered easily by noise, sweat, fatigue, etc. At the same time, there is still large room for the improvement in the recognition accuracy, especially in the gesture recognition of amputees by using EMG. To solve these problems, a wearable signal acquisition device based on pressure-based mechanomyogram (pMMG) is designed in this paper, which provides a signal source of high quality for gesture recognition. We proposed an improved multicalss neural fuzzy inference system (IMNFIS) by combining the full connection layer structure, typical sampling and uniformed regularization techniques in deep neural network, which significantly improved the generalization ability compared with the traditional adaptive neural fuzzy inference system (ANFIS). We recruited seven healthy subjects and an amputee subject, and then conducted an offline experiment in which eight algorithms are used. The proposed method got the highest average accuracy of 97.25% in the experiment of the disabled, and 98.18% in the experiment of the healthy. Compared with many reported gesture recognition researches in recent years, the method proposed in this paper achieves the better comprehensive performance.
-
电网频率是电力系统运行中最重要的关键参数之一, 系统频率的异常波动不仅影响电网的电能质量, 制约电力设备的正常使用, 严重的频率失衡甚至将导致并网设备脱网及系统崩溃, 对电力系统的稳定运行带来严重危害[1-2]. 因此, 系统频率的稳定控制是电力系统稳定不可或缺的重要组成. 电力系统频率的异常波动是由并网机组发电功率和用电负荷功率不平衡引起的, 因而, 电力系统能够根据系统频率的变化主动快速进行发用电控制以维持频率稳定就显得尤为重要[3]. 电力系统的一次调频就是当系统频率发生异常波动时, 自动迅速地作出响应, 改变并网机组的实时出力, 将系统频率稳定在一个合理安全的水平线上. 这对提高电能质量及电网频率控制水平, 迅速平息电网频率波动起到了重要作用[4-5].
电力系统是一个非线性时变强耦合的复杂系统, 其一次调频特性与每台并网运行的发电机调速器特性息息相关, 是每台并网运行机组调速器特性在当前电网工况下的综合特性[6]. 对于单台发电机组而言, 其一次调频特性与调速器自身的特性显著相关[7]. 而当机组并入电力系统时, 电力系统的一次调频特性不仅与机组自身的调速器特性有关, 也和所有并列运行的其他发电机组特性相关, 甚至还与机组所处的实时工况和机组间的相互作用有关[8]. 另一方面, 电网中各机组的发电类型、容量、设备结构(例如, 凝汽式机组的负荷响应速度明显高于背压式机组), 以及机组运行工况(例如, 滑压或定压、不同季节乃至一天内不同时间段的负荷率高低)等各类因素共同作用影响着电网的一次调频性能[9]. 所以, 在不同的电网运行工况下, 非常有必要对电网的一次调频性能进行在线估计, 使电网面对不断变化的负荷调节需求均能够具备足够的一次调频性能.
目前, 国内外对一次调频能力的研究偏向单台机组侧较多, 对电网一次调频性能的研究和计算比较分散, 整体性考虑电网一次调性能各个因素的综合性研究较少, 也未见实际的工程化应用. 此外, 由于很难、甚至无法对电网一次调频特性进行精确建模, 以至于无法有效估计电网真实的一次调频性能.
本文基于平行控制理论框架下由人工社会 (Artificial societies, A)、计算实验(Computational experiments, C)和平行执行 (Parallel execution, P) 构成的ACP理论体系, 构建了多源数据的电网一次调频性能平行计算平台. 本文提出的计算方法能够有效兼顾机组类型的静态特性和运行工况的动态特性, 并以平行执行方式完成人工估算系统与实际电力系统的滚动优化, 实现了电网一次调频性能的在线全面估计, 从而有效弥补其无法进行精确建模的不足, 有助于提高电网主动应对负荷突变的能力, 也是对电网一次调频性能管控的一种全新的可行性方式.
1. 平台整体框架设计
所谓平行系统, 是指由某一个自然的现实系统和对应的一个或多个虚拟或理想的人工系统所组成的共同系统[10]. 平行系统主要包括实际系统和人工系统, 其基本框架如图1所示[10].
设计平行系统的最主要目的是通过人工系统与实际系统的相互连接, 对二者之间的行为进行比对和分析, 完成对各自未来状况的“借鉴”和“预估”, 相应地调节各自的管理与控制方式, 达到实施有效解决方案以及学习和培训的目的[10-11].
人工系统试图尽可能地模拟实际系统, 对其行为进行预估, 从而为寻找对实际系统有效的解决方案或对当前方案进行改进提供依据[12]. 进一步, 通过观察实际系统与人工系统评估的状态之间的不同, 生成误差反馈信号, 对人工系统的评估方式或参数进行修正, 并重复进行新一轮的优化和评估[12].
目前, 平行系统理论已经在自动化码头[13]、智能交通[14]等领域中开展了应用研究, 但尚未用于现代电力系统. 现代电力系统是一个巨大的参数海量、广域分布、模型复杂的非线性复杂系统. 现阶段, 绝大多数电力仿真系统是基于物理模型进行的, 并没有深入地考虑人的行为、自然因素以及社会环境的影响, 这些因素限制了对复杂现代电力系统运行影响的评估, 尤其是无法定性、定量地对电力系统应对突发事件的能力问题进行分析评估[15].
为了克服采用传统仿真方法导致的上述局限, 本文将平行系统理论应用到电力系统的实际工程中, 改变传统的被动、离线的电力系统仿真方法, 建立人工电力平行计算平台, 进而实现管理与控制、实验与评估、学习与优化及故障诊断等功能. 以现代电网的海量数据为驱动, 电力平行计算平台为载体, 该平台可对正在运行的电力系统进行滚动式在线改进与不间断优化, 实现电网一次调频性能估算达到最优解, 为电力系统和复杂电网的能力优化与管理开辟一条新的途径[16].
基于上述平行系统理念, 本文设计了电网调频能力计算平台, 其整体框架如图2所示, 主要包括实际系统和人工平行系统两部分. 实际电力系统的频率控制本质上是一个以电力调控中心为中枢, 实现电源与负荷的电力供需平衡的闭环控制系统. 电力调控中心不仅需要满足区域内机组发电及负荷消纳的控制要求, 也肩负着电力系统实际运行产生的海量数据信息采集的重担. 人工平行系统中, 基于D5000电网调度系统的一次调频在线测试系统不仅能够在正常运行工况下测试控制区内的机组调频性能, 保证机组一次调频性能的测试结果更具真实性, 同时将历次测试数据存储、分析及评价, 实现机组一次调频的静态参数、动态性能的数据积累和量化管理. 将电力调控中心的海量数据和一次调频测试系统的特性数据应用于电网一次调频能力在线估计模型中, 通过实际电力系统和人工平行系统的彼此作用、相互修正, 既可以实现对实际电力系统的管理与控制,也能够完成对电网能力、行为与决策的评估与实验.
图 2 基于平行系统的计算平台整体框架Fig. 2 Overall framework of computing platform based on parallel systems2. 发电机组一次调频数据采集体系
当前, 电力系统已经发展成一个具有多源信息融合的综合性电力控制系统. 根据现有的SCADA(Supervisory control and data acquisition)、WAMS(Wide area measurement system)、OMS (Operations management system)等系统, 建立机组一次调频数据采集体系, 如图3所示. SCADA的一次调频在线测试系统能够在一次调频在线测试期间自动记录测试机组的有功PMU (Phasor measurement unit)数据, 整编录入历史库存储[17]. WAMS的机组性能考核系统能够不间断采集机组实际运行的一次调频动作、响应等广域量测数据, 实现机组转速不等率、积分电量、响应时间等一次调频性能指标的在线计算[18]. 基于OMS的网源协调信息管理系统, 通过设置统一的一次调频试验台账, 维护和存储机组的历次原始试验数据和参数, 统一存储于数据库中.
图 3 发电机组一次调频数据采集体系Fig. 3 Primary frequency control data acquisition system of generator sets2.1 机组一次调频在线测试系统
如图4所示, 在D5000系统中部署试验触发、结果展示、参数维护及统计分析功能. 试验触发功能可实现选定机组进行试验, 试验开始时将测试机组投入一次调频扰动测试模式, 调控中心向电厂侧下发调频测试开始信号和频率偏差信号, 机组根据接收的信号进行一次调频动作响应, 具体控制信息流程见图5所示.
图 4 机组一次调频在线测试与评价功能结构图Fig. 4 Functional structure diagram of unit primary frequency control online test and evaluation在一次调频在线测试期间, 系统自动记录测试机组的有功PMU数据, 整编录入D5000历史库, 永久保存, 供后期调频指标分析与数据展示使用[17]. 在结果展示功能中, 将频率扰动信号及电厂侧PMU装置上传的有功出力信号进行展示, 并根据有功出力曲线自动计算一次调频动作性能指标. 在浙江某厂#2机组进行一次调频在线测试系统联调试验时, 调控中心主站侧的操作及结果展示如图6所示, 发电厂机组一次调频动作响应曲线如图7所示.
图 6 机组一次调频在线测试主站侧界面图Fig. 6 Interface diagram of main station side of unit primary frequency control online test
图 7 浙江某厂#2机组一次调频在线测试动作曲线(75%负荷点 + 11 r/min转速偏差)Fig. 7 Action curves of primary frequency control online test of unit 2 of a power plant in Zhejiang (+11 r/min speed deviation at 75% load point)通过PMU实时采集高频度的功率、相位、相角等机组数据, 并利用PMU的GPS对时优势, 计算机组在一次调频在线测试中不同时刻、不同工况下的性能指标, 反向推算出机组实际的调频特性参数, 以此形成海量的机组一次调频特征数据, 实现机组一次调频静态参数、动态性能等特征数据的分析及管理[19].
2.2 网源协调信息管理系统
基于OMS的网源协调信息管理系统(如图8所示)能够实现设备台账、电气参数、涉网试验报告等技术资料的及时更新及有效管理, 打通厂网之间的涉网信息共享交互渠道, 更加科学、系统、高效地开展网源协调管理工作. 以各机组台账信息为基础, 实现对网内各并网机组数据信息的实时统计汇总; 以并网机组各类信息关键字为基础, 实现对特定类型、特定机组相关信息的快速查询和准确统计, 该系统所部署的功能模块如图9所示. 同时, 建立网源协调数据中心, 这是网源协调信息管理平台的核心内容, 可实现参数模型库、文档资料库、督办信息库及信息变更库的数据获取与整合、数据质量校核、合格率督办和数据查询等服务.
图 8 基于OMS的网源协调信息管理系统主界面Fig. 8 Main interface of grid power coordination information management system based on OMS
图 9 网源协调信息管理系统的功能模块Fig. 9 Functional modules of grid power coordination information management system在网源协调信息管理系统中, 通过设置统一的一次调频试验台账模板, 采集机组的历次试验数据, 以维护和存储机组的原始试验数据. 图10为浙江某百万千瓦机组A修后一次调频试验的数据台账.
图 10 浙江某厂#1机组A修后一次调频试验台账Fig. 10 Ledger of primary frequency control test after class a maintenance of unit 1 of a power plant in Zhejiang根据设置统一的数据台账, 可以追溯机组每一次试验工况, 分析不同工况下机组一次调频试验时的真实性能.
2.3 发电机组一次调频性能考核系统
发电机组一次调频性能考核系统是基于WAMS部署于D5000上的一次调频数据量化及评价系统[18]. 该系统通过PMU将现场采集的机组出力数据送至电网调控中心EMS (Energy management system)系统中, 并由EMS采集电网中的频率值及时间, 打包形成一个数据包存储于该系统, 以便进行量化处理.
机组一次调频贡献电量是机组一次调频多项指标的综合, 它反映机组对电网频率的实际贡献大小. 这一参数的引入使量化机组一次调频功能的动作力度和快速性成为可能.
在系统设计中, 机组对象被抽象封装成一个类, 每个机组都是该类的实例化. 类中数据为机组的实时信息, 包括机组基本参数(代号、装机容量、积分间隔、频率超前时间等)和积分电量以及其他计算辅助变量; 类中对数据的操作包含积分电量的计算和对数据库的读写操作等. 其中最为核心的函数为贡献电量积分函数, 用来计算和累积贡献电量, 以反映机组在实际运行时的一次调频动作性能.
一次调频性能考核系统通过采集各机组一次调频投退信号、有功功率和频率信息等广域量测信息, 实现机组转速不等率、积分电量、响应时间等一次调频性能指标的在线计算, 使系统能够考核机组是否贡献了满足要求的一次调频服务.
该系统依托于调度D5000系统, 能够不间断实时获取所有并网机组的一次调频动作性能. 据系统统计, 每个月一台机组的一次调频动作次数可达到近5000次(如表1所示), 大量原始数据的累积奠定了电网一次调频性能计算的数据基础.
表 1 机组一次调频月动作统计Table 1 Monthly action statistics of unit primary frequency control月份 强蛟厂 #3 机组 镇燃厂 #11 机组 动作总次数 正确动作数 动作总次数 正确动作数 8 2 491 2 389 408 381 9 4 085 3 948 1 464 1 369 10 4 965 4 760 2 956 2 822 3. 电网一次调频能力数据建模
目前, 国内外各类指标都是针对机组调频能力, 比如控制性能标准(Control performance standard, CPS)、投运率、贡献电量、调频效果等[20]. 现有的评价指标本质为事后评价, 不能预测电网总的调频裕度. 本文构建的数据平行驱动体系(如图11所示), 依靠电网运行的大数据, 通过不断演算、迭代优化, 形成调频能力的概率区间来衡量电网调频性能, 解决了电网一次调频能力的在线估计问题, 从而使电网的运行状态具有更高的可控性及透明度, 提高大电网应对负荷突变的能力.
图 11 电网一次调频能力估算技术路线图Fig. 11 Technical roadmap of grid primary frequency control capacity estimation3.1 同步电网一次调频能力建模
定义 1. PFCA (Primary frequency control ability)表示仅在一次调频作用情况下(即没有二次调频作用时), 某时间段内负荷变化与电网频率变化的比值[8], 可表示为
$$ {R_{{\rm{PFCA}}}} = \frac{{\left| {\text{电网负荷的变化}} \right|}}{{\left| {\text{电网频率的变化}} \right|}} $$ (1) 同步发电机并列运行的电网调频能力模型如图12所示. 在图12中,
${R_i}$ 为二次调频的给定值;${\delta _i}$ 为电网中第i台机组的不等率;${\alpha _i}$ 为电网中第i台机组装机容量/电力系统装机容量;${G_i}(s)$ 为第i台汽轮机的传递函数;${P_{NL}}(s)$ 为负载的标幺值;${\varphi _s}(s)$ 为电网转速变化的标幺值.${T_{\alpha \Sigma }} = \sum\nolimits_{i = 1}^M {{\alpha _i}} {T_{\alpha i}}$ 为电网的惯性时间常数;${T_{\alpha i}}$ 为第i台汽轮机的转子时间常数;$\beta _{\Sigma }$ 为电网负载频率特性系数.
图 12 M台同步发电机并列运行的电网调频能力模型图Fig. 12 Model diagram of power grid frequency regulation capacity of M synchronous generators operating in parallel$${\varphi _s}(s) = - \frac{1}{{{T_{\alpha \Sigma }}s + {\beta _{\Sigma} } + G(s)}}{P_{NL}}(s)$$ (2) 其中,
$$G(s) = \sum\limits_{i = 1}^M {\frac{{{\alpha _i}}}{{{\delta _i}}}} {G_i}(s)$$ (3) 各变量前加
$\Delta $ 作为静态标记, 静态下电网周波变化可表示为$$\Delta \varphi = - \mathop {\lim }\limits_{s \to 0} s\dfrac{1}{{{T_{\alpha \Sigma }}s + {\beta _{\Sigma} } + G(s)}}\Delta {P_{NL}}\dfrac{1}{s}$$ (4) 由于汽轮机传递函数可以表示为
$${G_i}(s) = \frac{{\prod\limits_{i = 1}^m {({T_i}s + 1)} }}{{\prod\limits_{i = 1}^n {({T_i}s + 1)} }}$$ (5) 故
$$\Delta \varphi = -\frac{ \Delta {P_{NL}}}{\sum\limits_{i = 1}^M {\frac{{{\alpha _i}}}{{{\delta _i}}}} + {\beta _{\Sigma} }}$$ (6) 由此得到, 在静态情况下电网的一次调频能力可表示为
$${R_{{\rm{PFCA}}}}(static) = {R_G} + {R_L}\qquad\qquad\qquad\;\;$$ (7) $${R_{{\rm{PFCA}}}}(static) = \frac{{\left| {\Delta {P_{NL}}} \right|}}{{\left| {\Delta \varphi } \right|}} = \sum\limits_{i = 1}^M {\frac{{{\alpha _i}}}{{{\delta _i}}}} + {\beta _{\sum} }$$ (8) ${R_G} = \sum\nolimits_{i = 1}^M ({{{{\alpha _i}}}/{{{\delta _i}}}})$ 定义为电网的等值调差率, 等效为所有挂网同步发电机组一次调频出力总和的最大值.${R_L} = {\beta _{\sum} }$ 为电网负载频率特性系数, 等效为电网负荷对频率变化的贡献能力.对于某台额定功率
${P_{ei}}$ , 不等率${\delta _i}$ 的机组而言, 其在$\Delta f$ 频率扰动下的一次调频能力的实际贡献为$${k_i} = \frac{{\Delta P_{{\rm{real}}}^ * }}{{\Delta {f^ * }}} = \frac{{{\eta _i}}}{{{\delta _i}}}\qquad\qquad\qquad\qquad\qquad\;\;$$ (9) $${\eta _i} = \frac{{\Delta {P_{{\rm{real}}}}}}{{\Delta {P_{{\rm{ expect}}}}}} =\frac{ \Delta {P_{{\rm{real}}}}}{\frac{{\Delta f \times {P_{ei}}}}{{{f_e} \times {\delta _i}}}} = \frac{{\Delta P_{{\rm{real}}}^ * }}{{\Delta {f^ * }}} \times {\delta _i}$$ (10) 显然, 在
$\Delta f$ 频率扰动下的同步电网中所有挂网运行机组的一次调频实际出力总和为$$ {R_{{\rm{PFCA}}}}(real) = \sum\limits_{i = 1}^M {{\alpha _i}\frac{{{\eta _i}}}{{{\delta _i}}}} = \sum\limits_{i = 1}^M {{\alpha _i}{k_i}} \le {R_G} = \sum\limits_{i = 1}^M {\dfrac{{{\alpha _i}}}{{{\delta _i}}}} $$ (11) 3.2 机组一次调频能力数据模型
从式(11)中可以发现电网一次调频能力主要取决于各挂网运行机组的转速不等率.
机组在不同功频下的转速不等率不尽相同, 将每台机组的转速不等率设定为常数势必难以准确描述机组的一次调频能力. 调度D5000系统具有各类丰富的一次调频数据(如图13所示), 其中, 机组性能考核系统实时采集的机组一次调频动作数据量最大, 每月达几千次. 而网源协调信息管理系统中存储的一次调频参数及试验数据, 则表征机组一次调频性能的理论出力和实际能力上下限值. 一次调频在线测试系统可以探测机组在不同工况下的大频差真实能力, 完善整个机组的一次调频数据集. 本文将采用极大似然估计及数值拟合等算法对调度D5000系统中的多源一次调频数据进行有效整合, 以获取机组真实的一次调频性能功频特性图谱.
图 13 发电机组一次调频动作历史数据库Fig. 13 Historical database of primary frequency control actions of generator sets可以从D5000系统中获取不同时刻机组一次调频能力的数据, 从而构造出机组转速不等率的样本空间. 我们用随机变量
${\delta _i}$ 表示单台机组在第i个时间段内的转速不等率. 利用D5000系统中不同时间的转速不等率数据, 可以构造用于描述转速不等率的随机变量${\delta _1},\cdots,{\delta _n}$ . 由于机组的短期特性通常不会发生变化, 可以认为${\delta _1},\cdots,{\delta _n}$ 独立同分布. 按照中心极限定理, 随机变量${\delta _1},\cdots,{\delta _n}$ 独立同分布, 当n很大时,$({{\sum\nolimits_{i = 1}^n {{\delta _i}} - n\mu }})/({{\sqrt n \sigma }})$ 近似服从标准正态分布${\rm{N}}({\rm{0,1}})$ . 设$X = \sum\nolimits_{i = 1}^n {{\delta _i}}$ , 则X服从正态分布${\rm{N}}(n\mu ,n{\sigma ^2})$ .本文通过极大似然估计来估计
$\mu $ 和${\sigma ^2}$ , 似然函数可以表示为$$ \begin{split} L(n\mu ,n{\sigma ^2}) =\;& \prod\limits_{i = 1}^N {\frac{1}{{\sqrt {2{\text{π}} } \sigma }}{{\rm{e}}^{ - \frac{{{{({X_i} - \mu )}^2}}}{{2{\sigma ^2}}}}}}= \\ & {(2{\text{π}} {\sigma ^2})^{ - \frac{N}{2}}}{{\rm{e}}^{ - \frac{1}{{2{\sigma ^2}}}\sum\limits_{i = 1}^N {{{({X_i} - \mu )}^2}} }} \end{split} $$ (12) 那么, 似然函数的对数为
$$ \begin{split} &\lg L(n\mu ,n{\sigma ^2}) = - \frac{N}{2}\lg (2{\text{π}} {\sigma ^2}) - \frac{1}{{2{\sigma ^2}}}\sum\limits_{i = 1}^N {{{({X_i} - \mu )}^2}} = \\ &\qquad - \frac{N}{2}\lg (2{\text{π}} ) - \frac{N}{2}\lg ({\sigma ^2}) - \frac{1}{{2{\sigma ^2}}}\sum\limits_{i = 1}^N {{{({X_i} - \mu )}^2}} \end{split} $$ (13) 则似然方程组为
$$\left\{ {\begin{aligned} & {\frac{{\partial \lg L(n\mu ,n{\sigma ^2})}}{{\partial \mu }} = \frac{1}{{{\sigma ^2}}}\sum\limits_{i = 1}^N {({X_i} - \mu ) = 0} } \qquad\;\;\,(14{\rm{a}})\\ & {\frac{{\partial \lg L(n\mu ,n{\sigma ^2})}}{{\partial {\sigma ^2}}} = - \frac{N}{{2{\sigma ^2}}} + \frac{1}{{2{\sigma ^4}}}\sum\limits_{i = 1}^N {{{({X_i} - \mu )}^2} = 0} } \end{aligned}} \right.\tag{14b}$$ 求解(14a), 得
$$n\hat \mu = \overline X = \frac{1}{N}\sum\limits_{i = 1}^N {{X_i}} $$ (15) 代入(14b), 得
$$n{\hat \sigma ^2} = \frac{1}{N}\sum\limits_{i = 1}^N {({X_i} - \overline X } {)^2}$$ (16) 那么,
$\mu $ 和${\sigma ^2}$ 的最大似然估计为$$ \hat \mu = \frac{{\rm{1}}}{n}\overline X , \;{\hat \sigma ^2} = \frac{{\rm{1}}}{n}\frac{1}{N}\sum\limits_{i = 1}^N {({X_i} - \overline X } {)^2} $$ (17) 本文用大写字母表示所有涉及的样本, 因为最大似然估计
$\hat \mu $ 和${\hat \sigma ^2}$ 均为统计量, 脱离具体的试验或观测, 它们都是随机的. 依据3套系统得到3个维度的真实观测数据, 采用上述的估计算法, 即可得到机组功频坐标下的一次调频性能估计值, 记为$K$ , 显然$K = \hat \mu $ , 从而得到功频坐标下的机组一次调频性能网格如表2所示, 对浙江某660万千瓦级超临界机组进行原始数据采集和数值迭代模拟, 可得到该机组的一次调频性能三维分布, 如图14所示, 其中, 横轴$f$ 和$P$ 分别为电网频率和机组实际出力值, 纵轴$K$ 表示为机组一次调频性能值.表 2 发电机组一次调频性能网格表Table 2 Grid table of primary frequency control performance of generator sets频率
功率f1 f2 Λ fN−1 fN P1 K11 K12 K1Λ K1(N−1) K1N P2 K21 K22 K2Λ K2(N−1) K2N Λ KΛ1 KΛ2 KΛΛ KΛ(N−1) KΛN PN−1 K(N−1)1 K(N−1)1 K(N−1)Λ K(N−1)(N−1) K(N−1)N PN KN1 K(N−1)1 KNΛ KN(N−1) KNN 
图 14 发电机组一次调频性能功频图谱Fig. 14 Power frequency spectrum of primary frequency control performance of generator sets4. 电网一次调频性能数据分析及估计
4.1 浙江电网一次调频性能分析
2018年10月22日10 : 36时, 浙江省全社会口径用电负荷39753.02 MW, 受电17849.91 MW (占全省用电负荷的52.9%), 发电21903.11 MW (占全省用电负荷的47.1%). 其中, 浙江省调统调机组共并网45台(不含新能源), 其中纯凝汽燃煤机组39台, 背压供热燃煤机组2台, 燃气机组4台, 并网机组平均负荷率约为58.25%, 总发电负荷16221.59 MW, 占当时全省用电负荷的40.8%.
据此, 按照式(11)进行计算, 得到浙江2018年10月22日10 : 36时省调统调并网机组的实际调频能力为12.2327, 根据当时的用电负荷进行折算后约为6 762.23 MW/Hz, 略高于浙江电网ACE系统中自然频率特性系数5 979 MW/Hz[21]. 推此及彼, 采用本文的方法可以得到浙江统调机组每天任意时刻的实际调频能力值, 如图15所示.
图 15 2018年10月22日浙江统调机组实际调频能力时刻图Fig. 15 Time map of actual primary frequency control performance of Zhejiang dispatching unit on October 22, 2018显然, 由图15可以明显得知, 电网的实际调频性能与并网机组台数、并网机组类型、机组平均负荷率等因素息息相关, 是一项电网拓扑结构的综合表征.
4.2 大频差时浙江电网一次调频容量分析
电网频率是描述电力系统发电与用电平衡关系的一个重要变量. 电网中用电负荷的随机特性使电网频率也成为随机变量. 电力负荷的随机变化部分是相互独立、不相关的[22]. 由中心极限定理可知, 连续的独立随机变量叠加后是正态分布的, 因此, 电力系统中负荷的随机性分量也是正态分布的. 在通常情况下, 电网频率的概率分布也就可近似认为是正态分布[23], 因此, 假设发生直流闭锁时的浙江电网系统最低频率近似正态分布, 并结合近三年的历次直流闭锁数据分析, 可得表3数据.
表 3 大频差时浙江电网一次调频数据分析Table 3 Data analysis of primary frequency control data of Zhejiang power grid in large frequency difference situations指标 均值 标准差 95 % 置信区间 最低频率 (Hz) 49.877 0.1022 [49.841, 49.914] 实际出力 (标幺值) 11.9216 6.5165 [9.5856, 14.2575] 出力限值 (%Pe) 3.94 0.90 [3.61, 4.27] 装机容量 (MW) 36 633.71 5 509.92 [34 658.58, 8 608.85] 装机总数 (台数) 68 12 [63, 72] 根据表3, 可知电网频率的最大变化差值为
$$ \Delta f=50-49.841=0.159\;{\rm{Hz}} $$ 最大概率功率缺额为
$$ {P_{{\rm{Loss}}}} = {\beta _{{\rm{ACE}}}} \times \frac{{\Delta f}}{{{f_N}}}=5\;979 \times 0.159=950.66 \;{\rm{MW}} $$ 按浙江历次直流闭锁数据, 可知在不少于63台同步发电机组, 总装机不低于34 658.58 MW的条件下, 浙江统调机组一次调频能力为9.5856 (95%以上), 折算后为6 644.46 MW/Hz, 在最大概率频差0.159 Hz下, 浙江电网可提供一次调频最大能力至少约为1 056.47 MW (3.1% MCR (Maximum continuous rating)), 能够满足功率缺额950.66 MW的要求, 也就是说, 浙江统调机组一次调频能力可以保证电网在直流闭锁等大频差工况下系统频率不低于49.841 Hz.
5. 结束语
机组的转速不等率与电网频率往往呈非线性关系, 采用固定的转速不等率来计算机组的一次调频能力误差较大. 本文基于数据平行驱动体系, 依靠电网运行的大数据, 建立了更为准确的电网一次调频模型, 对提升电网频率控制水平具有非常重要的意义.
-
图 8 MC_TS_UR融合算法在每名受试者数据集上的混淆矩阵
Fig. 8 The CM of the MC_TS_UR fusion algorithm applied to the datasets of every subject
图 9 4种基于MC_ANFIS的算法在训练过程中的分类误差随时间变化曲线
Fig. 9 The classification error changes curve of four MC_ANFIS based algorithms with time during the training process
表 1 参与手势识别实验的受试者信息
Table 1 Information of the subjects participating in the gesture recognition experiment
受试者 性别 年龄 身高 (cm) 体重 (kg) 腕围 (cm) 健康状况 Subject-1 男 25 180.4 72.4 18.8 健康 Subject-2 男 24 169.5 58.5 16.5 健康 Subject-3 男 56 164.6 61.2 15.8 手部截肢 Subject-4 男 25 172.3 62.8 17.9 健康 Subject-5 男 22 177.5 57.0 16.8 健康 Subject-6 男 26 166.6 65.7 18.4 健康 Subject-7 男 23 170.1 73.3 19.1 健康 Subject-8 男 25 175.5 66.9 17.1 健康 
下载: 导出CSV
表 2 6种手腕手势对应的肌肉信息
Table 2 Muscles information of the corresponding six gestures
手势 肌肉 作用 屈腕 尺侧腕屈肌 手腕屈曲和尺侧偏移 握拳 指浅屈肌 手指弯曲 尺侧倾 桡侧腕屈肌 手腕弯曲和径向偏移 伸腕 尺侧腕伸肌 手腕伸展和尺侧偏移 伸掌 指伸肌 手指伸展 桡侧倾 桡侧腕伸肌 手腕伸展和径向偏移
下载: 导出CSV
表 3 8种算法在健康人数据集上的离线实验结果
Table 3 The offline experiment results of eight algorithms on datasets of the normal
指标 SVM GBDT LDA TSK_GD_LSE MC MC_TS MC_UR MC_TS_UR ${\rm{RER}}$ 6.07% 7.82% 5.15% 5.26% 3.16% 2.52% 2.30% 1.82% ${\rm{BER}}$ 6.18% 8.74% 5.21% 5.35% 2.83% 2.41% 2.33% 1.77% $\kappa$ 0.9258 0.9018 0.9375 0.9358 0.9660 0.9711 0.9720 0.9787 $T_t$ 224.6 4.4 0.6 1121.9 796.1 886.5 734.7 310.2
下载: 导出CSV
表 4 8种算法在残疾人数据集上的离线实验结果
Table 4 The offline experiment results of eight algorithms on datasets of the disabled
指标 SVM GBDT LDA TSK_GD_LSE MC MC_TS MC_UR MC_TS_UR ${\rm{RER}}$ 5.94% 8.13% 4.46% 5.77% 4.64% 3.83% 3.77% 2.75% ${\rm{BER}}$ 6.10% 8.27% 4.48% 6.11% 4.72% 3.98% 3.65% 2.73% $\kappa$ 0.9268 0.9008 0.9462 0.9267 0.9434 0.9522 0.9562 0.9672 $T_t$ 173.0 5.3 0.7 1006.5 766.8 942.9 768.6 313.1
下载: 导出CSV
表 5 与近期同类研究工作文献的比较
Table 5 Comparison with similar research work literature
文献 传感器 实验对象 是否为公共数据集 分类算法 手势类别数 识别准确率 [25] 6 通道 pMMG 6 名健康人 否 Fuzzy logic 6 95.30% [26] 8 通道 FMG 10 名健康人 否 SVM 6 93.00% [27] 2 通道 sEMG 7 名健康人 否 SVM 4 95.00% [28] 4 通道 sEMG + 1 通道 IMU 10 名健康人 否 LDA 8 92.60% [29] 8 通道 sEMG 21 名健康人 否 LDA 6 94.70% [30] 8 通道 sEMG 8 名健康人 否 Hidden Markov model 6 94.20% Proposed 6 通道 pMMG + 1 通道 IMU 7 名健康人 否 MC_TS_UR 6 98.18% [31] 8 通道 sEMG 4 名残疾人 否 LDA 7 92.00% [32] 7 通道 sEMG 3 名残疾人 否 SVM 5 94.02% Proposed 6 通道 pMMG + 1 通道 IMU 1 名残疾人 否 MC_TS_UR 6 97.25%
下载: 导出CSV
-
[1] Liu H Y, Wang L H. Gesture recognition for human-robot collaboration: A review. International Journal of Industrial Ergonomics, 2018, 68(1): 355-367 [2] Ding J, Lin R Z, Lin Z Y. Service robot system with integration of wearable Myo armband for specialized hand gesture human–computer interfaces for people with disabilities with mobility problems. Computers Electrical Engineering, 2018, 69(1): 815-827 [3] 丁其川, 赵新刚, 李自由, 韩建达. 基于自更新混合分类模型的肌电运动识别方法. 自动化学报, 2019, 45(8): 1464-1474Ding Q C, Zhao X G, Li Z Y, Han J D. An EMG-motion recognition method with self-update hybrid classification model. Acta Automatica Sinica, 2019, 45(8): 1464-1474 [4] Zhang X, Chen X, Li Y, Lantz V, Wang K, Yang J. A framework for hand gesture recognition based on accelerometer and EMG sensors. IEEE Transactions on Systems, Man, and Cybernetics-Part A: Systems and Humans, 2011, 41(6): 1064-1076 doi: 10.1109/TSMCA.2011.2116004 [5] Duan T, Huang J, Xie Z, Wang L, Xiong C H. Continuous control of wrist-hand prosthesis by extracting independent sEMG signals from cross-talk muscle groups. In: Proceedings of the 2019 Chinese Control Conference. Guangzhou, China: IEEE, 2019. 4537−4542 [6] 李自由, 王丰焱, 赵新刚, 丁其川, 张道辉, 韩建达. 基于 Myo 旋转偏移估计与自适应校正的手势识别方法. 自动化学报, 2020, 46(9): 1896-1907Li Z Y, Wang F Y, Zhao X G, Ding Q C, Zhang D H, Han J D. The method for gestures recognition based on Myo rotation shifts estimation and adaptive correction. Acta Automatica Sinica, 2020, 46(9): 1896-1907 [7] Orizio C, Liberati D, Locatelli C, et al. Surface mechanomyogram reflects muscle fibres twitches summation. Journal of Biomechanics, 1996, 29(4): 475-481 doi: 10.1016/0021-9290(95)00063-1 [8] Akataki K, Mita K, Watakabe M, Itoh K. Mechanomyogram and force relationship during voluntary isometric ramp contractions of the biceps brachii muscle. European Journal of Applied Physiology, 2001, 84(1-2): 19-25 doi: 10.1007/s004210000321 [9] McIntosh J, Marzo A, Fraser M, Phillips C. Echoflex: Hand gesture recognition using ultrasound imaging. In: Proceedings of the 2017 CHI Conference on Human Factors in Computing Systems. Denver, USA: ACM, 2017. 1923−1934 [10] Liu M K, Lin Y T, Qiu Z W, Kuo C K, Wu C K. Hand Gesture Recognition by a MMG-Based Wearable Device. Journal of The Neurological Sciences, 2020, 20(24): 14703-14712 [11] Feng W, Xia C, Zhang Y, Yu J, Jiang W. Research on Chinese sign language recognition methods based on mechanomyogram signals analysis. In: Proceedings of the 4th International Conference on Signal and Image Processing. Wuxi, China: IEEE, 2019. 46−50 [12] Stokes M J, Dalton P A. Acoustic myographic activity increases linearly up to maximal voluntary isometric force in the human quadriceps muscle. Journal of The Neurological Sciences, 1991, 101(2): 163-167 doi: 10.1016/0022-510X(91)90041-5 [13] Ahsan M R, Ibrahimy M I, Khalifa O O. Electromygraphy (EMG) signal based hand gesture recognition using artificial neural network (ANN). In: Proceedings of the 4th International Conference on Mechatronics. Kuala Lumpur, Malaysia: IEEE, 2011. 1−6 [14] Zhang Z, Yang K, Qian J, Zhang L. Real-time surface emg pattern recognition for hand gestures based on an artificial neural network. Sensors, 2019, 19(14): 3170-3184 doi: 10.3390/s19143170 [15] Yao B, Hagras H, Alhaddad M J, Alghazzawi D. A fuzzy logic-based system for the automation of human behavior recognition using machine vision in intelligent environments. Soft Computing, 2015, 19(2): 499-506 doi: 10.1007/s00500-014-1270-4 [16] Hachaj T, Ogiela M R. Rule-based approach to recognizing human body poses and gestures in real time. Multimedia Systems, 2014, 20(1): 81-99 doi: 10.1007/s00530-013-0332-2 [17] Mufarroha F A, Utaminingrum F. Hand gesture recognition using adaptive network based fuzzy inference system and k-nearest neighbor. International Journal of Technology, 2017, 8(3): 559-567 doi: 10.14716/ijtech.v8i3.3146 [18] Khezri M, Jahed M. A neuro–fuzzy inference system for sEMG-based identification of hand motion commands. IEEE Transactions on Industrial Electronics, 2010, 58(5): 1952-60 [19] Hill A V. The Heat of Shortening and the Dynamic Constants of Muscle. Proceedings of the Royal Society of London, 1938, 126(843): 136-195 [20] Belyea A, Englehart K, Scheme E. FMG Versus EMG: A comparison of usability for real-time pattern recognition based control. IEEE Transactions on Biomedical Engineering, 2019, 66(11): 3098-3104 doi: 10.1109/TBME.2019.2900415 [21] Peng X Y, Li L, Wang F Y. Accelerating minibatch stochastic gradient descent using typicality sampling. IEEE Transactions on Neural Networks and Learning Systems, 2019, 31(11): 4649-4659 [22] Cui Y Q, Wu D R, Huang J. Optimize TSK fuzzy systems for classification problems: Mini-batch gradient descent with uniform regularization and batch normalization. IEEE Transactions on Fuzzy Systems, 2020, 28(12): 3065-3075 doi: 10.1109/TFUZZ.2020.2967282 [23] You H, Ma Z, Tang Y, Wang Y, Yan J, Ni M, Cen K, Huang Q. Comparison of ANN (MLP), ANFIS, SVM, and RF models for the online classification of heating value of burning municipal solid waste in circulating fluidized bed incinerators. Waste Management, 2017, 68(10): 186-197 [24] Shariati S, Haghighi M M. Comparison of anfis neural network with several other ANNs and support vector the machine for diagnosing hepatitis and thyroid diseases. In: Proceedings of the 2010 International Conference on Computer Information Systems and Industrial Management Applications. Krackow, Poland: IEEE, 2010. 596−599 [25] Jung P G, Lim G, Kim S, Kong K. A wearable gesture recognition device for detecting muscular activities based on air-pressure sensors. IEEE Transactions on Industrial Informatics, 2015, 11(2): 485-494 [26] Anvaripour M, Saif M. Hand gesture recognition using force myography of the forearm activities and optimized features. In: Proceedings of the 2018 IEEE International Conference on Industrial Technology. Lyon, France: IEEE, 2018. 187−192 [27] Tavakoli M, Benussi C, Lopes P A, et al. Robust hand gesture recognition with a double channel surface EMG wearable armband and SVM classifier. Biomedical Signal Processing and Control, 2018, 46(1): 121-130 [28] Jiang S, Lv B, Guo W, Zhang C, Wang H, Sheng X, Shull P B. Feasibility of wrist-worn, real-time hand, and surface gesture recognition via sEMG and IMU Sensing. IEEE Transactions on Industrial Informatics, 2017, 14(8): 3376-3385 [29] Botros F, Phinyomark A, Scheme E. EMG-based gesture recognition: Is it time to change focus from the forearm to the wrist?. IEEE Transactions on Industrial Informatics, 2022, 18(1), 174−184 [30] Shaabana A, Legere J, Li J, Zheng R, Mohrenschildt MV, Shedden JM. Portable electromyography: A case study on ballistic finger movement recognition. IEEE Sensors Journal, 2019, 19(16): 7043-55 doi: 10.1109/JSEN.2019.2908312 [31] Powell M A, Kaliki R R, Thakor N V. User training for pattern recognition-based myoelectric prostheses: Improving phantom limb movement consistency and distinguishability. IEEE Transactions on Neural Systems and Rehabilitation Engineering, 2014, 22(3): 522-532 doi: 10.1109/TNSRE.2013.2279737 [32] Kartsch V, Benatti S, Mancini M, Magno M, Benini L. Smart wearable wristband for EMG based gesture recognition powered by solar energy harvester. In: Proceedings of the 2018 IEEE International Symposium on Circuits and Systems. Florence, Italy: IEEE, 2018. 1−5 期刊类型引用(9)
1. 任姗. 基于动态手势识别算法的VSTi系统研究. 自动化与仪器仪表. 2025(01): 214-217+222 . 
百度学术2. 李婷婷，王靖，骆亚丽，刘红梅. 基于Leap Motion传感器的弹琴触键手势自动控制系统设计. 自动化与仪器仪表. 2025(02): 223-227 . 百度学术3. 邹灵果，张美花. 基于数理统计特征的人机交互图像手势识别. 黑龙江工业学院学报(综合版). 2024(01): 97-104 . 百度学术4. 边宝丽. 基于手势识别的幼儿游戏机器系统设计. 自动化与仪器仪表. 2024(03): 171-174+179 . 百度学术5. 闫颢月，王伟，田泽. 复杂环境下基于改进YOLOv5的手势识别方法. 计算机工程与应用. 2023(04): 224-234 . 百度学术6. 张琳钦. 基于支持向量机的人机交互媒体播放界面手势识别方法. 西安航空学院学报. 2023(03): 83-88 . 百度学术7. 盛博莹，侯进，李嘉新，党辉. 面向复杂交通场景的道路目标检测方法. 计算机工程与应用. 2023(15): 87-96 . 百度学术8. 李俊文，张红英，韩宾. 深层特征聚合引导的轻量级显著性目标检测. 计算机工程与应用. 2023(19): 122-129 . 百度学术9. 王鹏飞，黄汉明，王梦琪. 改进YOLOv5的复杂道路目标检测算法. 计算机工程与应用. 2022(17): 81-92 . 百度学术其他类型引用(4)
-
AAS-CN-2020-0901手势数据.zip
-
下载:
计量
- 文章访问数: 1483
- HTML全文浏览量: 817
- PDF下载量: 228
- 被引次数: 13
