CN107169557A - 一种对布谷鸟优化算法进行改进的方法 - Google Patents
一种对布谷鸟优化算法进行改进的方法 Download PDFInfo
- Publication number
- CN107169557A CN107169557A CN201710341246.9A CN201710341246A CN107169557A CN 107169557 A CN107169557 A CN 107169557A CN 201710341246 A CN201710341246 A CN 201710341246A CN 107169557 A CN107169557 A CN 107169557A
- Authority
- CN
- China
- Prior art keywords
- cuckoo
- optimized algorithm
- algorithm
- bird
- nest
- 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
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/004—Artificial life, i.e. computing arrangements simulating life
- G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
Abstract
本发明公开了一种对布谷鸟优化算法进行改进的方法,有效解决传统布谷鸟优化算法收敛精度低、迭代后期易陷入局部最优的缺陷。首先,通过动态自适应调整步长a和发现概率Pa,实现算法精细化搜索过程;其次,引入反向学习策略,增加种群多样性,提升算法迭代运行效率;最后,根据预先设置的停滞次数判定依据,启用多开端策略,跳出局部最优,进而得到最优解。根据本发明建立上述改进型布谷鸟优化算法的实施例,仿真结果表明该方法在收敛速度、收敛精度、收敛可靠性均得到一定程度的改善。
Description
技术领域
本发明涉及智能优化算法技术领域,具体是一种对布谷鸟优化算法进行改进的方法。
背景技术
在自然界中,布谷鸟寻找适合自己产卵的鸟巢位置是随机的或类似随机,为了模拟布谷鸟寻窝的方式,首先,需要假定以下3个理想的状态:
(1)每只布谷鸟一次只产一个卵,并随机选择鸟巢位置来孵化它;
(2)在随机选择的一组鸟巢中,最好的鸟巢将会被保留到下一代;
(3)所使用的宿主鸟巢的数量是固定的,而布谷鸟的鸟卵有一定的概率Pa∈[0,1]被宿主鸟发现,在这种情况下,宿主鸟就会把布谷鸟的鸟卵扔掉或者抛弃自己的巢到另外的地方重新建巢。在第三个规则当中,我们认为n个鸟巢中Pa部分由新的鸟巢(具有新的随机解)所代替。
对于一个最大值问题,一个解的适应度或质量是和它的目标函数值成正比的,这类似于其它只能算法如遗传算法,在具体到算法中,每一个巢的位置表示一个解,当要产生一个新解x(t+1)时,对于每个布谷鸟都要实施一次莱维飞行,其目的是要用新的和可能更优的解去代替不太好的解,每个布谷鸟寻巢的路径和位置更新公式如下:
式中,a>0是步长,它与所要解决问题的尺度有关,在算法中通常取a=1;表示点乘积;Levy(λ)为Lévy随机搜索路径。
一般来讲,布谷鸟搜索算法的随机运动就是马尔可夫链,它的下一个状态或位置只取决于当前位置(公式第一项)和转换概率(第二项)。表示点对点乘法,类似的运算在粒子群算法中也可以看到,但是这种通过莱维飞行产生的随机运动会更长,这使得它在探索解空间时会更加有效。从公式可以看到由于莱维飞行的随机运动,一些新解会产生在局部最优值附近,因此莱维飞行短步长加速了局部搜索。另外,由于莱维飞行产生的长步长,相当一部分新解会产生在距离局部最优值较远的地方,这就确保算法不会陷入局部最小值。
在Mantegna的算法中,基于莱维飞行的步长计算公式为:
式中,u、v服从正态分布;Γ为标准的Gamma函数,该概率分布的方差与均值都是无界的。
因此,布谷鸟寻找鸟巢的运动公式为:
现有的布谷鸟优化算法中种群多样性缺失,易陷入局部最优,不利于工程化应用。
发明内容
本发明的目的在于提供一种有效解决传统布谷鸟优化算法收敛精度低、迭代后期易陷入局部最优的缺陷的对布谷鸟优化算法进行改进的方法,以解决上述背景技术中提出的问题。
为实现上述目的,本发明提供如下技术方案:
一种对布谷鸟优化算法进行改进的方法,包括如下步骤:
步骤10:动态自适应参数a和Pa的获取;
步骤20:反向学习策略的引入;
步骤30:多开端策略的启用。
作为本发明进一步的方案:所述步骤10中参数a和Pa的获取过程如下:
式中,astart、aend分别表示a的初值和终值;Pastart、Paend分别表示Pa的初值和终值;t为当前迭代次数,Maxgen为最大迭代次数。
作为本发明进一步的方案:所述步骤20中反向学习策略的引入过程如下:
a)初始阶段:随机产生n个鸟巢的初始位置并依据反向点原理计算其对应的反向点找出当前最优鸟巢的位置
b)迭代阶段:根据随机数A值的大小,判断是否需要启动反向学习策略,根据适应度值的计算结果淘汰适应度值较差的个体。
作为本发明再进一步的方案:所述步骤30中多开端策略的启用过程如下:
根据最优解停止更新的次数S来判断是否需要初始化种群,如果S大于允许停滞次数Smax,则开始初始化种群,并将S重新置零,并保存个体最优和全局最优的信息。
与现有技术相比,本发明的有益效果是:
本发明建立了一种改进型布谷鸟优化算法,与现有布谷鸟优化算法相比,本发明方法充分考虑了算法收敛性和计算复杂性,设定对应的门槛条件来分别启用反向学习策略及多开端策略,具有更快的收敛速度,更高的收敛精度,在对复杂非线性问题进行求解时,该方法既能够有效提高计算精度,同时又能够加快计算速度。
附图说明
图1为本发明的流程示意图。
图2为本发明基于不同A值时的迭代过程示意图。
图3为本发明基于不同Smax值时的迭代过程示意图。
图4为本发明实施例中Schwefel函数的迭代过程示意图。
图5为本发明实施例中Sphere函数函数的迭代过程示意图。
具体实施方式
下面结合具体实施方式对本专利的技术方案作进一步详细地说明。
请参阅图1-5,一种对布谷鸟优化算法进行改进的方法,包括如下步骤:
步骤10:动态自适应参数a和Pa的获取;
步骤20:反向学习策略的引入;
步骤30:多开端策略的启用。
所述步骤10中参数a和Pa的获取过程如下:
式中,astart、aend分别表示a的初值和终值;Pastart、Paend分别表示Pa的初值和终值;t为当前迭代次数,Maxgen为最大迭代次数。
所述步骤20中反向学习策略的引入过程如下:
a)初始阶段:随机产生n个鸟巢的初始位置并依据反向点原理计算其对应的反向点找出当前最优鸟巢的位置
b)迭代阶段:根据随机数A值的大小,判断是否需要启动反向学习策略,根据适应度值的计算结果淘汰适应度值较差的个体。
所述步骤30中多开端策略的启用过程如下:
根据最优解停止更新的次数S来判断是否需要初始化种群,如果S大于允许停滞次数Smax,则开始初始化种群,并将S重新置零,并保存个体最优和全局最优的信息。
所述对布谷鸟优化算法进行改进的方法,具体包括如下步骤:
(1)对参数进行初始化操作,包括鸟巢的数量n、外来鸟卵发现概率范围[Pastart,Paend]、步长范围[astart,aend]、迭代次数Maxgen、当前迭代次数t和多开端策略允许迭代次数Smax等参数;
(2)求取种群中各布谷鸟个体的适应度值。随机产生n个鸟巢的初始位置并依据定义1计算其对应的反向点找出当前最优鸟巢的位置
(3)根据当前迭代次数t,对CSA中参数进行自适应调整;
(4)保留上代最优鸟巢的位置并利用Lévy-flights策略对鸟巢位置进行更新,得到一组新的鸟巢位置对这组新的鸟巢位置进行评估,与上代产生的鸟巢位置进行比较,让适应度值较好的鸟巢位置替换适应度值较差的鸟巢位置从而得到一组更优的鸟巢位置;
(5)通过概率Pa保留一部分当前代中不容易被发现的鸟巢位置,同时随机改变容易被发现的鸟巢位置,得到另一组新的鸟巢;与未改变之前的鸟巢位置进行比较,让适应度值较好的鸟巢位置替换适应度值较差的鸟巢位置从而得到下一组更优的鸟巢位置,对当前组鸟巢进行评估,找到当前组最优鸟巢位置;
(6)根据随机数A值的大小,判断是否需要启动反向学习策略,淘汰适应度值较差的个体;
(7)如果全局最优解满足更新判据,则执行t=t+1;否则,执行S=S+1,并根据是否达到停滞次数Smax,进而启动多开端策略,循环执行Step3~Step6;
(8)判断终止条件,若满足终止条件,则输出最优值,算法结束;否则,继续迭代执行Step3~Step7,直至算法终止条件被满足。
算法收敛性是衡量智能优化算法的一项重要指标。设智能优化算法为Γ,该算法第w次迭代结果为xw,第w+1次迭代结果为xw+1=Γ(xw,ξ),其中,ξ是算法Γ在整个迭代中已寻优过的解。
条件1:函数f(X)是搜索空间S内的连续函数,若f(Γ(xw,ξ))≤f(x),且ξ∈S,则有f(Γ(x,ξ))≤f(x)。
条件1是智能优化算法收敛的基础,可以保证算法在迭代中总可以产生优于当前个体的新个体。
条件2:对于任意B∈S,s.t.v[B]>0,则有
由条件2可知,当B在搜索空间S中满足v[B]>0时,智能优化算法Γ经过连续无穷次的搜索总能够寻找到B中的点。
定理1(改进型布谷鸟优化算法的收敛性):改进型布谷鸟优化算法能够以概率1收敛到全局最优解。
证明:在改进型布谷鸟优化算法中,布谷鸟鸟巢个体极值和全局极值根据式(7)进行更新。
以最小化的优化问题为例,布谷鸟优化算法在迭代过程中更新极值序列均为单调递减;同时,多开端策略是一种直接搜索方法,通过对当前全局最优极值点保留并优化能保证解进一步单调下降,即保证新解优于旧解时才对其进行替换,故改进型布谷鸟优化算法能够满足条件1。
当布谷鸟鸟巢全部收敛到全局极值点位置时,布谷鸟优化算法处于停滞状态。此时,全局极值点可能仅仅是局部最优点,因算法早熟收敛所引起。由于改进型布谷鸟优化算法中根据停滞次数Smax引入多开端策略,使原本聚集在局部最优点附近的鸟巢能够重新随机分布在整个搜索空间中。因此,当迭代次数趋于无穷大时,该算法总能够搜索到全局最优解,故改进型布谷鸟优化算法也能够满足条件2。
由上述证明可知,当迭代次数足够大时,改进型布谷鸟优化算法总能够以概率1收敛到全局最优解。
原理:本发明通过将参数动态调整、反向学习策略及多开端策略引入到布谷鸟优化算法中,其目的旨在有效解决传统布谷鸟优化算法种群多样性少、收敛精度低、迭代后期易陷入局部最优的缺陷。
为了验证上述改进型布谷鸟优化算法的正确性和有效性,本发明选取3个标准测试函数进行仿真验证,测试函数分别如表1所示,改进布谷鸟优化算法的参数如表2所示。
表1 标准测试函数
表2 改进型布谷鸟优化算法中参数设置
不同随机数A值和停滞次数Smax的选取,对于改进算法的迭代次数具有较大的影响。当A值由0.1~0.9之间变化时,采用C_OBL(布谷鸟优化算法+反向学习策略)对Rastrigin函数进行迭代寻优,迭代过程如图2所示;固定A=0.5,当Smax由50~250之间变化时,采用C_MOBL(布谷鸟优化算法+反向学习策略+多开端策略)对Rastrigin函数进行迭代寻优,迭代过程如图3所示。
由图2可知,当随机数A值选取过大或者过小时,迭代次数均较大,因此,本发明专利将A值设定为0.5。由图3可知,随着Smax的不断增大,迭代次数也不断递增,使得多开端策略逐渐失效,因此,本发明专利将Smax设定为50。
在固定进化次数条件下,分别采用Cuckoo(布谷鸟优化算法)、C_Pars(布谷鸟优化算法+参数动态自适应)、C_OBL、C_MOBL和C_MOBLPars(布谷鸟优化算法+反向学习策略+多开端策略+参数动态自适应)对函数f2~f3(M=10)进行迭代实验比较,结果如图4和图5所示。由图4和图5可知,对于上述2种不同的函数,本发明专利所提的C_MOBLPars均具有较高的收敛速度和收敛精度,说明反向学习策略、参数自适应调整和多开端策略的引入,更容易实现布谷鸟优化算法的全局收敛。
本发明建立了一种改进型布谷鸟优化算法,与现有布谷鸟优化算法相比,本发明方法充分考虑了算法收敛性和计算复杂性,设定对应的门槛条件来分别启用反向学习策略及多开端策略,具有更快的收敛速度,更高的收敛精度,在对复杂非线性问题进行求解时,该方法既能够有效提高计算精度,同时又能够加快计算速度。
上面对本专利的较佳实施方式作了详细说明,但是本专利并不限于上述实施方式,在本领域的普通技术人员所具备的知识范围内,还可以在不脱离本专利宗旨的前提下作出各种变化。
Claims (4)
1.一种对布谷鸟优化算法进行改进的方法,其特征在于,包括如下步骤:
步骤10:动态自适应参数a和Pa的获取;
步骤20:反向学习策略的引入;
步骤30:多开端策略的启用。
2.根据权利要求1所述的对布谷鸟优化算法进行改进的方法,其特征在于,所述步骤10中参数a和Pa的获取过程如下:
式中,astart、aend分别表示a的初值和终值;Pastart、Paend分别表示Pa的初值和终值;t为当前迭代次数,Maxgen为最大迭代次数。
3.根据权利要求1所述的对布谷鸟优化算法进行改进的方法,其特征在于,所述步骤20中反向学习策略的引入过程如下:
a)初始阶段:随机产生n个鸟巢的初始位置并依据反向点原理计算其对应的反向点找出当前最优鸟巢的位置
b)迭代阶段:根据随机数A值的大小,判断是否需要启动反向学习策略,根据适应度值的计算结果淘汰适应度值较差的个体。
4.根据权利要求1所述的对布谷鸟优化算法进行改进的方法,其特征在于,所述步骤30中多开端策略的启用过程如下:
根据最优解停止更新的次数S来判断是否需要初始化种群,如果S大于允许停滞次数Smax,则开始初始化种群,并将S重新置零,并保存个体最优和全局最优的信息。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710341246.9A CN107169557A (zh) | 2017-05-12 | 2017-05-12 | 一种对布谷鸟优化算法进行改进的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710341246.9A CN107169557A (zh) | 2017-05-12 | 2017-05-12 | 一种对布谷鸟优化算法进行改进的方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107169557A true CN107169557A (zh) | 2017-09-15 |
Family
ID=59815358
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710341246.9A Pending CN107169557A (zh) | 2017-05-12 | 2017-05-12 | 一种对布谷鸟优化算法进行改进的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107169557A (zh) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107747930A (zh) * | 2017-09-25 | 2018-03-02 | 华侨大学 | 一种基于万有引力加速布谷鸟算法的圆度误差评定方法 |
CN108387820A (zh) * | 2018-03-20 | 2018-08-10 | 东北电力大学 | 含分布式电源的配电网故障区段定位方法 |
CN108804384A (zh) * | 2018-06-21 | 2018-11-13 | 哈尔滨工业大学 | 一种用于最小区域球度评定的最优导向自适应搜索方法 |
CN108804390A (zh) * | 2018-06-21 | 2018-11-13 | 哈尔滨工业大学 | 一种基于改进布谷鸟搜索策略的最小区域球度评定方法 |
CN109115161A (zh) * | 2018-06-21 | 2019-01-01 | 哈尔滨工业大学 | 一种基于空间定向收缩及改进布谷鸟搜索算法的球度评定方法 |
CN109323677A (zh) * | 2018-08-21 | 2019-02-12 | 上海隧道工程有限公司 | 改进布谷鸟搜索算法的圆度误差评定算法 |
CN110555548A (zh) * | 2019-08-05 | 2019-12-10 | 三峡大学 | 基于数据挖掘原始误差校正的ics-elm超短期风电功率预测方法 |
CN111797899A (zh) * | 2020-06-04 | 2020-10-20 | 国网江西省电力有限公司电力科学研究院 | 一种低压台区kmeans聚类方法及系统 |
CN112100824A (zh) * | 2020-08-26 | 2020-12-18 | 西安工程大学 | 一种改进的布谷鸟算法及优化机器人结构参数的方法 |
CN113032921A (zh) * | 2021-03-16 | 2021-06-25 | 山东科技大学 | 基于并行自适应参数布谷鸟搜索与最低水平线的排样算法 |
CN115018221A (zh) * | 2022-08-10 | 2022-09-06 | 浙江浩普智能科技有限公司 | 基于改进布谷鸟搜索算法的锅炉负荷分配方法及系统 |
CN116051591A (zh) * | 2023-03-29 | 2023-05-02 | 长春工业大学 | 一种基于改进布谷鸟搜索算法的带钢图像阈值分割方法 |
-
2017
- 2017-05-12 CN CN201710341246.9A patent/CN107169557A/zh active Pending
Cited By (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107747930A (zh) * | 2017-09-25 | 2018-03-02 | 华侨大学 | 一种基于万有引力加速布谷鸟算法的圆度误差评定方法 |
CN107747930B (zh) * | 2017-09-25 | 2019-12-31 | 华侨大学 | 一种基于万有引力加速布谷鸟算法的圆度误差评定方法 |
CN108387820A (zh) * | 2018-03-20 | 2018-08-10 | 东北电力大学 | 含分布式电源的配电网故障区段定位方法 |
CN108804390B (zh) * | 2018-06-21 | 2020-02-07 | 哈尔滨工业大学 | 一种基于改进布谷鸟搜索策略的最小区域球度评定方法 |
CN108804384A (zh) * | 2018-06-21 | 2018-11-13 | 哈尔滨工业大学 | 一种用于最小区域球度评定的最优导向自适应搜索方法 |
CN108804390A (zh) * | 2018-06-21 | 2018-11-13 | 哈尔滨工业大学 | 一种基于改进布谷鸟搜索策略的最小区域球度评定方法 |
CN109115161A (zh) * | 2018-06-21 | 2019-01-01 | 哈尔滨工业大学 | 一种基于空间定向收缩及改进布谷鸟搜索算法的球度评定方法 |
CN109115161B (zh) * | 2018-06-21 | 2020-02-07 | 哈尔滨工业大学 | 一种基于空间定向收缩及改进布谷鸟搜索算法的球度评定方法 |
CN109323677B (zh) * | 2018-08-21 | 2020-08-11 | 上海隧道工程有限公司 | 改进布谷鸟搜索算法的圆度误差评定算法 |
CN109323677A (zh) * | 2018-08-21 | 2019-02-12 | 上海隧道工程有限公司 | 改进布谷鸟搜索算法的圆度误差评定算法 |
CN110555548A (zh) * | 2019-08-05 | 2019-12-10 | 三峡大学 | 基于数据挖掘原始误差校正的ics-elm超短期风电功率预测方法 |
CN111797899A (zh) * | 2020-06-04 | 2020-10-20 | 国网江西省电力有限公司电力科学研究院 | 一种低压台区kmeans聚类方法及系统 |
CN111797899B (zh) * | 2020-06-04 | 2023-11-07 | 国网江西省电力有限公司电力科学研究院 | 一种低压台区kmeans聚类方法及系统 |
CN112100824A (zh) * | 2020-08-26 | 2020-12-18 | 西安工程大学 | 一种改进的布谷鸟算法及优化机器人结构参数的方法 |
CN112100824B (zh) * | 2020-08-26 | 2024-02-27 | 西安工程大学 | 一种改进的布谷鸟算法及优化机器人结构参数的方法 |
CN113032921A (zh) * | 2021-03-16 | 2021-06-25 | 山东科技大学 | 基于并行自适应参数布谷鸟搜索与最低水平线的排样算法 |
CN113032921B (zh) * | 2021-03-16 | 2022-12-13 | 山东科技大学 | 基于并行自适应参数布谷鸟搜索与最低水平线的排样算法 |
CN115018221A (zh) * | 2022-08-10 | 2022-09-06 | 浙江浩普智能科技有限公司 | 基于改进布谷鸟搜索算法的锅炉负荷分配方法及系统 |
CN115018221B (zh) * | 2022-08-10 | 2022-11-11 | 浙江浩普智能科技有限公司 | 基于改进布谷鸟搜索算法的锅炉负荷分配方法及系统 |
CN116051591A (zh) * | 2023-03-29 | 2023-05-02 | 长春工业大学 | 一种基于改进布谷鸟搜索算法的带钢图像阈值分割方法 |
CN116051591B (zh) * | 2023-03-29 | 2023-06-16 | 长春工业大学 | 一种基于改进布谷鸟搜索算法的带钢图像阈值分割方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107169557A (zh) | 一种对布谷鸟优化算法进行改进的方法 | |
CN109116841B (zh) | 一种基于蚁群算法的路径规划平滑优化方法 | |
CN111024092B (zh) | 一种多约束条件下智能飞行器航迹快速规划方法 | |
CN108828436B (zh) | 基于混沌云自适应萤火虫算法的模拟电路故障诊断方法 | |
CN110262548A (zh) | 一种考虑抵达时间约束的无人机航迹规划方法 | |
CN114115362B (zh) | 一种基于双向apf-rrt*算法的无人机航迹规划方法 | |
CN105893694A (zh) | 一种基于重采样粒子群优化算法的复杂系统设计方法 | |
CN109931943B (zh) | 无人船舶全局路径规划方法及电子设备 | |
CN114460941B (zh) | 一种基于改进麻雀搜索算法的机器人路径规划方法及系统 | |
CN113296496B (zh) | 基于多采样点的引力自适应步长双向rrt路径规划方法 | |
CN107633105B (zh) | 一种基于改进混合蛙跳算法的四旋翼无人机参数辨识方法 | |
CN113296520A (zh) | 融合a*与改进灰狼算法的巡检机器人路径规划方法 | |
CN109800849A (zh) | 动态布谷鸟搜索算法 | |
CN115113628A (zh) | 一种基于改进灰狼算法的巡检机器人路径规划方法 | |
CN112465160A (zh) | 一种基于vr的车辆维修保养辅助系统 | |
CN115525038A (zh) | 一种基于联邦分层优化学习的设备故障诊断方法 | |
CN109472105A (zh) | 半导体产品良率上限分析方法 | |
CN112329327A (zh) | 一种硬件感知的液体状态机网络生成方法及系统 | |
CN111062462A (zh) | 基于差分进化算法的局部搜索和全局搜索融合方法及系统 | |
CN115019510A (zh) | 一种基于动态自适应生成对抗网络的交通数据修复方法 | |
CN108615069A (zh) | 一种基于改进的具有适应性的量子粒子群算法的优化计算方法 | |
CN111381605A (zh) | 一种应用于多无人机大范围海域水下多目标协同搜索方法 | |
CN109961129A (zh) | 一种基于改进粒子群的海上静止目标搜寻方案生成方法 | |
CN106650028A (zh) | 一种基于敏捷卫星设计参数的优化方法及系统 | |
CN105678844A (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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170915 |
|
RJ01 | Rejection of invention patent application after publication |