CN107528514A - Pmsm调速系统的近似离散快速终端滑模控制方法 - Google Patents
Pmsm调速系统的近似离散快速终端滑模控制方法 Download PDFInfo
- Publication number
- CN107528514A CN107528514A CN201710981166.XA CN201710981166A CN107528514A CN 107528514 A CN107528514 A CN 107528514A CN 201710981166 A CN201710981166 A CN 201710981166A CN 107528514 A CN107528514 A CN 107528514A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mtd
- msubsup
- mfrac
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000004364 calculation method Methods 0.000 claims abstract description 6
- 238000000034 method Methods 0.000 claims description 12
- 230000001105 regulatory effect Effects 0.000 claims description 8
- 230000001360 synchronised effect Effects 0.000 claims description 6
- 238000006243 chemical reaction Methods 0.000 claims description 5
- 230000001276 controlling effect Effects 0.000 claims description 5
- 230000004907 flux Effects 0.000 claims description 4
- 230000008859 change Effects 0.000 abstract description 6
- 230000001133 acceleration Effects 0.000 abstract description 3
- 230000010354 integration Effects 0.000 abstract description 2
- 230000003044 adaptive effect Effects 0.000 abstract 1
- 238000000605 extraction Methods 0.000 abstract 1
- 238000013459 approach Methods 0.000 description 4
- 238000013461 design Methods 0.000 description 3
- 230000007613 environmental effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000011217 control strategy Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000005312 nonlinear dynamic Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P21/00—Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
- H02P21/0003—Control strategies in general, e.g. linear type, e.g. P, PI, PID, using robust control
- H02P21/0007—Control strategies in general, e.g. linear type, e.g. P, PI, PID, using robust control using sliding mode control
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P21/00—Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
- H02P21/14—Estimation or adaptation of machine parameters, e.g. flux, current or voltage
- H02P21/18—Estimation of position or speed
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P21/00—Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
- H02P21/22—Current control, e.g. using a current control loop
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Control Of Ac Motors In General (AREA)
- Pharmaceuticals Containing Other Organic And Inorganic Compounds (AREA)
Abstract
本发明提供了一种PMSM调速系统的近似离散快速终端滑模控制方法，首先接收角速率主令，提取电机角速度反馈，计算速度误差的一阶导数；然后计算速度环的快速终端滑模面和调速系统的电流环指令；再利用角加速度反馈计算电流环控制量，采集线电流计算PMSM在dq坐标下的交直轴电流及其误差的一阶变化量，计算dq轴电流环的快速终端滑模面；最后计算的电流环dq轴电压控制量作为PMSM的逆park变换的输入来完成电机控制。本发明提升了随动调速系统的控制鲁棒性，对电机参数的变化和外部力矩干扰不敏感，增强了系统的环境适应能力；速度环和电流环控制项都采用积分对滑模变结构切换相进行平滑处理，避免了传统的线性滑模的控制抖振。
Description
技术领域
本发明属于火炮随动系统领域，主要涉及需要精确跟踪的随动驱动调速系统对参数变化和负载扰动不敏感的鲁棒控制方法。
背景技术
火炮随动系统尤其是高炮随动系统要求随动驱动调速系统具有快速响应指令的能力和对指令精确跟踪的能力。由于PI控制器控制策略结构简单，易于调节，传统的驱动调速系统的电流环和速度环控制器常采用PI控制器。电流环常采用时域的最佳模原则、频域的幅值裕度和相位裕度要求进行设计PI控制参数，速度环也常采用频域的理想中频宽原则设计合适的PI控制器参数。但是驱动调速系统的被控对象或者是随动系统的执行机构交流永磁同步电机(PMSM)是强耦合的非线性系统，且参数会随环境温度变化而变化，常常有难以预测的外部扰动和复杂的未建模状态，仅仅采用PI控制难以获得高性能。为了弥补PI经典控制理论对非线性系统调节能力不足的问题，国内外学者将模型参考自适应(MARS)、反步控制(back_step)、二自由度控制及滑模变结构控制等用于电机控制，以期望获得良好的静态和动态性能。但是有的算法在离散状态的研究和应用并不成熟。
滑模变结构控制由于对内部参数摄动和外部干扰具有较强的鲁棒性和较高的控制精度，且实现简单，故成为提高永磁同步电机控制系统的有效手段之一。传统滑模变结构控制常采用是线性滑模s＝CX，其中C为滑模系数，X为状态，如角速度、速度、加速度等。该类滑模控制系统稳定性分析简单、便捷，参数设计相对容易，滑模系数可以采用最优化控制法、极点配置法等进行设计，但线性滑模在非线性动态系统中的的收敛速度不够，无法使系统在有限时间内收敛到平衡点，对抑制不确定的鲁棒性难以发挥，故在复杂非线性系统中的应用则控制能力则明显不足，限制了在随动驱动调速系统中的应用。
发明内容
为了克服现有技术的不足，本发明提供一种PMSM调速系统的近似离散快速终端滑模控制方法，能够提升随动调速系统的控制性能和鲁棒性。
本发明解决其技术问题所采用的技术方案包括以下步骤：
(1)等待速度指令ω*(k)，速度指令到达后进入步骤(2)，k为控制的步数；
(2)提取电机角速度反馈ω(k)；
(3)计算速度误差的一阶导数deω(k)：
eω(k-1)＝ω*(k-1)-ω(k-1)
eωz(k-1)＝zω0(k-1)-eω(k-1)
zω0(k)＝zω0(k-1)+Ts(zω1(k-1)-λω1|eωz(k-1)|1/2sgn(eωz(k-1)))
zω1(k)＝zω1(k-1)+Tsλω2sgn(-λω1|eωz(k-1)|1/2sgn(eωz(k-1)))
deω(k)＝zω1(k)
其中，Ts为计算周期；eω(k-1)为速度误差；eωz(k-1)为速度误差的估计误差；zω0(k)为速度误差的估计值；zω1(k)为速度误差一阶导数deω(k)的估计值；λω1与λω2分别为微分估计器一阶、二阶增益；sgn(·)为符号函数；
(4)计算速度环的快速终端滑模面其中，αω与βω是快速终端滑模面sω(k)的系数，且αω＞0，βω＞0；qω与pω是奇数，pω＞qω＞0；
(5)计算调速系统的电流环指令
其中，和分别为电流指令的最大值和最小值；iqsw(k)为滑模积分控制分量；iqeq(k)为等效控制分量；J为电机轴与负载等效至电机轴的总转动惯量；np为电机的极对数；ψf为定子永磁磁链；ηω与kω分别为指数趋近率的系数，且ηω＞0，kω＞0；B为粘性摩擦系数；
(6)采集线电流ia(k),ib(k)，计算PMSM在dq坐标下的交直轴电流id(k),iq(k)：
(7)计算PMSM在dq坐标下的交直轴电流误差的一阶变化量ded(k),deq(k)
ed(k-1)＝-id(k-1)
edz(k-1)＝zd0(k-1)-ed(k-1)
zd0(k)＝zd0(k-1)+Ts(zd1(k-1)-λd1|edz(k-1)|1/2sgn(edz(k-1)))
zd1(k)＝zd1(k-1)+Tsλd2sgn(-λd1|edz(k-1)|1/2sgn(edz(k-1)))
ded(k)＝zd1(k)
eqz(k-1)＝zq0(k-1)-eq(k-1)
zq0(k)＝zq0(k-1)+Ts(zq1(k-1)-λq1|eqz(k-1)|1/2sgn(eqz(k-1)))
zq1(k)＝zq1(k-1)+Tsλq2sgn(-λq1|eqz(k-1)|1/2sgn(eqz(k-1)))
deq(k)＝zq1(k)
其中，ed(k-1)和eq(k-1)分别为dq轴电流误差；edz(k-1)和eqz(k-1)分别为dq轴电流误差的估计误差；zd0(k)和zq0(k)分别为dq轴电流误差的估计值；zd1(k)和zq1(k)分别为dq轴电流误差一阶导数ded(k)、deq(k)的估计值；λd1与λd2分别为d轴微分估计器一阶、二阶增益；λq1与λq2分别为d轴微分估计器一阶、二阶增益；
(8)计算dq轴电流环的快速终端滑模面sd(k),sq(k)
其中，αd与βd是d轴电流环快速终端滑模面sd(k)的系数，且αd＞0，βd＞0；qd与pd是奇数，pd＞qd＞0；αq与βq是q轴电流环快速终端滑模面sq(k)的系数，且αq＞0，βq＞0；qq与pq是奇数，pd＞qd＞0；
(9)计算的电流环d轴电压控制量ud(k)
ud(k)＝udeq(k)+udsw(k)
其中，umax和umin分别电压控制量ud(k)、uq(k)的最大值和最小值；L为PMSM电机dq轴交直轴电感；Rs为电机定子电阻；udsw(k)为滑模积分控制分量；udeq(k)为等效控制分量；ηd与kd分别为指数趋近率的系数，且ηd＞0，kd＞0；
(10)计算电流环q轴电压控制量uq(k)
uq(k)＝uqeq(k)+uqsw(k)
其中，uqsw(k)为滑模积分控制分量；uqeq(k)为等效控制分量；ηq与kq分别为指数趋近率的系数，且ηω＞0，kω＞0；
(11)将dq轴电压控制量ud(k),uq(k)作为PMSM的逆park变换的输入来完成电机的k步控制。
本发明的有益效果是：在电流环实现电流前馈的解耦控制，并对反电动势进行补偿，提升电流环的动态控制精度；提升了随动调速系统的控制鲁棒性，对电机参数的变化和外部力矩干扰不敏感，增强了系统的环境适应能力；速度环和电流环控制项都采用积分对滑模变结构切换相进行平滑处理，避免了传统的线性滑模的控制抖振。
附图说明
图1是本发明的控制原理图；
图2本发明的计算流程图。
具体实施方式
下面结合附图和实施例对本发明进一步说明，本发明包括但不仅限于下述实施例。
本发明能够解决火炮随动系统传统PI控制器的性能对执行机构PMSM参数(如定子永磁磁链ψf、定子电阻Rs，交直轴电感Ld,Lq随环境变化而变化，以及外部力矩干扰等鲁棒性不强，实现高性能的高控制效果具有很大难度的问题。针对火炮随动系统常采用的隐极电机(交直轴电感相等，Ld＝Lq＝L)，采用快速终端滑模变结构控制方法来设计速度环和电流环控制器，并采用一阶欧拉离散化来近似实现，以提升随动调速系统的控制性能和鲁棒性。其方案有如下步骤：
(1)判断速度指令ω*(k)是否到？是，进入步骤(2)，否则就等待，k为控制的步数；
(2)提取电机角速度反馈ω(k)；
(3)计算速度误差的一阶导数deω(k)：
eω(k-1)＝ω*(k-1)-ω(k-1)
eωz(k-1)＝zω0(k-1)-eω(k-1)
zω0(k)＝zω0(k-1)+Ts(zω1(k-1)-λω1|eωz(k-1)|1/2sgn(eωz(k-1)))
zω1(k)＝zω1(k-1)+Tsλω2sgn(-λω1|eωz(k-1)|1/2sgn(eωz(k-1)))
deω(k)＝zω1(k)
其中：Ts为计算周期；eω(k-1)为速度误差；eωz(k-1)为速度误差的估计误差；zω0(k)为速度误差的估计值；zω1(k)为速度误差一阶导数deω(k)的估计值；λω1与λω2分别为微分估计器一阶、二阶增益；sgn(·)为符号函数；
(4)计算速度环的快速终端滑模面sω(k)：
其中：αω与βω是快速终端滑模面sω(k)的系数，且αω＞0，βω＞0；qω与pω是奇数，pω＞qω＞0；
(5)计算调速系统的电流环指令
其中：和分别电流指令的最大值和最小值；iqsw(k)为滑模积分控制分量；iqeq(k)为等效控制分量；J为电机轴与负载等效至电机轴的总转动惯量；np为电机的极对数；ψf为定子永磁磁链；ηω与kω分别为指数趋近率的系数，且ηω＞0，kω＞0；B为粘性摩擦系数；
(6)采集线电流ia(k),ib(k)计算PMSM在dq坐标下的交直轴电流id(k),iq(k)
(7)计算PMSM在dq坐标下的交直轴电流误差的一阶变化量ded(k),deq(k)
ed(k-1)＝-id(k-1)
edz(k-1)＝zd0(k-1)-ed(k-1)
zd0(k)＝zd0(k-1)+Ts(zd1(k-1)-λd1|edz(k-1)|1/2sgn(edz(k-1)))
zd1(k)＝zd1(k-1)+Tsλd2sgn(-λd1|edz(k-1)|1/2sgn(edz(k-1)))
ded(k)＝zd1(k)
eqz(k-1)＝zq0(k-1)-eq(k-1)
zq0(k)＝zq0(k-1)+Ts(zq1(k-1)-λq1|eqz(k-1)|1/2sgn(eqz(k-1)))
zq1(k)＝zq1(k-1)+Tsλq2sgn(-λq1|eqz(k-1)|1/2sgn(eqz(k-1)))
deq(k)＝zq1(k)
其中：ed(k-1)和eq(k-1)分别为dq轴电流误差；edz(k-1)和eqz(k-1)分别为dq轴电流误差的估计误差；zd0(k)和zq0(k)分别为dq轴电流误差的估计值；zd1(k)和zq1(k)分别为dq轴电流误差一阶导数ded(k)、deq(k)的估计值；λd1与λd2分别为d轴微分估计器一阶、二阶增益；λq1与λq2分别为d轴微分估计器一阶、二阶增益；
(8)计算dq轴电流环的快速终端滑模面sd(k),sq(k)
其中：αd与βd是d轴电流环快速终端滑模面sd(k)的系数，且αd＞0，βd＞0；qd与pd是奇数，pd＞qd＞0；αq与βq是q轴电流环快速终端滑模面sq(k)的系数，且αq＞0，βq＞0；qq与pq是奇数，pd＞qd＞0；
(9)计算的电流环d轴电压控制量ud(k)
ud(k)＝udeq(k)+udsw(k)
其中：umax和umin分别电压控制量ud(k)、uq(k)的最大值和最小值；L为PMSM电机dq轴交直轴电感；Rs为电机定子电阻；udsw(k)为滑模积分控制分量；udeq(k)为等效控制分量；ηd与kd分别为指数趋近率的系数，且ηd＞0，kd＞0；
(10)计算电流环q轴电压控制量uq(k)
uq(k)＝uqeq(k)+uqsw(k)
其中：uqsw(k)为滑模积分控制分量；uqeq(k)为等效控制分量；ηq与kq分别为指数趋近率的系数，且ηω＞0，kω＞0；
(11)将dq轴电压控制量ud(k),uq(k)作为PMSM的逆park变换的输入来完成电机的k步控制。
本发明实施例的控制原理见图1。图中在一般的火炮随动隐极PMSM驱动驱动原理图，将传统的PI控制器采用快速终端滑模变结构控制器取代。该控制方法大致步骤：首先，接收角速率主令ω*(k)，提取电机角速度反馈ω(k)，计算速度误差的一阶导数deω(k)；然后，计算速度环的快速终端滑模面sω(k)，计算调速系统的电流环指令其次，利用角加速度反馈计算电流环控制量再次，采集线电流ia(k),ib(k)计算PMSM在dq坐标下的交直轴电流id(k),iq(k)，计算PMSM在dq坐标下的交直轴电流误差的一阶变化量ded(k),deq(k)；再再次，计算dq轴电流环的快速终端滑模面sd(k),sq(k)；最后，计算的电流环dq轴电压控制量ud(k),uq(k)，将dq轴电压控制量作为PMSM的逆park变换的输入来完成电机的k步控制。
根据隐极PMSM误差连续状态方程：
其中：d轴电流误差q轴电流误差角速率误差 分别为电流环dq轴电流主令；id,iq分别为电流环dq轴电流反馈；ω*,ω分别为角速度主令和角速度反馈；εd(Δd),εq(Δq)分别为dq轴状态方程为建模动态或扰动分量，且有上界，εω(TL,Δω)为速度误差方程含负载转矩和为建模误差，且有上界。分别取如下非线性快速终端滑模面：
分别对以上滑模面求导，分别采用指数趋近率其中η,k为指数趋近率系数，可得到电流环指令电流环dq轴电压控制ud,uq，可以验证分别满足李雅普诺夫稳定性要求，即
其中Vd,Vq,Vω分别为dq和速率滑模面能量函数。由于控制周期Ts为0.1ms，时间非常短，利用一阶欧拉离散化近似实现连续状态快速终端滑模控制是可行的。
实施该控制方法的火炮随动驱动调速系统，该计算平台是以DSP28335为核心的控制板。功率驱动采用IPM驱动，选用PMSM电机，母线电压为56VDC，极对数np＝3，转矩系数ψf＝0.368Wb，定子电阻Rs＝0.008Ω，交直轴电感L＝0.2mH，功率3Kw，可过载2倍，额定转速3500RPM，电机转子等效转动惯量J为0.0006329kg·m2。
图2为本发明控制方法的计算流程图，下面将结合流程图详述具体实施过程：
(1)判断速度指令ω*(k)是否到？是，进入步骤(2)，否则就等待，k为控制的步数；
(2)提取电机角速度反馈ω(k)；
(3)计算速度误差的一阶导数deω(k)：
eω(k-1)＝ω*(k-1)-ω(k-1)
eωz(k-1)＝zω0(k-1)-eω(k-1)
zω0(k)＝zω0(k-1)+Ts(zω1(k-1)-λω1|eωz(k-1)|2/3sgn(eωz(k-1)))
zω1(k)＝zω1(k-1)+Tsλω2sgn(-λω1|eωz(k-1)|1/2sgn(eωz(k-1)))
deω(k)＝zω1(k)
其中：计算周期Ts＝1ms；速度误差的估计值初值zω0(0)＝0；速度误差一阶导数deω(k)的估计值初值zω1(0)＝0；微分估计器一阶增益λω1＝5.84、二阶增益λω2＝15；
(4)计算速度环的快速终端滑模面sω(k)：
其中：是快速终端滑模面sω(k)的系数αω＝300，βω＝90；奇数qω＝3，奇数pω＝5；
(5)计算调速系统的电流环指令
其中：电流指令的最大值最小值电机轴与负载等效至电机轴的总转动惯量J为0.0006329；电机的极对数np＝3；定子永磁磁链ψf＝0.368；指数趋近率的系数分别为ηω＝2500，kω＝10000；B＝0.05487；
(6)采集线电流ia(k),ib(k)计算PMSM在dq坐标下的交直轴电流id(k),iq(k)
(7)计算PMSM在dq坐标下的交直轴电流误差的一阶变化量ded(k),deq(k)
ed(k-1)＝-id(k-1)
edz(k-1)＝zd0(k-1)-ed(k-1)
zd0(k)＝zd0(k-1)+Ts(zd1(k-1)-λd1|edz(k-1)|1/2sgn(edz(k-1)))
zd1(k)＝zd1(k-1)+Tsλd2sgn(-λd1|edz(k-1)|1/2sgn(edz(k-1)))
ded(k)＝zd1(k)
eqz(k-1)＝zq0(k-1)-eq(k-1)
zq0(k)＝zq0(k-1)+Ts(zq1(k-1)-λq1|eqz(k-1)|1/2sgn(eqz(k-1)))
zq1(k)＝zq1(k-1)+Tsλq2sgn(-λq1|eqz(k-1)|1/2sgn(eqz(k-1)))
deq(k)＝zq1(k)
其中：dq轴电流误差的估计值初值分别为zd0(k)＝0，zq0(k)＝0；dq轴电流误差一阶导数ded(k)、deq(k)的估计值初值zd1(0)＝0，zq1(k)＝0；d轴微分估计器一阶增益λd1＝6.52、二阶增益λd2＝24；分别为d轴微分估计器一阶增益λq1＝6.52，二阶增益λq2＝24；
(8)计算dq轴电流环的快速终端滑模面sd(k),sq(k)
其中：d轴电流环快速终端滑模面sd(k)的系数分别为αd＝450，βd＝150；奇数qd＝3，奇数pd＝5；q轴电流环快速终端滑模面sq(k)的系数分别为αq＝450，βq＝150；奇数qq＝3，奇数pq＝5；
(9)计算的电流环d轴电压控制量ud(k)
ud(k)＝udeq(k)+udsw(k)
其中：电压控制量的最大值umax＝32，最小值umin＝-32；PMSM电机dq轴交直轴电感L＝0.2mH；电机定子电阻Rs＝0.008；指数趋近率的系数分别为ηd＝9000，kd＝80000；
(10)计算电流环q轴电压控制量uq(k)
uq(k)＝uqeq(k)+uqsw(k)
其中：指数趋近率的系数分别为ηq＝9000，kq＝80000；
(11)将dq轴电压控制量ud(k),uq(k)作为PMSM的逆park变换的输入来完成电机的k步控制。
本发明所需要参数见表1。
表1本发明所需要参数取值表
Claims (1)
1.一种PMSM调速系统的近似离散快速终端滑模控制方法，其特征在于包括下述步骤：
(1)等待速度指令ω*(k)，速度指令到达后进入步骤(2)，k为控制的步数；
(2)提取电机角速度反馈ω(k)；
(3)计算速度误差的一阶导数deω(k)：
eω(k-1)＝ω*(k-1)-ω(k-1)
eωz(k-1)＝zω0(k-1)-eω(k-1)
zω0(k)＝zω0(k-1)+Ts(zω1(k-1)-λω1|eωz(k-1)|1/2sgn(eωz(k-1)))
zω1(k)＝zω1(k-1)+Tsλω2sgn(-λω1|eωz(k-1)|1/2sgn(eωz(k-1)))
deω(k)＝zω1(k)
其中，Ts为计算周期；eω(k-1)为速度误差；eωz(k-1)为速度误差的估计误差；zω0(k)为速度误差的估计值；zω1(k)为速度误差一阶导数deω(k)的估计值；λω1与λω2分别为微分估计器一阶、二阶增益；sgn(·)为符号函数；
(4)计算速度环的快速终端滑模面其中，αω与βω是快速终端滑模面sω(k)的系数，且αω＞0，βω＞0；qω与pω是奇数，pω＞qω＞0；
(5)计算调速系统的电流环指令
<mrow>
<msub>
<mi>i</mi>
<mrow>
<mi>q</mi>
<mi>s</mi>
<mi>w</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mn>1.5</mn>
<msub>
<mi>n</mi>
<mi>p</mi>
</msub>
<msub>
<mi>&psi;</mi>
<mi>f</mi>
</msub>
</mrow>
<mi>J</mi>
</mfrac>
<msub>
<mi>T</mi>
<mi>s</mi>
</msub>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>k</mi>
</munderover>
<mrow>
<mo>(</mo>
<msub>
<mi>&alpha;</mi>
<mi>&omega;</mi>
</msub>
<mi>d</mi>
<mi>e</mi>
<mo>(</mo>
<mi>j</mi>
<mo>)</mo>
<mo>+</mo>
<msub>
<mi>&beta;</mi>
<mi>&omega;</mi>
</msub>
<mfrac>
<msub>
<mi>q</mi>
<mi>&omega;</mi>
</msub>
<msub>
<mi>p</mi>
<mi>&omega;</mi>
</msub>
</mfrac>
<mi>e</mi>
<msup>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>q</mi>
<mi>&omega;</mi>
</msub>
<mo>-</mo>
<msub>
<mi>p</mi>
<mi>&omega;</mi>
</msub>
<mo>)</mo>
<mo>/</mo>
<msub>
<mi>p</mi>
<mi>&omega;</mi>
</msub>
</mrow>
</msup>
<msub>
<mi>de</mi>
<mi>&omega;</mi>
</msub>
<mo>(</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>&eta;</mi>
<mi>&omega;</mi>
</msub>
<mi>s</mi>
<mi>g</mi>
<mi>n</mi>
<mo>(</mo>
<mrow>
<msub>
<mi>s</mi>
<mi>&omega;</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
<mo>+</mo>
<msub>
<mi>k</mi>
<mi>&omega;</mi>
</msub>
<msub>
<mi>s</mi>
<mi>&omega;</mi>
</msub>
<mo>(</mo>
<mi>j</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mrow>
<msub>
<mi>i</mi>
<mrow>
<mi>q</mi>
<mi>e</mi>
<mi>q</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mi>J</mi>
<mrow>
<mn>1.5</mn>
<msub>
<mi>n</mi>
<mi>p</mi>
</msub>
<msub>
<mi>&psi;</mi>
<mi>f</mi>
</msub>
</mrow>
</mfrac>
<msup>
<mi>&omega;</mi>
<mo>*</mo>
</msup>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mfrac>
<mi>B</mi>
<mrow>
<mn>1.5</mn>
<msub>
<mi>n</mi>
<mi>p</mi>
</msub>
<msub>
<mi>&psi;</mi>
<mi>f</mi>
</msub>
</mrow>
</mfrac>
<mi>&omega;</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msubsup>
<mi>i</mi>
<mi>q</mi>
<mo>*</mo>
</msubsup>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>i</mi>
<mrow>
<mi>q</mi>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
<mo>*</mo>
</msubsup>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<msubsup>
<mi>i</mi>
<mi>q</mi>
<mo>*</mo>
</msubsup>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>&GreaterEqual;</mo>
<msubsup>
<mi>i</mi>
<mrow>
<mi>q</mi>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
<mo>*</mo>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>i</mi>
<mrow>
<mi>q</mi>
<mi>e</mi>
<mi>q</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>i</mi>
<mrow>
<mi>q</mi>
<mi>s</mi>
<mi>w</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<msubsup>
<mi>i</mi>
<mrow>
<mi>q</mi>
<mi>min</mi>
</mrow>
<mo>*</mo>
</msubsup>
<mo><</mo>
<msubsup>
<mi>i</mi>
<mi>q</mi>
<mo>*</mo>
</msubsup>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo><</mo>
<msubsup>
<mi>i</mi>
<mrow>
<mi>q</mi>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
<mo>*</mo>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>i</mi>
<mrow>
<mi>q</mi>
<mi>min</mi>
</mrow>
<mo>*</mo>
</msubsup>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<msubsup>
<mi>i</mi>
<mi>q</mi>
<mo>*</mo>
</msubsup>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>&le;</mo>
<msubsup>
<mi>i</mi>
<mrow>
<mi>q</mi>
<mi>min</mi>
</mrow>
<mo>*</mo>
</msubsup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中，和分别为电流指令的最大值和最小值；iqsw(k)为滑模积分控制分量；iqeq(k)为等效控制分量；J为电机轴与负载等效至电机轴的总转动惯量；np为电机的极对数；ψf为定子永磁磁链；ηω与kω分别为指数趋近率的系数，且ηω＞0，kω＞0；B为粘性摩擦系数；
(6)采集线电流ia(k),ib(k)，计算PMSM在dq坐标下的交直轴电流id(k),iq(k)：
<mrow>
<msub>
<mi>i</mi>
<mi>d</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>2</mn>
<mn>3</mn>
</mfrac>
<mo>&lsqb;</mo>
<msub>
<mi>i</mi>
<mi>a</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mi>cos</mi>
<mo>(</mo>
<mrow>
<msub>
<mi>&theta;</mi>
<mi>e</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
<mo>+</mo>
<msub>
<mi>i</mi>
<mi>b</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mi>cos</mi>
<mo>(</mo>
<mrow>
<msub>
<mi>&theta;</mi>
<mi>e</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mfrac>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
</mrow>
<mn>3</mn>
</mfrac>
</mrow>
<mo>)</mo>
<mo>+</mo>
<msub>
<mi>i</mi>
<mi>c</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mi>e</mi>
</msub>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
<mo>+</mo>
<mfrac>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
</mrow>
<mn>3</mn>
</mfrac>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<msub>
<mi>i</mi>
<mi>q</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mn>2</mn>
<mn>3</mn>
</mfrac>
<mo>&lsqb;</mo>
<msub>
<mi>i</mi>
<mi>a</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mi>e</mi>
</msub>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>i</mi>
<mi>b</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mi>e</mi>
</msub>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
<mo>-</mo>
<mfrac>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
</mrow>
<mn>3</mn>
</mfrac>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>i</mi>
<mi>c</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mi>e</mi>
</msub>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
<mo>+</mo>
<mfrac>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
</mrow>
<mn>3</mn>
</mfrac>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>;</mo>
</mrow>
(7)计算PMSM在dq坐标下的交直轴电流误差的一阶变化量ded(k),deq(k)
ed(k-1)＝-id(k-1)
edz(k-1)＝zd0(k-1)-ed(k-1)
zd0(k)＝zd0(k-1)+Ts(zd1(k-1)-λd1|edz(k-1)|1/2sgn(edz(k-1)))
zd1(k)＝zd1(k-1)+Tsλd2sgn(-λd1|edz(k-1)|1/2sgn(edz(k-1)))
ded(k)＝zd1(k)
<mrow>
<msub>
<mi>e</mi>
<mi>q</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msubsup>
<mi>i</mi>
<mi>q</mi>
<mo>*</mo>
</msubsup>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>i</mi>
<mi>q</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
eqz(k-1)＝zq0(k-1)-eq(k-1)
zq0(k)＝zq0(k-1)+Ts(zq1(k-1)-λq1|eqz(k-1)|1/2sgn(eqz(k-1)))
zq1(k)＝zq1(k-1)+Tsλq2sgn(-λq1|eqz(k-1)|1/2sgn(eqz(k-1)))
deq(k)＝zq1(k)
其中，ed(k-1)和eq(k-1)分别为dq轴电流误差；edz(k-1)和eqz(k-1)分别为dq轴电流误差的估计误差；zd0(k)和zq0(k)分别为dq轴电流误差的估计值；zd1(k)和zq1(k)分别为dq轴电流误差一阶导数ded(k)、deq(k)的估计值；λd1与λd2分别为d轴微分估计器一阶、二阶增益；λq1与λq2分别为d轴微分估计器一阶、二阶增益；
(8)计算dq轴电流环的快速终端滑模面sd(k),sq(k)
<mrow>
<msub>
<mi>s</mi>
<mi>d</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>de</mi>
<mi>d</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>&alpha;</mi>
<mi>d</mi>
</msub>
<msub>
<mi>e</mi>
<mi>d</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>&beta;</mi>
<mi>d</mi>
</msub>
<msubsup>
<mi>e</mi>
<mi>d</mi>
<mrow>
<msub>
<mi>q</mi>
<mi>d</mi>
</msub>
<mo>/</mo>
<msub>
<mi>p</mi>
<mi>d</mi>
</msub>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>s</mi>
<mi>q</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>de</mi>
<mi>q</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>&alpha;</mi>
<mi>q</mi>
</msub>
<msub>
<mi>e</mi>
<mi>q</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>&beta;</mi>
<mi>q</mi>
</msub>
<msubsup>
<mi>e</mi>
<mi>q</mi>
<mrow>
<msub>
<mi>q</mi>
<mi>q</mi>
</msub>
<mo>/</mo>
<msub>
<mi>p</mi>
<mi>q</mi>
</msub>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
其中，αd与βd是d轴电流环快速终端滑模面sd(k)的系数，且αd＞0，βd＞0；qd与pd是奇数，pd＞qd＞0；αq与βq是q轴电流环快速终端滑模面sq(k)的系数，且αq＞0，βq＞0；qq与pq是奇数，pd＞qd＞0；
(9)计算的电流环d轴电压控制量ud(k)
<mrow>
<msub>
<mi>u</mi>
<mrow>
<mi>d</mi>
<mi>s</mi>
<mi>w</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>LT</mi>
<mi>s</mi>
</msub>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>k</mi>
</munderover>
<mrow>
<mo>(</mo>
<msub>
<mi>&alpha;</mi>
<mi>d</mi>
</msub>
<mi>d</mi>
<mi>e</mi>
<mo>(</mo>
<mi>j</mi>
<mo>)</mo>
<mo>+</mo>
<msub>
<mi>&beta;</mi>
<mi>d</mi>
</msub>
<mfrac>
<msub>
<mi>q</mi>
<mi>d</mi>
</msub>
<msub>
<mi>p</mi>
<mi>d</mi>
</msub>
</mfrac>
<mi>e</mi>
<msup>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>q</mi>
<mi>d</mi>
</msub>
<mo>-</mo>
<msub>
<mi>p</mi>
<mi>d</mi>
</msub>
<mo>)</mo>
<mo>/</mo>
<msub>
<mi>p</mi>
<mi>d</mi>
</msub>
</mrow>
</msup>
<mi>d</mi>
<mi>e</mi>
<mo>(</mo>
<mi>j</mi>
<mo>)</mo>
<mo>+</mo>
<msub>
<mi>&eta;</mi>
<mi>d</mi>
</msub>
<mi>sgn</mi>
<mo>(</mo>
<mrow>
<msub>
<mi>s</mi>
<mi>d</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
<mo>+</mo>
<msub>
<mi>k</mi>
<mi>d</mi>
</msub>
<msub>
<mi>s</mi>
<mi>d</mi>
</msub>
<mo>(</mo>
<mi>j</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>u</mi>
<mrow>
<mi>d</mi>
<mi>e</mi>
<mi>q</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msubsup>
<mi>Li</mi>
<mi>d</mi>
<mo>*</mo>
</msubsup>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>Ln</mi>
<mi>p</mi>
</msub>
<mi>&omega;</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>i</mi>
<mi>q</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>R</mi>
<mi>s</mi>
</msub>
<msub>
<mi>i</mi>
<mi>d</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
ud(k)＝udeq(k)+udsw(k)
<mrow>
<msub>
<mi>u</mi>
<mi>d</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>u</mi>
<mi>max</mi>
</msub>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>u</mi>
<mi>d</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>&GreaterEqual;</mo>
<msub>
<mi>u</mi>
<mi>max</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>u</mi>
<mi>d</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>u</mi>
<mi>min</mi>
</msub>
<mo><</mo>
<msub>
<mi>u</mi>
<mi>d</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo><</mo>
<msub>
<mi>u</mi>
<mi>max</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>u</mi>
<mi>min</mi>
</msub>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>u</mi>
<mi>d</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>&le;</mo>
<msub>
<mi>u</mi>
<mi>min</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中，umax和umin分别电压控制量ud(k)、uq(k)的最大值和最小值；L为PMSM电机dq轴交直轴电感；Rs为电机定子电阻；udsw(k)为滑模积分控制分量；udeq(k)为等效控制分量；ηd与kd分别为指数趋近率的系数，且ηd＞0，kd＞0；
(10)计算电流环q轴电压控制量uq(k)
<mrow>
<msub>
<mi>u</mi>
<mrow>
<mi>q</mi>
<mi>s</mi>
<mi>w</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>LT</mi>
<mi>s</mi>
</msub>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>k</mi>
</munderover>
<mrow>
<mo>(</mo>
<msub>
<mi>&alpha;</mi>
<mi>q</mi>
</msub>
<mi>d</mi>
<mi>e</mi>
<mo>(</mo>
<mi>j</mi>
<mo>)</mo>
<mo>+</mo>
<msub>
<mi>&beta;</mi>
<mi>q</mi>
</msub>
<mfrac>
<msub>
<mi>q</mi>
<mi>q</mi>
</msub>
<msub>
<mi>p</mi>
<mi>q</mi>
</msub>
</mfrac>
<mi>e</mi>
<msup>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>q</mi>
<mi>q</mi>
</msub>
<mo>-</mo>
<msub>
<mi>p</mi>
<mi>q</mi>
</msub>
<mo>)</mo>
<mo>/</mo>
<msub>
<mi>p</mi>
<mi>q</mi>
</msub>
</mrow>
</msup>
<mi>d</mi>
<mi>e</mi>
<mo>(</mo>
<mi>j</mi>
<mo>)</mo>
<mo>+</mo>
<msub>
<mi>&eta;</mi>
<mi>q</mi>
</msub>
<mi>sgn</mi>
<mo>(</mo>
<mrow>
<msub>
<mi>s</mi>
<mi>q</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
<mo>+</mo>
<msub>
<mi>k</mi>
<mi>q</mi>
</msub>
<msub>
<mi>s</mi>
<mi>q</mi>
</msub>
<mo>(</mo>
<mi>j</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>u</mi>
<mrow>
<mi>q</mi>
<mi>e</mi>
<mi>q</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msubsup>
<mi>Li</mi>
<mi>q</mi>
<mo>*</mo>
</msubsup>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>Ln</mi>
<mi>p</mi>
</msub>
<mi>&omega;</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>i</mi>
<mi>d</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>R</mi>
<mi>s</mi>
</msub>
<msub>
<mi>i</mi>
<mi>q</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>n</mi>
<mi>p</mi>
</msub>
<msub>
<mi>&psi;</mi>
<mi>f</mi>
</msub>
<mi>&omega;</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
uq(k)＝uqeq(k)+uqsw(k)
<mrow>
<msub>
<mi>u</mi>
<mi>q</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>u</mi>
<mi>max</mi>
</msub>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>u</mi>
<mi>q</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>&GreaterEqual;</mo>
<msub>
<mi>u</mi>
<mi>max</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>u</mi>
<mi>q</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>u</mi>
<mi>min</mi>
</msub>
<mo><</mo>
<msub>
<mi>u</mi>
<mi>q</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo><</mo>
<msub>
<mi>u</mi>
<mi>max</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>u</mi>
<mi>min</mi>
</msub>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>u</mi>
<mi>q</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>&le;</mo>
<msub>
<mi>u</mi>
<mi>min</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中，uqsw(k)为滑模积分控制分量；uqeq(k)为等效控制分量；ηq与kq分别为指数趋近率的系数，且ηω＞0，kω＞0；
(11)将dq轴电压控制量ud(k),uq(k)作为PMSM的逆park变换的输入来完成电机的k步控制。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710981166.XA CN107528514B (zh) | 2017-10-20 | 2017-10-20 | Pmsm调速系统的近似离散快速终端滑模控制方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710981166.XA CN107528514B (zh) | 2017-10-20 | 2017-10-20 | Pmsm调速系统的近似离散快速终端滑模控制方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107528514A true CN107528514A (zh) | 2017-12-29 |
| CN107528514B CN107528514B (zh) | 2020-03-03 |
Family
ID=60685481
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710981166.XA Active CN107528514B (zh) | 2017-10-20 | 2017-10-20 | Pmsm调速系统的近似离散快速终端滑模控制方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107528514B (zh) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109458336A (zh) * | 2018-12-13 | 2019-03-12 | 青岛海尔空调器有限总公司 | 用于控制单转子压缩机转速的方法 |
| CN109560736A (zh) * | 2018-12-18 | 2019-04-02 | 东南大学 | 基于二阶终端滑模的永磁同步电机控制方法 |
| CN109946973B (zh) * | 2019-04-12 | 2021-05-04 | 东北大学 | 一种结合快速终端滑模和线性滑模的组合滑模控制方法 |
| CN113746392A (zh) * | 2021-09-17 | 2021-12-03 | 广州市香港科大霍英东研究院 | 一种基于二维系统的伺服电机电流控制方法及控制系统 |
| CN118399830A (zh) * | 2024-06-28 | 2024-07-26 | 广汽埃安新能源汽车股份有限公司 | 基于终端滑模的电机控制方法、装置、电子设备和存储介质 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102035456A (zh) * | 2010-12-14 | 2011-04-27 | 长春工业大学 | 基于终端滑模的永磁同步电机直接转矩控制系统 |
| CN103312244A (zh) * | 2013-06-18 | 2013-09-18 | 中南林业科技大学 | 基于分段式滑模变结构的无刷直流电机直接转矩控制方法 |
| CN105450123A (zh) * | 2015-11-25 | 2016-03-30 | 浙江工业大学 | 一种基于神经网络的永磁同步电机混沌系统快速终端滑模控制方法 |
| CN104270054B (zh) * | 2014-10-24 | 2016-09-21 | 哈尔滨工业大学 | 基于相对阶的永磁同步电机Anti-reset Windup平滑非奇异终端滑模控制方法 |
| CN106549616A (zh) * | 2016-10-28 | 2017-03-29 | 江苏大学 | 永磁同步电机的变指数系数趋近律滑模变结构控制方法 |
-
2017
- 2017-10-20 CN CN201710981166.XA patent/CN107528514B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102035456A (zh) * | 2010-12-14 | 2011-04-27 | 长春工业大学 | 基于终端滑模的永磁同步电机直接转矩控制系统 |
| CN103312244A (zh) * | 2013-06-18 | 2013-09-18 | 中南林业科技大学 | 基于分段式滑模变结构的无刷直流电机直接转矩控制方法 |
| CN104270054B (zh) * | 2014-10-24 | 2016-09-21 | 哈尔滨工业大学 | 基于相对阶的永磁同步电机Anti-reset Windup平滑非奇异终端滑模控制方法 |
| CN105450123A (zh) * | 2015-11-25 | 2016-03-30 | 浙江工业大学 | 一种基于神经网络的永磁同步电机混沌系统快速终端滑模控制方法 |
| CN106549616A (zh) * | 2016-10-28 | 2017-03-29 | 江苏大学 | 永磁同步电机的变指数系数趋近律滑模变结构控制方法 |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109458336A (zh) * | 2018-12-13 | 2019-03-12 | 青岛海尔空调器有限总公司 | 用于控制单转子压缩机转速的方法 |
| CN109560736A (zh) * | 2018-12-18 | 2019-04-02 | 东南大学 | 基于二阶终端滑模的永磁同步电机控制方法 |
| CN109946973B (zh) * | 2019-04-12 | 2021-05-04 | 东北大学 | 一种结合快速终端滑模和线性滑模的组合滑模控制方法 |
| CN113746392A (zh) * | 2021-09-17 | 2021-12-03 | 广州市香港科大霍英东研究院 | 一种基于二维系统的伺服电机电流控制方法及控制系统 |
| CN118399830A (zh) * | 2024-06-28 | 2024-07-26 | 广汽埃安新能源汽车股份有限公司 | 基于终端滑模的电机控制方法、装置、电子设备和存储介质 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107528514B (zh) | 2020-03-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107528514B (zh) | Pmsm调速系统的近似离散快速终端滑模控制方法 | |
| Liu et al. | Combined speed and current terminal sliding mode control with nonlinear disturbance observer for PMSM drive | |
| CN110289795B (zh) | 一种电动汽车用永磁同步电机控制系统及控制方法 | |
| CN106549616B (zh) | 永磁同步电机的变指数系数趋近律滑模变结构控制方法 | |
| CN106788086B (zh) | 考虑输入饱和的异步电机命令滤波有限时间模糊控制方法 | |
| CN106788046B (zh) | 永磁同步电机命令滤波有限时间模糊控制方法 | |
| Castaneda et al. | Discrete-time neural sliding-mode block control for a DC motor with controlled flux | |
| CN110190795B (zh) | 一种永磁同步电机级联式鲁棒预测电流控制方法 | |
| CN107577149B (zh) | 一种采用分数阶快速终端滑模控制的随动控制方法 | |
| CN112422004B (zh) | 一种永磁同步电机弱磁控制模式下的扰动抑制方法 | |
| CN108448991B (zh) | 基于零阶保持离散化模型的永磁电机电流预测控制方法 | |
| Liu et al. | Improved backstepping control with nonlinear disturbance observer for the speed control of permanent magnet synchronous motor | |
| CN110112971B (zh) | 一种基于有限时间动态面的异步电动机位置跟踪控制方法 | |
| CN108390597A (zh) | 带有扰动观测器的永磁同步电机非线性预测控制器设计 | |
| CN105811826A (zh) | 一种感应电机新型趋近律滑模控制方法 | |
| CN112769367B (zh) | 永磁直线电机数据驱动离散迭代积分滑模控制方法及系统 | |
| CN112187127B (zh) | 一种永磁同步电机控制方法 | |
| CN112737435B (zh) | 一种基于t-s模糊滑模控制的步进电机抗干扰系统 | |
| CN114280944B (zh) | 一种具有输出约束的pmsm系统有限时间动态面控制方法 | |
| CN106788053B (zh) | 基于观测器的电动车永磁同步电机系统误差补偿控制方法 | |
| Meng et al. | An EPCH control strategy for complex nonlinear systems with actuator saturation and disturbances | |
| CN109873583B (zh) | 基于状态受限的永磁同步电机模糊位置跟踪控制方法 | |
| Qu et al. | Sliding-mode anti-disturbance speed control of permanent magnet synchronous motor based on an advanced reaching law | |
| Feng et al. | Flux estimation of induction motors using high-order terminal sliding-mode observer | |
| CN114257130B (zh) | 一种适用于双永磁电机系统的基于二次型价值函数的预测速度同步控制方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |