张钦培 李 健 卢 阳 吴凌豪 杨 凯 孙佳伟
低载波比牵引系统的感应电机特征根离散化模型研究
张钦培1李 健1卢 阳1吴凌豪2杨 凯1孙佳伟3
(1. 强电磁工程与新技术国家重点实验室(华中科技大学电气与电子工程学院) 武汉 430074 2. 华中科技大学电气与电子工程学院新型电机与特种电磁装备教育部工程研究中心 武汉 430074 3. 中车大连电力牵引研发中心有限公司 大连 116052)
在大功率和高速电机驱动领域,电机控制系统将运行于低载波比工况。传统的一阶欧拉、二阶双线性等降阶离散化模型在低载波比下由于离散化误差过大,对应的状态观测将出现幅值和相位的稳态误差,严重时甚至出现发散不收敛现象。针对上述问题,该文提出了感应电机特征根离散化模型。通过构建感应电机的复矢量模型状态空间方程,将满秩的状态转移矩阵进行对角化,得到状态转移矩阵的精确离散化结果,该模型在低载波比时仍具有较高的离散化精度。同时,提出了一种基于伯德图的离散化误差定量分析方法,通过定量对比不同离散化模型和连续域模型之间观测变量的幅值和相位误差,从理论上证明了提出方法的优越性。最后,通过仿真和实验验证了上述感应电机特征根离散化模型在低载波比下均具有良好的稳态精度与暂态跟随性能。
感应电机 离散化模型 传递函数 伯德图 低载波比
在电机驱动系统中,电机离散模型的精度在状态观测器设计、无位置/无速度传感器算法、电流控制器设计等方面均具有重要意义。离散模型的准确度依赖于电磁关系的建模、参数的准确度、离散化误差的大小等几个重要方面。其中,参数的离线测量和在线辨识已有国内外学者对其进行了大量的研究工作,且目前仍在继续深入[1-5]。而在离散化方面,随着高速电机、大功率牵引电机的研究持续深入,离散化误差对控制系统的影响被越来越多学者所关注[6-8]。
在大功率牵引电机驱动系统中,当电机进入中高速区域时控制系统的开关频率与电频率的比值(“载波比”)将逐渐下降,随着高性能同步调制策略的应用,载波比将低至3,最后甚至进入方波工况。在低载波比工况下,传统降阶的离散化模型由于离散化误差过大,无法适用于低载波工况的控制系统[9]。文献[10]从状态矩阵的角度分析了一阶欧拉、二阶双线性变换的离散化误差,但未在低载波比工况下进行分析。文献[11]提出变坐标系的离散化方法,将定子方程在静止坐标系下离散化,转子方程在转子旋转坐标系下离散化,大大降低了离散化误差,然而在极低载波比下仍然无法满足高性能要求。文献[12]对比了一阶欧拉、二阶龙格库塔法、双线性变换离散化方法对无速度传感器算法的影响,分析了弱磁区域的无速度观测器稳态精度及动态性能。文献[13]研究了一种基于电流源的永磁同步电机预测转矩控制,但仅采用一阶欧拉离散化方法,没有在低载波比情况下研究其离散化误差。文献[14-15]提出在连续域下采用复数域模型进行电流控制器设计,该方法未考虑离散化误差的影响,在载波比较低时性能较差。因此,文献[16]提出直接在离散域下进行建模,提升了高速下的解耦性能,然而受限于模型精度的影响,极低载波比下的性能提升依旧困难。文献[17-18]研究了低载波比下永磁同步电机离散模型的近似误差,文献[8, 19]在此基础上研究了考虑转子位置补偿的离散模型,该模型对预测电流控制器性能有显著改善。文献[20]分析了感应电机混合型磁链观测器在离散域内的实现问题,通过考虑转子偏转角,提出了改进的数字实现方案。文献[21]提出了一种基于状态空间拆分重组的离散化方法,将状态矩阵拆分成动态系数矩阵和常量系数矩阵,仅对动态系数矩阵进行离散化,该方法在2 kHz的离散化频率下具有较高的精度,但并没有在更低的离散化频率下验证其离散化精度。文献[22]以定子电流和转子磁链作为状态变量构建感应电机的状态空间方程,实现了对状态转移矩阵的精确离散化,但是缺少在极低载波比工况下开展定量的离散化误差分析。
本文以感应电机为例,研究了连续域下感应电机的传递函数和极点分布轨迹,推导了零阶保持离散化数学模型,分析了传统包括二阶双线性变换在内的降阶离散化模型在低载波比工况下的局限性,提出了特征根离散化模型,该模型在低载波比下仍具有良好的稳态精度和暂态跟随性能。同时,提出了基于伯德图的离散化误差定量分析方法,通过推导不同离散模型的磁链观测域传递函数,绘制其伯德图,以连续模型的伯德图为评判标准,对比了传统的离散化模型和提出的特征根离散化模型在磁链观测幅值和相位上的误差大小。最后,仿真和实验验证了所提出方法的有效性。
其中
式(1)和式(2)分别为描述定子电流和转子磁链变化模态的复数微分方程,可将该系统看作一个单输入、双输出,定转子间具有强耦合特性的二阶系统。
其中
绘制由式(3)和式(4)化简后的系统频率响应在空载和满载时的伯德图,能更直观地分析输入输出变量的稳态特性。感应电机连续域传递函数的伯德图如图1所示,电机系统是一个幅值衰减剧烈,相位滞后严重的系统。以连续系统的伯德图为基准能精确地定量对比不同离散化模型对系统描述精度的影响,后文将以此为手段进行离散化误差的相关分析。
(a)s/s伯德图
(b)r/s伯德图
图1 感应电机连续域传递函数的伯德图
Fig.1 The Bode plot of induction machine in continous domain
传统的离散化方法可以分为基于数值积分的降阶离散化方法(也叫近似法)和输入响应不变的离散化方法两类。基于数值积分的降阶离散化是将与的无理关系近似地化为有理关系,即
式中,s为离散化频率。常用的近似法包括一阶前向欧拉、后向欧拉、双线性变换法、二阶龙格库塔法等。其中,一阶前向欧拉模型具有实现简单、计算量小的优点,而双线性变换在以上降阶的离散化方法中具有最小的离散化误差,本文以此两种典型离散化模型为参考,在后文分析其在低载波比下的局限性。基于数值积分的离散化方法等效于对状态量在一拍范围内的变化进行线性近似,变化的斜率示意图如图2所示。状态变量在t-1~t之间的变化斜率在一阶前向欧拉模型中为1,该斜率只与t-1时刻的状态量和输入电压有关。对后向欧拉和双线性变化,斜率则分别为2和(1+2)/2。具体而言,一阶前向欧拉离散模型fir()和双线性变换离散模型secd()分别为
式中,()为连续域下定子电压到转子磁链的传递函数r/s,如式(4)所示。当控制系统处于低载波比工况时,由于忽略了高阶项的影响,基于数值积分的降阶离散化方法会出现明显的离散化误差,对应的状态观测将出现幅值和相位的稳态误差。
图2 降阶离散模型状态变量变化斜率示意图
输入响应不变的离散化方法具有离散化前后连续系统和离散系统对输入信号的响应不变这一特性,通常可分为变换(冲击响应不变离散化)、零阶保持离散化(阶跃响应不变离散化)和一阶保持离散化(脉冲响应不变离散化)三种。考虑到实际的电压指令由电压型逆变器基于一个开关周期内的伏秒平衡等效作用于电机系统;而零阶保持离散化方法假定输入量在一拍内维持恒定,和实际电压的作用方式一致,能保证在采样处连续系统和离散系统的状态变量尽可能一致。因此,本文也将零阶保持离散化方法作为参考,感应电机的零阶保持离散化模型zoh()为
式中,[ ]为变换。
感应电机通过气隙磁场切割转子,在转子内感应出电流,从而产生电磁转矩。转差频率在暂态时为时变量,因此严格意义上传递函数无法准确描述这一时变系统。
为更准确地描述系统的暂态特性,采取状态空间描述的方式对系统进行建模[11],有
其中
假定电气时间常数远小于机械时间常数,r为定值,则状态方程式(9)为线性定常系统。在静止坐标系下,由线性控制理论,状态变量()的精确解为
其中,在两相静止ab参考坐标系下有
其中
图3 感应电机在静止坐标系和转子旋转参考系下的极点分布轨迹
其中
为了获得完全离散化的感应电机数学模型,需对离散状态方程中的系统矩阵和输入矩阵精确求解,有
该方程本质是对系统矩阵的精确求解。系统矩阵的定义如式(11)所示,是关于转子速度的非线性矩阵函数,传统降阶近似计算不可避免引入近似误差,且随着转速升高而不断增加,而通过定义式完全求解该系统矩阵十分困难。
其中
将上述中间参数代入式(21),可得状态矩阵指数和输入矩阵的精确解分别为
其中
因此,结合式(15)~式(18)、式(22)和式(23)最终获得感应电机特征根离散化数学模型为
以Frobenius范数来表征矩阵数字特征,通过状态矩阵指数的近似程度进行定量分析[21],传统的离散化误差分析方法为
以图1为参考标准,通过传统离散模型的域传递函数式(6)和式(7),零阶保持离散模型的传递函数式(8)和特征根离散模型的域传递函数式(25)。在满载和空载时,分别绘出离散化频率为1 kHz时的伯德图,如图4和图5所示。
图4 1 kHz离散化频率下满载时不同离散化模型的伯德图对比
在图4和图5中,随着电频率的上升,不同离散化模型的误差在满载与空载的情况下,表现出相同的变化规律。随着转速的上升,传统一阶离散化模型与二阶离散化模型在幅值和相位曲线和连续系统相比均出现较大的偏差,且随着频率增加逐渐恶化,同样无法适用于低载波比的场合,而提出的特征根离散化模型仍然和连续系统模型高度吻合,离散化误差较小。
基于Matlab/Simulink对本文的研究内容进行了仿真验证,仿真和实验所用感应电机参数见表1,逆变器允许的最大开关频率为500 Hz,采用非对称采样模式,故采样频率为1 kHz,所有观测器及控制算法均在PWM中断中进行。
图5 1 kHz离散化频率下空载时不同离散化模型的伯德图对比
表1 2.2 kW感应电机参数
Tab.1 The parameters of 2.2 kW induction machine
图6为电机升速过程中一阶欧拉离散化模型的观测波形。图中,a为a相电流,Dqe为转子磁场角度误差,e为电频率,在所有仿真结果中均表示此含义。在中高速时电流观测出现发散不收敛趋势。随着转速的上升,观测的转子磁场角度和真实角度的偏差越来越大。
图7为不同离散模型观测的转子磁场角度和真实转子磁场角度的误差对比波形,在载波比为5时,二阶双线性离散化模型转子磁场角度和真实值偏差高达20°,而零阶保持离散化和特征根离散化模型观测得到的转子磁场角度稳态误差较小,磁场角度误差均在1°左右。
图8为在高速3 000 r/min,载波比低至5时的稳态观测电流对比。传统的二阶双线性离散化模型观测电流幅值和相位均出现了较大的稳态误差,而零阶保持离散化模型和特征根离散化模型的观测电流和真实电流基本一致。
图6 一阶欧拉离散化模型在升速时的观测波形
图7 不同离散化模型观测磁链角度的误差大小对比
图8 不同离散模型稳态电流波形对比(3 000 r/min)
图9 不同离散化模型在转矩阶跃时的电流波形对比(2 160 r/min)
图10 零阶保持和特征根离散模型在降速时的电流波形对比
为进一步验证本文的分析结果,采用dSPACE公司MicroLabBox控制器对文中的离散化数学模型进行了实验验证,小功率异步电机实验平台如图11所示,包括被试电机及其驱动器和控制器、陪试电机、高压电源、低压电源、调试计算机、示波器、变频器等。电机运行在基于特征根离散化模型设计的无差拍预测控制策略下,开关频率为500 Hz,采用非对称采样方式,计算和采样频率为1 kHz。离散化模型计算均由PWM中断同步触发运行,观测器输出的电流、转子磁链、磁场角度、电频率等变量通过D-A转换输出,实际电流通过录波仪实时记录。
图11 小功率异步电机实验平台
图12为升速实验中一阶欧拉离散化模型在中高速下的观测波形。在频率为20 Hz左右,一阶离散化模型观测的电流fir出现明显观测误差,且随着转速的上升,观测的电流开始出现发散趋势,和仿真(见图6)基本一致,而特征根离散化模型的观测电流eigen观测误差较小。
图12 一阶欧拉离散化模型在中高速下的观测波形
图13为高速下,载波比为7时,不同观测器的观测电流稳态对比波形。图中,a为电流探头采样的电流信号,同步采样电流a_syn为A-D在PWM中断中采样的真实A相电流信号,再通过D-A从3个通道输出作为其他观测器的参考波形。其中,二阶双线性模型观测电流secd在低载波比下能维持稳定,但观测的电流幅值和相位出现了较大误差,提出的零阶保持观测电流zoh和特征根离散化模型观测电流eigen均能很好地跟踪实际同步采样电流。
图13 不同离散化模型的稳态电流波形(2 160 r/min)
图14为2 160 r/min时,转矩从1 N·m升至3 N·m再降至1 N·m时的暂态电流波形对比,转矩阶跃时,二阶双线性离散化模型的观测电流出现了明显的相位误差,而零阶保持模型和特征根离散化模型的观测电流也出现了一定的观测误差,与仿真结果图9一致,但相比于二阶双线性离散化模型具有更好的动态跟随性能。
(a)转矩突增 (b)转矩突减
图14 转矩阶跃时的观测波形对比(2 160 r/min)
Fig.14 Comparison of observed waveforms during torque step at 2 160 r/min
图15为三种离散模型的转子磁场角度在稳态下的对比波形,由于真实的磁场角度未知,各个模型的磁场角度通过互相印证进行说明。在2 400 r/min时,特征根离散化模型和零阶保持器模型所对应的磁场角度具有较好的一致性,而二阶双线性模型和其他模型相比存在较大的偏差,约为10°,和仿真图7在80 Hz时二阶双线性模型对应的角度误差结果基本一致。
图15 2 400 r/min时的转子磁场角度对比
本文对低载波比下的感应电机离散模型和离散化误差分析方法展开了研究,主要取得了以下成果:
1)提出了感应电机特征根离散化模型,该模型在全速度范围内均能实现高精度的转子磁链和定子电流观测,在载波比低至5时,磁链角度观测误差小于1°,相比于传统二阶双线性模型高达20°的角度偏差,大大提升了模型的稳态精度,同时仿真和实验也验证了该模型具有良好的稳态精度和暂态跟随性能。
2)推导了不同离散模型的域传递函数,提出了基于伯德图的离散误差定量分析方法。该方法克服了传统基于矩阵误差的方法无法分析离散误差对状态量的幅值和相位影响的缺点,通过对比不同离散模型和连续模型的伯德图,从理论上验证了传统方法的局限性和提出方法的优越性。
[1] Odhano S A, Pescetto P, Ali Awan H A, et al. Parameter identification and self-commissioning in AC motor drives: a technology status review[J]. IEEE Transactions on Power Electronics, 2019, 34(4): 3603-3614.
[2] Zhao Lihang, Huang Jin, Chen Jiahao, et al. A parallel speed and rotor time constant identification scheme for indirect field oriented induction motor drives[J]. IEEE Transactions on Power Electronics, 2016, 31(9): 6494-6503.
[3] 李婕, 杨淑英, 谢震, 等. 基于有效信息迭代快速粒子群优化算法的永磁同步电机参数在线辨识[J]. 电工技术学报, 2022, 37(18): 4604-4613.
Li Jie, Yang Shuying, Xie Zhen, et al. Online para- meter identification of permanent magnet syn- chronous motor based on fast particle swarm optimization algorithm with effective information iterated[J]. Transactions of China Electrotechnical Society, 2022, 37(18): 4604-4613.
[4] 徐伟, 董定昊, 葛健, 等. 基于在线参数辨识补偿的直线感应电机低开关频率模型预测控制策略[J]. 电工技术学报, 2022, 37(16): 4116-4133.
Xu Wei, Dong Dinghao, Ge Jian, et al. Low switching frequency model predictive control strategy based on online parameter identification compensation of linear induction motor for urban rail application[J]. Transa- ctions of China Electrotechnical Society, 2022, 37(16): 4116-4133.
[5] 黄科元, 周佳新, 刘思美, 等. 考虑逆变器非线性永磁同步电机高频注入电感辨识方法[J]. 电工技术学报, 2021, 36(8): 1607-1616.
Huang Keyuan, Zhou Jiaxin, Liu Simei, et al. Indu- ctance identification method of permanent magnet synchronous motor considering inverter nonlinearity based on high-frequency injection[J]. Transactions of China Electrotechnical Society, 2021, 36(8): 1607-1616.
[6] 刘亚静, 段超. 全数字自适应滤波器不同离散结构的性能对比分析[J]. 电工技术学报, 2021, 36(20): 4339-4349.
Liu Yajing, Duan Chao. Performance comparison and analysis of all-digital adaptive filter with different discrete methods[J]. Transactions of China Electro- technical Society, 2021, 36(20): 4339-4349.
[7] Xu Yang, Morito C, Lorenz R D. Accurate discrete- time modeling for improved torque control accuracy for induction machine drives at very low sampling- to-fundamental frequency ratios[J]. IEEE Transa- ctions on Transportation Electrification, 2020, 6(2): 668-678.
[8] 孙建业, 王志强, 谷鑫, 等. 高速低载波比下永磁同步电机预测电流控制[J]. 中国电机工程学报, 2020, 40(11): 3663-3672.
Sun Jianye, Wang Zhiqiang, Gu Xin, et al. Predictive current control of permanent magnet synchronous motor under high speed and low carrier ratio[J]. Proceedings of the CSEE, 2020, 40(11): 3663-3672.
[9] 李杰, 詹榕, 宋文祥. 感应电机低采样频率的磁链观测器离散化模型研究[J]. 电工技术学报, 2019, 34(15): 3136-3146.
Li Jie, Zhan Rong, Song Wenxiang. Improved discrete observer model of induction motor at low sampling frequency[J]. Transactions of China Electrotechnical Society, 2019, 34(15): 3136-3146.
[10] 罗慧. 感应电机全阶磁链观测器和转速估算方法研究[D]. 武汉: 华中科技大学, 2009.
[11] Diao Lijun, Sun Danan, Dong Kan, et al. Optimized design of discrete traction induction motor model at low-switching frequency[J]. IEEE Transactions on Power Electronics, 2013, 28(10): 4803-4810.
[12] Wang Bo, Zhao Yongzheng, Yu Yong, et al. Speed- sensorless induction machine control in the field- weakening region using discrete speed-adaptive full- order observer[J]. IEEE Transactions on Power Elec- tronics, 2016, 31(8): 5759-5773.
[13] 李昱, 郭宏, 平朝春, 等. 基于电流源变流器的永磁同步电机驱动系统全状态变量预测转矩控制[J]. 电工技术学报, 2021, 36(1): 15-26.
Li Yu, Guo Hong, Ping Zhaochun, et al. A full-state variable predictive torque control of current source converter fed permanent magnet synchronous motor drives[J]. Transactions of China Electrotechnical Society, 2021, 36(1): 15-26.
[14] Holtz J, Quan Juntao, Pontt J, et al. Design of fast and robust current regulators for high-power drives based on complex state variables[J]. IEEE Transactions on Industry Applications, 2004, 40(5): 1388-1397.
[15] 国敬, 范涛, 章回炫, 等. 高速低载波比下永磁同步电机电流环稳定性分析[J]. 中国电机工程学报, 2019, 39(24): 7336-7346, 7506.
Guo Jing, Fan Tao, Zhang Huixuan, et al. Stability analysis of permanent magnet synchronous motor current loop control at high speed and low carrier ratio[J]. Proceedings of the CSEE, 2019, 39(24): 7336-7346, 7506.
[16] Kim H, Degner M, Guerrero J M, et al. Discrete-time current regulator design for AC machine drives[C]// 2009 IEEE Energy Conversion Congress and Expo- sition, San Jose, CA, USA, 2009: 1317-1324.
[17] Jarzebowicz L. Errors of a linear current approxi- mation in high-speed PMSM drives[J]. IEEE Transa- ctions on Power Electronics, 2017, 32(11): 8254- 8257.
[18] Jarzebowicz L. Quasi-discrete modelling of PMSM phase currents in drives with low switching-to- fundamental frequency ratio[J]. IET Power Electro- nics, 2019, 12(12): 3280-3285.
[19] Dai Shangjian, Wang Jiabin, Sun Zhigang, et al. Deadbeat predictive current control for high-speed permanent magnet synchronous machine drives with low switching-to-fundamental frequency ratios[J]. IEEE Transactions on Industrial Electronics, 2022, 69(5): 4510-4521.
[20] West N T, Lorenz R D. Digital implementation of stator and rotor flux-linkage observers and a stator- current observer for deadbeat direct torque control of induction machines[J]. IEEE Transactions on Industry Applications, 2009, 45(2): 729-736.
[21] 赵雷廷, 刁利军, 董侃, 等. 基于状态空间拆分重组的牵引异步电机闭环离散全阶转子磁链观测器[J]. 电工技术学报, 2013, 28(10): 103-112.
Zhao Leiting, Diao Lijun, Dong Kan, et al. A novel discretized closed-loop full-order rotor flux observer for induction motor based on re-organization of state space[J]. Transactions of China Electrotechnical Society, 2013, 28(10): 103-112.
[22] Jean C A. Direct eigen control for induction machines and synchronous motors[M]. United Kingdom: John Wiley&Sons, 2013.
[23] 尾形克彦 (Katsuhiko Ogata). 离散时间控制系统[M]. 陈杰, 蔡涛, 张娟, 等译. 北京: 机械工业出版社, 2006.
Research on Discretization Model of Induction Motor for Low Switching-to-Fundamental Frequency Ratio Traction System
111213
(1. State Key Laboratory of Advanced Electromagnetic Engineering and Technology School of Electrical and Electronic Engineering Huazhong University of Science and Technology Wuhan 430074 China 2. Engineering Research Center of Novel Electrical Machines and Special Electromagnetic Equipment Ministry of Education School of Electrical and Electronic Engineering Huazhong University of Science and Technology Wuhan 430074 China 3. CRRC Dalian R&D Co. Ltd Dalian 116052 China)
In high-power and high-speed motor drives, the control system will operate in low switching- to-fundamental frequency ratio conditions. Due to the large discretization error, the traditional reduced-order discrete model cannot be applied to the control system. Therefore, this paper proposes an eigenvalue-based discrete induction motor model, which still has good steady-state accuracy and transient tracking performance under a low switching-to-fundamental frequency ratio. At the same time, a quantitative analysis method of discretization error based on the Bode diagram is proposed.
Firstly, the mathematical model of the induction motor is modeled in the continuous domain using state space description. The full-rank state matrix is diagonalized through the transformation matrix. The elements on the diagonal of the diagonal matrix are the characteristic roots of the corresponding stator and rotor voltage equations of the induction motor. Then, the exact solution of the system matrix can be obtained through the transformation matrix and the diagonal matrix, and the eigenvalue-based discrete model of the induction motor is derived. Moreover, by deducing the z-domain transfer function of flux linkage observation of different discrete models, the Bode diagram is drawn. Taking the Bode diagram of the continuous model as the evaluation standard, the traditional and the proposed discrete models’ errors in the magnitude and phase of flux linkage observation is compared.
The simulation and experimental results show that in terms of the steady-state observation accuracy, at the electric frequency of about 20 Hz, the observed current of the first-order discretization model has an obvious observation error. In contrast, the observed current of the eigenvalue-based discrete model has a small observation error. When the switching-to-fundamental frequency ratio is seven, the observed current of the second-order bilinear model can maintain stability. However, the observed current amplitude and phase have large errors, while the observed current of the characteristic root discretization model can track the actual synchronous sampling current. Regarding dynamic observation accuracy, in speed reduction, the zero-order holder discrete model has a specific amplitude and phase deviation in current observation. The maximum amplitude deviation is up to 30%, and the deviation disappears after the speed enters the steady state. The eigenvalue-based discrete model has good performance in both transient and steady speeds.
The following conclusions can be drawn: (1) The proposed discrete model of the characteristic root of the induction motor can achieve high-precision rotor flux and stator current observation in the full-speed range. When the switching-to-fundamental frequency ratio is as low as five, the observation error of the flux angle is less than 1°. (2) Compared with the zero-order holder discretization model, the proposed eigenvalue-based discrete model has higher current observation accuracy in speed reduction.
Induction machine, discretization method, transfer function, Bode plot, low switching-to- fundamental frequency ratio
10.19595/j.cnki.1000-6753.tces.222046
TM301.2
2022-10-28
2023-01-19
张钦培 男,1998年生,硕士研究生,主要研究方向为感应电机数学模型与控制策略。E-mail: zqp@hust.edu.cn
李 健 男,1982年生,研究员,博士生导师,主要研究方向为大功率牵引变流器控制方面的理论和技术开发。E-mail: jianli@hust.edu.cn(通信作者)
(编辑 崔文静)