CN106019949A - 一种自适应阶数分数阶模糊PIλ控制器方法 - Google Patents

一种自适应阶数分数阶模糊PIλ控制器方法 Download PDF

Info

Publication number
CN106019949A
CN106019949A CN201610633219.4A CN201610633219A CN106019949A CN 106019949 A CN106019949 A CN 106019949A CN 201610633219 A CN201610633219 A CN 201610633219A CN 106019949 A CN106019949 A CN 106019949A
Authority
CN
China
Prior art keywords
omega
controller
fuzzy
fractional order
order
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
Application number
CN201610633219.4A
Other languages
English (en)
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.)
Xian University of Technology
Original Assignee
Xian University of 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 Xian University of Technology filed Critical Xian University of Technology
Priority to CN201610633219.4A priority Critical patent/CN106019949A/zh
Publication of CN106019949A publication Critical patent/CN106019949A/zh
Pending legal-status Critical Current

Links

Classifications

    • 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

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)
  • Feedback Control In General (AREA)

Abstract

本发明公开了一种自适应阶数分数阶模糊PIλ控制器方法,具体按照以下步骤实施:步骤1:针对被控对象GP,确定分数阶PIλ控制器的三个控制参数;步骤2:利用模糊集合理论,设立模糊规则以及隶属度函数,将建立的模糊控制器加入到分数阶PIλ控制器,用于修正参数;步骤3:将经过模糊控制器修正的分数阶控制器连入控制系统中,模糊控制器根据输入实时修正PIλ控制器参数,达到需要的控制效果。本发明一种自适应阶数分数阶模糊PIλ控制器方法,响应速度快,超调量小,对扰动具有较好的鲁棒性。

Description

一种自适应阶数分数阶模糊PIλ控制器方法
技术领域
本发明属于自动控制原理技术领域,具体涉及一种自适应阶数分数阶模糊PIλ控制器方法。
背景技术
PI控制器,全称为比例积分控制器,其在工业控制中运用广泛,但对不良工况处理能力有限,近年来研究多集中在改进方面。
分数阶微积分可以认为是经典整数阶微积分的拓展,但它却有自己独特的逻辑和语法规则。常用的分数阶微积分定义有:Riemann-Liouville(RL)定义、Grunwald-Letnikov(GL)定义和Caputo定义。将PI控制器从整数阶拓展至分数阶,可以显著提高控制器的自由度,更好的控制被控系统。
利用模糊数学的基本思想和理论的控制方法。在传统的控制领域里,控制系统动态模式的精确与否是影响控制优劣的最主要关键,系统动态的信息越详细,则越能达到精确控制的目的。然而,对于复杂的系统,由于变量太多,往往难以正确的描述系统的动态,于是工程师便利用各种方法来简化系统动态,以达成控制的目的,但却不尽理想。换言之,传统的控制理论对于明确系统有强而有力的控制能力,但对于过于复杂或难以精确描述的系统,则显得无能为力了。因此便尝试着以模糊数学来处理这些控制问题。
发明内容
本发明的目的是提供一种自适应阶数分数阶模糊PIλ控制器方法,解决了现有分数阶PIλ控制器无法针对控制系统动态调整参数的问题。
本发明所采用的技术方案是,一种自适应阶数分数阶模糊PIλ控制器方法,具体按照以下步骤实施:
步骤1:针对被控对象GP,确定分数阶PIλ控制器的三个控制参数;
步骤2:利用模糊集合理论,设立模糊规则以及隶属度函数,将建立的模糊控制器加入到分数阶PIλ控制器,用于修正参数;
步骤3:将经过模糊控制器修正的分数阶控制器连入控制系统中,模糊控制器根据输入实时修正PIλ控制器参数,达到需要的控制效果。
本发明的特点还在于:
步骤1具体为:
步骤1.1:根据被控对象GP需要,确定幅值裕量Am和相位裕量φm
步骤1.2:根据幅值裕量和相位裕量来设计分数阶PIλ控制器以满足系统的实际要求,即:
φm=arg[Gc(jωc)Gp(jωc)]+π (1)
A m = 1 G c ( jω c ) G p ( jω c ) - - - ( 2 )
d { arg [ G c ( j ω ) G p ( j ω ) ] } d ω | ω = ω c = 0 - - - ( 3 )
控制器传递函数Gc为:
G c = K p + K i s λ - - - ( 4 )
将控制器传递函数Gc和被控对象传递函数GP改写为频域模式后代入公式(1)-公式(3)中,得到控制器的比例增益系数Kp、积分增益系数Ki及积分环节系数λ。
步骤1.2中sλ的频域表达式为:
s α = K Π k = - N N s + ω k ′ s + ω k - - - ( 5 )
其中,N为滤波器阶数,频域上限为ωb,下限为ωh,ωb·ωh=1,ωk、ωk′分别为极点和零点,K为增益系数:
K = ( ω h / ω b ) - α / 2 Π k = - N N ω k ′ ω k - - - ( 6 )
ω k = ω b ( ω h / ω b ) K + N + ( 1 + α ) / 2 ( 2 N + 1 ) - - - ( 7 )
ω k ′ = ω b ( ω h / ω b ) K + N + ( 1 - α ) / 2 ( 2 N + 1 ) - - - ( 8 ) .
步骤2具体为:
设误差e及误差变化率ec的模糊子集均为{NB,NM,NS,ZO,PS,PM,PB},其中NB,NM,NS,ZO,PS,PM,PB分别代表负大、负中、负小、零、正小、正中、正大;
根据被控对象GP,设置隶属度函数为:
&mu; F ( x ) = 0 , x < a x - a b - a , a &le; x < b c - x c - b , b &le; x < c 0 , x > c - - - ( 9 )
模糊规则:
比例增益系数Kp的模糊规则如表1:
表1
积分增益系数Ki的模糊规则如表2:
表2
使用模糊控制器模糊积分环节系数λ,积分环节系数λ的模糊规则如表
表3
隶属度函数的输入论域选择[-6,6],输出论域选择[-10,10]。
本发明的有益效果是:本发明一种自适应阶数分数阶模糊PIλ控制器方法,将模糊规则与分数阶PIλ控制器相结合。与传统的整数阶PI控制器、分数阶PIλ的控制器控制效果相比,自适应阶数分数阶模糊PIλ控制器,具有响应速度快,超调量小,对扰动具有较好的鲁棒性等优点。该方法可以有效的应用到含有非线性现象和多扰动输入的复杂控制系统中,提高系统的鲁棒性,达到更好的控制效果。
附图说明
图1是本发明自适应阶数分数阶模糊PIλ控制器方法中隶属度函数图;
图2是本发明自适应阶数分数阶模糊PIλ控制器方法中λ的阶跃响应图;
图3是本发明自适应阶数分数阶模糊PIλ控制器方法中λ模糊控制规则生成的控制面图;
图4是本发明自适应阶数分数阶模糊PIλ控制器方法与现有控制器的对比图。
具体实施方式
下面结合附图和具体实施方式对本发明进行详细说明。
本发明一种自适应阶数分数阶模糊PIλ控制器方法,具体按照以下步骤实施:
步骤1:针对被控对象GP,确定分数阶PIλ控制器的三个控制参数,具体为:
步骤1.1:根据被控对象GP需要,确定幅值裕量Am和相位裕量φm
步骤1.2:根据幅值裕量和相位裕量来设计分数阶PIλ控制器以满足系统的实际要求,即:
φm=arg[Gc(jωc)Gp(jωc)]+π (1)
A m = 1 G c ( j&omega; c ) G p ( j&omega; c ) - - - ( 2 )
d { arg &lsqb; G c ( j &omega; ) G p ( j &omega; ) &rsqb; } d &omega; | &omega; = &omega; c = 0 - - - ( 3 )
控制器传递函数Gc为:
G c = K p + K i s &lambda; - - - ( 4 )
将控制器传递函数Gc和被控对象传递函数GP改写为频域模式后代入公式(1)-公式(3)中,在设计控制器时,Gc、GP和期望的幅值裕量Am以及相位裕量φm都是已知,即可得到控制器的比例增益系数Kp、积分增益系数Ki及积分环节系数λ。
其中,为了实现拉普拉斯算子,我们采用Oustaloup算法在实际仿真中来近似分数阶次微积分,其频域的表达式如下:
s &alpha; = K &Pi; k = - N N s + &omega; k &prime; s + &omega; k - - - ( 5 )
其中,N为滤波器阶数,本发明中N统一取5,设定频域上下限为(ωhb),ωb·ωh=1,ωk、ωk′分别为极点和零点,K为增益系数:
K = ( &omega; h / &omega; b ) - &alpha; / 2 &Pi; k = - N N &omega; k &prime; &omega; k - - - ( 6 )
&omega; k = &omega; b ( &omega; h / &omega; b ) K + N + ( 1 + &alpha; ) / 2 ( 2 N + 1 ) - - - ( 7 )
&omega; k &prime; = &omega; b ( &omega; h / &omega; b ) K + N + ( 1 - &alpha; ) / 2 ( 2 N + 1 ) - - - ( 8 ) .
步骤2:利用模糊集合理论,设立模糊规则以及隶属度函数,将建立的模糊控制器加入到分数阶PIλ控制器,用于修正参数,具体为:
设误差e及误差变化率ec的模糊子集均为{NB,NM,NS,ZO,PS,PM,PB},其中NB、NM、NS、ZO、PS、PM、PB分别代表负大、负中、负小、零、正小、正中、正大;
根据被控对象GP,设置隶属度函数形式,本发明中选用的隶属函数为三角型,其形式如公式(9)所示,输入输出的模糊子集所对应的隶属度函数如图1所示(图1中所示为输入论域),输入论域选择[-6,6],输出论域选择[-10,10]:
&mu; F ( x ) = 0 , x < a x - a b - a , a &le; x < b c - x c - b , b &le; x < c 0 , x > c - - - ( 9 )
根据实际工程经验,以及被控对象对于干扰等复杂工况的要求得出Kp、Ki模糊规则如表1、表2,在模糊以上两个系数的基础上,本发明使用模糊控制器模糊积分环节系数λ,由自控原理知识可知,当λ升高会逐渐增大系统的超调量、增加振荡的时间以及拉长调节时间,但此时的静态误差相当小。如图2所示,当λ的值接近1时,系统收敛加快,稳定时间缩短,但当λ的值大于1时,振荡更为剧烈,稳定时间延长,甚至有发散的趋势。这说明λ的取值既不能太小,也不能过大,本发明制定的λ的模糊规则如表3:
比例增益系数Kp的模糊规则如表1:
表1
积分增益系数Ki的模糊规则如表2:
表2
使用模糊控制器模糊积分环节系数λ,积分环节系数λ的模糊规则如表3:
表3
在仿真软件中正确按照本发明中λ的模糊规则输入后,可以生成如图3所示的控制面图。
步骤3:将经过模糊控制器修正的分数阶控制器连入控制系统中,根据模糊控制器输入输出和系统输入输出对比后,决定比例因子Ke、Kec、Ku,使模糊控制器的输入输出与控制系统处于同一数量级:
Ke=n/xe,Kec=m/xec,Ku=y/l (10)
其中,xe、xec、y分别为误差、误差变化率、控制输出量的基本论域,n、m、l分别为误差、误差变化率、控制输出量的模糊集论域。
实施例
PWM整流器的电压外环传递函数为:
G p ( s ) = K 1 - T v s 1 + T i s
其中各参数的值为:
T i = 0.5 R L C = 0.495 T v = LP m 1.5 E m 2 = 0.000151 K = 3 R L E m 4 V d c = 9.06
按照步骤1叙述的方法得到的分数阶PIλ控制器如下式所示,按照参数搭建分数阶控制器模块。
G c ( s ) = 1.6 + 33.76 s 0.857
将整数阶PI控制器、传统分数阶PIλ控制器与本发明中的自适应阶数分数阶模糊PIλ控制器采用同一组参数(即由步骤1得出的参数,此时整数阶λ恒为1)进行对比仿真研究,对系统施加阶跃响应,如图4所示。
在同样一组参数下,自适应模糊分数阶PIλ控制器依靠模糊自整定参数的优势,在阶跃响应中具有较短的上升时间以及较快的稳定时间。VOFFPI(Variable OrderFractional Fuzzy PI,自适应阶数模糊分数阶PI)控制器超调1.665%,峰值时间为0.0143秒。FOPI(Fractional Order PI,分数阶PI)控制器超调1.981%,峰值时间为0.0156秒,VOFFPI各方面性能都优于分数阶PIλ控制器及OPI(Ordinary PI,整数阶PI)控制器。
系统稳定后,在0.25s时加入负载扰动,仿真结果数据如表4所示。
表4
在加入扰动后,三种控制器都能恢复稳定运行状态。其中VOFFPI在干扰波动后恢复快,鲁棒性强于其他两种控制器,分数阶控制器次之,整数阶控制器最差,对于干扰扰动过大,且恢复时间较慢。
本发明一种自适应阶数分数阶模糊PIλ控制器方法,将模糊规则与分数阶PIλ控制器相结合。与传统的整数阶PI控制器、分数阶PIλ的控制器控制效果相比,自适应阶数分数阶模糊PIλ控制器,具有响应速度快,超调量小,对扰动具有较好的鲁棒性等优点。该方法可以有效的应用到含有非线性现象和多扰动输入的复杂控制系统中,提高系统的鲁棒性,达到更好的控制效果。

Claims (5)

1.一种自适应阶数分数阶模糊PIλ控制器方法,其特征在于,具体按照以下步骤实施:
步骤1:针对被控对象GP,确定分数阶PIλ控制器的三个控制参数;
步骤2:利用模糊集合理论,设立模糊规则以及隶属度函数,将建立的模糊控制器加入到分数阶PIλ控制器,用于修正参数;
步骤3:将经过模糊控制器修正的分数阶控制器连入控制系统中,模糊控制器根据输入实时修正PIλ控制器参数,达到需要的控制效果。
2.根据权利要求1所述的一种自适应阶数分数阶模糊PIλ控制器方法,其特征在于,所述步骤1具体为:
步骤1.1:根据被控对象GP需要,确定幅值裕量Am和相位裕量φm
步骤1.2:根据幅值裕量和相位裕量来设计分数阶PIλ控制器以满足系统的实际要求,即:
φm=arg[Gc(jωc)Gp(jωc)]+π (1)
A m = 1 G c ( j&omega; c ) G p ( j&omega; c ) - - - ( 2 )
d { arg &lsqb; G c ( j &omega; ) G p ( j &omega; ) &rsqb; } d &omega; | &omega; = &omega; c = 0 - - - ( 3 )
控制器传递函数Gc为:
G c = K p + K i s &lambda; - - - ( 4 )
将控制器传递函数Gc和被控对象传递函数GP改写为频域模式后代入公式(1)-公式(3)中,得到控制器的比例增益系数Kp、积分增益系数Ki及积分环节系数λ。
3.根据权利要求2所述的一种自适应阶数分数阶模糊PIλ控制器方法,其特征在于,所述步骤1.2中sλ的频域表达式为:
s &alpha; = K &Pi; k = - N N s + &omega; k &prime; s + &omega; k - - - ( 5 )
其中,N为滤波器阶数,频域上限为ωb,下限为ωh,ωb·ωh=1,ωk、ω′k分别为极点和零点,K为增益系数:
K = ( &omega; h / &omega; b ) - &alpha; / 2 &Pi; k = - N N &omega; k &prime; &omega; k - - - ( 6 )
&omega; k = &omega; b ( &omega; h / &omega; b ) K + N + ( 1 + &alpha; ) / 2 ( 2 N + 1 ) - - - ( 7 )
&omega; k &prime; = &omega; b ( &omega; h / &omega; b ) K + N + ( 1 - &alpha; ) / 2 ( 2 N + 1 ) - - - ( 8 ) .
4.根据权利要求1所述的一种自适应阶数分数阶模糊PIλ控制器方法,其特征在于,所述步骤2具体为:
设误差e及误差变化率ec的模糊子集均为{NB,NM,NS,ZO,PS,PM,PB},其中NB,NM,NS,ZO,PS,PM,PB分别代表负大、负中、负小、零、正小、正中、正大;
根据被控对象GP,设置隶属度函数为:
&mu; F ( x ) = 0 , x < a x - a b - a , a &le; x < b c - x c - b , b &le; x < c 0 , x > c - - - ( 9 )
模糊规则:
比例增益系数Kp的模糊规则如表1:
表1
积分增益系数Ki的模糊规则如表2:
表2
使用模糊控制器模糊积分环节系数λ,积分环节系数λ的模糊规则如表3:
表3
5.根据权利要求4所述的一种自适应阶数分数阶模糊PIλ控制器方法,其特征在于,所述隶属度函数的输入论域选择[-6,6],输出论域选择[-10,10]。
CN201610633219.4A 2016-08-04 2016-08-04 一种自适应阶数分数阶模糊PIλ控制器方法 Pending CN106019949A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610633219.4A CN106019949A (zh) 2016-08-04 2016-08-04 一种自适应阶数分数阶模糊PIλ控制器方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610633219.4A CN106019949A (zh) 2016-08-04 2016-08-04 一种自适应阶数分数阶模糊PIλ控制器方法

Publications (1)

Publication Number Publication Date
CN106019949A true CN106019949A (zh) 2016-10-12

Family

ID=57133847

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610633219.4A Pending CN106019949A (zh) 2016-08-04 2016-08-04 一种自适应阶数分数阶模糊PIλ控制器方法

Country Status (1)

Country Link
CN (1) CN106019949A (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108988384A (zh) * 2018-07-24 2018-12-11 西安工业大学 基于分数阶pir的并网电流直流分量抑制方法
CN109039180A (zh) * 2018-08-02 2018-12-18 山东大学 双馈感应发电机并网过程的分数阶控制方法
CN109459928A (zh) * 2018-06-12 2019-03-12 陕西科技大学 模糊分数阶PIDμ控制器的DDS置换蒸煮温度控制方法
CN110118945A (zh) * 2019-04-22 2019-08-13 华南理工大学 一种能馈型直流电子负载系统及其自组织模糊控制方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104444647A (zh) * 2014-11-07 2015-03-25 沈阳工业大学 直线电梯的磁悬浮导向系统稳定运行的控制方法
CN104617850A (zh) * 2014-12-31 2015-05-13 广西科技大学 永磁同步电机的双闭环控制器及控制方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104444647A (zh) * 2014-11-07 2015-03-25 沈阳工业大学 直线电梯的磁悬浮导向系统稳定运行的控制方法
CN104617850A (zh) * 2014-12-31 2015-05-13 广西科技大学 永磁同步电机的双闭环控制器及控制方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
梁涛年: "《分数阶PID控制器及参数不确定分数阶系统稳定域分析》", 《中国博士学位论文全文数据库 信息科技辑》 *
王宏文 等: "《永磁直线同步电机的模糊自适应分数阶PIλ Dμ 控制》", 《机床与液压》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109459928A (zh) * 2018-06-12 2019-03-12 陕西科技大学 模糊分数阶PIDμ控制器的DDS置换蒸煮温度控制方法
CN109459928B (zh) * 2018-06-12 2021-10-08 陕西科技大学 模糊分数阶PIDμ控制器的DDS置换蒸煮温度控制方法
CN108988384A (zh) * 2018-07-24 2018-12-11 西安工业大学 基于分数阶pir的并网电流直流分量抑制方法
CN109039180A (zh) * 2018-08-02 2018-12-18 山东大学 双馈感应发电机并网过程的分数阶控制方法
CN109039180B (zh) * 2018-08-02 2020-01-17 山东大学 双馈感应发电机并网过程的分数阶控制方法
CN110118945A (zh) * 2019-04-22 2019-08-13 华南理工大学 一种能馈型直流电子负载系统及其自组织模糊控制方法
CN110118945B (zh) * 2019-04-22 2020-09-22 华南理工大学 一种能馈型直流电子负载系统及其自组织模糊控制方法

Similar Documents

Publication Publication Date Title
CN106019949A (zh) 一种自适应阶数分数阶模糊PIλ控制器方法
Li et al. Robust H∞ sliding mode observer design for a class of Takagi–Sugeno fuzzy descriptor systems with time-varying delay
Chen et al. Stochastic h/sub 2//h/sub/spl infin//control with state-dependent noise
Du et al. Recursive design of finite-time convergent observers for a class of time-varying nonlinear systems
Li et al. A heuristic approach to static output-feedback controller synthesis with restricted frequency-domain specifications
Wang et al. Direct adaptive neural control for stabilization of nonlinear time-delay systems
Doye et al. Robust fractional-order proportional-integral observer for synchronization of chaotic fractional-order systems
Lam et al. and Model Reduction for System Input with Sector Nonlinearities
Verma et al. Intelligent automatic generation control of two-area hydrothermal power system using ANN and fuzzy logic
Cui et al. Adaptive neural tracking control of nonlinear stochastic switched non-lower triangular systems with input saturation
CN108599649A (zh) Pmsm位置伺服系统高阶对象控制器设计及参数确定方法
CN106774379A (zh) 一种智能超螺旋强鲁棒姿态控制方法
Kolmanovksy et al. Discontinuous feedback stabilization of nonholonomic systems in extended power form
CN105892297A (zh) 一种自适应分数阶动态滑模控制算法
Kang et al. Finite-time command filter-based adaptive fuzzy tracking control for stochastic nonlinear induction motors systems with unknown backlash-like hysteresis
Trofino et al. Switching rule design for affine switched systems with ℋ∞ performance
CN103558761B (zh) 一种具有控制器输入饱和的非线性化学反应循环不确定时滞系统的控制方法
Yang et al. Decentralized adaptive funnel control of uncertain large-scale interconnected nonlinear system
CN107276474A (zh) 一种基于模糊滑模的同步发电机励磁控制方法
CN107102555A (zh) 一种镇定一阶惯性加纯滞后系统的线性自抗扰控制器设计方法
Lagrat et al. Fuzzy sliding mode PI controller for nonlinear systems
Sereni et al. Relaxed LMI conditions for the design of robust static output feedback controllers
Yanev et al. Multivariable system's parameters interaction and robust control design
CN109388065A (zh) 一种包含闭环参考模型的干扰观测器设计方法
Ashraf-Modarres et al. A new decentralized sliding mode observer design for interconnected power systems with unknown and time-varying delay

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication

Application publication date: 20161012

RJ01 Rejection of invention patent application after publication