KR101633187B1 - 레이더의 도플러 주파수 추정 장치 및 방법 - Google Patents
레이더의 도플러 주파수 추정 장치 및 방법 Download PDFInfo
- Publication number
- KR101633187B1 KR101633187B1 KR1020140161161A KR20140161161A KR101633187B1 KR 101633187 B1 KR101633187 B1 KR 101633187B1 KR 1020140161161 A KR1020140161161 A KR 1020140161161A KR 20140161161 A KR20140161161 A KR 20140161161A KR 101633187 B1 KR101633187 B1 KR 101633187B1
- Authority
- KR
- South Korea
- Prior art keywords
- doppler frequency
- estimating
- fractional part
- update time
- fractional
- Prior art date
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/50—Systems of measurement based on relative movement of target
- G01S13/52—Discriminating between fixed and moving objects or between objects moving at different speeds
- G01S13/522—Discriminating between fixed and moving objects or between objects moving at different speeds using transmissions of interrupted pulse modulated waves
- G01S13/524—Discriminating between fixed and moving objects or between objects moving at different speeds using transmissions of interrupted pulse modulated waves based upon the phase or frequency shift resulting from movement of objects, with reference to the transmitted signals, e.g. coherent MTi
- G01S13/53—Discriminating between fixed and moving objects or between objects moving at different speeds using transmissions of interrupted pulse modulated waves based upon the phase or frequency shift resulting from movement of objects, with reference to the transmitted signals, e.g. coherent MTi performing filtering on a single spectral line and associated with one or more range gates with a phase detector or a frequency mixer to extract the Doppler information, e.g. pulse Doppler radar
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
- G01R23/02—Arrangements for measuring frequency, e.g. pulse repetition rate; Arrangements for measuring period of current or voltage
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/06—Systems determining position data of a target
- G01S13/08—Systems for measuring distance only
- G01S13/32—Systems for measuring distance only using transmission of continuous waves, whether amplitude-, frequency-, or phase-modulated, or unmodulated
- G01S13/34—Systems for measuring distance only using transmission of continuous waves, whether amplitude-, frequency-, or phase-modulated, or unmodulated using transmission of continuous, frequency-modulated waves while heterodyning the received signal, or a signal derived therefrom, with a locally-generated signal related to the contemporaneously transmitted signal
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/50—Systems of measurement based on relative movement of target
- G01S13/505—Systems of measurement based on relative movement of target using Doppler effect for determining closest range to a target or corresponding time, e.g. miss-distance indicator
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/50—Systems of measurement based on relative movement of target
- G01S13/52—Discriminating between fixed and moving objects or between objects moving at different speeds
- G01S13/536—Discriminating between fixed and moving objects or between objects moving at different speeds using transmission of continuous unmodulated waves, amplitude-, frequency-, or phase-modulated waves
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
레이더의 도플러 주파수 추정 장치 및 방법이 개시된다.
레이더의 도플러 주파수 추정 방법은 오브젝트에 의해 반사된 레이더 신호의 업데이트 시간 별 도플러 주파수를 추정하는 단계; 및 상기 레이더 신호의 업데이트 시간별 수신 전력을 기초로 상기 업데이트 시간별 도플러 주파수에 가중평균(weighted average)을 적용하여 상기 레이더 신호의 최종 도플러 주파수를 추정하는 단계를 포함할 수 있다.
레이더의 도플러 주파수 추정 방법은 오브젝트에 의해 반사된 레이더 신호의 업데이트 시간 별 도플러 주파수를 추정하는 단계; 및 상기 레이더 신호의 업데이트 시간별 수신 전력을 기초로 상기 업데이트 시간별 도플러 주파수에 가중평균(weighted average)을 적용하여 상기 레이더 신호의 최종 도플러 주파수를 추정하는 단계를 포함할 수 있다.
Description
본 발명은 레이더의 도플러 주파수 추정 장치 및 방법에 관한 것으로, 보다 상세하게는 도플러 주파수를 이용하는 레이더에서 도플러 주파수의 추정 성능을 높이는 장치 및 방법에 관한 것이다.
FMCW (Frequency Modulation Continuous Wave) 레이더는 선형적 주파수 변조 레이더 신호를 송신하고, 송신한 레이더 신호가 탐지하고자 하는 오브젝트에 반사된 신호를 수신하며, 수신한 신호와 송신한 신호 간의 도플러 주파수를 추정하여 오브젝트의 거리 및 속도를 탐지할 수 있다.
종래의 FMCW 레이더는 도플러 주파수의 추정 성능을 높이기 위하여 MUSIC이나 ESPRIT와 같은 고 정밀 알고리즘을 사용하고 있다. 종래의 FMCW 레이더가 사용하는 고정밀 알고리즘은 고유값 분해(Eigen value decomposition) 함수를 이용해서 정확한 도플러 주파수를 추정할 수 있으나, 하드웨어 비용이 비싸다는 문제점이 있다.
그러나, 고정밀 알고리즘을 사용하기 위한 하드웨어의 고비용에 의하여 차량에 탑재기 위한 FMCW 레이더와 같은 저가형 FMCW 레이더를 구현하기 힘든 실정이다.
따라서, 고비용의 하드웨어 없이도 레이더에서 도플러 주파수의 추정 성능을 높일 수 있는 방법이 요청되고 있다.
본 발명은 업데이트 시간별로 도플러 주파수를 추정하고, 업데이트 시간별로 추정한 도플러 주파수에 가중 평균을 적용함으로써, 도플러 주파수의 추정 성능을 향상시키는 장치 및 방법을 제공할 수 있다.
본 발명의 일실시예에 따른 레이더의 도플러 주파수 추정 방법은 오브젝트에 의해 반사된 레이더 신호의 업데이트 시간 별 도플러 주파수를 추정하는 단계; 및 상기 레이더 신호의 업데이트 시간별 수신 전력을 기초로 상기 업데이트 시간별 도플러 주파수에 가중평균(weighted average)을 적용하여 상기 레이더 신호의 최종 도플러 주파수를 추정하는 단계를 포함할 수 있다.
본 발명의 일실시예에 따른 레이더의 도플러 주파수 추정 방법의 시간별 도플러 주파수를 추정하는 단계는, 이산 푸리에 변환한 상기 레이더 신호를 기초로 상기 시간별 도플러 주파수의 정수부(integer part)를 추정하는 단계; 상기 정수부를 이용하여 상기 시간별 도플러 주파수의 소수부(fractional part)를 추정하는 단계; 및 상기 정수부, 상기 소수부 및 샘플링 주파수에 기초하여 상기 레이더 신호의 업데이트 시간 동안 정규화된 도플러 주파수를 추정하는 단계를 포함할 수 있다.
본 발명의 일실시예에 따른 레이더의 도플러 주파수 추정 방법의 소수부를 추정하는 단계는, 상기 정수부를 이용하여 상기 시간별 도플러 주파수의 첫 번째 소수부를 추정하는 단계; 및 첫 번째 소수부, 또는 이전 단계에서 추정한 소수부와 상기 정수부를 이용하여 상기 시간별 도플러 주파수의 두 번째 이후의 나머지 소수부를 추정하는 단계를 포함하고, 상기 나머지 소수부를 추정하는 단계는, 상기 시간별 도플러 주파수의 소수부를 모두 추정할 때까지 반복하여 수행될 수 있다.
본 발명의 일실시예에 따른 레이더의 도플러 주파수 추정 방법의 첫 번째 소수부를 추정하는 단계; 상기 정수부를 이용하여 제1 이산 푸리에 변환 계수를 정의하는 단계; 및 상기 제1 이산 푸리에 변환 계수 및 레이더 신호를 이용하여 상기 제1 소수부를 추정하는 단계를 포함할 수 있다.
본 발명의 일실시예에 따른 레이더의 도플러 주파수 추정 방법의 최종 도플러 주파수를 추정하는 단계는, 상기 레이더 신호의 수신 전력, 및 상기 시간별 도플러 주파수의 소수부에 가중평균을 적용하여 최종 도플러 주파수의 소수부를 추정하는 단계; 및 상기 최종 도플러 주파수의 소수부, 상기 최종 도플러 주파수의 정수부, 및 샘플링 주파수에 기초하여 최종 도플러 주파수를 추정하는 단계를 포함할 수 있다.
본 발명의 일실시예에 따른 레이더의 도플러 주파수 추정 방법의 최종 도플러 주파수의 소수부를 추정하는 단계는, 상기 도플러 주파수의 소수부 및 상기 수신 전력과 상기 업데이트 시간 동안 누적된 수신 전력 간의 비율을 추정하는 단계; 및 상기 비율의 누적에 따라 최종 도플러 주파수의 소수부를 추정하는 단계를 포함할 수 있다.
본 발명의 일실시예에 따른 레이더의 도플러 주파수 추정 장치는 오브젝트에 의해 반사된 레이더 신호의 업데이트 시간 별 도플러 주파수를 추정하는 업데이트 시간 별 도플러 주파수를 추정부; 및 상기 레이더 신호의 업데이트 시간별 수신 전력을 기초로 상기 업데이트 시간별 도플러 주파수에 가중평균을 적용하여 상기 레이더 신호의 최종 도플러 주파수를 추정하는 최종 도플러 주파수 추정부를 포함할 수 있다.
본 발명의 일실시예에 따른 레이더의 도플러 주파수 추정 장치의 업데이트 시간 별 도플러 주파수를 추정부는, 이산 푸리에 변환한 상기 레이더 신호를 기초로 상기 시간별 도플러 주파수의 정수부를 추정하는 정수부 추정부; 상기 정수부를 이용하여 상기 시간별 도플러 주파수의 소수부를 추정하는 소수부 추정부; 및 상기 정수부, 상기 소수부 및 샘플링 주파수에 기초하여 상기 레이더 신호의 업데이트 시간 동안 정규화된 도플러 주파수를 추정하는 정규화된 도플러 주파수 추정부를 포함할 수 있다.
본 발명의 일실시예에 따른 레이더의 도플러 주파수 추정 장치의 소수부 추정부는, 상기 정수부를 이용하여 상기 시간별 도플러 주파수의 첫 번째 소수부를 추정하고, 첫 번째 소수부, 또는 이전 단계에서 추정한 소수부와 상기 정수부를 이용하여 상기 시간별 도플러 주파수의 두 번째 이후의 나머지 소수부를 추정하며, 상기 나머지 소수부를 추정하는 과정은 상기 시간별 도플러 주파수의 소수부를 모두 추정할 때까지 반복하여 수행될 수 있다.
본 발명의 일실시예에 따른 레이더의 도플러 주파수 추정 장치의 소수부 추정부는, 상기 정수부를 이용하여 제1 이산 푸리에 변환 계수를 정의하고, 상기 제1 이산 푸리에 변환 계수 및 레이더 신호를 이용하여 상기 제1 소수부를 추정할 수 있다.
본 발명의 일실시예에 따른 레이더의 도플러 주파수 추정 장치의 최종 도플러 주파수 추정부는, 상기 레이더 신호의 수신 전력, 및 상기 시간별 도플러 주파수의 소수부에 가중평균을 적용하여 최종 도플러 주파수의 소수부를 추정하는 소수부 추정부; 및 상기 최종 도플러 주파수의 소수부, 상기 최종 도플러 주파수의 정수부, 및 샘플링 주파수에 기초하여 최종 도플러 주파수를 추정하는 도플러 주파수 추정부를 포함할 수 있다.
본 발명의 일실시예에 따른 레이더의 도플러 주파수 추정 장치의 소수부 추정부는, 상기 도플러 주파수의 소수부 및 상기 수신 전력과 상기 업데이트 시간 동안 누적된 수신 전력 간의 비율을 추정하고, 상기 비율의 누적에 따라 최종 도플러 주파수의 소수부를 추정할 수 있다.
본 발명의 일실시예에 의하면, 업데이트 시간별로 도플러 주파수를 추정하고, 업데이트 시간별로 추정한 도플러 주파수에 가중 평균을 적용함으로써, 도플러 주파수의 추정 성능을 향상시킬 수 있다.
도 1은 본 발명의 일실시예에 따른 레이더 시스템을 나타내는 도면이다.
도 2는 본 발명의 일실시예에 따른 도플러 주파수 추정 장치의 구성을 나타내는 도면이다.
도 3은 본 발명의 일실시예에 따른 업데이트 시간별 도플러 주파수 추정부의 구성을 나타내는 도면이다.
도 4는 본 발명의 일실시예에 따른 최종 도플러 주파수 추정부의 구성을 나타내는 도면이다.
도 5는 본 발명의 일실시예에 따라 도플러 주파수를 추정하는 레이더 시스템의 동작 방법을 도시한 플로우차트이다.
도 6은 본 발명의 일실시예에 따른 업데이트 시간별 도플러 주파수 추정 과정을 도시한 플로우차트이다.
도 7은 본 발명의 일실시예에 따른 최종 주파수 추정 과정을 도시한 플로우차트이다.
도 2는 본 발명의 일실시예에 따른 도플러 주파수 추정 장치의 구성을 나타내는 도면이다.
도 3은 본 발명의 일실시예에 따른 업데이트 시간별 도플러 주파수 추정부의 구성을 나타내는 도면이다.
도 4는 본 발명의 일실시예에 따른 최종 도플러 주파수 추정부의 구성을 나타내는 도면이다.
도 5는 본 발명의 일실시예에 따라 도플러 주파수를 추정하는 레이더 시스템의 동작 방법을 도시한 플로우차트이다.
도 6은 본 발명의 일실시예에 따른 업데이트 시간별 도플러 주파수 추정 과정을 도시한 플로우차트이다.
도 7은 본 발명의 일실시예에 따른 최종 주파수 추정 과정을 도시한 플로우차트이다.
이하, 본 발명의 실시예를 첨부된 도면을 참조하여 상세하게 설명한다. 본 발명의 일실시예에 따른 도플러 주파수 추정 방법은 레이더 시스템에 포함된 도플러 주파수 추정 장치에 의해 수행될 수 있다.
도 1은 본 발명의 일실시예에 따른 레이더 시스템을 나타내는 도면이다.
도 1을 참고하면, 본 발명의 일실시예에 따른 레이더 시스템은 송신 장치(110), 수신 장치(120), 도플러 주파수 추정 장치(130), 및 오브젝트 탐지 장치(140)를 포함할 수 있다.
송신 장치(110)는 송신 안테나를 통해 레이더 신호(101)를 송신할 수 있다. 구체적으로, 송신 장치(110)는 주파수 변조 연속 파(Frequency Modulation Continuous Wave) 방식으로 생성된 레이더 신호(101)를 송신할 수 있다.
수신 장치(120)는 수신 안테나를 통해 레이더 신호(101)가 오브젝트(102)에 의해 반사된 레이더 신호(103)를 수신할 수 있다. 이때, 수신 장치(120)는 기 설정된 업데이트 시간이 되면, 이전 업데이트 시간부터 현재 업데이트 시간까지 수신한 레이더 신호(103)를 도플러 주파수 추정 장치(130)에 업데이트할 수 있다.
예를 들어, 수신 장치(120)가 수신한 전체 레이더 신호(103)인 s(t)는 수학식 1로 나타낼 수 있다.
이때, AT는 송신장치(110)가 송신한 레이더 신호(101)의 진폭이고, M은 수신 장치(120)가 레이더 신호(103)를 도플러 주파수 추정 장치(130)에 업데이트한 회수일 수 있다. 또한, N은 수신 장치(120)가 이전 업데이트 시간부터 현재 업데이트 시간까지 수신한 레이더 신호(103)에 포함된 전송 주기(PRI: Pulse Repetition Interval)의 개수일 수 있다. 그리고, tp는 시간 정규화 인자이고, TUP는 레이더 신호(103)의 N 전송 주기가 업데이트되는 시간일 수 있다. 또한, TPRI는 1 / fPRI일 수 있다. 이때, fPRI는 전송 주기의 주파수일 수 있다.
또한, j번째 업데이트 시간 동안 수신 장치(120)가 수신한 레이더 신호 rj(t)는 도플러 시프트(Doppler shift)에 따라 수학식 2와 같이 정리될 수 있다.
이때, sj(t)는 N개의 전송 주기를 포함하는 레이더 신호(103)이고, aj는 복소 진폭이며, fd,j는도플러 주파수일 수 있다. 또한, 는 지연 시간이고, 는 j 번째 업데이트 시간 동안 이동한 오브젝트에 의하여 영향을 받은 초기 위상일 수 있다. 그리고, W(t)는 분산 o2을 가진 복소 백색 부가 가우시안 잡음(AWGN: additive white Gaussian noise)일 수 있다.
또한, 레이더 신호(103)가 fs=1/Ts의 비율로 샘플링된 경우, 수신 장치(120)가 수신한 레이더 신호(103)의 N개의 샘플 r(nTs), n=0,...,(N-1)들은 도플러 주파수의 사인 성분으로 변환될 수 있다. 이때, Ts= TPRI일 수 있다. 또한, 샘플 rj는 수학식 3으로 나타낼 수 있다.
이때, fdn,j는 j 번째 업데이트 동안 정규화된 도플러 주파수이고, W[n]은 t=nTs인 경우, ω(t)의 샘플링 된 값일 수 있다. 또한, 기저 대역 신호는 S(t)의 피크 포인트에서 샘플링 되는 것으로 가정될 수 있다. 그리고, j번째 업데이트 동안 fdn,j = fd,j / fs = v/N인 경우, 정규화 도플러 주파수 fdn,j는 도플러 주파수 fdn로 정의될 수 있다. 이때, N은 샘플의 개수이고, v는 실수일 수 있다. 또한, 실수 v는 정수부인 kpj과 j 번째 업데이트 시간의 소수부인 δ로구성될 수 있다.
도플러 주파수 추정 장치(130)는 수신 장치(120)가 수신한 레이더 신호(103)의 도플러 주파수를 추정할 수 있다. 이때, 도플러 주파수 추정 장치(130)는 업데이트 시간 별로 레이더 신호(103)의 도플러 주파수를 추정하고, 업데이트 시간별 도플러 주파수에 가중 평균을 적용하여 레이더 신호의 최종 도플러 주파수를 추정할 수 있다.
도플러 주파수 추정 장치(130)의 구체적인 구성 및 동작은 이하 도 2를 참조하여 상세히 설명한다.
오브젝트 탐지 장치(140)는 도플러 주파수 추정 장치(130)가 추정한 레이더 신호(103)의 최종 도플러 주파수를 이용하여 오브젝트(102)의 방향 및 속도를 탐지할 수 있다.
레이더 시스템(100)은 수신한 레이더 신호(103)의 신호 대 잡음비(SNR: signal to noise ratio)에 따라 업데이트 시간별 도플러 주파수에 가중 평균을 적용함으로써, 도플러 주파수의 추정 성능을 향상시킬 수 있다.
도 2는 본 발명의 일실시예에 따른 도플러 주파수 추정 장치의 구성을 나타내는 도면이다.
도 2를 참고하면, 본 발명의 일실시예에 따른 도플러 주파수 추정 장치(130)는 업데이트 시간별 도플러 주파수 추정부(210), 및 최종 도플러 주파수 추정부(220)를 포함할 수 있다.
업데이트 시간별 도플러 주파수 추정부(210)는 오브젝트에 의해 반사된 레이더 신호의 업데이트 시간 별 도플러 주파수를 추정할 수 있다. 이때, 업데이트 시간별 도플러 주파수 추정부(210)는 반복 주파수 추정(IFE: Iterative Frequency Estimation)를 수행하여 업데이트 시간 별 도플러 주파수를 추정할 수 있다. 업데이트 시간별 도플러 주파수 추정부(210)의 구체적인 구성 및 동작은 도 3을 참조로 상세히 설명한다.
최종 도플러 주파수 추정부(220)는 레이더 신호의 업데이트 시간별 수신 전력을 기초로 업데이트 시간별 도플러 주파수 추정부(210)가 추정한 업데이트 시간별 도플러 주파수에 가중평균(weighted average)을 적용하여 레이더 신호의 최종 도플러 주파수를 추정할 수 있다. 레이더 신호의 신호 대 잡음비가 낮은 경우, 레이더 신호가 확산되므로, 도플러 주파수를 추정하기 어려울 수 있다. 따라서, 최종 도플러 주파수 추정부(220)는 시간 다이버시티(time diversity)를 고려하여 레이더 신호의 최종 도플러 주파수를 추정할 수 있다.
최종 도플러 주파수 추정부(220)의 구체적인 구성 및 동작은 도 4를 참조로 상세히 설명한다.
도 3은 본 발명의 일실시예에 따른 업데이트 시간별 도플러 주파수 추정부의 구성을 나타내는 도면이다.
도 3을 참고하면, 본 발명의 일실시예에 따른 업데이트 시간별 도플러 주파수 추정부(210)는 정수부 추정부(310), 소수부 추정부(320), 및 정규화된 도플러 주파수 추정부(330)를 포함할 수 있다.
정수부 추정부(310)는 수신 장치(120)가 수신한 레이더 신호를 기초로 업데이트 시간별 도플러 주파수의 정수부(integer part)를 추정할 수 있다.
이산 푸리에 변환(DFT: discrete Fourier transform)의 최대 진폭 계수에 대응하는 주파수는 최대우도 추정법(MLE: Maximum Likelihood Estimation)에서 추정한 주파수의 근사치로 선택될 수 있다. 그리고, 레이더 신호의 주파수는 업데이트 시간에 포함된 전송 주기와 샘플링 주파수 및 업데이트 시간별 도플러 주파수의 정수부로 추정될 수 있다. 예를 들어, 정수부 추정부(310)는 수학식 4를 이용하여 레이더 신호의 주파수 를 추정할 수 있다.
이때, N은 업데이트 시간 j에 포함되는 전송 주기이고, fs는 레이더 신호를 1/Ts비율로 샘플링한 샘플링 주파수일 수 있다. 그리고, kmax,j는 업데이트 시간별 도플러 주파수의 정수부일 수 있다. 또한, 정수부 추정부(310)는 수학식 5를 이용하여 업데이트 시간별 도플러 주파수의 정수부를 계산할 수 있다.
이때, yj[k]는 수신 장치(120)가 수신한 레이더 신호 rj[t]를 이산 푸리에 변환한 값일 수 있다. 예를 들어, yj[k]와 rj[t]를 수학식 6과 같이 정의할 수 있다.
이때, yj는 DFT (rj)일 수 있다. 또한, DFT (N)은 DFT 연산자의 N 포인트 복소수일 수 있다.
즉, 정수부 추정부(310)는 수신 장치(120)가 수신한 레이더 신호를 이산 푸리에 변환하여 yj[k]를 추정하고, yj[k]를 이용하여 업데이트 시간별 도플러 주파수의 정수부인 kmax,j를 추정할 수 있다.
소수부 추정부(320)는 정수부 추정부(310)가 추정한 업데이트 시간별 도플러 주파수의 정수부를 이용하여 업데이트 시간별 도플러 주파수의 소수부(fractional part)를 추정할 수 있다.
이때, 소수부 추정부(320)는 업데이트 시간별 도플러 주파수의 정수부를 이용하여 업데이트 시간별 도플러 주파수의 첫 번째 소수부를 추정할 수 있다. 그리고, 소수부 추정부(320)는 첫 번째 소수부, 또는 이전 단계에서 추정한 소수부와 정수부를 이용하여 업데이트 시간별 도플러 주파수의 두 번째 이후의 나머지 소수부를 순차적으로 추정할 수 있다.
먼저, 소수부 추정부(320)는 첫 번째 반복 계산법에 따라 레이더 신호의 주기도표(period gram)의 값을 이용하여 j번째 업데이트 시간별 도플러 주파수의 첫 번째 소수부 δj,1를 추정할 수 있다.
구체적으로, 소수부 추정부(320)는 업데이트 시간별 도플러 주파수의 정수부를 이용하여 제1 이산 푸리에 변환 계수를 정의하고, 제1 이산 푸리에 변환 계수 및 레이더 신호를 이용하여 업데이트 시간별 도플러 주파수의 첫 번째 소수부를 추정할 수 있다.
예를 들어, 소수부 추정부(320)는 수학식 7을 이용하여 첫 번째 소수부를 추정하기 위하여 수정한 제1 수정 DFT 계수 Xp를 정의할 수 있다.
그리고, 소수부 추정부(320)는 수학식 7과 같이 제1 수정 DFT 계수 Xp 에서 수신 장치(120)가 수신한 레이더 신호를 결합할 수 있다. 이때, 수신 장치(120)가 수신한 레이더 신호는 정현파 성분을 포함할 수 있다. 또한, 소수부 추정부(320)는 무소음 상태의 케이스에 필요한 작업을 수행할 수 있다. 예를 들어, 소수부 추정부(320)는 수학식 8을 이용하여 무소음 상태의 케이스에 필요한 작업의 수율(yield)을 추정할 수 있다.
다음으로, 소수부 추정부(320)는 첫 번째 소수부, 또는 이전 단계에서 추정한 소수부와 정수부를 이용하여 업데이트 시간별 도플러 주파수의 두 번째 이후의 나머지 소수부를 추정할 수 있다. 그리고, 소수부 추정부(320)는 시간별 도플러 주파수의 소수부를 모두 추정할 때까지 상기 과정을 반복하여 수행할 수 있다.
예를 들어, 소수부 추정부(320)는 업데이트 시간별 도플러 주파수의 첫 번째 소수부와 정수부 추정부(310)가 추정한 정수부에 두 번째 반복 계산법(second iterative scheme)을 이용하여 시간별 도플러 주파수의 두 번째 소수부를 추정할 수 있다. 이때, 소수부 추정부(320)는 q를 2로 설정하고, 수학식 9를 이용하여 나머지 소수부를 추정하기 위하여 수정한 제2 수정 DFT 계수 Xp를 정의할 수 있다.
그리고, 소수부 추정부(320)는 의 q를 3 내지 j번째 업데이트 시간별 도플러 주파수의 마지막 소수부의 값까지 순차적으로 증가하며, 이전 단계에서 추정한 소수부와 정수부를 이용하여 업데이트 시간별 도플러 주파수의 현재 단계 소수부를 추정할 수 있다.
정규화된 도플러 주파수 추정부(330)는 정수부 추정부(310)가 추정한 정수부와 소수부 추정부(320)가 추정한 소수부 및 샘플링 주파수에 기초하여 레이더 신호의 업데이트 시간 동안 정규화된 도플러 주파수를 추정할 수 있다. 예를 들어, 정규화된 도플러 주파수 추정부(330)는 수학식 11을 이용하여 j 번째 업데이트 시간 동안 정규화된 도플러 주파수를 추정할 수 있다.
이때, 는 j 번째 업데이트 동안의 도플러 주파수의 소수부일 수 있다. 또한, 이고, q = 2일 수 있다. 즉, 정규화된 도플러 주파수 추정부(330)는 수학식 9에서 추정한 를 이용하여 도플러 주파수의 소수부 를 획득할 수 있다.
도 4는 본 발명의 일실시예에 따른 최종 도플러 주파수 추정부의 구성을 나타내는 도면이다.
도 4를 참고하면, 본 발명의 일실시예에 따른 최종 도플러 주파수 추정부(220)는 수신 전력 추정부(410), 소수부 추정부(420) 및 도플러 주파수 추정부(430)를 포함할 수 있다.
수신 전력 추정부(410)는 업데이트 시간 별로 수신 장치(120)가 수신한 레이더 신호의 수신 전력을 추정할 수 있다. 예를 들어, 수신 전력 추정부(410)는 수학식 12를 사용하여 j 번째 업데이트 시간 동안 수신된 레이더 신호의 수신 전력 을 추정할 수 있다. 이때, 수신 전력은 수신한 레이더 신호의 복소 진폭의 파워일 수 있다.
소수부 추정부(420)는 수신 전력 추정부(410)가 추정한 레이더 신호의 수신 전력, 및 소수부(320)가 추정한 시간별 도플러 주파수의 소수부에 가중평균을 적용하여 최종 도플러 주파수의 소수부를 추정할 수 있다.
구체적으로, 소수부 추정부(420)는 도플러 주파수의 소수부 및 상기 수신 전력과 상기 업데이트 시간 동안 누적된 수신 전력 간의 비율을 추정하고, 상기 비율의 누적에 따라 최종 도플러 주파수의 소수부를 추정할 수 있다. 예를 들어, 소수부 추정부(420)는 수학식 13을 이용하여 최종 도플러 주파수의 소수부인 를 추정할 수 있다.
도플러 주파수 추정부(430)는 소수부 추정부(420)가 추정한 최종 도플러 주파수의 소수부, 최종 도플러 주파수의 정수부, 및 샘플링 주파수에 기초하여 최종 도플러 주파수를 추정할 수 있다. 예를 들어, 도플러 주파수 추정부(430)는 수학식 14를 이용하여 최종 도플러 주파수인 를 추정할 수 있다.
이때, kmax는 최종 도플러 주파수의 정수부이며, kmax=kmax,j일 수 있다.
도 5는 본 발명의 일실시예에 따라 도플러 주파수를 추정하는 레이더 시스템의 동작 방법을 도시한 플로우차트이다.
단계(510)에서 송신 장치(110)는 송신 안테나를 통해 레이더 신호(101)를 송신할 수 있다.
단계(520)에서 수신 장치(120)는 수신 안테나를 통해 단계(510)에서 송신한 레이더 신호(101)가 오브젝트(102)에 의해 반사된 레이더 신호(103)를 수신할 수 있다. 이때, 수신 장치(120)는 기 설정된 업데이트 시간이 되면, 이전 업데이트 시간부터 현재 업데이트 시간까지 수신한 레이더 신호(103)를 도플러 주파수 추정 장치(130)에 업데이트할 수 있다.
단계(530)에서 업데이트 시간별 도플러 주파수 추정부(210)는 단계(520)에서 수신한 레이더 신호의 업데이트 시간 별로 도플러 주파수를 추정할 수 있다. 이때, 업데이트 시간별 도플러 주파수 추정부(210)는 반복 주파수 추정(IFE)를 수행하여 업데이트 시간 별 도플러 주파수를 추정할 수 있다.
단계(540)에서 최종 도플러 주파수 추정부(220)는 단계(530)에서 추정한 레이더 신호의 업데이트 시간별 수신 전력을 기초로 업데이트 시간별 도플러 주파수 추정부(210)가 추정한 업데이트 시간별 도플러 주파수에 가중평균을 적용하여 레이더 신호의 최종 도플러 주파수를 추정할 수 있다.
단계(550)에서 오브젝트 탐지 장치(140)는 단계(540)에서 추정한 레이더 신호(103)의 최종 도플러 주파수를 이용하여 오브젝트(102)의 방향 및 속도를 탐지할 수 있다.
도 6은 본 발명의 일실시예에 따른 업데이트 시간별 도플러 주파수 추정 과정을 도시한 플로우차트이다. 이때, 도 6에 도시된 단계(610) 내지 단계(660)은 도 5의 단계(530)에 포함될 수 있다.
단계(610)에서 정수부 추정부(310)는 단계(520)에서 수신한 레이더 신호를 기초로 업데이트 시간별 도플러 주파수의 정수부를 추정할 수 있다.
이때, 정수부 추정부(310)는 수신 장치(120)가 수신한 레이더 신호를 이산 푸리에 변환하여 yj[k]를 추정하고, yj[k]를 이용하여 업데이트 시간별 도플러 주파수의 정수부인 kmax,j를 추정할 수 있다.
단계(620)에서 소수부 추정부(320)는 단계(610)에서 추정한 업데이트 시간별 도플러 주파수의 정수부를 이용하여 업데이트 시간별 도플러 주파수의 첫 번째 소수부를 추정할 수 있다.
단계(630)에서 소수부 추정부(320)는 단계(620)에서 추정한 첫 번째 소수부, 단계(610)에서 추정한 정수부를 이용하여 업데이트 시간별 도플러 주파수의 두 번째 소수부를 추정할 수 있다.
단계(630)에서 소수부 추정부(320)는 업데이트 시간별 도플러 주파수의 모든 소수부를 추정하였는지 여부를 확인할 수 있다. 업데이트 시간별 도플러 주파수의 소수부 중에서 추정하지 않은 소수부가 있는 경우, 소수부 추정부(320)는 단계(650)를 수행할 수 있다. 또한, 단계(630), 또는 단계(650)에서 업데이트 시간별 도플러 주파수의 모든 소수부가 추정된 경우, 소수부 추정부(320)는 단계(660)를 수행할 수 있다.
단계(650)에서 소수부 추정부(320)는 단계(630), 또는 이전에 수행된 단계(650)에서 추정한 소수부와 정수부를 이용하여 업데이트 시간별 도플러 주파수의 현재 단계의 소수부를 추정할 수 있다.
단계(660)에서 정규화된 도플러 주파수 추정부(330)는 단계(610)에서 추정한 정수부와 단계(630), 또는 단계(650)에서 추정한 소수부 및 샘플링 주파수에 기초하여 레이더 신호의 업데이트 시간 동안 정규화된 도플러 주파수를 추정할 수 있다.
도 7은 본 발명의 일실시예에 따른 최종 주파수 추정 과정을 도시한 플로우차트이다. 이때, 도 7에 도시된 단계(710) 내지 단계(730)은 도 5의 단계(540)에 포함될 수 있다.
단계(710)에서 수신 전력 추정부(410)는 업데이트 시간 별로 수신 장치(120)가 수신한 레이더 신호의 수신 전력을 추정할 수 있다.
단계(720)에서 소수부 추정부(420)는 단계(720)에서 추정한 레이더 신호의 수신 전력, 및 단계(530)에서 추정한 시간별 도플러 주파수의 소수부에 가중평균을 적용하여 최종 도플러 주파수의 소수부를 추정할 수 있다. 구체적으로, 소수부 추정부(420)는 도플러 주파수의 소수부 및 상기 수신 전력과 상기 업데이트 시간 동안 누적된 수신 전력 간의 비율을 추정하고, 상기 비율의 누적에 따라 최종 도플러 주파수의 소수부를 추정할 수 있다.
단계(730)에서 도플러 주파수 추정부(430)는 단계(720)에서 추정한 최종 도플러 주파수의 소수부, 최종 도플러 주파수의 정수부, 및 샘플링 주파수에 기초하여 최종 도플러 주파수를 추정할 수 있다.
본 발명은 업데이트 시간별로 도플러 주파수를 추정하고, 업데이트 시간별로 추정한 도플러 주파수에 가중 평균을 적용함으로써, 도플러 주파수의 추정 성능을 향상시킬 수 있다.
이상과 같이 본 발명은 비록 한정된 실시예와 도면에 의해 설명되었으나, 본 발명은 상기의 실시예에 한정되는 것은 아니며, 본 발명이 속하는 분야에서 통상의 지식을 가진 자라면 이러한 기재로부터 다양한 수정 및 변형이 가능하다.
그러므로, 본 발명의 범위는 설명된 실시예에 국한되어 정해져서는 아니 되며, 후술하는 특허청구범위뿐 아니라 이 특허청구범위와 균등한 것들에 의해 정해져야 한다.
110: 송신 장치
120: 수신 장치
130: 도플러 주파수 추정 장치
310: 정수부 추정부
320: 소수부 추정부
120: 수신 장치
130: 도플러 주파수 추정 장치
310: 정수부 추정부
320: 소수부 추정부
Claims (12)
- 오브젝트에 의해 반사된 레이더 신호의 업데이트 시간 별 도플러 주파수를 추정하는 단계; 및
상기 레이더 신호의 업데이트 시간별 수신 전력을 기초로 상기 업데이트 시간별 도플러 주파수에 가중평균(weighted average)을 적용하여 상기 레이더 신호의 최종 도플러 주파수를 추정하는 단계
를 포함하는 레이더의 도플러 주파수 추정 방법. - 제1항에 있어서,
상기 시간별 도플러 주파수를 추정하는 단계는,
이산 푸리에 변환한 상기 레이더 신호를 기초로 상기 시간별 도플러 주파수의 정수부(integer part)를 추정하는 단계;
상기 정수부를 이용하여 상기 시간별 도플러 주파수의 소수부(fractional part)를 추정하는 단계; 및
상기 정수부, 상기 소수부 및 샘플링 주파수에 기초하여 상기 레이더 신호의 업데이트 시간 동안 정규화된 도플러 주파수를 추정하는 단계
를 포함하는 레이더의 도플러 주파수 추정 방법. - 제2항에 있어서,
상기 소수부를 추정하는 단계는,
상기 정수부를 이용하여 상기 시간별 도플러 주파수의 첫 번째 소수부를 추정하는 단계; 및
첫 번째 소수부, 또는 이전 단계에서 추정한 소수부와 상기 정수부를 이용하여 상기 시간별 도플러 주파수의 두 번째 이후의 나머지 소수부를 추정하는 단계
를 포함하고,
상기 나머지 소수부를 추정하는 단계는,
상기 시간별 도플러 주파수의 소수부를 모두 추정할 때까지 반복하여 수행되는 레이더의 도플러 주파수 추정 방법. - 제3항에 있어서,
상기 첫 번째 소수부를 추정하는 단계;
상기 정수부를 이용하여 제1 이산 푸리에 변환 계수를 정의하는 단계; 및
상기 제1 이산 푸리에 변환 계수 및 레이더 신호를 이용하여 상기 제1 소수부를 추정하는 단계
를 포함하는 레이더의 도플러 주파수 추정 방법. - 제1항에 있어서,
상기 최종 도플러 주파수를 추정하는 단계는,
상기 도플러 주파수의 소수부와 상기 수신 전력의 곱을 계산하는 단계;
상기 곱과 상기 업데이트 시간 동안 누적된 수신 전력 간의 비율을 추정하는 단계;
상기 비율의 누적에 따라 최종 도플러 주파수의 소수부를 추정하는 단계; 및
상기 최종 도플러 주파수의 소수부, 상기 최종 도플러 주파수의 정수부, 및 샘플링 주파수에 기초하여 최종 도플러 주파수를 추정하는 단계
를 포함하는 레이더의 도플러 주파수 추정 방법. - 삭제
- 오브젝트에 의해 반사된 레이더 신호의 업데이트 시간 별 도플러 주파수를 추정하는 업데이트 시간 별 도플러 주파수를 추정부; 및
상기 레이더 신호의 업데이트 시간별 수신 전력을 기초로 상기 업데이트 시간별 도플러 주파수에 가중평균(weighted average)을 적용하여 상기 레이더 신호의 최종 도플러 주파수를 추정하는 최종 도플러 주파수 추정부
를 포함하는 레이더의 도플러 주파수 추정 장치. - 제7항에 있어서,
상기 업데이트 시간 별 도플러 주파수를 추정부는,
이산 푸리에 변환한 상기 레이더 신호를 기초로 상기 시간별 도플러 주파수의 정수부(integer part)를 추정하는 정수부 추정부;
상기 정수부를 이용하여 상기 시간별 도플러 주파수의 소수부(fractional part)를 추정하는 소수부 추정부; 및
상기 정수부, 상기 소수부 및 샘플링 주파수에 기초하여 상기 레이더 신호의 업데이트 시간 동안 정규화된 도플러 주파수를 추정하는 정규화된 도플러 주파수 추정부
를 포함하는 레이더의 도플러 주파수 추정 장치. - 제8항에 있어서,
상기 소수부 추정부는,
상기 정수부를 이용하여 상기 시간별 도플러 주파수의 첫 번째 소수부를 추정하고, 첫 번째 소수부, 또는 이전 단계에서 추정한 소수부와 상기 정수부를 이용하여 상기 시간별 도플러 주파수의 두 번째 이후의 나머지 소수부를 추정하며, 상기 나머지 소수부를 추정하는 과정은 상기 시간별 도플러 주파수의 소수부를 모두 추정할 때까지 반복하여 수행되는 레이더의 도플러 주파수 추정 장치. - 제9항에 있어서,
상기 소수부 추정부는,
상기 정수부를 이용하여 제1 이산 푸리에 변환 계수를 정의하고, 상기 제1 이산 푸리에 변환 계수 및 레이더 신호를 이용하여 상기 제1 소수부를 추정하는 레이더의 도플러 주파수 추정 장치. - 제7항에 있어서,
상기 최종 도플러 주파수 추정부는,
상기 도플러 주파수의 소수부와 상기 수신 전력의 곱을 계산하고, 상기 곱과 상기 업데이트 시간 동안 누적된 수신 전력 간의 비율을 추정하며, 상기 비율의 누적에 따라 최종 도플러 주파수의 소수부를 추정하는 소수부 추정부; 및
상기 최종 도플러 주파수의 소수부, 상기 최종 도플러 주파수의 정수부, 및 샘플링 주파수에 기초하여 최종 도플러 주파수를 추정하는 도플러 주파수 추정부
를 포함하는 레이더의 도플러 주파수 추정 장치.
- 삭제
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
KR1020140161161A KR101633187B1 (ko) | 2014-11-18 | 2014-11-18 | 레이더의 도플러 주파수 추정 장치 및 방법 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
KR1020140161161A KR101633187B1 (ko) | 2014-11-18 | 2014-11-18 | 레이더의 도플러 주파수 추정 장치 및 방법 |
Publications (2)
Publication Number | Publication Date |
---|---|
KR20160059359A KR20160059359A (ko) | 2016-05-26 |
KR101633187B1 true KR101633187B1 (ko) | 2016-06-23 |
Family
ID=56104741
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
KR1020140161161A KR101633187B1 (ko) | 2014-11-18 | 2014-11-18 | 레이더의 도플러 주파수 추정 장치 및 방법 |
Country Status (1)
Country | Link |
---|---|
KR (1) | KR101633187B1 (ko) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2019528458A (ja) * | 2016-06-21 | 2019-10-10 | シーアイエス4ユー | ドップラーレーダーを利用した人感センサー |
KR101967812B1 (ko) * | 2017-02-17 | 2019-04-11 | 옴니센서(주) | 가중치 적용 칼만 필터를 이용한 레이더 신호 처리 장치 및 그를 이용한 표적 검출 시스템 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2011069779A (ja) * | 2009-09-28 | 2011-04-07 | Mitsubishi Electric Corp | レーダ装置 |
JP2011191195A (ja) * | 2010-03-15 | 2011-09-29 | Honda Elesys Co Ltd | レーダ装置及びコンピュータプログラム |
JP2012018058A (ja) * | 2010-07-07 | 2012-01-26 | Honda Elesys Co Ltd | レーダ装置及びコンピュータプログラム |
KR101449869B1 (ko) * | 2013-10-10 | 2014-10-13 | 재단법인대구경북과학기술원 | 레이더 신호의 주파수 추정 장치 및 그 방법 |
-
2014
- 2014-11-18 KR KR1020140161161A patent/KR101633187B1/ko active IP Right Grant
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2011069779A (ja) * | 2009-09-28 | 2011-04-07 | Mitsubishi Electric Corp | レーダ装置 |
JP2011191195A (ja) * | 2010-03-15 | 2011-09-29 | Honda Elesys Co Ltd | レーダ装置及びコンピュータプログラム |
JP2012018058A (ja) * | 2010-07-07 | 2012-01-26 | Honda Elesys Co Ltd | レーダ装置及びコンピュータプログラム |
KR101449869B1 (ko) * | 2013-10-10 | 2014-10-13 | 재단법인대구경북과학기술원 | 레이더 신호의 주파수 추정 장치 및 그 방법 |
Also Published As
Publication number | Publication date |
---|---|
KR20160059359A (ko) | 2016-05-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10126417B2 (en) | Radar apparatus and object detecting method | |
KR20190096291A (ko) | 위상을 보정하는 레이더 감지 | |
Bechter et al. | Automotive radar interference mitigation by reconstruction and cancellation of interference component | |
US8619908B2 (en) | Wireless ranging system and related methods | |
KR101674747B1 (ko) | I/q 불균형을 보상하기 위한 레이더 신호 처리 장치 및 방법 | |
CN103777178A (zh) | 一种同步误差补偿方法、设备及系统 | |
US11125857B2 (en) | Moving object detection system and moving object detection method | |
US20180149693A1 (en) | Reflectometry method and device for diagnosing cables in use | |
KR101633187B1 (ko) | 레이더의 도플러 주파수 추정 장치 및 방법 | |
KR101900008B1 (ko) | 이동 표적 탐지 장치 및 방법 | |
EP3335058B1 (en) | Processing received radiation reflected from a target | |
KR102625211B1 (ko) | 무선 통신 신호의 수신 시각 산출 장치 및 방법 | |
KR101449869B1 (ko) | 레이더 신호의 주파수 추정 장치 및 그 방법 | |
CN104202103A (zh) | 一种欠采样的信号到达时差测定方法 | |
CN104320360B (zh) | 一种基于分数阶傅里叶变换的线性调频信号时延估计方法 | |
KR101584449B1 (ko) | 레이더 신호 처리 방법 및 장치 | |
US20070237213A1 (en) | Method for estimating one or more parameters of a ultra wideband signal and a receiver system for receiving ultra wideband signals | |
JP2010060318A (ja) | レーダ装置 | |
KR101446439B1 (ko) | Fmcw 레이더의 고정밀 주파수 추정 방법 및 고정밀 주파수 추정 장치 | |
Häfner et al. | High-resolution parameter estimation for chirp-sequence radar considering hardware impairments | |
KR100916640B1 (ko) | 무선 통신 기반 송,수신 장치 간 거리 추정 방법 | |
JP2012149951A (ja) | レーダ装置 | |
US8532203B1 (en) | Communication systems and methods using blind channel estimation based on cumulants | |
RU2824754C1 (ru) | Способ частотно-временной обработки сигналов от нескольких движущихся объектов | |
CN104217119A (zh) | 一种基于S-Method时频分解的机动目标自适应补偿检测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A201 | Request for examination | ||
E902 | Notification of reason for refusal | ||
E701 | Decision to grant or registration of patent right | ||
GRNT | Written decision to grant | ||
FPAY | Annual fee payment |
Payment date: 20190325 Year of fee payment: 4 |