CN107065548A - 基于安排过渡过程的双向dc‑dc变换器控制方法 - Google Patents
基于安排过渡过程的双向dc‑dc变换器控制方法 Download PDFInfo
- Publication number
- CN107065548A CN107065548A CN201710232174.4A CN201710232174A CN107065548A CN 107065548 A CN107065548 A CN 107065548A CN 201710232174 A CN201710232174 A CN 201710232174A CN 107065548 A CN107065548 A CN 107065548A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mfrac
- mtd
- msup
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
- G05B13/042—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Dc-Dc Converters (AREA)
Abstract
本发明公开了一种基于安排过渡过程的双向DC‑DC变换器控制方法,包括采用系统辨识非参数建模描述双向DC‑DC变换器动态特性,将变换器的阶跃响应特性曲线等效为线性二阶模型,计算需要辨识的系统参数;对原始PID控制器差分方程中的微分项和积分项采用后向差分近似,安排过渡过程;利用双向DC‑DC变换器的过渡过程的加速度函数,最终获得安排的过程函数。本发明在利用辨识法获得变换器二阶线性模型传递函数的基础上,使用安排过渡过程的PID控制器对双向DC‑DC变换器控制,安排过渡过程可以在变换器初始运行阶段消除输出电压电流的超调量,一定程度上解决了系统超调量与过渡时间之间的矛盾问题,将达到稳定的时间作为控制标准,实现系统无超调地达到稳定状态。
Description
技术领域
本发明涉及一种双向DC-DC变换器控制方法,尤其涉及一种基于安排过渡过程的双向DC-DC变换器控制方法。
背景技术
随着新能源产业的迅猛发展,储能系统的重要性日益显著。作为储能系统的关键部分,双向DC-DC变换器为系统的安全运行提供保障,其控制方式在电能储存过程中极其重要。双向DC-DC变换器常用控制方式为双闭环PID控制,但往往会导致控制器鲁棒性较弱,储能系统超调量与调节时间之间存在矛盾。
DC-DC变换器是储能系统的关键组成部分,其拓扑类型众多。为了保证储能系统的正常运行,对双向DC-DC变换器的控制也尤为重要。
发明内容
发明目的:针对以上问题,本发明提出一种基于安排过渡过程的双向DC-DC变换器控制方法。
技术方案:为实现本发明的目的,本发明所采用的技术方案是:一种基于安排过渡过程的双向DC-DC变换器控制方法,具体包括以下步骤:
(1)采用系统辨识非参数建模描述双向DC-DC变换器动态特性,将变换器的阶跃响应特性曲线等效为线性二阶模型,计算需要辨识的系统参数;
(2)对原始PID控制器差分方程中的微分项和积分项采用后向差分近似,安排过渡过程;
(3)利用双向DC-DC变换器的过渡过程的加速度函数,最终获得安排的过程函数。
有益效果:本发明在利用辨识法获得变换器二阶线性模型传递函数的基础上,使用安排过渡过程的PID控制器对双向DC-DC变换器控制,安排过渡过程可以在变换器初始运行阶段消除输出电压电流的超调量,一定程度上解决了系统超调量与过渡时间之间的矛盾问题。
本发明采用基于安排过渡过程的新型双闭环PID控制方式,使储能系统在初始运行阶段消除输出超调量,将达到稳定的时间作为控制标准,实现系统无超调地达到稳定状态。
附图说明
图1是双向DC-DC变换器框图;
图2是二阶传递函数零状态单位阶跃响应曲线图;
图3是双向DC-DC变换器传递函数matlab/simulink仿真曲线图;
图4是PID控制原理图;
图5是PID控制系统结构图;
图6是具有安排过渡过程的PID控制器模型图;
图7是安排过渡过程的单位阶跃响应效果图;其中,7(a)、7(b)为一对比组;7(c)、7(d)为一对比组;
图8是安排过渡过程的方波跟踪效果图;其中,8(a)为无过渡过程;8(b)为安排过渡过程;
图9是目标2A输出电流阶跃响应图;其中,9(a)为无过渡过程PID控制;9(b)为安排过渡过程PID控制。
具体实施方式
下面结合附图和实施例对本发明的技术方案作进一步的说明。
如图1所示是双向DC-DC变换器框图,其中,Udc为双向DC-DC变换器的输出电压,Id为锂电池充放电电流。根据双向DC-DC变换器在不同工作状态下控制对象的不同,采用双闭环控制策略,外环1控制输出电压,内环2控制充放电电流,如图虚线框所示,3为双向DC-DC变换器输出方向选择项。
采用系统辨识非参数建模描述双向DC-DC变换器动态特性,将变换器运行在某激励条件下的阶跃响应特性曲线近似等效为线性二阶模型,通常二阶线性模型传递函数如式(1):
根据线性系统理论,纯延时不会改变系统特征参数。为使系统辨识过程方便简单,将系统作归一化处理,使式(1)变为式(2):
其中,K为系统增益,ζ为系统阻尼系数,ωn为系统自然频率,τ为系统延时。式(2)经过拉普拉斯变换得到其在时域上的单位阶跃响应微分方程如式(3):
其中,
二阶传递函数零状态单位阶跃响应曲线如图2所示。结合图2,根据二阶常系数微分方程特性,对式(3)求导,令导函数等于零可以求得:
其中,Mn是t=tn时刻y(t)的极大值。根据式(4)、(5)可以求出:
通过求解方程组可得需要辨识的系统参数,即自然频率ωn与系统阻尼系数ζ:
利用搭建好的双向DC-DC变换器驱动控制平台,测得双向DC-DC变换器动态响应曲线,增益归一化后如图3中虚线所示,图中可以看出,双向DC-DC变换器在约35ms时趋于稳定,前两次的超调量M0、M1,以及分别对应的时间t0、t1,代入式(7),可以求得该激励条件下的系统特征参数ωn及ζ。通过matlab/simulink对双向DC-DC变换器传递函数进行仿真,仿真曲线如图3中实线所示,与实测曲线比较可以看到系统辨识的近似度较高,拟合效果好,可以满足绝大多数应用场合。
PID控制是基于被控变量的实际值与期望值的偏差,利用偏差量纠正系统的响应,执行调节控制,包括比例、积分、微分控制三个环节,控制原理如图4所示。因此根据PID控制原理可得出如下公式:
其中,Kp、Ki、Kd分别为比例增益、积分增益和微分增益。比例环节的作用是加快系统的响应速度,提高系统的调节精度。积分环节的主要作用是消除系统的稳态误差。微分环节的主要作用是在响应过程中抑制系统偏差向任何方向进行变化,起到提前预知偏差变化的效果。将式(8)中的模拟信号数字化后,可得下式:
将式(9)与式(10)相减,可得到增量式PID的控制增量Δu(k):
Δu(k)=u(k)-u(k-1)
=Kp(e(k)-e(k-1))+Kie(k)+
Kd(e(k)-2e(k-1)+e(k-2)) (11)
通过DSP芯片控制增量Δu(k),误动作影响小,不会产生积分失控,容易获得较好的速度调节效果。
PID控制系统结构如图5所示,在PID控制器中,微分环节的作用是抑制系统偏差向各个方向变化,能够有效改善系统的动态性能,但出现外界干扰时,会使得系统变得非常敏感。在双向DC-DC变换器实际运行时,电网电压及负载存在突变情况,变换器在电压电流控制过程中的微分项会将此视为干扰,产生较大的误差变化率,如果Kd设置不当则可能引起系统失稳。此外,双向DC-DC变换器自身的响应比较快,调节时间短,因此更加需要对Kd进行合理设计。结合电压电流采样电路,DSP通过采样值计算双向DC-DC变换器当前输出电压电流并将其作为反馈量,在CPU中执行控制算法产生新的PWM单元配置参数,输出相应的PWM驱动信号,从而保证双向DC-DC变换器工作状态稳定或完成双向DC-DC变换器工作状态的调整。
对原始PID控制器差分方程中的微分项和积分项采用后向差分近似的方法,差分方程可表示为:
其中,U(k)是PID控制器k时刻的输出;e是被控制量目标值减去实际值的差值;T是控制周期;Ti是积分时间;Td是微分时间;Kp是比例系数;Ki=KpT/Ti是积分系数;Kd=KpTd/T是微分系数;Up(k)是与当前差值成比例的比例控制量;Ui(k)是与差值累积成比例的积分控制量;Ud(k)是与前后差值变化率成比例的微分控制量。
对于偏差比例控制:
Up(k)=Kpe(k)=Kp(S(k)-R(k)) (13)
其中,S(k)是逐步接近最终目标值S0的一系列过渡目标值。通常,过渡目标值与实际值的差S(k)-R(k)要比最终目标值与实际值的差S0-R(k)小很多,因此,需要较大的比例系数获得足够的控制量。当被控量低于过渡目标值时,Up(k)表现为正向激励;当被控量高于过渡目标值时,Up(k)表现为负向激励。Up(k)能够使得被控量跟踪过渡过程中各个目标值变化而变化。而在不安排过渡过程时,被控量始终与目标值比较,Up(k)的切换并不及时。因此,安排过渡过程,提高了Up(k)正负控制切换的及时性。
对于偏差积分控制:
其中,ε是设置了积分分离门限,k是控制周期。计算第k个控制周期的积分控制量时,如果前几个周期(l=0,1,…,j-1)偏差大于积分分离门限,则积分量为0,当后几个周期(l=j,j+1,…,k)偏差小于积分分离门限时,进行积分量的累积。如果没有安排过渡过程,被控量始终与最终目标值比较,积分控制只在被控量与最终目标差值小于积分分离门限时开始产生抑制作用,积分项的控制作用并不明显。安排过渡过程之后,被控量在每一个控制周期内都与逐步接近最终目标值的过渡目标值比较,如果偏差小于积分分离门限,积分控制起作用,参数设计的合理积分控制有可能全程都有效果。因此在安排了过渡过程之后,积分分离门限可以参考系统允许的稳态误差,设置小一些,积分系数也可以小一些。这样整个积分控制可以设计得比较温和。一般地动态误差主要靠比例项控制,静差则由积分控制消除。
对于偏差微分控制:
Ud(k)=Kd(e(k)-e(k-1))
=Kd[(S(k)-R(k))-(S(k-1)-R(k-1))] (15)
通常,微分控制总是抑制过程变化的趋势,即当没有安排过渡过程时,Ud(k)与R(k)的变化趋势是相反的:
Ud(k)=Kd[R(k-1)-R(k)] (16)
而如果安排了过渡过程,式(16)可以改写为:
Ud(k)=Kd(e(k)-e(k-1))
=Kd[(R(k)-R(k-1))-(S(k)-S(k-1))]=ΔR-ΔS (17)
式(17)中的偏差量由两部分组成,第一部分是实际被控量的增量ΔR,第二部分是过渡目标的增量ΔS。当第一部分增量ΔR与第二部分增量ΔS相等时,说明被控量的变化等于安排的目标变化量,此时微分控制为零;当ΔR大于ΔS时,说明被控量的变化大于安排的目标变化量,此时微分控制表现为反向抑制作用;当ΔS大于ΔR时,说明被控量的变化小于安排的目标变化量,此时微分控制表现为正向推进作用。可见,安排过渡过程后,微分控制更为精细,它不再是简单地抑制过程的变化,而是既能使被控量接近目标值,又能使被控量远离目标值。
从上述安排过渡过程的影响可以看出,PID参数的选定范围扩大了。比例系数可以选择临界震荡比例系数,大一些或者小一些都满足要求;积分系数和微分系数可以由系统允许的稳态误差及控制周期确定,即期望在多少个控制周期内消除很小的误差,与实验法、经验法获取PID参数规律比较,更为简单快捷。
对于双向DC-DC变换器的二阶线性模型传递函数,过渡过程的加速度函数V2(t):
其中,v0为设定值,T表示过渡过程时间。则速度函数V1(t)为加速度函数在时间t内的积分:
最终获得安排的过程函数V(t):
一般地,加速度函数V2(t)可以按照如下规则选取:在区间(0,T)的前一部分取正,后一部分取负,正负面积相等即可。由于V1(t)在(0,T)上是为正,因此安排的过程V(t)是从零开始单调递增至v0的函数,没有超调。
根据双向DC-DC变换器系统的驱动时间和所能承受快慢的能力,选取过渡时间T,v0则为单位阶跃函数。利用Matlab/simulink建立双向DC-DC变换器单位阶跃响应安排过渡过程后的PID控制器模型,如图6所示,MATLAB FUN中即为安排的过渡过程,其输入为系统仿真时间,输出值与单位阶跃函数的乘积是不断接近最终目标值的过渡目标值。
在同一PID参数控制下,仿真曲线对比如下图7所示,7(a)、7(b)为一对比组,7(c)、7(d)为一对比组。通过对比在图中可以看出,安排过渡过程之后,被控量一直跟随逐步接近最终目标值的过渡目标值,因此系统在过渡时间内没有超调量,可以选取较大的比例增益系数,使系统快速响应,缩短系统达到稳定状态的时间。这样一来,系统超调量大小与过渡时间长短这一矛盾得到了解决,只需要将达到稳定的时间作为控制标准,系统可以无超调的达到稳定状态。
同样的方法和步骤,首先选取过渡时间T,v0选择为单位方波函数,通过Matlab仿真,验证安排过渡过程对于双向DC-DC变换器开断的优化效果。在同一PID参数控制下,图8(a)是没有安排过渡过程的仿真曲线,从图中可以看出,在变换器开断阶段,存在一定的超调量,即输出电流切换的过渡时间内,电流不稳定;图8(b)是在0s、1s、2s三个时间点安排过渡过程后的仿真曲线,从图中可以看出,安排过渡过程后,变换器输出不再始终跟踪最终目标值,而是在安排的过渡时间内跟踪不断接近最终目标值的过渡目标值,因此在变换器开断阶段,输出没有产生超调的现象,而未安排过渡过程的其余时间,则存在超调量。
在仿真的基础上,利用双向DC-DC变换器的驱动平台,配合控制器及采样电路对变换器的恒流控制进行实际测试。通过设置驱动PWM波的频率和占空比,确定变换器输出电流。实验方案为变换器输出电流由1A向2A的阶跃响应,实验结果如图9(a)、(b)所示。从图中可以看出,安排过渡过程无法像仿真那样基本消除变换器输出电流的振荡,但可以在初始运行阶段,消除输出电流的超调量。
Claims (4)
1.一种基于安排过渡过程的双向DC-DC变换器控制方法,其特征在于:具体包括以下步骤:
(1)采用系统辨识非参数建模描述双向DC-DC变换器动态特性,将变换器的阶跃响应特性曲线等效为线性二阶模型,计算需要辨识的系统参数;
(2)对原始PID控制器差分方程中的微分项和积分项采用后向差分近似,安排过渡过程;
(3)利用双向DC-DC变换器的过渡过程的加速度函数,最终获得安排的过程函数。
2.根据权利要求1所述的基于安排过渡过程的双向DC-DC变换器控制方法,其特征在于:所述步骤(1)具体包括:
二阶线性模型传递函数如式:
将系统作归一化处理,变为:
其中,K为系统增益,ζ为系统阻尼系数,ωn为系统自然频率,τ为系统延时;经过拉普拉斯变换得到时域单位阶跃响应微分方程:
其中,
求导,令导函数等于零可得:
<mrow>
<msub>
<mi>t</mi>
<mi>n</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
<mi>&pi;</mi>
</mrow>
<mrow>
<msqrt>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msup>
<mi>&zeta;</mi>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<msub>
<mi>&omega;</mi>
<mi>n</mi>
</msub>
</mrow>
</mfrac>
<mo>,</mo>
<mi>n</mi>
<mo>=</mo>
<mn>0</mn>
<mo>,</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>...</mo>
<mo>...</mo>
</mrow>
<mrow>
<msub>
<mi>M</mi>
<mi>n</mi>
</msub>
<mo>=</mo>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>&zeta;</mi>
<mfrac>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
<mi>&pi;</mi>
</mrow>
<msqrt>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msup>
<mi>&zeta;</mi>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mfrac>
</mrow>
</msup>
<mo>,</mo>
<mi>n</mi>
<mo>=</mo>
<mn>0</mn>
<mo>,</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>...</mo>
<mo>...</mo>
</mrow>
其中,Mn是t=tn时刻y(t)的极大值;可以求出:
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mfrac>
<msub>
<mi>M</mi>
<mi>n</mi>
</msub>
<msub>
<mi>M</mi>
<mrow>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
</mfrac>
<mo>=</mo>
<msup>
<mi>e</mi>
<mfrac>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
<mi>&zeta;</mi>
</mrow>
<msqrt>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msup>
<mi>&zeta;</mi>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mfrac>
</msup>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>t</mi>
<mrow>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>t</mi>
<mi>n</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
</mrow>
<mrow>
<msub>
<mi>&omega;</mi>
<mi>n</mi>
</msub>
<msqrt>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msup>
<mi>&zeta;</mi>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
通过求解方程组可得需要辨识的系统参数,即自然频率ωn与系统阻尼系数ζ:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>&zeta;</mi>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<msqrt>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>&pi;</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mi>n</mi>
<mfrac>
<msub>
<mi>M</mi>
<mi>n</mi>
</msub>
<msub>
<mi>M</mi>
<mrow>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mfrac>
<mi>l</mi>
<mi>n</mi>
<mfrac>
<msub>
<mi>M</mi>
<mi>n</mi>
</msub>
<msub>
<mi>M</mi>
<mrow>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&omega;</mi>
<mi>n</mi>
</msub>
<mo>=</mo>
<mfrac>
<msqrt>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>&pi;</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mi>n</mi>
<mfrac>
<msub>
<mi>M</mi>
<mi>n</mi>
</msub>
<msub>
<mi>M</mi>
<mrow>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mrow>
<msub>
<mi>t</mi>
<mrow>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>t</mi>
<mi>n</mi>
</msub>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>.</mo>
</mrow>
1
3.根据权利要求2所述的基于安排过渡过程的双向DC-DC变换器控制方法,其特征在于:所述步骤(2)具体包括:
PID控制原理:
<mrow>
<mi>u</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>K</mi>
<mi>p</mi>
</msub>
<mi>e</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>K</mi>
<mi>i</mi>
</msub>
<msubsup>
<mo>&Integral;</mo>
<mi>o</mi>
<mi>t</mi>
</msubsup>
<mi>e</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>t</mi>
<mo>+</mo>
<msub>
<mi>K</mi>
<mi>d</mi>
</msub>
<mfrac>
<mrow>
<mi>d</mi>
<mi>e</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mfrac>
</mrow>
其中,Kp、Ki、Kd分别为比例增益、积分增益和微分增益;将式中模拟信号数字化后,可得下式:
<mrow>
<mi>u</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>K</mi>
<mi>p</mi>
</msub>
<mi>e</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>K</mi>
<mi>i</mi>
</msub>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>k</mi>
</msubsup>
<mi>e</mi>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>K</mi>
<mi>d</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>e</mi>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
<mo>-</mo>
<mi>e</mi>
<mo>(</mo>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>u</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>K</mi>
<mi>p</mi>
</msub>
<mi>e</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>K</mi>
<mi>i</mi>
</msub>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mi>e</mi>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>K</mi>
<mi>d</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>e</mi>
<mo>(</mo>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
<mo>-</mo>
<mi>e</mi>
<mo>(</mo>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>2</mn>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
将两式相减,可得到增量式PID的控制增量Δu(k):
Δu(k)=u(k)-u(k-1)
=Kp(e(k)-e(k-1))+Kie(k)+
Kd(e(k)-2e(k-1)+e(k-2))
PID控制器的差分方程为:
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>U</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>k</mi>
<mi>p</mi>
</msub>
<mi>e</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>K</mi>
<mi>p</mi>
</msub>
<mi>T</mi>
</mrow>
<msub>
<mi>T</mi>
<mi>i</mi>
</msub>
</mfrac>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>l</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>k</mi>
</munderover>
<mi>e</mi>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>K</mi>
<mi>p</mi>
</msub>
<msub>
<mi>T</mi>
<mi>d</mi>
</msub>
</mrow>
<mi>T</mi>
</mfrac>
<mrow>
<mo>(</mo>
<mi>e</mi>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
<mo>-</mo>
<mi>e</mi>
<mo>(</mo>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<msub>
<mi>K</mi>
<mi>p</mi>
</msub>
<mi>e</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>K</mi>
<mi>i</mi>
</msub>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>l</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>k</mi>
</munderover>
<mi>e</mi>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>K</mi>
<mi>d</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>e</mi>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
<mo>-</mo>
<mi>e</mi>
<mo>(</mo>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<msub>
<mi>U</mi>
<mi>p</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>U</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>U</mi>
<mi>d</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中,U(k)是PID控制器k时刻的输出;e是被控制量目标值减去实际值的差值;T是控制周期;Ti是积分时间;Td是微分时间;Kp是比例系数;Ki=KpT/Ti是积分系数;Kd=KpTd/T是微分系数;Up(k)是与当前差值成比例的比例控制量;Ui(k)是与差值累积成比例的积分控制量;Ud(k)是与前后差值变化率成比例的微分控制量;
对于偏差比例控制:
Up(k)=Kpe(k)=Kp(S(k)-R(k))
其中,S(k)是逐步接近最终目标值S0的一系列过渡目标值;
对于偏差积分控制:
其中,ε是设置积分分离门限,k是控制周期;
对于偏差微分控制:
Ud(k)=Kd(e(k)-e(k-1))
=Kd[(S(k)-R(k))-(S(k-1)-R(k-1))]
Ud(k)与R(k)的变化趋势是相反的:
Ud(k)=Kd[R(k-1)-R(k)]
安排过渡过程,为:
Ud(k)=Kd(e(k)-e(k-1))
=Kd[(R(k)-R(k-1))-(S(k)-S(k-1))]=ΔR-ΔS。
4.根据权利要求3所述的基于安排过渡过程的双向DC-DC变换器控制方法,其特征在于:所述步骤(3)具体包括:
双向DC-DC变换器的过渡过程的加速度函数V2(t):
<mrow>
<msub>
<mi>V</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mn>4</mn>
<mfrac>
<msub>
<mi>v</mi>
<mn>0</mn>
</msub>
<msup>
<mi>T</mi>
<mn>2</mn>
</msup>
</mfrac>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>t</mi>
<mo>&le;</mo>
<mfrac>
<mi>T</mi>
<mn>2</mn>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mn>4</mn>
<mfrac>
<msub>
<mi>v</mi>
<mn>0</mn>
</msub>
<msup>
<mi>T</mi>
<mn>2</mn>
</msup>
</mfrac>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mfrac>
<mi>T</mi>
<mn>2</mn>
</mfrac>
<mo><</mo>
<mi>t</mi>
<mo>&le;</mo>
<mi>T</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>0</mn>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>t</mi>
<mo><</mo>
<mi>T</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中,v0为设定值,T表示过渡过程时间;则速度函数V1(t):
<mrow>
<msub>
<mi>V</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msubsup>
<mo>&Integral;</mo>
<mn>0</mn>
<mi>t</mi>
</msubsup>
<msub>
<mi>V</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>&tau;</mi>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mn>4</mn>
<mfrac>
<msub>
<mi>v</mi>
<mn>0</mn>
</msub>
<msup>
<mi>T</mi>
<mn>2</mn>
</msup>
</mfrac>
<mi>t</mi>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>t</mi>
<mo>&le;</mo>
<mfrac>
<mi>T</mi>
<mn>2</mn>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mn>4</mn>
<mfrac>
<msub>
<mi>v</mi>
<mn>0</mn>
</msub>
<msup>
<mi>T</mi>
<mn>2</mn>
</msup>
</mfrac>
<mi>t</mi>
<mo>+</mo>
<mn>4</mn>
<mfrac>
<msub>
<mi>v</mi>
<mn>0</mn>
</msub>
<mi>T</mi>
</mfrac>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mfrac>
<mi>T</mi>
<mn>2</mn>
</mfrac>
<mo><</mo>
<mi>t</mi>
<mo>&le;</mo>
<mi>T</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>0</mn>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>t</mi>
<mo><</mo>
<mi>T</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
最终获得安排的过程函数V(t):
<mrow>
<mi>V</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msubsup>
<mo>&Integral;</mo>
<mn>0</mn>
<mi>t</mi>
</msubsup>
<msub>
<mi>V</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>&tau;</mi>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mn>4</mn>
<mfrac>
<msub>
<mi>v</mi>
<mn>0</mn>
</msub>
<msup>
<mi>T</mi>
<mn>2</mn>
</msup>
</mfrac>
<msup>
<mi>t</mi>
<mn>2</mn>
</msup>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>t</mi>
<mo>&le;</mo>
<mfrac>
<mi>T</mi>
<mn>2</mn>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mi>2</mi>
<mfrac>
<msub>
<mi>v</mi>
<mn>0</mn>
</msub>
<msup>
<mi>T</mi>
<mn>2</mn>
</msup>
</mfrac>
<msup>
<mi>t</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<mn>4</mn>
<mfrac>
<msub>
<mi>v</mi>
<mn>0</mn>
</msub>
<mi>T</mi>
</mfrac>
<mi>t</mi>
<mo>-</mo>
<msub>
<mi>v</mi>
<mn>0</mn>
</msub>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mfrac>
<mi>T</mi>
<mn>2</mn>
</mfrac>
<mo><</mo>
<mi>t</mi>
<mo>&le;</mo>
<mi>T</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>v</mi>
<mn>0</mn>
</msub>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>t</mi>
<mo><</mo>
<mi>T</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>.</mo>
</mrow>
3
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710232174.4A CN107065548A (zh) | 2017-04-11 | 2017-04-11 | 基于安排过渡过程的双向dc‑dc变换器控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710232174.4A CN107065548A (zh) | 2017-04-11 | 2017-04-11 | 基于安排过渡过程的双向dc‑dc变换器控制方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107065548A true CN107065548A (zh) | 2017-08-18 |
Family
ID=59602685
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710232174.4A Pending CN107065548A (zh) | 2017-04-11 | 2017-04-11 | 基于安排过渡过程的双向dc‑dc变换器控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107065548A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109711034A (zh) * | 2018-12-24 | 2019-05-03 | 西北工业大学 | 一种飞机起落架健康监测方法 |
CN111240192A (zh) * | 2018-11-28 | 2020-06-05 | 中国科学院沈阳自动化研究所 | 一种基于平滑切换函数的过渡过程目标值控制方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1812246A (zh) * | 2006-01-04 | 2006-08-02 | 天津大学 | 非接触液体媒质超声波电机频率跟踪转速控制系统 |
CN102694484A (zh) * | 2012-01-14 | 2012-09-26 | 河南科技大学 | 一种超声波电机混沌控制系统及控制方法 |
US8538572B2 (en) * | 2009-06-30 | 2013-09-17 | Lam Research Corporation | Methods for constructing an optimal endpoint algorithm |
CN103378724A (zh) * | 2013-07-25 | 2013-10-30 | 重庆大学 | DC-DC buck变换器高阶滑模控制方法 |
-
2017
- 2017-04-11 CN CN201710232174.4A patent/CN107065548A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1812246A (zh) * | 2006-01-04 | 2006-08-02 | 天津大学 | 非接触液体媒质超声波电机频率跟踪转速控制系统 |
US8538572B2 (en) * | 2009-06-30 | 2013-09-17 | Lam Research Corporation | Methods for constructing an optimal endpoint algorithm |
CN102694484A (zh) * | 2012-01-14 | 2012-09-26 | 河南科技大学 | 一种超声波电机混沌控制系统及控制方法 |
CN103378724A (zh) * | 2013-07-25 | 2013-10-30 | 重庆大学 | DC-DC buck变换器高阶滑模控制方法 |
Non-Patent Citations (3)
Title |
---|
吴笑天: "基于DDS技术的超声波电机驱动控制研究", 《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》 * |
郭赣等: "基于过渡过程安排的 PID 控制器设计", 《宇航学报》 * |
黄焕袍等: "安排过渡过程是提高闭环系统"鲁棒性、适应性和稳定性"的一种有效方法", 《控制理论与应用》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111240192A (zh) * | 2018-11-28 | 2020-06-05 | 中国科学院沈阳自动化研究所 | 一种基于平滑切换函数的过渡过程目标值控制方法 |
CN111240192B (zh) * | 2018-11-28 | 2022-02-01 | 中国科学院沈阳自动化研究所 | 一种基于平滑切换函数的过渡过程目标值控制方法 |
CN109711034A (zh) * | 2018-12-24 | 2019-05-03 | 西北工业大学 | 一种飞机起落架健康监测方法 |
CN109711034B (zh) * | 2018-12-24 | 2022-04-15 | 西北工业大学 | 一种飞机起落架健康监测方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Maiti et al. | An adaptive speed sensorless induction motor drive with artificial neural network for stability enhancement | |
Oliva et al. | Digital control of a voltage-mode synchronous buck converter | |
Efe | Integral sliding mode control of a quadrotor with fractional order reaching dynamics | |
Ahmad et al. | Unified disturbance-estimation-based control and equivalence with IMC and PID: Case study on a DC–DC boost converter | |
CN106950839A (zh) | 一种自抗扰控制结构的构造方法 | |
Papadopoulos et al. | Explicit analytical PID tuning rules for the design of type-III control loops | |
Aharon et al. | Uncertainty and disturbance estimator-based controllers design under finite control bandwidth constraint | |
CN105186863A (zh) | 一种基于Buck变换器的连续非奇异终端滑模控制方法 | |
CN107065548A (zh) | 基于安排过渡过程的双向dc‑dc变换器控制方法 | |
Hassan et al. | Brushless DC motor speed control using PID controller, fuzzy controller, and neuro fuzzy controller | |
Ramos-Paja et al. | Sliding-mode control of a CuK converter for voltage regulation of a dc-bus | |
Linares-Flores et al. | Sliding mode-delta modulation GPI control of a DC motor through a buck converter | |
Sabatta et al. | Super capacitor management using a Split-Pi symmetrical bi-directional DC-DC power converter with feed-forward gain control | |
CN106094508A (zh) | 基于δ算子的数字控制开关稳压电源的电压补偿器设计方法 | |
Viji et al. | Improved delta operator based discrete sliding mode fuzzy controller for buck converter | |
Ali et al. | Improved internal model control technique for position control of AC servo motors | |
Ghosh et al. | Study on chaos and bifurcation in DC-DC flyback converter | |
Yi et al. | Adaptive backstepping sliding mode nonlinear control for buck DC/DC switched power converter | |
CN112615546B (zh) | Llc谐振变换器高性能控制方法及系统 | |
Wang et al. | Robust smooth-trajectory control of nonlinear servo systems based on neural networks | |
Guedri et al. | Design of single closed loop control for chopper fed DC motor drive using IMC principles | |
Vladimir et al. | Three-loop control system of energy storage device in the frequency-controlled electric drive | |
Kim et al. | A DC Motor Speed Control System with Disturbance Rejection and Noise Reduction. | |
Abbas et al. | Design and FPGA implementation of 1-degree-of-freedom discrete PID controller for power switching converter | |
Sankeswari et al. | Sliding mode control for performance improvement of linear systems |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170818 |
|
RJ01 | Rejection of invention patent application after publication |