TWI697688B - Frequency modulated continuous wave processing device - Google Patents
Frequency modulated continuous wave processing device Download PDFInfo
- Publication number
- TWI697688B TWI697688B TW108130214A TW108130214A TWI697688B TW I697688 B TWI697688 B TW I697688B TW 108130214 A TW108130214 A TW 108130214A TW 108130214 A TW108130214 A TW 108130214A TW I697688 B TWI697688 B TW I697688B
- Authority
- TW
- Taiwan
- Prior art keywords
- continuous wave
- frequency domain
- signals
- frequency
- azimuth
- Prior art date
Links
- 230000001427 coherent effect Effects 0.000 claims abstract description 16
- 230000035559 beat frequency Effects 0.000 claims description 42
- 238000000034 method Methods 0.000 claims description 8
- 238000010586 diagram Methods 0.000 description 5
- 230000005540 biological transmission Effects 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 238000001514 detection method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000012615 high-resolution technique Methods 0.000 description 1
- 235000015096 spirit Nutrition 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
Description
本發明係關於一種處理裝置,且特別關於一種調頻連續波(FMCW)處理裝置。 The present invention relates to a processing device, and more particularly to a frequency modulated continuous wave (FMCW) processing device.
調頻連續波(FMCW)雷達是一種很成熟的技術。各家企業發展雷達技術主要著重於雷達資料的處理,主要包括雷達目標探測、目標參數估測、目標識別,以及更後端的雷達資料運用,其中參數包含距離、速度與方位。在目標參數獲取上,方位角一般而言是比較困難的,熟知的高解析度技術有多訊號分類(Multiple Signal Classification,MUSIC)演算法、基於旋轉不變技術的信號參數估計(estimation of signal parameter via rotational invariant techniques,ESPRIT)演算法、最小絕對值收斂和選擇算子(LASSO)演算法以及各種演算法。 Frequency modulated continuous wave (FMCW) radar is a very mature technology. The development of radar technology by various companies mainly focuses on the processing of radar data, which mainly includes radar target detection, target parameter estimation, target recognition, and more back-end radar data application. The parameters include distance, speed and azimuth. In terms of target parameter acquisition, the azimuth angle is generally more difficult. The well-known high-resolution techniques include the Multiple Signal Classification (MUSIC) algorithm and the estimation of signal parameter based on the rotation invariant technology. Via rotational invariant techniques, ESPRIT) algorithm, minimum absolute value convergence and selection operator (LASSO) algorithm, and various algorithms.
調頻連續波雷達的一個重要應用是車用雷達,獲取的資料特徵是複雜、訊噪比(SNR)不佳、以及變動快速。上述技術在這樣的環境難以達到高效能、低運算複雜度與高即時性。此外,在多用戶使用環境下,雷達之間的互相干擾將會是個嚴重問題,習知抗干擾技術包括錯開使用頻段、錯開發射時間、跳頻與加入展頻編碼等。錯開使用頻段或發射時間無須對雷達本身多作改變,但需要事先協調各用戶。而加入跳頻機制或展頻編碼等則會大幅增加雷達訊號處理的複雜度。 An important application of FM continuous wave radar is automotive radar. The characteristics of the acquired data are complex, poor signal-to-noise ratio (SNR), and rapid changes. The above-mentioned technology is difficult to achieve high performance, low computational complexity and high real-time in such an environment. In addition, in a multi-user environment, mutual interference between radars will be a serious problem. Conventional anti-interference techniques include staggering the use of frequency bands, staggering the transmission time, frequency hopping, and adding spread spectrum coding. Staggering the use of frequency bands or transmission time does not require much change to the radar itself, but it is necessary to coordinate users in advance. The addition of frequency hopping mechanism or spread spectrum coding will greatly increase the complexity of radar signal processing.
因此,本發明係在針對上述的困擾,提出一種調頻連續波處理裝 置,以解決習知所產生的問題。 Therefore, the present invention proposes a FM continuous wave processing device for the above-mentioned problems. Set up to solve the problems caused by conventional knowledge.
本發明的主要目的,在於提供一種調頻連續波(FMCW)處理裝置,其係擬合複數個時變模型之參數,以有效提高訊雜比,並精準計算出目標物相對調頻連續波雷達之距離、速度或方位角,並有效降低資料擬合的複雜度,且提高資料處理的即時性。 The main purpose of the present invention is to provide a frequency modulated continuous wave (FMCW) processing device, which fits the parameters of multiple time-varying models to effectively improve the signal-to-noise ratio and accurately calculate the distance of the target relative to the frequency modulated continuous wave radar , Speed or azimuth angle, and effectively reduce the complexity of data fitting, and improve the real-time data processing.
本發明的另一目的,在於提供一種調頻連續波處理裝置,其係利用分群演算法對資料區分為複數個群組,並對每一個群組進行後續處理,然後將每一群組處理的結果融合起來,以有效降低複雜資料對方位估測帶來的干擾,並提高資料處理的即時性。 Another object of the present invention is to provide a FM continuous wave processing device, which uses a grouping algorithm to classify data into a plurality of groups, performs subsequent processing on each group, and then combines the results of each group processing Combine them to effectively reduce the interference caused by the position estimation of complex data and improve the real-time data processing.
本發明的再一目的,在於提供一種調頻連續波處理裝置,其係降低掃頻訊號在相干處理時間區間內出現的時間與機率,以降低多用戶干擾的機率。 Another object of the present invention is to provide a FM continuous wave processing device, which reduces the time and probability of the frequency sweep signal appearing in the coherent processing time interval, so as to reduce the probability of multi-user interference.
為達上述目的,本發明提供一種調頻連續波(FMCW)處理裝置,其包括一調頻連續波雷達與一調頻連續波處理器。調頻連續波處理器電性連接該調頻連續波雷達。調頻連續波雷達具有複數個天線,所有天線之前方有複數目標物。調頻連續波處理器電性連接調頻連續波雷達,調頻連續波處理器控制調頻連續波雷達之天線之其中一者在一相干處理時間區間內(coherent processing time interval)向所有目標物依序發射複數個掃頻訊號(chirp signals),所有目標物反射所有掃頻訊號,以形成複數個回波訊號。調頻連續波處理器透過調頻連續波雷達之所有天線接收所有回波訊號,並解調所有回波訊號,且對所有回波訊號在快時間維度(fast time dimension)、慢時間維度(slow time dimension)或空間維度(space dimension)進行傅立葉轉換(Fourier transform),以分別產生複數個距離頻域訊號、複數個速度頻域訊號或複數個方位頻域訊號,距離頻域訊號、速度 頻域訊號或方位頻域訊號為具有複數個波峰之波訊號,所有波峰之數量與所有目標物之數量相同。調頻連續波處理器利用所有距離頻域訊號、所有速度頻域訊號或所有方位頻域訊號之所有波峰建立複數個時變模型,所有時變模型之數量與所有波峰之數量相同,調頻連續波處理器根據所有時變模型計算出所有目標物相對調頻連續波雷達之距離、速度或方位角。 To achieve the above objective, the present invention provides a frequency modulated continuous wave (FMCW) processing device, which includes a frequency modulated continuous wave radar and a frequency modulated continuous wave processor. The FM continuous wave processor is electrically connected to the FM continuous wave radar. The FM continuous wave radar has multiple antennas, and there are multiple targets in front of all the antennas. The FM continuous wave processor is electrically connected to the FM continuous wave radar, and the FM continuous wave processor controls one of the antennas of the FM continuous wave radar to transmit complex numbers to all targets in a coherent processing time interval. A chirp signal, all the targets reflect all the chirp signals to form a plurality of echo signals. The FM continuous wave processor receives all the echo signals through all the antennas of the FM continuous wave radar, and demodulates all the echo signals, and for all the echo signals in the fast time dimension and slow time dimension ) Or space dimension to perform Fourier transform to generate multiple distance frequency domain signals, multiple speed frequency domain signals or multiple azimuth frequency domain signals, distance frequency domain signals, speed The frequency domain signal or the azimuth frequency domain signal is a wave signal with multiple crests, and the number of all crests is the same as the number of all targets. The FM continuous wave processor uses all the peaks of all distance frequency domain signals, all speed frequency domain signals or all azimuth frequency domain signals to create multiple time-varying models. The number of all time-varying models is the same as the number of all peaks. FM continuous wave processing The device calculates the distance, speed or azimuth angle of all targets relative to the FM continuous wave radar based on all time-varying models.
在本發明之一實施例中,調頻連續波處理器利用分群演算法(clustering algorithm)對同一距離頻域訊號進行分群,以形成複數個子距離頻域訊號,所有子距離頻域訊號分別對應複數個不同距離範圍,調頻連續波處理器依據所有距離範圍區分所有目標物為不同群組,並根據此不同群組分別對所有回波訊號在慢時間維度或空間維度進行傅立葉轉換,以分別產生所有速度頻域訊號或所有方位頻域訊號。 In an embodiment of the present invention, the FM continuous wave processor uses a clustering algorithm to cluster the same distance frequency domain signal to form a plurality of sub-distance frequency domain signals, and all the sub-distance frequency domain signals correspond to a plurality of For different distance ranges, the FM continuous wave processor distinguishes all targets into different groups according to all distance ranges, and performs Fourier transform on all echo signals in the slow time dimension or space dimension according to this different group, to generate all speeds respectively Frequency domain signal or all azimuth frequency domain signal.
在本發明之一實施例中,調頻連續波處理器計算出同一目標物對應之每一掃頻訊號及其對應之回波訊號之拍頻(beat frequency),並依據拍頻、動態規劃演算法(dynamic programming algorithm)與所有距離範圍取得所有距離頻域訊號之波峰對應之距離。 In an embodiment of the present invention, the FM continuous wave processor calculates the beat frequency of each sweep signal and its corresponding echo signal corresponding to the same target, and based on the beat frequency, the dynamic programming algorithm ( dynamic programming algorithm) Get the distance corresponding to the peak of all distance frequency domain signals with all distance ranges.
在本發明之一實施例中,調頻連續波處理器利用分群演算法對同一速度頻域訊號進行分群,以形成複數個子速度頻域訊號,所有子速度頻域訊號分別對應複數個不同速度範圍,調頻連續波處理器依據所有速度範圍區分所有目標物為不同群組,並根據此不同群組分別對所有回波訊號在快時間維度或空間維度進行傅立葉轉換,以分別產生所有距離頻域訊號或所有方位頻域訊號。 In an embodiment of the present invention, the FM continuous wave processor uses a grouping algorithm to group the same speed frequency domain signal to form a plurality of sub-speed frequency domain signals, and all the sub-speed frequency domain signals correspond to a plurality of different speed ranges, respectively. The FM continuous wave processor distinguishes all targets into different groups according to all speed ranges, and performs Fourier transforms on all echo signals in the fast time dimension or space dimension according to the different groups to generate all distance frequency domain signals or All azimuth frequency domain signals.
在本發明之一實施例中,調頻連續波處理器計算出同一目標物對應之每一掃頻訊號及其對應之回波訊號之拍頻,並依據拍頻、動態規劃演算法與所有速度範圍取得所有速度頻域訊號之波峰對應之速度。 In one embodiment of the present invention, the FM continuous wave processor calculates the beat frequency of each sweep signal and its corresponding echo signal corresponding to the same target, and obtains it according to the beat frequency, dynamic programming algorithm and all speed ranges The speed corresponding to the peak of all speed frequency domain signals.
在本發明之一實施例中,調頻連續波處理器利用分群演算法對同 一方位頻域訊號進行分群,以形成複數個子方位頻域訊號,所有子方位頻域訊號分別對應複數個不同方位範圍,調頻連續波處理器依據所有方位範圍區分所有目標物為不同群組,並根據此不同群組分別對所有回波訊號在快時間維度或慢時間維度進行傅立葉轉換,以分別產生所有距離頻域訊號或所有速度頻域訊號。 In an embodiment of the present invention, the FM continuous wave processor uses a grouping algorithm to compare the same The one-sided frequency domain signals are grouped to form a plurality of sub-azimuth frequency domain signals. All the sub-azimuth frequency domain signals correspond to a plurality of different azimuth ranges. The FM continuous wave processor distinguishes all targets into different groups according to all the azimuth ranges, and According to the different groups, Fourier transform is performed on all the echo signals in the fast time dimension or the slow time dimension to generate all distance frequency domain signals or all speed frequency domain signals respectively.
在本發明之一實施例中,調頻連續波處理器計算出同一目標物對應之每一掃頻訊號及其對應之回波訊號之拍頻,並依據拍頻、動態規劃演算法與所有方位範圍取得所有方位頻域訊號之波峰對應之方位角。 In an embodiment of the present invention, the FM continuous wave processor calculates the beat frequency of each sweep signal and its corresponding echo signal corresponding to the same target, and obtains it according to the beat frequency, dynamic programming algorithm and all azimuth ranges The azimuth angle corresponding to the peak of all azimuth frequency domain signals.
在本發明之一實施例中,分群演算法為k平均(k-means)演算法或具有雜訊之密度基礎空間分群應用法(DBSCAN,Density-based Spatial Clustering of Applications with Noise)。 In one embodiment of the present invention, the clustering algorithm is k-means algorithm or Density-based Spatial Clustering of Applications with Noise (DBSCAN).
在本發明之一實施例中,調頻連續波處理器利用一常數誤警率(Constant False-Alarm Rate,CFAR)偵測出所有波峰。 In an embodiment of the present invention, the FM continuous wave processor uses a Constant False-Alarm Rate (CFAR) to detect all peaks.
在本發明之一實施例中,調頻連續波處理器計算出同一目標物對應之每一掃頻訊號及其對應之回波訊號之拍頻,並根據拍頻對所有回波訊號在快時間維度上進行傅立葉轉換,以產生所有距離頻域訊號。 In an embodiment of the present invention, the FM continuous wave processor calculates the beat frequency of each sweep signal and its corresponding echo signal corresponding to the same target, and performs a fast time dimension on all echo signals according to the beat frequency Perform Fourier transform to generate all range frequency domain signals.
在本發明之一實施例中,調頻連續波處理器計算出同一目標物對應之每一掃頻訊號及其對應之回波訊號之拍頻,並計算出同一目標物對應之拍頻之差值,且根據差值對所有回波訊號在慢時間維度上進行傅立葉轉換,以產生所有速度頻域訊號。 In an embodiment of the present invention, the FM continuous wave processor calculates the beat frequency of each sweep signal and its corresponding echo signal corresponding to the same target, and calculates the difference between the beat frequencies corresponding to the same target, And according to the difference, all echo signals are Fourier transformed in the slow time dimension to generate all speed frequency domain signals.
在本發明之一實施例中,調頻連續波處理器取得所有天線所接收到同一目標物對應之回波訊號之時間點,並計算出此時間點之時間差,且利用時間差對所有回波訊號在空間維度上進行傅立葉轉換,以產生所有方位頻域訊號。 In an embodiment of the present invention, the FM continuous wave processor obtains the time point of the echo signal corresponding to the same target object received by all antennas, calculates the time difference at this time point, and uses the time difference to calculate the time difference between all echo signals. Fourier transformation is performed in the spatial dimension to generate all azimuth frequency domain signals.
在本發明之一實施例中,時變模型係以n次多項式(n-th degree polynomial)表示,n大於或等於1。 In an embodiment of the present invention, the time-varying model is represented by an n-th degree polynomial, and n is greater than or equal to 1.
在本發明之一實施例中,調頻連續波處理器存有對應不同方位角之複數個方位頻率響應訊號,並依據最小絕對值收斂和選擇算子(LASSO)演算法比對所有方位頻率響應訊號與所有方位頻域訊號,以取得所有方位頻域訊號之波峰對應之方位角。 In an embodiment of the present invention, the FM continuous wave processor stores a plurality of azimuth frequency response signals corresponding to different azimuth angles, and compares all azimuth frequency response signals according to the minimum absolute value convergence and selection operator (LASSO) algorithm With all azimuth frequency domain signals, to obtain the azimuth angles corresponding to the peaks of all azimuth frequency domain signals.
在本發明之一實施例中,調頻連續波處理器存有對應不同速度之複數個速度頻率響應訊號,並依據最小絕對值收斂和選擇算子(LASSO)演算法比對所有速度頻率響應訊號與所有速度頻域訊號,以取得所有速度頻域訊號之波峰對應之速度。 In an embodiment of the present invention, the FM continuous wave processor stores a plurality of speed frequency response signals corresponding to different speeds, and compares all speed frequency response signals with the minimum absolute value convergence and selection operator (LASSO) algorithm All speed frequency domain signals to obtain the speed corresponding to the peak of all speed frequency domain signals.
在本發明之一實施例中,調頻連續波處理器存有對應不同距離之複數個距離頻率響應訊號,並依據最小絕對值收斂和選擇算子(LASSO)演算法比對所有距離頻率響應訊號與所有距離頻域訊號,以取得所有距離頻域訊號之波峰對應之距離。 In one embodiment of the present invention, the FM continuous wave processor stores a plurality of distance frequency response signals corresponding to different distances, and compares all distance frequency response signals with the minimum absolute value convergence and selection operator (LASSO) algorithm All distance frequency domain signals to obtain the distance corresponding to the peak of all distance frequency domain signals.
在本發明之一實施例中,相干處理時間區間內之所有掃頻訊號之對應時間區間互相錯開,且該對應時間區間隨機分佈,相干處理時間區間為1~10毫秒(ms),所有對應時間區間佔相干處理時間區間小於或等於20%,又大於0。 In an embodiment of the present invention, the corresponding time intervals of all the frequency sweep signals in the coherent processing time interval are staggered with each other, and the corresponding time intervals are randomly distributed. The coherent processing time interval is 1-10 milliseconds (ms), and all corresponding times The interval accounts for less than or equal to 20% of the coherent processing time interval and greater than 0.
茲為使 貴審查委員對本發明的結構特徵及所達成的功效更有進一步的瞭解與認識,謹佐以較佳的實施例圖及配合詳細的說明,說明如後: In order to enable your reviewers to have a better understanding and understanding of the structural features of the present invention and the effects achieved, I would like to provide a better embodiment diagram and detailed descriptions. The description is as follows:
10:調頻連續波雷達 10: FM continuous wave radar
12:調頻連續波處理器 12: FM continuous wave processor
14:天線 14: Antenna
16:目標物 16: target
16’:目標物 16’: Target
16”:目標物 16": target
第1圖為本發明之調頻連續波(FMCW)處理裝置之一實施例之電路方塊圖。 Figure 1 is a circuit block diagram of an embodiment of the frequency modulated continuous wave (FMCW) processing device of the present invention.
第2圖為本發明之掃頻訊號之一實施例之波形圖。 Figure 2 is a waveform diagram of an embodiment of the frequency sweep signal of the present invention.
第3圖為本發明之距離頻域訊號、速度頻域訊號或方位頻域訊號之波形圖。 Figure 3 is a waveform diagram of the distance frequency domain signal, speed frequency domain signal, or azimuth frequency domain signal of the present invention.
第4圖為本發明之掃頻訊號之另一實施例之波形圖。 Figure 4 is a waveform diagram of another embodiment of the frequency sweep signal of the present invention.
於下文中關於“一個實施例”或“一實施例”之描述係指關於至少一實施例內所相關連之一特定元件、結構或特徵。因此,於下文中多處所出現之“一個實施例”或“一實施例”之多個描述並非針對同一實施例。再者,於一或多個實施例中之特定構件、結構與特徵可依照一適當方式而結合。 The following description of "one embodiment" or "an embodiment" refers to at least one specific element, structure, or feature related to the embodiment. Therefore, multiple descriptions of "one embodiment" or "an embodiment" appearing in various places below are not directed to the same embodiment. Furthermore, specific components, structures, and features in one or more embodiments can be combined in an appropriate manner.
以下請參閱第1圖、第2圖與第3圖,並介紹本發明之調頻連續波(FMCW)處理裝置,其包括一調頻連續波雷達10與一調頻連續波處理器12。調頻連續波處理器12電性連接該調頻連續波雷達10。調頻連續波雷達10可以作為車用雷達、監視雷達或合成孔徑雷達(synthetic aperture radar,SAR)。調頻連續波雷達10具有複數個天線14,所有天線14之前方有複數目標物16、16’與16”,在此實施例中,以三個目標物16、16’與16”為例。調頻連續波處理器12電性連接調頻連續波雷達10,調頻連續波處理器12控制調頻連續波雷達10之天線14之其中一者在一相干處理時間區間T0-T1內(coherent processing time interval)向所有目標物16、16’與16”依序發射複數個掃頻訊號(chirp signals),如第2圖所示,一條斜實線代表一個掃頻訊號。所有目標物16、16’與16”反射所有掃頻訊號,以形成複數個回波訊號,回波訊號以斜虛線表示。調頻連續波處理器12透過調頻連續波雷達10之所有天線14接收所有回波訊號,並解調所有回波訊號,且對所有回波訊號在快時間維度(fast time dimension)、慢時間維度(slow time dimension)或空間維度(space dimension)進行傅立葉轉換(Fourier transform),以分別產生複數個距離頻域訊號、複數個速度頻域訊號或複數個方位頻域訊號,每
一距離頻域訊號、每一速度頻域訊號或每一方位頻域訊號為具有複數個波峰A、B與C之波訊號,每一距離頻域訊號、每一速度頻域訊號或每一方位頻域訊號之所有波峰A、B與C之數量與所有目標物16之數量相同,且波峰A、B與C分別對應目標物16、16’與16”。由於波訊號會受雜訊干擾,故波訊號並非所有波峰都能夠代表目標物16、16’與16”,只有比較明顯的波峰A、B與C可以代表目標物16、16’與16”,因此,如第2圖所示,調頻連續波處理器12利用一常數誤警率(Constant False-Alarm Rate,CFAR)偵測出波峰A、B與C。調頻連續波處理器12利用所有距離頻域訊號、所有速度頻域訊號或所有方位頻域訊號之所有波峰A、B與C建立複數個時變(time varying)模型,所有時變模型之數量與所有波峰A、B與C之數量相同。一個時變模型代表一個目標物16、16’或16”的運動行為,因為目標物16、16’或16”的運動行為可以用多項式表示,故時變模型係以n次多項式(n-th degree polynomial)表示,n大於或等於1。調頻連續波處理器12根據所有時變模型計算出所有目標物16、16’與16”相對調頻連續波雷達10之距離、速度或方位角。
Please refer to FIG. 1, FIG. 2 and FIG. 3 below to introduce the frequency modulated continuous wave (FMCW) processing device of the present invention, which includes a frequency modulated
在本發明之一實施例中,調頻連續波處理器12計算出同一目標物16、16’或16”對應之每一掃頻訊號及其對應之回波訊號之拍頻(beat frequency)Fb,並根據拍頻Fb對所有回波訊號在快時間維度上進行傅立葉轉換,以產生所有距離頻域訊號。調頻連續波處理器12利用分群演算法(clustering algorithm)對同一距離頻域訊號進行分群,以形成複數個子距離頻域訊號,所有子距離頻域訊號分別對應複數個不同距離範圍,調頻連續波處理器12依據所有距離範圍區分所有目標物16、16’與16”為不同群組,並根據此不同群組分別對所有回波訊號在慢時間維度或空間維度進行傅立葉轉換,以分別產生所有速度頻域訊號或所有方位頻域訊號。在本發明之一實施例中,分群演算法可為k平均(k-means)演算法或具有雜訊之密度基礎空間分群應用法(DBSCAN,
Density-based Spatial Clustering of Applications with Noise)。以第2圖為例,同一距離頻域訊號分為二個子距離頻域訊號,其中一距離頻域訊號對應一距離範圍,另一距離頻域訊號對應另一距離範圍。二不同距離範圍對應不同的目標物16、16’與16”,其中一距離範圍對應一目標物16,另一距離範圍對應二目標物16’與16”。因此,調頻連續波處理器12可以先對目標物16之所有回波訊號在慢時間維度或空間維度進行傅立葉轉換,以分別產生所有速度頻域訊號或所有方位頻域訊號。接著,調頻連續波處理器12再對目標物16’與16”之所有回波訊號在慢時間維度或空間維度進行傅立葉轉換,以分別產生所有速度頻域訊號或所有方位頻域訊號。本發明利用分群演算法對資料區分為複數個群組,並對每一個群組進行後續處理,然後將每一群組處理的結果融合起來,以有效降低複雜資料對方位估測帶來的干擾,並提高資料處理的即時性。
In one embodiment of the present invention, the FM
在本發明之一實施例中,調頻連續波處理器12計算出同一目標物16、16’或16”對應之每一掃頻訊號及其對應之回波訊號之拍頻Fb,並計算出同一目標物16、16’或16”對應之拍頻Fb之差值,且根據此差值對所有回波訊號在慢時間維度上進行傅立葉轉換,以產生所有速度頻域訊號。調頻連續波處理器12利用分群演算法對同一速度頻域訊號進行分群,以形成複數個子速度頻域訊號,所有子速度頻域訊號分別對應複數個不同速度範圍。調頻連續波處理器12依據所有速度範圍區分所有目標物16、16’與16”為不同群組,並根據此不同群組分別對所有回波訊號在快時間維度或空間維度進行傅立葉轉換,以分別產生所有距離頻域訊號或所有方位頻域訊號。在本發明之一實施例中,分群演算法可為k平均(k-means)演算法或具有雜訊之密度基礎空間分群應用法(DBSCAN,Density-based Spatial Clustering of Applications with Noise)。以第2圖為例,同一速度頻域訊號分為二個子速度頻域訊號,其中一速度頻域訊號對應一速度範圍,另一速度頻域訊號對應另一速度範圍。二不同速度範圍對應不同的目標物
16、16’與16”,其中一速度範圍對應一目標物16,另一速度範圍對應二目標物16’與16”。因此,調頻連續波處理器12可以先對目標物16之所有回波訊號在快時間維度或空間維度進行傅立葉轉換,以分別產生所有距離頻域訊號或所有方位頻域訊號。接著,調頻連續波處理器12再對目標物16’與16”之所有回波訊號在快時間維度或空間維度進行傅立葉轉換,以分別產生所有距離頻域訊號或所有方位頻域訊號。本發明利用分群演算法對資料區分為複數個群組,並對每一個群組進行後續處理,然後將每一群組處理的結果融合起來,以有效降低複雜資料對方位估測帶來的干擾,並提高資料處理的即時性。
In an embodiment of the present invention, the FM
在本發明之一實施例中,調頻連續波處理器12取得所有天線14所接收到同一目標物16、16’或16”對應之回波訊號之時間點,並計算出所有時間點之時間差,且利用所有時間差對所有回波訊號在空間維度上進行傅立葉轉換,以產生所有方位頻域訊號。調頻連續波處理器12利用分群演算法對同一方位頻域訊號進行分群,以形成複數個子方位頻域訊號,所有子方位頻域訊號分別對應複數個不同方位範圍。調頻連續波處理器12依據所有方位範圍區分所有目標物16、16’與16”為不同群組,並根據此不同群組分別對所有回波訊號在快時間維度或慢時間維度進行傅立葉轉換,以分別產生所有距離頻域訊號或所有速度頻域訊號。在本發明之一實施例中,分群演算法可為k平均(k-means)演算法或具有雜訊之密度基礎空間分群應用法(DBSCAN,Density-based Spatial Clustering of Applications with Noise)。以第2圖為例,同一方位頻域訊號分為二個子方位頻域訊號,其中一方位頻域訊號對應一方位範圍,另一方位頻域訊號對應另一方位範圍。二不同方位範圍對應不同的目標物16、16’與16”,其中一方位範圍對應一目標物16,另一方位範圍對應二目標物16’與16”。因此,調頻連續波處理器12可以先對目標物16之所有回波訊號在快時間維度或慢時間維度進行傅立葉轉換,以分別產生所有距離頻域訊號或所有速度頻域訊號。接著,
調頻連續波處理器12再對目標物16’與16”之所有回波訊號在快時間維度或慢時間維度進行傅立葉轉換,以分別產生所有距離頻域訊號或所有速度頻域訊號。本發明利用分群演算法對資料區分為複數個群組,並對每一個群組進行後續處理,然後將每一群組處理的結果融合起來,以有效提高資料處理的即時性。
In an embodiment of the present invention, the FM
以下介紹本發明取得波峰A、B與C對應之距離、速度與方位角的方式。若未取得同一目標物16、16’或16”對應之每一掃頻訊號及其對應之回波訊號之拍頻Fb時,則調頻連續波處理器12計算出同一目標物16、16’或16”對應之每一掃頻訊號及其對應之回波訊號之拍頻Fb。若已取得同一目標物16、16’或16”對應之每一掃頻訊號及其對應之回波訊號之拍頻Fb時,則調頻連續波處理器12可直接使用同一目標物16、16’或16”對應之每一掃頻訊號及其對應之回波訊號之拍頻Fb。同一目標物16、16’或16”對應之每一掃頻訊號及其對應之回波訊號之拍頻Fb可能呈特定關係,例如為線性關係,則調頻連續波處理器12依據同一目標物16、16’或16”對應之每一掃頻訊號及其對應之回波訊號之拍頻Fb、動態規劃演算法(dynamic programming algorithm)與所有不同距離範圍取得所有距離頻域訊號之波峰A、B與C對應之距離。除了上述方式之外,調頻連續波處理器12可預存有對應不同距離之複數個距離頻率響應訊號,並依據最小絕對值收斂和選擇算子(LASSO)演算法比對所有距離頻率響應訊號與距離頻域訊號,以取得距離頻域訊號之波峰對應之距離。假設1公里對應的距離頻率響應訊號為,是方位角,d是天線之間的間距,f c 是中心頻率,c則是光速。若具波峰A的距離頻域訊號與此距離頻率響應訊號相同時,代表此距離頻域訊號之波峰A所對應的距離為1公里。依此類推,調頻連續波處理器12取得所有距離頻域訊號之波峰A、B與C所對應的距離。
The following describes the method of the present invention to obtain the distances, speeds and azimuths corresponding to the peaks A, B and C. If the beat frequency Fb of each sweep signal and its corresponding echo signal corresponding to the
調頻連續波處理器12將所有波峰A對應之距離代入n次多項式中可以取得目標物16在距離上之時變模型及其參數,將所有波峰B對應之距離代入
n次多項式中可以取得目標物16’在距離上之時變模型及其參數,將所有波峰C對應之距離代入n次多項式中可以取得目標物16”在距離上之時變模型及其參數。調頻連續波處理器12根據所有時變模型計算出所有目標物16、16’與16”相對調頻連續波雷達10之距離。本發明擬合複數個時變模型之參數,以有效提高訊雜比,並精準計算出目標物16、16’與16”相對調頻連續波雷達10之距離,並有效降低資料擬合的複雜度,且提高資料處理的即時性。
The FM
若未取得同一目標物16、16’或16”對應之每一掃頻訊號及其對應之回波訊號之拍頻Fb時,則調頻連續波處理器12計算出同一目標物16、16’或16”對應之每一掃頻訊號及其對應之回波訊號之拍頻Fb。若已取得同一目標物16、16’或16”對應之每一掃頻訊號及其對應之回波訊號之拍頻Fb時,則調頻連續波處理器12可直接使用同一目標物16、16’或16”對應之每一掃頻訊號及其對應之回波訊號之拍頻Fb。同一目標物16、16’或16”對應之每一掃頻訊號及其對應之回波訊號之拍頻Fb可能呈特定關係,例如為線性關係,則調頻連續波處理器12依據同一目標物16、16’或16”對應之每一掃頻訊號及其對應之回波訊號之拍頻Fb、動態規劃演算法與所有不同速度範圍取得所有速度頻域訊號之波峰A、B與C對應之速度。除了上述方式之外,調頻連續波處理器12可預存有對應不同速度之複數個速度頻率響應訊號,並依據最小絕對值收斂和選擇算子(LASSO)演算法比對所有速度頻率響應訊號與速度頻域訊號,以取得速度頻域訊號之波峰對應之速度。假設1公里/分對應的速度頻率響應訊號為,是方位角,d是天線之間的間距,f c 是中心頻率,c則是光速。若具波峰A的速度頻域訊號與此速度頻率響應訊號相同時,代表此速度頻域訊號之波峰A所對應的速度為1公里/分。依此類推,調頻連續波處理器12取得所有速度頻域訊號之波峰A、B與C所對應的速度。
If the beat frequency Fb of each sweep signal and its corresponding echo signal corresponding to the
調頻連續波處理器12將所有波峰A對應之速度代入n次多項式中
可以取得目標物16在速度上之時變模型及其參數,將所有波峰B對應之速度代入n次多項式中可以取得目標物16’在速度上之時變模型及其參數,將所有波峰C對應之速度代入n次多項式中可以取得目標物16”在速度上之時變模型及其參數。調頻連續波處理器12根據所有時變模型計算出所有目標物16、16’與16”相對調頻連續波雷達10之速度。本發明擬合複數個時變模型之參數,以有效提高訊雜比,並精準計算出目標物16、16’與16”相對調頻連續波雷達10之速度,並有效降低資料擬合的複雜度,且提高資料處理的即時性。
The FM
若未取得同一目標物16、16’或16”對應之每一掃頻訊號及其對應之回波訊號之拍頻Fb時,則調頻連續波處理器12計算出同一目標物16、16’或16”對應之每一掃頻訊號及其對應之回波訊號之拍頻Fb。若已取得同一目標物16、16’或16”對應之每一掃頻訊號及其對應之回波訊號之拍頻Fb時,則調頻連續波處理器12可直接使用同一目標物16、16’或16”對應之每一掃頻訊號及其對應之回波訊號之拍頻Fb。同一目標物16、16’或16”對應之每一掃頻訊號及其對應之回波訊號之拍頻Fb可能呈特定關係,例如為線性關係,則調頻連續波處理器12依據同一目標物16、16’或16”對應之每一掃頻訊號及其對應之回波訊號之拍頻Fb、動態規劃演算法與所有不同方位範圍取得所有方位頻域訊號之波峰A、B與C對應之方位。除了上述方式之外,調頻連續波處理器12可預存有對應不同方位之複數個方位頻率響應訊號,並依據最小絕對值收斂和選擇算子(LASSO)演算法比對所有方位頻率響應訊號與方位頻域訊號,以取得方位頻域訊號之波峰對應之方位角。假設15度對應的方位頻率響應訊號為,是方位角,d是天線之間的間距,f c 是中心頻率,c則是光速。若具波峰A的方位頻域訊號與此方位頻率響應訊號相同時,代表此方位頻域訊號之波峰A所對應的方位為15度。依此類推,調頻連續波處理器12取得所有方位頻域訊號之波峰A、B與C所對應的方位角。
If the beat frequency Fb of each sweep signal and its corresponding echo signal corresponding to the
調頻連續波處理器12將所有波峰A對應之方位角代入n次多項式中可以取得目標物16在方位上之時變模型及其參數,將所有波峰B對應之方位角代入n次多項式中可以取得目標物16’在方位上之時變模型及其參數,將所有波峰C對應之方位角代入n次多項式中可以取得目標物16”在方位上之時變模型及其參數。調頻連續波處理器12根據所有時變模型計算出所有目標物16、16’與16”相對調頻連續波雷達10之方位角。本發明擬合複數個時變模型之參數,以有效提高訊雜比,並精準計算出目標物16、16’與16”相對調頻連續波雷達10之方位角,並有效降低資料擬合的複雜度,且提高資料處理的即時性。
The FM
在本發明之一實施例中,如第4圖所示,一條斜實線代表一個掃頻訊號,回波訊號以斜虛線表示,相干處理時間區間T0-T1內之所有掃頻訊號之對應時間區間互相錯開,並避免重疊,且所有對應時間區間隨機分佈,相干處理時間區間為1~10毫秒(ms),所有對應時間區間佔相干處理時間區間T0-T1小於或等於20%,又大於0。本發明降低掃頻訊號在相干處理時間區間T0-T1內出現的時間與機率,以降低後續處理的干擾,同時降低多用戶干擾的機率,提升估測距離、速度或方位角之準確性。 In an embodiment of the present invention, as shown in Figure 4, an oblique solid line represents a sweep signal, and the echo signal is represented by an oblique dashed line. The corresponding time of all sweep signals in the coherent processing time interval T0-T1 The intervals are staggered to avoid overlapping, and all corresponding time intervals are randomly distributed. The coherent processing time interval is 1-10 milliseconds (ms). All corresponding time intervals account for less than or equal to 20% of the coherent processing time interval T0-T1, and greater than 0 . The present invention reduces the time and probability of the frequency sweep signal appearing in the coherent processing time interval T0-T1 to reduce subsequent processing interference, while reducing the probability of multi-user interference, and improving the accuracy of estimated distance, speed or azimuth.
綜上所述,本發明擬合複數個時變模型之參數,以有效提高訊雜比,並精準計算出目標物相對調頻連續波雷達之距離、速度或方位角,降低資料擬合的複雜度,且提高資料處理的即時性。 In summary, the present invention fits the parameters of multiple time-varying models to effectively improve the signal-to-noise ratio, and accurately calculates the distance, speed or azimuth of the target relative to the FM continuous wave radar, reducing the complexity of data fitting , And improve the real-time data processing.
以上所述者,僅為本發明一較佳實施例而已,並非用來限定本發明實施之範圍,故舉凡依本發明申請專利範圍所述之形狀、構造、特徵及精神所為之均等變化與修飾,均應包括於本發明之申請專利範圍內。 The above is only a preferred embodiment of the present invention, and is not used to limit the scope of implementation of the present invention. Therefore, all the shapes, structures, features and spirits described in the scope of the patent application of the present invention are equally changed and modified. , Should be included in the scope of patent application of the present invention.
10:調頻連續波雷達 10: FM continuous wave radar
12:調頻連續波處理器 12: FM continuous wave processor
14:天線 14: Antenna
16:目標物 16: target
16’:目標物 16’: Target
16”:目標物 16": target
Claims (17)
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| TW108130214A TWI697688B (en) | 2019-08-23 | 2019-08-23 | Frequency modulated continuous wave processing device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| TW108130214A TWI697688B (en) | 2019-08-23 | 2019-08-23 | Frequency modulated continuous wave processing device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| TWI697688B true TWI697688B (en) | 2020-07-01 |
| TW202109077A TW202109077A (en) | 2021-03-01 |
Family
ID=72601812
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| TW108130214A TWI697688B (en) | 2019-08-23 | 2019-08-23 | Frequency modulated continuous wave processing device |
Country Status (1)
| Country | Link |
|---|---|
| TW (1) | TWI697688B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114076941A (en) * | 2020-08-21 | 2022-02-22 | 上海禾赛科技有限公司 | Method, radar, and computer-readable storage medium for detection using frequency-modulated continuous waves |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20130009807A1 (en) * | 2011-07-07 | 2013-01-10 | Lamb Brian M | Apparatus & Method for Short Dwell Inverse Synthetic Aperture Radar (ISAR) Imaging of Turning Moving Vehicles |
| CN106646447A (en) * | 2017-01-18 | 2017-05-10 | 武汉雷博合创电子技术有限公司 | Detection method for radar target long-time accumulation based on linear frequency modulation continuous wave |
| TWI586987B (en) * | 2016-12-22 | 2017-06-11 | Nat Chung-Shan Inst Of Science And Tech | Signal processing device for continuous wave radar sensing system |
| TW201723529A (en) * | 2015-12-28 | 2017-07-01 | Panasonic Ip Man Co Ltd | Sensor and faucet device using the sensor comprising a detection unit and a processing unit with a correction unit and acalculation unit |
| EP3425419A1 (en) * | 2017-07-05 | 2019-01-09 | Stichting IMEC Nederland | A method and a system for localization and monitoring of living beings |
-
2019
- 2019-08-23 TW TW108130214A patent/TWI697688B/en active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20130009807A1 (en) * | 2011-07-07 | 2013-01-10 | Lamb Brian M | Apparatus & Method for Short Dwell Inverse Synthetic Aperture Radar (ISAR) Imaging of Turning Moving Vehicles |
| TW201723529A (en) * | 2015-12-28 | 2017-07-01 | Panasonic Ip Man Co Ltd | Sensor and faucet device using the sensor comprising a detection unit and a processing unit with a correction unit and acalculation unit |
| TWI586987B (en) * | 2016-12-22 | 2017-06-11 | Nat Chung-Shan Inst Of Science And Tech | Signal processing device for continuous wave radar sensing system |
| CN106646447A (en) * | 2017-01-18 | 2017-05-10 | 武汉雷博合创电子技术有限公司 | Detection method for radar target long-time accumulation based on linear frequency modulation continuous wave |
| EP3425419A1 (en) * | 2017-07-05 | 2019-01-09 | Stichting IMEC Nederland | A method and a system for localization and monitoring of living beings |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114076941A (en) * | 2020-08-21 | 2022-02-22 | 上海禾赛科技有限公司 | Method, radar, and computer-readable storage medium for detection using frequency-modulated continuous waves |
Also Published As
| Publication number | Publication date |
|---|---|
| TW202109077A (en) | 2021-03-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US11693104B2 (en) | Radar angular ambiguity resolution | |
| US11988736B2 (en) | Systems and methods for virtual aperture radar tracking | |
| CN112630768B (en) | Noise reduction method for improving frequency modulation continuous wave radar target detection | |
| CN114488109B (en) | Low complexity super resolution technique for object detection in frequency modulated continuous wave radar | |
| CN109975807B (en) | A dimension-reduced subspace angle measurement method suitable for millimeter-wave vehicle radar | |
| US10845886B2 (en) | Coherent multi-look radar processing | |
| CN109324322B (en) | A Direction Finding and Target Recognition Method Based on Passive Phased Array Antenna | |
| Wagner et al. | Modification of DBSCAN and application to range/Doppler/DoA measurements for pedestrian recognition with an automotive radar system | |
| CN112241003B (en) | Method and system for object detection | |
| CN110632587A (en) | Weak moving object monitoring method based on rapid FMCW radar | |
| JP2015514971A (en) | Target detection method | |
| Wang et al. | An antideceptive jamming method for multistatic synthetic aperture radar based on collaborative localization and spatial suppression | |
| CN114114153B (en) | A multi-sound source localization method, system, microphone array and terminal device | |
| Lim et al. | Radar sensor-based estimation of vehicle orientation for autonomous driving | |
| JP2012194051A (en) | Electronic scanning radar device, received wave direction estimation method, and program | |
| Cherniakov et al. | Forward scattering micro sensor for vehicle classification | |
| TWI697688B (en) | Frequency modulated continuous wave processing device | |
| JP7224292B2 (en) | Radar device and automobile equipped with it | |
| Kaiser et al. | Complex-valued neural networks for millimeter wave fmcw-radar angle estimations | |
| CN106291450B (en) | Direction finding method and device for a single-transmission and dual-reception passive radar | |
| US10845475B2 (en) | Method of measuring azimuth of radar target | |
| CN1063268C (en) | Radar apparatus | |
| CN115616562A (en) | Apparatus and method for detecting target using radar | |
| CN119104990A (en) | Interference target detection method, device, radar equipment and storage medium | |
| Ciattaglia et al. | mmWave radar features extraction of drones for machine learning classification |