CN103001567A - Multi-mode inverse model identification method for speed regulating system of six-phase permanent magnet synchronous linear motor - Google Patents
Multi-mode inverse model identification method for speed regulating system of six-phase permanent magnet synchronous linear motor Download PDFInfo
- Publication number
- CN103001567A CN103001567A CN2012104536057A CN201210453605A CN103001567A CN 103001567 A CN103001567 A CN 103001567A CN 2012104536057 A CN2012104536057 A CN 2012104536057A CN 201210453605 A CN201210453605 A CN 201210453605A CN 103001567 A CN103001567 A CN 103001567A
- Authority
- CN
- China
- Prior art keywords
- lssvm
- inverse
- phase permanent
- local
- output
- 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
Links
- 230000001360 synchronised effect Effects 0.000 title claims abstract description 34
- 238000000034 method Methods 0.000 title claims abstract description 20
- 230000001105 regulatory effect Effects 0.000 title abstract 4
- 238000012549 training Methods 0.000 claims description 6
- 238000002790 cross-validation Methods 0.000 claims description 3
- 238000001914 filtration Methods 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- 230000000737 periodic effect Effects 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 238000013507 mapping Methods 0.000 abstract description 3
- 238000012843 least square support vector machine Methods 0.000 abstract 2
- 230000001276 controlling effect Effects 0.000 abstract 1
- 238000013178 mathematical model Methods 0.000 abstract 1
- 230000001537 neural effect Effects 0.000 description 6
- 238000013528 artificial neural network Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000013459 approach Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000004907 flux Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000006978 adaptation Effects 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000011217 control strategy Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 230000002441 reversible effect Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
Images
Abstract
The invention discloses a multi-mode inverse model identification method for a speed regulating system of a six-phase permanent magnet synchronous linear motor. The multi-mode inverse model identification method includes that the speed regulating system of the six-phase permanent magnet synchronous linear motor is formed by process for controlling a d shaft component i<d> of primary current of the six-phase permanent magnet synchronous linear motor so that the i<d> is equal to zero, and the reversibility of a mathematical model of the speed regulating system is verified; and input space and output space of the system are divided according to a multi-mode idea, nonlinear mapping in local inverse models of the system is fit by a least square support vector machine (LSSVM), the various local inverse models are outputted in a weighted manner, and an offline initial inverse model is built. Weights of the local inverse models are adaptively adjusted by an improved recursive least square (RLS) algorithm according to deviation of output of the inverse model and input of the system, and the initial inverse model is corrected in an online manner, so that the initial inverse model can adapt to change of an object, the identification precision is improved, and the convergence rate is increased.
Description
Technical field
The present invention relates to a kind of multimode identification of inverse model method of six-phase permanent-magnet linear synchronous motor governing system, be applicable to the technical field of electric drive control.
Background technology
The multi-phase permanent linear synchronous motor has plurality of advantages as a kind of New-type electric machine: 1) high-power with the low-voltage device realization, be particularly suitable for obtaining high pressure but the powerful occasion of needs output; 2) significantly reduce motor harmonic loss, reduce motor torque ripple, improve system effectiveness and stability.Therefore the multi-phase permanent linear synchronous motor is widely used in logistics system, industrial equipment, information and automated system, traffic and the many-side such as civilian, military.
The multi-phase permanent linear synchronous motor is a multivariable close coupling non linear system, and owing to adopt direct drive mode, load disturbance, the disturbance of ripple thrust, frictional force disturbance and other uncertain disturbances can directly act on motor, have a strong impact on the control precision of motor.Traditional PID control strategy is difficult to realize the control requirement of its increasingly stringent.Neural network inverse control is as a kind of feedback linearization method of intelligence, the nonlinear function that combines the height of the clear physics conception of method of inverse, understandable characteristics and neural net approaches the advantage strong with adaptive ability, the decoupling zero of non linear system is controlled have preferably effect.Be applied to the multi-phase permanent linear synchronous motor, can realize the high accuracy control of motor.
The Neural network inverse control method is that the system's inversion model that builds is connected with original system, thereby the non linear system of complexity is approximately pseudo-linear system, and the design additional controller is realized its high performance control.Therefore, the structure of system's inversion model is most important, and increasing document is studied it, as: neural network generalized inverse model building method, but the limit of the pseudo-linear hybrid system of arbitrary disposition gained; The online adjustment strategy of system neural network inversion model is constructed first RBF neural net off-line training and is obtained initial parameter, adjusts online its network weight by the improved BP algorithm again, has strengthened the generalization ability of neural net.
Yet present research all is the static non linear mappings that approach system's inversion model with a neural net, considers that system condition is complicated, only is not enough to accurate descriptive system at the inversion model of different conditions with a network; And, although literature research the has been arranged on-line identification method of system's inversion model, because it needs to adjust in real time each weights and the threshold value of neural net inside, performance difficulty.
Summary of the invention
The present invention is take the six-phase permanent-magnet linear synchronous motor as object, purpose provides a kind of multimode identification of inverse model method that is applicable to multi-phase permanent linear synchronous motor governing system, so that adopt the multi-phase permanent linear synchronous motor of inverse system control to realize more high performance control.
The technical solution used in the present invention is: the multimode identification of inverse model method of six-phase permanent-magnet linear synchronous motor governing system comprises the steps:
1) adopts primary current d axle component i
dThe six-phase permanent-magnet linear synchronous motor governing system of=0 control mode operation, with square wave at random as input data u
qEncourage, and to input data u
q, output data ω
rSample; To the sampled data smothing filtering, ask for derivative and periodic sampling, obtain sample data
2) according to motor speed ω
rLow speed, middling speed and high velocity between, sample data is divided into three different subspaces; With sample data
As input, { u
qAs output, the LSSVM off-line training is carried out in each subspace; Select Gaussian function as the LSSVM kernel function, adopt cross-validation method commonly used to obtain suitable regularization parameter γ and nuclear parameter σ, obtain respectively three LSSVM;
3) LSSVM that each is trained adds respectively two integrators and consists of local LSSVM inverse system; To local LSSVM inverse system weights initialize, with addition after three local inverse system weightings, obtain the output of the initial inversion model of six-phase permanent-magnet linear synchronous motor governing system respectively
4) utilize variable forgetting factor function lambda (k) to improving in the RLS algorithm, the RLS algorithm is improved:
In the formula, y is the output of modified RLS algorithm; X=[x
1, L, x
m]
TM input for modified RLS algorithm; W=[w
1, L, w
m]
TWeights for m input; D is desired output; E is desired output and the actual deviation of exporting; K is gain; P is the inverse matrix of k; Int () is bracket function; ρ is responsive gain; λ
MinMinimum value for forgetting factor.When deviation e is tending towards infinite, λ=λ
Min, so that algorithm can comparatively fast be followed the tracks of the local trend of non-stationary signal; E is tending towards at 0 o'clock, and λ=1 is to reduce parameter estimating error.
5) local LSSVM inverse system is regarded respectively as the input x=[x of modified RLS algorithm
1, x
2, x
3]
T, the weights of corresponding local LSSVM inverse system are regarded respectively the weight w of modified RLS algorithm=[w as
1, w
2, w
3]
T, the output of six-phase permanent-magnet linear synchronous motor governing system inversion model
Be the output y of modified RLS algorithm;
6) according to six-phase permanent-magnet linear synchronous motor governing system input data u
qExport with six-phase permanent-magnet linear synchronous motor governing system inversion model
Deviation e, adjust online the weights of local LSSVM inverse system, so that obtain more accurately the output of six-phase permanent-magnet linear synchronous motor governing system inversion model after the weighting of local LSSVM inverse system
The invention has the beneficial effects as follows:
1) with the inner a plurality of weights of weights adjustment replacement neural net of LSSVM partial model and the adjustment of threshold value, greatly reduces the online difficulty of adjusting of inversion model, improved the real-time of modeling; And LSSVM identification precision and computational efficiency are higher, are easy to engineering and use.
2) improved RLS Weight number adaptively adjustment algorithm can be adjusted forgetting factor in real time according to the size of Identification Errors, taken into account Model Distinguish initial convergence speed, the time become follow-up control and identification precision.
3) multi-model LSSVM adopts and decomposes synthetic idea about modeling, only needs the input and output sample set of system, has preferably local fit precision and generalization ability, can obtain gratifying identification of inverse model effect.
Description of drawings
The LSSVM51 that it is good that Fig. 1 is off-line training adds the schematic diagram that two integrators consist of local LSSVM inverse system 61.
Fig. 2 is the schematic diagram of the initial inversion model of addition constructing system after local LSSVM inverse system 61, local LSSVM inverse system 62 and local LSSVM inverse system 63 weightings.
Fig. 3 is the multimode inversion model on-line identification schematic diagram of six-phase permanent-magnet linear synchronous motor governing system 1.
Embodiment
Below in conjunction with accompanying drawing the present invention is described in detail, implementation divides following 8 steps:
During 1, by secondary flux linkage orientation, the six-phase permanent-magnet linear synchronous motor adopts primary current d axle component i
d=0 control method, and combine with load and to consist of six-phase permanent-magnet linear synchronous motor governing system 1.System simplification is the non linear system of the single output of single input, and input variable is primary voltage q axle component u
q2, output variable is motor speed ω
r4.Select primary current q axle component i
q3 and motor speed ω
r4 is quantity of state, and the Mathematical Modeling that obtains is the second order state equation.
2, adopt the Interactor algorithm to output variable ω
r4 ask local derviation, until the aobvious input variable u that contains
q2.According to the inverse system principle as can be known, the inverse system of six-phase permanent-magnet linear synchronous motor governing system 1 correspondence exists, and can determine that the input variable of its inverse system is motor speed ω
r4 second dervative
Output variable is primary voltage q axle component u
q2, and expression formula is
3, six-phase permanent-magnet linear synchronous motor governing system 1 is adopted primary current d axle component i
d=0 control mode operation, with the at random square wave of realistic range of operation as input variable u
q2 come excitation system, and to input data u
q2, output data ω
r4 sample.To the sampled data smothing filtering, ask for derivative and periodic sampling, obtain sample data
4, according to motor speed ω
rBetween 4 low speed, middling speed and high velocity, sample data is divided into three different subspaces.With
As input, { u
qAs output, the LSSVM off-line training is carried out in each subspace.Select Gaussian function as the LSSVM kernel function, adopt cross-validation method commonly used to obtain suitable regularization parameter γ and nuclear parameter σ, obtain respectively LSSVM51, LSSVM52 and LSSVM53.
5, the LSSVM51 that off-line training is good, LSSVM52 and LSSVM53 add two integrators, consist of respectively local LSSVM inverse system 61(as shown in Figure 1), local LSSVM inverse system 62 and local LSSVM inverse system 63.The structure of the initial inversion model of system as shown in Figure 2, respectively to the weight w of local LSSVM inverse system 61, local LSSVM inverse system 62 and local LSSVM inverse system 63
171, w
272 and w
373 initializes with addition after three local inverse system weightings, obtain the output of the initial inversion model of system
6, utilize variable forgetting factor function lambda (k) to improving in the RLS algorithm, the mathematic(al) representation of the RLS algorithm that is improved:
In the formula, y is the output of modified RLS algorithm; X=[x
1, L, x
m]
TM input for modified RLS algorithm; W=[w
1, L, w
m]
TWeights for m input; D is desired output; E is desired output and the actual deviation of exporting; K is gain; P is the inverse matrix of k; Int () is bracket function; ρ is responsive gain; λ
MinMinimum value for forgetting factor.When deviation e is tending towards infinite, λ=λ
Min, so that algorithm can comparatively fast be followed the tracks of the local trend of non-stationary signal; E is tending towards at 0 o'clock, and λ=1 is to reduce parameter estimating error.
7, the input x=[x that local LSSVM inverse system 61, local LSSVM inverse system 62 and local LSSVM inverse system 63 is regarded respectively as modified RLS algorithm
1, x
2, x
3]
T, corresponding weight w
171, w
272 and w
373 weight w of regarding respectively modified RLS algorithm as=[w
1, w
2, w
3]
T, the output of system's inversion model
Be the output y of modified RLS algorithm.
8, as shown in Figure 3, according to the input u of system
q2 with the output of system inversion model
Deviation e9, adjust online weight w
171, w
272 and w
373, so that obtain more accurately system's inversion model output after local LSSVM inverse system 61, local LSSVM inverse system 62 and local LSSVM inverse system 63 weightings
When the present invention pressed secondary flux linkage orientation, the six-phase permanent-magnet linear synchronous motor adopted primary current d axle component i
d=0 control method consists of six-phase permanent-magnet linear synchronous motor governing system, and checking governing system Mathematical Modeling is reversible.According to multi-model thought, system's input-output space is divided, utilize the Nonlinear Mapping in the local inversion model of least square method supporting vector machine (LSSVM) match system, with each local inversion model weighting output, set up the initial inversion model of off-line.According to the deviation of inversion model output with system's input, adjust local inversion model weights by improved recurrence least square (RLS) algorithm self adaptation, to initial inversion model on-line correction, make the variation of its adaption object, improve identification precision and convergence rate.
Claims (1)
1. the multimode identification of inverse model method of a six-phase permanent-magnet linear synchronous motor governing system is characterized in that comprising the steps:
1) adopts primary current d axle component i
d=0 control mode operation six-phase permanent-magnet linear synchronous motor governing system, with square wave at random as input data u
qEncourage, and to input data u
q, output data ω
rSample; To the sampled data smothing filtering, ask for derivative and periodic sampling, obtain sample data
2) according to motor speed ω
rLow speed, middling speed and high velocity between, sample data is divided into three different subspaces; With sample data
As input, { u
qAs output, the LSSVM off-line training is carried out in each subspace; Select Gaussian function as the LSSVM kernel function, adopt cross-validation method commonly used to obtain suitable regularization parameter γ and nuclear parameter σ, obtain respectively three LSSVM that train;
3) LSSVM that each is trained adds respectively two integrators and consists of local LSSVM inverse system; To local LSSVM inverse system weights initialize, with addition after three local inverse system weightings, obtain the output of the initial inversion model of six-phase permanent-magnet linear synchronous motor governing system respectively
4) utilize variable forgetting factor function lambda (k) to improving in the RLS algorithm, the RLS algorithm is improved:
In the formula, y is the output of modified RLS algorithm; X=[x
1, x
2, x
3]
TM input for modified RLS algorithm; W=[w
1, w
2, w
3]
TM the weights of inputting for modified RLS algorithm; D is desired output; E is desired output and the actual deviation of exporting; K is gain; P is the inverse matrix of k; Int () is bracket function; ρ is responsive gain; λ
MinMinimum value for forgetting factor: when deviation e is tending towards infinite, λ=λ
MinE is tending towards at 0 o'clock, λ=1;
5) local LSSVM inverse system is regarded respectively as the input x=[x of modified RLS algorithm
1, x
2, x
3]
T, the weights of corresponding local LSSVM inverse system are regarded respectively the weight w of modified RLS algorithm=[w as
1, w
2, w
3]
T, the output of six-phase permanent-magnet linear synchronous motor governing system inversion model
Be the output y of modified RLS algorithm;
6) according to six-phase permanent-magnet linear synchronous motor governing system input data u
qExport with six-phase permanent-magnet linear synchronous motor governing system inversion model
Deviation e, adjust online the weights of local LSSVM inverse system, so that obtain more accurately the output of six-phase permanent-magnet linear synchronous motor governing system inversion model after the weighting of local LSSVM inverse system
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210453605.7A CN103001567B (en) | 2012-11-13 | 2012-11-13 | Multi-mode inverse model identification method for speed regulating system of six-phase permanent magnet synchronous linear motor |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210453605.7A CN103001567B (en) | 2012-11-13 | 2012-11-13 | Multi-mode inverse model identification method for speed regulating system of six-phase permanent magnet synchronous linear motor |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103001567A true CN103001567A (en) | 2013-03-27 |
CN103001567B CN103001567B (en) | 2015-06-17 |
Family
ID=47929795
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210453605.7A Expired - Fee Related CN103001567B (en) | 2012-11-13 | 2012-11-13 | Multi-mode inverse model identification method for speed regulating system of six-phase permanent magnet synchronous linear motor |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103001567B (en) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105425587A (en) * | 2015-11-16 | 2016-03-23 | 北京理工大学 | Hysteresis nonlinear motor identification and control method |
CN106354695A (en) * | 2016-08-22 | 2017-01-25 | 北京理工大学 | Output-only linear time-varying structure modal parameter identification method |
CN108233817A (en) * | 2018-01-15 | 2018-06-29 | 中国人民解放军海军工程大学 | A kind of six phase line inductance electromotor Energy Chain handover control systems and method |
CN111564995A (en) * | 2020-05-25 | 2020-08-21 | 华中科技大学 | Linear oscillation motor control method based on self-adaptive full-order displacement observer |
CN112152529A (en) * | 2020-09-28 | 2020-12-29 | 长沙贝士德电气科技有限公司 | Maximum thrust control method and system for permanent magnet linear motor |
CN112415392A (en) * | 2020-11-03 | 2021-02-26 | 珠海格力电器股份有限公司 | Method for determining forgetting factor, electronic equipment, storage medium and device |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102769424A (en) * | 2012-06-18 | 2012-11-07 | 江苏大学 | Support vector machine inverse-based bearingless asynchronous motor control system |
-
2012
- 2012-11-13 CN CN201210453605.7A patent/CN103001567B/en not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102769424A (en) * | 2012-06-18 | 2012-11-07 | 江苏大学 | Support vector machine inverse-based bearingless asynchronous motor control system |
Non-Patent Citations (3)
Title |
---|
CHAOXU MU ET AL: "Inverse System Identification of Nonlinear Systems Using Least Square Support Vector Machine Based on FCM Clustering", 《NEURAL NETWORKS,2008.IJCNN 2008.(IEEE WORLD CONGRESS ON COMPUTATIONAL INTELLIGENCE)》 * |
张懿: "两电机变频调速系统的最小二乘支持向量机逆控制", 《中国博士学位论文全文数据库》 * |
沈曙光等: "最小支持向量机在系统逆动力学辨识与控制中的应用", 《中国电机工程学报》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105425587A (en) * | 2015-11-16 | 2016-03-23 | 北京理工大学 | Hysteresis nonlinear motor identification and control method |
CN106354695A (en) * | 2016-08-22 | 2017-01-25 | 北京理工大学 | Output-only linear time-varying structure modal parameter identification method |
CN106354695B (en) * | 2016-08-22 | 2019-09-17 | 北京理工大学 | One kind only exporting linear Time variable structure Modal Parameters Identification |
CN108233817A (en) * | 2018-01-15 | 2018-06-29 | 中国人民解放军海军工程大学 | A kind of six phase line inductance electromotor Energy Chain handover control systems and method |
CN108233817B (en) * | 2018-01-15 | 2020-04-17 | 中国人民解放军海军工程大学 | Six-phase linear induction motor energy chain switching control system and method |
CN111564995A (en) * | 2020-05-25 | 2020-08-21 | 华中科技大学 | Linear oscillation motor control method based on self-adaptive full-order displacement observer |
CN112152529A (en) * | 2020-09-28 | 2020-12-29 | 长沙贝士德电气科技有限公司 | Maximum thrust control method and system for permanent magnet linear motor |
CN112152529B (en) * | 2020-09-28 | 2022-08-12 | 长沙贝士德电气科技有限公司 | Maximum thrust control method and system for permanent magnet linear motor |
CN112415392A (en) * | 2020-11-03 | 2021-02-26 | 珠海格力电器股份有限公司 | Method for determining forgetting factor, electronic equipment, storage medium and device |
Also Published As
Publication number | Publication date |
---|---|
CN103001567B (en) | 2015-06-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103001567B (en) | Multi-mode inverse model identification method for speed regulating system of six-phase permanent magnet synchronous linear motor | |
CN102497156B (en) | Neural-network self-correcting control method of permanent magnet synchronous motor speed loop | |
CN104595106B (en) | Wind-power generating variable pitch control method based on intensified learning compensation | |
CN102662323B (en) | Adoptive sliding mode control method and adoptive sliding mode control system of wind power generation variable-pitch actuator | |
CN107179682A (en) | A kind of passive type load simulator and Surplus Moment suppressing method | |
CN105888971B (en) | A kind of large scale wind power machine blade active load shedding control system and method | |
CN102411305A (en) | Design method of comprehensive disturbance rejection control system for single-rotor wing helicopter/turboshaft engine | |
CN107203137A (en) | The non-linear heuristic Adaptive PID Control method of pump-storage generator governing system gain | |
CN103557037A (en) | Method for controlling rotating speed of steam turbine on basis of self-adaptive inverse control | |
CN104199294A (en) | Motor servo system bilateral neural network friction compensation and limited time coordination control method | |
CN111273544B (en) | Radar pitching motion control method based on prediction RBF feedforward compensation type fuzzy PID | |
CN107800343A (en) | The design method of asynchronous machine automatic disturbance rejection controller | |
CN102664580B (en) | Mixed smart control method of power system multi-stabilizer | |
CN104135205A (en) | Control method for maximum torque current rate of induction motor | |
CN112523944A (en) | Self-adaptive dynamic surface control method for variable pitch system of wind driven generator | |
US9154069B2 (en) | Integrated generation control system | |
CN102269971A (en) | Self-adapting servo controller based on model tracking | |
Wang et al. | The identification of electric load simulator for gun control systems based on variable-structure WNN with adaptive differential evolution | |
CN113659895A (en) | Permanent magnet synchronous motor all-state constraint finite time control method based on instruction filtering | |
CN110323741B (en) | Low-frequency oscillation analysis and suppression method for power system based on wave recording curve | |
CN109884881A (en) | A kind of design for surely taking aim at servo controller based on nonlinear PID controller technology | |
Sarkar et al. | Real time swing up and stabilization of rotary inverted pendulum system | |
CN107959453A (en) | A kind of improved MRAS speed observation procedure | |
CN102969948A (en) | Two-motor variable-frequency speed control system considering efficiency optimization and decoupling control method thereof | |
CN105099263A (en) | Ultrasonic motor position and speed control system for generalized regression neural network-based identifier |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20170619 Address after: 226681, Haian, Jiangsu County, Nantong Province, Mo Zhen, Yan Macun Patentee after: JIANGSU FEISI ALUMINIUM INDUSTRY CO.,LTD. Address before: 212003 Zhenjiang City, Jiangsu province dream Creek Road, No. 2 Patentee before: Jiangsu University of Science and Technology |
|
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150617 |