-
摘要: 以小型多轴无人机为代表的“低慢小”目标, 通常难以被常规手段探测, 而此类目标又会严重威胁某些重要设施. 因此对该类目标的识别已经成为一个亟待解决的重要问题. 本文基于目标运动特征, 提出了一种无人机目标识别方法, 并揭示了二阶运动参量以及重力方向运动参量是无人机识别过程中的关键参数. 该方法首先提取候选目标的多阶运动参量, 建立梯度提升树(Gradient boosting decision tree, GBDT)和门控制循环单元(Gate recurrent unit, GRU)记忆神经网络分别完成短时和长期识别, 然后融合表观特征识别结果得到最终判别结果. 此外, 本文还建立了一个综合多尺度无人机数据集(Multi-scale UAV dataset, MUD), 本文所提出的方法在该数据集上相对于传统基于运动特征的方法, 其识别精度(Average precision, AP)提升103%, 融合方法提升26%.Abstract: Due to the features of low, slow and small aircraft, such as quadrotors, it is a challenging and urgent problem to detect UAVs (Unmanned aerial vehicles) in the wild. Different from the past literatures directly using deep learning method, this paper exploits motion features by extracting multi-order kinematic parameters such as velocity, accelerate, angular velocity, angular velocity vectors and it is exposed that 2nd order and gravity direction motion parameters are key motion patterns for UAV detection. By building GBDT (Gradient boosting decision tree) and GRU (Gate recurrent unit) network, it comes out with a short-term and a long-term detection result, respectively. This recognition process integrates appearance detection result into motion detection result and obtains the final determination. The experimental results achieve state-of-the-art result, with a 103% increase on the precision index AP (Average precision) with respect to the previous work and a 26% increase for hybrid method.
-
Key words:
- Quadrotors /
- object detection /
- motion feature /
- fusion method
-
电网频率是电力系统运行中最重要的关键参数之一, 系统频率的异常波动不仅影响电网的电能质量, 制约电力设备的正常使用, 严重的频率失衡甚至将导致并网设备脱网及系统崩溃, 对电力系统的稳定运行带来严重危害[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. 发电机组一次调频数据采集体系
当前, 电力系统已经发展成一个具有多源信息融合的综合性电力控制系统. 根据现有的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的网源协调信息管理系统, 通过设置统一的一次调频试验台账, 维护和存储机组的历次原始试验数据和参数, 统一存储于数据库中.
2.1 机组一次调频在线测试系统
如图4所示, 在D5000系统中部署试验触发、结果展示、参数维护及统计分析功能. 试验触发功能可实现选定机组进行试验, 试验开始时将测试机组投入一次调频扰动测试模式, 调控中心向电厂侧下发调频测试开始信号和频率偏差信号, 机组根据接收的信号进行一次调频动作响应, 具体控制信息流程见图5所示.
在一次调频在线测试期间, 系统自动记录测试机组的有功PMU数据, 整编录入D5000历史库, 永久保存, 供后期调频指标分析与数据展示使用[17]. 在结果展示功能中, 将频率扰动信号及电厂侧PMU装置上传的有功出力信号进行展示, 并根据有功出力曲线自动计算一次调频动作性能指标. 在浙江某厂#2机组进行一次调频在线测试系统联调试验时, 调控中心主站侧的操作及结果展示如图6所示, 发电厂机组一次调频动作响应曲线如图7所示.
通过PMU实时采集高频度的功率、相位、相角等机组数据, 并利用PMU的GPS对时优势, 计算机组在一次调频在线测试中不同时刻、不同工况下的性能指标, 反向推算出机组实际的调频特性参数, 以此形成海量的机组一次调频特征数据, 实现机组一次调频静态参数、动态性能等特征数据的分析及管理[19].
2.2 网源协调信息管理系统
基于OMS的网源协调信息管理系统(如图8所示)能够实现设备台账、电气参数、涉网试验报告等技术资料的及时更新及有效管理, 打通厂网之间的涉网信息共享交互渠道, 更加科学、系统、高效地开展网源协调管理工作. 以各机组台账信息为基础, 实现对网内各并网机组数据信息的实时统计汇总; 以并网机组各类信息关键字为基础, 实现对特定类型、特定机组相关信息的快速查询和准确统计, 该系统所部署的功能模块如图9所示. 同时, 建立网源协调数据中心, 这是网源协调信息管理平台的核心内容, 可实现参数模型库、文档资料库、督办信息库及信息变更库的数据获取与整合、数据质量校核、合格率督办和数据查询等服务.
在网源协调信息管理系统中, 通过设置统一的一次调频试验台账模板, 采集机组的历次试验数据, 以维护和存储机组的原始试验数据. 图10为浙江某百万千瓦机组A修后一次调频试验的数据台账.
根据设置统一的数据台账, 可以追溯机组每一次试验工况, 分析不同工况下机组一次调频试验时的真实性能.
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所示), 依靠电网运行的大数据, 通过不断演算、迭代优化, 形成调频能力的概率区间来衡量电网调频性能, 解决了电网一次调频能力的在线估计问题, 从而使电网的运行状态具有更高的可控性及透明度, 提高大电网应对负荷突变的能力.
3.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 }$ 为电网负载频率特性系数.$${\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系统中的多源一次调频数据进行有效整合, 以获取机组真实的一次调频性能功频特性图谱.
可以从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 4. 电网一次调频性能数据分析及估计
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可以明显得知, 电网的实际调频性能与并网机组台数、并网机组类型、机组平均负荷率等因素息息相关, 是一项电网拓扑结构的综合表征.
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. 结束语
机组的转速不等率与电网频率往往呈非线性关系, 采用固定的转速不等率来计算机组的一次调频能力误差较大. 本文基于数据平行驱动体系, 依靠电网运行的大数据, 建立了更为准确的电网一次调频模型, 对提升电网频率控制水平具有非常重要的意义.
-
表 1 本文所采集数据与其他运动目标数据集的对比
Table 1 Comparison of different datasets for moving objects
表 2 MUD数据集采集设备说明
Table 2 Main equipment for acquisition of multi-scale UAV dataset (MUD)
设备 参数 精度 相机 SONY A7 ILCE-7M2, $6\,000 \times 4\,000$像素FE 24 ~ 240 mm, F 3.5 ~ 6.3 — GPS GPS/GLONASS双模 垂直$\pm 0.5\;{\rm{m} },$ 水平$ \pm 1.5\;{\rm{m}}$ 激光测距仪 SKIL Xact 0530, 0 ~ 80 m $ \pm 0.2\;{\rm{mm}}$ 表 3 运动目标区域提取算法性能对比
Table 3 Comparison between performance of different motion ROIs
表 4 不同深度估计方法误差对比
Table 4 Error of different depth estimation methods
方法 探测范围 绝对误差 平方误差 均方根误差 $\delta < 1.25$ $\delta < {1.25^2}$ $\delta < {1.25^3}$ DORN[51] 0 ~ 100 m 0.103 0.321 9.014 0.832 0.875 0.922 GeoNet[63] 0 ~ 100 m 0.280 2.813 14.312 0.817 0.849 0.895 双目视觉[64] 0 ~ 100 m 0.062 1.210 0.821 0.573 0.642 0.692 激光测距 0 ~ 200 m 0.041 2.452 1.206 0.875 0.932 0.961 DORN[51] 200 ~ 500 m 0.216 1.152 13.021 0.672 0.711 0.748 GeoNet[63] 200 ~ 500 m 0.398 5.813 18.312 0.617 0.649 0.696 双目视觉[64] 200 ~ 500 m 0.786 5.210 25.821 0.493 0.532 0.562 激光测距 200 ~ 500 m 0.078 3.152 2.611 0.891 0.918 0.935 参数 说明 ${\boldsymbol v}$ 速度 ${\boldsymbol a}$ 加速度 ${\boldsymbol \omega}$ 角速度 ${\boldsymbol \alpha }$ 角加速度 X 轴分量方向 与图像平面坐标系中 u 轴方向保持一致 Y 轴分量方向 与图像平面坐标系中 v 轴方向保持一致 Z 轴分量方向 铅垂向上 表 6 运动参量的决策树模型识别结果混淆矩阵
Table 6 Confusion matrix of MokiP by using GDBT
真实值
预测值旋翼无人机 鸟类 行人 车辆 其他物体 旋翼无人机 0.67 0.25 0.02 0.01 0.12 鸟类 0.21 0.58 0.01 0.00 0.10 行人 0.01 0.02 0.75 0.06 0.09 车辆 0.01 0.00 0.10 0.80 0.08 其他物体 0.10 0.15 0.12 0.13 0.61 表 7 不同识别方法性能指标对比表
Table 7 Comparison of performance indexes for different detection method
表 8 运动参量的性质对无人机识别的影响表
Table 8 Impact of the parameter properties on UAV detection
参量贡献度$D$ 平动参量 旋转参量 总贡献度 一阶参量 7.2% 20.1% 27.3% 二阶参量 34.1% 38.6% 72.7% 总贡献度 41.3% 58.7% 1 表 9 运动参量的方向对无人机识别的影响表
Table 9 Impact of the parameter direction on UAV detection
参量贡献度$D$ 沿 X 轴方向 沿 Y 轴方向 沿 Z 轴方向 总贡献度 平动参量 8.3% 8.8% 24.2% 41.3% 旋转参量 18.7% 18.8% 22.2% 58.7% 总贡献度 27.0% 27.6% 46.4% 1 -
[1] 李菠, 孟立凡, 李晶, 刘春美, 黄广炎. 低空慢速小目标探测与定位技术研究. 中国测试, 2016, 42(12): 64-69 doi: 10.11857/j.issn.1674-5124.2016.12.014Li Bo, Meng Li-Fan, Li Jing, Liu Chun-Mei, Huang Guang-Yan. Research on detecting and locating technology of LSS-UAV. China Measurement & Test, 2016, 42(12): 64-69 doi: 10.11857/j.issn.1674-5124.2016.12.014 [2] Wang Z H, Lin X P, Xiang X Y, Blasch E, Pham K, Chen G S, et al. An airborne low SWaP-C UAS sense and avoid system. In: Proceedings of SPIE 9838, Sensors and Systems for Space Applications IX. Baltimore, USA: SPIE, 2016. 98380C [3] Busset J, Perrodin F, Wellig P, Ott B, Heutschi K, Rühl T, et al. Detection and tracking of drones using advanced acoustic cameras. In: Proceedings of SPIE 9647, Unmanned/Unattended Sensors and Sensor Networks XI; and Advanced Free-Space Optical Communication Techniques and Applications. Toulouse, France: SPIE, 2015. 96470F [4] Mezei J, Fiaska V, Molnár A. Drone sound detection. In: Proceedings of the 16th IEEE International Symposium on Computational Intelligence and Informatics (CINTI). Budapest, Hungary: IEEE, 2015. 333−338 [5] 张号逵, 李映, 姜晔楠. 深度学习在高光谱图像分类领域的研究现状与展望. 自动化学报, 2018, 44(6): 961-977Zhang Hao-Kui, Li Ying, Jiang Ye-Nan. Deep learning for hyperspectral imagery classification: The state of the art and prospects. Acta Automatica Sinica, 2018, 44(6): 961-977 [6] 贺霖, 潘泉, 邸, 李远清. 高光谱图像高维多尺度自回归有监督检测. 自动化学报, 2009, 35(5): 509-518He Lin, Pan Quan, Di Wei, Li Yuan-Qing. Supervised detection for hyperspectral imagery based on high-dimensional multiscale autoregression. Acta Automatica Sinica, 2009, 35(5): 509-518 [7] 叶钰, 王正, 梁超, 韩镇, 陈军, 胡瑞敏. 多源数据行人重识别研究综述. 自动化学报, 2020, 46(9): 1869-1884Ye Yu, Wang Zheng, Liang Chao, Han Zhen, Chen Jun, Hu Rui-Min. A survey on multi-source person re-identification. Acta Automatica Sinica, 2020, 46(9): 1869-1884 [8] Zhao J F, Feng H J, Xu Z H, Li Q, Peng H. Real-time automatic small target detection using saliency extraction and morphological theory. Optics & Laser Technology, 2013, 47: 268-277 [9] Nguyen P, Ravindranatha M, Nguyen A, Han R, Vu T. Investigating cost-effective RF-based detection of drones. In: Proceedings of the 2nd Workshop on Micro Aerial Vehicle Networks, Systems, and Applications for Civilian Use. Singapore: Association for Computing Machinery, 2016. 17−22 [10] Drozdowicz J, Wielgo M, Samczynski P, Kulpa K, Krzonkalla J, Mordzonek M, et al. 35 GHz FMCW drone detection system. In: Proceedings of the 17th International Radar Symposium. Krakow, Poland: IEEE, 2016. 1−4 [11] Felzenszwalb P F, Girshick R B, McAllester D, Ramanan D. Object detection with discriminatively trained part-based models. IEEE Transactions on Pattern Analysis and Machine Intelligence, 2010, 32(9): 1627-1645 doi: 10.1109/TPAMI.2009.167 [12] Ren S Q, He K M, Girshick R, Sun J. Faster R-CNN: Towards real-time object detection with region proposal networks. IEEE Transactions on Pattern Analysis and Machine Intelligence, 2017, 39(6): 1137-1149 doi: 10.1109/TPAMI.2016.2577031 [13] Dollar P, Tu Z W, Perona P, Belongie S. Integral channel features. In: Proceedings of the 2009 British Machine Vision Conference. London, UK: BMVA Press, 2009. 91.1−91.11 [14] Aker C, Kalkan S. Using deep networks for drone detection. In: Proceedings of the 14th IEEE International Conference on Advanced Video and Signal Based Surveillance (AVSS). Lecce, Italy: IEEE, 2017. 1−6 [15] Rozantsev A, Lepetit V, Fua P. Flying objects detection from a single moving camera. In: Proceedings of the 2015 IEEE Conference on Computer Vision and Pattern Recognition (CVPR). Boston, USA: IEEE, 2015. 4128−4136 [16] Coluccia A, Ghenescu M, Piatrik T, De Cubber G, Schumann A, Sommer L, et al. Drone-vs-Bird detection challenge at IEEE AVSS2017. In: Proceedings of the 14th IEEE International Conference on Advanced Video and Signal Based Surveillance (AVSS). Lecce, Italy: IEEE, 2017. 1−6 [17] Schumann A, Sommer L, Klatte J, Schuchert T, Beyerer J. Deep cross-domain flying object classification for robust UAV detection. In: Proceedings of the 14th IEEE International Conference on Advanced Video and Signal Based Surveillance (AVSS). Lecce, Italy: IEEE, 2017. 1−6 [18] Sommer L, Schumann A, Müller T, Schuchert T, Beyerer J. Flying object detection for automatic UAV recognition. In: Proceedings of the 14th IEEE International Conference on Advanced Video and Signal Based Surveillance (AVSS). Lecce, Italy: IEEE, 2017. 1−6 [19] Sapkota K R, Roelofsen S, Rozantsev A, Lepetit V, Gillet D, Fua P, et al. Vision-based unmanned aerial vehicle detection and tracking for sense and avoid systems. In: Proceedings of the 2016 IEEE/RSJ International Conference on Intelligent Robots and Systems (IROS). Daejeon, Korea: IEEE, 2016. 1556−1561 [20] Carrio A, Vemprala S, Ripoll A, Saripall S, Campoy P. Drone detection using depth maps. In: Proceedings of the 2018 IEEE/RSJ International Conference on Intelligent Robots and Systems (IROS). Madrid, Spain: IEEE, 2018. 1034−1037 [21] Carrio A, Tordesillas J, Vemprala S, Saripalli S, Campoy P, How J P. Onboard detection and localization of drones using depth maps. IEEE Access, 2020, 8: 30480-30490 doi: 10.1109/ACCESS.2020.2971938 [22] Ganti S R, Kim Y. Implementation of detection and tracking mechanism for small UAS. In: Proceedings of the 2016 International Conference on Unmanned Aircraft Systems (ICUAS). Arlington, USA: IEEE, 2016. 1254−1260 [23] Farhadi M, Amandi R. Drone detection using combined motion and shape features. In: IEEE International Workshop on Small-Drone Surveillance Detection and Counteraction Techniques. Lecce, Italy: IEEE, 2017. 1−6 [24] Alom M Z, Hasan M, Yakopcic C, Taha T M, Asari V K. Improved inception-residual convolutional neural network for object recognition. Neural Computing and Applications, 2020, 32(1): 279-293 doi: 10.1007/s00521-018-3627-6 [25] Saqib M, Khan S D, Sharma N, Blumenstein M. A study on detecting drones using deep convolutional neural networks. In: Proceedings of the 14th IEEE International Conference on Advanced Video and Signal Based Surveillance (AVSS). Lecce, Italy: IEEE, 2017. 1−5 [26] Liu W, Anguelov D, Erhan D, Szegedy C, Reed S, Fu C Y, et al. SSD: Single shot MultiBox detector. In: Proceedings of the 14th European Conference on Computer Vision - ECCV 2016. Amsterdam, The Netherlands: Springer, 2016. 21−37. [27] 周卫祥, 孙德宝, 彭嘉雄. 红外图像序列运动小目标检测的预处理算法研究. 国防科技大学学报, 1999, 21(5): 60-63Zhou Wei-Xiang, Sun De-Bao, Peng Jia-Xiong. The study of preprocessing algorithm of small moving target detection in infrared image sequences. Journal of National University of Defense Technology, 1999, 21(5): 60-63 [28] Wu Y W, Sui Y, Wang G H. Vision-based real-time aerial object localization and tracking for UAV sensing system. IEEE Access, 2017, 5: 23969-23978 doi: 10.1109/ACCESS.2017.2764419 [29] Lv P Y, Lin C Q, Sun S L. Dim small moving target detection and tracking method based on spatial-temporal joint processing model. Infrared Physics & Technology, 2019, 102: Article No. 102973 [30] Van Droogenbroeck M, Paquot O. Background subtraction: Experiments and improvements for ViBe. In: Proceedings of the 2012 IEEE Computer Society Conference on Computer Vision and Pattern Recognition Workshops. Providence, USA: IEEE, 2012. 32−37 [31] Zamalieva D, Yilmaz A. Background subtraction for the moving camera: A geometric approach. Computer Vision and Image Understanding, 2014, 127: 73-85 doi: 10.1016/j.cviu.2014.06.007 [32] Sun Y F, Liu G, Xie L. MaxFlow: A convolutional neural network based optical flow algorithm for large displacement estimation. In: Proceedings of the 17th International Symposium on Distributed Computing and Applications for Business Engineering and Science (DCABES). Wuxi, China: IEEE, 2018. 119−122 [33] Dosovitskiy A, Fischer P, Ilg E, Häusser P, Hazirbas C, Golkov V, et al. FlowNet: Learning optical flow with convolutional networks. In: Proceedings of the 2015 International Conference on Computer Vision (ICCV). Santiago, Chile: IEEE, 2015. 2758−2766 [34] Ilg E, Mayer N, Saikia T, Keuper M, Dosovitskiy A, Brox T. FlowNet 2.0: Evolution of optical flow estimation with deep networks. In: Proceedings of the 2017 IEEE Conference on Computer Vision and Pattern Recognition (CVPR). Honolulu, USA: IEEE, 2017. 1647−1655 [35] 陈鑫, 魏海军, 吴敏, 曹卫华. 基于高斯回归的连续空间多智能体跟踪学习. 自动化学报, 2013, 39(12): 2021-2031Chen Xin, Wei Hai-Jun, Wu Min, Cao Wei-Hua. Tracking learning based on gaussian regression for multi-agent systems in continuous space. Acta Automatica Sinica, 2013, 39(12): 2021-2031 [36] Shi S N, Shui P L. Detection of low-velocity and floating small targets in sea clutter via income-reference particle filters. Signal Processing, 2018, 148: 78-90 doi: 10.1016/j.sigpro.2018.02.005 [37] Kang K, Li H S, Yan J J, Zeng X Y, Yang B, Xiao T, et al. T-CNN: Tubelets with convolutional neural networks for object detection from videos. IEEE Transactions on Circuits and Systems for Video Technology, 2018, 28(10): 2896-2907 doi: 10.1109/TCSVT.2017.2736553 [38] Zhu X Z, Xiong Y W, Dai J F, Yuan L, Wei Y C. Deep feature flow for video recognition. In: Proceedings of the 2017 IEEE Conference on Computer Vision and Pattern Recognition (CVPR). Honolulu, USA: IEEE, 2017. 4141−4150 [39] Zhu X Z, Wang Y J, Dai J F, Yuan L, Wei Y C. Flow-guided feature aggregation for video object detection. In: Proceedings of the 2017 IEEE International Conference on Computer Vision (ICCV). Venice, Italy: IEEE, 2017. 408−417 [40] Bertasius G, Torresani L, Shi J B. Object detection in video with spatiotemporal sampling networks. In: Proceedings of the 15th European Conference on Computer Vision. Munich, Germany: Springer, 2018. 342−357 [41] Luo H, Huang L C, Shen H, Li Y, Huang C, Wang X G. Object detection in video with spatial-temporal context aggregation. arXiv: 1907.04988, 2019. [42] Xiao F Y, Lee Y J. Video object detection with an aligned spatial-temporal memory. In: Proceedings of the 15th European Conference on Computer Vision. Munich, Germany: Springer, 2018. 494−510. [43] Chen X Y, Yu J Z, Wu Z X. Temporally identity-aware SSD with attentional LSTM. IEEE Transactions on Cybernetics, 2020, 50(6): 2674-2686. doi: 10.1109/TCYB.2019.2894261 [44] Shi X J, Chen Z R, Wang H, Yeung D Y, Wong W K, Woo W C. Convolutional LSTM network: A machine learning approach for precipitation nowcasting. In: Proceedings of the 28th International Conference on Neural Information Processing Systems. Montreal, Canada: MIT Press, 2015. 802−810. [45] 高雪琴, 刘刚, 肖刚, Bavirisetti D P, 史凯磊. 基于FPDE的红外与可见光图像融合算法. 自动化学报, 2020, 46(4): 796-804.Gao Xue-Qin, Liu Gang, Xiao Gang, Bavirisetti Durga Prasad, Shi Kai-Lei. Fusion Algorithm of Infrared and Visible Images Based on FPDE. Acta Automatica Sinica, 2020, 46(4): 796-804 [46] Bluche T, Messina R. Gated convolutional recurrent neural networks for multilingual handwriting recognition. In: Proceedings of the 14th IAPR International Conference on Document Analysis and Recognition (ICDAR). Kyoto, Japan: IEEE, 2017. 646−651 [47] Deng J, Dong W, Socher R, Li L J, Li K, Li F F. ImageNet: A large-scale hierarchical image database. In: Proceedings of the 2009 IEEE Conference on Computer Vision and Pattern Recognition. Miami, USA: IEEE, 2009. 248−255 [48] Son J, Jung I, Park K, Han B. Tracking-by-segmentation with online gradient boosting decision tree. In: Proceedings of the 2015 IEEE International Conference on Computer Vision (ICCV). Santiago, Chile: IEEE, 2015. 3056−3064 [49] Wang F, Jiang M Q, Qian C, Yang S, Li C, Zhang H G, et al. Residual attention network for image classification. In: Proceedings of the 2017 IEEE Conference on Computer Vision and Pattern Recognition (CVPR). Honolulu, USA: IEEE, 2017. 6450−6458 [50] 张秀伟, 张艳宁, 郭哲, 赵静, 仝小敏. 可见光-热红外视频运动目标融合检测的研究进展及展望. 红外与毫米波学报, 2011, 30(4): 354-360Zhang Xiu-Wei, Zhang Yan-Ning, Guo Zhe, Zhao Jing, Tong Xiao-Min. Advances and perspective on motion detection fusion in visual and thermal framework. Journal of Infrared and Millimeter Waves, 2011, 30(4): 354-360 [51] Fu H, Gong M M, Wang C H, Batmanghelich K, Tao D C. Deep ordinal regression network for monocular depth estimation. In: Proceedings of the 2018 IEEE/CVF Conference on Computer Vision and Pattern Recognition. Salt Lake City, USA: IEEE, 2018. 2002−2011 [52] Dragon R, van Gool L. Ground plane estimation using a hidden Markov model. In: Proceedings of the 2014 IEEE Conference on Computer Vision and Pattern Recognition. Columbus, USA: IEEE, 2014. 4026−4033 [53] Rublee E, Rabaud V, Konolige K, Bradski G. ORB: An efficient alternative to SIFT or SURF. In: Proceedings of the 2011 International Conference on Computer Vision. Barcelona. Spain: IEEE, 2011. 2564−2571 [54] Lepetit V, Moreno-Noguer F, Fua P. EPnP: An accurate O(n) solution to the PnP problem. International Journal of Computer Vision, 2009, 81(2): 155-166 doi: 10.1007/s11263-008-0152-6 [55] Geiger A, Lenz P, Urtasun R. Are we ready for autonomous driving? The KITTI vision benchmark suite. In: Proceedings of the 2012 IEEE Conference on Computer Vision and Pattern Recognition. Providence, USA: IEEE, 2012. 3354−3361 [56] Bagautdinov T, Fleuret F, Fua P. Probability occupancy maps for occluded depth images. In: Proceedings of the 2015 IEEE Conference on Computer Vision and Pattern Recognition (CVPR). Boston, USA: IEEE, 2015. 2829−2837 [57] Kranstauber B, Cameron A, Weinzerl R, Fountain T, Tilak S, Wikelski M, et al. The Movebank data model for animal tracking. Environmental Modelling & Software, 2011, 26(6): 834-835 [58] Belongie S, Perona P, Van Horn G, Branson S. NABirds dataset: Download it now! [Online], available: https://dl.allaboutbirds.org/nabirds, March 30, 2020 [59] Xiang Y, Mottaghi R, Savarese S. Beyond PASCAL: A benchmark for 3D object detection in the wild. In: Proceedings of the 2014 IEEE Winter Conference on Applications of Computer Vision. Steamboat Springs, USA: IEEE, 2014. 75−82 [60] Silberman N, Hoiem D, Kohli P, Fergus R. Indoor segmentation and support inference from RGBD images. In: Proceedings of the 12th European conference on Computer Vision. Florence, Italy: Springer, 2012. 746−760 [61] Lucas B D, Kanade T. An iterative image registration technique with an application to stereo vision. In: Proceedings of the 7th International Joint Conference on Artificial Intelligence. Vancouver, Canada: Morgan Kaufmann, 1981. 674−679 [62] Yazdian-Dehkordi M, Rojhani O R, Azimifar Z. Visual target tracking in occlusion condition: A GM-PHD-based approach. In: Proceedings of the 16th CSI International Symposium on Artificial Intelligence and Signal Processing (AISP 2012). Shiraz, Iran: IEEE, 2012. 538−541 [63] Yin Z C, Shi J P. GeoNet: Unsupervised learning of dense depth, optical flow and camera pose. In: Proceedings of the 2018 IEEE/CVF Conference on Computer Vision and Pattern Recognition. Salt Lake City, USA: IEEE, 2018. 1983−1992 [64] Mur-Artal R, Tardós J D. ORB-SLAM2: An open-source SLAM system for monocular, stereo, and RGB-D cameras. IEEE Transactions on Robotics, 2017, 33(5): 1255-1262 doi: 10.1109/TRO.2017.2705103 期刊类型引用(5)
1. 曹强,党晓圆,高俊,南亮生,吴胜磊,曾庆欣,张杰. 模型预测控制理论的研究与应用. 自动化应用. 2024(09): 8-14+18 . 百度学术
2. 马乐乐,刘向杰,高福荣. 基于知识迁移的数据驱动迭代学习模型预测控制. 中国科学:信息科学. 2024(07): 1752-1774 . 百度学术
3. 孔小兵,王文文,刘向杰. 风电场分布式经济模型预测控制. 中国科学:信息科学. 2024(09): 2200-2223 . 百度学术
4. 何思宇,鲁浩林,冯朝,肖晓晖. 基于平行度误差补偿的多通道双晶单色器运动控制. 光学精密工程. 2023(09): 1335-1346 . 百度学术
5. 申永鹏,袁小芳,赵素娜,孟步敏,王耀南. 智能网联电动汽车节能优化控制研究进展与展望. 自动化学报. 2023(12): 2437-2456 . 本站查看
其他类型引用(23)
-