CN102590837B - Method for calibrating phase errors of global positioning system (GPS) signal source - Google Patents
Method for calibrating phase errors of global positioning system (GPS) signal source Download PDFInfo
- Publication number
- CN102590837B CN102590837B CN201110394219.0A CN201110394219A CN102590837B CN 102590837 B CN102590837 B CN 102590837B CN 201110394219 A CN201110394219 A CN 201110394219A CN 102590837 B CN102590837 B CN 102590837B
- Authority
- CN
- China
- Prior art keywords
- phase
- frequency control
- control word
- correct
- point format
- 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
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
The invention relates to the technical field of simulative generation of GPS signal source, in particular to a method for calibrating phase errors of a GPS signal source. The technical scheme of the method includes calculating a frequency control word of a floating point format by utilizing an input Doppler frequency parameter; subjecting the frequency control word to a fix-point process and calculating an overall error produced within a calibrating time internal; controlling phase accumulation by using the frequency control word after the fix-point process; if the calibrating time is reached, adding the overall error into the accumulation result, storing and outputting the final result, and storing and outputting the accumulation result directly otherwise; and if a parameter updated time interval is reached, recalculating the frequency control word of the floating point format according to the Doppler frequency at the moment and repeating the above steps, and proceeding phase accumulation otherwise. According to the technical scheme of the method for calibrating the phase errors of the GPS signal source, the high accuracy of floating point accomplishment and simplicity and low cost of fixed point accomplishment are balanced and combined, and the finite word-length effect in the fixed point accomplishment process can be calibrated simultaneously.
Description
Technical field
The present invention relates to the technical field that simulation produces the gps signal source, be specifically related to a kind of method for the phase error of proofreading and correct the gps signal source.
Background technology
The gps signal source can provide repeatably simulating signal for every research, emulation and the test in satellite navigation field, contributes to the propelling of research work.Simultaneously, in the environment that is not easy to erect a television antenna, as the area of remote mountain areas and awful weather, the gps signal source also can guarantee the normal operation of navigation terminal, greatly reduces the construction cost and the maintenance cost that erect a television antenna.
In the simulation of gps signal source produces the process of gps signal, need to calculate the frequency control word of intermediate frequency carrier, the frequency control word of C/A code and the frequency control word of navigation message, and then with these three frequency control words, go control phase cumulative respectively, obtain respectively the phase place of intermediate frequency carrier, C/A code and navigation message.
Calculate the said frequencies control word and carry out phase-accumulated method and mainly comprise the steps:
Step 1: utilize the frequency control word of Doppler frequency calculation of parameter intermediate frequency carrier, the frequency control word of C/A code and the frequency control word of navigation message;
Step 2: according to sampling interval, the frequency control word that utilizes step 1 to obtain carries out phase-accumulated, obtains the phase place of intermediate frequency carrier, C/A code and navigation message.
Wherein, utilize the formula of the frequency control word FCW_Carrier_Float of Doppler frequency calculation of parameter intermediate frequency carrier to be:
F
imean IF-FRE, f
dmean Doppler frequency, f
smean the signal sampling frequency.
Utilize the formula of the frequency control word FCW_Code_Float of Doppler frequency calculation of parameter C/A code to be:
F
cfor the code check of C/A code, f
tfor the gps signal carrier frequency.
The formula that utilizes the Doppler frequency calculation of parameter to calculate the frequency control word FCW_Data_Float of navigation message is:
F
dmean the navigation message code check.
Step 2 is being carried out with frequency control word when phase-accumulated, if Doppler frequency used calculates for the first time for this satellite, should start from the initial phase of gps signal to add up, otherwise just continue cumulative on the cumulative phase-accumulated value basis obtained in last time.
As can be seen here, for the phase place of the intermediate frequency carrier, C/A code and the navigation message that calculate the gps signal source, the initial phase of Doppler frequency and this gps satellite signal is its basic input parameter, and these two parameters are generally that other modules by the gps signal source generate and export to the phase calculation module.
During the specific implementation above-mentioned steps, both can adopt the digital platform of fixed point on a digital platform, also can adopt the digital platform of floating-point.We the analysis found that, no matter adopt merely which kind of platform to realize, all there is following difficult point: if adopt the fixed point platform to realize, mean while with a time-limited binary string, meaning intermediate frequency carrier frequency control word, C/A code frequency control word and navigation message frequency control word, this just exists and means error, between the value that its fixed point mode means and ideal value, has deviation that is:.There is gap between the gps signal that this makes simulation produce and real signal.And, along with the increase of time, this gap may become greatly gradually, causes performance further to worsen.If mean and deal with data, can greatly reduce the error that finite word length effect causes, but this will increase the difficulty of realization, the cost of hardware platform and the power consumption of whole system on the floating-point platform.
From domestic and international disclosed document, we do not find the technical scheme that can effectively solve above-mentioned two aspect technical matterss simultaneously.
Summary of the invention
A kind of method for the phase error of proofreading and correct the gps signal source that the present invention proposes is characterized in that in turn including the following steps:
Step 1: the word length of setting digital platform is L
w, intermediate frequency carrier, C/A code and the navigation message signals phase error of largest tolerable separately are Max_Phase_Err, parameter interval T update time
u, the time interval T that finite word length effect is proofreaied and correct
e, the sampling time interval T of signal
s;
Step 2: utilize the Doppler frequency parameter of current time, calculate the frequency control word of floating-point format;
Step 3: it is L that the frequency control word of the floating-point format that step 2 is obtained converts word length to
wthe frequency control word of fixed point format, record the error that single conversion produces, and calculate the time interval T proofreaied and correct at finite word length effect
ethe summation of the error of interior generation, then preserve this sum of the deviations with fixed point format;
Step 4: the frequency control word of the fixed point format obtained by step 3, with the sampling time interval T of signal
sfor the time beat, the phase place of signal is added up, then judge whether the time of cumulative process reaches the time interval T that finite word length effect is proofreaied and correct
eif, reached this interval, first by this phase-accumulated value, add the sum of the deviations that step 3 obtains, obtain the phase error correction value, by this phase error correction value pair
preserve after delivery and export, then entering step 5, if do not reach, directly that phase-accumulated value is right
preserve and export this result after delivery, then enter step 5;
Step 5: judge whether to reach parameter interval T update time
uif, reach, return to step 2; If do not reach, return to step 4.
Parameter interval T update time
ube greater than or equal to the time interval T that finite word length effect is proofreaied and correct
e.
Parameter interval T update time
uthe time interval T that finite word length effect is proofreaied and correct
eintegral multiple.
Described frequency control word is in intermediate frequency carrier frequency control word, C/A code frequency control word and navigation message frequency control word.
The time interval T that finite word length effect is proofreaied and correct
eadopt following method to determine:
1) for each described frequency control word, calculate its needed minimum cumulative frequency Min_N=Max_Phase_Err/Max_ Δ FCW according to the described maximum phase error Max_Phase_Err that can tolerate respectively, wherein, the absolute value that Max_ Δ FCW is the maximum error that may cause while in described step 3, the frequency control word of floating-point format being converted to the frequency control word of fixed point format, for adopting the round up conversion method of mode of decimal place, Max_ Δ FCW=0.5, for adopting the conversion method of directly removing decimal place, Max_ Δ FCW=1;
2) for each described frequency control word, calculate respectively T
e=T
s* Min_N, wherein T
ssampling time interval for described signal;
3) get 2) three T obtaining of step
ein minimum value, be the time interval T that finite word length effect is proofreaied and correct
e.
It is L that described step 3 converts the frequency control word of floating-point format to word length
wthe step of fixed point format comprise:
5) by 1) value that obtains of step carry out round numbers computing that decimal place rounds up or by 1) value that obtains of step directly removes the rounding operation of decimal place.
The described time interval T proofreaied and correct at finite word length effect of step 3
eerror * T that the summation of the error of interior generation=described single conversion produces
e/ T
s.
It is that following formula is carried out in circulation: the frequency control word of phase-accumulated value=phase-accumulated value+fixed point format that described step 3 obtains that the phase place of described step 4 pair signal is carried out cumulative method, if be to calculate for the first time its Doppler frequency parameter for working as this gps satellite of front simulation, the prima facies place value that initial value of phase-accumulated value is signal, otherwise the initial value of phase-accumulated value is the phase-accumulated value preserved in described step 4.
Technical scheme of the present invention has been carried out good compromise on floating-point realization and fixed-point implementation, combine the accuracy of floating-point realization and the simplicity of fixed-point implementation, considered that floating-point format turns the error of introducing due to finite word length effect in the fixed point format process simultaneously, and it has been proofreaied and correct.Before the step 3 of the technical program, the calculated rate of parameter is lower (to be equaled parameter and upgrades interval T
u), use floating-point to realize, to calculate more accurately the value of flat rate control word according to the Doppler frequency of current time satellite-signal; After step 3, calculated rate is higher (is mainly by sampling interval T
scarry out phase-accumulated calculating), use fixed-point implementation, more the generation satellite-signal of Simple fast.Simultaneously, also calculate the accumulation of phase error that finite word length effect causes, and upgraded interval T in the parameter of fixed-point implementation part
eproofread and correct at place, has guaranteed the accuracy of fixed-point implementation.
The accompanying drawing explanation
The overview flow chart of a kind of method for the phase error of proofreading and correct the gps signal source that accompanying drawing 1 is that the present invention proposes.
Embodiment
Below illustrate the embodiment of technical solution of the present invention.
Step 1: initialization specifically comprises: setup parameter comprises gps signal carrier frequency f
t=1575.42MHz, C/A code check f
c=1.023MHz, navigation message code check f
d=50Hz, sample frequency f
s=16.384MHz, sampling interval T
s=1/16.384us, intermediate frequency f
i=5.12MHz, the word length L of digital platform
w=32bit,, the maximum phase error tolerance value that three frequency control word errors cause is equal 8192, adopts the method rounded up to realize fixed point, i.e. Max_ Δ FCW=0.5, the finite word length effect obtained is thus proofreaied and correct interval T
e=1ms, proofread and correct the number of data points N in interval
e=T
e/ T
s=16384, in addition, setting parameter upgrades interval T
u=2ms, N
u=T
u/ T
e=2.
Step 2: the frequency control word that calculates floating-point format according to the Doppler frequency of current time
(a) calculate the frequency control word FCW_Carrier_Float of intermediate frequency carrier,
(b) calculate the frequency control word FCW_Code_Float of C/A code,
(c) calculate the frequency control word FCW_Data_Float of navigation message,
F wherein
dfor Doppler frequency.
Step 3: frequency control word is carried out to the fixed point processing, and record the error that finite word length effect causes, specifically comprise the steps:
Step 3-1: the IF-FRE control word FCW_Carrier_Float to floating-point format is multiplied by
after carry out, to the rounding up of decimal place, obtaining the IF-FRE control word FCW_Carrier of fixed point format,
FCW_Carrier=round{FCW_Carrier_Float×2
32}
Calculating is proofreaied and correct interval T at finite word length effect
ein, the intermediate frequency carrier accumulation of phase error delta Phase_Carrier caused by finite word length effect, still mean by fixed point format,
ΔPhase_Carrier=round{(FCW_Carrier_Float×2
32-FCW_Carrier)×16384}
Step 3-2: the C/A code frequency control word FCW_Code_Float to floating-point format is multiplied by
after carry out, to the rounding up of decimal place, obtaining the C/A code frequency control word FCW_Code of fixed point format,
FCW_Code=round{FCW_Code_Float×2
32}
Calculate finite word length effect and proofread and correct interval T
ein, the C/A code phase cumulative errors Δ Phase_Code caused by finite word length effect, still mean by fixed point format,
ΔPhase_Code=round{(FCW_Code_Float×2
32-FCW_Code)×16384}
Step 3-3: the navigation message frequency control word FCW_Data_Float to floating-point format is multiplied by
after carry out, to the rounding up of decimal place, obtaining the navigation message frequency control word FCW_Data of fixed point format,
FCW_Data=round{FCW_Data_Float×2
32}
Calculate finite word length effect and proofread and correct interval T
ein, the navigation message accumulation of phase error delta Phase_Data caused by finite word length effect, still mean by fixed point format,
ΔPhase_Data=round{(FCW_Data_Float×2
32-FCW_Data)×16384}
Round{} computing in above-mentioned formula, the round up operation of expression to decimal place, except adopting the mode that decimal place is rounded up to complete floating-point format to the conversion of fixed point format, can also directly decimal place be given up, the way that only retains integer-bit converts the frequency control word of fixed point format to, if adopt this method, and Max_ Δ FCW=1, other parameters relevant to Max_ Δ FCW, also should make corresponding change.
Step 4: utilize the frequency control word of fixed point format to carry out phase-accumulated, specifically comprise the steps:
Step 4-1: utilize the frequency control word of intermediate frequency carrier to carry out phase-accumulated,
Phase_Carrier=Phase_Carrier+FCW_Carrier
Step 4-2: utilize the frequency control word of C/A code to carry out phase-accumulated,
Phase_Code=Phase_Code+FCW_Code
Step 4-3: utilize the frequency control word of navigation message to carry out phase-accumulated,
Phase_Data=Phase_Data+FCW_Data
Before the beginning that phase place added up with frequency control word, whether the Doppler frequency that needs judgement to adopt calculates for the first time for this satellite, if so, and the first phase that initial value of the phase place in top three cumulative formula is this satellite-signal; If not, the initial value of the phase place in top three cumulative formula is the last cumulative phase-accumulated value preserved later.
Step 5: if arrive finite word length effect, proofread and correct interval T
e, utilize the cumulative errors amount to be proofreaied and correct finite word length effect, and will proofread and correct later accumulation result and preserve, specifically comprise the steps:
Step 5-1: proofread and correct the intermediate frequency carrier phase place,
Phase_Carrier=Phase_Carrier+ΔPhase_Carrier
Step 5-2: proofread and correct the C/A code phase,
Phase_Code=Phase_Code+ΔPhase_Code
Step 5-3: proofread and correct the navigation message phase place,
Phase_Data=Phase_Data+ΔPhase_Data
Do not proofread and correct interval T if reach finite word length effect
e, directly accumulation result is preserved.
Step 6: three phase-accumulated values that step 5 is preserved are by 2
32carry out delivery output, specifically comprise the steps:
Step 6-1: the intermediate frequency carrier phase is carried out to delivery,
Phase_Carrier=mod{Phase_Carrier,2
32}
Step 6-2: the C/A code phase is carried out to delivery,
Phase_Code=mod{Phase_Code,2
32}
Step 6-3: the navigation message phase place is carried out to delivery,
Phase_Data=mod{Phase_Data,2
32}
In the fixed-point implementation process, delivery is exactly simply to get low level, to 2
32delivery is exactly to take out simply low 32 of phase-accumulated value to get final product.
Step 7: judge whether to have reached parameter and upgrade interval T
uif, reach, return to step 2, otherwise return to step 4.
The step 2 of this specific embodiment adopts floating system to realize, step 3,4,5,6 adopts the fixed point mode to realize, from top step, describes and can find out, has only reached parameter and has upgraded interval T
uthe time, just can require the floating-point platform to recalculate each frequency control word, and parameter is upgraded interval T
uvalue generally all than finite word length effect, proofread and correct interval T
egrow, this has just well taken into account precision and system cost, makes the floating-point operation amount of expensive and high power consumption fewer, and a large amount of computings is placed on the fixed point platform and realizes.Therefore, the technical program is when specific implementation, and suggestion realizes that by floating-point part and fixed-point implementation partly separate, be placed on different platforms, as adopted the pattern of " DSP+FPGA ", DSP is responsible for the calculating that floating-point is relevant, FPGA is responsible for phase-accumulated and tables look-up, and jointly completes the task of simulating GPS signal.
Claims (6)
1. the method for the phase error of proofreading and correct the gps signal source is characterized in that in turn including the following steps:
Step 1: the word length of setting digital platform is L
w, intermediate frequency carrier, C/A code and the navigation message signals phase error of largest tolerable separately are Max_Phase_Err, parameter interval T update time
u, the time interval T that finite word length effect is proofreaied and correct
e, the sampling time interval T of signal
s;
Step 2: utilize the Doppler frequency parameter of current time, calculate the frequency control word of floating-point format;
Step 3: it is L that the frequency control word of the floating-point format that step 2 is obtained converts word length to
wthe frequency control word of fixed point format, record the error that single conversion produces, and calculate the time interval T proofreaied and correct at finite word length effect
ethe summation of the error of interior generation, then preserve this sum of the deviations with fixed point format;
Step 4: the frequency control word of the fixed point format obtained by step 3, with the sampling time interval T of signal
sfor the time beat, the phase place of signal is added up, then judge whether the time of cumulative process reaches the time interval T that finite word length effect is proofreaied and correct
eif, reached this interval, first by this phase-accumulated value, add the sum of the deviations that step 3 obtains, obtain the phase error correction value, by this phase error correction value pair
preserve after delivery and export, then entering step 5, if do not reach, directly that phase-accumulated value is right
preserve and export this result after delivery, then enter step 5;
Step 5: judge whether to reach parameter interval T update time
uif, reach, return to step 2; If do not reach, return to step 4;
Wherein, described frequency control word is in intermediate frequency carrier frequency control word, C/A code frequency control word and navigation message frequency control word; The time interval Te that finite word length effect is proofreaied and correct adopts following method to determine:
1) for each described frequency control word, calculate its needed minimum cumulative frequency Min_N=Max_Phase_Err/Max_ Δ FCW according to the described maximum phase error Max_Phase_Err that can tolerate respectively, wherein, the absolute value that Max_ Δ FCW is the maximum error that may cause while in described step 3, the frequency control word of floating-point format being converted to the frequency control word of fixed point format, for adopting the round up conversion method of mode of decimal place, Max_ Δ FCW=0.5, for adopting the conversion method of directly removing decimal place, Max_ Δ FCW=1;
2) for each described frequency control word, calculate respectively Te=Ts * Min_N, the sampling time interval that wherein Ts is described signal;
3) get 2) minimum value in three Te obtaining of step, be the time interval Te that finite word length effect is proofreaied and correct.
2. a kind of method for the phase error of proofreading and correct the gps signal source according to claim 1, is characterized in that: parameter interval T update time
ube greater than or equal to the time interval T that finite word length effect is proofreaied and correct
e.
3. a kind of method for the phase error of proofreading and correct the gps signal source according to claim 1, is characterized in that: parameter interval T update time
uthe time interval T that finite word length effect is proofreaied and correct
eintegral multiple.
4. a kind of method for the phase error of proofreading and correct the gps signal source according to claim 1, it is characterized in that: it is L that described step 3 converts the frequency control word of floating-point format to word length
wthe step of fixed point format comprise:
2) by 1) value that obtains of step carry out round numbers computing that decimal place rounds up or by 1) value that obtains of step directly removes the rounding operation of decimal place.
5. a kind of method for the phase error of proofreading and correct the gps signal source according to claim 1, is characterized in that: relate to the rapid 3 described time interval T that proofread and correct at finite word length effect
eerror * T that the summation of the error of interior generation=described single conversion produces
e/ T
s.
6. a kind of method for the phase error of proofreading and correct the gps signal source according to claim 1, it is characterized in that: it is that following formula is carried out in circulation: the frequency control word of phase-accumulated value=phase-accumulated value+fixed point format that described step 3 obtains that the phase place of described step 4 pair signal is carried out cumulative method, if be to calculate for the first time its Doppler frequency parameter for working as this gps satellite of front simulation, the prima facies place value that initial value of phase-accumulated value is signal, otherwise the initial value of phase-accumulated value is the phase-accumulated value preserved in described step 4.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110394219.0A CN102590837B (en) | 2011-11-22 | 2011-11-22 | Method for calibrating phase errors of global positioning system (GPS) signal source |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110394219.0A CN102590837B (en) | 2011-11-22 | 2011-11-22 | Method for calibrating phase errors of global positioning system (GPS) signal source |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102590837A CN102590837A (en) | 2012-07-18 |
CN102590837B true CN102590837B (en) | 2014-01-08 |
Family
ID=46479747
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110394219.0A Expired - Fee Related CN102590837B (en) | 2011-11-22 | 2011-11-22 | Method for calibrating phase errors of global positioning system (GPS) signal source |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102590837B (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113132027B (en) * | 2019-12-30 | 2023-02-10 | 江西联智集成电路有限公司 | Method and apparatus for correcting operating frequency of radio transmitter |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4445118A (en) * | 1981-05-22 | 1984-04-24 | The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration | Navigation system and method |
US5194871A (en) * | 1982-03-01 | 1993-03-16 | Western Atlas International, Inc. | System for simultaneously deriving position information from a plurality of satellite transmissions |
GB2245445A (en) * | 1990-06-18 | 1992-01-02 | Philips Electronic Associated | Method of and apparatus for obtaining vehicle heading information |
US5418538A (en) * | 1993-05-21 | 1995-05-23 | Trimble Navigation Limited | Rapid satellite signal acquisition in a satellite positioning system |
-
2011
- 2011-11-22 CN CN201110394219.0A patent/CN102590837B/en not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN102590837A (en) | 2012-07-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9261597B2 (en) | Real-time and post-processed orbit determination and positioning | |
CN108732597B (en) | Method and system for establishing time reference of multi-satellite navigation system | |
CN103323888B (en) | Method for eliminating delay errors of troposphere of GNSS atmospheric probing data | |
CN102197317B (en) | Satellite navigation device | |
CN101920762B (en) | Ship dynamic positioning method based on real-time correction of noise matrix | |
CN101261317A (en) | High dynamic state multi- mode satellite navigation signal source analogue method and its device | |
CN103163533B (en) | Seamless fusion expression and correction method of global navigation satellite system (GNSS) global and regional ionospheric delay | |
RU2010146231A (en) | METHOD AND DEVICE FOR DETERMINING POSITION USING SPS ORBIT HYBRID DATA | |
WO2010049564A3 (en) | Optimized-cost method for computer-assisted calculation of the aerodynamic forces in an aircraft | |
CN107367739A (en) | A kind of GLONASS civil navigation signals simulation production method and device | |
CN106093967A (en) | The ionosphere delay method for solving that a kind of pseudorange phase place is comprehensive | |
CN107544082A (en) | The step modeling of Big Dipper IGSO/MEO satellite pseudorange codes deviation one | |
CN103777518A (en) | Satellite navigation timing system and method based on improved least square method | |
CN102590837B (en) | Method for calibrating phase errors of global positioning system (GPS) signal source | |
CN103278826A (en) | Beidou B1 frequency point intermediate frequency signal simulation method | |
CN117388881A (en) | Method and system for tracing satellite-borne atomic clock of low-orbit satellite to UTC (k) | |
CN110398764A (en) | A kind of Beidou Navigation System two-band fusion and positioning method and device | |
CN114286286A (en) | Time synchronization method, apparatus, medium, and program product | |
CN111398994B (en) | Method and device for positioning and time service of medium-orbit communication satellite | |
CN103869345B (en) | For associating time synchronized and the localization method of Beidou satellite navigation and positioning system | |
CN107505637B (en) | Hot start positioning method based on relative pseudo-range | |
CN109738912B (en) | Method for realizing fixed point time service based on GNSS satellite signals | |
CN111399022A (en) | Satellite navigation and inertial navigation ultra-tight combined simulation positioning method and terminal | |
CN115603847A (en) | Subnanosecond time synchronization method and device based on Beidou short message and PPP-B2B service | |
CN113341686B (en) | Single-satellite multi-calendar timing method and device, computer equipment and storage medium |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20140108 Termination date: 20161122 |