CN104267601A - 二重随机跳变系统基于观测器的有限短时间控制方法 - Google Patents

二重随机跳变系统基于观测器的有限短时间控制方法 Download PDF

Info

Publication number
CN104267601A
CN104267601A CN201410488294.7A CN201410488294A CN104267601A CN 104267601 A CN104267601 A CN 104267601A CN 201410488294 A CN201410488294 A CN 201410488294A CN 104267601 A CN104267601 A CN 104267601A
Authority
CN
China
Prior art keywords
observer
short time
double random
limited short
random jump
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
CN201410488294.7A
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.)
Jiangnan University
Original Assignee
Jiangnan University
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 Jiangnan University filed Critical Jiangnan University
Priority to CN201410488294.7A priority Critical patent/CN104267601A/zh
Publication of CN104267601A publication Critical patent/CN104267601A/zh
Pending legal-status Critical Current

Links

Landscapes

  • Feedback Control In General (AREA)

Abstract

二重随机跳变系统基于观测器的有限时间控制方法,涉及多重随机跳变系统的多重概率描述以及随机系统在有限短时间内的暂态性能控制。首先,用高斯概率密度函数来描述二重随机跳变系统的概率随机分布特性;其次,通过放宽对系统李雅普诺夫能量函数在采样时刻严格递减的条件,定义新的二重随机跳变过程在各模态下有限短时间稳定性定义;再次,针对系统的能量有界干扰设计基于观测器的控制器,使得闭环系统的过程轨迹在平衡点的一定范围内受限运动,从而放弃对渐进稳定的要求。本发明针对生产实际要求的短时间工作系统,考虑二重随机跳变现象以及系统状态不可测的实际情形,设计基于观测器的有限时间控制方法。

Description

二重随机跳变系统基于观测器的有限短时间控制方法
技术领域
本发明涉及一种随机跳变系统基于观测器的有限时间控制方法,特别是针对二重随机跳变过程基于观测器的有限短时间控制方法,该方法可用于生化系统、网络系统、机器人系统、通信系统、经济系统、航空航天等领域。
背景技术
很多实际系统,人们更感兴趣的常常是其能否在有限短时间内满足暂态要求。为此,Dorato于1961年提出了有限短时间稳定的概念,进而分析了系统的有限时间控制问题。
另一方面很多实际过程,系统的状态常常不可测量,基于观测器的控制方法得到了学者的广泛关注。但二重随机跳变系统基于观测器的有限短时间控制问题仍未解决。本发明针对实际工程过程中存在的短时间工作系统以及二重随机跳变过程,考虑状态不可测情形,提供一种基于观测器的有限短时间控制方法,使得系统状态轨迹在平衡点的一定范围内受限运动。
发明内容
本发明要解决的技术问题是:针对实际工程过程中存在的短时间工作系统以及二重随机跳变过程,用高斯过程来描述跳变系统的二重随机跳变特性,考虑外部干扰能量有界,提供一种基于观测器的有限短时间控制方法,使得系统状态轨迹在平衡点的一定范围内受限运动,从时间角度为降低一般渐进稳定的工程保守性提供了思路。
本发明的技术解决方案为:首先,针对二重随机跳变系统,用高斯随机分布来描述跳变系统的二重随机跳变特性;其次,从能量角度,通过允许系统的李雅普诺夫能量函数在采样时刻递增,定义跳变系统各模态下新的有限短时间稳定性定义;再次,假设外部干扰能量有界,设计基于观测器的控制器,使得二重跳变下的各模态状态轨迹在平衡点的一定范围内受限运动并具有干扰抑制能力,具体步骤如下:
(1)二重随机跳变过程描述:
a.对二重随机跳变过程构造跳变系统模型;
b.用高斯随机分布来描述随机跳变过程的二重跳变特性;
c.建立二重随机跳变过程转移概率密度函数矩阵;
(2)二重随机特性下各子系统有限短时间稳定性定义:
a.给出各模态下的初始状态的受限空间;
b.根据初始空间,结合轨迹要求,通过新的李雅普诺夫能量函数条件,定义二重随机跳变过程在各模态下有限短时间稳定性定义;
(3)基于观测器的有限短时间控制器设计:
a.针对步骤(1)中构造的系统,设计全阶状态观测器;
b.设计基于观测器的控制器,并将其代入步骤(1)中构造的系统,得到闭环误差动态系统;
c.选取李雅普诺夫函数,并定义有限短时间抗干扰性能指标;
d.基于李雅普诺夫稳定性定理及有限短时间抗干扰控制方法,利用已经获得的期望跳变转移概率矩阵,设计基于观测器的控制器,获取使闭环系统有限时间稳定并满足干扰抑制性能的控制器存在条件;本发明针对实际工程应用中普遍存在的二重随机跳变系统以及短时间工作系统,首次设计基于观测器的控制器,使得闭环跳变系统有限短时间稳定的同时还要具有干扰抑制能力,与现有技术相比的优点在于:
1.本发明用高斯概率密度函数对二重随机跳变过程进行描述,具有实际意义。
2.本发明通过允许李雅普诺夫能量函数在采样时刻递增,定义了新的二重随机跳变过程确保各子系统有限短时间稳定的定义。
3.本发明利用线性矩阵不等式技术设计基于观测器的控制器,一方面不仅计算简单,便捷可行,另一方面,不仅能够使系统的状态轨迹在平衡点的一定范围内受限运动,而且能针对所有频段的外部干扰信号具有干扰抑制能力。
4.本发明设计的基于观测器的控制器还可以应用到不确定系统、时滞系统等复杂工业过程,具有普适性。
附图说明
图1二重随机跳变过程模态图
图2系统状态轨迹图
图3系统观测误差图
具体实施方式
下面结合附图所示实施例,对本发明作进一步详细描述。
需要强调的是,本发明涉及的技术并不仅适用于下面提及的例子,这些技术可以被用于任何适用的随机跳变控制系统。
本发明基于观测器的二重随机跳变系统有限短时间控制方法,包括以下步骤:
(1)二重随机跳变过程描述
(2)二重随机特性下各子系统有限短时间稳定性定义
(3)基于观测器的有限短时间控制器设计
(4)仿真实验验证
下面介绍具体步骤:
(1)二重随机跳变过程描述
考虑如下一类离散跳变系统:
x(k+1)=A(rk)x(k)+B(rk)u(k)+Bw(rk)w(k)
z(k)=C(rk)x(k)+D(rk)u(k)+Dw(rk)w(k)
y(k)=E(rk)x(k)
x(k)=x0,rk=r0,k=0
其中,x(k)∈Rn是系统的状态向量;u(k)∈Rm是系统的控制向量;是外部扰动信号;z(k)∈Rl是系统的被控输出;A(rk),B(rk),Bw(rk),C(rk),D(rk),Dw(rk)和E(rk)分别为已知的与模态rk相关的适当维数的系数矩阵,其中rk表示系统的模态,为在有限集合Μ={1,2,…,s}中随时间k取值的二重随机跳变过程,其跳变转移概率定义如下:
π ij ( ξ k ) = P r ( r k = j | r k - 1 = i , k )
式中表示从模态i跳变到模态j的转移概率。为了方便起见,当rk=i时,分别用Ai,Bi,Bwi,Ci,Di,Dwi及Ei表征A(rk),B(rk),Bw(rk),C(rk),D(rk),Dw(rk)和E(rk)。
不失一般性,用高斯随机分布{ξk,k∈K}来描述二重随机跳变过程的连续时变特性,其受限高斯概率密度函数表征为:
p ( π ij ( ξ k ) ) = 1 σ ij f ( π ij ( ξ k ) - μ ij σ ij ) F ( 1 - μ ij σ ij ) - F ( 0 - μ ij σ ij )
式中f(·)为高斯概率密度函数的标准分布;F(·)为f(·)的累积分布函数,μij和σij为转移概率矩阵中各元素的高斯概率密度函数的均值和方差信息。基于上述描述,转移概率密度函数矩阵可表达为下式:
其中 n ( μ ij , σ ij ) = p ( π ij ( ξ k ) ) 的受限高斯概率密度函数。
(2)二重随机特性下各子系统有限短时间稳定性定义:
定义:设u(k)=0以及w(k)=0,被控二重随机跳变系统的各子系统是关于(c1 c2 N Ri)有限时间稳定的,其中c1是初始空间,c2是受限空间,且满足c1<c2,Ri>0,N为要求的时间常数,如果下列条件成立:
x 0 T R i x 0 &le; c 1 &DoubleRightArrow; x k T R i x k < c 2 , &ForAll; k &Element; { 1,2 , . . . , N }
对于上述定义,如果考虑系统受到外部干扰的影响,并假设干扰信号能量有界,则对于u(k)=0,被控系统的各子系统是关于(c1 c2 N Ri d)有限时间有界的,其中d为未知输入信号的上界,如果对于满足能量有界的干扰信号,均有上述条件成立。
同样,如果对于上述定义,使用状态反馈控制,则被控系统的各子系统是关于(c1 c2 N Ri d)有限时间可镇定的,如果对于满足能量有界的干扰信号,均有上述条件成立。
上述有限时间稳定的定义与Lyapunov意义下的渐进稳定是两个不同的概念,两者并无直接的关联,系统Lyapunov意义下的稳定并不能确保有限时间稳定;同样,有限时间稳定也不能保证Lyapunov渐进稳定。
(3)基于观测器的有限短时间控制器设计:
a.针对步骤(1)中构造的系统,设计全阶状态观测器;
x &OverBar; k + 1 = A i x &OverBar; k + B i u k + H i , &xi; k ( y k - y &OverBar; k ) y &OverBar; k = E i x &OverBar; k u k = K i , &xi; k x &OverBar; k
其中为观测器状态,为观测器估计输出,分别为待求的观测器和控制器参数。
b.定义观测误差为: e k = x k - x &OverBar; k , 状态变量为: x ~ k = x k T e k T T
则可得如下闭环误差动态系统:
x ~ ( k + 1 ) = A ~ i x ~ ( k ) + B ~ wi w ( k ) z ( k ) = C ~ i x ~ ( k ) + D wi w ( k )
其中 A ~ i = A i + B i K i , &xi; k - B i K i , &xi; k 0 A i - E i H i , &xi; k , B ~ wi = B wi B wi , C ~ i = C i + D i K i , &xi; k - D i K i , &xi; k .
c.选取李雅普诺夫泛函 V ( x ~ k , r k , &xi; k ) = x ~ k T P ~ i , &xi; k x ~ k , 以及抗干扰性能指标 J = &Delta; E { &Sigma; k = 0 &infin; [ z k T z k - &gamma; 2 w k T w k ] } .
d.基于李雅普诺夫稳定性定理及有限短时间H控制方法,结合Markov跳变理论,利用已经获得的期望跳变转移概率矩阵,获取使闭环系统有限短时间稳定并满足H性能的充分条件;
结论:如果对于给定标量α≥0及γ>0,闭环控制系统有限时间有界,且对所有满足能量有界条件的外部干扰 w ( k ) &Element; l 2 p 0 + &infin; , 满足性能指标 | | z | | 2 2 < &gamma; 2 | | w | | 2 2 , 如果存在对称正定矩阵以及矩阵 使得如下的矩阵不等式成立
&lambda; 2 &alpha; k c 1 2 + &alpha; k - 1 &gamma; 2 h 2 - &lambda; 1 c 2 2 < 0
G &lambda; 1 < P ~ i , &xi; k < G &lambda; 2
其中
M = &alpha;E ( &pi; i 1 ( &xi; k ) ) A ~ i T &alpha;E ( &pi; i 2 ( &xi; k ) ) A ~ i T . . . &alpha;E ( &pi; is ( &xi; k ) ) A ~ i T T
N = &alpha;E ( &pi; i 1 ( &xi; k ) ) B ~ &omega;i T &alpha;E ( &pi; i 2 ( &xi; k ) ) B ~ &omega;i T . . . &alpha;E ( &pi; is ( &xi; k ) ) B ~ &omega;i T T
A ~ i = A i + B i K i , &xi; k - B i K i , &xi; k 0 A i - H i , &xi; k E i = A i 0 n &times; n 0 n &times; n A i + B i 0 n &times; m K i , &xi; k I n &times; n - I n &times; n + 0 n &times; n - I n &times; n H i , &xi; k 0 p &times; n E i = &Theta; 1 + &Theta; 2 K i , &xi; k &Theta; 3 + &Theta; 4 H i , &xi; k &Theta; 5
B ~ &omega;i = B &omega;i B &omega;i = I n &times; n I n &times; n B &omega;i = &Theta; 6 B &omega;i
&Theta; 1 = A i 0 n &times; n 0 n &times; n A i , &Theta; 2 = B i 0 n &times; m , Θ3=[In×n -In×n],
&Theta; 4 = 0 n &times; n - I n &times; n , Θ5=[0p×n Ei], &Theta; 6 = I n &times; n I n &times; n
那么,闭环跳变系统有限时间稳定且满足H性能指标。
(4)仿真实验验证
为了验证系统状态不可测情形下,所提基于观测器H控制算法在有限时间内能将系统的状态轨迹限制在给定范围内,并具有一定的干扰抑制能力,本例考虑如下三模态参数的跳变系统描述:
A 1 = 0.88 - 0.05 0.4 - 0 . 72 , A 2 = 1 - 0.03 0 . 5 - 1.1 , A 3 = - 0.8 0 . 16 0.8 - 0 . 64 , B 1 = 2 1 , B 2 = 1 - 1 , B 3 = 1 1 B w 1 = 0.4 0.5 , B w 2 = 0.2 0 . 6 , B w 3 = 0.1 0 . 3 , C1=[-0.5 0.4],C2=[-0.3 0.1],C3=[-0.2 0.4]D1=0.9,D2=-0.6,D2=0.4,Dw1=0.5,Dw2=-0.5,Dw3=0.2
同时,反映三个模态之间跳变关系的转移概率矩阵的高斯概率密度函数分布给定如下:
N = n ( 0.2,0.05 ) n ( 0.3,0.05 ) n ( 0.5,0.05 ) n ( 0.2,0.05 ) n ( 0.4,0.05 ) n ( 0.4,0.05 ) n ( 0.2,0.05 ) n ( 0.6,0.05 ) n ( 0.2,0.05 )
先求取如下期望转移概率矩阵:
&Pi; = 0.2456 0.3049 0.4495 0.2471 0.3764 0.3764 0.2418 0.5165 0.2418
令γ=0.8,c1=1,R=I2,N=6,h2=10,α=0.4,根据结论中的不等式,利用MATLAB LMI Toolbox进行求解,得到各模态下的有限时间H控制器增益和观测器参数分别为:
K 1 , &pi; ij k = - 0.0846 - 0.0888 , K 2 , &pi; ij k = - 0 . 1476 - 0 . 1449 , K 3 , &pi; ij k = 0 . 3973 - 0 . 2276 H 1 , &pi; ij k = 1.5519 - 0.2538 , H 2 , &pi; ij k = 0.1787 - 0 . 4563 , H 3 , &pi; ij k = 0 . 1546 - 0 . 1683
假设系统的初始状态x0=[0.2 0.4]T,外界干扰输入为w(k)=0.1exp(-0.1k)sin(0.01πk),将设计的有限时间H控制器和观测器应用于步骤(1)的系统进行仿真实验,得到系统模态跳变图,状态轨迹图和观测误差图分别如图1,图2和图3所示。从图2中可以看出,在一定外界干扰下,经控制器作用后闭环控制系统的状态轨迹在有限时间内限定在指定范围c2内,同时从图3中可以看出观测器误差也在有限时间内收敛到平衡点。

Claims (2)

1.二重随机跳变系统基于观测器的有限短时间控制方法,其特征在于,所述方法包括以下步骤: 
(1)二重随机跳变过程描述 
考虑如下一类离散跳变系统: 
x(k+1)=A(rk)x(k)+B(rk)u(k)+Bw(rk)w(k) 
z(k)=C(rk)x(k)+D(rk)u(k)+Dw(rk)w(k) 
y(k)=E(rk)x(k) 
x(k)=x0,rk=r0,k=0 
其中,x(k)∈Rn是系统的状态向量;u(k)∈Rm是系统的控制向量;是外部扰动信号;z(k)∈Rl是系统的被控输出;A(rk),B(rk),Bw(rk),C(rk),D(rk),Dw(rk)和E(rk)分别为已知的与模态rk相关的适当维数的系数矩阵,其中rk表示系统的模态,为在有限集合Μ={1,2,…,s}中随时间k取值的二重随机跳变过程,其跳变转移概率定义如下: 
式中表示从模态i跳变到模态j的转移概率。为了方便起见,当rk=i时,分别用Ai,Bi,Bwi,Ci,Di,Dwi及Ei表征A(rk),B(rk),Bw(rk),C(rk),D(rk),Dw(rk)和E(rk)。 
不失一般性,用高斯随机分布{ξk,k∈K}来描述二重随机跳变过程的连续时变特性,其受限高斯概率密度函数表征为: 
式中f(·)为高斯概率密度函数的标准分布;F(·)为f(·)的累积分布函数,μij和σij为转移概率矩阵中各元素的高斯概率密度函数的均值和方差信息。基于上述描述,转移概率密度函数矩阵可表达为下式: 
其中的受限高斯概率密度函数。 
(2)二重随机特性下各子系统有限短时间稳定性定义: 
定义:设u(k)=0以及w(k)=0,被控二重随机跳变系统的各子系统是关于(c1 c2 N Ri)有限时间稳定的,其中c1是初始空间,c2是受限空间,且满足c1<c2,Ri>0,N为要求的时间常数,如果下列条件成立: 
对于上述定义,如果考虑系统受到外部干扰的影响,并假设干扰信号能量有界,则对于u(k)=0,被控系统的各子系统是关于(c1 c2 N Ri d)有限时间有界的,其中d为未知输入信号的上界,如果对于满足能量有界的干扰信号,均有上述条件成立。 
同样,如果对于上述定义,使用状态反馈控制,则被控系统的各子系统是关于(c1 c2 N Ri d)有限时间可镇定的,如果对于满足能量有界的干扰信号,均有上述条件成立。 
上述有限时间稳定的定义与Lyapunov意义下的渐进稳定是两个不同的概念,两者并无直接的关联,系统Lyapunov意义下的稳定并不能确保有限时间稳定;同样,有限时间稳定也不能保证Lyapunov渐进稳定。 
(3)基于观测器的有限短时间控制器设计: 
a.针对步骤(1)中构造的系统,设计全阶状态观测器; 
其中为观测器状态,为观测器估计输出,分别为待求的观测器和控制器参数。 
b.定义观测误差为:状态变量为:
则可得如下闭环误差动态系统: 
其中
c.选取李雅普诺夫泛函以及抗干扰性能指标
d.基于李雅普诺夫稳定性定理及有限短时间H控制方法,结合Markov跳变理论,利用已经获得的期望跳变转移概率矩阵,获取使闭环系统有限短时间稳定并满足H性能的充分条件。 
2.根据权利要求1所述的二重随机跳变系统基于观测器的有限短时间控制方法,其特征在于,设计观测器对系统的状态进行估计,再利用估计的状态设计控制器,确保系统有限短时间稳定的同时对所有频段的外部干扰信号具有一定的干扰抑制能力。 
CN201410488294.7A 2014-09-23 2014-09-23 二重随机跳变系统基于观测器的有限短时间控制方法 Pending CN104267601A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410488294.7A CN104267601A (zh) 2014-09-23 2014-09-23 二重随机跳变系统基于观测器的有限短时间控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410488294.7A CN104267601A (zh) 2014-09-23 2014-09-23 二重随机跳变系统基于观测器的有限短时间控制方法

Publications (1)

Publication Number Publication Date
CN104267601A true CN104267601A (zh) 2015-01-07

Family

ID=52159129

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410488294.7A Pending CN104267601A (zh) 2014-09-23 2014-09-23 二重随机跳变系统基于观测器的有限短时间控制方法

Country Status (1)

Country Link
CN (1) CN104267601A (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107065545A (zh) * 2017-04-01 2017-08-18 同济大学 基于马尔科夫跳变的分布式事件触发滤波系统及设计方法
CN109149604A (zh) * 2018-09-19 2019-01-04 湖南工业大学 一种保证电力系统稳定性的采样周期优化方法
CN109991848A (zh) * 2019-03-28 2019-07-09 鲁东大学 抗干扰控制方法和装置
CN112198802A (zh) * 2020-11-19 2021-01-08 南京工程学院 一种有限模态系统的有限时间抗扰异步估计方法

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107065545A (zh) * 2017-04-01 2017-08-18 同济大学 基于马尔科夫跳变的分布式事件触发滤波系统及设计方法
CN107065545B (zh) * 2017-04-01 2020-03-24 同济大学 基于马尔科夫跳变的分布式事件触发滤波系统及设计方法
CN109149604A (zh) * 2018-09-19 2019-01-04 湖南工业大学 一种保证电力系统稳定性的采样周期优化方法
CN109991848A (zh) * 2019-03-28 2019-07-09 鲁东大学 抗干扰控制方法和装置
CN109991848B (zh) * 2019-03-28 2022-02-01 鲁东大学 抗干扰控制方法和装置
CN112198802A (zh) * 2020-11-19 2021-01-08 南京工程学院 一种有限模态系统的有限时间抗扰异步估计方法
CN112198802B (zh) * 2020-11-19 2022-04-29 南京工程学院 一种有限模态系统的有限时间抗扰异步估计方法

Similar Documents

Publication Publication Date Title
Sui et al. Finite-time filter decentralized control for nonstrict-feedback nonlinear large-scale systems
Boltyanski et al. The robust maximum principle: theory and applications
Li et al. Observer-based adaptive fuzzy tracking control of MIMO stochastic nonlinear systems with unknown control directions and unknown dead zones
He et al. Finite-time $ H_ {\infty} $ fuzzy control of nonlinear jump systems with time delays via dynamic observer-based state feedback
CN103279031B (zh) 一种不确定多智能体系统的鲁棒趋同控制方法
Qiu et al. Output feedback guaranteed cost control for networked control systems with random packet dropouts and time delays in forward and feedback communication links
CN104267601A (zh) 二重随机跳变系统基于观测器的有限短时间控制方法
CN110376882A (zh) 基于有限时间扩张状态观测器的预定性能控制方法
CN105068420A (zh) 一种区间约束的非仿射不确定系统自适应控制方法
CN104298126A (zh) 非齐次Markov跳变控制系统的观测器设计方法
Yu et al. Quantized output feedback control of networked control systems with packet dropout
CN105372994A (zh) 一种基于鲁棒自适应最优控制的飞行员操纵行为描述方法
CN108564279A (zh) 一种考虑认知的生产线手工工位人因复杂性测度方法
CN106054606A (zh) 基于级联观测器的无模型控制方法
Munoz-Vázquez et al. Attitude tracking control of a quadrotor based on absolutely continuous fractional integral sliding modes
Fernández-González et al. Mixed discrete-continuous heuristic generative planning based on flow tubes
CN104238364A (zh) 二重随机跳变系统的有限短时间控制方法
CN104460336A (zh) 非齐次Markov跳变系统的干扰抑制方法
Wang et al. Finite-time consensus of nonlinear delayed multi-agent system via state-constraint impulsive control under switching topologies
Mukaidani Stackelberg strategy for discrete-time stochastic system and its application to H 2/H∞ control
Linsenmayer et al. Nonlinear event-triggered platooning control with exponential convergence
Shi et al. Neural adaptive intermittent output feedback control for autonomous underwater vehicles with full-state quantitative designs
Liu et al. Dissipativity-based asynchronous repetitive control for networked Markovian jump systems: 2-D system approach
Mohan et al. Performance comparison of LQR and ANFIS controller for stabilizing double inverted pendulum system
CN105824239A (zh) 带注视差异约束的平滑追踪眼动自适应跟踪控制器及方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
WD01 Invention patent application deemed withdrawn after publication
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20150107