CN106605179B - 预测值整形系统及其方法以及控制系统和控制方法 - Google Patents

预测值整形系统及其方法以及控制系统和控制方法 Download PDF

Info

Publication number
CN106605179B
CN106605179B CN201580045459.6A CN201580045459A CN106605179B CN 106605179 B CN106605179 B CN 106605179B CN 201580045459 A CN201580045459 A CN 201580045459A CN 106605179 B CN106605179 B CN 106605179B
Authority
CN
China
Prior art keywords
value
control
predicted value
algorithm
target
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.)
Active
Application number
CN201580045459.6A
Other languages
English (en)
Other versions
CN106605179A (zh
Inventor
东俊一
南裕树
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
National Research Institute For Science And Technology
Original Assignee
National Research Institute For Science And Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by National Research Institute For Science And Technology filed Critical National Research Institute For Science And Technology
Publication of CN106605179A publication Critical patent/CN106605179A/zh
Application granted granted Critical
Publication of CN106605179B publication Critical patent/CN106605179B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B11/00Automatic controllers
    • G05B11/01Automatic controllers electric
    • G05B11/36Automatic controllers electric with provision for obtaining particular characteristics, e.g. proportional, integral, differential
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/04Adaptive 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
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/04Adaptive 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/048Adaptive 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 using a predictor

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Health & Medical Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Evolutionary Computation (AREA)
  • Medical Informatics (AREA)
  • Software Systems (AREA)
  • Feedback Control In General (AREA)

Abstract

本发明提供通过预测值的整形来计算更高精度的控制值的预测值整形系统。计算用于对控制对象(12)进行控制的控制值(v)的预测调节器(20)具备:实绩值获得部(22),其获得控制对象(12)的过去的目标值即实绩值(r(t‑1));预测值获得部(21),其获得预测控制对象(12)的目标值而得到的预测值(r^(t));以及控制值计算部(23),其通过将实绩值(r(t‑1))和预测值(r^(t))应用于预测值整形算法(G)修正预测值(r^(t)),来计算出用于对控制对象(12)进行控制的控制值(v(t))。预测值整形算法(G)使用控制对象(12)的控制模型(P)的参数。

Description

预测值整形系统及其方法以及控制系统和控制方法
关联的申请
本申请主张在2014年10月10日在日本国提交的专利申请号2014-209483的权利,并通过引用将该申请的内容纳入本申请。
技术领域
本发明涉及计算用于对控制对象进行控制的控制值的预测值整形系统和包含该预测值整形系统的控制系统、预测值整形方法和包含该预测值整形方法的控制方法、以及预测值整形程序。
背景技术
存在对被控制的对象(控制对象)使用控制器来进行前馈控制的前馈控制系统。在理想的前馈控制系统中,向系统赋予目标值作为用于对控制对象进行控制的控制值,并从控制对象输出与目标值相应的物理量。
图12是表示理想的前馈控制系统的构成的方框图。通过向包含控制器和控制对象的系统101的控制器赋予目标值r,控制器利用依照目标值r的驱动信号对控制对象进行驱动,从而从系统101得到与目标值r相应的输出yr
例如,对于发电系统,包含发电机及其控制器的发电系统成为系统101,构成上述的前馈控制系统。若作为目标值r向对发电机的发电量进行控制的控制器输入了电力需要量,则从发电机输出与电力需要大致相等的量的电力yr
在实际情况中,很多情况下,应向系统赋予的目标值(即,当前时间点需要的控制对象的输出量)是未知的,不得不将利用规定的预测模型预测目标值而得到的预测值赋予给系统。
例如,在上述的发电系统的情况下,不能准确地得到当前时间点的电力需要量(实际需要的电力量),因此,需要向发电系统赋予利用规定的预测模型得到的预测值。
图13是表示利用预测值来进行控制的控制系统的构成的方框图。向系统102赋予利用规定的预测模型得到的预测值
数式1:
Figure GDA0002013960160000011
此外,下面将
数式2:
Figure GDA0002013960160000021
也表示为“r^”。从系统102得到与该预测值r^相应的输出yr^。
然而,在利用预测值进行控制的情况下,与预测误差(即,预测值与实际的目标值之差)相应地,也在控制对象的输出中产生误差。因此,已知有为了使预测误差的影响不影响系统的输出而设置了预测调节器(预测值整形器)的前馈控制系统,该预测调节器将过去的时间点的目标值(即,在过去的时间点实际需要的控制值)作为实绩值来利用,对要赋予给系统的预测值进行修正(整形)。在此,实绩值不是过去的输出值,而是实际需要的控制值,实绩值与输出值独立。
图14是表示设置了预测调节器的控制系统的构成的方框图。向预测调节器201,赋予预测值r^和过去的时间点的目标值(实绩值)r。预测调节器201基于这些预测值r^和实绩值r计算出控制值v,并向系统103输出。在系统103中,得到与控制值v相应的输出yv
例如,在上述的发电系统的情况下,能够在预测调节器201中,例如使用几分钟前的实绩值r(即,几分钟前实际需要的控制值)对将天气、气温等应用于预测模型而得到的电力需要量的预测值r^进行修正,并向系统103赋予。利用这样的前馈控制系统进行同时同量的控制。
现有技术文献
专利文献
专利文献1:日本特开2006-288151号公报
发明内容
发明要解决的问题
本发明的目的在于提供通过预测值的整形来计算出精度更高的控制值的预测值整形系统、控制值计算方法、及控制值计算程序。
解决问题的方案
本发明一形态的预测值整形系统对预测用于对控制对象进行控制的目标值而得到的预测值进行修正,其具有如下构成,具备:预测值获得部,其获得所述预测值;实绩值获得部,其获得所述控制对象的过去的目标值即实绩值;以及控制值计算部,其通过将所述实绩值和所述预测值应用于算法修正所述预测值,来计算出用于对所述控制对象进行控制的控制值,所述算法使用所述控制对象的控制模型的参数。
根据该构成,利用使用了控制对象的控制模型的参数的算法来修正预测值,因此,能够进行适于控制对象的修正,能够得到使系统的输出更接近在系统用目标值被控制时的系统的输出的控制值。
在上述的预测值整形系统中,所述算法可以是以使过去的所述预测值与所述实绩值之间的差即过去的预测误差对所述控制对象的输出带来的影响较小的方式修正所述预测值的算法。
根据该构成,能够提供以使过去的预测值与实绩值之间的差(过去的预测误差)对控制对象的输出带来的影响较小的方式修正预测值的算法。
本发明另一形态的预测值整形系统对预测用于对控制对象进行控制的目标值而得到的预测值进行修正,其具有以下构成,具备:预测值获得部,其获得所述预测值;实绩值获得部,其获得所述控制对象的过去的目标值即实绩值;以及控制值计算部,其通过将所述实绩值和所述预测值应用于算法修正所述预测值,来计算出用于对所述控制对象进行控制的控制值,所述算法是以使过去的所述预测值与所述实绩值之间的差即过去的预测误差对所述控制对象的输出带来的影响较小的方式修正所述预测值的算法。
根据该构成,以使过去的预测值与实绩值之间的差(过去的预测误差)对控制对象的输出带来的影响较小的方式修正预测值,因此,能够进行适于控制对象的修正,能够得到使系统的输出更接近在系统用目标值被控制时的系统的输出的控制值。
在上述的预测值整形系统中,所述影响可以是应用了没有所述预测误差的情况下的目标值时的输出与实际的输出之间的输出差。
根据该构成,能够以使理想的输出(即,应用了没有预测误差的情况下的目标值时的输出)与实际的输出之间的输出差为较小的方式修正预测值。
上述的预测值整形系统可以还具备评价部,该评价部使用表现所述输出差的评价函数来评价所述影响的大小。
根据该构成,评价部评价过去的预测误差对所述控制对象的输出带来的影响的大小,因此,能够在影响变大的情况下检测出这种情况并采取重新计算算法等应对措施。
上述的预测值整形系统可以还具备算法决定部,该算法决定部决定所述算法。
根据该构成,能够动态地决定用于修正预测值的算法。
上述的预测值整形系统可以还具备评价部,该评价部评价所述影响的大小,所述算法决定部可以在由所述评价部评价出的所述影响的大小比规定的大小大时,决定所述算法。
根据该构成,能够在过去的所述预测值与所述实绩值之间的差即过去的预测误差对所述控制对象的输出带来的影响变大时,重新决定算法。
在上述的预测值整形系统中,所述算法可以是使所述影响为最小的算法。
根据该构成,能够得到使过去的预测误差对所述控制对象的输出带来的影响即输出差为最小的算法作为用于修正预测值的算法。
在上述的预测值整形系统中,所述算法可以是使所述影响比规定的值小的算法。
根据该构成,能够在过去的预测误差对所述控制对象的输出带来的影响即输出差变得比某规定值小时,得到用于修正预测值的算法。
在上述的预测值整形系统中,所述算法可以是判定为所述影响收敛了时的算法。
根据该构成,能够在过去的预测误差对所述控制对象的输出带来的影响即输出差收敛了时,得到用于修正预测值的算法。
在上述的预测值整形系统中,将所述控制模型的所述参数设为AP、BP、CP,所述控制模型可以通过式(1)表示。
数式3:
Figure GDA0002013960160000041
在此,v(t)是时刻t的所述控制值,y(t)是时刻t的所述控制对象的输出。另外,所述算法可以通过式(2)表示。
数式4:
Figure GDA0002013960160000042
在此,r^(t)是时刻t的所述预测值,r^(t-1)是与时刻t相比过去的所述预测值,r(t-1)是所述实绩值,其中,
数式5:
Figure GDA0002013960160000051
Figure GDA0002013960160000052
Figure GDA0002013960160000053
Figure GDA0002013960160000054
根据该构成,能够优化用于修正预测值的算法。
本发明一形态的控制系统具有以下构成,具备:上述的任意一项的预测值整形系统;所述控制对象;以及控制器,其按照由所述预测值整形系统计算出的所述控制值对所述控制对象进行控制。
根据该构成,利用使用了控制对象的控制模型的参数的算法来修正预测值,因此,也能够进行适于控制对象的修正,或者,以使过去的预测值与实绩值之间的差(过去的预测误差)对控制对象的输出带来的影响较小的方式修正预测值,因此,也能够进行适于控制对象的修正,能够得到使系统的输出更接近在系统用目标值被控制时的系统的输出的控制值。
本发明一形态的控制值计算方法是对预测用于对控制对象进行控制的目标值而得到的预测值进行修正的预测值整形方法,其具有以下构成,包括:预测值获得步骤,获得所述预测值;实绩值获得步骤,获得所述控制对象的过去的目标值即实绩值;以及控制值计算步骤,通过将所述实绩值和所述预测值应用于算法修正所述预测值,来计算出用于对所述控制对象进行控制的控制值,所述算法使用所述控制对象的控制模型的参数。
根据该构成,利用使用了控制对象的控制模型的参数的算法来修正预测值,因此,也能够进行适于控制对象的修正,能够得到使系统的输出更接近在系统用目标值被控制时的系统的输出的控制值。
本发明另一形态的控制值计算方法是对预测用于对控制对象进行控制的目标值而得到的预测值进行修正的预测值整形方法,其具有以下构成,包括:预测值获得步骤,获得所述预测值;实绩值获得步骤,获得所述控制对象的过去的目标值即实绩值;以及控制值计算步骤,通过将所述实绩值和所述预测值应用于算法修正所述预测值,来计算出用于对所述控制对象进行控制的控制值,所述算法是以使过去的所述预测值与所述实绩值之间的差即过去的预测误差对所述控制对象的输出带来的影响较小的方式修正所述预测值的算法。
根据该构成,以使过去的预测值与实绩值之间的差(过去的预测误差)对控制对象的输出带来的影响较小的方式修正预测值,因此,也能够进行适于控制对象的修正,能够得到使系统的输出更接近在系统用目标值被控制时的系统的输出的控制值。
本发明一形态的控制方法具有以下构成,包括:控制值计算步骤,利用上述的任意形态的预测值整形方法来计算控制值;以及控制步骤,利用通过所述控制值计算步骤计算出的所述控制值对所述控制对象进行控制。
根据该构成,利用使用了控制对象的控制模型的参数的算法来修正预测值,因此,也能够进行适于控制对象的修正,或者,以使过去的预测值与实绩值之间的差(过去的预测误差)对控制对象的输出带来的影响较小的方式修正预测值,因此,也能够进行适于控制对象的修正,能够得到使系统的输出更接近在系统用目标值被控制时的系统的输出的控制值。
本发明一形态的控制值计算程序具有以下构成,使对预测用于对控制对象进行控制的目标值而得到的预测值进行修正的预测值整形装置的计算机执行以下步骤:预测值获得步骤,获得所述预测值;实绩值获得步骤,获得所述控制对象的过去的目标值即实绩值;以及控制值计算步骤,通过将所述实绩值和所述预测值应用于算法修正所述预测值,来计算出用于对所述控制对象进行控制的控制值,所述算法使用所述控制对象的控制模型的参数。
根据该构成,利用使用了控制对象的控制模型的参数的算法来修正预测值,因此,也能够进行适于控制对象的修正,能够得到使系统的输出更接近在系统用目标值被控制时的系统的输出的控制值。
本发明另一形态的预测值整形程序具有以下构成,使对预测用于对控制对象进行控制的目标值而得到的预测值进行修正的预测值整形装置的计算机执行以下步骤:预测值获得步骤,获得所述预测值;实绩值获得步骤,获得控制对象的过去的目标值即实绩值;以及控制值计算步骤,通过将所述实绩值和所述预测值应用于算法修正所述预测值,来计算出用于对所述控制对象进行控制的控制值,所述算法是以使过去的所述预测值与所述实绩值之间的差即过去的预测误差对所述控制对象的输出带来的影响较小的方式修正所述预测值的算法。
根据该构成,以使过去的预测值与实绩值之间的差(过去的预测误差)对控制对象的输出带来的影响较小的方式修正预测值,因此,也能够进行适于控制对象的修正,能够得到使系统的输出更接近在系统用目标值被控制时的系统的输出的控制值。
发明效果
根据本发明,利用使用了控制对象的控制模型的参数的算法来修正(整形)预测值,因此,能够进行适于控制对象的修正,能够得到使系统的输出更接近在系统用目标值被控制时的系统的输出的控制值。
如以下说明的那样,本发明存在其他的形态。因此,该发明公开意图提供本发明的一部分,而不是意图限制在此描述要求的发明的范围。
附图简要说明
图1是表示本发明的实施方式涉及的包含预测值整形系统的控制系统的构成的方框图。
图2是表示本发明的实施方式涉及的前馈控制系统的构成的方框图。
图3是用于对本发明的实施方式涉及的评价函数进行说明的方框图。
图4是用于对本发明的实施方式涉及的评价函数进行说明的曲线图。
图5是表示本发明的实施方式的数值例中的目标值的曲线图。
图6是表示本发明的实施方式的数值例中的预测值的曲线图。
图7是表示使用预测值但不使用预测调节器的情况下的比较例的曲线图。
图8是表示本发明的实施方式涉及的使用预测调节器修正了预测值的情况下的数值例的曲线图。
图9是用于对本发明的实施方式涉及的前馈控制系统的效果进行说明的图。
图10是表示实现本发明的实施方式涉及的前馈控制系统的顺序的流程图。
图11是表示本发明的实施方式涉及的动态地决定预测值整形算法的预测值整形系统的构成的方框图。
图12是表示以往的理想的前馈的控制系统的构成的方框图。
图13是表示以往的利用预测值来进行控制的控制系统的构成的方框图。
图14是表示设置了预测调节器的控制系统的构成的方框图。
具体实施方式
下面,参照附图对本发明的实施方式的预测值整形系统进行说明。此外,下面说明的实施方式示出实施本发明时的一例,而不是将本发明限定于下面说明的具体构成。在实施本发明时,可以适宜采用根据实施方式的具体构成。
图1是表示本发明的实施方式涉及的含有预测值整形系统的控制系统的构成的方框图。控制系统100是具备预测调节器(预测值整形系统)20和系统10的前馈控制系统。系统10是具备控制对象12和对其进行控制的控制器11的反馈控制系统。即,在控制系统100中,预测调节器20与反馈控制系统10串联连接,并作为整体构成了前馈控制系统。控制器11接受控制对象12的输出值的反馈,并基于此输出值的反馈进行用于使控制对象12的输出成为从预测调节器20指定的控制值的反馈控制。具体而言,控制器11为了对控制对象12进行控制,而将作为控制信号的驱动信号向控制对象12输出。控制对象12按照从控制器11接收的驱动信号进行驱动,其结果,输出与驱动信号相应的物理量。
本实施方式中,用于对控制对象12进行控制的目标值是未知的。作为那样的控制对象12,例如有发电机。从发电机输出(发电)的电力具有以下特点:难以以电的形式直接储存,另外若供给低于需求则导致停电。因此,以符合时刻变动的电力需要量,产生确保了某种程度的余量的电力量来发电的方式,进行发电量控制。然而,若该余量过大则造成能源的浪費,因此,期望准确地估计需要量来进行同时同量的发电量控制。因此,本实施方式的控制值计算技术优选地适用于发电量控制。
作为预测值整形系统的预测调节器20具备:预测值获得部21、实绩值获得部22、控制值计算部23及输出部24。对于作为预测值整形系统的预测调节器20的构成要素即预测值获得部21、实绩值获得部22、控制值计算部23及输出部24,可以将一部分或全部分散地配置在网络上,也可以将一部分或全部配备在同一装置(预测值整形装置)中。
预测值获得部21获得通过利用规定的预测模型对控制对象12的目标值进行预测而得到的预测值r^。实绩值获得部22获得控制对象12的过去的目标值即实绩值r。控制值计算部23通过将实绩值r和预测值r^应用于规定的预测值整形算法修正(整形)预测值r^,来计算出用于对控制对象12进行控制的控制值v。输出部24将由控制值计算部23计算出的控制值v向系统10的控制器11输出。
若使用预测调节器20利用实绩值r对预测值r^进行了修正,则提高预测的精度,得到使系统10的输出接近系统10用目标值被控制时的系统10的输出的值得预测值r^。然而,预测调节器中的预测值整形算法不使用控制对象的信息,因此,没有针对个别的控制对象进行优化。换言之,预测调节器中的预测值整形算法只不过是对什么样的控制对象都拥有一定程度的效果而能够适用的通用的算法。
因此,本实施方式的预测调节器20使用基于系统10中的反馈控制的控制模型(下面,也称为“控制对象12的控制模型”。)而得到的预测值整形算法,来进行预测值的修正。下面,具体地对这样的预测值整形算法进行说明。
1.符号的准备
在说明预测值整形算法之前,首先,如以下那样定义各符号。
数式6:
Figure GDA0002013960160000091
分别表示实数的集合、正实数的集合、及自然数的集合。
将由
数式7:
Figure GDA0002013960160000092
Figure GDA0002013960160000093
表示的向量x和矩阵M的∞范数分别定义为数式8:
||x||:=maxi|xi|
Figure GDA0002013960160000094
将离散时间信号e的∞范数设为数式9:
Figure GDA0002013960160000095
对于由
数式10:
M:={Mij}
表示的矩阵M,设为,abs(M)为由矩阵M的各要素的绝对值构成的矩阵(设为abs(M)={Mij})。
2.问题设定
图2是表示本实施方式的前馈控制系统的构成的方框图。图2表示了由采用预测值整形算法G的预测调节器20、和利用控制模型P进行反馈控制的系统10构成的前馈控制系统Σ。在此,控制模型P是离散时间线性模型。另外,“z-1”表示获得一定时刻前的信息的操作。
控制模型P由下式(1)表示。
数式11:
Figure GDA0002013960160000101
在此,v(t)是预测调节器20输出的控制值,yv(t)是系统10的输出。另外,(Ap、Bp、Cp)是控制模型P的参数,根据输入输出数据使用以往的系统识别方法(例如,最小二乘法)决定(Ap、Bp、Cp)。
另外,
数式12:
Figure GDA0002013960160000102
表示离散时刻,
数式13:
Figure GDA0002013960160000103
表示状态,
数式14:
Figure GDA0002013960160000104
表示输入,
数式15:
Figure GDA0002013960160000105
表示输出,
数式16:
Figure GDA0002013960160000106
是常数矩阵。将AP的特征值设为存在于复平面上的单位圆内(设P为稳定)。另外,对于初始状态,根据
数式17:
Figure GDA0002013960160000107
给出
数式18:
xP(0)=xP0
另一方面,预测调节器20的预测值整形算法G由下面的式(2)表示。
数式19:
Figure GDA0002013960160000111
在此,r(t-1)是实绩值(一定时刻前的目标值),r^(t-1)是利用规定的预测模型而得到的一定时刻前的预测值,r^(t)是当前的预测值。另外,v(t)是当前的控制值(对预测值进行整形而得到的值)。另外,(A,B,C,D)是预测值整形算法的设计参数。
另外,
数式20:
Figure GDA0002013960160000112
是G的状态(ξ(0)=0),
数式21:
Figure GDA0002013960160000113
是输入,
数式22:
Figure GDA0002013960160000114
是输出,
数式23:
Figure GDA0002013960160000115
是常数矩阵。
信号r是实际的目标值。不能获得当前的值r(t),但是能够获得实绩值r(t-1)。另一方面,信号r^是预测值,因此,当前的值r^(t)和过去的值r^(t-1)中的任意一项都能够利用。对于该r和r^,设以下的假定(A1)。
数式24:
Figure GDA0002013960160000116
该假定(A1)意味着,预测误差为数式25:
Figure GDA0002013960160000117
以下,
即,使用一些估计方法以Δ以下的精度得到预测值r。
在此,为了得到控制值计算部23计算出使系统10的输出接近系统10用目标值被控制时的系统10的输出的值的控制值v的预测值整形算法,考虑图3、图4及式(3)所示的评价函数J(G)。图3是用于对评价函数J(G)进行说明的方框图,图4是用于对评价函数J(G)进行说明的曲线图。评价函数J(G)是在前馈控制系统Σ中,任意地给出控制对象12的控制模型
数式26:
Figure GDA0002013960160000121
时的评价函数,由下式给出。
数式27:
Figure GDA0002013960160000122
在此,yv是利用由本实施方式的预测调节器20计算出的控制值v进行了控制时的系统10的输出,yr是利用目标值r进行了控制时的系统10的输出。
如图3、图4及式(3)所示,评价函数J(G)表示本实施方式的前馈控制系统Σ的输出yv与能够利用目标值r的情况下的理想的前馈控制系统ΣI的输出yr之间的最大的差、即赋予了最差的预测值时的最大的输出差。换言之,评价函数J(G)是最差情形的评价指标。因此,使评价函数J(G)最小化的预测值整形算法G为最优的预测值整形算法。在下面,对求出使评价函数J(G)最小化的预测值整形算法G的、预测调节器的最优设计进行说明。
3.预测调节器的最优设计
3-1.性能分析
在给出使评价函数J(G)的值为较小的预测调节器20的预测值整形算法G的解之前,对已给出的预测值整形算法G的性能进行分析。然后,基于其结果,给出最优的预测值整形算法G。
首先,将式(1)和式(2)合并,则成为数式28:
Figure GDA0002013960160000123
因此,若设
数式29:
Figure GDA0002013960160000124
则前馈控制系统Σ能够如下面那样表示。
数式30:
Figure GDA0002013960160000131
接着,在时刻
数式31:
Figure GDA0002013960160000132
的Σ的输出yv与ΣI的输出yr之差以
数式32:
Figure GDA0002013960160000133
表示。
在此,设
数式33:
Figure GDA0002013960160000134
根据假定(A1),
数式34:
||w(t)||≤Δ。
此时,若注意到w(t)=0(t≦0),则下式成立,
数式35:
Figure GDA0002013960160000141
由此,
数式36:
||w(t)||≤Δ,
因此,得到下式。
数式37:
Figure GDA0002013960160000142
而且,存在上式的等号成立的预测值r^,因此,下式成立
数式38:
Figure GDA0002013960160000143
根据该式(3),能够预先估计由于预测误差Δ而对系统10的输出产生怎样程度的影响。并且,式(3)的第1项的值仅由P的参数决定,第2项的值由控制模型P和预测值整形算法G的参数决定。因此,由于全部的项是非负的值,所以,若以使第2项的值尽可能地小的方式决定预测值整形算法G,则能够得到使评价函数J(G)的值为较小的预测值整形算法G。
3-2.最优预测调节器
能够得到使式(3)的第2项为零的预测值整形算法。实际上,最优的预测值整形算法G为
数式39:
Figure GDA0002013960160000151
这时的评价函数的值为
数式40:
Figure GDA0002013960160000152
该式(4)是G的性能界限。
这样,本实施方式的预测调节器20的预测值整形算法G由系统10的控制模型P的参数被赋予特征。即,预测值整形算法G是针对控制模型P而决定的订制型的预测值整形算法。由此,能够准确地估计过去产生的预测误差r^(t-1)-r(t-1)对由控制模型P表示的系统10的输出带来的影响,并且,能够进行尽可能使该影响较小那样的修正。
另外,根据式(4)可知,若预测误差Δ变得较小,则J(G)的值变得较小,也就是,输出误差变得较小。并且,越是
数式41:
||CpBp||
小的控制模型P,则J(G)越小。如上述那样,控制模型P是离散时间系统,因此,将连续时间系统Pc离散化时的采样周期越短,则
数式42:
||CpBp||
越小。据此可知,采样周期越短,则越能够提高预测调节器20的性能。
4.数值例
下面,示出关于由下面的控制模型P控制的前馈控制系统(最小相位系统)的数值例。
数式43:
Figure GDA0002013960160000153
其中,初始值为数式44:
xP0:=[0.50]T
图5是表示目标值的曲线图,图6是表示针对图5的目标值的预测值的曲线图。作为目标值,使用r(t)=0.5sin(0.1πt)+0.5cos(0.03πt)。设图6所示的预测值的预测误差为Δ=1.6429。
图7是表示使用预测值r^但未使用预测调节器20的情况下的结果(控制对象12的输出)yr^的曲线图,图8是表示,使用采用被优化的预测值整形算法G的预测调节器20修正预测值r^来计算出控制值v,利用该控制值v对控制对象12进行控制而得到的结果(控制对象12的输出)yv的曲线图。在图7和图8中,虚线yr是使用了目标值r时的理想的输出。此外,J(G)=0.0049,
数式45:
||yv-yr||=0.0049
。参照图8,在利用被优化的预测值整形算法G进行了预测值r^的修正的情况下,系统10的输出yv与利用目标值r进行了控制的情况下的系统10的输出yr几乎一致。
如以上那样,本实施方式的前馈控制系统100中,预测调节器20的控制值计算部23中采用的预测值整形算法G使用控制对象12的控制模型P的参数,因此,能够与控制对象12相应地修正预测值。
图9是用于对本实施方式的前馈控制系统100的效果进行说明的图。本实施方式的预测调节器20的控制值计算部23中采用的预测值整形算法G能够准确地估计过去产生的预测误差(即,过去的预测值r^(t-1)与实绩值r(t-1)之间的差r^(t-1)-r(t-1))对控制对象12的输出带来的影响,能够以使该影响为较小的方式,来修正当前时间点的预测值r^(t)。
在预测调节器20中,使用预测值r^(t)、过去的预测值r^(t-1)及实绩值(过去的目标值)r(t-1),利用预测值整形算法G滤波在过去产生的预测误差r^(t-1)-r(t-1),来修正预测值r^(t)计算出控制值v(t),此时,作为该控制值v(t),能够得到使系统10的输出接近系统10用目标值r(t)被控制时的系统10的输出的值的控制值v(t),能够使系统10的输出的误差
数式46:
||yv-yr||
较小。此外,如图9所示,即使使用预测调节器20,也不能改善下一步骤的影响,但是,能够使其以后的影响较小。
下面,对实现本实施方式的前馈控制系统100的顺序进行说明。图10是表示实现本实施方式的前馈控制系统100的顺序的流程图。首先,通过系统识别,来求系统10的控制模型P的参数(Ap、Bp、Cp)(步骤S101)。
接着,使用控制模型P的参数(Ap、Bp、Cp),求预测值整形算法G的参数(A、B、C、D)(步骤S102)。然后,使用预测值整形算法G的参数(A、B、C、D),将采用了预测值整形算法G的预测调节器20作为程序装在计算机中(步骤S103)。
上述的实施方式中,进行预测调节器的最优设计,求出使评价函数J(G)最小化的预测值整形算法G,但是作为预测值整形算法G,不是一定要采用使评价函数J(G)为最小的算法。例如,在利用基于数学公式的变形的解析方法、一般的优化的方法(梯度法、牛顿法等)、或称为试探法的试探性的方法来求使评价函数J(G)为最小的预测值整形算法G的情况下,也可以采用评价函数J(G)变成比规定的阈值小时的预测值整形算法G。或者,在利用上述的梯度法、牛顿法等优化方法求预测值整形算法G的情况下,也可以采用判定为评价函数J(G)收敛了时的预测值整形算法G。即,若将预测值整形算法G的候补更新为G1、G2、…,则评价函数J(G)的值也如J(G1)、J(G2)、…进行变化,但是,在这时的J(G)的值的变化量的绝对值(|J(Gi)-J(G(i+1))|)变为ε以下(ε是事先决定的较小的正的值)的阶段视为收敛,也可以采用这时的预测值整形算法G。
另外,上述的实施方式中,作为评价函数J(G),使用实际的前馈控制系统Σ的输出yv与能够利用目标值r的情况下的理想的前馈控制系统ΣI的输出yr之间的最大的差、即赋予了最差的预测值时的最大的输出差,但是,成为该评价函数J(G)的输出差可以是将曲线图的线yv-yr的接近0的接近程度量化后的任意的函数,例如,也可以是输出yvと输出yr之间的差的平方范数。
上述的实施方式中,说明了使用求出预测值整形算法并采用其而得到的预测调节器20的情况,但是,预测值整形系统也可以追加用于动态地决定预测值整形算法的构成。即,预测值整形系统除了上述的预测调节器20的构成以外,还可以具备评价部和算法决定部。
图11是表示动态地决定预测值整形算法的预测值整形系统的构成的方框图。预测值整形系统20′除了上述的预测调节器20的构成以外,还具备评价部25和算法决定部26。即使在预测值整形系统20′中,可以将一部分或全部的构成分散地配置在网络上,也可以将一部分或全部的构成配备在同一装置(预测值整形装置)中。
评价部25输入实际的前馈控制系统Σ的输出yv,且对表现与使用了目标值r的理想的前馈控制系统ΣI的输出yr之间的输出差的评价函数J(G)的大小进行评价,由此评价过去的预测误差对控制对象12的输出带来的影响的大小。评价部25在评价函数J(G)比规定的阈值大的情况下,判断为需要重新决定预测值整形算法,并使算法决定部26重新决定预测值整形算法。
算法决定部26在从评价部25接受到预测值整形算法的重新决定的指示的情况下,如上述说明的那样,求使评价函数J(G)最小化的预测值整形算法G。此外,在该情况下,不是一定要使评价函数J(G)为最小,也可以将使评价函数J(G)比规定的阈值小的G作为预测值整形算法,还可以将判断为评价函数J(G)收敛了时的G作为预测值整形算法。
(应用例)
上述中,对本实施方式的前馈控制系统100能够应用于预测电力需要量的发电系统的情况进行了说明,但是本发明的控制系统不限于发电系统,也能够应用于其他的系统。在系统的目标值是未知的情况下,能够有效地应用本发明。
作为目标值是未知的系统,除了上述的发电系统以外,例如有太阳能发电系统中的进行日射量预测的系统、使用了响应性差的传感器的系统(例如,控制发动机的系统)、利用网络通信发送传感器信息的系统、及将传感器信息进行处理并使用的系统(例如,如机器人控制或汽车的自动驾驶控制那样,在所获得的图像的处理上花费时间的系统),也能够将本发明应用于这些系统。
另外,在人工胰脏系统中也能应用本发明。在人工胰脏系统中,利用在皮下组织安装的传感器进行CGM(持续血糖测定)。CGM中测定组织间液中的糖浓度,糖浓度与血糖值存在一定的相关性,但是已知糖浓度实际上比血糖值延迟10~15分钟。
另外,例如,将控制对象设为供应链,将控制目标设为仓库的库存管理,在物流领域也能够应用本发明。在该情况下,作为目标值乃至预测值,可以使用满意度,作为实绩信息(过去的值),可以使用基于客户调查的满意度。
另外,例如,将控制对象设为胰岛素用药,将控制目标设为血糖值的控制,在医疗领域也能够应用本发明。在这种情况下,作为目标值乃至预测值,可以使用血糖值,作为实绩信息(过去的值),可以使用基于体液测量的10分钟前的血糖值。
另外,例如,将控制对象设为汽车的辅助驾驶乃至自动驾驶,将控制目标设为线路跟踪控制,在ITS领域也能够应用本发明。在该情况下,作为目标值乃至预测值,可以使用环境信息(障碍物、信号等信息),作为实绩信息(过去的值),可以使用花时间而得到的详细的信息。
另外,例如,将控制对象设为植物的计划生产,将控制目标设为生长环境的控制,在农业领域也能够应用本发明。在该情况下,作为目标值乃至预测值,可以使用植物的生长状况,作为实绩信息(过去的值),可以使用实际的生长状况的信息。
另外,例如,将控制对象设为空调控制,将控制目标设为运转时间和运转台数的决定,在建筑领域也能够应用本发明。在该情况下,作为目标值乃至预测值,可以使用环境信息(温度、天气等信息),作为实绩信息(过去的值),可以使用实际的环境信息。
以上对在当前时间点考虑到的本发明的适宜的实施方式进行了说明,但是,针对本实施方式能够进行多种的变形,并且,附加的权利要求包含本发明的真实的精神和范围内的全部的变形。
工业实用性
本发明利用使用了控制对象的控制模型的参数的算法来修正(整形)预测值,因此具有能够进行适于控制对象的修正、能够得到更接近理想输出的输出之类的效果,作为计算用于对控制对象进行控制的控制值的预测值整形系统等是有用的。
符号说明
100 前馈控制系统(控制系统)
10 系统
11 控制器
12 控制对象
20 预测调节器(预测值整形系统)
21 预测值获得部
22 实绩值获得部
23 控制值计算部
24 输出部
25 评价部
26 算法决定部

Claims (10)

1.一种预测值整形系统,其对预测用于对控制对象进行控制的目标值而得到的预测值进行修正,其特征在于,
具备:
预测值获得部,其获得所述预测值;
实绩值获得部,其获得所述控制对象的过去的目标值即实绩值;以及
控制值计算部,其通过将所述实绩值和所述预测值应用于算法修正所述预测值,来计算出用于对所述控制对象进行控制的控制值,
所述算法使用所述控制对象的控制模型的参数,并且/或者所述算法是以使过去的所述预测值与所述实绩值之间的差即过去的预测误差对所述控制对象的输出带来的影响较小的方式修正所述预测值的算法,其中所述影响是应用了没有所述预测误差的情况下的目标值时的输出与实际的输出之间的输出差,
在所述算法使用所述控制对象的控制模型的参数的情况下,所述控制模型的所述参数是AP、BP、CP
所述控制模型通过式(1)表示,
数式1:
Figure FDA0002522053510000011
在此,v(t)是时刻t的所述控制值,y(t)是时刻t的所述控制对象的输出,所述算法通过式(2)表示,
数式2:
Figure FDA0002522053510000012
在此,r^(t)是时刻t的所述预测值,r^(t-1)是与时刻t相比过去的所述预测值,r(t-1)是所述实绩值,其中,
数式3:
Figure FDA0002522053510000021
Figure FDA0002522053510000022
Figure FDA0002522053510000023
Figure FDA0002522053510000024
2.如权利要求1所述的预测值整形系统,其特征在于,
所述算法是以使过去的所述预测值与所述实绩值之间的差即过去的预测误差对所述控制对象的输出带来的影响较小的方式修正所述预测值的算法,还具备评价部,该评价部使用表现所述输出差的评价函数来评价所述影响的大小。
3.如权利要求1所述的预测值整形系统,其特征在于,
还具备算法决定部,该算法决定部决定所述算法。
4.如权利要求3所述的预测值整形系统,其特征在于,
所述算法是以使过去的所述预测值与所述实绩值之间的差即过去的预测误差对所述控制对象的输出带来的影响较小的方式修正所述预测值的算法,还具备评价部,该评价部评价所述影响的大小,
所述算法决定部在由所述评价部评价出的所述影响的大小比规定的大小大时,决定所述算法。
5.如权利要求1所述的预测值整形系统,其特征在于,
所述算法是以使过去的所述预测值与所述实绩值之间的差即过去的预测误差对所述控制对象的输出带来的影响较小的方式修正所述预测值的算法,并且所述算法是使所述影响为最小的算法。
6.如权利要求1所述的预测值整形系统,其特征在于,
所述算法是以使过去的所述预测值与所述实绩值之间的差即过去的预测误差对所述控制对象的输出带来的影响较小的方式修正所述预测值的算法,并且所述算法是使所述影响比规定的值小的算法。
7.如权利要求1所述的预测值整形系统,其特征在于,
所述算法是以使过去的所述预测值与所述实绩值之间的差即过去的预测误差对所述控制对象的输出带来的影响较小的方式修正所述预测值的算法,并且所述算法是判定为所述影响收敛了时的算法。
8.一种利用预测值整形的控制系统,其特征在于,
具备:
权利要求1至7中任意一项记载的预测值整形系统,其对预测用于对控制对象进行控制的目标值而得到的预测值进行修正;
所述控制对象;以及
控制器,其按照由所述预测值整形系统计算出的控制值对所述控制对象进行控制。
9.一种预测值整形方法,其对预测用于对控制对象进行控制的目标值而得到的预测值进行修正,其特征在于,
包括:
预测值获得步骤,获得所述预测值;
实绩值获得步骤,获得所述控制对象的过去的目标值即实绩值;以及
控制值计算步骤,通过将所述实绩值和所述预测值应用于算法修正所述预测值,来计算出用于对所述控制对象进行控制的控制值,
所述算法使用所述控制对象的控制模型的参数,并且/或者所述算法是以使过去的所述预测值与所述实绩值之间的差即过去的预测误差对所述控制对象的输出带来的影响较小的方式修正所述预测值的算法,其中所述影响是应用了没有所述预测误差的情况下的目标值时的输出与实际的输出之间的输出差,
在所述算法使用所述控制对象的控制模型的参数的情况下,所述控制模型的所述参数是AP、BP、CP
所述控制模型通过式(1)表示,
数式1:
Figure FDA0002522053510000031
在此,v(t)是时刻t的所述控制值,y(t)是时刻t的所述控制对象的输出,所述算法通过式(2)表示,
数式2:
Figure FDA0002522053510000032
在此,r^(t)是时刻t的所述预测值,r^(t-1)是与时刻t相比过去的所述预测值,r(t-1)是所述实绩值,其中,
数式3:
Figure FDA0002522053510000041
Figure FDA0002522053510000042
Figure FDA0002522053510000043
Figure FDA0002522053510000044
10.一种利用预测值整形的控制方法,其特征在于,
包括:
控制值计算步骤,利用权利要求9中记载的预测值整形方法来计算控制值;以及
控制步骤,利用通过所述控制值计算步骤计算出的所述控制值对所述控制对象进行控制。
CN201580045459.6A 2014-10-10 2015-09-14 预测值整形系统及其方法以及控制系统和控制方法 Active CN106605179B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP2014209483 2014-10-10
JP2014-209483 2014-10-10
PCT/JP2015/004670 WO2016056176A1 (ja) 2014-10-10 2015-09-14 予測値整形システム、制御システム、予測値整形方法、制御方法、及び予測値整形プログラム

Publications (2)

Publication Number Publication Date
CN106605179A CN106605179A (zh) 2017-04-26
CN106605179B true CN106605179B (zh) 2020-09-25

Family

ID=55652816

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201580045459.6A Active CN106605179B (zh) 2014-10-10 2015-09-14 预测值整形系统及其方法以及控制系统和控制方法

Country Status (5)

Country Link
US (1) US10372089B2 (zh)
EP (1) EP3206094B1 (zh)
JP (1) JP6632538B2 (zh)
CN (1) CN106605179B (zh)
WO (1) WO2016056176A1 (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6859633B2 (ja) * 2016-09-05 2021-04-14 オムロン株式会社 モデル予測制御装置、モデル予測制御装置の制御方法、情報処理プログラム、および記録媒体
JP2018041150A (ja) * 2016-09-05 2018-03-15 オムロン株式会社 モデル予測制御装置、モデル予測制御装置の制御方法、情報処理プログラム、および記録媒体
JP6946861B2 (ja) * 2017-08-29 2021-10-13 トヨタ自動車株式会社 自動運転評価装置及び自動運転評価方法
JP7151546B2 (ja) * 2019-02-25 2022-10-12 オムロン株式会社 制御装置、制御方法、及び制御プログラム
US11740268B2 (en) * 2020-05-01 2023-08-29 University Of South Carolina Predictive zero phase filtering

Family Cites Families (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4539633A (en) * 1982-06-16 1985-09-03 Tokyo Shibaura Denki Kabushiki Kaisha Digital PID process control apparatus
JPH06102939A (ja) * 1992-09-22 1994-04-15 Kazuyuki Aihara 上水需要量予測装置
JP3641918B2 (ja) * 1997-10-09 2005-04-27 富士電機システムズ株式会社 電力需要量予測値補正方法
JP2000066731A (ja) * 1998-08-26 2000-03-03 Mitsubishi Electric Corp 上水道システムの運用制御システム及び運用制御装置
US6760716B1 (en) * 2000-06-08 2004-07-06 Fisher-Rosemount Systems, Inc. Adaptive predictive model in a process control system
US6616759B2 (en) * 2001-09-06 2003-09-09 Hitachi, Ltd. Method of monitoring and/or controlling a semiconductor manufacturing apparatus and a system therefor
JP2004295226A (ja) * 2003-03-25 2004-10-21 Matsushita Electric Works Ltd 需要量予測支援システム及びそのプログラム並びにそのプログラムを記録したコンピュータで読み取り可能な記録媒体
JP2006288151A (ja) 2005-04-04 2006-10-19 Kansai Electric Power Co Inc:The 電力需給制御システム、電力需給制御方法、および、記録媒体
EP2299338B1 (en) * 2008-07-07 2013-05-29 Honda Motor Co., Ltd. Controller
CN101634837A (zh) * 2009-08-17 2010-01-27 浙江大学 空分装置制氩系统氮塞的防控方法
CN102444784B (zh) * 2011-11-14 2013-11-20 上海金自天正信息技术有限公司 基于动态矩阵控制的钢铁企业蒸汽管网压力控制系统

Also Published As

Publication number Publication date
EP3206094A4 (en) 2018-08-15
WO2016056176A1 (ja) 2016-04-14
JP6632538B2 (ja) 2020-01-22
US20170235285A1 (en) 2017-08-17
CN106605179A (zh) 2017-04-26
US10372089B2 (en) 2019-08-06
JPWO2016056176A1 (ja) 2017-07-27
EP3206094A1 (en) 2017-08-16
EP3206094B1 (en) 2021-06-09

Similar Documents

Publication Publication Date Title
CN106605179B (zh) 预测值整形系统及其方法以及控制系统和控制方法
Noel et al. Control of a nonlinear liquid level system using a new artificial neural network based reinforcement learning approach
US11294339B2 (en) Control device, control method, and non-transitory recording medium
US20220326664A1 (en) Improved machine learning for technical systems
CN110688920B (zh) 一种无人驾驶控制方法、装置及服务器
US6801810B1 (en) Method and device for state estimation
CN113874865A (zh) 借助于贝叶斯优化方法确定技术系统的调节策略的模型参数的方法和装置
EP3200038A1 (en) Model evaluation device, model evaluation method, and program recording medium
CN114761966A (zh) 用于以轨迹为中心的基于模型的强化学习的稳健优化的系统及方法
US20200151545A1 (en) Update of attenuation coefficient for a model corresponding to time-series input data
Dutka et al. Optimized discrete-time state dependent Riccati equation regulator
US20040181300A1 (en) Methods, apparatus and computer program products for adaptively controlling a system by combining recursive system identification with generalized predictive control
KR102326733B1 (ko) 엑츄에이터 조절 시스템을 작동시키기 위한 방법 및 장치, 컴퓨터 프로그램 및 기계 판독가능한 저장 매체
Lakshmanan et al. Estimation and model predictive control of non-linear batch processes using linear parameter varying models
Zong et al. Reinforced iLQR: A sample-efficient robot locomotion learning
US11231694B2 (en) Information processing device, information processing method, and non-transitory recording medium
Cheng et al. DiffTune $^+ $: Hyperparameter-Free Auto-Tuning using Auto-Differentiation
CN106033189A (zh) 飞行机器人位姿神经网络预测控制器
Viswanadhapalli et al. Deep reinforcement learning with reward shaping for tracking control and vibration suppression of flexible link manipulator
US20160321087A1 (en) Methods of increasing processing speed in a processing system that performs a nonlinear optimization routine
Coen et al. Cruise control using model predictive control with constraints
Zhang et al. Controller Design for Time-Varying Sampling, Co-Regulated Systems
Fang et al. An integrated model predictive control scheme with disturbance preview
CN113985900B (zh) 一种四旋翼无人机姿态动态特性模型、辨识方法及自适应柔化预测控制方法
US11631030B2 (en) Learning with moment estimation using different time constants

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