CN102710349A - Data selection-based frequency spectrum sensing method used under pulse interference environment - Google Patents
Data selection-based frequency spectrum sensing method used under pulse interference environment Download PDFInfo
- Publication number
- CN102710349A CN102710349A CN2012101819317A CN201210181931A CN102710349A CN 102710349 A CN102710349 A CN 102710349A CN 2012101819317 A CN2012101819317 A CN 2012101819317A CN 201210181931 A CN201210181931 A CN 201210181931A CN 102710349 A CN102710349 A CN 102710349A
- Authority
- CN
- China
- Prior art keywords
- absolute value
- centrifugal pump
- signal
- time
- domain sampling
- 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
Abstract
The invention discloses a data selection-based frequency spectrum sensing method used under a pulse interference environment, which comprises the steps as follows: firstly, an antenna in a cognitive radio system repeatedly executes time-domain sampling aiming at a signal in a monitoring channel; secondly, an absolute value discrete value of the signal of each time domain sampling point is computed, and each discrete value in a first absolute value discrete value sequence formed by all discrete values is compared with the average value of all the discrete values, when two or more than two continuous discrete values in the first absolute value discrete value sequence are large than the average value of the discrete values by three times, and the discrete values with interference pulse are removed; thirdly, the inspection statistical quantity of a second absolute value discrete value sequence formed by all discrete values without pulse interference is computed; and lastly, the sizes of the inspection statistical quantity and a decision threshold are compared, and whether the monitoring channel is idle or not is judged according to a comparison result. The method disclosed by the invention effectively solves a problem that the false alarm probability is increased owing to existing pulse interference, and improves the utilization ratio of frequency spectrum resources.
Description
Technical field
The present invention relates to the frequency spectrum perception technology in a kind of cognitive radio system, especially relate to the frequency spectrum sensing method of selecting based on data under a kind of impulse disturbances environment.
Background technology
Along with the appearance of increasing radio communication service, situation more and more in short supply has appearred in limited frequency spectrum resources.And in reality; The physics frequency spectrum resource that is assigned to specific communication service in a large number is in idle condition on many times, space; This shows that it is not to be lacked by the physics frequency spectrum resource to cause that present frequency spectrum resource lacks phenomenon, but since fixing frequency spectrum resource distribution policy cause.The proposition of cognitive radio provides a feasible thinking for improving the availability of frequency spectrum; Promptly effectively utilize the mandate frequency range that is in idle condition; But before utilizing these mandate frequency ranges that are in idle condition; Which must effectively detect earlier authorize frequency range idle condition to occur, so frequency spectrum perception is one of key technology in the cognitive radio.
Existing frequency spectrum sensing method has energy measuring method, covariance matrix detection method, feature detection method etc., these frequency spectrum sensing methods all are to consider under the stationary noise environment, directly utilizes all sampled datas to carry out frequency spectrum detection.But under the situation of the impulse disturbances that exists thunder and lightning, automobile sparking etc. to cause, the false alarm probability of existing frequency spectrum sensing method can sharply increase, and then makes frequency spectrum resource utilization rate seriously low.
Summary of the invention
Technical problem to be solved by this invention provides the frequency spectrum sensing method of selecting based on data under a kind of impulse disturbances environment, and it can effectively be reduced in the false alarm probability under the impulse disturbances environment, improves frequency spectrum resource utilization rate.
The present invention solves the problems of the technologies described above the technical scheme that is adopted: the frequency spectrum sensing method of selecting based on data under a kind of impulse disturbances environment is characterized in that may further comprise the steps: at first the antenna in the cognitive radio system carries out repeatedly time-domain sampling to the signal in the supervisory channel and obtains the time-domain sampling signal that the signal by a plurality of time-domain sampling points constitutes; Calculate the absolute value centrifugal pump of the signal of each time-domain sampling point then; And the size of the mean value of the absolute value centrifugal pump of the signal of each the absolute value centrifugal pump in the first absolute value centrifugal pump sequence that relatively constitutes by the absolute value centrifugal pump of the signal of all time-domain sampling points and all time-domain sampling points, when continuous two and above absolute value centrifugal pump in the first absolute value centrifugal pump sequence are rejected the absolute value centrifugal pump that there is impulse disturbances in these during greater than 3 times of the mean value of the absolute value centrifugal pump of the signal of all time-domain sampling points; Then calculate the test statistics that does not have the second absolute value centrifugal pump sequence that the absolute value centrifugal pump of impulse disturbances constitutes by all; The size of last comparing check statistic and decision threshold judges according to comparative result whether this supervisory channel is in idle condition.
It comprises following concrete steps:
1. the antenna in the cognitive radio system carries out time-domain sampling K time to the signal in the supervisory channel, obtains the time-domain sampling signal that the signal by K time-domain sampling point constitutes, and the signal of k time-domain sampling point is designated as x (k), wherein k ∈ [1, K];
2. the signal of each time-domain sampling point in the time-domain sampling signal is asked absolute value operation respectively; Obtain the absolute value centrifugal pump of the signal of each the time-domain sampling point in the time-domain sampling signal; The absolute value centrifugal pump of the signal x (k) of k time-domain sampling point is designated as y (k); Y (k)=| x (k) |, wherein, symbol " || " is the modulo operation symbol;
3. calculate the mean value of absolute value centrifugal pump of the signal of all the time-domain sampling points in the time-domain sampling signal, be designated as m
y,
4. the absolute value centrifugal pump with the signal of all time-domain sampling points constitutes one first absolute value centrifugal pump sequence, then each absolute value centrifugal pump and m in the first absolute value centrifugal pump sequence relatively
ySize, if continuous two and above absolute value centrifugal pump are all greater than m in the first absolute value centrifugal pump sequence
y3 times, then judge in these absolute value centrifugal pumps to have impulse disturbances, and reject the absolute value centrifugal pump that there is impulse disturbances in these, if do not exist continuous two and above absolute value centrifugal pump in the first absolute value centrifugal pump sequence all greater than m
y3 times, judge that then all there is not impulse disturbances in all the absolute value centrifugal pumps in the first absolute value centrifugal pump sequence;
5. calculate the test statistics that does not have the second absolute value centrifugal pump sequence that the absolute value centrifugal pump of impulse disturbances constitutes by all; Be designated as T;
wherein; N ∈ [1, N], N represent the number of the absolute value centrifugal pump in the second absolute value centrifugal pump sequence; N absolute value centrifugal pump in z (n) the expression second absolute value centrifugal pump sequence, symbol " || " is the modulo operation symbol;
6. according to false alarm probability P
FCalculate decision threshold, be designated as λ
T,
Wherein, σ
2The power of expression noise, Q
-1(P
F) expression Q (P
F) counter-transformation function, Q (P
F) be the tail area function of standardized normal distribution probability density,
T is the calculus variable;
7. comparing check statistic T and decision threshold λ
TSize, if T>=λ
T, judge that then this supervisory channel is in busy condition, if T<λ
T, judge that then this supervisory channel is in idle condition.
Compared with prior art; The invention has the advantages that through utilizing the comparison of time-domain sampling signal amplitude, reject the sampled point that has impulse disturbances, utilize remaining sampled point then; Whether the supervisory channel of judging the current perception cycle through the size of judging inspection statistics value and decision threshold is in idle condition; Realize frequency spectrum perception, the inventive method has effectively overcome the existing frequency spectrum sensing method problem that false alarm probability increases when having impulse disturbances, thereby has improved frequency spectrum resource utilization rate.
Description of drawings
Fig. 1 is the FB(flow block) of frequency spectrum sensing method of the present invention;
Fig. 2 is under different impulse disturbances and the noise power ratio, the comparison sketch map of the false alarm probability of existing energy measuring method and the inventive method.
Embodiment
Embodiment describes in further detail the present invention below in conjunction with accompanying drawing.
The frequency spectrum sensing method of selecting based on data under a kind of impulse disturbances environment that the present invention proposes, its FB(flow block) is as shown in Figure 1.It may further comprise the steps: at first the antenna in the cognitive radio system carries out repeatedly time-domain sampling to the signal in the supervisory channel and obtains the time-domain sampling signal that the signal by a plurality of time-domain sampling points constitutes; Calculate the absolute value centrifugal pump of the signal of each time-domain sampling point then; And the size of the mean value of the absolute value centrifugal pump of the signal of each the absolute value centrifugal pump in the first absolute value centrifugal pump sequence that relatively constitutes by the absolute value centrifugal pump of the signal of all time-domain sampling points and all time-domain sampling points, when continuous two and above absolute value centrifugal pump in the first absolute value centrifugal pump sequence are rejected the absolute value centrifugal pump that there is impulse disturbances in these during greater than 3 times of the mean value of the absolute value centrifugal pump of the signal of all time-domain sampling points; Then calculate the test statistics that does not have the second absolute value centrifugal pump sequence that the absolute value centrifugal pump of impulse disturbances constitutes by all; The size of last comparing check statistic and decision threshold judges according to comparative result whether this supervisory channel is in idle condition.Specifically comprise following concrete steps:
1. the antenna in the cognitive radio system carries out time-domain sampling K time to the signal in the supervisory channel, obtains the time-domain sampling signal that the signal by K time-domain sampling point constitutes, and the signal of k time-domain sampling point is designated as x (k), wherein k ∈ [1, K].In actual mechanical process, generally the signal in the supervisory channel is carried out minimum hundreds of time time-domain sampling.
2. the signal of each time-domain sampling point in the time-domain sampling signal is asked absolute value operation respectively; Obtain the absolute value centrifugal pump of the signal of each the time-domain sampling point in the time-domain sampling signal; The absolute value centrifugal pump of the signal x (k) of k time-domain sampling point is designated as y (k); Y (k)=| x (k) |, wherein, symbol " || " is the modulo operation symbol.
3. calculate the mean value of absolute value centrifugal pump of the signal of all the time-domain sampling points in the time-domain sampling signal, be designated as m
y,
4. the absolute value centrifugal pump with the signal of all time-domain sampling points constitutes one first absolute value centrifugal pump sequence, then each absolute value centrifugal pump and m in the first absolute value centrifugal pump sequence relatively
ySize, if continuous two and above absolute value centrifugal pump are all greater than m in the first absolute value centrifugal pump sequence
y3 times, then judge in these absolute value centrifugal pumps to have impulse disturbances, and reject the absolute value centrifugal pump that there is impulse disturbances in these, if do not exist continuous two and above absolute value centrifugal pump in the first absolute value centrifugal pump sequence all greater than m
y3 times, judge that then all there is not impulse disturbances in all the absolute value centrifugal pumps in the first absolute value centrifugal pump sequence.
5. calculate the test statistics that does not have the second absolute value centrifugal pump sequence that the absolute value centrifugal pump of impulse disturbances constitutes by all; Be designated as T;
wherein; N ∈ [1, N], N represent the number of the absolute value centrifugal pump in the second absolute value centrifugal pump sequence; N absolute value centrifugal pump in z (n) the expression second absolute value centrifugal pump sequence, symbol " || " is the modulo operation symbol.At this; The second absolute value centrifugal pump sequence that constitutes comprises two kinds of situation; First kind of situation is to reject in the first absolute value centrifugal pump sequence to have the new absolute value centrifugal pump sequence that forms after the absolute value centrifugal pump of impulse disturbances, and the number of the absolute value centrifugal pump that comprises in this new absolute value centrifugal pump sequence is less than the number of the absolute value centrifugal pump that comprises in the first absolute value centrifugal pump sequence; Second kind of situation is to need not to reject any the absolute value centrifugal pump in the first absolute value centrifugal pump sequence, and the absolute value centrifugal pump that promptly comprises in the second absolute value centrifugal pump sequence is original absolute value centrifugal pump in the first absolute value centrifugal pump sequence.
6. according to false alarm probability P
FCalculate decision threshold, be designated as λ
T,
Wherein, σ
2The power of expression noise, Q
-1(P
F) expression Q (P
F) counter-transformation function, Q (P
F) be the tail area function of standardized normal distribution probability density,
T is the calculus variable.
7. comparing check statistic T and decision threshold λ
TSize, if T>=λ
T, judge that then this supervisory channel is in busy condition, if T<λ
T, judge that then this supervisory channel is in idle condition.
Through feasibility and the validity of following emulation to further specify frequency spectrum sensing method of the present invention.
Suppose the signal time-domain sampling number of times K=1000 in each perception cycle of cognitive radio system, wherein have impulsive noise in the signal of 30 continuous sampling points, preset false alarm probability is P
f=0.1.Fig. 2 has provided under different impulse disturbances and the noise power ratio, and existing energy measuring method and the inventive method are through the comparison of the false alarm probability of 10000 Monte Carlo simulation acquisitions.As can be seen from Figure 2; When impulse disturbances and noise power are relatively lower; Existing energy measuring method can realize the false alarm probability preset, and along with the increase of impulse disturbances power, it departs from preset false alarm probability can become big gradually; Even reach 1, and the inventive method can make the false alarm probability of emulation maintain all the time near the preset false alarm probability.
Claims (2)
1. the frequency spectrum sensing method of selecting based on data under the impulse disturbances environment is characterized in that may further comprise the steps: at first the antenna in the cognitive radio system carries out repeatedly time-domain sampling to the signal in the supervisory channel and obtains the time-domain sampling signal that the signal by a plurality of time-domain sampling points constitutes; Calculate the absolute value centrifugal pump of the signal of each time-domain sampling point then; And the size of the mean value of the absolute value centrifugal pump of the signal of each the absolute value centrifugal pump in the first absolute value centrifugal pump sequence that relatively constitutes by the absolute value centrifugal pump of the signal of all time-domain sampling points and all time-domain sampling points, when continuous two and above absolute value centrifugal pump in the first absolute value centrifugal pump sequence are rejected the absolute value centrifugal pump that there is impulse disturbances in these during greater than 3 times of the mean value of the absolute value centrifugal pump of the signal of all time-domain sampling points; Then calculate the test statistics that does not have the second absolute value centrifugal pump sequence that the absolute value centrifugal pump of impulse disturbances constitutes by all; The size of last comparing check statistic and decision threshold judges according to comparative result whether this supervisory channel is in idle condition.
2. the frequency spectrum sensing method of selecting based on data under a kind of impulse disturbances environment according to claim 1 is characterized in that comprising following concrete steps:
1. the antenna in the cognitive radio system carries out time-domain sampling K time to the signal in the supervisory channel, obtains the time-domain sampling signal that the signal by K time-domain sampling point constitutes, and the signal of k time-domain sampling point is designated as x (k), wherein k ∈ [1, K];
2. the signal of each time-domain sampling point in the time-domain sampling signal is asked absolute value operation respectively; Obtain the absolute value centrifugal pump of the signal of each the time-domain sampling point in the time-domain sampling signal; The absolute value centrifugal pump of the signal x (k) of k time-domain sampling point is designated as y (k); Y (k)=| x (k) |, wherein, symbol " || " is the modulo operation symbol;
3. calculate the mean value of absolute value centrifugal pump of the signal of all the time-domain sampling points in the time-domain sampling signal, be designated as m
y,
4. the absolute value centrifugal pump with the signal of all time-domain sampling points constitutes one first absolute value centrifugal pump sequence, then each absolute value centrifugal pump and m in the first absolute value centrifugal pump sequence relatively
ySize, if continuous two and above absolute value centrifugal pump are all greater than m in the first absolute value centrifugal pump sequence
y3 times, then judge in these absolute value centrifugal pumps to have impulse disturbances, and reject the absolute value centrifugal pump that there is impulse disturbances in these, if do not exist continuous two and above absolute value centrifugal pump in the first absolute value centrifugal pump sequence all greater than m
y3 times, judge that then all there is not impulse disturbances in all the absolute value centrifugal pumps in the first absolute value centrifugal pump sequence;
5. calculate the test statistics that does not have the second absolute value centrifugal pump sequence that the absolute value centrifugal pump of impulse disturbances constitutes by all; Be designated as T;
wherein; N ∈ [1, N], N represent the number of the absolute value centrifugal pump in the second absolute value centrifugal pump sequence; N absolute value centrifugal pump in z (n) the expression second absolute value centrifugal pump sequence, symbol " || " is the modulo operation symbol;
6. according to false alarm probability P
FCalculate decision threshold, be designated as λ
T,
Wherein, σ
2The power of expression noise, Q
-1(P
F) expression Q (P
F) counter-transformation function, Q (P
F) be the tail area function of standardized normal distribution probability density,
T is the calculus variable;
7. comparing check statistic T and decision threshold λ
TSize, if T>=λ
T, judge that then this supervisory channel is in busy condition, if T<λ
T, judge that then this supervisory channel is in idle condition.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210181931.7A CN102710349B (en) | 2012-05-31 | 2012-05-31 | Data selection-based frequency spectrum sensing method used under pulse interference environment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210181931.7A CN102710349B (en) | 2012-05-31 | 2012-05-31 | Data selection-based frequency spectrum sensing method used under pulse interference environment |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN102710349A true CN102710349A (en) | 2012-10-03 |
| CN102710349B CN102710349B (en) | 2014-09-10 |
Family
ID=46902931
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201210181931.7A Expired - Fee Related CN102710349B (en) | 2012-05-31 | 2012-05-31 | Data selection-based frequency spectrum sensing method used under pulse interference environment |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN102710349B (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103346847A (en) * | 2013-06-21 | 2013-10-09 | 哈尔滨工业大学 | United compressing spectrum sensing method based on iteration attack detection |
| CN104994046A (en) * | 2015-07-14 | 2015-10-21 | 宁波大学 | Interframe spectrum sensing method in cognitive radio system |
| CN105375996A (en) * | 2015-10-12 | 2016-03-02 | 桂林电子科技大学 | Frequency spectrum sensing method based on sequence statistics in impulsive noise environment |
| CN117496582A (en) * | 2023-12-28 | 2024-02-02 | 苏州元脑智能科技有限公司 | Face recognition model training method and device, electronic equipment and storage medium |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104703190A (en) * | 2015-03-16 | 2015-06-10 | 哈尔滨工业大学深圳研究生院 | Combined space-time frequency spectrum sharing method in cognitive radio system |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20090323766A1 (en) * | 2006-03-16 | 2009-12-31 | The Boeing Company | Method and device of peak detection in preamble synchronization for direct sequence spread spectrum communication |
| CN101630983A (en) * | 2009-07-30 | 2010-01-20 | 哈尔滨工业大学 | Idle frequency spectrum detecting method by using cyclic spectrum statistic value in cognitive radio |
| CN101848044A (en) * | 2010-05-20 | 2010-09-29 | 北京邮电大学 | Low power consumption time domain and frequency domain double threshold combined energy detection algorithm |
| US7860197B2 (en) * | 2006-09-29 | 2010-12-28 | Samsung Electro-Mechanics | Spectrum-sensing algorithms and methods |
| CN101944961A (en) * | 2010-09-03 | 2011-01-12 | 电子科技大学 | Double threshold cooperative sensing method in cognitive wireless network |
-
2012
- 2012-05-31 CN CN201210181931.7A patent/CN102710349B/en not_active Expired - Fee Related
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20090323766A1 (en) * | 2006-03-16 | 2009-12-31 | The Boeing Company | Method and device of peak detection in preamble synchronization for direct sequence spread spectrum communication |
| US7860197B2 (en) * | 2006-09-29 | 2010-12-28 | Samsung Electro-Mechanics | Spectrum-sensing algorithms and methods |
| CN101630983A (en) * | 2009-07-30 | 2010-01-20 | 哈尔滨工业大学 | Idle frequency spectrum detecting method by using cyclic spectrum statistic value in cognitive radio |
| CN101848044A (en) * | 2010-05-20 | 2010-09-29 | 北京邮电大学 | Low power consumption time domain and frequency domain double threshold combined energy detection algorithm |
| CN101944961A (en) * | 2010-09-03 | 2011-01-12 | 电子科技大学 | Double threshold cooperative sensing method in cognitive wireless network |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103346847A (en) * | 2013-06-21 | 2013-10-09 | 哈尔滨工业大学 | United compressing spectrum sensing method based on iteration attack detection |
| CN103346847B (en) * | 2013-06-21 | 2015-02-04 | 哈尔滨工业大学 | United compressing spectrum sensing method based on iteration attack detection |
| CN104994046A (en) * | 2015-07-14 | 2015-10-21 | 宁波大学 | Interframe spectrum sensing method in cognitive radio system |
| CN104994046B (en) * | 2015-07-14 | 2018-01-05 | 宁波大学 | A kind of interframe frequency spectrum sensing method in cognitive radio system |
| CN105375996A (en) * | 2015-10-12 | 2016-03-02 | 桂林电子科技大学 | Frequency spectrum sensing method based on sequence statistics in impulsive noise environment |
| CN105375996B (en) * | 2015-10-12 | 2017-10-03 | 桂林电子科技大学 | The frequency spectrum sensing method based on order statistics under impulse noise environment |
| CN117496582A (en) * | 2023-12-28 | 2024-02-02 | 苏州元脑智能科技有限公司 | Face recognition model training method and device, electronic equipment and storage medium |
Also Published As
| Publication number | Publication date |
|---|---|
| CN102710349B (en) | 2014-09-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102710349B (en) | Data selection-based frequency spectrum sensing method used under pulse interference environment | |
| KR101689795B1 (en) | Method and system for detecting vulnerabilities of communication protocol software | |
| CN102324959B (en) | Frequency spectrum sensing method based on multi-aerial system covariance matrix | |
| CN102271022B (en) | Spectrum sensing method based on maximum generalized characteristic value | |
| CN104483602A (en) | Local discharge signal identification method and device | |
| CN104535905A (en) | Partial discharge diagnosis method based on naive bayesian classification | |
| CN102710345B (en) | Cognition radio frequency spectrum sensing method based on multi-antenna Friedman inspection | |
| Honghao et al. | Spectrum anomalies autonomous detection in cognitive radio using hidden markov models | |
| CN107071788B (en) | Spectrum sensing method and device in cognitive wireless network | |
| CN106170139B (en) | A kind of frequency spectrum detecting method and system | |
| CN107276691B (en) | Method and device for determining direction angle of antenna | |
| CN102590647A (en) | Audible noise measuring method for high voltage DC transmission line in corona cage | |
| CN102457853A (en) | Method and device capable of dividing cell clusters | |
| CN108680798B (en) | Lightning monitoring and early warning method and system | |
| CN104361529A (en) | Reliability detecting and evaluating method of power distribution system on basis of cloud model | |
| CN102075263B (en) | Segmented energy processing-based truncated sequential test method | |
| CN104821852A (en) | Frequency spectrum sensing method based on multi-antenna instantaneous power | |
| CN103346847B (en) | United compressing spectrum sensing method based on iteration attack detection | |
| CN104954089B (en) | A kind of frequency spectrum sensing method being compared based on multiple antennas instantaneous power | |
| CN104363064A (en) | Cooperative spectrum sensing method based on preference users | |
| CN106452627B (en) | A kind of noise power estimation method and device for broader frequency spectrum perception | |
| CN106936527B (en) | Multi-user cooperative spectrum sensing method based on double thresholds | |
| CN104320209B (en) | A kind of frequency spectrum sensing method based on the test of fitness of fot | |
| KR101704812B1 (en) | Apparatus and method to estimate noise and interference power in a communication system | |
| KR101542659B1 (en) | Determining method for presence of primary user's signal using spectrum sensing and spectrum sensing apparatus |
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: 20140910 Termination date: 20170531 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |