CN101545778B - Initial alignment method of electric control compass based on nonlinear filtering - Google Patents
Initial alignment method of electric control compass based on nonlinear filtering Download PDFInfo
- Publication number
- CN101545778B CN101545778B CN2008102368756A CN200810236875A CN101545778B CN 101545778 B CN101545778 B CN 101545778B CN 2008102368756 A CN2008102368756 A CN 2008102368756A CN 200810236875 A CN200810236875 A CN 200810236875A CN 101545778 B CN101545778 B CN 101545778B
- Authority
- CN
- China
- Prior art keywords
- gamma
- clj
- error
- delta
- initial alignment
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Gyroscopes (AREA)
Abstract
The present invention provides an initial alignment method of an electric control compass based on nonlinear filtering, which relates to the technical problems of the initial alignment of the electric control compass, and the like. A UKF nonlinear error model for initial alignment of the electric control compass is established according to the characteristics of the gyroscope and the acceleration in the electric control compass. The method comprises the following steps: establishing a nonlinear initial alignment model of the electric control compass; and then carrying out the UKF filtering initial alignment of the electric control compass. An indirect output correcting filter structure is adopted by the UKF filtering of the electric control compass. The UKF optimal estimation state equation adopted in the electric control compass is x(k+1) = xk + g (xk) delta t + B omegak, and the observation noise equation is zk = h(xk) + vk. The initial alignment method can effectively solve the nonlinear problem of the large misalignment angle in initial alignment, thereby shortening the time of initial alignment and improving the performance of the whole system.
Description
Technical field
The present invention relates to the technical matterss such as initial alignment in the equipment (following general designation " gyroscopic navigation equipment ") such as electromagnet control gyrocompass, platform compass, inertial navigation system, specifically be meant a kind of electromagnet control gyrocompass initial alignment method based on nonlinear filtering.
Background technology
Gyroscopic navigation equipment at first is in and prepares and the initial alignment state before startup enters operate as normal.Initial alignment is to set up the necessary starting condition of navigation equipment, as the initial attitude angle error of initial position, initial velocity, platform etc.The speed and the position of naval vessels obtain through integral and calculating, and this will provide starting condition V
X0, V
Y0,
λ
0Because start at harbour, initial V
X0, V
Y0Be zero,
λ
0Can have known measurement point to provide, perhaps provide with radio position-finder etc., but precision meet the demands should be within tens meters.Initial alignment can utilize the responsive error angle signal of the inertance element of gyroscopic navigation equipment own that system is calibrated, and designs a rational initial alignment scheme, and reduce error to satisfy request for utilization be very necessary as far as possible.The time of alignment precision and aligning is two important technology indexs when carrying out initial alignment.The initial alignment precision influences the performance of inertial navigation system, and the initial alignment time tag quick-reaction capability (QRC), therefore requires initial alignment precision height, aligning time weak point, and is promptly smart and fast.In order to reach this requirement, the initial alignment error model and the corresponding initial alignment Filtering Estimation method of realistic system be arranged.
Always there is nonlinear with different degree in real system, some system can be similar to and regard linear system as, but most systems then can not only be described with linear differential equation, nonlinear factor wherein can not be ignored, or for analysis integrated result better, must use the nonlinear mathematical model of reflection real system, the stochastic non-linear system filtering problem is arisen at the historic moment.
Have a large amount of random noises inevitably in the actual gyroscopic navigation equipment, traditional filtering method such as the high pass in the frequency domain, low pass, band are logical is powerless to this.The scholar had proposed to compose the Wiener filtering method of eliminating random noise by rated output afterwards.Though the Wiener filtering method is effective to random white noise, the calculated amount of this method is big, can't be suitable in real time complicated slightly system.Nineteen sixty, the Kalman has proposed the Kalman filtering algorithm based on state-space method and matrix operation, it provides strong tool for the optimal filtering of time domain recursion, it adopts time domain recursive algorithm, must storage time measurement amount in the process, thereby Kalman filtering is used widely in every field such as navigation, communication, remote sensing.But Kalman filtering can only be handled linear system.For solving the nonlinear problem in the engineering reality, the scholar had proposed improved Kalman filtering algorithm EKF (Extended KalmanFilter EKF) again afterwards, but existed not enough.Up to the present, though EKF is had numerous improving one's methods, block EKF, iteration EKF etc. as high-order.But the linearization error of EKF has reduced the accuracy of model, prolongation in time, precision is difficult to guarantee, sometimes even disperse, these defectives still are difficult to overcome, so the various countries navigation expert of circle is making great efforts to study new more effective nonlinear noise reduction processing method always.For this reason, utilize UKF method, open up new technological approaches for solving the navigation field nonlinear problem based on the statistics nonlinear transformation.
Summary of the invention
The objective of the invention is in order to overcome the deficiency of above-mentioned background, propose a kind of electromagnet control gyrocompass initial alignment method, improve initial alignment precision and the practical operation effect of rapid alignment time to reach based on nonlinear filtering.
At above-mentioned purpose, method proposed by the invention is set up the UKF nonlinearity erron model of electromagnet control gyrocompass initial alignment according to the characteristic of gyro in the electromagnet control gyrocompass and acceleration, obtain the optimal estimation of system's random noise effectively, guarantee the precision of compensation, shorten start-up time, its concrete grammar is:
1) foundation of the non-linear initial alignment model of electromagnet control gyrocompass
Electromagnet control gyrocompass adopts the fixing northern formula inertial navigation machinery layout that refers to.Complete error model will be considered the error model and the acceleration error model of gyro.Gyro error includes three calibration factor errors, the alignment errors on three axles on the axle; Contain three misalignment errors of three calibration factor sum of errors on the axle in the same accelerometer error.INS errors generally has: 3 platform misalignments (α, beta, gamma); 3 velocity error (δ V
x, δ V
y, δ V
z); Longitude, latitude, 3 site errors of height; Consider that gyro, all error model factors of accelerometer can reach 21 dimension state variables.21 dimension error states do not allow on computing velocity and complexity for electromagnet control gyrocompass.The inventor drops to 15 dimensions through dimension-reduction treatment by 21 dimensions, and filtering can reduce more than 40% computing time like this, and precision property is suitable.So, be necessary the error model of electromagnet control gyrocompass is simplified processing.Because of electromagnet control gyrocompass is a shipborne equipment, vertically to acceleration, speed, height can be ignored, and only need to calculate east orientation and north orientation acceleration error, and three-dimensional speed and position also are simplified to two dimensional surface.The initial alignment time is not long, and gyro error and North-East Bound acceleration all are reduced to white noise at random.According to above-mentioned, the error state of electromagnet control gyrocompass is defined as the vector of one 5 dimension: x=[α β γ δ V
xδ V
y]
TAccording to physical construction layout in the electromagnet control gyrocompass, can obtain α, β, γ, the δ V of system alignment
x, δ V
yThe nonlinearity erron model:
Following formula is rewritten as following form:
In the formula: x=[α β γ δ V
xδ V
y]
T, w=[w
xw
yw
zw
δ Vxw
δ Vy]
TIt is the random white noise of gyro and accelerometer; B=I
5 * 5R wherein
eIt is earth radius; ω
IeIt is earth rotational angular velocity; α, β, γ are that the platform error angle is three projections on the axle at navigation coordinate; δ V
x, δ V
yBe respectively latitude and North-East Bound velocity error; ε
x, ε
y, ε
zBe the projection of gyroscope constant value error in local geographic coordinate system oxyz,
Be that horizontal accelerometer often is worth biasing.
G () is non-linear, and it is carried out discretize, and establishing sampling interval is Δ t, and X (t+ Δ t) is located to carry out Taylor expansion at X (t), obtains:
Make X (t+ Δ t)=x
K+1, X (t)=x
k, ω (t)=ω
kAccording to accuracy requirement, ignore second order and the above a small amount of of second order of launching progression, the discrete equation that obtains system is:
x
k+1=x
k+g(x
k)Δt+Bω
k (4)
Formula (4) is the non-linear initial alignment model of electromagnet control gyrocompass.
2) the UKF filtering initial alignment of electromagnet control gyrocompass
● state equation:
Indirect output calibration filter construction, its structural drawing such as Fig. 4 are adopted in the UKF filtering of electromagnet control gyrocompass.The UKF optimal estimation state equation that adopts in the electromagnet control gyrocompass is the depression of order error model of above-mentioned electromagnet control gyrocompass, sees formula (4).
● the error observation equation
Because the terrestrial magnetic field very easily is subjected on every side ferromagnetic magnetic interference over the ground, thereby Magnetic Sensor can produce than mistake, so precision is not high; And the geomagnetic observation of Magnetic Sensor is also to be a kind of autonomous observation, external information is relied on lack, and does not also need externally to carry out any communication interaction simultaneously, thereby also is used widely in the combination of navigational system in recent years.In the compass of being discussed, we have installed additional and the Magnetic Sensor measurement mechanism, the amount that can observe has course and cradle angle in length and breadth, the volume of Magnetic Sensor is little, as long as it is close with SINS as far as possible when installing, thereby the armed lever effect between Magnetic Sensor and SINS is very little, with respect to the bigger observational error of Magnetic Sensor, can ignore.
If the attitude of Magnetic Sensor output is respectively pitching R
Clj, rolling P
Clj, course H
Clj, the attitude angle of compass output is R
INS, P
INS, H
INSSo can be similar to and think and obtain being observed of attitude error angle
α=R
clj-R
INS
β=P
clj-P
INS (5)
γ=H
clj-H
INS
Formula (5) is out of shape, and makes Z=[R
CljP
CljH
Clj]
T, have
W wherein
CljFor the error observation noise of Magnetic Sensor, can think the zero-mean white noise in the short time of aligning.To measure the equation brief note is
z
k=h(x
k)+v
k (7)
V wherein
kBe observation noise.
● the UKF filtering algorithm of discrete recursion
At the error equation of the top electromagnet control gyrocompass of setting up and the system of observation noise equation (4) and (7) definition, be rewritten as following formula:
x(k+1)=f[x(k),u(k),w(k)] (8)
z(k)=h[x(k),v(k)] (9)
In UKF filtering, owing to have noise item, need expand dimension to state and handle, make x
a=[x
Tw
Tv
T]
T, concrete UKF filtering flow process is as follows:
The state starting condition is
For k ∈ 1,2 .... ∞, calculate the sigma point rebuild:
The predictive equation of UKF is:
The renewal equation of UKF is:
X wherein
a=[x
Tw
Tv
T]
T, χ
a=[(χ
x)
T(χ
w)
T(χ
v)
T]
T, the ξ list of values sigma point sampling method in the formula (11), the i.e. method of UT conversion.Weighting coefficient W wherein
iCalculating formula be
Wherein
τ is that a little positive number (is taken as 1>τ 〉=1e usually
-4), determined the distribution of sigma point around x.λ=τ
2(L+k)-L is a scaling factor.K also is a scaling factor, often is taken as 0 or 3-L, and η is that η gets to the added parameter of the distribution of needle-like attitude vector x,
Be the i row of matrix square root.
Superiority of the present invention is:
Electromagnet control gyrocompass initial alignment test findings shows, utilize error model that system derives and observation equation to carry out UKF filtering initial alignment after, when reaching same accuracy requirement, the time of alignment of orientation is shortened in the 200s by original 2879s; Simultaneously the longitudinal and transverse misalignment that shakes also under same precision the initial alignment time go up original 2400s respectively and 2700s shortens to 220s and 180s; In addition, the speed of system is dispersed also with the accumulation of acceleration stochastic error and is improved greatly, and in 1500s, electromagnet control gyrocompass autoregistration speed diffuses to 14m/s as seen from Figure 5, and after adopting UKF filtering initial alignment, speed only diffuses to about 4m/s.
Analyze α, β, γ and speed (the North-East Bound velocity and) curve of its initial alignment testing experiment, can obtain to draw a conclusion:
(1) the initial alignment nonlinearity erron model of electromagnet control gyrocompass system derivation and UKF combination initial alignment algorithm can effectively solve non-linear that big misalignment exists in the initial alignment, thereby shorten the initial alignment time, and the total system performance is improved;
(2) electromagnet control gyrocompass adopts UKF initial alignment technology, can effectively suppress conventional static initial alignment down speed disperse phenomenon;
The initial alignment nonlinearity erron model that electromagnet control gyrocompass adopts, UKF combination initial alignment algorithm has important reference for other engineering practical applications such as optimal estimation of the navigational state error of big nonlinear problem that exists in other system and navigational system.
Description of drawings
Fig. 1 is the electromagnet control gyrocompass combined filter align structures block diagram based on UKF
Fig. 2 compares from the γ of initial alignment and UKF combination initial alignment in the embodiment
Fig. 3 compares from the β of initial alignment and UKF combination initial alignment in the embodiment
Fig. 4 compares from the α of initial alignment and UKF combination initial alignment in the embodiment
Fig. 5 be in the embodiment from the velocity ratio of initial alignment and UKF combination initial alignment
Embodiment
Describe performance of the present invention in detail below in conjunction with accompanying drawing, but they do not constitute the qualification to patent of the present invention, only do for example.Simultaneously by illustrating that advantage of the present invention will become clear more and understanding easily.
Embodiments of the invention:
According to the state equation and the observation equation of electromagnet control gyrocompass, according to discrete stepping type (9)-(12) of UKF, system carries out initial alignment again, and the scheme flow process of its realization is seen accompanying drawing 1.In system's initial alignment test, initial value is chosen:
Angular unit is: °, speed unit is: m/s; Test is carried out in Chongqing factory, and its latitude is 30 °.By the error statistics characteristic of electromagnet control gyrocompass, calculate P according to the method for introducing in " Kalman filtering and integrated navigation principle " book
0, getting each state-noise variance intensity more respectively is 5% of initial value.Given initial value
And P
0, electromagnet control gyrocompass utilizes C/C++ in Control Software aforementioned error model and UKF to be realized, τ value 0.5, and η gets 2, and the filtering cycle is 0.1s, is respectively 20Hz and 10Hz to the sampling rate of electromagnet control gyrocompass and Magnetic Sensor.
In the static initial alignment test that the electromagnet control gyrocompass system carries out, the installation position benchmark of system is aimed at by optic theodolite optical and is guaranteed to aim at north orientation.System is carried out independent autoregistration and UKF combined filter aim at and compare, can obtain the result of accompanying drawing 2,3,4,5.In the experimental test, often be worth and be corrected after alignment error is measured by measuring method.In the test result of the autoregistration orientation misalignment γ of system and the UKF combined filter initial alignment of electromagnet control gyrocompass, orientation misalignment γ test result is seen Fig. 2.
Claims (1)
1. electromagnet control gyrocompass initial alignment method based on nonlinear filtering, it is characterized in that: the non-linear initial alignment model of electromagnet control gyrocompass, at first, set up electromagnet control gyrocompass and adopt the fixing northern formula inertial navigation machinery layout that refers to, obtain comprising the error model of 21 dimension state variables of gyro error, accelerometer error, platform misalignment, velocity error, all factors of site error, carry out dimension-reduction treatment to 15 dimension then, get α, β, γ, δ V
x, δ V
yAs the error state variable, the nonlinearity erron model that obtains system alignment is as follows:
Following formula is rewritten as following form:
In the formula: x=[α β γ δ V
xδ V
y]
T, w=[w
xw
yw
zw
δ Vxw
δ Vy]
TIt is the random white noise of gyro and accelerometer; B=I
5 * 5R wherein
eIt is earth radius; ω
IeIt is earth rotational angular velocity; α, β, γ are that the platform error angle is three projections on the axle at navigation coordinate; δ V
x, δ V
yBe respectively latitude and North-East Bound velocity error; ε
x, ε
y, ε
zBe the projection of gyroscope constant value error in local geographic coordinate system oxyz,
Be that horizontal accelerometer often is worth biasing;
G () is non-linear, and it is carried out discretize, and establishing sampling interval is Δ t, and X (t+ Δ t) is located to carry out Taylor expansion at X (t), obtains
Make X (t+ Δ t)=x
K+1, X (t)=x
k, ω (t)=ω
k, according to accuracy requirement, ignore second order and the above a small amount of of second order of launching progression, the discrete equation that obtains system is:
x
k+1=x
k+g(x
k)Δt+Bω
k (4)
Formula (4) is the error model of electromagnet control gyrocompass;
Build UKF filter state equation, the output calibration filter construction is adopted in electromagnet control gyrocompass UKF filtering, and state equation is the nonlinearity erron model equation (4) of electromagnet control gyrocompass;
Build UKF filtering measurement equation, the attitude of establishing Magnetic Sensor output is respectively pitching R
Clj, rolling P
Clj, course H
Clj, the attitude angle of compass output is R
INS, P
INS, H
INSSo, can be similar to and think and obtain being observed of attitude error angle
α=R
clj-R
INS
β=P
clj-P
INS (5)
γ=H
clj-H
INS
Formula (5) is out of shape, and makes z=[R
CljP
CljH
Clj]
T, have
W wherein
CljFor the error observation noise of Magnetic Sensor, can think the zero-mean white noise in the short time of aligning; To measure the equation brief note is
z
k=h(x
k)+v
k (7)
V wherein
kBe observation noise.
Press the error model equation of electromagnet control gyrocompass and the system of observation noise equation (4) and (7) definition, be rewritten as following formula:
x(k+1)=f[x(k),u(k),w(k)] (8)
z(k)=h[x(k),v(k)] (9)
In UKF filtering, owing to have noise item, need expand dimension to state and handle, make x
a=[x
Tw
Tv
T]
T, concrete UKF filtering flow process is as follows:
The state starting condition is
For k ∈ 1,2 .... ∞, calculate the sigma point rebuild:
The predictive equation of UKF is:
(12)
The renewal equation of UKF is:
(13)
X wherein
a=[x
Tw
Tv
T]
T, χ
a=[(χ
x)
T(χ
w)
T(χ
v)
T]
T, the ζ list of values sigma point sampling method in the formula (11), the i.e. method of UT conversion; Weighting coefficient W wherein
iCalculating formula be
(14)
Wherein
τ is that a little positive number (is taken as 1>τ 〉=1e usually
-4), determined the sigma point to exist
Distribution on every side, λ=τ
2(L+k)-and L is a scaling factor, and k also is a scaling factor, often is taken as 0 or 3-L, and η is that η gets 2 to the added parameter of the distribution of needle-like attitude vector x,
Be the i row of matrix square root.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008102368756A CN101545778B (en) | 2008-12-17 | 2008-12-17 | Initial alignment method of electric control compass based on nonlinear filtering |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008102368756A CN101545778B (en) | 2008-12-17 | 2008-12-17 | Initial alignment method of electric control compass based on nonlinear filtering |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101545778A CN101545778A (en) | 2009-09-30 |
CN101545778B true CN101545778B (en) | 2011-11-16 |
Family
ID=41193048
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2008102368756A Expired - Fee Related CN101545778B (en) | 2008-12-17 | 2008-12-17 | Initial alignment method of electric control compass based on nonlinear filtering |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101545778B (en) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102175267B (en) * | 2011-03-04 | 2013-05-08 | 中国人民解放军第二炮兵工程学院 | High-precision compensation method for horizontal angle of electro-optic theodolite |
CN102305636B (en) * | 2011-08-18 | 2013-02-06 | 江苏科技大学 | Rapid alignment method based on nonlinear initial alignment model |
CN102495831B (en) * | 2011-11-17 | 2014-08-06 | 西北工业大学 | Quaternion Hermitian approximate output method based on angular velocities for aircraft during extreme flight |
CN104103902B (en) * | 2014-07-23 | 2016-03-16 | 武汉虹信通信技术有限责任公司 | Based on the alignment methods of the point-to-point of compass and gradiometer |
CN106597498B (en) * | 2017-01-18 | 2020-04-24 | 哈尔滨工业大学 | Space-time deviation joint calibration method for multi-sensor fusion system |
CN111366143A (en) * | 2018-12-26 | 2020-07-03 | 北京信息科技大学 | Combined polar region compass device capable of automatically positioning and orienting |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101183004A (en) * | 2007-12-03 | 2008-05-21 | 哈尔滨工程大学 | Method for online real-time removing oscillation error of optical fibre gyroscope SINS system |
CN101246023A (en) * | 2008-03-21 | 2008-08-20 | 哈尔滨工程大学 | Closed-loop calibration method of micro-mechanical gyroscope inertial measuring component |
-
2008
- 2008-12-17 CN CN2008102368756A patent/CN101545778B/en not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101183004A (en) * | 2007-12-03 | 2008-05-21 | 哈尔滨工程大学 | Method for online real-time removing oscillation error of optical fibre gyroscope SINS system |
CN101246023A (en) * | 2008-03-21 | 2008-08-20 | 哈尔滨工程大学 | Closed-loop calibration method of micro-mechanical gyroscope inertial measuring component |
Also Published As
Publication number | Publication date |
---|---|
CN101545778A (en) | 2009-09-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107525503B (en) | Adaptive cascade Kalman filtering method based on combination of dual-antenna GPS and MIMU | |
CN103575299B (en) | Utilize dual-axis rotation inertial navigation system alignment and the error correcting method of External Observation information | |
CN103900565B (en) | A kind of inertial navigation system attitude acquisition method based on differential GPS | |
Wu et al. | Velocity/position integration formula part I: Application to in-flight coarse alignment | |
CN101706284B (en) | Method for increasing position precision of optical fiber gyro strap-down inertial navigation system used by ship | |
CN102486377B (en) | Method for acquiring initial course attitude of fiber optic gyro strapdown inertial navigation system | |
CN101476894B (en) | Vehicle-mounted SINS/GPS combined navigation system performance reinforcement method | |
CN103323026B (en) | The attitude reference estimation of deviation of star sensor and useful load and modification method | |
CN104344837B (en) | Speed observation-based redundant inertial navigation system accelerometer system level calibration method | |
CN102169184B (en) | Method and device for measuring installation misalignment angle of double-antenna GPS (Global Position System) in integrated navigation system | |
CN101949703B (en) | Strapdown inertial/satellite combined navigation filtering method | |
CN101706287B (en) | Rotating strapdown system on-site proving method based on digital high-passing filtering | |
CN103674030B (en) | The deviation of plumb line dynamic measurement device kept based on astronomical attitude reference and method | |
CN101545778B (en) | Initial alignment method of electric control compass based on nonlinear filtering | |
CN103344260B (en) | Based on the strapdown inertial navitation system (SINS) Initial Alignment of Large Azimuth Misalignment On method of RBCKF | |
CN102928858B (en) | GNSS (Global Navigation Satellite System) single-point dynamic positioning method based on improved expanded Kalman filtering | |
CN103278163A (en) | Nonlinear-model-based SINS/DVL (strapdown inertial navigation system/doppler velocity log) integrated navigation method | |
CN104344836B (en) | Posture observation-based redundant inertial navigation system fiber-optic gyroscope system level calibration method | |
CN103076025B (en) | A kind of optical fibre gyro constant error scaling method based on two solver | |
CN104019828A (en) | On-line calibration method for lever arm effect error of inertial navigation system in high dynamic environment | |
CN103245359A (en) | Method for calibrating fixed errors of inertial sensor in inertial navigation system in real time | |
CN103822633A (en) | Low-cost attitude estimation method based on second-order measurement update | |
CN103217699B (en) | Integrated navigation system recursion optimizing initial-alignment method based on polarization information | |
CN104374401A (en) | Compensating method of gravity disturbance in strapdown inertial navigation initial alignment | |
CN106643806B (en) | A kind of inertial navigation system alignment precision appraisal procedure |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20111116 Termination date: 20141217 |
|
EXPY | Termination of patent right or utility model |