CN105553442A - Network-based H-infinity filter information processing method for Lipschitz nonlinear system - Google Patents

Network-based H-infinity filter information processing method for Lipschitz nonlinear system Download PDF

Info

Publication number
CN105553442A
CN105553442A CN201510967000.3A CN201510967000A CN105553442A CN 105553442 A CN105553442 A CN 105553442A CN 201510967000 A CN201510967000 A CN 201510967000A CN 105553442 A CN105553442 A CN 105553442A
Authority
CN
China
Prior art keywords
filter
omega
lipschitz
network
estimated signal
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.)
Granted
Application number
CN201510967000.3A
Other languages
Chinese (zh)
Other versions
CN105553442B (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.)
JIANGSU EMPOWER INTELLIGENT TECHNOLOGY Co.,Ltd.
Original Assignee
Nantong 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 Nantong University filed Critical Nantong University
Priority to CN201510967000.3A priority Critical patent/CN105553442B/en
Publication of CN105553442A publication Critical patent/CN105553442A/en
Application granted granted Critical
Publication of CN105553442B publication Critical patent/CN105553442B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03HIMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
    • H03H17/00Networks using digital techniques
    • H03H17/02Frequency selective networks
    • H03H17/0248Filters characterised by a particular frequency response or filtering method

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Mathematical Physics (AREA)
  • Compression, Expansion, Code Conversion, And Decoders (AREA)

Abstract

The invention discloses a network-based H-infinity filter information processing method for a Lipschitz nonlinear system. The method comprises the steps of: 1) analyzing a networked Lipschitz nonlinear system interfered by noise whose statistical characteristics are difficult to determine; and performing quantization processing on a measuring output and an estimated signal of the system and transmitting the measuring output and the estimated signal to the filter over the network; 2) after the filter receives the data, separating the measuring output and the estimated signal; and discussing issues of the H-infinity filter for the Lipschitz nonlinear system based on the data transfer process described by the Markov chain; and 3) establishing a new networked Lipschitz nonlinear system and a filter model based on the state transfer probability information of the Markov chain, and estimating the system state based on the newly established filter model. According to the invention, issues, including quantization and packet loss, in the design of H-infinity filters for Lipschitz nonlinear systems are studied, an information processing method based on the filter is provided, and estimation of the filter for the system state is achieved.

Description

The H of network Lipschitz non linear system ∞filter information processing method
Technical field
The present invention relates to a kind of network control system, especially relate to the H of network Lipschitz non linear system filter information processing method.
Background technology
Along with the development of network technology, communication network application is in the controls more and more extensive.This control system by communication network formation closed loop is called network control system or network control system.In network control system, transducer, controller and actuator utilize network to carry out information exchange, therefore, network control system have convenient for installation and maintenance, flexibility is high and can realize the advantages such as Long-distance Control.
Because network control system introduces communication network in control loop, therefore create the problems such as data quantification and data-bag lost.If when the state of networked system can not be surveyed and be subject to noise disturbance, the state of designing filter to system is needed to estimate.When the noise be subject to by estimating system is white noise, can designs and adopt estimation error variance to estimate system mode as the Kalman filter of performance metric and linear filter.
At present, network design of filter mainly linear system, there is fan-shaped nonlinear system to unify the non linear system described by T-S fuzzy model.Lipschitz non linear system is the important non linear system of a class, such as, the robot system that there is trigonometric function item is Lipschitz non linear system, and the non linear system of some complexity can be described as linear system and the nonlinear terms sum meeting Lipschitz condition.But existing network design of filter Problems existing has: one, by estimated signal by wire instead of by Internet Transmission to long-range filter, the system of adding realizes cost, is unfavorable for system installation and maintenance.Two, main consideration network transfer delay and data-bag lost, and make great efforts distribution data of description transmitting procedure with shellfish, be not suitable for the situation of data by wireless network transmissions.Therefore, when Lipschitz non linear system be subject to statistical property be difficult to the noise disturbance determined time, measurement exports and is undertaken quantizing by estimated signal and by Internet Transmission to long-range filter by the present invention, and by Markov chain data of description transmitting procedure, devises network H filter.
Summary of the invention
The object of the invention is the deficiency existed for prior art, propose a kind of H of network Lipschitz non linear system filter information processing method.
The technical scheme that the present invention takes is as follows:
The H of network Lipschitz non linear system filter information processing method, comprises the steps:
1) first the networking Lipschitz non linear system being subject to the noise jamming that statistical property is difficult to determine is analyzed; And the measurement of system is exported and carried out quantification treatment by estimated signal, be sent to filter end by network;
2) after filter receives data, measurement is exported and be separated with by estimated signal; And based on the data transmission procedure that Markov chain describes, the H of Lipschitz non linear system is discussed filter problem;
3) based on the state transition probability information of Markov chain, set up novel networking Lipschitz non linear system and filter model, and estimate according to the state of newly-established filter model to system.
Described step 1) in, the state equation that Lipschitz non linear system is set up is as follows
x(k+1)=Ax(k)+Ff(k,x(k))+Bw(k)
y(k)=Cx(k)+Gg(k,x(k))+Dw(k)(1)
z(k)=Lx(k)
Wherein: x (k) ∈ R nstate vector, y (k) ∈ R pmeasure to export, z (k) ∈ R qby estimated signal, w (k) ∈ R mspace L 2[0, ∞) on noise signal, A, F, B, C, G, D and L are the coefficient matrixes of known corresponding dimension, f (k, x (k)) and g (k, x (k)) be the Nonlinear Vector function of satisfied following Lipschitz condition: f (k, 0)=0, g (k, 0)=0, || f (k, x (k)) ||≤|| F 1x (k) ||, || g (k, x (k)) ||≤|| G 1x (k) ||, wherein F 1, G 1it is known respective dimension matrix number.
Described step 1) in, the measurement of system to be exported and carried out the method for quantification treatment by estimated signal specific as follows:
Measure and export y (k) and after quantizing, broken into a packet by Internet Transmission to filter node by the element of estimated signal z (k), first quantize it before transmission, quantizer is:
Q ( v ) = u i , 1 1 + &delta; u i < v &le; 1 1 - &delta; u i , v > 0 0 , v = 0 - Q ( - v ) , v < 0 - - - ( 2 )
Wherein, the quantization level collection U={ ± u of system i, u iiu 0, i=± 1, ± 2 ... ∪ { ± u 0∪ { 0}, 0< ρ <1, u 0>0, ρ is the quantization resolution of quantizer Q (v).
Described step 2) in, the data transmission procedure method that Markov chain describes is as follows:
σ (k)=0 represents that data are not lost by during Internet Transmission; σ (k)=1 represents that data are lost by producing during Internet Transmission; State transition probability matrix P=[the p of Markov chain ij], wherein p ij>=0, work as p 01+ p 10when=1, transmitting procedure for shellfish make great efforts distribution;
When data are by network normal transmission, then after filter receives data, measurement exported and be separated with by estimated signal; Filter input y f(k)=Q (y (k)), filter receive by estimated signal z c(k)=Q (z (k)); The quantizer of system adopts logarithmic quantization device, therefore y f(k)=(I+H (k)) y (k), z c(k)=(I+H (k)) z (k), wherein I is the unit matrix of corresponding dimension, and H (k) is uncertain matrix, and meets || H (k) ||≤δ;
When data occur to lose in transmitting procedure, filter does not receive the measurement that network transmits and exports and by estimated signal, now suppose that filter input keeps the value of previous moment constant, i.e. y f(k)=y f(k-1), filter does not receive by estimated signal, i.e. z c(k)=0, what such filter inputted and received can be turned to by estimated signal:
y f(k)=(1-σ(k))(I+H(k))y(k)+σ(k)y f(k-1)(3)
z c(k)=(1-σ(k))(I+H(k))z(k)(4)。
Described step 2) H of Lipschitz non linear system is discussed the method of filter problem is specific as follows:
For output z (k) of estimator (1), design the pattern-dependent filter of following form:
x f(k+1)=A fσ(k)x f(k)+B fσ(k)y f(k)
(5)
z f(k)=C fσ(k)x f(k)
Wherein: x f(k) ∈ R nfilter status vector, y f(k) ∈ R pfilter input, z f(k) ∈ R qestimated signal, A f σ (k), B f σ (k)and C f σ (k), { 0,1} is filter parameter undetermined to σ (k)=i ∈;
Definition vector x ~ ( k ) = x T ( k ) x f T ( k ) y f T ( k ) T , η (k)=[f t(k, x (k)) g t(k, x (k))] t, evaluated error e (k)=z c(k)-z fk (), by formula (1), formula (3), formula (4) and formula (5), can obtain following filtering error system:
x ~ ( k + 1 ) = A ~ &sigma; ( k ) x ~ ( k ) + F ~ &sigma; ( k ) &eta; ( k ) + B ~ &sigma; ( k ) w ( k ) e ( k ) = C ~ &sigma; ( k ) x ~ ( k ) , &sigma; ( k ) = i &Element; { 0 , 1 } - - - ( 6 )
Wherein: A ~ 0 = A ^ 0 + D ^ H ( k ) E 1 , F ~ 0 = F ^ 0 + D ^ H ( k ) E 2 , B ~ 0 = B ^ 0 + D ^ H ( k ) E 3 , C ~ 0 = C ^ + H ( k ) E 4 , A ^ 0 = A 0 0 B f 0 A f 0 0 C 0 0 , F ^ 0 = F 0 0 B f 0 G 0 G , B ^ 0 = B B f 0 D D , D ^ = 0 B f 0 I , E 1=[C00],E 2=[0G],E 3=D, C ^ = L - C f 0 0 , E 4=[L00], A ~ 1 = A 0 0 0 A f 1 B f 1 0 0 I , F ~ 1 = I 0 0 0 0 0 , B ~ 1 = B 0 0 , C ~ 1 = 0 - C f 1 0 .
Described step 3), set up novel networking Lipschitz non linear system and filter model specific as follows:
For given positive number γ and quantization resolution ρ, if there is symmetric positive definite matrix P i, M i, R iand matrix N i, V i, B fi, i=0,1, positive scalar ε 1, ε 2and ε 3, following linear MATRIX INEQUALITIES (7) and (8) are set up, then filtering error system (6) formula is Stochastic stable and has H performance γ; Wherein filter parameter
&Pi; 11 * * * * * * * &Pi; 21 &Pi; 22 * * * * * * &Pi; 31 &Pi; 3 2 &Pi; 33 * * * * * &Pi; 41 &Pi; 4 2 &Pi; 43 &Pi; 44 * * * * &Pi; 51 &Pi; 5 2 &Pi; 53 0 &Pi; 55 * * * &Pi; 61 0 0 0 0 &Pi; 66 * * 0 0 0 &Pi; 74 &Pi; 75 0 &Pi; 77 * 0 0 0 0 0 &Pi; 86 0 &Pi; 88 < 0 - - - ( 7 )
&Omega; 11 * * * * * 0 &Omega; 22 * * * * 0 0 &Omega; 33 * * * &Omega; 41 &Omega; 42 &Omega; 43 &Omega; 44 * * &Omega; 51 &Omega; 52 &Omega; 53 0 &Omega; 55 * &Omega; 61 0 0 0 0 &Omega; 66 < 0 - - - ( 8 )
Wherein: &Pi; 11 = d i a g - P 0 + &epsiv; 1 ( F 1 T F 1 + G 1 T G 1 ) + &epsiv; 2 &delta; 2 C T C + &epsiv; 3 &delta; 2 L T L , - M 0 , - R 0 ,
&Pi; 21 = 0 0 0 &epsiv; 2 &delta; 2 G T C 0 0 , Π 22=diag{-ε 1I,-ε 1I+ε 2δ 2G TG},
Π 31=[ε 2δ 2D TC00],Π 32=[0ε 2δ 2D TG],Π 33=-γ 2I+ε 2δ 2D TD,
&Pi; 41 = P 0 A 0 0 B f 0 C N 0 0 R 0 C 0 0 , &Pi; 42 = P 0 F 0 0 B f 0 G 0 R 0 G , &Pi; 43 = P 0 B B f 0 D R 0 D ,
&Pi; 44 = d i a g - p 00 - 1 P 0 , - p 00 - 1 M 0 , - p 00 - 1 R 0 , , &Pi; 51 = P 1 A 0 0 B f 0 C N 0 0 R 1 C 0 0 ,
&Pi; 52 = P 1 F 0 0 B f 0 G 0 R 1 G , &Pi; 53 = P 1 B B f 0 D R 1 D , &Pi; 55 = d i a g - p 01 - 1 P 1 , - p 01 - 1 M 1 , - p 01 - 1 R 1 ,
Π 61=[L-V 00],Π 66=-I, &Pi; 74 = 0 B f 0 T I , Π 75=Π 74,Π 77=-ε 2I,
Π 86=I,Π 88=-ε 3I, &Omega; 11 = d i a g - P 1 + &epsiv; 1 ( F 1 T F 1 + G 1 T G 1 ) , - M 1 , - R 1 , Ω 22=-ε 1I,Ω 33=-γ 2I,
&Omega; 41 = P 0 A 0 0 0 N 1 B f 1 0 0 R 0 , &Omega; 42 = P 0 F 0 0 0 0 0 , &Omega; 43 = P 0 B 0 0 ,
&Omega; 44 = d i a g - p 10 - 1 P 0 , - p 10 - 1 M 0 , - p 10 - 1 R 0 , &Omega; 51 = P 1 A 0 0 0 N 1 B f 1 0 0 R 1 , &Omega; 52 = P 1 F 0 0 0 0 0 , &Omega; 53 = P 1 B 0 0 ,
&Omega; 55 = d i a g - p 11 - 1 P 1 , - p 11 - 1 M 1 , - p 11 - 1 R 1 , Ω 61=[0-V 10],Ω 66=-I。
Advantage of the present invention and remarkable result specific as follows:
(1) measurement of Lipschitz non linear system exported and undertaken quantizing by estimated signal and by Internet Transmission to long-range filter, devise network H filter, compared to existing technology, what reduce filtering system realizes cost, improve convenience and the flexibility of system installation and maintenance, even if data are lost by existing in network transmission process, system is subject to the noise jamming that statistical property is difficult to determine, estimated signal also can be followed the tracks of by estimated signal preferably.
(2) the present invention adopts Markov chain data of description transmitting procedure, makes great efforts distribution data of description transmitting procedure and has more generality, be particularly useful for the situation of data by wireless network transmissions compared to employing shellfish.
(3) the present invention utilizes Lyapunov Functional Approach, have studied the H that a class has the Lipschitz non linear system of quantification and data-bag lost filter design problem, and give the information processing method based on this filter, achieve the estimation of filter to system mode.
Accompanying drawing explanation
Fig. 1 is network non linear system filter structure figure;
Fig. 2 is data transmission procedure figure;
Fig. 3 is the drawing for estimate of filter to system mode.
Embodiment
Below in conjunction with accompanying drawing, technical scheme of the present invention is described in detail.
As shown in Figure 1, Figure 2 and Figure 3, be respectively network non linear system filter structure figure, data transmission procedure figure and filter to the drawing for estimate of system mode, the present invention is based on the H of the Lipschitz non linear system of network filter information processing method, comprises the steps:
1) first the networking Lipschitz non linear system being subject to the noise jamming that statistical property is difficult to determine is analyzed; And the measurement of system is exported and carried out quantification treatment by estimated signal, be sent to filter end by network;
2) after filter receives data, measurement is exported and be separated with by estimated signal; And based on the data transmission procedure that Markov chain describes, the H of Lipschitz non linear system is discussed filter problem;
3) based on the state transition probability information of Markov chain, set up novel networking Lipschitz non linear system and filter model, and estimate according to the state of newly-established filter model to system.
Wherein, step 1) in be subject to the noise jamming that statistical property is difficult to determine networking Lipschitz non linear system set up state equation such as formula (1).The measurement of system is exported and quantized by estimated signal and the method for process of pack as follows: measure and export y (k) and after quantification, broken into a packet by Internet Transmission to filter node by the element of estimated signal z (k).First quantize it before transmission, quantizer is formula (2).
Step 2) in Markov chain describe data transmission procedure method as follows: represent that data are not lost by during Internet Transmission with σ (k)=0; σ (k)=1 represents that data are lost by producing during Internet Transmission.State transition probability matrix P=[the p of Markov chain ij], wherein p i j = Pr &sigma; ( k + 1 ) = j | &sigma; ( k ) = i &ForAll; i , j &Element; { 0 , 1 } , p ij≥0, &Sigma; j = 0 1 p i j = 1. Work as p 01+ p 10when=1, transmitting procedure for shellfish make great efforts distribution, therefore have more generality with Markov chain data of description transmitting procedure, be especially applicable to data of description pass through wireless network transmissions situation.
When data are by network normal transmission, then after filter receives data, measurement exported and be separated with by estimated signal.Filter input y f(k)=Q (y (k)), filter receive by estimated signal z c(k)=Q (z (k)).Quantizer due to system have employed logarithmic quantization device, therefore y f(k)=(I+H (k)) y (k), z c(k)=(I+H (k)) z (k), wherein || H (k) ||≤δ.When data occur to lose in transmitting procedure, filter does not receive the measurement that network transmits and exports and by estimated signal, now suppose that filter input keeps the value of previous moment constant, i.e. y f(k)=y f(k-1), filter does not receive by estimated signal, i.e. z c(k)=0.The formula (3) that can be turned to by estimated signal of such filter input and reception and formula (4)
For output z (k) of estimator (1), design the pattern-dependent filter such as formula (5).Definition vector x ~ ( k ) = x T ( k ) x f T ( k ) y f T ( k ) T , η (k)=[f t(k, x (k)) g t(k, x (k))] t, evaluated error e (k)=z c(k)-z fk (), by formula (1), formula (3), formula (4) and formula (5), can obtain filtering error systematic (6).
Step 3) set up novel networking Lipschitz non linear system and filter model specific as follows: for given positive number γ and quantization resolution ρ, if there is symmetric positive definite matrix P i, M i, R i, matrix N i, V i, B fi, i=0,1, positive scalar ε 1, ε 2and ε 3, LMI (7) and (8) are set up, then filtering error system (6) is Stochastic stable and has H performance γ; Wherein filter parameter
The given matrix W with corresponding dimension of lemma, D and E, wherein W is symmetrical matrix.F is met for all tthe matrix F (k) of (k) F (k)≤I, W+DF (k) E+E tf t(k) D tthe necessary and sufficient condition that <0 sets up there is ε >0, makes W+ ε DD t+ ε -1e te<0.
Select following Lyapunov functional:
V &lsqb; x ~ ( k ) , &sigma; ( k ) = i &rsqb; = x ~ T ( k ) S i x ~ ( k ) = x T ( k ) P i x ( k ) + x f T ( k ) Q i x f ( k ) + y f T ( k - 1 ) R i y f ( k - 1 )
When w (k)=0, difference along system (6) is:
&Delta; V ( k ) = E { V &lsqb; x ~ ( k + 1 ) , &sigma; ( k + 1 ) &rsqb; | x ~ ( k ) , &sigma; ( k ) = i } - V &lsqb; x ~ ( k ) , &sigma; ( k ) = i &rsqb; = x ~ T ( k + 1 ) &Sigma; j = 0 1 p i j S i x ~ ( k + 1 ) - x ~ T ( k ) S i x ~ ( k ) = x ~ T ( k ) A ~ i T &Sigma; j = 0 1 p i j S j A ~ i x ~ ( k ) + 2 x ~ T ( k ) A ~ i T &Sigma; j = 0 1 p i j S j F ~ i &eta; ( k ) + &eta; T ( k ) F ~ i T &Sigma; j = 0 1 p i j S j F ~ i &eta; ( k ) - x ~ T ( k ) S i x ~ ( k )
Due to
&eta; T ( k ) &eta; ( k ) = f T ( k , x ( k ) ) f ( k , x ( k ) ) + g T ( k , x ( k ) ) g ( k , x ( k ) ) &le; x T ( k ) ( F 1 T F 1 + G 1 T G 1 ) x ( k ) , So there is ε 1>0, makes &epsiv; 1 x T ( k ) ( F 1 T F 1 + G 1 T G 1 ) x ( k ) - &epsiv; 1 &eta; T ( k ) &eta; ( k ) &GreaterEqual; 0 , Then
&Delta; V ( k ) &le; x ~ T ( k ) A ~ i T &Sigma; j = 0 1 p i j S j A ~ i x ~ ( k ) + 2 x ~ T ( k ) A ~ i T &Sigma; j = 0 1 p i j S j F ~ i &eta; ( k ) + &eta; T ( k ) F ~ i T &Sigma; j = 0 1 p i j S j F ~ i &eta; ( k ) - x ~ T ( k ) S i x ~ ( k ) + &epsiv; 1 x T ( k ) ( F 1 T F 1 + G 1 T G 1 ) x ( k ) - &epsiv; 1 &eta; T ( k ) &eta; ( k ) = &xi; T ( k ) &Theta; i &xi; ( k )
Wherein &xi; ( k ) = x ~ T ( k ) &eta; T ( k ) T , &Theta; i = A ~ i T &Sigma; j = 0 1 p i j S j A ~ i - S i + W A ~ i T &Sigma; j = 0 1 p i j S j F ~ F ~ i T &Sigma; j = 0 1 p i j S j A ~ i F ~ i T &Sigma; j = 0 1 p i j S j F ~ i - &epsiv; 1 I , W = I 0 0 T &epsiv; 1 ( F 1 T F 1 + G 1 T G 1 ) I 0 0 .
If n 0=A f0m 0, V 0=C f0m 0, N 1=A f1m 1, V 1=C f1m 1,
(7) formula both sides are multiplied by respectively d i a g { I , Q 0 , I , I , I , I , P 0 - 1 , I , R 0 - 1 , P 1 - 1 , I , R 1 - 1 , I , I , I } :
- S 0 + W + &epsiv; 2 &delta; 2 E 1 T E 1 + &epsiv; 3 &delta; 2 E 4 T E 4 * * * * * * * &epsiv; 2 &delta; 2 E 2 T E 1 - &epsiv; 1 I + &epsiv; 2 &delta; 2 E 2 T E 2 * * * * * * &epsiv; 2 &delta; 2 E 3 T E 1 &epsiv; 2 &delta; 2 E 3 T E 2 - &gamma; 2 I + &epsiv; 2 &delta; 2 E 3 T E 3 * * * * * A ^ 0 F ^ 0 B ^ 0 - p 00 - 1 S 0 - 1 * * * * A ^ 0 F ^ 0 B ^ 0 0 - p 01 - 1 S 1 - 1 * * * C ^ 0 0 0 0 - I * * 0 0 0 D ^ T D ^ T 0 - &epsiv; 2 I * 0 0 0 0 0 I 0 - &epsiv; 3 I < 0
Namely
- S 0 + W + &epsiv; 2 &delta; 2 E 1 T E 1 * * * * * * &epsiv; 2 &delta; 2 E 2 T E 1 - &epsiv; 1 I + &epsiv; 2 &delta; 2 E 2 T E 2 * * * * * &epsiv; 2 &delta; 2 E 3 T E 1 &epsiv; 2 &delta; 2 E 3 T E 2 - &gamma; 2 I + &epsiv; 2 &delta; 2 E 3 T E 3 * * * * A ^ 0 F ^ 0 B ^ 0 - p 00 - 1 S 0 - 1 * * * A ^ 0 F ^ 0 B ^ 0 0 - p 01 - 1 S 1 - 1 * * C ^ 0 0 0 0 - I * 0 0 0 D ^ T D ^ T 0 - &epsiv; 3 I
+ &epsiv; 3 - 1 0 0 0 0 0 I 0 0 0 0 0 0 I 0 + &epsiv; 3 &delta; 2 E 4 T 0 0 0 0 0 0 E 4 0 0 0 0 0 0 < 0
Will substitute into, according to lemma, above formula can turn to:
- S 0 + W + &epsiv; 2 &delta; 2 E 1 T E 1 * * * * * * &epsiv; 2 &delta; 2 E 2 T E 1 - &epsiv; 1 I + &epsiv; 2 &delta; 2 E 2 T E 2 * * * * * &epsiv; 2 &delta; 2 E 3 T E 1 &epsiv; 2 &delta; 2 E 3 T E 2 - &gamma; 2 I + &epsiv; 2 &delta; 2 E 3 T E 3 * * * * A ^ 0 F ^ 0 B ^ 0 - p 00 - 1 S 0 - 1 * * * A ^ 0 F ^ 0 B ^ 0 0 - p 01 - 1 S 1 - 1 * * C ^ 0 0 0 0 0 - I * 0 0 0 D ^ T D ^ T 0 - &epsiv; 3 I < 0 Namely
- S 0 + W * * * * * 0 - &epsiv; 1 I * * * * 0 0 - &gamma; 2 I * * * A ^ 0 F ^ 0 B ^ 0 - p 00 - 1 S 0 - 1 * * A ^ 0 F ^ 0 B ^ 0 0 - p 01 - 1 S 1 - 1 * C ^ 0 0 0 0 0 - I + &epsiv; 2 - 1 0 0 0 D ^ D ^ 0 0 0 0 D ^ T D ^ T 0
+ &epsiv; 2 &delta; 2 E 1 T E 2 T E 3 T 0 0 0 E 1 E 2 E 3 0 0 0 < 0
Will A ~ 0 = A ^ 0 + D ^ H ( k ) E 1 , F ~ 0 = F ^ 0 + D ^ H ( k ) E 2 , B ~ 0 = B ^ 0 + D ^ H ( k ) E 3 Substitute into, according to lemma, above formula can turn to:
- S 0 + W * * * * * 0 - &epsiv; 1 I * * * * 0 0 - &gamma; 2 I * * * A ~ 0 F ~ 0 B ~ 0 - p 00 - 1 S 0 - 1 * * A ~ 0 F ~ 0 B ~ 0 0 - p 01 - 1 S 1 - 1 * C ~ 0 0 0 0 0 - I < 0
Mend lemma by Schur, above formula can turn to:
&Xi; 0 = A ~ 0 T &Sigma; j = 0 1 p i j S j A ~ 0 - S 0 + W + C ~ 0 T C ~ 0 * * F ~ 0 T &Sigma; j = 0 1 p i j S j A ~ 0 F ~ 0 T &Sigma; j = 0 1 p i j S j F ~ 0 - &epsiv; 1 I * F ~ 0 T &Sigma; j = 0 1 p i j S j A ~ 0 B ~ 0 T &Sigma; j = 0 1 p i j S j F ~ 0 B ~ 0 T &Sigma; j = 0 1 p i j S j B ~ 0 - &gamma; 2 I < 0
Can be obtained fom the above equation:
&Theta; 0 = A ~ 0 T &Sigma; j = 0 1 p i j S j A ~ 0 - S 0 + W * F ~ 0 T &Sigma; j = 0 1 p i j S j A ~ 0 F ~ 0 T &Sigma; j = 0 1 p i j S j A ~ 0 - &epsiv; 1 I < 0
In like manner, by (8) Shi Ke get:
&Xi; 1 = A ~ 1 T &Sigma; j = 0 1 p i j S j A ~ 1 - S 1 + W + C ~ 1 T C ~ 1 * * F ~ 1 T &Sigma; j = 0 1 p i j S j A ~ 1 F ~ 1 T &Sigma; j = 0 1 p i j S j F ~ 1 - &epsiv; 1 I * F ~ 1 T &Sigma; j = 0 1 p i j S j A ~ 1 B ~ 1 T &Sigma; j = 0 1 p i j S j F ~ 1 B ~ 1 T &Sigma; j = 0 1 p i j S j B ~ 1 - &gamma; 2 I < 0
Can obtain further:
&Theta; 1 = A ~ 1 T &Sigma; j = 0 1 p i j S j A ~ 1 - S 1 + W * F ~ 1 T &Sigma; j = 0 1 p i j S j A ~ 1 F ~ 1 T &Sigma; j = 0 1 p i j S j F ~ 1 - &epsiv; 1 I < 0.
Therefore,
E { V &lsqb; x ~ ( k + 1 ) , &sigma; ( k + 1 ) &rsqb; | x ~ ( k ) , &sigma; ( k ) = i } - V &lsqb; x ~ ( k ) , &sigma; ( k ) = i &rsqb; &le; &xi; T ( k ) &Theta; i &xi; ( k ) &le; &beta;&xi; T ( k ) &xi; ( k ) < 0
Wherein for Θ ieigenvalue of maximum.
When N >=1,
E { V &lsqb; x ~ ( N + 1 ) , &sigma; ( N + 1 ) &rsqb; } - V ( x ~ 0 , &sigma; 0 ) &le; &Sigma; k = 0 N &beta; E &lsqb; &xi; T ( k ) &xi; ( k ) &rsqb; &le; &Sigma; k = 0 N &beta; E &lsqb; x ~ T ( k ) x ~ ( k ) &rsqb;
Above formula can turn to:
&Sigma; k = 0 N E &lsqb; x ~ T ( k ) x ~ ( k ) &rsqb; &le; &beta; - 1 E { V &lsqb; x ~ ( N + 1 ) , &sigma; ( N + 1 ) &rsqb; } - &beta; - 1 V ( x ~ 0 , &sigma; 0 )
As N → ∞, can obtain:
&Sigma; k = 0 &infin; E &lsqb; x ~ T ( k ) x ~ ( k ) &rsqb; < - &beta; - 1 V ( x ~ 0 , &sigma; 0 ) < &infin;
Therefore, filtering error system (6) is Stochastic stable.
Prove that filtering error system (6) has H below performance γ.Under zero initial condition, definition
J N = E { &Sigma; k = 0 N &lsqb; e T ( k ) e ( k ) - &gamma; 2 w T ( k ) w ( k ) &rsqb; }
Above formula can be expressed as:
Wherein
As N → ∞, can obtain:
J &infin; = E { &Sigma; k = 0 &infin; &lsqb; z T ( k ) z ( k ) - &gamma; 2 w T ( k ) w ( k ) &rsqb; } < 0
Therefore, filtering error system (6) has H performance γ.
Below according to example, verify validity of the present invention and superiority.
In FIG, the such as Lipschitz non linear system shown in formula (1), parameter is as follows:
A = 0.2 0.3 0.1 0.1 , F = 0.3 0.4 0.4 0.2 , B = 0.4 0.1 , C=[0.10.1], G=[0.10.2], D=0.3, L=[11], nonlinear function f (k, x (k))=0.1sin (x (k)), g (k, x (k))=0.1sin (x (k)), quantization resolution ρ=0.5 of quantizer Q (v), H performance γ=1.8, data transmission procedure is as shown in Figure 2, its state transition probability matrix P = 0.8 0.2 0.7 0.3 , Noise signal w ( k ) = 0.5 , 5 < k &le; 15 - 0.5 , 55 < k &le; 65 0 , o t h e r w i s e .
Utilize the LMI tool box in MATLAB to solve LMI formula (7) and formula (8) filter parameter can be obtained be:
A f 0 = 0.1106 0.3278 0.1357 0.4021 , B f 0 = 2.6164 4.6841 , C f0=[0.10130.1001],
A f 1 = 0.1203 0.0319 0.1583 0.2794 , B f 1 = 4.8483 8.8104 , C f1=[0.11950.0962]。

Claims (6)

1. the H of network Lipschitz non linear system filter information processing method, is characterized in that, comprises the steps:
1) first the networking Lipschitz non linear system being subject to the noise jamming that statistical property is difficult to determine is analyzed; And the measurement of system is exported and carried out quantification treatment by estimated signal, be sent to filter end by network;
2) after filter receives data, measurement is exported and be separated with by estimated signal; And based on the data transmission procedure that Markov chain describes, the H of Lipschitz non linear system is discussed filter problem;
3) based on the state transition probability information of Markov chain, set up novel networking Lipschitz non linear system and filter model, and estimate according to the state of newly-established filter model to system.
2. the H of network Lipschitz non linear system according to claim 1 design of filter, is characterized in that, described step 1) in, the state equation that Lipschitz non linear system is set up is as follows
x(k+1)=Ax(k)+Ff(k,x(k))+Bw(k)
y(k)=Cx(k)+Gg(k,x(k))+Dw(k)(1)
z(k)=Lx(k)
Wherein: x (k) ∈ R nstate vector, y (k) ∈ R pmeasure to export, z (k) ∈ R qby estimated signal, w (k) ∈ R mspace L 2[0, ∞) on noise signal, A, F, B, C, G, D and L are the coefficient matrixes of known corresponding dimension, f (k, x (k)) and g (k, x (k)) be the Nonlinear Vector function of satisfied following Lipschitz condition: f (k, 0)=0, g (k, 0)=0, || f (k, x (k)) ||≤|| F 1x (k) ||, || g (k, x (k)) ||≤|| G 1x (k) ||, wherein F 1, G 1it is known respective dimension matrix number.
3. the H of network Lipschitz non linear system according to claim 2 design of filter, is characterized in that, described step 1) in, the measurement of system to be exported and carried out the method for quantification treatment by estimated signal specific as follows:
Measure and export y (k) and after quantizing, broken into a packet by Internet Transmission to filter node by the element of estimated signal z (k), first quantize it before transmission, quantizer is:
Q ( v ) = u i , 1 1 + &delta; u i < v &le; 1 1 - &delta; u i , v > 0 0 , v = 0 - Q ( - v ) , v < 0 - - - ( 2 )
Wherein, the quantization level collection U={ ± u of system i, u iiu 0, i=± 1, ± 2 ... ∪ { ± u 0∪ { 0}, 0< ρ <1, u 0>0, ρ is the quantization resolution of quantizer Q (v).
4. the H of network Lipschitz non linear system according to claim 3 design of filter, is characterized in that, described step 2) in, the data transmission procedure method that Markov chain describes is as follows:
σ (k)=0 represents that data are not lost by during Internet Transmission; σ (k)=1 represents that data are lost by producing during Internet Transmission; State transition probability matrix P=[the p of Markov chain ij], wherein p i j = Pr { &sigma; ( k + 1 ) = j | &sigma; ( k ) = i } &ForAll; i , j &Element; { 0 , 1 } , P ij>=0, work as p 01+ p 10when=1, transmitting procedure for shellfish make great efforts distribution;
When data are by network normal transmission, then after filter receives data, measurement exported and be separated with by estimated signal; Filter input y f(k)=Q (y (k)), filter receive by estimated signal z c(k)=Q (z (k)); The quantizer of system adopts logarithmic quantization device, therefore y f(k)=(I+H (k)) y (k), z c(k)=(I+H (k)) z (k), wherein I is the unit matrix of corresponding dimension, and H (k) is uncertain matrix, and meets || H (k) ||≤δ;
When data occur to lose in transmitting procedure, filter does not receive the measurement that network transmits and exports and by estimated signal, now suppose that filter input keeps the value of previous moment constant, i.e. y f(k)=y f(k-1), filter does not receive by estimated signal, i.e. z c(k)=0, what such filter inputted and received can be turned to by estimated signal:
y f(k)=(1-σ(k))(I+H(k))y(k)+σ(k)y f(k-1)(3)
z c(k)=(1-σ(k))(I+H(k))z(k)(4)。
5. the H of network Lipschitz non linear system according to claim 4 design of filter, is characterized in that, described step 2) H of Lipschitz non linear system is discussed the method of filter problem is specific as follows:
For output z (k) of estimator (1), design the pattern-dependent filter of following form:
x f(k+1)=A fσ(k)x f(k)+B fσ(k)y f(k)
z f(k)=C fσ(k)x f(k)(5)
Wherein: x f(k) ∈ R nfilter status vector, y f(k) ∈ R pfilter input, z f(k) ∈ R qestimated signal, A f σ (k), B f σ (k)and C f σ (k), { 0,1} is filter parameter undetermined to σ (k)=i ∈;
Definition vector x ~ ( k ) = x T ( k ) x f T ( k ) y f T ( k ) T , &eta; ( k ) = f T ( k , x ( k ) ) g T ( k , x ( k ) ) T , Evaluated error e (k)=z c(k)-z fk (), by formula (1), formula (3), formula (4) and formula (5), can obtain following filtering error system:
x ~ ( k + 1 ) = A ~ &sigma; ( k ) x ~ ( k ) + F ~ &sigma; ( k ) &eta; ( k ) + B ~ &sigma; ( k ) w ( k ) e ( k ) = C ~ &sigma; ( k ) x ~ ( k ) , &sigma; ( k ) = i &Element; { 0 , 1 } - - - ( 6 )
Wherein: A ~ 0 = A ^ 0 + D ^ H ( k ) E 1 , F ~ 0 = F ^ 0 + D ^ H ( k ) E 2 , B ~ 0 = B ^ 0 + D ^ H ( k ) E 3 , C ~ 0 = C ^ + H ( k ) E 4 ,
A ^ 0 = A 0 0 B f 0 A f 0 0 C 0 0 , F ^ 0 = F 0 0 B f 0 G 0 G , B ^ 0 = B B f 0 D D , D ^ = 0 B f 0 I , E 1=[C00],
E 2=[0G],E 3=D, C ^ = L - C f 0 0 , E 4=[L00], A ~ 1 = A 0 0 0 A f 1 B f 1 0 0 I , F ~ 1 = I 0 0 0 0 0 ,
B ~ 1 = B 0 0 , C ^ 1 = 0 - C f 1 0 .
6. the H of network Lipschitz non linear system according to claim 5 design of filter, is characterized in that, described step 3), set up novel networking Lipschitz non linear system and filter model specific as follows:
For given positive number γ and quantization resolution ρ, if there is symmetric positive definite matrix P i, M i, R iand matrix N i, V i, B fi, i=0,1, positive scalar ε 1, ε 2and ε 3, following linear MATRIX INEQUALITIES (7) and (8) are set up, then filtering error system (6) formula is Stochastic stable and has H performance γ; Wherein filter parameter
&Pi; 11 * * * * * * * &Pi; 2 1 &Pi; 22 * * * * * * &Pi; 3 1 &Pi; 32 &Pi; 33 * * * * * &Pi; 4 1 &Pi; 42 &Pi; 43 &Pi; 44 * * * * &Pi; 5 1 &Pi; 52 &Pi; 53 0 &Pi; 55 * * * &Pi; 6 1 0 0 0 0 &Pi; 66 * * 0 0 0 &Pi; 74 &Pi; 75 0 &Pi; 77 * 0 0 0 0 0 &Pi; 86 0 &Pi; 88 < 0 - - - ( 7 )
&Omega; 11 * * * * * 0 &Omega; 22 * * * * 0 0 &Omega; 33 * * * &Omega; 41 &Omega; 42 &Omega; 43 &Omega; 44 * * &Omega; 51 &Omega; 52 &Omega; 53 0 &Omega; 55 * &Omega; 61 0 0 0 0 &Omega; 66 < 0 - - - ( 8 )
Wherein: &Pi; 11 = d i a g - P 0 + &epsiv; 1 ( F 1 T F 1 + G 1 T G 1 ) + &epsiv; 2 &delta; 2 C T C + &epsiv; 3 &delta; 2 L T L , - M 0 , - R 0 ,
&Pi; 21 = 0 0 0 &epsiv; 2 &delta; 2 G T C 0 0 , Π 22=diag{-ε 1I,-ε 1I+ε 2δ 2G TG},
Π 31=[ε 2δ 2D TC00],Π 32=[0ε 2δ 2D TG],Π 33=-γ 2I+ε 2δ 2D TD,
&Pi; 41 = P 0 A 0 0 B f 0 C N 0 0 R 0 C 0 0 , &Pi; 4 2 = P 0 F 0 0 B f 0 G 0 R 0 G , &Pi; 4 3 = P 0 B B f 0 D R 0 D ,
&Pi; 44 = d i a g - p 00 - 1 P 0 , - p 00 - 1 M 0 , - p 00 - 1 R 0 , &Pi; 51 = P 1 A 0 0 B f 0 C N 0 0 R 1 C 0 0 ,
&Pi; 52 = P 1 F 0 0 B f 0 G 0 R 1 G , &Pi; 53 = P 1 B B f 0 D R 1 D , &Pi; 55 = d i a g - p 01 - 1 P 1 , - p 01 - 1 M 1 , - p 01 - 1 R 1 ,
Π 61=[L-V 00],Π 66=-I, &Pi; 74 = 0 B f 0 T I , Π 75=Π 74,Π 77=-ε 2I,
Π 86=I,Π 88=-ε 3I, &Omega; 11 = d i a g - P 1 + &epsiv; 1 ( F 1 T F 1 + G 1 T G 1 ) , - M 1 , - R 1 , Ω 22=-ε 1I,Ω 33=-γ 2I,
&Omega; 41 = P 0 A 0 0 0 N 1 B f 1 0 0 R 0 , &Omega; 42 = P 0 F 0 0 0 0 0 , &Omega; 43 = P 0 B 0 0 ,
&Omega; 44 = d i a g - p 10 - 1 P 0 , - p 10 - 1 M 0 , - p 10 - 1 R 0 , &Omega; 51 = P 1 A 0 0 0 N 1 B f 1 0 0 R 1 , &Omega; 52 = P 1 F 0 0 0 0 0 , &Omega; 53 = P 1 B 0 0 ,
&Omega; 55 = d i a g - p 11 - 1 P 1 , - p 11 - 1 M 1 , - p 11 - 1 R 1 , Ω 61=[0-V 10],Ω 66=-I。
CN201510967000.3A 2015-12-18 2015-12-18 The H of network-based Lipschitz nonlinear system∞Filter information processing method Active CN105553442B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510967000.3A CN105553442B (en) 2015-12-18 2015-12-18 The H of network-based Lipschitz nonlinear system∞Filter information processing method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510967000.3A CN105553442B (en) 2015-12-18 2015-12-18 The H of network-based Lipschitz nonlinear system∞Filter information processing method

Publications (2)

Publication Number Publication Date
CN105553442A true CN105553442A (en) 2016-05-04
CN105553442B CN105553442B (en) 2018-12-25

Family

ID=55832416

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510967000.3A Active CN105553442B (en) 2015-12-18 2015-12-18 The H of network-based Lipschitz nonlinear system∞Filter information processing method

Country Status (1)

Country Link
CN (1) CN105553442B (en)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105988368A (en) * 2016-07-27 2016-10-05 江南大学 Fault-tolerant control method for networked control system with time-varying delay
CN106227700A (en) * 2016-08-15 2016-12-14 江南大学 A kind of uncatalyzed coking dissipation filtering method of network control system
CN108847828A (en) * 2018-07-23 2018-11-20 哈尔滨理工大学 A kind of non-linear event triggering filtering method with stochastic modeling error
CN109976160A (en) * 2019-04-23 2019-07-05 哈尔滨理工大学 A kind of H ∞ control system of unilateral Lipschitz nonlinear systems with delay
CN111897227A (en) * 2020-08-25 2020-11-06 南通大学 Non-fragile optimal control method for specified attenuation rate of Lipschitz nonlinear system
CN112015087A (en) * 2020-08-25 2020-12-01 南通大学 Non-fragile optimal control method of Lipschitz nonlinear system

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6574336B1 (en) * 1996-06-19 2003-06-03 Nokia Telecommunications Oy Echo suppressor and non-linear processor of echo canceller
CN103701433A (en) * 2013-12-29 2014-04-02 哈尔滨理工大学 Quantitative filtering method of time-varying target tracking system under condition with multiple measurement loss

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6574336B1 (en) * 1996-06-19 2003-06-03 Nokia Telecommunications Oy Echo suppressor and non-linear processor of echo canceller
CN103701433A (en) * 2013-12-29 2014-04-02 哈尔滨理工大学 Quantitative filtering method of time-varying target tracking system under condition with multiple measurement loss

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
楼斌: "基于Markov模型的无线网络量化系统的H∞滤波研究", 《中国优秀硕士学位论文全文数据库》 *

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105988368A (en) * 2016-07-27 2016-10-05 江南大学 Fault-tolerant control method for networked control system with time-varying delay
CN105988368B (en) * 2016-07-27 2019-05-07 江南大学 A kind of fault tolerant control method of the network control system with time-vary delay system
CN106227700A (en) * 2016-08-15 2016-12-14 江南大学 A kind of uncatalyzed coking dissipation filtering method of network control system
CN106227700B (en) * 2016-08-15 2019-01-08 江南大学 A kind of uncatalyzed coking dissipation filtering method of network control system
CN108847828A (en) * 2018-07-23 2018-11-20 哈尔滨理工大学 A kind of non-linear event triggering filtering method with stochastic modeling error
CN108847828B (en) * 2018-07-23 2022-02-11 哈尔滨理工大学 Nonlinear event trigger filtering method with random modeling error
CN109976160A (en) * 2019-04-23 2019-07-05 哈尔滨理工大学 A kind of H ∞ control system of unilateral Lipschitz nonlinear systems with delay
CN111897227A (en) * 2020-08-25 2020-11-06 南通大学 Non-fragile optimal control method for specified attenuation rate of Lipschitz nonlinear system
CN112015087A (en) * 2020-08-25 2020-12-01 南通大学 Non-fragile optimal control method of Lipschitz nonlinear system
CN112015087B (en) * 2020-08-25 2022-05-24 南通大学 Non-fragile optimal control method of Lipschitz nonlinear system
CN111897227B (en) * 2020-08-25 2022-05-24 南通大学 Non-fragile optimal control method for specified attenuation rate of Lipschitz nonlinear system

Also Published As

Publication number Publication date
CN105553442B (en) 2018-12-25

Similar Documents

Publication Publication Date Title
CN105553442A (en) Network-based H-infinity filter information processing method for Lipschitz nonlinear system
CN111629380B (en) Dynamic resource allocation method for high concurrency multi-service industrial 5G network
CN101232180B (en) Power distribution system load obscurity model building device and method
CN106529479A (en) Non-fragile dissipative filtering method of nonlinear networked control system
CN104168131A (en) Flow generation method of power dispatching exchange network based on multicast communication
CN106383442B (en) A kind of H of networking Linear Parameter-Varying Systems∞Control method
CN109558436B (en) Airport flight delay cause and effect relationship mining method based on transfer entropy
CN112327810B (en) Fault estimation method of dynamic event triggered transmission Markov jump system
CN106227700A (en) A kind of uncatalyzed coking dissipation filtering method of network control system
CN105606381A (en) Distributed filtering network system and design method
CN102651641B (en) L2-L infinity filtering information processing method based on logarithmic quantization for network control system
CN107124158A (en) Wireless sensor network filtering information processing system and method based on logarithmic quantization
Sekine et al. LACSLE: lightweight and adaptive compressed sensing based on deep learning for edge devices
CN102664843A (en) Method and device for single-frequency interference detection and elimination
CN113950018B (en) Asynchronous multi-sensor network system and global ellipsoid state estimation method
CN102299766B (en) Joint optimization method for dimensionality reduction and quantification of communication signal for object state estimation
CN104331630A (en) State estimation and data fusion method for multi-rate observation data
CN116520814B (en) Equipment fault prediction method and device based on federal learning under cloud edge cooperative architecture
CN106454892A (en) Set-valued Kalman filtering method based on quantification and events
CN111046697A (en) Adaptive modulation signal identification method based on fuzzy logic system
Wu et al. Short-message communication Lossy data compression algorithm for BeiDou-3 satellite information transmission
CN111030644B (en) Finite time dissipation filtering method of nonlinear networked control system
Guo et al. System identification with binary-valued output observations under either-or communication and data packet dropout
Chakravorty et al. Distortion-transmission trade-off in real-time transmission of Gauss-Markov sources
CN108629096B (en) Non-fragile quantization dissipation filtering method of nonlinear networked system

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right

Effective date of registration: 20201224

Address after: Room 504, building 6, Tianan Digital City, 99 Shennan Road, Gangzha District, Nantong City, Jiangsu Province, 226000

Patentee after: JIANGSU EMPOWER INTELLIGENT TECHNOLOGY Co.,Ltd.

Address before: 226019 Jiangsu Province, Nantong City Chongchuan District sik Road No. 9

Patentee before: NANTONG University

TR01 Transfer of patent right