CN100472223C - Anti-RF interference method for high frequency radar - Google Patents

Anti-RF interference method for high frequency radar Download PDF

Info

Publication number
CN100472223C
CN100472223C CNB200610018265XA CN200610018265A CN100472223C CN 100472223 C CN100472223 C CN 100472223C CN B200610018265X A CNB200610018265X A CN B200610018265XA CN 200610018265 A CN200610018265 A CN 200610018265A CN 100472223 C CN100472223 C CN 100472223C
Authority
CN
China
Prior art keywords
signal
interference
radar
frequency
spectrum
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
Application number
CNB200610018265XA
Other languages
Chinese (zh)
Other versions
CN1804655A (en
Inventor
周浩
文必洋
吴世才
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Wuhan University WHU
Original Assignee
Wuhan University WHU
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Wuhan University WHU filed Critical Wuhan University WHU
Priority to CNB200610018265XA priority Critical patent/CN100472223C/en
Publication of CN1804655A publication Critical patent/CN1804655A/en
Application granted granted Critical
Publication of CN100472223C publication Critical patent/CN100472223C/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

The invention relates to a high frequency radar anti-radio frequency interference method, which detects the interference signal on the radar collected far distance element accumulating Doppler spectrum by the setting signal-to-noise ratio gate and adopts different scan frequency impulse phase modulating modes to reduce the bad effect of the detecting result by the radio frequency interference distribution received by the high frequency radar.

Description

Anti-RF interference method for high frequency radar
Technical field
The present invention relates to the method for the anti-Radio frequency interference (RFI) of a kind of higher-frequency radar.
Background technology
High frequency (3-30MHz) radars are because of its unique over the horizon and round-the-clock detectivity, and obtained broad research and be successfully applied to during the ocean surface kinetic parameter measures, but and moving targets such as real-time detection naval vessel, aircraft and guided missile.The Radio frequency interference (RFI) of generations such as intensive short wave communication and broadcasting station can have a strong impact on the detection performance of radar in the high band, and this is the operating subject matter of higher-frequency radar.In existing higher-frequency radar anti-interference method, the adaptive beam that adopts forms anti-interference method and depends on large phased-array usually, and system is comparatively complicated, and is not suitable for small-bore widebeam radar; And from the time domain of signal or the anti-interference method of frequency domain character, as based on detection-amputation-interpolation method of disturbing instantaneous feature with based on the undesired signal subspace orthographic projection of interference distance distribution characteristics etc., all need the received signal that contains Radio frequency interference (RFI) is carried out Filtering Processing, thereby inevitably can to useful signal cause to a certain degree loss and distortion (as signal spectrum main lobe broadening, secondary lobe raises), and when having a plurality of interference, this influence is then even more serious.
Summary of the invention
The purpose of this invention is to provide a kind of anti-RF interference method for high frequency radar, avoid the interference of existing filter processing method to suppress not exclusively to reach to adverse effects such as signal generation distortion, obtain interfere information more accurately, recover disturbance spectrum fully, make the useful signal randomization, reduce to disturb and the influencing each other of useful signal, realize the independent analysis of the two, for higher-frequency radar provides a kind of simple and direct effective and cheap radio frequency interference resistant method, to improve the especially detection performance of portable high frequency radar system of higher-frequency radar.
Technical scheme of the present invention is: a kind of anti-RF interference method for high frequency radar is characterized in that:
On the remote unit accumulation doppler spectral of radar collection, detect undesired signal by default signal-noise ratio threshold;
Different scanning frequency pulse phase modulation methods is adopted in the distribution of the Radio frequency interference (RFI) that receives according to higher-frequency radar, realizes anti-interference:
When disturbing when frequency domain be intensive distribution, adopt at random pm mode with in disturbance spectrum albefaction to the average substrate;
When disturbing when frequency domain is sparse distribution, adopt linear pm mode that disturbance spectrum is moved outside the interested spectrum district, recover concealed signal;
When disturb do not exist or its influence little during to admissible degree, adopt the constant phase mode, signal is not made any extra process.
Aforesaid method is characterized in that: modulation phase sequence is compensated, recover disturbance spectrum fully, make the useful signal randomization, reduce to disturb and the influencing each other of useful signal, realize the independent analysis of the two.
In interference distance-doppler spectral that higher-frequency radar receives, the formation of interference stripes is actually owing to the chirp waveform of radar emission in the relevant processing time is in full accord, thereby the formation of disturbance spectrum is similar to a kind of process of Radio Frequency Interfere being carried out extraction with radar pulse repetition frequency, disturbs to have obtained coherent accumulation.Based on this, the present invention considers the radar anti-interference problem along another kind of thinking, and promptly the chirp to radar emission carries out phase modulation (PM), with the spectrum form of change interference, thereby realizes anti-interference.Concentrated Radio frequency interference (RFI) at facing usually in the higher-frequency radar can adopt linear pm mode, and disturbance spectrum is moved to interested frequency range (as wave single order and second order spectrum district); And at expansion profile Radio frequency interference (RFI), then can adopt the mode of phase modulation at random, with disturbance spectrum randomization (albefaction) on a uniform substrate level.Because radar adopts coherent demodulation, local oscillation signal with transmit consistent, therefore phase modulation without any influence, but for radar provides convenient, antijamming measure cheaply, is particularly useful for being used for the portable high frequency groundwave radar of extra large attitude monitoring to useful echoed signal.Especially, adopt another benefit of pulse-phase modulation method to be, if carry out phase compensation to received signal, then can the former disturbance spectrum of reconstruct, and with useful echo spectrum randomization, thereby can obtain more accurate interference parameter estimation.Useful signal and undesired signal can separate processes, and influencing each other to each other reaches minimum.And do not exist or its intensity is lower than when allowing level when disturbing, can adopt constant initial phase frequency modulation on pulse again, signal is not made extra process.The present invention can be effectively applied in the anti-interference process of higher-frequency radar flexibly, thereby has improved the work efficiency and the detection performance of radar greatly.
Advantage of the present invention is: need not to be provided with any auxiliary antenna, make anti-interference more cheap; Do not have complicated signal processing, calculated amount is very little, and disposal route is very convenient; Disturbing number more after a little while, by linear phase modulation can move out disturbance spectrum fully and to useful signal without any loss, thereby avoided the interference of existing filter processing method to suppress to produce adverse effects such as distortion not exclusively and to signal; On the accumulation doppler spectral in the remote unit, detect and disturb, can obtain interfere information more accurately; Modulation phase sequence is compensated, can recover disturbance spectrum fully, and make the useful signal randomization, reduce influencing each other of interference and useful signal, realize the independent analysis of the two, and the accurate analysis of disturbing helps further optimizing the detection performance of radar system.
Description of drawings
Fig. 1, the higher-frequency radar fundamental diagram of the embodiment of the invention.Wherein, 1 receiving antenna, 2 local oscillation signals, 3 separate range conversion, 4 Doppler-shifts
Fig. 2, the accumulation doppler spectral of concentrating profile to disturb.
Fig. 3, the accumulation doppler spectral that the expansion profile disturbs.
Fig. 4, means for anti-jamming principle of work block diagram.Wherein, 11 transmit, 12 local oscillation signals, and 13 received signals, 14 target analysis of spectrum modules, 15 disturbance spectrum analysis modules, 16 pm modes are selected module, 17 follow-up signal processing modules, * represents conjugation.
Embodiment
Below in conjunction with drawings and Examples, the present invention is done more detailed explanation.
The principle of work of linear frequency modulation system higher-frequency radar and the feature of disturbance spectrum simply are described below:
The principle of work block diagram of high-frequency ground wave radar as shown in Figure 1.Radar adopts the linear frequency modulation system, in each frequency sweep cycle to demodulation after the signal discrete Fourier transform (DFT) of carrying out low-pass filtering, sampling and fast time domain obtain the distance spectrum of this frequency sweep cycle, wherein each spectrum is o'clock corresponding to a spectrum value sample on the distance element, in time a plurality of distance spectrum sequences are carried out slow time domain discrete Fourier transform at the coherent accumulation of being made up of a plurality of frequency sweep cycles and obtain distance-Doppler's two-dimensional spectrum (as shown in Figure 1, wherein 1 is receiving antenna, 2 is local oscillation signal, 3 for separating range conversion, and 4 is Doppler-shift).Distance-doppler spectral by signal on a plurality of receiving antenna passages just can extract extra large attitude information and moving target information.
Linear FM radar transmit and local oscillation signal can be expressed as
g ( t ) = e j 2 &pi; ( f 0 t + 1 2 kt 2 ) , - T / 2 &le; t < T / 2 - - - ( 1 )
F wherein 0, k and T be respectively radar center frequency, chirp slope and frequency sweep cycle, t is fast time variable.
The spectrum form that amplitude modulation (AM) the type communication that following analysis of high frequency radar faces is usually disturbed.N (n=0 ..., N-1) in the individual frequency sweep cycle, interference can be expressed as
Figure C200610018265D00052
S wherein c(t ') and s m(t ') is respectively carrier signal and modulation signal, a, f lWith
Figure C200610018265D00061
Be respectively amplitude, frequency and the first phase of carrier wave.T '=nT+t is common time variable.
Mixing output undesired signal is s (t ')=s i(t ') g (t).The interior distance spectrum of n frequency sweep cycle that forms by baseband low pass filters H (f) back is
S n ( f ) = F [ s ( t &prime; ) ] H ( f ) = F [ s i ( t &prime; ) g ( t ) ] H ( f ) = F [ s c ( t &prime; ) g ( t ) + s m ( t &prime; ) g ( t ) ] H ( f )
(3)
= F [ s c ( t &prime; ) g ( t ) ] H ( f ) + F [ s m ( t &prime; ) g ( t ) ] H ( f ) = def S c n ( f ) + S m n ( f ) ,
F[wherein] Fourier transform of expression t when fast, subscript " expression frequency sweep cycle sequence number.Here disturbance spectrum is made up of definite carrier signal spectrum and Stochastic Modulation signal spectrum two parts.
Have by formula (2)
Figure C200610018265D00064
Therefore the interfering carrier spectrum is
S c n ( f ) = F [ s c ( t &prime; ) g ( t ) ] H ( f ) = { F [ s c ( t ) g ( t ) ] H ( f ) } e j 2 &pi; f l nT = S c 0 ( f ) e j 2 &pi; f l nT . - - - ( 5 )
This shows the characteristic of interfering carrier spectrum, it is a sinusoidal signal in slow time domain (is variable with n), thereby presents the line spectrum form in the Doppler territory, and this form and range-independence, and disturbance spectrum will appear on all distance elements.Definitely, interfering carrier is rendered as the bright wisp band that is parallel to distance axis on distance-doppler spectral, and its frequency location is
f D = mod ( f l , f r ) , mod ( f l , f r ) &le; f r 2 mod ( f l , f r ) - f r , mod ( f l , f r ) > f r 2 , - - - ( 6 )
F wherein r=1/T is the pulse repetition rate of radar.When disturbing phase place to exist fluctuation when (as owing to carry the ship motion or ionosphere changes), disturb then broadening and become strips to some extent of spectral line.
See the modulation signal part of interference again.Usually the modulation signal bandwidth is 5-10kHz, and its correlation time (inverse that is about bandwidth) is less than 1ms, and frequency sweep cycle is generally hundreds of ms, so can be considered the independent noise signal at different frequency sweep cycle modulation signals.After Doppler's processing, the power ratio of carrier wave and modulation signal will improve, and its gain equals scanning frequency pulse and counts N.Under this meaning, the modulation composition of interference only improves the generation contribution to the noise floor of echoed signal, and carrier wave then becomes the principal ingredient of disturbing the radar detection performance, therefore, and the work that the inhibition of interfering carrier composition is normally even more important.
As the above analysis, to present line spectrum be because obtained coherent accumulation in different frequency sweep cycles to interfering carrier.If each scanning frequency pulse parameter is all identical, will be equivalent to radar pulse repetition frequency the processing of interfering carrier so extraction will be carried out in interference, interfering frequency also will fold in the doppler analysis bandwidth on the corresponding position.Yet, if the parameter of scanning frequency pulse is time dependent, so just can carry out the transformation of certain mode to disturbance spectrum, avoid or reduce its interference echo spectrum interested district.As in extra large attitude is surveyed, interested is positive and negative single order Bragg spectrum of wave and second order spectrum district, and usually when extracting ocean current and unrestrained high information, generally only utilize the higher side spectrum district of signal to noise ratio (S/N ratio), thereby can consider disturbance spectrum moved outside this useful spectrum district and remove its adverse effect.Investigate the disturbance spectrum forming process, as can be known transponder pulse is carried out phase modulation, change the initial phase of each scanning frequency pulse, just can realize that disturbance spectrum moves.Suppose that the disturbance spectrum amount of moving is Δ f, then the respective phase modulation sequence is α (n)=2 π Δ fnT, and the transponder pulse behind the phase modulation is
g m(t′)=g(t)e ja(n), (7)
Use g m(t ') replaces g (t) substitution formula (3), and utilizes the character of Fourier transform, has
S m n ( f ) = F [ s i ( t &prime; ) g m ( t &prime; ) ] H ( f ) = F [ s i ( t &prime; ) g ( t ) e j 2 &pi;f nT ] H ( f ) = S n ( f - &Delta;f ) . - - ( 8 )
By selecting the suitable amount of moving, just can remove the influence of interference effectively, this linear phase modulation method can be effectively applied to disturb daytime the less situation about concentrating when distributing that is.Yet, at some more under the condition of severe, as between the lights or be intensive extended pattern disturbing evening and distribute, disturb the application of moving method then to be restricted, can adopt the method for random phase modulation this moment, with disturbance spectrum randomization (or albefaction), thus the inhibition that realization is disturbed.This moment, the phase modulation (PM) sequence α (n) in formula (7) elected any random series as, and in formula (8)
Figure C200610018265D00072
Also by correspondingly randomization, because the spread spectrum interfering carrier is attenuated N doubly, letter (number) do (disturbing) than corresponding acquisition gain N then.
If to adopting the received signal behind the phase-modulated pulse at random to carry out phase compensation, the signal times that is about to n frequency sweep cycle is with e -ja (n), then can recover former disturbance spectrum, realize independent analysis to useful signal and undesired signal.Useful signal is because of being subjected to the random phase modulation by albefaction after phase compensation, and jam-to-signal obtains gain N, helps extracting more accurately interfere information, and this information can be used for subsequent treatment to obtain more excellent system's detection performance.
Key of the present invention is to detect undesired signal in the accumulation doppler spectral in remote unit, and it is anti-interference to realize to select for use suitable manner to carry out phase modulation according to the distribution characteristics of disturbing then.
In typically containing the echo distance-Doppler spectrogram of Radio frequency interference (RFI), useful signal is concentrated and is distributed in the nearer distance element section, intensity increases with distance and reduces, and Radio frequency interference (RFI) is the bright wisp zonal arrangement that is parallel to distance axis in distance-doppler spectral, its intensity is comparatively steady on different distance, therefore, and by rational radar parameter design, a larger distance section can be provided, make not contain useful signal on this distance segment and only contain undesired signal.
Selection obtains disturbance spectrum than carrying out the doppler spectral accumulation between far region.Typical concentrated profile disturbance spectrum wherein disturbs the unimodal form to separate to occur as shown in Figure 2.Typical expansion profile disturbance spectrum disturbs the spectrum district form that then distributes with expansion to occur as shown in Figure 3.By default signal-noise ratio threshold, can detect the existence and the frequency separation thereof of Radio frequency interference (RFI), and judge the distribution pattern of disturbing by interfering frequency length of an interval degree.
After judging the interference profile type, select corresponding pm mode to carry out the setting of scanning frequency pulse initial phase.For pm mode at random, with first phase α (n) (n=0 ... N-1) be made as random series arbitrarily, then it is made as α (n)=2 π Δ fnT for linear pm mode,, then establishes it and be arbitrarily normal value sequence if do not have interference or interference strength less than permissible range.
Modulation phase sequence is compensated, and the signal times that is about to n frequency sweep cycle is with e -ja (n), recover disturbance spectrum fully, obtain the accurate information of interference, for follow-up further processing.
The anti-interference work block diagram of phase modulation as shown in Figure 4.Wherein 11 is to transmit, the 12nd, and local oscillation signal, the 13rd, received signal, the 14th, target analysis of spectrum module, the 15th, the disturbance spectrum analysis module, the 16th, pm mode is selected module, and the 17th, the follow-up signal processing module, * represents conjugation.

Claims (2)

1, a kind of anti-RF interference method for high frequency radar is characterized in that:
On the remote unit accumulation doppler spectral of radar collection, detect undesired signal by default signal-noise ratio threshold;
Different scanning frequency pulse phase modulation methods is adopted in the distribution of the Radio frequency interference (RFI) that receives according to higher-frequency radar, realizes anti-interference:
When disturbing when frequency domain be intensive distribution, adopt at random pm mode with in disturbance spectrum albefaction to the average substrate;
When disturbing when frequency domain is sparse distribution, adopt linear pm mode that disturbance spectrum is moved outside the interested spectrum district, recover concealed signal;
When disturb do not exist or its influence little during to admissible degree, adopt the constant phase mode, signal is not made any extra process.
2, the method for claim 1 is characterized in that: modulation phase sequence is compensated, recover disturbance spectrum fully, make the useful signal randomization, reduce to disturb and the influencing each other of useful signal, realize the independent analysis of the two.
CNB200610018265XA 2006-01-20 2006-01-20 Anti-RF interference method for high frequency radar Expired - Fee Related CN100472223C (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CNB200610018265XA CN100472223C (en) 2006-01-20 2006-01-20 Anti-RF interference method for high frequency radar

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CNB200610018265XA CN100472223C (en) 2006-01-20 2006-01-20 Anti-RF interference method for high frequency radar

Publications (2)

Publication Number Publication Date
CN1804655A CN1804655A (en) 2006-07-19
CN100472223C true CN100472223C (en) 2009-03-25

Family

ID=36866739

Family Applications (1)

Application Number Title Priority Date Filing Date
CNB200610018265XA Expired - Fee Related CN100472223C (en) 2006-01-20 2006-01-20 Anti-RF interference method for high frequency radar

Country Status (1)

Country Link
CN (1) CN100472223C (en)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1924607B (en) * 2006-09-14 2010-05-12 武汉大学 Method for controlling multiple working conditions high frequency radar receiver
CN101089653B (en) * 2007-07-20 2011-03-09 西安理工大学 Short-range frequency-modulation continuous wave FMCW radar anti-interference method
CN101334469B (en) * 2008-08-04 2011-01-12 北京理工大学 Wind profile radar clutter suppression method based on fraction order Fourier transform
CN101581782B (en) * 2009-06-15 2012-05-23 武汉大学 Method for inhibiting ionospheric clutter in portable high frequency groundwave radar
CN102243300B (en) * 2011-04-15 2012-10-10 中国人民解放军国防科学技术大学 Method for radio frequency interference suppression and error correction of low-frequency synthetic aperture radar
CN103954944A (en) * 2014-05-14 2014-07-30 武汉大学 Radio-frequency interference suppression method of high-frequency ground wave radar
WO2016165281A1 (en) * 2015-04-16 2016-10-20 深圳大学 Radar communication method and system
EP3098623A1 (en) 2015-05-25 2016-11-30 Autoliv Development AB A vehicle radar system
US11002829B2 (en) * 2016-04-15 2021-05-11 Mediatek Inc. Radar interference mitigation method and apparatus
CN106680790B (en) * 2016-12-29 2019-05-17 中国电子科技集团公司第五十四研究所 A kind of pseudo-random code binary-phase modulation interference detection method
CN113267751A (en) * 2021-06-29 2021-08-17 珠海上富电技股份有限公司 Anti-interference method for vehicle-mounted millimeter wave radar

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN2441147Y (en) * 2000-10-20 2001-08-01 武汉大学 Ocean echo wave simulgating and marking device
CN2443543Y (en) * 2000-10-20 2001-08-15 武汉大学 Low-level wideband electronic phase sifter

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN2441147Y (en) * 2000-10-20 2001-08-01 武汉大学 Ocean echo wave simulgating and marking device
CN2443543Y (en) * 2000-10-20 2001-08-15 武汉大学 Low-level wideband electronic phase sifter

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
利用时频分析抑制射频干扰. 许瑜昀,文必洋,吴雄斌.无线电工程,第35卷第11期. 2005 *
基于短时傅立叶变换的高频地波雷达射频干扰抑制. 许瑜昀,文必洋,吴雄斌,罗郁杨.武汉大学学报(理学版),第50卷第5期. 2004 *
应用时频分析进行高频雷达射频干扰抑制. 周浩,文必洋,吴世才,刘晓峰.电子学报,第32卷第9期. 2004 *
高频海态雷达软件无线电接收机的设计. 万显荣,杨子杰,文必洋,柯亨玉.数据采集与处理,第19卷第4期. 2004 *
高频线性调频雷达的高精度目标瞬时距离估计. 周浩,文必洋,吴世才,马志刚.电子学报,第33卷第4期. 2005 *
高频雷达在射频干扰下的目标探测研究. 周浩,文必洋,吴世才.现代雷达,第26卷第9期. 2004 *
高频雷达接收机模拟前端的设计与实现. 张国军,文必洋,吴雄斌,白立云.无线电工程,第34卷第9期. 2004 *

Also Published As

Publication number Publication date
CN1804655A (en) 2006-07-19

Similar Documents

Publication Publication Date Title
CN100472223C (en) Anti-RF interference method for high frequency radar
US7151483B2 (en) System and method for concurrent operation of multiple radar or active sonar systems on a common frequency
CN103116155A (en) Homotype radar same frequency interference suppression method used for ship formation condition
US9075138B2 (en) Efficient pulse Doppler radar with no blind ranges, range ambiguities, blind speeds, or Doppler ambiguities
CN110109117B (en) Satellite-borne synthetic aperture radar convolution modulation interference method based on frequency control array
CN103728593B (en) A kind of method realizing ground wave OTHR simultaneously multifrequency transmitting/receiving
CN113093123B (en) Jammer for resisting pulse Doppler radar and interference method thereof
CN101825707A (en) Monopulse angular measurement method based on Keystone transformation and coherent integration
CN103954944A (en) Radio-frequency interference suppression method of high-frequency ground wave radar
CN102890272A (en) Method for processing millimeter wave cloud radar signal
CN108414966A (en) A kind of wideband correlation direction-finding system and method based on time-modulation
CN104914429B (en) Target indication radar system capable of adaptively selecting waveform according to target distance
CN110398718B (en) Radio frequency interference suppression method based on FRFT frequency estimation subspace
CN114063023A (en) Device and method for interfering SAR radar
US7081850B2 (en) Coherent detection of ultra wideband waveforms
CN114355310A (en) Multi-mode agile waveform generation and processing method of cognitive radar
Xianrong et al. Adaptive cochannel interference suppression based on subarrays for HFSWR
CN112881982B (en) Method for restraining dense forwarding interference by frequency agile radar
Tian et al. Fully digital multi‐frequency compact high‐frequency radar system for sea surface remote sensing
Shi et al. A novel ionospheric oblique-incidence sounding network consisting of the ionospheric oblique backscatter sounder and the parasitic oblique-incidence sounder
Wu et al. A novel range detection method for 60GHz LFMCW radar
CN102998659B (en) Doppler frequency spectrum shaping method and system based on interpulse modulation
Hakobyan et al. Sweep-based spectrum sensing method for interference-aware cognitive automotive radar
Dzvonkovskaya et al. Software-improved range resolution for oceanographic HF FMCW radar
RU2596229C1 (en) Method for increasing range resolution of radar station

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
C17 Cessation of patent right
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20090325

Termination date: 20110120