CN100525109C - Direct digital frequency synthesizer of shaped frequency noise - Google Patents

Direct digital frequency synthesizer of shaped frequency noise Download PDF

Info

Publication number
CN100525109C
CN100525109C CNB2005101171751A CN200510117175A CN100525109C CN 100525109 C CN100525109 C CN 100525109C CN B2005101171751 A CNB2005101171751 A CN B2005101171751A CN 200510117175 A CN200510117175 A CN 200510117175A CN 100525109 C CN100525109 C CN 100525109C
Authority
CN
China
Prior art keywords
frequency
shaping
accumulator
output
sine
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
CNB2005101171751A
Other languages
Chinese (zh)
Other versions
CN1960182A (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.)
Institute of Semiconductors of CAS
Original Assignee
Institute of Semiconductors of CAS
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 Institute of Semiconductors of CAS filed Critical Institute of Semiconductors of CAS
Priority to CNB2005101171751A priority Critical patent/CN100525109C/en
Publication of CN1960182A publication Critical patent/CN1960182A/en
Application granted granted Critical
Publication of CN100525109C publication Critical patent/CN100525109C/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Analogue/Digital Conversion (AREA)

Abstract

The invention relates to a frequency noise reforming direct digital frequency synthesizer, which can raise the dynamic feature of the output signal. It consists of the frequency reformer inner-inserter, the phase accumulator, and the sine or cosine inquiry table (SCIT), the D/A converter and the low-pass filter. The output of the reformer connects to the accumulator to complete accumulation. The accumulator produces the related phase value, which is sent to SCIT. The output of SCIT goes through the D/A converter and the low-pass filter. The final output is a low noise sine or cosine wave.

Description

A kind of direct digital frequency synthesizer of shaped frequency noise
Technical field
The present invention relates to the frequency synthesizer technical field, particularly a kind of direct digital frequency synthesizer of shaped frequency noise that is used to improve the output signal dynamic characteristic.
Background technology
Direct Digital Frequency Synthesizers (DDFS) is a kind of important low-cost high resolution frequency synthetic method.Different with the synthesizer based on phase-locked loop, DDFS and digital CMOS process are compatible fully, often are embedded into base band ASIC transceiver.Though DDFS has come out much years, because the restriction of Digital Logical Circuits and digital to analog converter performance can only be employed with narrow band frequency synthetic.Along with the improvement of Digital Logical Circuits and digital to analog converter technology, DDFS has become the synthetic main flow of wideband frequency recently.As shown in Figure 1, traditional DDFS comprises a phase accumulator 30 that generates phase value based on incoming frequency word 10.DDFS uses the mode that searches sine or cosine question blank 40 phase value is converted to sinusoidal magnitude value, and its restricted width is in the resolution of digital to analog converter 50.The anti-spike filter 60 that is added in digital to analog converter 50 back will filter out the glitch noise in data conversion process.
But the structures shape of DDFS itself have a large amount of clutters in its output signal, the source of these clutters comprises: as phase-accumulated 20 truncation noise of sine table address, be stored in the sample magnitude truncation noise in the sine table 40, the noise of digital to analog converter 50 and the phase noise of clock.Wherein the influence that brings with phase truncation again is the most serious.Because the periodicity of DDFS itself, the process that has caused phase truncation also are periodic, the phase truncation noise of Yin Ruing is very obvious on some Frequency point thus.
The method of multiple inhibition phase truncation noise has been proposed so far.The most direct scheme is exactly the length that increases the phase place word, but can cause the rapid increase of synthesizer hardware.And along with the increase of phase place word, the figure place of phase accumulator 30 also increases thereupon, will influence the speed that phase accumulator 30 adds up like this.
Summary of the invention
The objective of the invention is to, a kind of direct digital frequency synthesizer of shaped frequency noise (DDFS) that is used to improve the output signal dynamic characteristic is provided, eliminated because the quantizing noise that phase truncation caused.This method is by the MATLAB simulating, verifying, and realizes on chip.The test result of experiment chip shows: compare with the DDFS that does not add frequency noise shaping interpolater, the dynamic range of no parasitic signal (SFDR) has improved 10dB, and chip area does not have to increase substantially.This invention also can reduce the ROM area of DDFS under the situation that does not reduce SFDR, and ROM accounts for the very most of of whole DDFS area.
A kind of direct digital frequency synthesizer of shaped frequency noise of the present invention, comprise: frequency shaping device interpolater (20), phase accumulator (30), sine or cosine question blank (40), digital to analog converter (50) and low pass filter (60), wherein, frequency shaping device interpolater (20) comprises shaping accumulator (21) and delay circuit (22), outside N bit frequency word is input to the frequency shaping accumulator (21) in the frequency shaping device interpolater (20), and the output of shaping accumulator (21) is divided into as the high p position of phase accumulator input with as intercepted low N-p position; Wherein, the output of low N-p position is connected to delay circuit (22), returns again through delay circuit (22) computing output and is connected to shaping accumulator (21), revises the high p position as the phase accumulator input, thereby realizes the shaping of noise; The output of shaping is connected to phase accumulator (30) and adds up, corresponding phase value adds up out, these phase value input sine or cosine question blanks (40), digital to analog converter (50) is passed through in the output of sine or cosine question blank (40) again and low pass filter (60) is finally exported low noise sine or cosine wave;
Wherein the transfer function of delay circuit is 1-(1-Z -1) n, wherein n is the exponent number of frequency noise shaping interpolater.
1, the synthesizer step-length is f Clk/ 2 nIn order to obtain higher frequency resolution, use the accumulator of longer figure place usually, for example 16.Yet big accumulator needs big with 2 nQuestion blank for the address.The ROM question blank has just occupied the major part of whole DDFS area like this.In order to reduce the area of ROM question blank, phase value was blocked a part usually before being used to the question blank address.This blocks the introducing quantizing noise, and the noise that this quantizing noise can be described as a linearity is added on the sinusoidal wave phase place.Signal to noise ratio snr by the phase truncation introducing p=6.02p-3.992dB, wherein p is the width of phase place word.
2, the output valve width of ROM question blank is subjected to the restriction of the figure place of digital to analog converter.Only with the byte of finite width represent sine wave must introduce quantizing noise be added to output signal on.Similar to the phase truncation noise, because the signal to noise ratio snr that limited output amplitude width is introduced A=6.02D+1.76dB, wherein D is the resolution of digital to analog converter.
Consider because phase truncation e pThe quantizing noise and the amplitude that produce are blocked (limited ROM value width) e A, suppose the phase quantization noise much smaller than phase place, this complete DDFS output can be expressed as:
A out = A sin ( 2 πWi 2 n + e p ( i ) ) + Ae A ( i ) ≈ A sin ( 2 πWi 2 n ) + Ae p ( i ) cos ( 2 πWi 2 n ) + Ae A ( i )
Therefore as can be seen, phase noise is output the integration of signal and has modulated on amplitude.Above-mentioned equation provides simple DDFS output and quantization noise model.
Overlapping for avoiding in data conversion process occurring, the frequency that is synthesized is less than the clock frequency of DDFS.Therefore have over-sampling in DDFS, this just provides possibility for noise shaping.Noise shaping can be transferred to high-frequency region to phase noise, makes things convenient for low pass filter filters out.As shown in Figure 2, the numerical frequency word that accumulator produces corresponding frequency synthesis to being used to of receiving is summed into the digital value of a N position, high p position in the digital value of this N position is as the phase value of sine or cosine, and low N-p position turns back to accumulator again through delay circuit and adds up.The transfer function of delay circuit is 1-(1-Z -1) n, wherein n is the exponent number of frequency noise shaping interpolater.Process is hanged down the delay computing of N-p position and is returned again and add up, and revises the high p position as the phase value of sine or cosine.At this moment DDFS output can be expressed as:
A out = A sin ( 2 πWi 2 n + e p ( i ) ( 1 - z - 1 ) k ) + Ae A ( i )
≈ A sin ( 2 πWi 2 n ) + Ae p ( i ) ( 1 - z - 1 ) k cos ( 2 πWi 2 n ) + Ae A ( i )
Phase noise has been added a high pass filter as can be seen from the above equation.This algorithm greatly reduces near the phase truncation noise the output frequency.Phase noise after the high pass shaping can be by the low pass filter filters out behind the digital to analog converter.Therefore the noise of being introduced by phase truncation is greatly diminished.
Shaping is exported high p position and is connected to sine or cosine table through the frequency shaping interpolater, converts digital sine or cosine amplitude to by sine or cosine table.This digital sine or cosine amplitude convert analog sine or cosine signal to through digital to analog converter and low pass filter again.
Description of drawings
For further specifying technology contents of the present invention, be described in detail as follows below in conjunction with example and accompanying drawing, wherein:
Fig. 1 is the structured flowchart of conventional digital direct synthesizer;
Fig. 2 is the structured flowchart that the present invention has increased the digital direct frequency synthesizer of frequency shaping interpolater;
Fig. 3 is the instance graph of one 4 order frequency shaping interpolater of design.
Fig. 4 is the instance graph of one 5 order frequency shaping interpolater of design.
Fig. 5 is the chip photo in conjunction with a Direct Digital Frequency Synthesizers of the structural design of Fig. 2 and Fig. 3.
Fig. 6 is integrated and the contrast figure of the spectral characteristic of the frequency synthesis of integrated frequency shaping not.
Embodiment
According to Fig. 1, the output of a desirable Direct Digital Frequency Synthesizers can be expressed as:
A out = A sin ( 2 πWt 2 n T clk ) = A sin ( 2 πWi 2 n )
Wherein A is the full amplitude output valve of digital to analog converter, time variable t=T ClkI.The output here is desirable sine wave, yet owing to certainly existing noise jamming in the following reason output wave signal.
The present invention is a kind of direct digital frequency synthesizer of shaped frequency noise, as shown in Figure 2, comprise: frequency shaping device interpolater (20), phase accumulator (30), sine or cosine question blank (40), digital to analog converter (50) and low pass filter (60), wherein, frequency shaping device interpolater (20) comprises shaping accumulator (21) and delay circuit (22), outside N bit frequency word is input to the frequency shaping accumulator (21) in the frequency shaping device interpolater (20), and the output of frequency shaping accumulator (21) is divided into as the high p position of phase accumulator input with as intercepted low N-p position; Wherein, the output of low N-p position is connected to delay circuit (22), and carrying out transfer function through delay circuit (22) is 1-(1-Z -1) nComputing, wherein n is the exponent number of noise shaping interpolater, delay circuit (22) computing output is returned again and is connected to shaping accumulator (21), revises the high p position as the phase accumulator input, thereby realizes the shaping of noise; The output of shaping is connected to phase accumulator (30) and adds up, corresponding phase value adds up out, these phase value input sine or cosine question blanks (40), digital to analog converter (50) is passed through in the output of sine or cosine question blank (40) again and low pass filter (60) is finally exported low noise sine or cosine wave.
The shaping accumulator (30) of interpolater is 16 bit accumulators, and wherein most-significant byte output is as the phase value of sine or cosine, and least-significant byte turns back to accumulator again through delay circuit and adds up.Shaping is exported high p position and is connected to sine or cosine table through the frequency shaping interpolater, converts digital sine or cosine amplitude to by sine or cosine table (40).This digital sine or cosine amplitude pass through digital to analog converter (50) again and low pass filter (60) converts analog sine or cosine signal to.
Described direct digital frequency synthesizer of shaped frequency noise, sine or cosine table are input 8 bit address, the sine or the cosine table of 12 amplitudes of output.
Described direct digital frequency synthesizer of shaped frequency noise, digital to analog converter (50) are 12 figure place weighted-voltage D/A converters.
Described direct digital frequency synthesizer of shaped frequency noise, low pass filter (60) are anti-spike low pass filter.
Described direct digital frequency synthesizer of shaped frequency noise, the shaping accumulator that is used for the frequency noise shaping interpolater of Direct Digital Frequency Synthesizers is 16 bit accumulators, wherein most-significant byte output is as the input of phase accumulator, and least-significant byte turns back to the shaping accumulator again through delay circuit and adds up.
By above analysis, two kinds of frequency shaping interpolaters of following specific design.In order to realize high-speed applications, adopted single cycle frequency shaping interpolater.Fig. 3 and Fig. 4 are requirement, 4 rank of design and the 5 order frequency shaping interpolaters according to different shaping degree.This structure has realized blocking (N-p) LSB (low-value) with equation 1-(1-z -1) kHigh pass noise shift " A " 24 expression (N-p) LSB among Fig. 3 and Fig. 4.P MSB (high-value) is used as the address of sine look up table 30.Because multiplier becomes the bottleneck of area and speed probably, shifting function is used to replace multiplier in conjunction with add operation.In Fig. 3 and Fig. 4 "<<" 25 representatives are to shifting left, then "<<2 " 25 representatives promptly are equivalent to take advantage of 4 to moving to left two.The frequency noise shaping interpolater of Fig. 3 and Fig. 4 adopts 16 inputs to add up, and 8 are blocked.
Now be that example specifies and how to realize blocking (N-p) LSB (low-value) with equation 1-(1-z with Fig. 3 -1) kHigh pass noise shift.Noise accumulator 22 adds up the output of the output valve of phase accumulator 20 and delay circuit 23 in Fig. 3, and output is divided into two parts, as the high p position of the phase value of sine or cosine with as intercepted low N-p position.Wherein low N-p position is through Z -1Postpone to obtain exporting AZ -1, AZ -1Be divided into three the tunnel.Wherein one road bit manipulation that moves to left promptly takes advantage of 2 operations to obtain 2AZ -1Output is connected to adder 26; Another road promptly takes advantage of 4 operations to obtain 4AZ through moving to left two bit manipulations -1Output, this output divide two-way wherein one the tunnel also to be connected to adder 26, export through the addition of adder and obtain 6AZ -1Output; Third Road postpones Z -2The AZ that attains the Way -3Output, this output also divides two-way wherein one the road to be connected to adder 28.The output 6AZ of adder 26 -1Also be divided into two-way, wherein one the tunnel also be connected to adder 28, adder 28 is 6AZ like this -1And AZ -3Addition obtains A (6Z -1+ Z -3) output, Z is passed through in this output again -1Postpone to obtain exporting A (6Z -2+ Z -4).AZ -3Another road of output promptly takes advantage of 4 operations to obtain 4AZ through moving to left two bit manipulations -3Output also is connected to adder 27, and 4AZ -1Another road of output also is connected to adder 27, so adder 27 is 4AZ -1And AZ -3Addition obtains 4A (Z -1+ Z -3) output.A (6Z -2+ Z -4) output gets complement code and be connected to adder 29 later on, 4A (Z -1+ Z -3) output also be connected to adder 29, obtain exporting A (4Z at last -1-6Z -2+ 4Z -3-Z -4), i.e. A (1-(1-Z -1) 4), then the shaping accumulator is returned in output again, has realized at last blocking (N-p) LSB (low-value) with equation 1-(1-z -1) kHigh pass noise shift.
Fig. 5 is the chip photo in conjunction with a direct digital frequency synthesizer of shaped frequency noise of the structural design of Fig. 2 and Fig. 3.As can be seen from the figure, the chip of this Direct Digital Frequency Synthesizers mainly comprises three parts: accumulator and 4 order frequency shaping interpolaters 70, sine look up table 60 and digital to analog converter 80.In accumulator and 4 order frequency shaping interpolaters, 70 modules, a control end is arranged, integrated or not integrated frequency shaping when being used for controlling frequency synthesis, be used for convenient contrast integrated with the spectral characteristic of the frequency synthesis of integrated frequency shaping not.
Fig. 6 is integrated and the not contrast of the spectral characteristic of the frequency synthesis of integrated frequency shaping, and obviously integrated as can be seen frequency shaping is better than not integrated spectral characteristic, and dynamic (SFDR) that promptly do not have parasitic signal wants big.The frequency spectrum of synthetic sine wave wherein, figure a does not add frequency shaping interpolater figure b and has added the frequency shaping interpolater.

Claims (7)

1, a kind of direct digital frequency synthesizer of shaped frequency noise, comprise: frequency shaping device interpolater (20), phase accumulator (30), sine or cosine question blank (40), digital to analog converter (50) and low pass filter (60), wherein, frequency shaping device interpolater (20) comprises shaping accumulator (21) and delay circuit (22), outside N bit frequency word is input to the frequency shaping accumulator (21) in the frequency shaping device interpolater (20), and the output of shaping accumulator (21) is divided into as the high p position of phase accumulator input with as intercepted low N-p position; Wherein, the output of low N-p position is connected to delay circuit (22), returns again through delay circuit (22) computing output and is connected to shaping accumulator (21), revises the high p position as the phase accumulator input, thereby realizes the shaping of noise; The output of shaping is connected to phase accumulator (30) and adds up, corresponding phase value adds up out, these phase value input sine or cosine question blanks (40), digital to analog converter (50) is passed through in the output of sine or cosine question blank (40) again and low pass filter (60) is finally exported low noise sine or cosine wave;
Wherein the transfer function of delay circuit is 1-(1-Z -1) n, wherein n is the exponent number of frequency noise shaping interpolater.
2, direct digital frequency synthesizer of shaped frequency noise according to claim 1 is characterized in that, is used for the frequency noise shaping interpolater of Direct Digital Frequency Synthesizers.
3, direct digital frequency synthesizer of shaped frequency noise according to claim 1, it is characterized in that, the shaping accumulator (30) that is used for the frequency noise shaping interpolater of Direct Digital Frequency Synthesizers is 16 bit accumulators, wherein most-significant byte output is as the input of phase accumulator, and least-significant byte turns back to the shaping accumulator again through delay circuit and adds up.
4, direct digital frequency synthesizer of shaped frequency noise according to claim 1 is characterized in that, phase-accumulated output phase value, and this phase value is connected to sine or cosine table.
5, direct digital frequency synthesizer of shaped frequency noise according to claim 1 is characterized in that, sine or cosine table (40) are input 8 bit address, the sine or the cosine table of 12 amplitudes of output.
6, direct digital frequency synthesizer of shaped frequency noise according to claim 1 is characterized in that, digital to analog converter (50) is 12 figure place weighted-voltage D/A converters.
7, direct digital frequency synthesizer of shaped frequency noise according to claim 1 is characterized in that, low pass filter (60) is anti-spike low pass filter.
CNB2005101171751A 2005-11-01 2005-11-01 Direct digital frequency synthesizer of shaped frequency noise Expired - Fee Related CN100525109C (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CNB2005101171751A CN100525109C (en) 2005-11-01 2005-11-01 Direct digital frequency synthesizer of shaped frequency noise

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CNB2005101171751A CN100525109C (en) 2005-11-01 2005-11-01 Direct digital frequency synthesizer of shaped frequency noise

Publications (2)

Publication Number Publication Date
CN1960182A CN1960182A (en) 2007-05-09
CN100525109C true CN100525109C (en) 2009-08-05

Family

ID=38071698

Family Applications (1)

Application Number Title Priority Date Filing Date
CNB2005101171751A Expired - Fee Related CN100525109C (en) 2005-11-01 2005-11-01 Direct digital frequency synthesizer of shaped frequency noise

Country Status (1)

Country Link
CN (1) CN100525109C (en)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN100520672C (en) * 2007-09-28 2009-07-29 电子科技大学 DDS signal source amplitude-frequency characteristic compensation method and related DDS signal source
CN101615051B (en) * 2008-06-25 2011-06-29 中国科学院半导体研究所 Compact direct digital frequency synthesizer used for system on chip
CN103095297B (en) * 2013-01-06 2015-03-25 中国电子科技集团公司第十研究所 Method for generating accurate frequency by direct digital frequency synthesizer
CN104734701A (en) * 2015-04-08 2015-06-24 中国科学院光电技术研究所 Low-stray DDS single-frequency signal generator
CN107782934A (en) * 2016-08-26 2018-03-09 江苏银佳电子设备有限公司 A kind of AC constant-current source for being used to measure accumulator internal resistance
CN108334157B (en) * 2017-01-20 2020-04-21 成都鼎桥通信技术有限公司 Carrier signal generation method and device

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4328554A (en) * 1980-07-03 1982-05-04 The United States Of America As Represented By The Secretary Of The Navy Programmable frequency synthesizer (PFS)
CN1076282A (en) * 1992-01-31 1993-09-15 莫托罗拉公司 Utilize the temperature compensation of Direct Digital synthetic crystal benchmark
US5430764A (en) * 1993-09-01 1995-07-04 Grand Valley State University Direct digital frequency synthesizer
US5986483A (en) * 1997-10-02 1999-11-16 National Science Council Direct digital frequency systhesizer
CN1469547A (en) * 2003-06-10 2004-01-21 湘潭师范学院 High-precision optional waveform generator based on FPGA
CN1472882A (en) * 2003-06-24 2004-02-04 深圳市东方汉华软件技术有限公司 Method and circuit structure for temperature compensation and frequency calibration of oscillation frequency source

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4328554A (en) * 1980-07-03 1982-05-04 The United States Of America As Represented By The Secretary Of The Navy Programmable frequency synthesizer (PFS)
CN1076282A (en) * 1992-01-31 1993-09-15 莫托罗拉公司 Utilize the temperature compensation of Direct Digital synthetic crystal benchmark
US5430764A (en) * 1993-09-01 1995-07-04 Grand Valley State University Direct digital frequency synthesizer
US5986483A (en) * 1997-10-02 1999-11-16 National Science Council Direct digital frequency systhesizer
CN1469547A (en) * 2003-06-10 2004-01-21 湘潭师范学院 High-precision optional waveform generator based on FPGA
CN1472882A (en) * 2003-06-24 2004-02-04 深圳市东方汉华软件技术有限公司 Method and circuit structure for temperature compensation and frequency calibration of oscillation frequency source

Also Published As

Publication number Publication date
CN1960182A (en) 2007-05-09

Similar Documents

Publication Publication Date Title
CN100525109C (en) Direct digital frequency synthesizer of shaped frequency noise
CN100568153C (en) But a kind of synchronous DDS device of binary channels phase-modulation and amplitude-modulation
CN105718404B (en) A kind of square-wave generator and method based on FPGA
CN100392976C (en) High order sigmatriangle noise shaping interpolator for direct digital frequency synthesis
JPH0548003B2 (en)
CN1333948A (en) Method and system for calibrating analog-to-digital conversion
CN1447935B (en) Method for generating clock signal and clock generator
CN101335509B (en) Method and digital control oscillator for sinusoidal and cosine signal generation
CN110488228B (en) Linear frequency modulation signal generation method and device and storage medium
WO2012093386A2 (en) Method and system for signal synthesis
CN103488245A (en) Phase-amplitude conversion method and device in DDS (Direct Digital Synthesizer)
CN104601148A (en) Digital voltage ramp generator
CN101615051B (en) Compact direct digital frequency synthesizer used for system on chip
CN102025375B (en) Analogue-to-digital converter and digital calibration circuit thereof
US8319563B2 (en) Digitally controlled oscillator
CN201509196U (en) Intermediate-frequency transceiver chip
CN104113333A (en) Direct digital frequency synthesizer
CN101854172B (en) Numerical control oscillator parallel design method based on two-dimensional sine table
JP4536946B2 (en) Frequency synthesizer
CN1783701A (en) High order sigma delta noise shaping direct digital frequency synthesizer
CN100557946C (en) Low power consumption non-ROM searching list phase amplitude converter
CN102006067B (en) DDS (Digital Display Scope) circuit structure with waveform correction ROM (Read Only Memory)
CN112771783B (en) Asynchronous sampling architecture and chip
CN103873409B (en) A kind of modulator, signal generator and method generating π/4-DQPSK modulated signal
CN102420611A (en) Sampling rate conversion method and device of digital signal

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: 20090805

Termination date: 20101101