CN201607527U - Intermediate-frequency direct-sequence spread spectrum receiver - Google Patents

Intermediate-frequency direct-sequence spread spectrum receiver Download PDF

Info

Publication number
CN201607527U
CN201607527U CN2009202780708U CN200920278070U CN201607527U CN 201607527 U CN201607527 U CN 201607527U CN 2009202780708 U CN2009202780708 U CN 2009202780708U CN 200920278070 U CN200920278070 U CN 200920278070U CN 201607527 U CN201607527 U CN 201607527U
Authority
CN
China
Prior art keywords
branch road
module
tracking loop
output
branch
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
CN2009202780708U
Other languages
Chinese (zh)
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.)
Beihang University
Original Assignee
Beihang University
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 Beihang University filed Critical Beihang University
Priority to CN2009202780708U priority Critical patent/CN201607527U/en
Application granted granted Critical
Publication of CN201607527U publication Critical patent/CN201607527U/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Abstract

The utility model discloses an intermediate-frequency direct-sequence spread spectrum receiver used for satellite ranging, comprising 37 parts, such as a front-end A/D (10), an FFT module (42), a local PN code generator (80), a correlator (81), an automatic threshold computation module (70) and other parts. The connection relation is as follows: the outputs of the front-end A/D (10) and a carrier tracking loop NCO (93) are respectively connected with an in-phase channel multiplier (20) and a quadrature channel multiplier (30), and the inputs thereof enters an in-phase channel FIR low-pass filter (21) and a quadrature channel FIR low-pass filter (31); on one hand, the result is output to an integral clearing device (40), and then sent to an FFT module (42), a channel 1 local PN code storage ROM (51) and a channel 2 local PN code storage ROM (61), next, the output enters a channel 1 plural multiplier (50) and channel 2 plural multiplier (60), then the output enters a channel 1 root-mean-square module (53) and a channel 2 root-mean-square module (63), and finally the output is transmitted to the automatic threshold computation module (70) and an acquisition and adjustment module (71) for code acquisition; and on the other hand, the result is output to the correlator (81) and the local PN code generator (80) for code tracking. The output of the correlator (81) is transmitted to a carrier tracking loop frequency discriminator/ phase discriminator (91) simultaneously, and then transmitted to a carrier tracking loop circuit filter (91), and the output of the carrier tracking loop circuit filter (91) enters a carrier tracking loop NCO (93) for carrier tracking.

Description

A kind of intermediate frequency direct sequence spread spectrum receiver
(1) technical field
The utility model relates to a kind of band spread receiver, relates in particular to a kind of intermediate frequency direct sequence spread spectrum receiver that is used for satellite ranging, and this invention belongs to technical field of satellite communication.
(2) background technology
Development along with spationautics, the survey of deep space strategy of China is formally launched, and bear first fruit, in October, 2007, the Chang'e I moon probing satellite succeeds in sending up, and has passed a large amount of moon telemetries back, for utilizing, China surveys space, utilize space to step the first step, simultaneously, the success of moon exploration will be established technical foundation for China carries out survey of deep space.
In the survey of deep space process, it is significant constantly to understand the residing exact position of satellite, location to satellite not only can make land station understand the duty of satellite, residing position constantly, and can provide foundation for change rail, the attitude adjustment of satellite.The principle of range finding is to measure the propagation delay time of radio signal to measured target, thereby calculates and target range.The satellite ranging system of China mainly is divided into two kinds, sidetone (continuant) range measurement system and pseudo-code (spread spectrum) range measurement system at present.Continuant range finding is main by sending the continuous distance measuring signal of a series of different frequencies, the time-delay of measuring each frequency signal that receives, and highest frequency component is called as the range finding clock, has determined the precision of side-tone ranging, and other component is used for fuzzy distance solution.The pseudo-random code ranging signal is by sending specific spread spectrum code sequence, and the phase place time-delay of the spread spectrum code sequence that instrumented satellite returns can obtain the distance between land station and the satellite.Compare with side-tone ranging, pseudo-random code ranging has the distance accuracy height, no fuzzy distance is big, antijamming capability is strong, and anti-fading ability is strong, and the anti-multipath interference performance is strong, good confidentiality, have advantages such as CDMA ability, pseudo-random code ranging is the important means of satellite observing and controlling, and China's satellite ranging mode just more and more adopts the pseudo-random code ranging mode at present.
The precision of pseudo-random code ranging is by the spreading rate decision of spreading code, and the high more distance measuring precision of spreading rate is high more, and the no fuzzy distance of range finding is by the length decision of ranging code.Ranging code is long more, and the distance range that can measure is big more.In order to improve distance accuracy, increase maximum measuring distance, need to improve the spreading rate of spreading code, and select longer spreading code, but spreading rate is high more, requirement to hardware is high more, the difficulty that realizes is big more, and the spreading code cycle is long more, and the difficulty of catching is big more, the time of catching is long more, therefore design a kind of can be significant to the band spread receiver that high spreading rate long code is caught fast.
(3) summary of the invention
1, purpose: the purpose of this utility model provides a kind of intermediate frequency direct sequence spread spectrum receiver that is used for satellite ranging.This band spread receiver can be caught fast to the long code of high spreading rate, and has higher distance accuracy.
2, technical scheme: as shown in Figure 1, the utility model is a kind of intermediate frequency direct sequence spread spectrum receiver that is used for satellite ranging, and composition comprises: the front end analog-to-digital conversion module (is called for short front-end A/D) 10, in-phase branch multiplier 20, in-phase branch finite impulse response low-pass filter (being called for short in-phase branch FIR low-pass filter) 21, quadrature branch multiplier 30, quadrature branch finite impulse response low-pass filter (being called for short quadrature branch FIR low-pass filter) 31, integration zero clearing device 40, buffer 41, fast Fourier transform module (being called for short the FFT module) 42, branch road 1 complex multiplier 50, branch road 1 local pseudo-random code storage ROM (read-only memory) (being called for short branch road 1 local PN sign indicating number storage ROM) 51, branch road 1 inverse fast fourier transform module (being called for short branch road 1IFFT module) 52, branch road 1 mean square root module 53, branch road 2 complex multipliers 60, branch road 2 local pseudo-random code storage ROM (read-only memory) (being called for short branch road 2 local PN sign indicating number storage ROM) 61, branch road 2 inverse fast fourier transform modules (being called for short branch road 2IFFT module) 62, branch road 2 mean square root module 63, automatic threshold computing module 70, catch judge module 71, pseudo-random code reference position computing module (being called for short PN sign indicating number reference position computing module) 72, carrier frequency adjustment module 73, local pseudo-random code maker (being called for short local PN code generator) 80, correlator 81, code tracking loop mean square root module 82, code tracking loop phase detector 83, code tracking loop loop filter 84, code tracking loop totalizer 85, code tracking loop digital controlled oscillator (being called for short code tracking loop NCO) 86, code tracking loop Doppler effect correction module 87, carrier tracking loop frequency discrimination/phase detector 90, carrier tracking loop loop filter 91, carrier tracking loop totalizer 92 and carrier tracking loop digital controlled oscillator (being called for short carrier tracking loop NCO) 93.
Wherein code tracking loop loop filter 84 is identical with carrier tracking loop loop filter 91 inner structures, has only the parameter difference, as shown in Figure 2: comprise straight-through branch amplifier 100, loop filter totalizer 101, integration branch amplifier 110, integration branch road totalizer 111 and integration branch road delay unit 112.
Above-mentioned all constituents, except that front-end A/D10 used the ready-made product of putting on the shelf, remainder was all realized in programmable gate array (FPGA).
Described front-end A/D10 carries out bandpass sampling with constant sampling rate to the 70M intermediate frequency spread-spectrum signal after nursing one's health, and the simulating signal of importing is transformed to digital signal, and front-end A/D10 uses ready-made product.
Carrier tracking loop NCO 93 uses direct frequency synthesizing algorithm (being called for short the DDS algorithm) to realize, be responsible for producing fixedly local carrier of the two-way identical with nominal intermediate frequency behind front-end A/D10 bandpass sampling, 90 ° of the phase phasic differences of two-way carrier wave, the signal behind the bandpass sampling of the output of carrier tracking loop NCO 93 and front-end A/D10 output is as the input of in-phase branch multiplier 20 and quadrature branch multiplier 30; In-phase branch multiplier 20 and quadrature branch multiplier 30 use the inner IP kernel of FPGA to realize, use as low-converter, input signal is down-converted to zero intermediate frequency, and in-phase branch multiplier 20 enters identical in-phase branch FIR low-pass filter 21 of structure and quadrature branch FIR low-pass filter 31 respectively with the result of quadrature branch multiplier 30 outputs; In-phase branch FIR low-pass filter 21 and quadrature branch FIR low-pass filter 31 adopt the FIR Structure Filter, use the inner IP kernel of FPGA to realize, be responsible for the signal after the down coversion is carried out filtering, frequency multiplication component and out-of-band noise after the filtering down coversion drop to the noise power in the baseband signal lower.
Result after the low-pass filtering exports to integration zero clearing device 40 on the one hand, carries out sign indicating number and catches, and exports to correlator 81 on the other hand and is used for code tracking.Integration zero clearing device 40 carries out the integration zero clearing respectively to in-phase branch FIR low-pass filter 21 and quadrature branch FIR low-pass filter 31 filtered results under the local PN sign indicating number clock control that code tracking loop NCO86 produces.The result of integration zero clearing device 40 outputs enters buffer 41, and the signal to homophase and quadrature branch carries out buffer memory simultaneously in buffer 41, and to adapt to the requirement of 42 pairs of input data rates of follow-up FFT module, buffer 41 uses the inner IP kernel of FPGA to realize.The homophase of exporting behind buffer 41 buffer memorys and the result of quadrature branch as the real part and the imaginary part input of FFT module 42, carry out Fast Fourier Transform (FFT) respectively in FFT module 42, FFT module 42 uses the inner IP kernel of FPGA to realize.The real part of FFT module 42 and imaginary part output are as one tunnel input of branch road 1 complex multiplier 50 and branch road 2 complex multipliers 60, and branch road 1 complex multiplier 50 and branch road 2 complex multipliers 60 use the inner IP kernel of FPGA to realize.Another road input of branch road 1 complex multiplier 50 and branch road 2 complex multipliers 60 is provided by branch road 1 local PN sign indicating number storage ROM51 and branch road 2 local PN sign indicating number storage ROM61 respectively.The output of FFT module 42 successively with local PN sign indicating number storage ROM51 of branch road 1 and branch road 2 local PN sign indicating numbers storage ROM61 in the local PN sign indicating number the stored result after through the FFT conversion carry out complex multiplication.Local PN sign indicating number storage ROM51 of branch road 1 and branch road 2 local PN sign indicating numbers storage ROM61 use the inner IP kernels of FPGA to realize, the data of inside solidification are the result after code phase differs half the local PN sign indicating number of the two-way sequence FFT conversion of FFT length.Result after branch road 1 complex multiplier 50 and branch road 2 complex multipliers 60 calculate sends into branch road 1IFFT module 52 respectively and branch road 2IFFT module 62 is carried out inverse fast fourier transform, and branch road 1IFFT module 52 and branch road 2IFFT module 62 use the inner IP kernel of FPGA to realize.Branch road 1 mean square root module 53 is sent in the output of branch road 1IFFT module 52 and branch road 2IFFT module 62 respectively and branch road 2 mean square root module 63 are calculated corresponding root-mean-square value, and the JPL approximate data is adopted in the calculating of root-mean-square value, promptly utilizes formula a 2 + b 2 ≈ Max ( abs ( a ) , abs ( b ) ) + 1 2 ( Min ( abs ( a ) , abs ( b ) ) ) Calculate root-mean-square value, calculate after the root-mean-square value, automatic threshold computing module 70 is sent in the output of branch road 1 mean square root module 53 and branch road 2 mean square root module 63 simultaneously, the state computation of relevant peaks goes out suitable threshold value according to this moment, output and branch road 1 mean square root module 53 with automatic threshold computing module 70, the output of branch road 2 mean square root module 63 is sent into simultaneously and is caught judge module 71, catches thresholding and branch road 1 mean square root module 53 of judge module 71 by relatively automatic threshold module 70 outputs, the correlation peak of the output of branch road 2 mean square root module 63 judges whether band spread receiver catches.Catch the output of judge module 71 and send into PN sign indicating number reference position computing module 72 and carrier frequency adjustment module 73 respectively.
When catching judge module 71 when judging that the spread spectrum answering machines are not caught, band spread receiver is operated in trapped state, and carrier frequency adjustment module 73 is being searched for all after dates of a sign indicating number, can adjust the frequency separation of output, thereby adjust the frequency of carrier tracking loop NCO 93 outputs, until finish catch till.PN sign indicating number reference position computing module 72 not outputs this moment.When catching judge module 71 when judging that the spread spectrum answering machines have been caught, will make band spread receiver change tracking mode over to by trapped state.At this moment, carrier frequency adjustment module 73 can stop to adjust, simultaneously code tracking loop Doppler effect correction module 87 is sent in output this moment, code tracking loop Doppler effect correction module 87 is according to the proportionate relationship between spreading rate and the carrier frequency, produce corresponding PN sign indicating number Doppler effect correction component, be used to adjust the sign indicating number clock of code tracking loop NCO 86 outputs, the compensating for doppler effect is to the influence of spreading code.And PN sign indicating number reference position computing module 72 will calculate the reference position of local PN by the position of when locking relevant peaks, and result calculated is sent into local PN code generator 80, adjust the phase place of the local PN sequence that local PN code generator 80 produces.
Local PN code generator 80 produce leading (PN-), on schedule (PN), (PN+) three tunnel local PN sign indicating numbers lag behind, this three tunnel local PN sign indicating number is the phase difference of half chip phase each other, and three tunnel local PN sign indicating numbers enter correlator 81 backs simultaneously and carry out related operation with the output of in-phase branch FIR low-pass filter 21 and quadrature branch FIR low-pass filter 31.Three road correlated results of output are sent into code tracking loop mean square root module 82 simultaneously and are calculated root mean square.Wherein correlated results need be sent into carrier tracking loop frequency discrimination/phase detector 90 on schedule, carries out being used for carrier track behind frequency discrimination, the phase demodulation.The result of code tracking loop mean square root module 82 outputs enters code tracking loop phase detector 83 phase demodulations, phase-demodulating principle be relatively leading relevant result and lag correlation the result, the phase place of determining code tracking loop NCO86 still lags behind in advance, and the result of phase demodulation sends into code tracking loop loop filter 84.The main effect of code tracking loop loop filter 84 is the high fdrequency components in the filtering error signal; and provide the memory of a short-term for phase-locked loop; when loop because when instantaneous noise and losing lock; can guarantee loop lock-on signal again rapidly; code tracking loop loop filter 84 uses desirable firstorder filter; structure as shown in Figure 2; by two branch roads: straight-through branch road and integration branch road are formed; straight-through branch road only contains a straight-through branch amplifier 100; the multiple that input signal is amplified appointment gets final product; the integration branch road comprises integration branch amplifier 110; integration branch road unit delay unit 112 and integration branch road totalizer 111 are formed; the signal of input can enter the integration branch road when entering straight-through branch road; input is by the amplifier amplification back of integration branch road with through the results added after the time-delay of integration branch road delay unit; result's one side after the addition is as the input of integration branch road delay unit; on the other hand as the output of integration branch road; pass through 101 additions of loop filter totalizer with the output of straight-through branch road; result after two branch road additions is as the output of code tracking loop loop filter 84; 84 outputs of code tracking loop loop filter are sent into code tracking loop totalizer 85 simultaneously with the output of code tracking loop Doppler effect correction module 87 and are carried out sum operation; result after the addition sends into code tracking loop NCO 86; adjust output PN sign indicating number clock frequency, the local PN sign indicating number clock of feasible output and the sign indicating number clock synchronization that receives.
The correlated results on schedule of correlator 81 outputs is sent into carrier tracking loop frequency discrimination/phase detector 90, carries out frequency discrimination, phase demodulation.Carrier tracking loop frequency discrimination/phase detector 90 uses cross product frequency discrimination/phase demodulation algorithms, calculates the frequency difference between input signal carrier wave and the local carrier and differs.The frequency difference that calculates and differ incoming carrier tracking loop loop filter 91, the structure of carrier tracking loop loop filter 91 and code tracking loop loop filter 84 are identical, but coefficient is different.The output of the output of carrier tracking loop loop filter 91 and carrier frequency adjustment module 73 is incoming carrier tracking loop totalizer 92 together, the tracking loop of incoming carrier as a result NCO 93 after the summation follows the tracks of the carrier wave of input signal, finishes despreading, demodulation to the input spread-spectrum signal.
3, advantage and effect: as can be seen from the above description, this intermediate frequency direct sequence spread spectrum receiver has following characteristics: adopt FFT part correlation arresting structure; Utilize ROM to store local FFT result calculated and input signal carries out related operation; One road FFT calculates the correlation of two-way simultaneously.The advantage that this structure is brought is as follows:
(1) adopt FFT part correlation arresting structure, whenever carry out the relevant peaks that a FFT computing can calculate the FFT length points, compare the method for matched filtering, significantly reduce resource consumption, the speed of catching is faster.
(2) utilize ROM to store local FFT result calculated and input signal carries out related operation, need not carry out the FFT computing, reduced local PN FFT conversion module, significantly reduced the consumption of resource local PN sign indicating number.Simultaneously can guarantee the requirement of catching under different carrier-to-noise ratios according to the length of the degree of depth of ROM and FFT.
(3) utilize one road FFT to calculate the correlation of two-way simultaneously, compare with independent F FT and calculate the two-way relevant peaks, have identical capture time, but reduced the resource consumption of one road FFT computing module, improved resource utilization.
(4) description of drawings
Fig. 1 the utility model intermediate frequency direct sequence spread spectrum receiver structural representation;
Fig. 2 the utility model loop filter structure synoptic diagram;
Symbol description is as follows among the figure:
10 front-end A/D; 20 in-phase branch multipliers;
21 in-phase branch FIR low-pass filters;
30 quadrature branch multipliers; 31 quadrature branch FIR low-pass filters;
40 integration zero clearing devices; 41 buffers; The 42FFT module;
50 branch roads, 1 complex multiplier; 51 branch roads, 1 local PN sign indicating number storage ROM;
52 branch road 1IFFT modules; 53 branch roads, 1 mean square root module;
60 branch roads, 2 complex multipliers; 61 branch roads, 2 local PN sign indicating number storage ROM;
62 branch road 2IFFT modules; 63 branch roads, 2 mean square root module;
70 automatic threshold computing modules; 71 catch judge module;
72PN sign indicating number reference position computing module; 73 carrier frequency adjustment module;
80 local PN code generators; 81 correlators; 82 code tracking loop mean square root module
83 code tracking loop phase detectors; 84 code tracking loop loop filters;
85 code tracking loop totalizers; 86 code tracking loop NCO;
87 code tracking loop Doppler effect correction modules; 90 carrier tracking loop frequency discrimination/phase detectors;
91 carrier tracking loop loop filters; 92 carrier tracking loop totalizers
93 carrier tracking loop NCO
100 straight-through branch amplifiers;
101 loop filter totalizers; 110 integration branch amplifiers;
111 integration branch road totalizers; 112 integration branch road delay units.
(5) embodiment
As shown in Figure 1, a kind of intermediate frequency direct sequence spread spectrum receiver that is used for satellite ranging of the utility model, composition comprises: the front end analog-to-digital conversion module (is called for short front-end A/D) 10, in-phase branch multiplier 20, in-phase branch finite impulse response low-pass filter (being called for short in-phase branch FIR low-pass filter) 21, quadrature branch multiplier 30, quadrature branch finite impulse response low-pass filter (being called for short quadrature branch FIR low-pass filter) 31, integration zero clearing device 40, buffer 41, fast Fourier transform module (being called for short the FFT module) 42, branch road 1 complex multiplier 50, branch road 1 local pseudo-random code storage ROM (read-only memory) (being called for short branch road 1 local PN sign indicating number storage ROM) 51, branch road 1 inverse fast fourier transform module (being called for short branch road 1IFFT module) 52, branch road 1 mean square root module 53, branch road 2 complex multipliers 60, branch road 2 local pseudo-random code storage ROM (read-only memory) (being called for short branch road 2 local PN sign indicating number storage ROM) 61, branch road 2 inverse fast fourier transform modules (being called for short branch road 2IFFT module) 62, branch road 2 mean square root module 63, automatic threshold computing module 70, catch judge module 71, pseudo-random code reference position computing module (being called for short PN sign indicating number reference position computing module) 72, carrier frequency adjustment module 73, local pseudo-random code maker (being called for short local PN code generator) 80, correlator 81, code tracking loop mean square root module 82, code tracking loop phase detector 83, code tracking loop loop filter 84, code tracking loop totalizer 85, code tracking loop digital controlled oscillator (being called for short code tracking loop NCO) 86, code tracking loop Doppler effect correction module 87, carrier tracking loop frequency discrimination/phase detector 90, carrier tracking loop loop filter 91, carrier tracking loop totalizer 92 and carrier tracking loop digital controlled oscillator (being called for short carrier tracking loop NCO) 93.
Wherein code tracking loop loop filter 84 is identical with carrier tracking loop loop filter 91 inner structures, has only the parameter difference, as shown in Figure 2: comprise straight-through branch amplifier 100, loop filter totalizer 101, integration branch amplifier 110, integration branch road totalizer 111 and integration branch road delay unit 112.
Above-mentioned all constituents, except that front-end A/D10 used the ready-made product of putting on the shelf, remainder was all realized in programmable gate array (FPGA).
Described front-end A/D10 carries out bandpass sampling with constant sampling rate to the 70M intermediate frequency spread-spectrum signal after nursing one's health, and the simulating signal of importing is transformed to digital signal, and front-end A/D10 uses ready-made product.
Carrier tracking loop NCO 93 uses direct frequency synthesizing algorithm (being called for short the DDS algorithm) to realize, be responsible for producing fixedly local carrier of the two-way identical with nominal intermediate frequency behind front-end A/D10 bandpass sampling, 90 ° of the phase phasic differences of two-way carrier wave, the signal behind the bandpass sampling of the output of carrier tracking loop NCO 93 and front-end A/D10 output is as the input of in-phase branch multiplier 20 and quadrature branch multiplier 30; In-phase branch multiplier 20 and quadrature branch multiplier 30 use the inner IP kernel of FPGA to realize, use as low-converter, input signal is down-converted to zero intermediate frequency, and in-phase branch multiplier 20 enters identical in-phase branch FIR low-pass filter 21 of structure and quadrature branch FIR low-pass filter 31 respectively with the result of quadrature branch multiplier 30 outputs; In-phase branch FIR low-pass filter 21 and quadrature branch FIR low-pass filter 31 adopt the FIR Structure Filter, use the inner IP kernel of FPGA to realize, be responsible for the signal after the down coversion is carried out filtering, frequency multiplication component and out-of-band noise after the filtering down coversion drop to the noise power in the baseband signal lower.
Result after the low-pass filtering exports to integration zero clearing device 40 on the one hand, carries out sign indicating number and catches, and exports to correlator 81 on the other hand and is used for code tracking.Integration zero clearing device 40 carries out the integration zero clearing respectively to in-phase branch FIR low-pass filter 21 and quadrature branch FIR low-pass filter 31 filtered results under the local PN sign indicating number clock control that code tracking loop NCO86 produces.The result of integration zero clearing device 40 outputs enters buffer 41, and the signal to homophase and quadrature branch carries out buffer memory simultaneously in buffer 41, and to adapt to the requirement of 42 pairs of input data rates of follow-up FFT module, buffer 41 uses the inner IP kernel of FPGA to realize.The homophase of exporting behind buffer 41 buffer memorys and the result of quadrature branch as the real part and the imaginary part input of FFT module 42, carry out Fast Fourier Transform (FFT) respectively in FFT module 42, FFT module 42 uses the inner IP kernel of FPGA to realize.The real part of FFT module 42 and imaginary part output are as one tunnel input of branch road 1 complex multiplier 50 and branch road 2 complex multipliers 60, and branch road 1 complex multiplier 50 and branch road 2 complex multipliers 60 use the inner IP kernel of FPGA to realize.Another road input of branch road 1 complex multiplier 50 and branch road 2 complex multipliers 60 is provided by branch road 1 local PN sign indicating number storage ROM51 and branch road 2 local PN sign indicating number storage ROM61 respectively.The output of FFT module 42 successively with local PN sign indicating number storage ROM51 of branch road 1 and branch road 2 local PN sign indicating numbers storage ROM61 in the local PN sign indicating number the stored result after through the FFT conversion carry out complex multiplication.Local PN sign indicating number storage ROM51 of branch road 1 and branch road 2 local PN sign indicating numbers storage ROM61 use the inner IP kernels of FPGA to realize, the data of inside solidification are the result after code phase differs half the local PN sign indicating number of the two-way sequence FFT conversion of FFT length.Result after branch road 1 complex multiplier 50 and branch road 2 complex multipliers 60 calculate sends into branch road 1IFFT module 52 respectively and branch road 2IFFT module 62 is carried out inverse fast fourier transform, and branch road 1IFFT module 52 and branch road 2IFFT module 62 use the inner IP kernel of FPGA to realize.Branch road 1 mean square root module 53 is sent in the output of branch road 1IFFT module 52 and branch road 2IFFT module 62 respectively and branch road 2 mean square root module 63 are calculated corresponding root-mean-square value, and the JPL approximate data is adopted in the calculating of root-mean-square value, promptly utilizes formula a 2 + b 2 ≈ Max ( abs ( a ) , abs ( b ) ) + 1 2 ( Min ( abs ( a ) , abs ( b ) ) ) Calculate root-mean-square value, calculate after the root-mean-square value, automatic threshold computing module 70 is sent in the output of branch road 1 mean square root module 53 and branch road 2 mean square root module 63 simultaneously, the state computation of relevant peaks goes out suitable threshold value according to this moment, output and branch road 1 mean square root module 53 with automatic threshold computing module 70, the output of branch road 2 mean square root module 63 is sent into simultaneously and is caught judge module 71, catches thresholding and branch road 1 mean square root module 53 of judge module 71 by relatively automatic threshold module 70 outputs, the correlation peak of the output of branch road 2 mean square root module 63 judges whether band spread receiver catches.Catch the output of judge module 71 and send into PN sign indicating number reference position computing module 72 and carrier frequency adjustment module 73 respectively.
When catching judge module 71 when judging that the spread spectrum answering machines are not caught, band spread receiver is operated in trapped state, and carrier frequency adjustment module 73 is being searched for all after dates of a sign indicating number, can adjust the frequency separation of output, thereby adjust the frequency of carrier tracking loop NCO 93 outputs, until finish catch till.PN sign indicating number reference position computing module 72 not outputs this moment.When catching judge module 71 when judging that the spread spectrum answering machines have been caught, will make band spread receiver change tracking mode over to by trapped state.At this moment, carrier frequency adjustment module 73 can stop to adjust, simultaneously code tracking loop Doppler effect correction module 87 is sent in output this moment, code tracking loop Doppler effect correction module 87 is according to the proportionate relationship between spreading rate and the carrier frequency, produce corresponding PN sign indicating number Doppler effect correction component, be used to adjust the sign indicating number clock of code tracking loop NCO 86 outputs, the compensating for doppler effect is to the influence of spreading code.And PN sign indicating number reference position computing module 72 will calculate the reference position of local PN by the position of when locking relevant peaks, and result calculated is sent into local PN code generator 80, adjust the phase place of the local PN sequence that local PN code generator 80 produces.
Local PN code generator 80 produce leading (PN-), on schedule (PN), (PN+) three tunnel local PN sign indicating numbers lag behind, this three tunnel local PN sign indicating number is the phase difference of half chip phase each other, and three tunnel local PN sign indicating numbers enter correlator 81 backs simultaneously and carry out related operation with the output of in-phase branch FIR low-pass filter 21 and quadrature branch FIR low-pass filter 31.Three road correlated results of output are sent into code tracking loop mean square root module 82 simultaneously and are calculated root mean square.Wherein correlated results need be sent into carrier tracking loop frequency discrimination/phase detector 90 on schedule, carries out being used for carrier track behind frequency discrimination, the phase demodulation.The result of code tracking loop mean square root module 82 outputs enters code tracking loop phase detector 83 phase demodulations, phase-demodulating principle be relatively leading relevant result and lag correlation the result, the phase place of determining code tracking loop NCO86 still lags behind in advance, and the result of phase demodulation sends into code tracking loop loop filter 84.The main effect of code tracking loop loop filter 84 is the high fdrequency components in the filtering error signal; and provide the memory of a short-term for phase-locked loop; when loop because when instantaneous noise and losing lock; can guarantee loop lock-on signal again rapidly; code tracking loop loop filter 84 uses desirable firstorder filter; structure as shown in Figure 2; by two branch roads: straight-through branch road and integration branch road are formed; straight-through branch road only contains a straight-through branch amplifier 100; the multiple that input signal is amplified appointment gets final product; the integration branch road comprises integration branch amplifier 110; integration branch road unit delay unit 112 and integration branch road totalizer 111 are formed; the signal of input can enter the integration branch road when entering straight-through branch road; input is by the amplifier amplification back of integration branch road with through the results added after the time-delay of integration branch road delay unit; result's one side after the addition is as the input of integration branch road delay unit; on the other hand as the output of integration branch road; pass through 101 additions of loop filter totalizer with the output of straight-through branch road; result after two branch road additions is as the output of code tracking loop loop filter 84; 84 outputs of code tracking loop loop filter are sent into code tracking loop totalizer 85 simultaneously with the output of code tracking loop Doppler effect correction module 87 and are carried out sum operation; result after the addition sends into code tracking loop NCO 86; adjust output PN sign indicating number clock frequency, the local PN sign indicating number clock of feasible output and the sign indicating number clock synchronization that receives.
The correlated results on schedule of correlator 81 outputs is sent into carrier tracking loop frequency discrimination/phase detector 90, carries out frequency discrimination, phase demodulation.Carrier tracking loop frequency discrimination/phase detector 90 uses cross product frequency discrimination/phase demodulation algorithms, calculates the frequency difference between input signal carrier wave and the local carrier and differs.The frequency difference that calculates and differ incoming carrier tracking loop loop filter 91, the structure of carrier tracking loop loop filter 91 and code tracking loop loop filter 84 are identical, but coefficient is different.The output of the output of carrier tracking loop loop filter 91 and carrier frequency adjustment module 73 is incoming carrier tracking loop totalizer 92 together, the tracking loop of incoming carrier as a result NCO 93 after the summation follows the tracks of the carrier wave of input signal, finishes despreading, demodulation to the input spread-spectrum signal.
The intermediate frequency direct sequence spread spectrum receiver mainly utilizes the algorithm of FFT and line correlation to finish the catching of spreading code, and principle is: the discrete Fourier transformation of simple crosscorrelation sequence z (n) is:
Z ( k ) = Σ n = 0 N - 1 Σ m = 0 N - 1 x ( m ) y ( n + m ) e - j 2 πkn / N
= Σ m = 0 N - 1 x ( m ) e j 2 πkn / N Σ n = 0 N - 1 y ( n + m ) e - j 2 πkn / N
= X * ( k ) Y ( k )
Promptly the cross correlation value of two sequences can be carries out doing after complex conjugate multiplies each other contrary FFT conversion and obtains by calculating two sequence Fourier transform results.
Also line correlation is a kind of part correlation, need that at first the signal after sampling, down coversion, the filtering is carried out the integration zero clearing and reduce sampling rate, to simplify requirement to hardware, then the data after the integration zero clearing are carried out the FFT computing, with FFT computing output result in this locality the FFT value with reference to the PN sequence carry out complex multiplication, carry out contrary FFT computing afterwards, output is the multiple correlation result who receives PN signal and local PN signal.
The multiple correlation peak is carried out square law envelope detection (asking mould), and the thresholding at the line correlation peak of going forward side by side is judged to determine whether to obtain correct catching.If all do not obtain correct catching in the cycle at a spreading code, then show because the influence of frequency deviation has exceeded capture range, need to control the centre frequency of carrier tracking loop NCO 93 to next carrier frequency by carrier frequency adjustment module 73, carry out sign indicating number in cycle at a spreading code equally and catch, till obtaining correct catching.
Catch finish after, by catching judge module 71, receiver can calculate the reference position of spreading code, by this reference position, local PN maker 80 can generate and receive the local PN sign indicating number of PN sign indicating number near-synchronous, and phase differential between the two is within one half chip period.Local PN maker 80 produces on schedule simultaneously, lead and lag three road PN sign indicating numbers.
Tracking to the PN sign indicating number is to finish by digital delay phase-locked loop leading, that hysteresis serial correlator is constituted.

Claims (1)

1. intermediate frequency direct sequence spread spectrum receiver, it is characterized in that: its composition comprises: the front end analog-to-digital conversion module is front-end A/D (10), in-phase branch multiplier (20), in-phase branch finite impulse response low-pass filter is an in-phase branch FIR low-pass filter (21), quadrature branch multiplier (30), quadrature branch finite impulse response low-pass filter is a quadrature branch FIR low-pass filter (31), integration zero clearing device (40), buffer (41), fast Fourier transform module is FFT module (42), branch road 1 complex multiplier (50), branch road 1 local pseudo-random code storage ROM (read-only memory) is branch road 1 local PN sign indicating number storage ROM (51), branch road 1 inverse fast fourier transform module is a branch road 1IFFT module (52), branch road 1 mean square root module (53), branch road 2 complex multipliers (60), branch road 2 local pseudo-random code storage ROM (read-only memory) are branch road 2 local PN sign indicating number storage ROM (61), branch road 2 inverse fast fourier transform modules are branch road 2IFFT module (62), branch road 2 mean square root module (63), automatic threshold computing module (70), catch judge module (71), pseudo-random code reference position computing module is a PN sign indicating number reference position computing module (72), carrier frequency adjustment module (73), the local PN code generator (80) of local pseudo-random code maker, correlator (81), code tracking loop mean square root module (82), code tracking loop phase detector (83), code tracking loop loop filter (84), code tracking loop totalizer (85), the code tracking loop digital controlled oscillator is code tracking loop NCO (86), code tracking loop Doppler effect correction module (87), carrier tracking loop frequency discrimination/phase detector (90), carrier tracking loop loop filter (91), carrier tracking loop totalizer (92) and carrier tracking loop digital controlled oscillator are carrier tracking loop NCO (93);
Wherein code tracking loop loop filter (84) is identical with carrier tracking loop loop filter (91) inner structure, have only the parameter difference, comprise straight-through branch amplifier (100), loop filter totalizer (101), integration branch amplifier (110), integration branch road totalizer (111) and integration branch road delay unit (112);
Above-mentioned all constituents, except that front-end A/D (10) used the ready-made product of putting on the shelf, remainder was to realize among the FPGA at programmable gate array all;
Described front-end A/D (10) carries out bandpass sampling with constant sampling rate to the 70M intermediate frequency spread-spectrum signal after nursing one's health, and the simulating signal of importing is transformed to digital signal;
Signal behind the bandpass sampling of the output of carrier tracking loop NCO 93 and front-end A/D (10) output is as the input of in-phase branch multiplier (20) and quadrature branch multiplier (30); In-phase branch multiplier (20) enters identical in-phase branch FIR low-pass filter (21) of structure and quadrature branch FIR low-pass filter (31) respectively with the result of quadrature branch multiplier (30) output;
Result after the low-pass filtering exports to integration zero clearing device (40) on the one hand, carries out sign indicating number and catches, and exports to correlator (81) on the other hand and is used for code tracking; Integration zero clearing device (40) carries out the integration zero clearing respectively to in-phase branch FIR low-pass filter (21) and the filtered result of quadrature branch FIR low-pass filter (31) under the local PN sign indicating number clock control that code tracking loop NCO (86) produces; The result of integration zero clearing device (40) output enters buffer (41); The homophase of exporting behind buffer (41) buffer memory and the result of quadrature branch import as the real part and the imaginary part of FFT module (42) respectively; The real part of FFT module (42) and imaginary part output are as one tunnel input of branch road 1 complex multiplier (50) and branch road 2 complex multipliers (60); Another road input of branch road 1 complex multiplier (50) and branch road 2 complex multipliers (60) is provided by branch road 1 local PN sign indicating number storage ROM (51) and branch road 2 local PN sign indicating number storage ROM (61) respectively; The output of FFT module (42) successively with local PN sign indicating number storage ROM (51) of branch road 1 and branch road 2 local PN sign indicating numbers storage ROM (61) in result after the local PN sign indicating number process FFT conversion of storage carry out complex multiplication; Result after branch road 1 complex multiplier (50) and branch road 2 complex multipliers (60) calculate sends into branch road 1IFFT module (52) respectively and branch road 2IFFT module (62) is carried out inverse fast fourier transform; Branch road 1 mean square root module (53) is sent in the output of branch road 1IFFT module (52) and branch road 2IFFT module (62) respectively and branch road 2 mean square root module (63) are calculated corresponding root-mean-square value, utilizes formula a 2 + b 2 ≈ Max ( abs ( a ) , abs ( b ) ) + 1 2 ( Min ( abs ( a ) , abs ( b ) ) ) Calculate after the root-mean-square value, automatic threshold computing module (70) is sent in the output of branch road 1 mean square root module (53) and branch road 2 mean square root module (63) simultaneously, the output of the output of automatic threshold computing module (70) and branch road 1 mean square root module (53), branch road 2 mean square root module (63) is sent into simultaneously caught judge module (71); Catch the output of judge module (71) and send into PN sign indicating number reference position computing module (72) and carrier frequency adjustment module (73) respectively;
When catching judge module (71) when judging that the spread spectrum answering machine is not caught, carrier frequency adjustment module (73) is being searched for all after dates of a sign indicating number, can adjust the frequency separation of output, thereby adjust the frequency of carrier tracking loop NCO (93) output, until finish catch till; PN sign indicating number reference position computing module this moment (72) is output not; When catching judge module (71) when judging that the spread spectrum answering machine has been caught, will make band spread receiver change tracking mode over to by trapped state; At this moment, carrier frequency adjustment module (73) can stop to adjust, and simultaneously code tracking loop Doppler effect correction module (87) is sent in output this moment, is used to adjust the sign indicating number clock of code tracking loop NCO (86) output; And PN sign indicating number reference position computing module (72) will calculate the reference position of local PN by the position of when locking relevant peaks, and result calculated is sent into local PN code generator (80);
Local PN code generator (80) produces in advance, on schedule, the three tunnel local PN sign indicating numbers that lag behind, three tunnel local PN sign indicating numbers enter correlator (81) back simultaneously and carry out related operation with the output of in-phase branch FIR low-pass filter (21) and quadrature branch FIR low-pass filter (31); Three road correlated results of output are sent into code tracking loop mean square root module (82) simultaneously and are calculated root mean square; Wherein correlated results need be sent into carrier tracking loop frequency discrimination/phase detector (90) on schedule; The result of code tracking loop mean square root module (82) output enters code tracking loop phase detector (83) phase demodulation, and the result of phase demodulation sends into code tracking loop loop filter (84); Code tracking loop loop filter (84) uses desirable firstorder filter, structure as shown in Figure 2, by two branch roads: straight-through branch road and integration branch road are formed, straight-through branch road only contains a straight-through branch amplifier (100), the integration branch road comprises integration branch amplifier (110), integration branch road unit delay unit (112) and integration branch road totalizer (111) are formed, the signal of input can enter the integration branch road when entering straight-through branch road, after input is amplified by integration branch amplifier (110) and through the results added after integration branch road delay unit (112) time-delay, result's one side after the addition is as the input of integration branch road delay unit (112), on the other hand as the output of integration branch road, pass through loop filter totalizer (101) addition with the output of straight-through branch road, result after two branch road additions is as the output of code tracking loop loop filter (84), code tracking loop loop filter (84) output is sent into code tracking loop totalizer (85) simultaneously with the output of code tracking loop Doppler effect correction module 87 and is carried out sum operation, result after the addition sends into code tracking loop NCO (86), adjust output PN sign indicating number clock frequency, the local PN sign indicating number clock of feasible output and the sign indicating number clock synchronization that receives;
The correlated results on schedule of correlator (81) output is sent into carrier tracking loop frequency discrimination/phase detector (90); The frequency difference that calculates and differ incoming carrier tracking loop loop filter (91), the structure of carrier tracking loop loop filter (91) and code tracking loop loop filter (84) are identical, but coefficient is different; The output of the output of carrier tracking loop loop filter (91) and carrier frequency adjustment module (73) is incoming carrier tracking loop totalizer (92) together, the tracking loop of incoming carrier as a result NCO (93) after the summation follows the tracks of the carrier wave of input signal, finishes despreading, demodulation to the input spread-spectrum signal.
CN2009202780708U 2009-12-14 2009-12-14 Intermediate-frequency direct-sequence spread spectrum receiver Expired - Fee Related CN201607527U (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2009202780708U CN201607527U (en) 2009-12-14 2009-12-14 Intermediate-frequency direct-sequence spread spectrum receiver

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2009202780708U CN201607527U (en) 2009-12-14 2009-12-14 Intermediate-frequency direct-sequence spread spectrum receiver

Publications (1)

Publication Number Publication Date
CN201607527U true CN201607527U (en) 2010-10-13

Family

ID=42952232

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2009202780708U Expired - Fee Related CN201607527U (en) 2009-12-14 2009-12-14 Intermediate-frequency direct-sequence spread spectrum receiver

Country Status (1)

Country Link
CN (1) CN201607527U (en)

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102096081A (en) * 2010-12-14 2011-06-15 东莞市泰斗微电子科技有限公司 Carrier wave and pseudo-random code stripping circuit
CN102195669A (en) * 2011-04-28 2011-09-21 马飒飒 Dual path 32-bit M-sequence digital correlator
CN102684737A (en) * 2011-03-10 2012-09-19 北京化工大学 Method for realizing multi-user spread spectrum broadcasting station based on parallel interference cancellation algorithm
CN102801434A (en) * 2012-07-20 2012-11-28 西安空间无线电技术研究所 Satellite-borne measurement and control receiver
CN103001663A (en) * 2011-09-08 2013-03-27 中国航天科工集团第三研究院第八三五七研究所 Novel wireless tracking method for spread spectrum communication (SSC) carrier waves
CN103023527A (en) * 2012-11-23 2013-04-03 陕西理工学院 Signal time-delay estimation method for direct sequence spread-spectrum ranging system and implementation thereof
CN103490746A (en) * 2012-06-13 2014-01-01 美国亚德诺半导体公司 Apparatus and methods for notch filtering
CN105301611A (en) * 2014-12-26 2016-02-03 上海华测导航技术股份有限公司 Baseband part structure for carrying out eight-frequency point processing on satellite signals
CN109407121A (en) * 2018-11-15 2019-03-01 北京遥感设备研究所 A kind of compatible acquisition and tracking device of configurable navigation signal
CN109586761A (en) * 2018-11-09 2019-04-05 北京航天万源科技有限公司 A kind of trace demodulation method of high dynamic spread spectrum signal
CN111786671A (en) * 2020-06-01 2020-10-16 南京熊猫电子股份有限公司 Method for enhancing discrete in-phase quadrature loop
CN116016072A (en) * 2023-03-22 2023-04-25 天津讯联科技有限公司 Zero intermediate frequency structure low-complexity MSK quadrature demodulation device and demodulation method thereof

Cited By (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102096081B (en) * 2010-12-14 2012-10-03 东莞市泰斗微电子科技有限公司 Carrier wave and pseudo-random code stripping circuit
CN102096081A (en) * 2010-12-14 2011-06-15 东莞市泰斗微电子科技有限公司 Carrier wave and pseudo-random code stripping circuit
CN102684737B (en) * 2011-03-10 2014-06-18 北京化工大学 Method for realizing multi-user spread spectrum broadcasting station based on parallel interference cancellation algorithm
CN102684737A (en) * 2011-03-10 2012-09-19 北京化工大学 Method for realizing multi-user spread spectrum broadcasting station based on parallel interference cancellation algorithm
CN102195669A (en) * 2011-04-28 2011-09-21 马飒飒 Dual path 32-bit M-sequence digital correlator
CN103001663A (en) * 2011-09-08 2013-03-27 中国航天科工集团第三研究院第八三五七研究所 Novel wireless tracking method for spread spectrum communication (SSC) carrier waves
CN103490746B (en) * 2012-06-13 2016-05-18 美国亚德诺半导体公司 For the apparatus and method of notch filter
CN103490746A (en) * 2012-06-13 2014-01-01 美国亚德诺半导体公司 Apparatus and methods for notch filtering
CN102801434A (en) * 2012-07-20 2012-11-28 西安空间无线电技术研究所 Satellite-borne measurement and control receiver
CN102801434B (en) * 2012-07-20 2014-08-27 西安空间无线电技术研究所 Satellite-borne measurement and control receiver
CN103023527A (en) * 2012-11-23 2013-04-03 陕西理工学院 Signal time-delay estimation method for direct sequence spread-spectrum ranging system and implementation thereof
CN105301611A (en) * 2014-12-26 2016-02-03 上海华测导航技术股份有限公司 Baseband part structure for carrying out eight-frequency point processing on satellite signals
CN109586761A (en) * 2018-11-09 2019-04-05 北京航天万源科技有限公司 A kind of trace demodulation method of high dynamic spread spectrum signal
CN109407121A (en) * 2018-11-15 2019-03-01 北京遥感设备研究所 A kind of compatible acquisition and tracking device of configurable navigation signal
CN109407121B (en) * 2018-11-15 2022-07-15 北京遥感设备研究所 Configurable navigation signal compatible capturing and tracking device
CN111786671A (en) * 2020-06-01 2020-10-16 南京熊猫电子股份有限公司 Method for enhancing discrete in-phase quadrature loop
CN111786671B (en) * 2020-06-01 2024-02-09 南京熊猫电子股份有限公司 Enhancement method of discrete type in-phase orthogonal loop
CN116016072A (en) * 2023-03-22 2023-04-25 天津讯联科技有限公司 Zero intermediate frequency structure low-complexity MSK quadrature demodulation device and demodulation method thereof

Similar Documents

Publication Publication Date Title
CN101726746B (en) Intermediate frequency direct sequence spread spectrum receiver for satellite ranging
CN201607527U (en) Intermediate-frequency direct-sequence spread spectrum receiver
CN113238261B (en) Signal capturing and tracking system of low-orbit satellite spread spectrum communication system
CN101777933B (en) Generation and capture system of encrypted frame hopping spread spectrum signal of air fleet link
CN102608626B (en) High-sensitivity satellite navigation signal capturing method and system
CN101132191B (en) Baseband signal processing method for GNSS receiver
CN101666869B (en) Method and device for secondary capturing weak satellite navigation signals
CN102594393A (en) Universal pseudocode synchronization system of comprehensive satellite baseband equipment
CN102162852B (en) Method and system for capturing weak GNSS (Global Navigation Satellite System) signal under condition of large-scale frequency deviation
CN101030787B (en) Method and apparatus for estimating signal noise ratio of frequency-amplifying signal
CN102116866B (en) Method and device for tracking global positioning system precision (GPS P) and/or Y code signal of full-cycle carrier
CN102841360B (en) Weak signal quick capture method on basis of part correlation
CN102901973B (en) Beidou satellite-based method for fast capturing signals in real time
CN102890280B (en) Compatibility capturing method of multi-mode GNSS (Global Navigation Satellite System) combination receiver
CN101261318A (en) High dynamic state spread-spectrum precise distance measurement receiving machine
CN204101732U (en) A kind of GNSS New System signal capture device
CN103278825A (en) Method for determining satellite navigation signal quality evaluation parameters
CN104360357A (en) Quick Beidou satellite signal capturing method and system based on circulation mode
CN104536016A (en) GNSS new-system signal capturing device and method
CN101458321A (en) GPS/CAPS dual mode software receiver based on PC machine
CN109412644B (en) Doppler frequency estimation method for direct sequence spread spectrum MSK signal
CN101216549B (en) Medium-high frequency wave spread-spectrum navigation system distance observed quantity extraction method
CN104570017A (en) Avoiding and stripping method for Neumann-Hoffman codes in navigation messages of Beidou navigation satellite system D1
CN111399006B (en) High-sensitivity GNSS carrier tracking loop optimization method
CN102854516A (en) Method and system for estimating carrier-to-noise ratio in GNSS (Global Navigation Satellite System) receiver

Legal Events

Date Code Title Description
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: 20101013

Termination date: 20121214