CN109738697A - A kind of frequency measurement method based on the correction of finite point discrete spectrum - Google Patents
A kind of frequency measurement method based on the correction of finite point discrete spectrum Download PDFInfo
- Publication number
- CN109738697A CN109738697A CN201910065988.2A CN201910065988A CN109738697A CN 109738697 A CN109738697 A CN 109738697A CN 201910065988 A CN201910065988 A CN 201910065988A CN 109738697 A CN109738697 A CN 109738697A
- Authority
- CN
- China
- Prior art keywords
- frequency
- correction
- signal
- error curve
- fitting
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Landscapes
- Measuring Frequencies, Analyzing Spectra (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The present invention relates to signal parameter fields of measurement, disclose a kind of frequency measurement method based on the correction of finite point discrete spectrum.Include: to inject intermediate-freuqncy signal in frequency measurement circuit, record the intermediate-freuqncy signal frequency point injected every time, the frequency for acquiring a correcting circuit every time carries position preliminary information, and is converted to the frequency values of measurement；According to the frequency values of the frequency values of measurement and injection, frequency error curve is obtained；Frequency error curve is fitted using sectional linear fitting algorithm analytical expression；According to the frequency error curve after fitting, position is carried to frequency and carries out secondary correction.Technical solution of the present invention is on the basis of a Spectrum Correction, the error curve based on fitting, carries position to frequency and carries out secondary correction, has further reduced the error range of spectrum position estimation；Technical solution of the present invention substantially reduces the time of measuring to signal carrier frequency information using parallel, streamlined Processing Algorithm, has the processing capability in real time to multi signal.
Description
Technical field
The present invention relates to signal parameter fields of measurement, especially a kind of frequency measurement based on the correction of finite point discrete spectrum
Method.
Background technique
When directlying adopt finite point discrete spectrum analysis method and measuring to the carrier frequency value of actual signal, work as sampling
It is maximum between direct measurement frequency and actual signal carrier frequency when frequency and measured signal carrier frequency be not at integral multiple relationship
It can produce a frequency resolution (sample frequency is divided by limited points) error.It is answered to frequency measurement accuracy is more demanding
With under environment, need to directly adopt finite point discrete spectrum analysis result be corrected, further decrease worst error range
To meet High Precision Frequency requirement.Currently, specifically including that interpolation to the bearing calibration of finite point discrete spectrum measurement result
Method, frequency spectrum refinement method, power enhanced.The width that interpolation method passes through time big spectral line and maximum spectral line in search discrete spectrum result
Ratio is spent, and then calculates the deviation of maximum spectral line Yu actual signal frequency, using the deviation to the frequency at maximum spectral line
It is corrected with phase, when deviation is smaller, interpolation anisotropy is caused to cause large error.Frequency spectrum refinement method passes through shifting
Frequently, part increase differentiates the modes such as points and improves frequency measurement accuracy to a certain degree, but processing complexity is high, equivalent analysis point
Number is more, real-time is poor.Power enhanced carries out discrete spectrum using the characteristic that the energy barycenter that various windows are composed approaches coordinate origin
Analysis, frequency estimation accuracy is related with signals selected modulated window and signal-to-noise ratio, poor universality.
Summary of the invention
The technical problems to be solved by the present invention are: said frequencies measurement there are aiming at the problem that, provide one kind and be based on
The frequency measurement method of finite point discrete spectrum correction.
The technical solution adopted by the invention is as follows: a kind of frequency measurement method based on the correction of finite point discrete spectrum, tool
Body includes following procedure:
Step 1, intermediate-freuqncy signal is injected in frequency measurement circuit, records the intermediate-freuqncy signal frequency point injected every time, and acquisition is every
The frequency of a secondary correcting circuit carries position preliminary information, and is converted to the frequency values of measurement；
Step 2, according to the frequency values of the frequency values of measurement and injection, frequency error curve is obtained；
Step 3, frequency error curve is fitted using sectional linear fitting algorithm analytical expression；
Step 4, according to the frequency error curve after fitting, position is carried to frequency and carries out secondary correction.
Further, the detailed process of the step 1 are as follows: step 11, intermediate-freuqncy signal is digital intermediate frequency letter by AD conversion
Number, then parallel FFT processor is sent into after carrying out time-domain windowed modulation；Step 12, fft processor is according to the flowing water signal segment time
Update discrete complex signal spectrum output；Step 13, each frequency point is sought according to the complex signal of output spectrum using amplitude fitting circuit
The approximate range, is judged the rough location of current demand signal by the approximate range value for distinguishing channel after Threshold detection；Step 14, it obtains
Rough location guidance peak search circuit be focused the preliminary information to determine local maxima position of spectral line.
Further, in the step 3, the frequency error curve after fitting is in working frequency range, the signal-to-noise ratio item of restriction
Meet the frequency accuracy requirement of setting under part.
It further, further include the verification process of frequency error curve after fitting: by analytical expression in the step 3
Error curve after fitting is added in the realization logic an of correcting circuit；Again by different capacity, different frequent points, different letters
It makes an uproar the intermediate-freuqncy signal injected frequency measuring circuit of ratio, repeats step 1~step 3, the error curve of adjustment analytical expression fitting,
Until the frequency measurement circuit in working frequency range, limit power bracket in, limit signal-to-noise ratio under the conditions of meet frequency measurement essence
Degree requires.
Further, the detailed process of the step 4 are as follows: the peak search circuit exports multiple position of spectral line simultaneously
Preliminary information, the preliminary information of the position of spectral line are transferred to position correction unit, and position correction unit is according to the spectral line of output
The error curve of position and the fitting being obtained ahead of time is corrected.
Further, the preliminary information of the position of spectral line is also transmitted to position synchronization unit, and position synchronization unit is completed
The timing synchronization of data processing, is then output to position correction unit.
Compared with prior art, by adopting the above technical scheme have the beneficial effect that technical solution of the present invention in primary frequency
On the basis of spectrum correction, the error curve based on fitting carries position to frequency and carries out secondary correction, further reduced spectrum position
The error range of estimation；Technical solution of the present invention is substantially reduced using parallel, streamlined Processing Algorithm believes signal carrier frequency
The time of measuring of breath, improve serialization, intermittent discrete Fourier transform technology bring discrete spectrum time of withdrawing information it is long,
The problem of real-time difference, has the processing capability in real time to multi signal；Technical solution of the present invention can be used for communication, radar,
In the electronic equipments such as electronic reconnaissance, navigation, the quick and precisely measurement to various signal carrier frequency is realized.
Detailed description of the invention
Fig. 1 is the flow diagram of the frequency measurement method based on the correction of finite point discrete spectrum.
Specific embodiment
The present invention is described further with reference to the accompanying drawing.
As shown in Figure 1, the frequency measurement method of the present embodiment includes correction course twice:
Primary correction: injecting intermediate-freuqncy signal in frequency measurement circuit, records the intermediate-freuqncy signal frequency point injected every time, acquisition
The frequency of a correcting circuit carries position preliminary information every time, and is converted to the frequency values of measurement；Its intermediate frequency carries position preliminary information
It generates after peak value searching process；After primary correction, band segment still remains biggish position deviation.
Secondary correction: according to the frequency values of the frequency values of measurement and injection, frequency error curve is obtained；Using piecewise linearity
Fitting algorithm is fitted frequency error curve with analytical expression；According to the frequency error curve after fitting, position is carried to frequency and is carried out
Secondary correction.The process of above-mentioned secondary correction is distributed quasi- according to the measurement error under the conditions of time domain, frequency domain and different signal-to-noise ratio
The secondary correction error curve for adapting to majority of case is closed out, the error range of spectrum position estimation has been further reduced.
Embodiment 1: the detailed process once corrected: the detailed process of the step 1 are as follows: intermediate-freuqncy signal is by AD conversion
Digital medium-frequency signal, then parallel FFT processor is sent into after carrying out time-domain windowed modulation；Fft processor is according to flowing water signal segment
Time updates discrete complex signal spectrum output；Each frequency discrimination is sought according to the complex signal of output spectrum using amplitude fitting circuit to believe
The approximate range, is judged the rough location of current demand signal by the approximate range value in road after Threshold detection；The rough location of acquisition
Guidance peak search circuit is focused the preliminary information to determine local maxima position of spectral line.
Embodiment 2: the frequency error curve after fitting is in working frequency range, satisfaction setting under the conditions of the signal-to-noise ratio of restriction
Frequency accuracy requirement.
After once fitting, the frequency accuracy of error curve is not necessarily met the requirements, so also needing to frequency after fitting
Rate error curve is verified: the error curve after analytical expression is fitted is added to the realization logic an of correcting circuit
In；Again by different capacity, different frequent points (the intermediate-freuqncy signal frequency point re-injected will be in operating frequency range), different noises
The intermediate-freuqncy signal injected frequency measuring circuit of ratio repeats step 1~step 3, adjusts the error curve of analytical expression fitting, directly
Meet frequency-measurement accuracy in working frequency range, in the power bracket of restriction, under the conditions of the signal-to-noise ratio of restriction to the frequency measurement circuit
It is required that.
Embodiment 3: the error curve obtained based on embodiment 2 carries out secondary correction on the basis of primary correction.It is described
Peak search circuit exports the preliminary information of multiple position of spectral line simultaneously, and the preliminary information of the position of spectral line is transferred to position school
Positive unit, position correction unit are corrected according to the position of spectral line of output and the error curve for the fitting being obtained ahead of time.
The preliminary information of the position of spectral line is also transmitted to position synchronization unit, and position synchronization unit completes data processing
Timing synchronization is then output to position correction unit.Position synchronization unit ensure that the timing synchronization of data processing.
The invention is not limited to specific embodiments above-mentioned.The present invention, which expands to, any in the present specification to be disclosed
New feature or any new combination, and disclose any new method or process the step of or any new combination.If this
Field technical staff is altered or modified not departing from the unsubstantiality that spirit of the invention is done, should belong to power of the present invention
The claimed range of benefit.
Claims (6)
1. a kind of frequency measurement method based on the correction of finite point discrete spectrum, which is characterized in that specifically include following procedure:
Step 1, it injects intermediate-freuqncy signal in frequency measurement circuit, records the intermediate-freuqncy signal frequency point injected every time, acquisition every time one
The frequency of secondary correcting circuit carries position preliminary information, and is converted to the frequency values of measurement；
Step 2, according to the frequency values of the frequency values of measurement and injection, frequency error curve is obtained；
Step 3, frequency error curve is fitted using sectional linear fitting algorithm analytical expression；
Step 4, according to the frequency error curve after fitting, position is carried to frequency and carries out secondary correction.
2. the frequency measurement method as described in claim 1 based on the correction of finite point discrete spectrum, which is characterized in that the step
Rapid 1 detailed process are as follows: step 11, intermediate-freuqncy signal is digital medium-frequency signal by AD conversion, then after carrying out time-domain windowed modulation
It is sent into parallel FFT processor；Step 12, fft processor updates discrete complex signal spectrum output according to the flowing water signal segment time；Step
Rapid 13, the approximate range value of each frequency discrimination channel is sought according to the complex signal of output spectrum using amplitude fitting circuit, by this
Approximate range judges the rough location of current demand signal after Threshold detection；Step 14, the rough location of acquisition guides peak value searching
Circuit is focused the preliminary information to determine local maxima position of spectral line.
3. the frequency measurement method as claimed in claim 2 based on the correction of finite point discrete spectrum, which is characterized in that the step
In rapid 3, the frequency accuracy that the frequency error curve after fitting meets setting in working frequency range, under the conditions of the signal-to-noise ratio of restriction is wanted
It asks.
4. the frequency measurement method as claimed in claim 3 based on the correction of finite point discrete spectrum, which is characterized in that the step
In rapid 3, further include the verification process of frequency error curve after fitting: the error curve after analytical expression is fitted is added to one
In the realization logic of secondary correcting circuit；The intermediate-freuqncy signal injected frequency of different capacity, different frequent points, different signal-to-noise ratio is surveyed again
Circuit is measured, step 1~step 3, the error curve of adjustment analytical expression fitting, until the frequency measurement circuit is working are repeated
Meet frequency-measurement accuracy requirement in frequency range, in the power bracket of restriction, under the conditions of the signal-to-noise ratio of restriction.
5. the frequency measurement method as claimed in claim 4 based on the correction of finite point discrete spectrum, which is characterized in that the step
Rapid 4 detailed process are as follows: the peak search circuit exports the preliminary information of multiple position of spectral line simultaneously, the position of spectral line
Preliminary information is transferred to position correction unit, and position correction unit is according to the position of spectral line of output and the fitting being obtained ahead of time
Error curve is corrected.
6. the frequency measurement method as claimed in claim 5 based on the correction of finite point discrete spectrum, which is characterized in that the spectrum
The preliminary information of line position is also transmitted to position synchronization unit, and position synchronization unit completes the timing synchronization of data processing, then
It exports to position correction unit.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910065988.2A CN109738697B (en) | 2019-01-24 | 2019-01-24 | Frequency measurement method based on finite point discrete spectrum correction |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910065988.2A CN109738697B (en) | 2019-01-24 | 2019-01-24 | Frequency measurement method based on finite point discrete spectrum correction |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109738697A true CN109738697A (en) | 2019-05-10 |
| CN109738697B CN109738697B (en) | 2021-02-26 |
Family
ID=66365837
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910065988.2A Active CN109738697B (en) | 2019-01-24 | 2019-01-24 | Frequency measurement method based on finite point discrete spectrum correction |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109738697B (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112162152A (en) * | 2020-08-31 | 2021-01-01 | 南京亿杰明信息技术有限公司 | Sine wave coherent pulse train signal frequency estimation method based on phase straight line fitting |
| CN112198365A (en) * | 2020-09-29 | 2021-01-08 | 中国船舶重工集团公司第七二四研究所 | Method for realizing real-time frequency spectrum monitoring equipment |
| CN113075450A (en) * | 2021-02-22 | 2021-07-06 | 中国电子科技集团公司第二十九研究所 | Method for analyzing radio frequency and intermediate frequency of broadband frequency compressed signal based on sampling rate |
| CN113884761A (en) * | 2021-09-13 | 2022-01-04 | 中国电子科技集团公司第二十九研究所 | Method, device, equipment and medium for improving frequency measurement precision |
| CN116879627A (en) * | 2023-09-04 | 2023-10-13 | 中国电子科技集团公司第二十九研究所 | Nanosecond non-coherent narrow pulse sequence frequency measurement system |
Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101334431A (en) * | 2007-12-20 | 2008-12-31 | 复旦大学 | Electric network harmonic frequency spectrum interpolation correction analytical method |
| CN101738601A (en) * | 2009-11-27 | 2010-06-16 | 西安电子科技大学 | System and method for measuring speed of locomotive based on radar near field echo power spectrum characteristics |
| JP2011153984A (en) * | 2010-01-28 | 2011-08-11 | Anritsu Corp | Apparatus and method for observing waveform |
| CN102185627A (en) * | 2011-05-06 | 2011-09-14 | 京信通信系统（中国）有限公司 | Frequency calibrating method and device under TD-SCDMA mode |
| CN103454497A (en) * | 2013-09-10 | 2013-12-18 | 南京理工大学 | Phase difference measuring method based on improved windowing discrete Fourier transform |
| CN103454495A (en) * | 2013-09-13 | 2013-12-18 | 电子科技大学 | Self-adaptive high-precision fast spectral analysis method |
| CN103944606A (en) * | 2014-02-28 | 2014-07-23 | 电子科技大学 | Self-adaptation frequency hopping pattern generation method |
| WO2014203708A1 (en) * | 2013-06-17 | 2014-12-24 | アルプス電気株式会社 | Signal frequency calculation method |
| CN105704081A (en) * | 2016-01-12 | 2016-06-22 | 中国电子科技集团公司第十研究所 | Optimization and demodulation method of multiple offset carrier wave amplitude modulation signals |
| CN105842536A (en) * | 2016-03-30 | 2016-08-10 | 国家电网公司 | Improved ZoomFFT and ApFFT algorithm integrated method for analyzing dense harmonics and inter-harmonics |
| CN106053937A (en) * | 2016-06-21 | 2016-10-26 | 电子科技大学 | Fundamental wave frequency measurement method based on FFT (Fast Fourier Transform) + FT (Fourier Transform) |
-
2019
- 2019-01-24 CN CN201910065988.2A patent/CN109738697B/en active Active
Patent Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101334431A (en) * | 2007-12-20 | 2008-12-31 | 复旦大学 | Electric network harmonic frequency spectrum interpolation correction analytical method |
| CN101738601A (en) * | 2009-11-27 | 2010-06-16 | 西安电子科技大学 | System and method for measuring speed of locomotive based on radar near field echo power spectrum characteristics |
| JP2011153984A (en) * | 2010-01-28 | 2011-08-11 | Anritsu Corp | Apparatus and method for observing waveform |
| CN102185627A (en) * | 2011-05-06 | 2011-09-14 | 京信通信系统（中国）有限公司 | Frequency calibrating method and device under TD-SCDMA mode |
| WO2014203708A1 (en) * | 2013-06-17 | 2014-12-24 | アルプス電気株式会社 | Signal frequency calculation method |
| CN103454497A (en) * | 2013-09-10 | 2013-12-18 | 南京理工大学 | Phase difference measuring method based on improved windowing discrete Fourier transform |
| CN103454495A (en) * | 2013-09-13 | 2013-12-18 | 电子科技大学 | Self-adaptive high-precision fast spectral analysis method |
| CN103944606A (en) * | 2014-02-28 | 2014-07-23 | 电子科技大学 | Self-adaptation frequency hopping pattern generation method |
| CN105704081A (en) * | 2016-01-12 | 2016-06-22 | 中国电子科技集团公司第十研究所 | Optimization and demodulation method of multiple offset carrier wave amplitude modulation signals |
| CN105842536A (en) * | 2016-03-30 | 2016-08-10 | 国家电网公司 | Improved ZoomFFT and ApFFT algorithm integrated method for analyzing dense harmonics and inter-harmonics |
| CN106053937A (en) * | 2016-06-21 | 2016-10-26 | 电子科技大学 | Fundamental wave frequency measurement method based on FFT (Fast Fourier Transform) + FT (Fourier Transform) |
Non-Patent Citations (7)
| Title |
|---|
| HONG CHEN 等: "Automatic correcting method of discrete spectrum with rectangular window", 《2011 4TH INTERNATIONAL CONGRESS ON IMAGE AND SIGNAL PROCESSING》 * |
| QIN SHUREN 等: "CORRECTION METHOD OF DISCRETE SPECTRUM BY EXTRACTING FREQUENCY", 《2007 IEEE INTERNATIONAL CONFERENCE ON SIGNAL PROCESSING AND COMMUNICATIONS》 * |
| ZHANG JIE: "Research of High Precision Frequency Measure Algorithm Based on LabVIEW", 《THE EIGHTH INTERNATIONAL CONFERENCE ON ELECTRONIC MEASUREMENT AND INSTRUMENTS》 * |
| 丁康 等: "FFT+FT离散频谱校正法参数估计精度", 《机械工程学报》 * |
| 朱小勇 等: "离散频谱校正方法的综合比较", 《信号处理》 * |
| 谢胜 等: "基于FFT并二次修正的Rife频率估计算法", 《探测与控制学报》 * |
| 郭振聚 等: "基于离散频谱校正提高数字化交流测量精度", 《计量与测试技术》 * |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112162152A (en) * | 2020-08-31 | 2021-01-01 | 南京亿杰明信息技术有限公司 | Sine wave coherent pulse train signal frequency estimation method based on phase straight line fitting |
| CN112162152B (en) * | 2020-08-31 | 2024-01-26 | 南京亿杰明信息技术有限公司 | Sine wave coherent pulse train signal frequency estimation method based on phase straight line fitting |
| CN112198365A (en) * | 2020-09-29 | 2021-01-08 | 中国船舶重工集团公司第七二四研究所 | Method for realizing real-time frequency spectrum monitoring equipment |
| CN113075450A (en) * | 2021-02-22 | 2021-07-06 | 中国电子科技集团公司第二十九研究所 | Method for analyzing radio frequency and intermediate frequency of broadband frequency compressed signal based on sampling rate |
| CN113884761A (en) * | 2021-09-13 | 2022-01-04 | 中国电子科技集团公司第二十九研究所 | Method, device, equipment and medium for improving frequency measurement precision |
| CN113884761B (en) * | 2021-09-13 | 2023-04-25 | 中国电子科技集团公司第二十九研究所 | Method, device, equipment and medium for improving frequency measurement precision |
| CN116879627A (en) * | 2023-09-04 | 2023-10-13 | 中国电子科技集团公司第二十九研究所 | Nanosecond non-coherent narrow pulse sequence frequency measurement system |
| CN116879627B (en) * | 2023-09-04 | 2023-11-21 | 中国电子科技集团公司第二十九研究所 | Nanosecond non-coherent narrow pulse sequence frequency measurement system |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109738697B (en) | 2021-02-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109738697A (en) | A kind of frequency measurement method based on the correction of finite point discrete spectrum | |
| CN109617631B (en) | Adaptive receiving method of reconnaissance system based on digital channelized instantaneous parameter measurement | |
| CN110518936B (en) | High dynamic Doppler fast capturing method and communication system for hypersonic aircraft | |
| CN110068727B (en) | Single-frequency signal frequency estimation method based on Candan-Rife comprehensive interpolation | |
| CN110161478B (en) | Waveform design method based on clutter power spectral density self-optimization | |
| CN105652085B (en) | A kind of fundamental wave and harmonic parameters method of estimation based on the incomplete S-transformation of improvement | |
| US9063227B2 (en) | Positioning apparatus and signal processing method thereof | |
| CN105223481A (en) | Based on the partial-discharge ultrahigh-frequency signal initial time defining method of difference energy function | |
| CN116892898B (en) | Track error detection method, device and system for agricultural machinery | |
| CN104330782A (en) | Time domain and modulation domain parameter combined measuring method of triangular frequency-modulation pulse signals | |
| Wei et al. | Doppler estimation based on dual-HFM signal and speed spectrum scanning | |
| CN104158605A (en) | PWM-code-based method capable of realizing self-adaption simultaneous wireless information and power transfer (SWIPT) frequency | |
| CN102495383B (en) | Automatic phase correction method | |
| EP3373540A1 (en) | Frequency offset processing method and system | |
| Yang et al. | Polynomial fitting and interpolation method in tdoa estimation of sensors network | |
| CN110048741A (en) | A kind of method for parameter estimation of the Frequency Hopping Signal based on Short-Time Fractional Fourier Transform | |
| CN105738698B (en) | A kind of harmonic parameters algorithm for estimating based on centre frequency shift | |
| CN109521269B (en) | Amplitude modulation signal digital frequency measurement method | |
| CN104133110B (en) | Method for achieving cursor for frequency hopping signal time-frequency characteristic measurement | |
| CN110673223A (en) | SIP data acquisition and processing method without synchronous current acquisition and transmission and frequency domain measurement | |
| CN111611686A (en) | Detection method for communication signal time-frequency domain | |
| CN108020818A (en) | A kind of sine pulse signal distance measuring method under noise based on sliding DFT | |
| CN105223552A (en) | Based on the disturbance identification method of radio spectrum line identification | |
| CN112883787B (en) | Short sample low-frequency sinusoidal signal parameter estimation method based on spectrum matching | |
| CN108241144B (en) | FMCW radar waveform modulation method and device |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |