TWI762243B - 頻率調變連續波雷達、數位訊號處理方法與表徵資訊偵測方法 - Google Patents
頻率調變連續波雷達、數位訊號處理方法與表徵資訊偵測方法 Download PDFInfo
- Publication number
- TWI762243B TWI762243B TW110109640A TW110109640A TWI762243B TW I762243 B TWI762243 B TW I762243B TW 110109640 A TW110109640 A TW 110109640A TW 110109640 A TW110109640 A TW 110109640A TW I762243 B TWI762243 B TW I762243B
- Authority
- TW
- Taiwan
- Prior art keywords
- signal
- signals
- frequency
- digital
- output signal
- Prior art date
Links
- 238000012545 processing Methods 0.000 title claims abstract description 48
- 238000012512 characterization method Methods 0.000 title claims abstract description 9
- 238000001514 detection method Methods 0.000 title claims description 111
- 238000000034 method Methods 0.000 title abstract description 20
- 238000003672 processing method Methods 0.000 claims description 24
- 238000010801 machine learning Methods 0.000 claims description 17
- 238000004458 analytical method Methods 0.000 claims description 4
- 230000004044 response Effects 0.000 claims description 3
- 230000007613 environmental effect Effects 0.000 abstract description 16
- 238000010586 diagram Methods 0.000 description 38
- 230000008569 process Effects 0.000 description 14
- 239000011159 matrix material Substances 0.000 description 13
- 230000033001 locomotion Effects 0.000 description 9
- 230000003068 static effect Effects 0.000 description 8
- 230000005540 biological transmission Effects 0.000 description 6
- 238000009826 distribution Methods 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 4
- 230000008859 change Effects 0.000 description 4
- 238000012937 correction Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000036544 posture Effects 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 238000012549 training Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 238000003491 array Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 230000029058 respiratory gaseous exchange Effects 0.000 description 1
- 230000036387 respiratory rate Effects 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/35—Details of non-pulse systems
- G01S7/352—Receivers
- G01S7/354—Extracting wanted echo-signals
-
- 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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/28—Details of pulse systems
-
- 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/56—Discriminating between fixed and moving objects or between objects moving at different speeds for presence detection
-
- 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/58—Velocity or trajectory determination systems; Sense-of-movement determination systems
- G01S13/583—Velocity or trajectory determination systems; Sense-of-movement determination systems using transmission of continuous unmodulated waves, amplitude-, frequency-, or phase-modulated waves and based upon the Doppler effect resulting from movement of targets
- G01S13/584—Velocity or trajectory determination systems; Sense-of-movement determination systems using transmission of continuous unmodulated waves, amplitude-, frequency-, or phase-modulated waves and based upon the Doppler effect resulting from movement of targets adapted for simultaneous range and velocity measurements
-
- 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/88—Radar or analogous systems specially adapted for specific applications
-
- 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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
- G01S7/415—Identification of targets based on measurements of movement associated with the target
-
- 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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
- G01S7/417—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section involving the use of neural networks
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
-
- 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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/35—Details of non-pulse systems
- G01S7/352—Receivers
- G01S7/356—Receivers involving particularities of FFT processing
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Artificial Intelligence (AREA)
- Evolutionary Computation (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
一種數位訊號處理方法,係疊加接收端收取雷達訊號所獲得的數位訊號,而可減少環境干擾造成的雜訊。因此,依據疊加後所得到的輸出訊號,可對偵測目標量測得精準的表徵資訊。
Description
本發明關於雷達訊號處理技術,特別是一種頻率調變連續波雷達、數位訊號處理方法與表徵資訊偵測方法。
雷達技術可用於許多室外場域應用,例如測速、測距。然而,在室內場域中,具有許多環境干擾因素,造成量測結果不準確。例如:地板、牆壁、櫥櫃造成信號反射干擾,窗廉擺動、風扇轉動等物體移動也會造成信號擾動。
有鑑於此,依據一些實施例,頻率調變連續波雷達包括發射單元及接收單元。發射單元用以發送複數啁啾訊號。接收單元接收反射的啁啾訊號,產生對應啁啾訊號的複數數位訊號。接收單元再疊加數位訊號而得到輸出訊號,並依據輸出訊號計算偵測目標的表徵資訊。
依據一些實施例,接收單元將位於同一訊框中的啁啾訊號所對應的數位訊號疊加而得到輸出訊號。
依據一些實施例,接收單元將相鄰的複數個訊框中的啁啾訊號所對應的數位訊號疊加而得到輸出訊號。
依據一些實施例,接收單元先對位數位訊號再對疊加。
依據一些實施例,發射單元包括至少一發射天線。接收單元疊加來自同一發射天線的啁啾訊號所對應的數位訊號。
依據一些實施例,接收單元包括複數接收天線。接收單元將其中一個接收天線所收到的來自同一發射天線的啁啾訊號所對應的數位訊號進行疊加而得到輸出訊號。
依據一些實施例,接收單元包括複數接收天線。接收單元疊加此些接收天線所收到的來自同一發射天線的啁啾訊號所對應的數位訊號而得到輸出訊號。
依據一些實施例,發射單元包括複數發射天線。接收單元疊加來自此些發射天線中的至少二者的啁啾訊號所對應的數位訊號。
依據一些實施例,接收單元包括複數接收天線。接收單元將其中一個接收天線所收到的來自此些發射天線中的至少二者的啁啾訊號所對應的數位訊號進行疊加而得到輸出訊號。
依據一些實施例,接收單元包括複數接收天線。接收單元疊加此些接收天線所收到的來自此些發射天線中的至少二者的啁啾訊號所對應的數位訊號而得到輸出訊號。
依據一些實施例,接收單元還根據輸出訊號產生統計資訊,並依據統計資訊修正輸出訊號。
依據一些實施例,統計資訊為絕對平方和、絕對最大值、標準差或變異數。
依據一些實施例,接收單元根據輸出訊號及機器學習模型判斷偵測區域中是否存在偵測目標。
依據一些實施例,接收單元於判斷到偵測目標存在時,依據輸出訊號計算偵測目標的表徵資訊。
依據一些實施例,數位訊號處理方法由訊號處理裝置中之處理器執行,包括:疊加對應於都卜勒雷達的接收端收到的複數啁啾訊號的複數數位訊號而得到輸出訊號;及依據輸出訊號計算偵測目標的表徵資訊。
依據一些實施例,所疊加的啁啾訊號所對應的數位訊號是位於同一訊框中或位於相鄰的複數個訊框中。
依據一些實施例,在疊加數位訊號之前,先對位待疊加的數位訊號。
依據一些實施例,數位訊號處理方法還包括:根據輸出訊號產生統計資訊;及依據統計資訊修正輸出訊號。
依據一些實施例,數位訊號處理方法還包括:根據輸出訊號及機器學習模型判斷偵測區域中是否存在偵測目標。
依據一些實施例,依據該輸出訊號計算偵測目標的表徵資訊的步驟是於判斷到偵測目標存在時執行。
依據一些實施例,表徵資訊偵測方法包括:接收對應於都卜勒雷達的複數數位偵測訊號;對數位偵測訊號進行頻域分析,以得到複數頻域偵測訊號;依據頻域偵測訊號產生複數統計資訊;依據統計資訊修正頻域偵測訊號;根據修正後的頻域偵測訊號及機器學習模型判斷偵測區域中是否存在偵測目標;及響應於存在偵測目標,依據修正後的頻域偵測訊號計算偵測目標的表徵資訊。
綜上所述,依據一些實施例的頻率調變連續波雷達、數位訊號處理方法與表徵資訊偵測方法,能解決靜態與動態環境干擾造成的訊號不佳的問題,並可減少資料處理量而可加快處理速度。依據一些實施例的頻率調變連續波雷達、數位訊號處理方法與表徵資訊偵測方法,能識別偵測目標存在與否,以增進處理效率並過濾誤判錯誤。
參照圖1,係為依據一些實施例的頻率調變連續波(Frequency Modulated Continuous Wave,FMCW)雷達10的使用狀態示意圖。頻率調變連續波雷達10可對一偵測區域40進行偵測。偵測結果可供計算偵測目標50的一種或多種資訊,例如:距離、方位、移動速度、生理資訊(如心跳、呼吸)等。偵測目標50可以是例如但不限於生物體(如人體)。
參照圖2,係為依據一些實施例的頻率調變連續波雷達10的方塊示意圖。頻率調變連續波雷達10包括發射單元20及接收單元30。發射單元20發送雷達訊號,接收單元30接收反射的雷達訊號。
參照圖3,圖3為例示雷達訊號的示意圖,上半部呈現雷達訊號的振幅對時間的變化,下半部呈現雷達訊號的頻率對時間的變化。發射單元20發送的雷達訊號包括複數啁啾(chirp)訊號SC。為了圖式清晰,圖3僅呈現一個啁啾訊號SC。在此,啁啾訊號SC為線性調頻脈衝訊號,指頻率隨時間以線性方式增加的正弦波。在一些實施例中,啁啾訊號SC的頻率是以非線性方式增加。為了方便說明,後續以線性方式來說明。如圖3所示,在一持續時間Tc(如40微秒)內,啁啾訊號SC根據一斜率S由一起始頻率(如77GHz)線性增加至一終止頻率(如81GHz)。起始頻率與終止頻率可選自毫米波頻段(即30GHz至300GHz)。起始頻率與終止頻率之差為脈衝帶寬B。
合併參照圖4及圖5。圖4為依據一些實施例的頻率調變連續波雷達10的細部方塊示意圖。圖5為例示發送與接收的雷達訊號的示意圖。發射單元20包括發射天線22和雷達發射器24。雷達發射器24包括訊號合成器,用以產生啁啾訊號Ct,並經由發射天線22發射。接收單元30包括接收天線32、雷達接收器34及處理單元36。接收天線32接收反射的雷達訊號(啁啾訊號Cr)。啁啾訊號Cr可視為啁啾訊號Ct的延遲版本。雷達接收器34包括混頻器、低通濾波器及類比數位轉換器。混頻器將來自雷達發射器24的啁啾訊號Ct和接收到的啁啾訊號Cr耦合,可產生兩啁啾訊號Ct、Cr的頻率之和以及頻率之差等兩種耦合訊號。低通濾波器將耦合後的訊號進行低通濾波而獲得兩啁啾訊號Ct、Cr的頻率之差的耦合訊號,於後稱「中頻(Intermediate Frequency)訊號SI」。類比數位轉換器將中頻訊號SI轉換為數位訊號,以利處理單元36進行數位訊號處理。處理單元36可以例如是中央處理單元(Central Processing Unit,CPU)、圖形處理器(Graphics Processing Unit,GPU),或是其他可程式化之一般用途或特殊用途的微處理器(Microprocessor)、數位訊號處理器(Digital Signal Processor,DSP)、可程式化控制器、特殊應用積體電路(Application Specific Integrated Circuits,ASIC)、可程式化邏輯裝置(Programmable Logic Device,PLD)或其他類似裝置、晶片、積體電路及其組合。
在本揭露之另一實施例中,頻率調變連續波雷達10更包括一傳輸模組,連接至處理單元36。該傳輸模組用以將處理單元36進行數位訊號處理之結果傳送至另一端之邊緣裝置或是雲端伺服器。
在本揭露之另一實施例中,頻率調變連續波雷達10的處理單元36僅對來自類比數位轉換器的數位訊號進行部分處理,而透過頻率調變連續波雷達10的一傳輸模組將部分處理的結果傳送至另一端之邊緣裝置或是雲端伺服器上進行後續數位訊號處理及運算。
在本揭露之另一實施例中,頻率調變連續波雷達10的處理單元36不對來自類比數位轉換器的數位訊號進行任何處理,而透過頻率調變連續波雷達10的一傳輸模組直接將來自類比數位轉換器的數位訊號傳送至另一端之邊緣裝置或是雲端伺服器上進行數位訊號處理及運算。
參照圖5,中頻訊號SI的頻率
可表示為式1,S為斜率,τ為雷達訊號發送至接收之間的延遲時間。因此,τ可表示為式2,d為頻率調變連續波雷達10之發射天線與待測對象間的距離,c為光速。將式2代入式1可獲得式3。由式3可以知道,中頻訊號SI的頻率
隱含有距離資訊(即頻率調變連續波雷達10與偵測目標50之間的距離)。
參照圖6,係為依據一些實施例的訊號處理示意圖。在此,將啁啾訊號SC分別依序編號為C1、C2、C3、…、Cn,n為正整數。雷達接收器34將所收到的對應各啁啾訊號C1~Cn的中頻訊號SI轉換為數位訊號SD(分別表示為D1、D2…、Dn,n為正整數)。各數位訊號SD的數值可表示為一維陣列(橫列(Row)矩陣)。將此些橫列矩陣依序縱向排列而可成為二維陣列A1。可以理解的是,也可以將數位訊號SD排列成直行(Column)矩陣並依序橫向排列,同樣可以得到二維陣列。二維陣列A1的值代表訊號強度(振幅)。二維陣列A1的直行的索引值(Index)對應於啁啾訊號SC(數位訊號SD)的次序。二維陣列A1的橫列的索引值具有時間的含意,亦即二維陣列A1的橫列矩陣為時域訊號。
處理單元36對二維陣列A1(即數位訊號SD)的各個橫列矩陣執行快速傅立葉轉換(Fast Fourier Transform,FFT)(後稱「距離傅立葉轉換」)而可取得頻域訊號SF(分別表示為F1、F2…、Fn,n為正整數),即二維陣列A2。因此,二維陣列A2的橫列矩陣相當於頻譜分布。如前述,中頻訊號SI的頻率隱含有距離資訊。亦即,二維陣列A2的橫列的索引值具有距離含意。二維陣列A2的值代表頻譜上各頻率的強度,可呈現相距頻率調變連續波雷達10不同距離所反射的雷達訊號強度。如圖6所示,二維陣列A2中的填色框為峰值處(即數值超過一閾值),表示在此頻率對應距離處有物體。從峰值處的頻率可以計算出頻率調變連續波雷達10與偵測目標50之間的距離。
由於各啁啾訊號SC之間的間隔時間很短(例如為數十微秒)。相對而言,反射此些啁啾訊號SC的同一物體的位置基本不變。因此,各個頻域訊號SF具有對應於同一距離的填色框,而呈現一直行填色框。如圖6所示,於此例為兩直行填色框。雖然各頻域訊號SF在同一直行出現峰值,些微的運動變化無法在頻率上看到明顯改變,然而相位成分卻會造成明顯的影響。二維陣列A2的直行索引值對應於頻域訊號SF(啁啾訊號SC)的次序,意味著為時間次序。因此,二維陣列A2的各個直行矩陣可視為時域訊號。處理單元36對二維陣列A2的各個直行矩陣執行快速傅立葉轉換(後稱「都卜勒傅立葉轉換」),而可分別取得相位頻域訊號SQ(分別表示為Q1、Q2、…、Qm,m為正整數),即二維陣列A3。因此,二維陣列A3的直行矩陣相當於相位頻譜分布。中頻訊號SI的相位
可表示為式4,經代入式2,可表示為式5,
為波長。根據式5,可以推導出式6,v為速度,
為兩相鄰的啁啾訊號(Cn-1、Cn)之間的相位差,
為兩相鄰的啁啾訊號SC之間的時間差。由式6可以知道,中頻訊號SI的相位隱含有運動資訊(速度)。因此,二維陣列A3的直行的索引值具有速度含意。從相位頻域訊號SQ可以計算出偵測目標50的移動速度或是週期性運動的頻率,而可得到偵測目標50的表徵資訊(如運動資訊、生理資訊(如呼吸頻率、心跳頻率))。如圖6的二維陣列A3所示為例,在兩個直行分別有二個填色框和三個填色框,表示在兩個相距頻率調變連續波雷達10的不同距離處,具有至少二個偵測目標50,各個填色框對應的速度(頻率)即為相應的偵測目標50的表徵資訊。
在一些實施例中,處理單元36可不對整個二維陣列A2執行快速傅立葉轉換,而僅針對頻域訊號SF的同一峰值處(於此為表示為填色框的兩個直行矩陣)執行快速傅立葉轉換,以減少運算量並節省計算時間。
由上述說明可以理解,透過對收到的啁啾訊號Cr進行訊號耦合、類比數位轉換與數位訊號處理,可以得到距離資訊及表徵資訊。然而,在室內環境下,雷達訊號容易受到偵測區域40之外的其他因素影響而造成資訊誤判。例如:訊號經牆壁或地板反射的靜態環境干擾、訊號受到其他移動物體(如電扇、窗簾)的擾動的動態環境干擾。因此,在取得數位訊號SD之後,且在進行前述距離傅立葉轉換與都卜勒傅立葉轉換等訊號處理之前,可先執行下述之數位訊號處理方法,以解決環境干擾造成的訊號不佳的問題。
合併參照圖7及圖8。圖7為依據一些實施例的訊號處理裝置60的方塊示意圖。圖8為依據一實施例的數位訊號處理方法的流程圖。訊號處理裝置60包含處理器61及儲存裝置62。儲存裝置62為電腦可讀取儲存媒體,供儲存供處理器61執行的程式63,以執行數位訊號處理方法。於一些實施例中,訊號處理裝置60為都卜勒雷達(如前述頻率調變連續波雷達10),處理器61為前述處理單元36。於一些實施例中,都卜勒雷達是連續波(Continuous Wave,CW)雷達或超寬頻(Ultra-wideband,UWB)雷達。於一些實施例中,訊號處理裝置60為邊緣裝置或是雲端伺服器,亦即頻率調變連續波雷達10獲得數位訊號SD之後,數位訊號SD將被傳送至邊緣裝置或雲端伺服器,由邊緣裝置或雲端伺服器進行數位訊號處理。
合併參照圖7至圖9,圖9為例示訊號疊加的示意圖。都卜勒雷達(頻率調變連續波雷達10)的接收端將所收到的對應各啁啾訊號C1~Cn的中頻訊號SI轉換為數位訊號P1~Pn。在流程S200中,疊加對應於啁啾訊號SC的數位訊號SP,而得到輸出訊號SE(如圖9所示,分別表示為E1~Ek,k為正整數)。以圖9所示之例來說明,數位訊號SP中每兩相鄰者進行疊加而成為輸出訊號SE。例如,數位訊號P1、P2疊加為輸出訊號E1,數位訊號P3、P4疊加為輸出訊號E2。透過疊加訊號可將環境造成的隨機雜訊平均化,而可提高訊雜比(Signal-to-Noise Ratio,SNR)。此外,訊號疊加後,可將多個一維陣列合併為一個一維陣列,將可大幅減少資料處理量,而可加快處理速度。在此,雖是以相鄰的兩個數位訊號SP進行疊加為例,然而本發明實施例非限於此,亦可以是將相鄰的兩個以上的數位訊號SP進行疊加。
在一些實施例中,相互疊加的數位訊號SP對應於同一訊框(Frame)中的複數啁啾訊號SC。在另一些實施例中,相互疊加的數位訊號SP對應於相鄰的複數訊框中的複數啁啾訊號SC。所述相鄰的複數訊框可以是兩個以上的訊框。
流程S300:依據疊加後的輸出訊號SE,可執行前述的距離傅立葉轉換和都卜勒傅立葉轉換,而可計算偵測目標50的距離資訊及表徵資訊(運動資訊及生理資訊)等資訊。
前述說明雖以發射單元20具有一個發射天線22,接收單元30具有一個接收天線32的情形為例。然而,在一些實施例中,發射單元20可具有多個發射天線22。相似地,在一些實施例中,接收單元30可具有多個接收天線32。
參照圖10,係為例示啁啾訊號SC的示意圖。在此,是以發射單元20具有二個發射天線22的情形為例。於此,啁啾訊號SC分別標示為Tx1及Tx2,分別表示由第一個發射天線22和第二個發射天線22發出的啁啾訊號SC。雷達訊號定義有複數訊框M(分別表示為M1~Mp,p為正整數)。每個訊框M中包括複數啁啾訊號SC。在此,是以各個發射天線22交替發送啁啾訊號Tx1、Tx2,共發送四個啁啾訊號SC為例。訊框M的期間可例如為20毫秒。在每一訊框M的期間內,處理器61依據各訊框M內的啁啾訊號SC執行前述數位訊號處理方法。也就是說,處理器61將位於同一訊框M中的啁啾訊號SC所對應的數位訊號SP疊加而得到輸出訊號SE,並執行流程S300。例如,在訊框M1期間內,依據訊框M1內二個啁啾訊號Tx1及二個啁啾訊號Tx2執行一次數位訊號處理方法;在訊框M2期間內,依據訊框M2內二個啁啾訊號Tx1及二個啁啾訊號Tx2執行一次數位訊號處理方法。以下說明不同的訊號疊加方式。
參照圖11,係為例示接收圖10訊號的示意圖。在此,是以具有四個接收天線Rx1、Rx2、Rx3、Rx4的接收單元30接收圖10所示的雷達訊號的情形為例。啁啾訊號Tx1、Tx2經過各個接收天線Rx1、Rx2、Rx3、Rx4接收,並由雷達接收器34轉換為數位訊號SP(於此僅以一個訊框M的四個數位訊號P1~P4來說明)。第一個數位訊號P1對應於第一個發射天線22發送的第一個啁啾訊號Tx1;第二個數位訊號P2對應於第二個發射天線22發送的第一個啁啾訊號Tx2;第三個數位訊號P3對應於第一個發射天線22發送的第二個啁啾訊號Tx1;第四個數位訊號P4對應於第二個發射天線22發送的第二個啁啾訊號Tx2。
參照圖12,係為另一例示訊號疊加的示意圖。處理器61可以疊加各接收天線Rx1~Rx4分別接收的來自同一發射天線22的啁啾訊號SC所對應的數位訊號SP。舉例來說,第一個數位訊號P1與第三個數位訊號P3均對應於來自第一個發射天線22的啁啾訊號Tx1,兩者疊加而成為疊加訊號SG(分別表示為G1~G4)。處理器61進一步將疊加訊號SG中對應各接收天線Rx1~Rx4的訊號(G1~G4)疊加,而得到輸出訊號SE。也就是說,處理器61疊加接收天線Rx1~Rx4所收到的來自同一發射天線22的啁啾訊號SC所對應的數位訊號SP而得到輸出訊號SE。在一些實施例中,處理器61是疊加來自其他發射天線22的啁啾訊號SC(如第二個發射天線22的啁啾訊號Tx2)所對應的數位訊號SP。
在一些實施例中,處理器61不將對應各接收天線Rx1~Rx4的疊加訊號SG進行疊加,而是選擇對應於其中一個接收天線Rx1~Rx4的疊加訊號G1、G2、G3或G4做為輸出訊號SE。也就是說,選擇此些接收天線Rx1~Rx4的其中一者所收到的來自同一發射天線22的啁啾訊號SC所對應的數位訊號SP進行疊加來獲得輸出訊號SE。
參照圖13,係為又一例示訊號疊加的示意圖。有別於如圖12所示僅疊加來自同一發射天線22的啁啾訊號SC所對應的數位訊號SP,處理器61也可以疊加來自不同的發射天線22(即至少兩個發射天線22)的啁啾訊號SC所對應的數位訊號SP。舉例來說,雖然數位訊號P1~P4來自不同的發射天線22,處理器61將各接收天線Rx1~Rx4分別接收的啁啾訊號SC所對應的數位訊號P1~P4疊加而成為疊加訊號SG(分別表示為G1~G4)。處理器61進一步將疊加訊號SG中對應各接收天線Rx1~Rx4的訊號(G1~G4)疊加,而得到輸出訊號SE。也就是說,處理器61疊加接收天線Rx1~Rx4所收到的來自不同發射天線22(即至少兩個發射天線22)的啁啾訊號SC所對應的數位訊號SP而得到輸出訊號SE。
參照圖14及圖15。圖14為例示在靜態環境中沒有進行訊號疊加的數位訊號SP進行距離傅立葉轉換後的示意圖。圖15為例示在靜態環境中數位訊號SP經過如圖13所示之訊號疊加並進行距離傅立葉轉換後的示意圖。橫軸為經過距離傅立葉轉換後的距離資訊,縱軸為訊號強度,豎軸為時間。可以看到經過訊號疊加後,左側的靜態環境干擾雜訊明顯減少了。
在一些實施例中,處理器61不將對應各接收天線Rx1~Rx4的疊加訊號SG進行疊加,而是選擇對應於其中一個接收天線Rx1~Rx4的疊加訊號G1、G2、G3或G4做為輸出訊號SE。也就是說,選擇此些接收天線Rx1~Rx4的其中一者所收到的來自不同發射天線22(即至少兩個發射天線22)的啁啾訊號SC所對應的數位訊號SP進行疊加來獲得輸出訊號SE。
在一些實施例中,雖前述是以疊加位於同一訊框M中的啁啾訊號SC來說明,然而處理器61也可以是將相鄰的兩個訊框M中的啁啾訊號SC所對應的數位訊號SP疊加而得到輸出訊號SE。在一些實施例中,處理器61也可以是將相鄰的至少三個訊框M中的啁啾訊號SC所對應的數位訊號SP疊加而得到輸出訊號SE。如前所述,相鄰至少兩訊框M中的數位訊號SP的疊加方式,可以是將來自同一發射天線22的啁啾訊號SC所對應的數位訊號SP進行疊加;或者是將來自不同發射天線22的啁啾訊號SC所對應的數位訊號SP進行疊加,於此不重複贅述。如前所述,可以選擇疊加後的訊號(疊加訊號SG)中的其中一個做為輸出訊號SE;也可以是將疊加訊號SG疊加而得到輸出訊號SE,於此不重複贅述。於此,由於是疊加至少兩相鄰訊框M中的啁啾訊號SC所對應的數位訊號SP,因而每至少兩個訊框M執行一次數位訊號處理方法。
參照圖16,係為依據另一實施例的數位訊號處理方法的流程圖。與圖8之差異在於,在流程S200之前,還先執行流程S100,對位待疊加的數位訊號SP。藉此,可修正此些數位訊號SP之間的相位延遲。在一些實施例中,若待疊加的數位訊號SP之間的相位延遲在容許範圍內,可不執行流程S100。例如,若處理器61是疊加位於同一訊框M的啁啾訊號SC,其相位延遲一般在容許範圍內,而可不執行流程S100。若處理器61是疊加相鄰的至少兩個訊框M中的啁啾訊號SC,且其相位延遲超出容許範圍,則在流程S200之前,還先執行流程S100。
參照圖17,係為依據一實施例的計算偵測目標的資訊的流程圖。流程S300除包括前述的距離傅立葉轉換(步驟S310)、都卜勒傅立葉轉換(步驟S370)和計算偵測目標50的資訊(步驟S390)之外,還包括步驟S330及S350。
合併參照圖17及圖18,圖18為例示統計資訊的產生示意圖。經過前述步驟S310而得到頻域訊號SF(橫軸為經過距離傅立葉轉換後的距離資訊,縱軸為訊號強度,豎軸為時間,如以二維表示即相當於圖6所示之二維陣列A2)之後,進入步驟S330。在步驟S330中,處理器61可先根據頻域訊號SF產生統計資訊SV。具體來說,係依據時間軸對頻域訊號SF進行統計資訊SV計算。透過統計資訊SV,可以分析訊號特性,以利於區別環境中的動態干擾與目標訊號源。統計資訊SV可例如為絕對平方和、絕對最大值、標準差(Standard Deviation)或變異數(Variance)。詳言之,處理器61是對於二維陣列A2的各個直行矩陣CL計算統計資訊SV,亦即對於每一個距離的時間序列訊號強度進行統計。例如,每一個直行矩陣CL都計算出一個標準差p。如圖18所示,係顯示按照此些標準差p所描繪出的分布曲線W,係能作為遮罩,以供後續步驟S350修正訊號。
在步驟S350中,處理器61依據統計資訊SV修正訊號。也就是說,依據統計資訊SV正規化頻域訊號SF。例如,對於每個距離的時間序列上的訊號乘上對應的統計資訊SV(如標準差p),即對於二維陣列A2的各個直行矩陣CL分別乘上對應各個直行矩陣CL對應的統計資訊SV(如標準差p)。如圖18所示,係顯示修正後的頻域訊號SF’。於後,便可依據修正後的頻域訊號SF’來計算更精確的偵測目標50的距離資訊、表徵資訊(如運動資訊、生理資訊)等資訊(步驟S390)。
合併參照圖19A~圖19E及圖20A~圖20E。圖19A~圖19E為例示不同環境干擾的頻域訊號SF的示意圖。圖20A~圖20E為例示不同環境干擾的修正後頻域訊號SF’的示意圖。圖19A及圖20A為風扇轉動但不擺頭的環境干擾。圖19B及圖20B為風扇轉動且擺頭的環境干擾。圖19C及圖20C為布簾靜止的環境干擾。圖19D及圖20D為布簾擺動的環境干擾。圖19E及圖20E為其他人移動的環境干擾。清楚可見,透過前述的訊號修正,確實可消除動態環境干擾。
合併參照圖21A及圖21B。圖21A為例示偵測目標50側躺的頻域訊號SF的示意圖。圖21B為例示偵測目標50側躺的修正後頻域訊號SF’的示意圖。透過前述的訊號修正,還可消除偵測目標50的不同姿勢(如側躺)所造成的床反射訊號干擾,因而對於不同的睡姿所造成的影響也有很好的抑制效果。
合併參照圖7及圖22。圖22為依據另一實施例的計算偵測目標的資訊的流程圖。儲存裝置62還儲存有機器學習模型64,係經過不同情形下所獲得的前述修正後訊號進行訓練。所述不同情形包括:各種環境條件(如前述風扇轉動、布簾擺動等)下,偵測區域40中存在偵測目標50時;以及各種環境條件(如前述風扇轉動、布簾擺動等)下,偵測區域40中不存在偵測目標50。圖22與圖17之差異在於,更包括步驟S360,係根據修正後的訊號,利用機器學習模型64判斷偵測區域40中是否存在偵測目標50。若存在偵測目標50,則繼續進行步驟S370、S390,以計算偵測目標50的距離資訊、表徵資訊(如運動資訊、生理資訊)等資訊;若否,則結束流程。透過機器學習技術,可學習偵測目標50存在與否的特徵,藉以判斷偵測目標50是否存在。藉此,除了可在不存在偵測目標50的情形下節省運算資源外,還可過濾誤判錯誤(例如避免發生不存在偵測目標50卻計算出表徵資訊的情形)。
在一些實施例中,當欲訓練機器學習模型64與利用機器學習模型64進行判斷時,可預先對於修正後的訊號進行特徵擷取處理,再將擷取到的特徵輸入至機器學習模型64。所述特徵擷取處理可以是統計直方圖、計算平均值、標準差、變異數、偏度(Skewness)、峰度(Kurtosis)等。以統計直方圖為例,可以將修正後的頻域訊號SF’的二維陣列A2數值歸一化,並依照歸一化後的數值統計各個訊號強度區間的元素數量。合併參照圖23及圖24所示,係分別例示不存在偵測目標50及存在偵測目標50的直方統計圖,以分為10個訊號強度區間為例,可以看到兩種情形的直方圖具有不同的分佈。
參照圖25,係為依據一些實施例的表徵資訊偵測方法的流程圖。所述表徵資訊偵測方法可由前述訊號處理裝置60的處理器61執行。在步驟S410中,接收對應於都卜勒雷達的複數數位偵測訊號(即前述數位訊號SP)。在步驟S420中,對數位偵測訊號進行頻域分析,以得到複數頻域偵測訊號(即前述頻域訊號SF)。在步驟S430中,依據頻域偵測訊號產生複數統計資訊(即前述統計資訊SV)。在步驟S440中,依據統計資訊SV修正頻域偵測訊號。具體地,步驟S440係依據統計資訊SV正規化頻域偵測訊號。統計資訊SV係頻域偵測訊號對應於一統計週期下的絕對平方和、絕對最大值、標準差、或是變異數。在步驟S450中,根據修正後的頻域偵測訊號(即前述修正後頻域訊號SF’)及機器學習模型64判斷偵測區域40中是否存在偵測目標50。在步驟460中,響應於存在偵測目標50,依據修正後的頻域偵測訊號SF’計算偵測目標50的表徵資訊。各步驟之相關說明已於前詳細說明,於此恕不重複贅述。
在一些實施例中,機器學習模型64係儲存在邊緣裝置或是雲端伺服器之中,頻率調變連續波雷達10將前述擷取到的特徵透過其傳輸模組傳送給該邊緣裝置或是該雲端伺服器,以進行後續之機器模型訓練或是偵測判斷。或是,頻率調變連續波雷達10將前述修正後訊號透過其傳輸模組傳送給該邊緣裝置或是該雲端伺服器,以進行後續之數位訊號處理、機器模型訓練、及/或是偵測判斷。
綜上所述,依據一些實施例的頻率調變連續波雷達與數位訊號處理方法,能解決靜態與動態環境干擾造成的訊號不佳的問題,並可減少資料處理量而可加快處理速度。依據一些實施例的頻率調變連續波雷達與數位訊號處理方法,能識別偵測目標50存在與否,以增進處理效率並過濾誤判錯誤。
10:頻率調變連續波雷達
20:發射單元
22:發射天線
24:雷達發射器
30:接收單元
32:接收天線
34:雷達接收器
36:處理單元
40:偵測區域
50:偵測目標
60:訊號處理裝置
61:處理器
62:儲存裝置
63:程式
64:機器學習模型
M,M1,M2,Mp:訊框
A1,A2,A3:二維陣列
B:脈衝帶寬
CL:直行矩陣
Ct,Cr:啁啾訊號
C1,C2,C3,Cn:啁啾訊號
D1,D2,D3,D4,Dn:數位訊號
d:傳輸距離
E1,E2,Ek:輸出訊號
F1,F2,Fn:頻域訊號
G1,G2,G3,G4:疊加訊號
Rx1,Rx2,Rx3,Rx4:接收天線
P1,P2,P3,P4,Pn:數位訊號
p:標準差
Q1,Q2,Qm:相位頻域訊號
S:斜率
SC:啁啾訊號
SD:數位訊號
SE:輸出訊號
SG:疊加訊號
SI:中頻訊號
SP:數位訊號
SF,SF’:頻域訊號
SQ:相位頻域訊號
SV:統計資訊
W:分布曲線
S100:對位待疊加的數位訊號
S200:疊加啁啾訊號所對應的數位訊號而得到輸出訊號
S300:依據輸出訊號計算偵測目標的資訊
S310:距離傅立葉轉換
S330:產生統計資訊
S350:依據統計資訊修正訊號
S360:根據機器學習模型判斷偵測區域中是否存在偵測目標
S370:都卜勒傅立葉轉換
S390:計算偵測目標的資訊
S410:接收對應於都卜勒雷達的複數數位偵測訊號
S420:對數位偵測訊號進行頻域分析,以得到複數頻域偵測訊號
S430:依據頻域偵測訊號產生複數統計資訊
S440:依據統計資訊修正頻域偵測訊號
S450:根據修正後的頻域偵測訊號及機器學習模型判斷偵測區域中是否存在偵測目標
S460:響應於存在偵測目標,依據修正後的頻域偵測訊號計算偵測目標的表徵資訊
Tx1,Tx2:啁啾訊號
Tc:持續時間
τ:延遲時間
[圖1]為依據一些實施例的頻率調變連續波雷達的使用狀態示意圖;
[圖2]為依據一些實施例的頻率調變連續波雷達的方塊示意圖;
[圖3]為例示雷達訊號的示意圖;
[圖4]為依據一些實施例的頻率調變連續波雷達的細部方塊示意圖;
[圖5]為例示發送與接收的雷達訊號的示意圖;
[圖6]為依據一些實施例的訊號處理示意圖;
[圖7]為依據一些實施例的訊號處理裝置的方塊示意圖;
[圖8]為依據一實施例的數位訊號處理方法的流程圖;
[圖9]為例示訊號疊加的示意圖;
[圖10]為例示啁啾訊號的示意圖;
[圖11]為例示接收圖10訊號的示意圖;
[圖12]為另一例示訊號疊加的示意圖;
[圖13]為又一例示訊號疊加的示意圖;
[圖14]為例示在靜態環境中沒有進行訊號疊加的數位訊號進行距離傅立葉轉換後的示意圖;
[圖15]為例示在靜態環境中數位訊號經過如圖13所示之訊號疊加並進行距離傅立葉轉換後的示意圖;
[圖16]為依據另一實施例的數位訊號處理方法的流程圖;
[圖17]為依據一實施例的計算偵測目標的資訊的流程圖;
[圖18]為例示統計資訊的產生示意圖;
[圖19A~圖19E]為例示不同環境干擾的頻域訊號的示意圖;
[圖20A~圖20E]為例示不同環境干擾的修正後頻域訊號的示意圖;
[圖21A]為例示偵測目標側躺的頻域訊號的示意圖;
[圖21B]為例示偵測目標側躺的修正後頻域訊號的示意圖;
[圖22]為依據另一實施例的計算偵測目標的資訊的流程圖;
[圖23]為例示不存在偵測目標的直方統計圖;
[圖24]為例示存在偵測目標的直方統計圖;及
[圖25]為依據一些實施例的表徵資訊偵測方法的流程圖。
S200:疊加啁啾訊號所對應的數位訊號而得到輸出訊號
S300:依據輸出訊號計算偵測目標的資訊
Claims (22)
- 一種頻率調變連續波雷達,包括:一發射單元,發送複數啁啾訊號;及一接收單元,接收反射的該些啁啾訊號,產生對應該些啁啾訊號的複數數位訊號,對每兩相鄰之該些數位訊號進行疊加而得到一輸出訊號,並依據該輸出訊號計算一偵測目標的表徵資訊。
- 如請求項1所述之頻率調變連續波雷達,其中該接收單元將位於同一訊框中的該些啁啾訊號所對應的該些數位訊號疊加而得到該輸出訊號。
- 如請求項1所述之頻率調變連續波雷達,其中該接收單元將相鄰的複數個訊框中的該些啁啾訊號所對應的該些數位訊號疊加而得到該輸出訊號。
- 如請求項3所述之頻率調變連續波雷達,其中該接收單元先對位該些數位訊號再對其疊加。
- 如請求項1所述之頻率調變連續波雷達,其中該發射單元包括至少一發射天線,該接收單元疊加來自同一該發射天線的該些啁啾訊號所對應的該些數位訊號。
- 如請求項5所述之頻率調變連續波雷達,其中該接收單元包括複數接收天線,該接收單元將其中一個該接收天線所收到的來自同一該發射天線的該些啁啾訊號所對應的該些數位訊號進行疊加而得到該輸出訊號。
- 如請求項5所述之頻率調變連續波雷達,其中該接收單元包括複數接收天線,該接收單元疊加該些接收天線所收到的來自同一該發射天線的該些啁啾訊號所對應的該些數位訊號而得到該輸出訊號。
- 如請求項1所述之頻率調變連續波雷達,其中該發射單元包括複數發射天線,該接收單元疊加來自該些發射天線中的至少二者的該些啁啾訊號所對應的該些數位訊號。
- 如請求項8所述之頻率調變連續波雷達,其中該接收單元包括複數接收天線,該接收單元將其中一個該接收天線所收到的來自該些發射天線中的至少二者的該些啁啾訊號所對應的該些數位訊號進行疊加而得到該輸出訊號。
- 如請求項8所述之頻率調變連續波雷達,其中該接收單元包括複數接收天線,該接收單元疊加該些接收天線所收到的來自該些發射天線中的至少二者的該些啁啾訊號所對應的該些數位訊號而得到該輸出訊號。
- 如請求項1至10中任一項所述之頻率調變連續波雷達,其中該接收單元還根據該輸出訊號產生一統計資訊,並以該統計資訊為遮罩修正該輸出訊號。
- 如請求項11所述之頻率調變連續波雷達,其中該統計資訊為絕對平方和、絕對最大值、標準差或變異數。
- 如請求項11所述之頻率調變連續波雷達,其中該接收單元根據該輸出訊號及一機器學習模型判斷一偵測區域中是否存在一偵測目標。
- 如請求項13所述之頻率調變連續波雷達,其中該接收單元於判斷到該偵測目標存在時,依據該輸出訊號計算該偵測目標的表徵資訊。
- 一種數位訊號處理方法,由一訊號處理裝置中之一處理器執行,包括:疊加對應於一都卜勒雷達的一接收端收到的複數啁啾訊號的複數數位訊號而得到一輸出訊號,該疊加係對每兩相鄰之該些數位訊號進行疊加;及依據該輸出訊號計算一偵測目標的表徵資訊。
- 如請求項15所述之數位訊號處理方法,其中所疊加的該些啁啾訊號所對應的該些數位訊號是位於同一訊框中或位於相鄰的複數個訊框中。
- 如請求項15所述之數位訊號處理方法,其中在該疊加對應於該都卜勒雷達的該接收端收到的該些啁啾訊號的該些數位訊號而得到該輸出訊號的步驟之前更包括:對位待疊加的該些數位訊號。
- 如請求項15所述之數位訊號處理方法,更包括:根據該輸出訊號產生一統計資訊;及以該統計資訊為遮罩修正該輸出訊號。
- 如請求項15或18所述之數位訊號處理方法,更包括:根據該輸出訊號及一機器學習模型判斷一偵測區域中是否存在一偵測目標。
- 如請求項19所述之數位訊號處理方法,其中該依據該輸出訊號計算該偵測目標的表徵資訊的步驟是於判斷到該偵測目標存在時執行。
- 一種表徵資訊偵測方法,包括:接收對應於一都卜勒雷達的複數數位偵測訊號;對該些數位偵測訊號進行頻域分析,以得到複數頻域偵測訊號;依據該些頻域偵測訊號產生複數統計資訊;以該些統計資訊為遮罩修正該些頻域偵測訊號;根據修正後的該些頻域偵測訊號及一機器學習模型判斷一偵測區域中是否存在一偵測目標;及響應於存在該偵測目標,依據修正後的該些頻域偵測訊號計算該偵測目標的表徵資訊。
- 如請求項21所述之表徵資訊偵測方法,其中依據該些統計資訊修正該些頻域偵測訊號更包括:依據該些統計資訊正規化該些頻域偵測訊號,其中該些統計資訊係該些頻域偵測訊號對應於一統計週期下的絕對平方和、絕對最大值、標準差、或是變異數。
Priority Applications (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
TW110109640A TWI762243B (zh) | 2021-03-17 | 2021-03-17 | 頻率調變連續波雷達、數位訊號處理方法與表徵資訊偵測方法 |
CN202110458353.6A CN115113143A (zh) | 2021-03-17 | 2021-04-27 | 频率调制连续波雷达、数字信号处理方法与表征信息检测方法 |
US17/349,264 US12135388B2 (en) | 2021-03-17 | 2021-06-16 | FMCW radar, method for processing digital signals, and characterization information detection method |
EP21190065.9A EP4060369A1 (en) | 2021-03-17 | 2021-08-06 | Fmcw radar, method for processing digital signals, and characterization information detection method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
TW110109640A TWI762243B (zh) | 2021-03-17 | 2021-03-17 | 頻率調變連續波雷達、數位訊號處理方法與表徵資訊偵測方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
TWI762243B true TWI762243B (zh) | 2022-04-21 |
TW202238167A TW202238167A (zh) | 2022-10-01 |
Family
ID=77249721
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
TW110109640A TWI762243B (zh) | 2021-03-17 | 2021-03-17 | 頻率調變連續波雷達、數位訊號處理方法與表徵資訊偵測方法 |
Country Status (4)
Country | Link |
---|---|
US (1) | US12135388B2 (zh) |
EP (1) | EP4060369A1 (zh) |
CN (1) | CN115113143A (zh) |
TW (1) | TWI762243B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1277065B1 (en) * | 2000-03-09 | 2009-12-02 | Tele-IP Limited | Acoustic sounding |
TW201643463A (zh) * | 2015-04-07 | 2016-12-16 | 光電波有限公司 | 小型光達系統 |
US20200011968A1 (en) * | 2017-03-03 | 2020-01-09 | Iee International Electronics & Engineering S.A. | Method and system for obtaining an adaptive angle-doppler ambiguity function in mimo radars |
CN111295596A (zh) * | 2019-02-28 | 2020-06-16 | 深圳市大疆创新科技有限公司 | 毫米波雷达的测角方法、设备及存储介质 |
TW202102979A (zh) * | 2019-05-22 | 2021-01-16 | 美商塔切爾實驗室公司 | 毫米波陣列 |
Family Cites Families (37)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3877010A (en) * | 1967-10-26 | 1975-04-08 | Hughes Aircraft Co | Adaptable moving target processor |
GB2072986B (en) * | 1980-03-26 | 1984-06-13 | Standard Telephones Cables Ltd | Phased array signal processing |
US7477758B2 (en) * | 1992-05-05 | 2009-01-13 | Automotive Technologies International, Inc. | System and method for detecting objects in vehicular compartments |
US6809681B1 (en) * | 1992-08-25 | 2004-10-26 | Raytheon Company | Random-modulation radar signal-induced interference cancellation method and apparatus |
US5390133A (en) * | 1992-09-30 | 1995-02-14 | Martin Marietta Corporation | Image processor for target detection and tracking |
JP2005520160A (ja) * | 2002-03-13 | 2005-07-07 | レイセオン・カナダ・リミテッド | レーダのスペクトル発生のためのシステムおよび方法 |
US6603424B1 (en) * | 2002-07-31 | 2003-08-05 | The Boeing Company | System, method and computer program product for reducing errors in synthetic aperture radar signals |
ATE519126T1 (de) * | 2004-06-24 | 2011-08-15 | Bae Systems Plc | Verbesserungen in bezug auf die geschwindigkeitsextraktion |
EP1672382A1 (de) * | 2004-12-18 | 2006-06-21 | Leica Geosystems AG | Einkanal-Heterodyn -Distanzmessverfahren |
IL170726A (en) * | 2005-09-07 | 2011-12-29 | Camero Tech Ltd | Signal acquisition system and method for ultra-wideband (uwb) radar |
US20100152600A1 (en) * | 2008-04-03 | 2010-06-17 | Kai Sensors, Inc. | Non-contact physiologic motion sensors and methods for use |
US8649446B2 (en) * | 2010-10-29 | 2014-02-11 | Texas Instruments Incorporated | System and method for channel classification |
US8742935B2 (en) * | 2011-06-30 | 2014-06-03 | General Electric Company | Radar based systems and methods for detecting a fallen person |
US8717230B1 (en) * | 2012-02-15 | 2014-05-06 | Lockheed Martin Corporation | Digital clutter spatial nulling |
US20200064444A1 (en) * | 2015-07-17 | 2020-02-27 | Origin Wireless, Inc. | Method, apparatus, and system for human identification based on human radio biometric information |
LU92331B1 (en) | 2013-12-10 | 2015-06-11 | Iee Sarl | Radar sensor with frequency dependent beam steering |
JP6940414B2 (ja) * | 2015-04-20 | 2021-09-29 | レスメッド センサー テクノロジーズ リミテッド | 特性信号から人間の検出及び識別 |
US11125866B2 (en) * | 2015-06-04 | 2021-09-21 | Chikayoshi Sumi | Measurement and imaging instruments and beamforming method |
US10037671B2 (en) * | 2015-06-29 | 2018-07-31 | Echocare Technologies Ltd. | Human respiration feature extraction in personal emergency response systems and methods |
US10078131B2 (en) | 2015-09-15 | 2018-09-18 | Texas Instruments Incorporated | Method and apparatus for FMCW radar processing |
US10514770B2 (en) * | 2016-06-17 | 2019-12-24 | Texas Instruments Incorporated | Hidden Markov model-based gesture recognition with FMCW radar |
US11412937B2 (en) | 2017-03-29 | 2022-08-16 | Texas Instruments Incorporated | Multi-person vital signs monitoring using millimeter wave (mm-wave) signals |
EP3425419B1 (en) | 2017-07-05 | 2024-01-03 | Stichting IMEC Nederland | A method and a system for localization and monitoring of living beings |
US11204647B2 (en) * | 2017-09-19 | 2021-12-21 | Texas Instruments Incorporated | System and method for radar gesture recognition |
DE102018123383A1 (de) * | 2017-10-13 | 2019-04-18 | Infineon Technologies Ag | Radarerfassung mit Störungsunterdrückung |
US11119186B2 (en) | 2017-12-07 | 2021-09-14 | Texas Instruments Incorporated | Radar processing chain for frequency-modulated continuous wave radar systems |
US10962637B2 (en) | 2018-06-25 | 2021-03-30 | Texas Instruments Incorporated | Radar data processing using neural network classifier and confidence metrics |
DE102018121987A1 (de) | 2018-09-10 | 2020-03-12 | Infineon Technologies Ag | Frequenzmoduliertes Dauerstrich-Radarsystem |
US11125869B2 (en) | 2018-10-16 | 2021-09-21 | Infineon Technologies Ag | Estimating angle of human target using mmWave radar |
CN111580086B (zh) * | 2019-02-19 | 2023-08-25 | 富士通株式会社 | 生命体检测方法、检测装置和电子设备 |
TWI712392B (zh) | 2019-02-27 | 2020-12-11 | 和碩聯合科技股份有限公司 | 多目標生命徵象偵測系統及方法 |
WO2020226638A1 (en) | 2019-05-08 | 2020-11-12 | Google Llc | Sleep tracking and vital sign monitoring using low power radio waves |
US11422230B2 (en) | 2019-09-06 | 2022-08-23 | Infineon Technologies Ag | System and method for receiving a radar signal |
CN111142102B (zh) | 2019-12-26 | 2022-07-19 | 华为技术有限公司 | 一种呼吸数据计算方法以及相关设备 |
US11385344B2 (en) * | 2020-03-20 | 2022-07-12 | Aptiv Technologies Limited | Frequency-modulated continuous-wave (FMCW) radar-based detection of living objects |
US11415670B2 (en) * | 2020-03-20 | 2022-08-16 | Aptiv Technologies Limited | Object classification using low-level radar data |
US11774553B2 (en) * | 2020-06-18 | 2023-10-03 | Infineon Technologies Ag | Parametric CNN for radar processing |
-
2021
- 2021-03-17 TW TW110109640A patent/TWI762243B/zh active
- 2021-04-27 CN CN202110458353.6A patent/CN115113143A/zh active Pending
- 2021-06-16 US US17/349,264 patent/US12135388B2/en active Active
- 2021-08-06 EP EP21190065.9A patent/EP4060369A1/en active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1277065B1 (en) * | 2000-03-09 | 2009-12-02 | Tele-IP Limited | Acoustic sounding |
TW201643463A (zh) * | 2015-04-07 | 2016-12-16 | 光電波有限公司 | 小型光達系統 |
US20200011968A1 (en) * | 2017-03-03 | 2020-01-09 | Iee International Electronics & Engineering S.A. | Method and system for obtaining an adaptive angle-doppler ambiguity function in mimo radars |
CN111295596A (zh) * | 2019-02-28 | 2020-06-16 | 深圳市大疆创新科技有限公司 | 毫米波雷达的测角方法、设备及存储介质 |
TW202102979A (zh) * | 2019-05-22 | 2021-01-16 | 美商塔切爾實驗室公司 | 毫米波陣列 |
Also Published As
Publication number | Publication date |
---|---|
CN115113143A (zh) | 2022-09-27 |
TW202238167A (zh) | 2022-10-01 |
US20220299598A1 (en) | 2022-09-22 |
EP4060369A1 (en) | 2022-09-21 |
US12135388B2 (en) | 2024-11-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111352102B (zh) | 一种基于调频连续波雷达的多目标个数检测方法及装置 | |
US10795012B2 (en) | System and method for human behavior modelling and power control using a millimeter-wave radar sensor | |
EP3869222A1 (en) | Radar vital signal tracking using a kalman filter | |
JP6270901B2 (ja) | Fmcwレーダ装置 | |
US9897685B2 (en) | Method for detecting interference in a received signal of a radar sensor, driver assistance device and motor vehicle | |
CN108652601A (zh) | 一种基于调频连续波毫米波雷达的睡眠监控方法、装置和雷达系统 | |
EP3342330B1 (en) | Phasor approach to signal to noise ratio measurement evaluation of physiological signals | |
CN109932700B (zh) | 一种多普勒速度的解模糊方法及装置 | |
EP4187275A1 (en) | Cfar phased array pre-processing using noncoherent and coherent integration in automotive radar systems | |
JP6995258B2 (ja) | レーダ信号処理器、レーダセンサシステム及び信号処理方法 | |
WO2013150992A1 (en) | Methods and apparatus for ultrasound imaging | |
CN111505631A (zh) | 基于lfmcw雷达的心率估计算法 | |
JP4729765B2 (ja) | パルスドプラ計測装置、その方法及びそのプログラム | |
CN110320509A (zh) | 一种毫米波雷达目标检测方法及系统 | |
CN117310670B (zh) | 基于超声波雷达的测量方法、装置、车载终端和存储介质 | |
Chen et al. | Ultrawideband synthetic aperture radar for respiratory motion detection | |
TWI762243B (zh) | 頻率調變連續波雷達、數位訊號處理方法與表徵資訊偵測方法 | |
CN115561728A (zh) | 一种fmcw雷达的单帧速度解模糊方法和装置 | |
CN107783099A (zh) | 基于组合波形的旋翼无人机短距离防撞系统信号处理系统及方法 | |
US11876582B2 (en) | Hybrid beam steering radar | |
JP7140568B2 (ja) | 到来方向推定装置及び到来方向推定方法 | |
CN103154767B (zh) | 距离测量装置 | |
TWI586987B (zh) | Signal processing device for continuous wave radar sensing system | |
KR102381262B1 (ko) | 레이더 시스템 및 그에서 수행되는 생체 신호 탐지방법 | |
JP2013113723A (ja) | レーダ装置 |