KR101628154B1 - Multiple target tracking method using received signal strengths - Google Patents

Multiple target tracking method using received signal strengths Download PDF

Info

Publication number
KR101628154B1
KR101628154B1 KR1020150030943A KR20150030943A KR101628154B1 KR 101628154 B1 KR101628154 B1 KR 101628154B1 KR 1020150030943 A KR1020150030943 A KR 1020150030943A KR 20150030943 A KR20150030943 A KR 20150030943A KR 101628154 B1 KR101628154 B1 KR 101628154B1
Authority
KR
South Korea
Prior art keywords
target
data
received signal
signal strength
probability
Prior art date
Application number
KR1020150030943A
Other languages
Korean (ko)
Inventor
김수진
정영헌
양동원
김성준
윤주홍
Original Assignee
국방과학연구소
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 국방과학연구소 filed Critical 국방과학연구소
Priority to KR1020150030943A priority Critical patent/KR101628154B1/en
Application granted granted Critical
Publication of KR101628154B1 publication Critical patent/KR101628154B1/en

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems 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/66Radar-tracking systems; Analogous systems
    • G01S13/72Radar-tracking systems; Analogous systems for two-dimensional tracking, e.g. combination of angle and range tracking, track-while-scan radar
    • G01S13/723Radar-tracking systems; Analogous systems for two-dimensional tracking, e.g. combination of angle and range tracking, track-while-scan radar by using numerical data
    • G01S13/726Multiple target tracking
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems 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/02Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
    • G01S13/50Systems of measurement based on relative movement of target
    • G01S13/52Discriminating between fixed and moving objects or between objects moving at different speeds
    • G01S13/522Discriminating between fixed and moving objects or between objects moving at different speeds using transmissions of interrupted pulse modulated waves
    • G01S13/524Discriminating 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/53Discriminating 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
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/28Details of pulse systems

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

The present invention relates to a multiple target tracking method using a received signal strength indicator (RSSI), and more specifically, to a multiple target tracking method using an RSSI, which generates RSSI data from a plurality of tracking signals detected by a radar, extracts only a tracking signal present within a certain effective distance among tracking signals in which the RSSI data is generated, generates the data related probability of the extracted tracking signal by using the RSSI data, and obtains a state vector computation value and an error covariance value for a target by using a tracking filter to which the data related probability is applied, thus updating the track of the target. According to the present invention, different RSSI data by target is applied to Joint Probability Data-Association (JPDA) or Cheap Joint Probability Data-Association (CJPDA), thus improving the performance of a data related technique.

Description

수신 신호 세기를 이용한 다중 표적 추적 방법{MULTIPLE TARGET TRACKING METHOD USING RECEIVED SIGNAL STRENGTHS}[0001] MULTIPLE TARGET TRACKING METHOD USING RECEIVED SIGNAL STRENGTHS [0002]

본 발명은 수신 신호 세기를 이용한 다중 표적 추적 방법에 관한 것으로, 더욱 상세하게는, 레이더가 탐지한 다수의 표적 신호로부터 수신 신호 세기 데이터를 생성하고 상기 수신 신호 세기 데이터가 생성된 표적 신호 중 소정의 유효 거리 이내에 존재하는 표적 신호만을 추출한 후, 상기 수신 신호 세기 데이터를 이용하여 상기 추출된 표적 신호의 데이터 연관 확률을 생성하고, 상기 데이터 연관 확률이 적용된 추적 필터를 이용하여 상기 표적에 대한 상태 벡터 연산치 및 오차공분산을 구함으로써 상기 표적의 항적을 갱신하는 수신 신호 세기를 이용한 다중 표적 추적 방법에 관한 것이다.
The present invention relates to a multi-target tracking method using received signal strengths, and more particularly, to a multi-target tracking method using received signal strengths to generate received signal strength data from a plurality of target signals detected by a radar, The method includes generating a data association probability of the extracted target signal using the received signal strength data after extracting only a target signal existing within an effective distance, performing a state vector calculation on the target using the tracking filter to which the data association probability is applied And a multi-target tracking method using the received signal strength to update the trajectory of the target by obtaining a covariance and an error covariance.

다중 표적 추적 기술은, 레이더에 의해 탐지된 다수의 표적을 동시에 추적하는 기술이다. 상기 다중 표적 추적 기술에서는, 다수의 이동하는 표적에 대한 항적 오할당률을 줄이고, 다수의 이동하는 표적에 대한 항적 유지율을 높임으로써, 실시간으로 다수의 표적을 정확하게 추적을 가능하게 하기 위하여 고정밀 추적 기술이 필요하다.Multi-target tracking technology is a technique that simultaneously tracks multiple targets detected by a radar. In the multi-target tracking technique, high-precision tracking technology is used to reduce the rate of trafic perturbation for a large number of moving targets and increase the retention rate of a large number of moving targets, thereby enabling accurate tracking of a large number of targets in real time need.

도 1 은 종래 기술에 의한 다중 표적 추적 방법을 보인 흐름도이다.1 is a flowchart illustrating a conventional multi-target tracking method.

도 1 을 참조하면, 레이더에 의해 탐지된 각 표적들에 대한 탐지된 결과는 각 표적의 거리, 속도 및 가속도 등의 정보를 통하여 연산될 수 있다. 그리고, 상기 각 표적들이 이동하는 항적들은, 상기 각 표적들에 대한 탐지 결과를 기반으로, 각 표적들의 상태 벡터 측정치와 항적의 관계를 정의해주는 데이터 연관 기법 및 상기 각 표적들의 항적을 시간에 따라 개선해 나가는 필터 기법을 통하여 수행될 수 있다. Referring to FIG. 1, the detected results for each of the targets detected by the radar can be calculated through information such as distance, speed, and acceleration of each target. The wakes in which each of the targets are moving include a data associating technique for defining a relationship between a state vector measurement of each target and a wake based on the detection result for each of the targets and improving the wake of each of the targets over time Outgoing filter technique.

따라서, 다중 표적 추적 기술의 성능은, 상기 데이터 연관 기법에 따라 다를 수 있다. 상기 데이터 연관 기법에는, 빠른 연산이 가능한 최근접 연관 기법(NN : Nearest Neighbor), 복잡도는 높지만 정확한 추적 성능을 보이는 결합 확률 데이터 연관 기법(JPDA : Joint Probability Data-Association) 및 상기 결합 확률 데이터 연관 기법(JPDA)의 복잡도를 줄이기 위해 부분적 이벤트만을 고려한 단순 결합 확률 데이터 연관 기법(CJPDA : Cheap Joint Probability Data-Association)이 사용될 수 있다. Thus, the performance of the multiple target tracking technique may vary depending on the data association technique. The data association technique includes a nearest neighbor technique (NN) capable of fast operation, a joint probability data association technique (JPDA) having a high complexity but showing accurate tracking performance, and the joint probability data association technique (JJPDA: Cheap Joint Probability Data-Association (CJPDA)) that considers only partial events can be used to reduce the complexity of JPDA.

한편, 다중 표적 추적 기술의 성능을 향상시키기 위하여, 각 표적의 거리 정보뿐만이 아니라, 각 표적으로부터 획득된 표적의 신호 세기 정보를 혼용하여, 데이터 연관 기법에 적용하는 기술이 개발되어 있다. 그러나, 표적의 종류 및 표적의 자세 특성에 따라서 변화하는 표적의 신호 세기 정보가 정확하게 반영되지 않는 경우, 데이터 연관 기법에서 표적의 신호 세기 정보를 사용하는 기술은 다중 표적 추적 기술의 성능을 악화시킬 수 있다. 즉, 표적의 수신 신호 세기 데이터를 다중 표적 추적 방법에 적용하는 경우 표적별 평균신호대 잡음비가 같다는 가정 하에 표적의 항적을 갱신하여 왔으므로, 표적의 신호 세기 정보가 제대로 반영되지 못하여 다중 표적 추적 장치의 성능이 악화될 수 있는 것이다. On the other hand, in order to improve the performance of the multi-target tracking technique, a technique of applying the signal intensity information of the target obtained from each target, not only the distance information of each target, to the data association technique has been developed. However, if the signal strength information of a target that changes according to the type of target and posture characteristics of the target is not accurately reflected, a technique using the signal strength information of the target in the data association technique may deteriorate the performance of the multi-target tracking technique have. That is, when the received signal strength data of the target is applied to the multi-target tracking method, since the target's wake-up rate has been updated under the assumption that the average signal-to-noise ratio per target is the same, The performance may deteriorate.

따라서, 다중 표적 추적 기술의 성능을 향상시키기 위하여 표적의 신호 세기 정보를 이용하는 경우, 표적에 따른 적응형 신호 세기 정보를 적용해야하는 필요가 있다.
Therefore, when the signal strength information of the target is used to improve the performance of the multi-target tracking technique, adaptive signal strength information according to the target needs to be applied.

KR 1402206 B1KR 1402206 B1 KR 1001559 B1KR 1001559 B1

1. Donald Lerro and Yaakov Bar-Shalom, "Automated Tracking with Target Amplitude Information," 2875-28801. Donald Lerro and Yaakov Bar-Shalom, "Automated Tracking with Target Amplitude Information," 2875-2880

본 발명은 상기와 같은 문제점을 해소하기 위하여 제안된 것으로, 각 표적별로 상이한 수신 신호 세기 데이터를 결합 확률 데이터 연관 기법 또는 단순 결합 확률 데이터 연관 기법에 적용하여 데이터 연관 기법의 성능을 향상시킴으로써 보다 정확한 수신 신호 세기를 이용한 다중 표적 추적 방법을 제공하는데 그 목적이 있다.
The present invention has been proposed in order to solve the above problems. By improving the performance of a data association technique by applying different received signal strength data for each target to a joint probability data association technique or a simple joint probability data association technique, And to provide a multi-target tracking method using signal strength.

상기 목적을 달성하기 위하여, 본 발명에 의한 수신 신호 세기를 이용한 다중 표적 추적 방법의 일 실시 예는, 레이더가 탐지한 표적 신호로부터 표적의 위치 데이터를 추출하는 표적 위치 추출단계; 추적 필터를 이용하여 상기 추출된 표적의 위치 데이터로부터 상기 표적의 항적 및 상기 표적에 대한 상태 벡터 예측치를 생성하는 항적 생성단계; 상기 추출된 표적의 위치 데이터로부터 표적의 거리, 표적의 속도 및 표적으로부터의 수신 신호 세기 데이터를 포함하는 표적 데이터를 생성하는 표적 데이터 생성단계; 상기 표적 데이터가 생성된 표적 신호들 중 상기 수신 신호 세기 데이터를 이용하여 생성된 표적의 최종 항적 위치로부터 소정의 유효 거리 이내에 존재하는 표적 신호를 추출하는 게이트닝 단계; 상기 표적으로부터의 수신 신호 세기 데이터를 이용하여 상기 생성된 표적의 최종 항적 위치와 추출된 표적 신호의 위치 사이의 데이터 연관 확률을 생성하는 데이터 연관단계; 및 상기 표적의 데이터 연관 확률을 추적 필터에 입력하여 상기 생성된 표적의 항적을 갱신하는 항적 갱신단계; 를 포함할 수 있다.In order to achieve the above object, a multi-target tracking method using the received signal strength according to the present invention includes: a target position extracting step of extracting target position data from a target signal detected by a radar; Generating a trajectory of the target and a state vector prediction for the target from the extracted positional data of the target using a tracking filter; A target data generation step of generating target data including the distance of the target, the speed of the target, and the received signal strength data from the target, from the extracted positional data of the target; A gating step of extracting a target signal existing within a predetermined effective distance from a final navigation position of a target generated using the received signal strength data of the target signals from which the target data is generated; A data associating step of generating a data associating probability between the final navigation position of the generated target and the position of the extracted target signal using the received signal strength data from the target; And a navigation updating step of inputting a data association probability of the target into a tracking filter to update the generated map of the generated target; . ≪ / RTI >

이때, 상기 표적 데이터 생성단계에서, 상기 표적으로부터의 수신 신호 세기 데이터는 표적의 종류, 표적의 자세에 따라 달라질 수 있다.At this time, in the target data generation step, the received signal strength data from the target may vary depending on the type of target and the posture of the target.

또한, 상기 데이터 연관단계는, 상기 생성된 표적의 최종 항적 위치와 상기 추출된 표적 신호의 위치로부터 연산된 상태 벡터 측정치(xj(k+1))와 상기 생성된 표적에 대한 상태 벡터 예측치(xt(k+1)) 차이인 이노베이션 값(vj,t(k+1))을 연산하고 상기 이노베이션 값(vj,t(k+1))으로부터 유효 행렬(Ω)을 생성하는 유효 행렬 생성단계; 및 상기 유효 행렬(Ω)로부터 발생 가능 행렬(θ)을 도출하고 상기 발생 가능 행렬(θ)을 이용하여 상기 연산된 상태 벡터 측정치(xj(k+1))가 상기 표적에 연관될 데이터 연관 확률(βj,t)을 생성하는 연관 확률 생성단계; 를 포함할 수 있다.The data associating step may further include a step of calculating a state vector measurement value (x j (k + 1)) calculated from a final trajectory position of the generated target and a position of the extracted target signal and a state vector predicted value computing a x t (k + 1)) of Innovation values (v j, t (k + 1) differ) and effective to produce an effective matrix (Ω) from the innovation value (v j, t (k + 1) ) A matrix generation step; And means for deriving a generator matrix from the valid matrix Ω and using the generator matrix to calculate the computed state vector measurements x j (k + 1) An association probability generation step of generating a probability (? J, t ); . ≪ / RTI >

여기서, 상기 연관 확률 생성단계에서는, 상기 도출된 발생 가능 행렬(θ)을 결합 확률 데이터 연관 기법(JPDA : Joint Probability Data-Association)에 적용하여, 상태 벡터 측정치(xj(k+1))가 상기 표적에 연관될 데이터 연관 확률(βj,t)을 생성할 수 있다. 이때, 상기 결합 확률 데이터 연관 기법(JPDA)을 적용한 데이터 연관 확률(βj,t)은 다음 식에 의해 산출될 수 있다.Here, in the association probability generation step, the derived probability matrix? Is applied to a joint probability data association (JPDA) to calculate a state vector measurement value x j (k + 1) ( J, t ) to be associated with the target. At this time, the data association probability ( j, t ) to which the joint probability data association technique (JPDA) is applied can be calculated by the following equation.

Figure 112015021898211-pat00001
Figure 112015021898211-pat00001

Figure 112015021898211-pat00002
Figure 112015021898211-pat00002

Figure 112015021898211-pat00003
Figure 112015021898211-pat00003

또한, 상기 연관 확률 생성단계에서는, 상기 도출된 발생 가능 행렬(θ)을 단순 결합 확률 데이터 연관 기법(CJPDA : Cheap Joint Probability Data-Association)에 적용하여, 상기 연산된 상태 벡터(Zj(k+1))가 상기 표적에 연관될 데이터 연관 확률(βj,t)을 생성할 수 있다. 이때, 상기 단순 결합 확률 데이터 연관 기법(CJPDA)을 적용한 데이터 연관 확률(βj,t)은 다음 식에 의해 산출될 수 있다.Also, in the association probability generation step, the derived generated matrix? Is applied to a simple joint probability data association (CJPDA), and the calculated state vector Z j (k + 1)) may generate a data association probability (? J, t ) to be associated with the target. At this time, the data association probability ( j, t ) to which the simple association probability data association technique (CJPDA) is applied can be calculated by the following equation.

Figure 112015021898211-pat00004
Figure 112015021898211-pat00004

Figure 112015021898211-pat00005
And
Figure 112015021898211-pat00005

특히, 상기 연관 확률 생성단계에서, 상기 수신 신호 세기(p1 η(aj t))는 다음 식에 의해 산출될 수 있다.In particular, in the association probability generation step, the received signal strength (p 1 ? (A j t )) can be calculated by the following equation.

Figure 112015021898211-pat00006
Figure 112015021898211-pat00006
And

Figure 112015021898211-pat00007
Figure 112015021898211-pat00007

또한, 상기 항적 갱신단계에서, 상기 표적의 항적 갱신은 상기 표적의 데이터 연관 확률을 상기 추적 필터의 필터 게인 가중치로 입력하여 표적의 상태 벡터 연산치(x(k+1)) 및 오차공분산 행렬(P(k+1))를 연산함으로써, 상기 생성된 표적의 항적을 갱신할 수 있다. 이때, 상기 표적의 상태 벡터 연산치(xt(k+1)) 및 오차공분산 행렬(Pt(k+1))은 다음 식에 의해 산출될 수 있다.(X (k + 1)) and an error covariance matrix (x (k)) by inputting the data association probability of the target as a filter gain weight of the tracking filter, P (k + 1)), thereby updating the wake of the generated target. At this time, the state vector calculation value (x t (k + 1)) and the error covariance matrix (P t (k + 1)) of the target can be calculated by the following equation.

Figure 112016019039665-pat00039
Figure 112016019039665-pat00039
And

Figure 112016019039665-pat00040
Figure 112016019039665-pat00040

본 발명에 의한 수신 신호 세기를 이용한 다중 표적 추적 방법에 의하면, 결합 확률 데이터 연관 기법에 표적의 거리 데이터뿐만 아니라 표적별 수신 신호 세기 데이터를 적용하여 서로 다른 표적의 크기 및 자세에 따라 달라지는 수신 신호 세기를 이용하므로, 다중 표적 추정 방법에서 사용되는 데이터 연관 기법의 성능을 향상시킬 수 있다. According to the multi-target tracking method using the received signal strength according to the present invention, not only the distance data of the target but also the received signal strength data of each target are applied to the combined probability data association technique, The performance of the data association technique used in the multi-target estimation method can be improved.

또한, 본 발명에 의한 수신 신호 세기를 이용한 다중 표적 추적 방법에 의하면, 표적에 대한 수신 신호 세기 데이터로 다수 개 누적된 수신 신호 세기 데이터의 평균을 사용함으로써 수신 신호 세기의 편차로 인해 발생할 수 있는 데이터 연관 확률의 오차를 감소시킬 수 있다.
In addition, according to the multi-target tracking method using the received signal strength according to the present invention, by using the average of the plurality of accumulated received signal strength data as the received signal strength data for the target, data The error of the association probability can be reduced.

도 1 은 종래 기술에 의한 다중 표적 추적 방법을 보인 흐름도.
도 2 는 본 발명에 의한 수신 신호 세기를 이용한 다중 표적 추적 장치를 보인 블록도.
도 3 은 본 발명에 의한 수신 신호 세기를 이용한 다중 표적 추적 방법을 보인 흐름도.
도 4 는 본 발명에 의한 수신 신호 세기를 이용한 다중 표적 추적 방법 중 유효 행렬 생성단계에서 생성된 유효 행렬(Ω)의 일예.
도 5는 본 발명에 의한 수신 신호 세기를 이용한 다중 표적 추적 방법 중 연관 확률 생성단계에서 도출된 발생 가능 행렬(θ)의 일예들.
BRIEF DESCRIPTION OF THE DRAWINGS Figure 1 is a flow chart illustrating a conventional multi-target tracking method.
2 is a block diagram illustrating a multi-target tracking apparatus using received signal strength according to the present invention.
3 is a flowchart illustrating a multi-target tracking method using received signal strength according to the present invention.
4 is an example of an effective matrix (Ω) generated in the effective matrix generation step among the multiple target tracking methods using the received signal strength according to the present invention.
FIG. 5 illustrates examples of a generable matrix θ derived from an association probability generation step among multiple target tracking methods using received signal strengths according to the present invention.

본 명세서에서 사용되는 기술적 용어는 단지 특정한 실시 예를 설명하기 위하여 사용된 것으로, 본 발명을 한정하려는 의도가 아님을 유의해야 한다. 또한, 본 명세서에서 사용되는 기술적 용어는 본 명세서에서 특별히 다른 의미로 정의되지 않는 한, 본 명세서에 개시된 기술 분야에서 통상의 지식을 가진 자에 의해 일반적으로 이해되는 의미로 해석되어야 하며, 과도하게 포괄적이거나, 과도하게 축소된 의미로 해석되지 않아야 한다. 또한, 첨부된 도면은 본 명세서에 개시된 기술의 사상을 쉽게 이해할 수 있도록 하기 위한 것일 뿐, 첨부된 도면에 의해 그 기술의 사상이 제한되는 것은 아니며, 본 명세서에 개시된 기술의 사상 및 기술 범위에 포함되는 모든 변경·균등물 내지 대체물을 포함하는 것으로 이해되어야 한다.
It is noted that the technical terms used herein are used only to describe specific embodiments and are not intended to limit the invention. Also, the technical terms used herein should be interpreted in a sense that is generally understood by those skilled in the art to which the present disclosure relates, unless otherwise specifically defined in the present specification, Or shall not be construed to mean excessively reduced. It is to be understood that both the foregoing general description and the following detailed description are exemplary and explanatory and are intended to provide further explanation of the invention as claimed. But is to be understood as including all modifications, equivalents, and alternatives falling within the scope of the appended claims.

이하에서는, 첨부된 도면을 참조하여 본 발명에 의한 수신 신호 세기를 이용한 다중 표적 추적 방법의 일 실시 예를 보다 상세하게 설명한다.Hereinafter, an embodiment of a multi-target tracking method using received signal strength according to the present invention will be described in detail with reference to the accompanying drawings.

우선, 본 발명에 의한 수신 신호 세기를 이용한 다중 표적 추적 방법에 사용되는 다중 표적 추적 장치를 간단히 설명한다.First, a multi-target tracking device used in a multi-target tracking method using received signal strength according to the present invention will be briefly described.

도 2 는 본 발명에 의한 수신 신호 세기를 이용한 다중 표적 추적 장치를 보인 블록도이다.2 is a block diagram illustrating a multi-target tracking apparatus using received signal strength according to the present invention.

도 2 를 참조하면, 본 발명에 의한 수신 신호 세기를 이용한 다중 표적 추적 장치(1)의 일 실시 예는, 표적 탐지부(1000), 게이트닝부(2000), 데이터 연관부(3000) 및 항적 데이터 관리부(4000)을 포함할 수 있다. 상기 표적 탐지부(1000)는, 레이더(1100)가 탐지한 표적 신호(100)를 획득하고 상기 표적 신호(100)로부터 표적의 위치 데이터(110)를 추출한다. 또한, 상기 표적 탐지부(1000)는, 상기 표적의 위치 데이터(110)로부터 표적의 거리(310), 표적의 속도(320) 및 표적으로부터의 수신 신호 세기(330) 데이터를 포함하는 표적 데이터(300)를 생성한다. 상기 게이트닝부(2000)는, 상기 표적 데이터(300)가 생성된 표적 신호(100) 중 소정의 유효 거리 이내에 존재하는 표적 신호(100)를 추출한다. 여기서, 상기 표적 데이터(300)가 생성된 표적 신호(100) 중 소정의 유효 거리 이내에 존재하는 표적 신호(100)를 추출하는 것을 게이트닝이라고 표현할 수 있다. 상기 데이터 연관부(3000)는, 상기 표적(10)으로부터의 수신 신호 세기 데이터(330)를 이용하여 상기 생성된 표적의 최종 항적 위치(211)와 추출된 표적 신호의 위치(340) 사이의 데이터 연관 확률(700)을 생성한다. 상기 항적 관리부(4000)는, 추적 필터(4100)를 이용하여 표적의 항적(200)을 생성 또는 갱신한다. 더욱 상세하게는, 상기 항적 관리부(4000)는, 추적 필터를 이용하여 상기 추출된 표적의 위치 데이터(110)로부터 상기 표적이 이동한 이동 항적(210)을 생성하거나, 상기 데이터 연관 확률(700)이 입력된 추적 필터(4100)를 이용하여 상기 생성된 표적의 항적을 갱신할 수 있다.
2, an embodiment of the multi-target tracking apparatus 1 using the received signal strength according to the present invention includes a target detection unit 1000, a gate unit 2000, a data association unit 3000, And may include a management unit 4000. The target detection unit 1000 acquires the target signal 100 detected by the radar 1100 and extracts the positional data 110 of the target from the target signal 100. The target detection unit 1000 further includes target data 310 including target distance 310, target speed 320 and received signal strength 330 data from the target from the target position data 110 300). The gate unit 2000 extracts a target signal 100 within a predetermined effective distance from the target signal 100 in which the target data 300 is generated. Here, the extraction of the target signal 100 existing within the predetermined effective distance among the target signals 100 in which the target data 300 is generated can be expressed as gating. The data associating unit 3000 uses the received signal strength data 330 from the target 10 to calculate the data between the final navigation position 211 of the generated target and the position 340 of the extracted target signal And generates an association probability 700. The trace manager 4000 uses the trace filter 4100 to generate or update the trajectory 200 of the target. More specifically, the trace management unit 4000 generates a movement trajectory 210 in which the target moved from the extracted positional data 110 of the target using a tracking filter, The input tracing filter 4100 can be used to update the generated trajectory of the target.

도 3 은 본 발명에 의한 수신 신호 세기를 이용한 다중 표적 추적 방법을 보인 흐름도이고, 도 4 는 본 발명에 의한 수신 신호 세기를 이용한 다중 표적 추적 방법 중 유효 행렬 생성단계에서 생성된 유효 행렬(Ω)을 보인 예시이며, 도 5 는 본 발명에 의한 수신 신호 세기를 이용한 다중 표적 추적 방법 중 연관 확률 생성단게에서 도출된 발생 가능 행렬(θ)을 보인 예시들이다.FIG. 3 is a flowchart illustrating a multi-target tracking method using the received signal strength according to an embodiment of the present invention. FIG. 4 is a flowchart illustrating a method of tracking an effective matrix (Ω) FIG. 5 is a graph showing examples of a generable matrix (?) Derived from an association probability generation step among multiple target tracking methods using received signal strength according to the present invention.

도 3 을 참조하면, 본 발명에 의한 수신 신호 세기를 이용한 다중 표적 추적 방법의 일실시예는, 표적 위치 추출단계(S310), 항적 생성단계(S320), 표적 데이터 생성단계(S321), 게이트닝 단계(S330), 데이터 연관단계(S350) 및 항적갱신단계 (S340,S360), 항적 제거 및 관리 단계(S370) 등을 포함할 수 있다.Referring to FIG. 3, an embodiment of a multi-target tracking method using received signal strength according to the present invention includes a target position extracting step S310, a snapshot generating step S320, a target data generating step S321, Step S330, a data association step S350, and a navigation update step S340, S360, a wake removal and management step S370, and the like.

상기 표적 위치 추출단계(S310)에서는, 레이더(1100)가 탐지한 표적 신호(100)로부터 표적의 위치 데이터(110)를 추출한다. 더욱 상세하게는, 상기 표적 위치 추출단계(S310)에서는, 상기 표적 탐지부(1000)가 이동하는 다수의 표적(10)을 탐지하는 레이더(1100)로부터 상기 이동하는 다수의 표적(10)에 대한 표적 신호(100)를 획득할 수 있다. 그 후, 상기 표적 탐지부(1000)는 상기 표적 신호(100)를 변환하여 표적의 위치 데이터(110)를 추출할 수 있다. 즉, 상기 표적 위치 추출단계(S310)에서는, 상기 표적 탐지부(1000)가 상기 이동하는 다수의 표적(10)에 대한 신호(100)를 탐지하고 상기 탐지된 표적 신호(100)를 표적의 방위각 데이터 또는 표적의 도플러 데이터로 변환할 수 있으며, 상기 변환된 표적의 방위각 데이터 또는 표적의 도플러 데이터를 다시 표적의 위치 데이터(110)로 변환할 수 있다. 여기서, 상기 레이더(1100)로부터 획득하는 표적 신호(100)는 다수 개의 표적(10)에 대한 표적 신호일 수 있다. 물론, 상기 표적 위치 추출 단계(S310)에 앞서 항적의 초기화가 수행된다.
In the target position extracting step S310, the position data 110 of the target is extracted from the target signal 100 detected by the radar 1100. More specifically, in the target location extraction step (S310), the target detection unit (1000) searches for a plurality of moving targets (10) from a radar (1100) The target signal 100 can be obtained. The target detection unit 1000 may then convert the target signal 100 to extract the target location data 110. That is, in the target position extraction step S310, the target detection unit 1000 detects the signal 100 for the moving plurality of targets 10 and outputs the detected target signal 100 to the target azimuth angle Data or target Doppler data, and the azimuth data of the transformed target or the Doppler data of the target can be converted into the target position data 110 again. Here, the target signal 100 obtained from the radar 1100 may be a target signal for a plurality of targets 10. Of course, prior to the target position extraction step (S310), initialization of the wake is performed.

상기 항적 생성단계(S320)에서는, 추적 필터(4100)를 이용하여 상기 추출된 표적의 위치 데이터(110)로부터 상기 표적의 항적(200) 및 상기 표적에 대한 상태 벡터 예측치를 생성한다. 더욱 상세하게는, 상기 항적 생성단계(S320)에서는, 상기 항적 관리부(4000)의 추적 필터(4100)를 이용하여 상기 추출된 각 표적의 위치 데이터(110)로부터 상기 추출된 각 표적이 이동한 이동 항적(210) 및 이동할 항적을 의미하는 상기 표적에 대한 상태 벡터 예측치를 생성할 수 있다. In the trail generating step S320, a trail filter 4100 is used to generate the trajectory 200 of the target and the state vector predictor for the target from the extracted position data 110 of the target. More specifically, in the generation step (S320), the tracing filter (4100) of the navigation controller (4000) is used to move the extracted target from the position data (110) It is possible to generate the state vector predictions for the target, which means the wake 210 and the wake to move.

즉, 상기 항적 생성단계(S320)에서는, 상기 추출된 표적의 위치 데이터(110)로부터 k 번째 시간까지 표적 t 에 대한 상태 벡터(221)를 연산하여 상기 표적(10)이 이동한 항적(210)을 생성할 수 있다. 또한, 상기 항적 생성단계(S320)에서는, 상기 표적의 위치 데이터(110)로부터, 상기 추적 필터(4100)를 이용하여, 상기 표적이 이동할 항적을 의미하는 k+1 번째 시간에서 t 번째 표적에 대한 상태 벡터 예측치 및 오차 공분산 행렬을 생성할 수 있다. That is, in the trail generation step S320, the state vector 221 for the target t is calculated from the extracted positional data 110 of the target to the kth time, and the trajectory 210, Can be generated. In the trail generating step S320, the tracing filter 4100 is used to calculate, from the target position data 110, the trajectory of the t-th target in the (k + 1) A state vector prediction value and an error covariance matrix.

여기서, k+1 번째 시간의 t 번째 표적에 대한 상태 벡터 예측치를 구하기 위한 k 번째 시간 동안 t 번째 표적까지의 거리의 집합은 k 번째 시간까지 t 번째 표적의 위치 데이터(110)로부터 연산될 수 있다. 여기서, 상기 표적에 대한 상태 벡터 예측치는 k+1 번째 시간의 t 번째 표적의 위치와 k 번째 시간의 t 번째 표적의 위치로부터 연산된 상태 벡터를 의미할 수 있으며, 상기 추적 필터(4100)는, 상기 표적의 상태 벡터 예측치 및 오차공분산 행렬과 상기 표적의 수신 신호 세기 데이터(330)를 반영한 상태 벡터 연산치 및 오차공분산을 도출할 수 있는 칼만 필터일 수 있다.Here, the set of distances up to the t-th target during the k-th time to obtain the state vector predictions for the t-th target in the (k + 1) -th time can be calculated from the position data 110 of the t- . Here, the state vector prediction value for the target may be a state vector calculated from the position of the t-th target in the (k + 1) -th time and the position of the t-th target in the k-th time. The tracking filter 4100, A Kalman filter that can derive a state vector calculation value and an error covariance reflecting the target state vector prediction value and the error covariance matrix and the target received signal strength data 330.

이때, 상기 표적의 위치 데이터(110)로부터 k 번째 시간까지 t 번째 표적에 대한 거리의 집합이 연산되는 경우, K+1 번째 시간의 t 번째 표적에 대한 상태 벡터에 대한 예측치(xt(k+1)) 및 K+1 번째 시간의 t 번째 표적의 상태 벡터에 대한 오차 공분산 행렬(Pt(k+1))는 다음 식,In this case, when the set of distances for the t-th target from the target position data 110 to the k-th time is calculated, the predicted value x t (k + 1) for the state vector for the t- 1)) and the error covariance matrix P t (k + 1) for the state vector of the t-th target of the (K + 1)

Figure 112015021898211-pat00010
Figure 112015021898211-pat00010

Figure 112015021898211-pat00011
Figure 112015021898211-pat00011

에 의해 산출될 수 있다.Lt; / RTI >

여기서, 상기 xt(k|k)는 k 번째 시간의 t 번째 표적의 상태 벡터로 xt(k)=[zX,k,zY,k]로 표현되는 k 번째 시간의 t 번째 표적의 위치 정보일 수 있다. 또한, 상기 F(k) 는 표적의 동적 운동 모델로 상태 변환 행렬이고, 상기 P(k|k) 는 k 번째 시간에서의 표적의 오차공분산 행렬이며, 상기 F'(k)는 상기 상태 변환 행렬의 역행렬이고, 상기 Q(k) 는 표적의 동적 운동 모델에서 과정 잡음을 나타내는 공분산 행렬이다.Here, x t (k | k) is a state vector of the t-th target of the k-th time, which is expressed as x t (k) = [z X, k , z Y, k ] Location information. (K) is an error covariance matrix of the target at the k-th time, and F '(k) is a state transformation matrix of the target dynamic motion model, (K) is a covariance matrix representing the process noise in the target dynamic motion model.

즉, 상기 [수학식 1] 및 [수학식 2]를 통하여, 상기 항적 생성단계(S320)에서는, 상기 추출된 표적의 위치 데이터(110)로부터 도출된 K 번째 시간까지의 t 번째 표적의 상태 벡터(xt(k))를 이용하여 K+1 번째 시간까지의 t 번째 표적의 상태 벡터(xt(k+1)) 예측치 및 오차 공분산 행렬(Pt(k+1))을 연산할 수 있다.That is, through the above-mentioned [Expression 1] and [Equation 2], in the trail generation step S320, the state vector of the t-th target from the extracted positional data 110 of the target up to the K- (x t (k)) to calculate a k + 1 th of the time t the second target to the state vector (x t (k + 1) ) predicted value and the error covariance matrix (P t (k + 1) ) by using the have.

상기 표적 데이터 생성단계(S321)에서는, 레이더 센서(1100)를 통해 상기 추출된 표적의 위치 데이터(110)로부터 표적까지의 거리, 표적의 속도 및 표적으로부터의 수신 신호 세기 데이터를 포함하는 표적 데이터(300)를 생성한다. 더욱 상세하게는, 상기 표적 데이터 생성단계(S321)에서는, 상기 이동하는 다수의 표적(10)에 대한 표적 신호(100)의 표적의 위치 데이터(110)로부터 상기 생성된 표적의 항적(200)을 갱신하는데 이용되는 표적까지의 거리 데이터(310), 표적의 속도 데이터(320) 및 표적으로부터의 수신 신호 세기 데이터(330)를 생성할 수 있다. In the target data generation step S321, the target data including the distance from the position data 110 of the extracted target to the target, the speed of the target, and the received signal strength data from the target through the radar sensor 1100 300). More specifically, in the target data generation step (S321), the generated wakefulness (200) of the generated target from the position data (110) of the target of the target signal (100) The distance data 310 to the target used to update, the velocity data 320 of the target, and the received signal strength data 330 from the target.

즉, 상기 표적 데이터 생성단계(S321)에서는, 상기 추출된 표적의 위치 데이터(110)로부터 상기 이동하는 다수의 표적(10)에 대한 각 표적까지의 거리(310), 각 표적의 속도(320) 및 각 표적으로부터의 수신 신호 세기 데이터(330)를 생성하는 것이다. 여기서, 상기 표적까지의 거리 데이터(310)와 상기 t 번째 표적까지의 거리는 동일할 수 있다. 또한, 여기서, 상기 표적으로부터의 수신 신호 세기 데이터(330)는 표적의 종류, 표적의 자세에 따라 달라질 수 있다.That is, in the target data generation step S321, the distance 310 from the extracted target position data 110 to each target of the moving target 10, the velocity 320 of each target, And the received signal strength data 330 from each target. Here, the distance data 310 to the target and the distance to the t-th target may be the same. Here, the received signal strength data 330 from the target may vary depending on the type of the target and the posture of the target.

상기 게이트닝 단계(S330)에서는, 상기 표적 데이터(300)가 생성된 표적 신호(100) 중 상기 생성된 표적의 최종 항적 위치(211)로부터 소정의 유효 거리 이내에 존재하는 표적 신호(100)를 추출한다. 더욱 상세하게는, 상기 게이트닝 단계(S330)에서는, 상기 K+1 번째 시간에서의 t 번째 표적의 위치를 추정하기 위하여, 상기 표적 데이터(300)가 생성된 표적 신호(100) 중 k 번째 시간에서의 t 번째 표적의 최종 항적 위치(211)를 중심으로 소정의 반경 이내에 존재하는 표적 신호(100)만을 추출할 수 있다. In the gatewaying step S330, the target signal 100 existing within a predetermined effective distance from the final navigation position 211 of the generated target among the generated target signals 100 is extracted from the target data 300 do. More specifically, in the gateing step S330, in order to estimate the position of the t-th target in the (K + 1) -th time, the target data 300 is generated at the k- It is possible to extract only the target signal 100 existing within a predetermined radius centering on the final trajectory position 211 of the t-th target in FIG.

즉, 상기 게이트닝 단계(S330)에서는, 레이더 센서(1100)에 의해 탐지되어 생성된 표적의 과거 항적(210)으로부터 상기 표적(10)이 이동할 다음 위치를 추정하기 위하여, 레이더에 의해 탐지된 표적 신호(100) 중 상기 표적의 과거 항적(210)에서 표적(10)의 최종 위치(211)를 중심으로 소정의 반경 이내에 존재하는 표적 신호(100)만을 추출하는 것이다. 이때, 상기 레이더에 의해 탐지된 표적 신호(100)는 표적(10)의 실제 개수보다 많을 수 있으며, 예를 들어, 2 개의 이동하는 표적(10)에 대하여 10 개의 표적 신호(100)가 레이더로 탐지될 수 있다. That is, in the gatewaying step S330, in order to estimate the next position at which the target 10 moves from the past warnings 210 of the target detected and generated by the radar sensor 1100, Only the target signal 100 existing within a predetermined radius around the final position 211 of the target 10 in the past trajectory 210 of the target 100 is extracted. At this time, the target signal 100 detected by the radar may be larger than the actual number of the target 10, for example, ten target signals 100 for two moving targets 10 may be radar Can be detected.

더욱 상세하게는, 상기 게이트닝 단계(S330)에서는, 레이더 센서(1100)에 의해 획득된 k 번째 시간까지 하나 이상의 표적 신호(100)로부터 생성된 하나 이상의 표적의 항적(210)으로부터 각 표적이 이동할 다음 위치를 추정하기 위하여, k+1 번째 시간에서 레이더에 의해 탐지된 다수의 표적 신호(100) 중 상기 각 표적의 과거 항적(210)에서의 각 표적(10)의 최종 위치(211), 즉, k 번째 시간에서의 각 표적의 위치(211)를 중심으로 소정의 반경 이내에 존재하는 표적 신호(100)만을 각각 추출하는 것이다. More specifically, in the gating step S330, each target moves from the wake of one or more targets 210 generated from one or more target signals 100 until the kth time obtained by the radar sensor 1100 In order to estimate the next position, the final position 211 of each target 10 in the past trail 210 of each target among the plurality of target signals 100 detected by the radar at the (k + 1) , and extracts only the target signal 100 existing within a predetermined radius around the position 211 of each target in the k-th time.

여기서, 상기 표적 데이터(300)가 생성된 표적의 신호(100) 중 소정의 유효 거리 이내에 존재하는 표적 신호(100)를 추출하는 것을 게이트닝이라고 표현할 수 있다. 또한, 여기서, 상기 소정의 유효 거리는, k+1 번째 시간에서 상기 표적이 이동할 수 있는 위치를 제한하는 거리로 상수일 수 있다.Here, the extraction of the target signal 100 existing within the predetermined effective distance among the signals 100 of the target in which the target data 300 is generated can be expressed as gating. Here, the predetermined effective distance may be a constant that limits a position at which the target can move in a (k + 1) -th time.

이때, k+1 번째 시간에서 레이더에 의해 탐지된 다수의 표적 신호(100) 중 소정의 유효 거리 이내에 존재하는 표적 신호(100)의 개수를 m 개라 하는 경우, k+1 번째 시간에서 소정의 유효 거리 이내에 존재하는 t 번째 표적의 상태 벡터(xt(k+1))(222)는 다음 식,In this case, when the number of target signals 100 existing within a predetermined effective distance among the plurality of target signals 100 detected by the radar at the (k + 1) th time is m, The state vector (x t (k + 1)) 222 of the t-th target existing within the distance is given by the following equation,

Figure 112015021898211-pat00012
Figure 112015021898211-pat00012

Figure 112015021898211-pat00013
Figure 112015021898211-pat00013

으로 표현될 수 있다.. ≪ / RTI >

여기서, 상기 Zt m(k+1)는 m 번째 표적 신호에서 k+1 번째 시간의 t 번째 표적에 대한 위치정보(zX,k+1 , zY,k+1) 및 속도 정보(z'X,k+1,z'Y,k+1)를 의미하고, 상기 Zt a(k+1)는 m 번째 표적 신호에서 k+1 번째 시간의 t 번째 표적으로부터 수신된 신호 세기 데이터(330)를 의미한다. 이때, 상기 수신 신호 세기 데이터( Zt a(k+1))(330)로부터 t 번째 표적의 종류, 표적의 크기 및 표적의 자세 정보를 유추할 수 있으며, 상기 데이터 연관단계(S350)에서는 상기 수신 신호 세기 정보( Zt a(k+1))(330)를 이용하여 서로 다른 표적(10) 사이에 상이하게 나타나는 수신 신호 세기의 분포 및 평균 신호대 잡음비를 통해 표적의 데이터 연관 확률의 신뢰성을 높일 수 있다.Here, the Z t m (k + 1) is the m-th in the target signal position information for the k + 1 t-th sign of the second time (z X, k + 1, z Y, k + 1) and speed information (z 'X, k + 1, z ' Y, wherein Z t means the k + 1), and a (k + 1) is the signal strength data received from the k + 1 t-th sign of the second time in an m-th target signal ( 330). At this time, the type of the t-th target, the size of the target, and the attitude information of the target can be deduced from the received signal strength data (Z t a (k + 1)) 330, The reliability of the data association probability of the target is obtained through the distribution of the received signal strengths and the average signal-to-noise ratio, which are different between the different targets 10, using the received signal strength information Z t a (k + .

상기 데이터 연관단계(S350)에서는, 상기 표적(10)으로부터의 수신 신호 세기 데이터(330)를 이용하여 상기 생성된 표적의 최종 항적 위치(211)와 상기 소정의 유효 거리 이내에 존재하는 표적 신호 위치(340) 사이의 데이터 연관 확률을 생성한다. 즉, 상기 데이터 연관단계(S350)에서는, 상기 각 표적(10)으로부터의 수신 신호 세기 데이터(330)를 이용하여 상기 생성된 각 표적의 최종 항적 위치(211)와 상기 소정의 유효 거리 이내에 존재하는 m 개의 표적 신호 위치(340) 사이의 데이터 연관 확률을 생성할 수 있다. 더욱 상세하게는, 상기 데이터 연관단계(S350)는, 유효 행렬 생성단계(미도시) 및 연관 확률 생성단계(미도시)를 포함할 수 있다.In the data associating step S350, the target signal position (211) existing within the predetermined effective distance from the final navigation position (211) of the generated target using the received signal strength data 330 from the target (10) 340). ≪ / RTI > That is, in the data associating step S350, it is determined whether or not there is a difference between the final navigation position 211 of each generated target using the received signal strength data 330 from each of the targets 10, m < / RTI > More specifically, the data association step S350 may include an effective matrix generation step (not shown) and an association probability generation step (not shown).

상기 유효 행렬 생성단계에서는, 상기 생성된 표적의 최종 항적 위치(211)와 상기 m 개의 추출된 표적 신호의 위치(340)로부터 연산된 상태 벡터 측정치(xj(k+1))와 상기 생성된 표적에 대한 상태 벡터 예측치(xt(k+1)) 차이인 이노베이션 값(vj,t(k+1))을 연산하고 상기 이노베이션 값(vj,t(k+1))으로 유효 행렬(Ω)을 생성한다. 즉, 상기 유효 행렬 생성단계(S510)에서는, 상기 생성된 각 표적의 최종 항적 위치(211)와 상기 m 개의 추출된 표적 신호의 위치(340)로부터 각각 연산된 상태 벡터(xj(k+1))와 상기 생성된 각 표적에 대한 상태 벡터 예측치(xt(k+1)) 사이의 차이를 의미하는 이노베이션 값(vj,t(k+1))을 연산하고 상기 이노베이션 값(vj,t(k+1))으로부터 유효 행렬(Ω)(500)을 생성한다.In the valid matrix generation step, the state vector measurement value x j (k + 1) calculated from the final navigation position 211 of the generated target and the position 340 of the m extracted target signals, The innovation value v j, t (k + 1)) is calculated by computing an innovation value v j, t (k + 1) which is the difference between the state vector prediction value x t (?). That is, in the effective matrix generation step (S510), the state vector (x j (k + 1)) calculated from the final navigation position 211 of each generated target and the position 340 of the m extracted target signals, )) and the generated status for each target vector predicted values (Innovation value indicating a difference between x t (k + 1)), (v j, t (k + 1)) for operation and the innovation value (v j , t (k + 1)).

이때, 상기 유효 행렬(Ω)(500)은 다음 식,At this time, the effective matrix (?) (500)

Figure 112015021898211-pat00014
Figure 112015021898211-pat00014

에 의해 산출될 수 있다.Lt; / RTI >

여기서, 상기 wj,t 는 상기 유효 행렬의 엘리먼트로써 j 번째 추출된 표적 신호의 위치가, k 번째 시간에서 t 번째 표적을 중심으로, 소정의 유효 범위 내에 존재하면 "1" 이고 그렇지 않을 경우 "0" 인 엘리먼트이고, 상기 m 은 상기 소정의 유효 범위 내에 존재하는 표적 신호(100)의 개수이며, 상기 T 는 항적을 생성한 표적의 개수이고, 상기 r 는 소정의 유효 범위 값을 의미이다.Here, w j, t denotes "1" if the position of the j-th extracted target signal as an element of the valid matrix is within a predetermined effective range centered on the t-th target in the k-th time, 0 ", where m is the number of target signals 100 present within the predetermined effective range, T is the number of targets that generated the trajectory, and r is a predetermined effective range value.

상기 연관 확률 생성단계(미도시)에서는, 상기 유효 행렬(Ω)(500)로부터 생성 가능한 발생 가능 행렬(θ)(510)(520)(530)(540)(550)(560)(570)(580)을 도출하고 상기 발생 가능 행렬(θ)을 이용하여 상기 연산된 상태 벡터(xj(k+1))가 상기 표적에 연관될 데이터 연관 확률(βj,t)(700)을 생성한다. 510, 520, 530, 540, 550, 560, 570 that can be generated from the effective matrix (?) 500 in the association probability generation step (not shown) (? J, t ) (700) to be associated with the target by using the generated state vector (x j (k + 1)) using the generator matrix do.

더욱 상세하게는, 도 4 및 도 5 를 참조하면, 상기 연관 확률 생성단계(미도시)에서, 2 개의 표적(t) 및 3 개의 표적 신호(j)가 소정의 유효 범위 이내에 존재하는 경우 유효 행렬(Ω)(500)은 도 4 와 같고, 상기 유효 행렬(Ω)(500)로부터 발생 가능 행렬(θ)(510)(520)(530)(540)(550)(560)(570)(580)은 도 5와 같이 도출될 수 있다. 이때, 상기 발생 가능 행렬을 도출하는 방법은, 하나의 표적(t)에는 하나의 표적 신호(j)만이 존재하고 하나의 표적 신호(j)는 하나의 표적(t)에 대해서만 존재하도록 상기 발생 가능 행렬을 도출하는 것이다.4 and 5, in the association probability generation step (not shown), when two targets t and three target signals j are within a predetermined effective range, 5 500, 530, 550, 560, 570, 530, 530, 530, 530, 530, 530, 580 may be derived as shown in FIG. At this time, the method for deriving the generable matrix is such that only one target signal (j) exists in one target (t) and one target signal (j) exists only in one target To derive a matrix.

또한, 상기 연관 확률 생성단계(미도시)에서는, 상기 도출된 발생 가능 행렬(θ)(510)(520)(530)(540)(550)(560)(570)(580)을 결합 확률 데이터 연관 기법(JPDA : Joint Probability Data-Association)에 적용하여, 상기 연산된 상태 벡터(xj(k+1))(400)가 상기 표적에 연관될 데이터 연관 확률(βj,t)(700)을 생성할 수 있다.510, 520, 530, 540, 550, 560, 570, and 580 of the derived probability matrix (not shown) (? J, t ) (700) to be associated with the target by applying the computed state vector (x j (k + 1)) (400) to the Joint Probability Data- Can be generated.

이때, 결합 확률 데이터 연관 기법을 적용한 상기 데이터 연관 확률(βj,t)(700)은 다음 식, At this time, the data association probability (β j, t ) (700), to which the joint probability data association technique is applied,

Figure 112015021898211-pat00015
Figure 112015021898211-pat00015

Figure 112015021898211-pat00016
Figure 112015021898211-pat00016

에 의해 산출될 수 있다.Lt; / RTI >

여기서, 상기 c 는 정규화 상수(normalized coefficent)이고, 상기 V 는 유효 영역 구간으로 πr|S(k)|1/2로 정의되며, r은 게이트 영역을 결정하는 게이트닝 임계치 상수를 의미하며, 상기 φ 는 현 상태에서 신호 개수로 여기서는 m-1 이고, 상기 N[vj,t(k+1);0;S(k+1)] 는 vj,t 의 평균이 0 이고 분산이 S(k+1) 인 가우시안 분포이며, 상기 τj(θ), δt 는 모든 발생 가능 행렬(θ)에서의 t 번째 표적(트랙)이 연관된 경우와 j 번째 추출된 표적 신호가 연관된 경우에만 1 값을 갖는 변수이고, Pt D 는 t 번째 표적(트랙)의 검출 확룔로 사전에 정의된 시스템 변수이다. 그리고, 상기 p1 η(at j) 는 j 번째 추출된 표적 신호의 수신 신호 세기이고, 상기 at j 는 PDF(Probability density function)이며, 상기 p0 η(at j) 는 노이즈 신호의 수신 신호 세기를 의미한다.Here, c is a normalized coefficient, V is a valid region interval, and πr | S (k) | 1/2 , where r is the gateing threshold constant for determining the gate region, and φ is the number of signals in the current state, m-1, where N [v j, t (k + 1); 0; S (k + 1) ] is a v j, which is a Gaussian distribution is zero mean and variance of the t S (k + 1), wherein τ j (θ), δ t from all generated available matrix (θ) Is a variable having a value of 1 only when the t-th target (track) of the t-th target (track) is associated with the j-th extracted target signal and P t D is a previously defined system variable defined as the detection probability of the t-th target (track). In addition, the p 1 η of (a t j) is the j-th and the received signal strength of the extracted target signal, wherein a t j is a PDF (Probability density function), the p 0 η (a t j) is the noise signal Means the received signal strength.

이때, 상기 j 번째 추출된 표적 신호의 수신 신호 세기(p1 η(at j)) 및 노이즈 신호의 수신 신호 세기(p0 η(at j))는 다음 식,At this time, the received signal strength 1 p (a t j)) and the received signal strength of the noise signal (p 0 η (a t j )) of the j-th extracted the target signal is obtained using the equation,

Figure 112015021898211-pat00017
Figure 112015021898211-pat00017

에 의해 산출될 수 있다.Lt; / RTI >

여기서, 상기 d는 표적 신호의 평균 신호대 잡음비를 나타내고 FMCW(Frequency Modulated Continous Wave) 레이더로부터 받은 수신 신호 세기에 비례하는 값을 가진다.
Here, d represents the average signal-to-noise ratio of the target signal and has a value proportional to the received signal strength received from the FMCW (Frequency Modulated Continuous Wave) radar.

본 발명은 표적별 신호 세기를 다르게 적용하는 것이므로, [수학식 10] 및 [수학식 11] 을 상기 [수학식 8]에 적용하게 되면 표적별 상이한 신호대 잡음비(λi)가 상기 데이터 연관 확률(βj,t)(700)에 적용되어, 표적으로부터의 수신 신호 세기 데이터가 반영된 데이터 연관 확률(βj,t)(700)을 획득할 수 있다.Since the present invention applies different target signal strengths, applying Equation (10) and Equation (11) to Equation (8) results in different signal-to-noise ratios (λ i ) β j, t ) 700 to obtain a data association probability β j, t 700 reflecting the received signal strength data from the target.

또한, 상기 연관 확률 생성단계(S350)에서는, 상기 발생가능 행렬(θ)(510)(520)(530)(540)(550)(560)(570)(580)을 단순 결합 확률 데이터 연관 기법(CJPDA : Cheap Joint Probability Data-Association)에 적용하여, 상기 연산된 상태 벡터(xj(k+1))가 상기 표적에 연관될 데이터 연관 확률(βj,t)을 생성할 수 있다.510, 520, 530, 540, 550, 560, 570, and 580 may be combined with a simple joint probability data association technique ( J j, t ) to be associated with the target by applying the computed state vector x j (k + 1) to the CJPDA (Cheap Joint Probability Data Association).

이때, 단순 결합 확률 데이터 연관 기법을 적용한 상기 데이터 연관 확률(βj,t)은 다음 식, At this time, the data association probability (? J, t ) to which the simple association probability data association technique is applied is expressed by the following equation

Figure 112015021898211-pat00018
Figure 112015021898211-pat00018

에 의해 산출될 수 있다.Lt; / RTI >

여기서, 상기 λj 는 신호 대 잡음비(Amplitude likelihood ratio)로서

Figure 112015021898211-pat00019
이고, 상기 p1 η(aj t)는 j 번째 추출된 표적 신호의 수신 신호 세기이며, 상기 p0 η(at j) 는 노이즈 신호의 수신 신호 세기이다.Here, the above-mentioned? J is an amplitude likelihood ratio
Figure 112015021898211-pat00019
And wherein η 1 p (a j t) is the received signal strength of the j-th extracted the target signal, the η 0 p (a t j) is the received signal strength of the noise signal.

여기서, 상기 수신 신호 세기 (p1 η(aj t))는 다음 식,Here, the received signal strength p 1 ? (A j ? T )

Figure 112015021898211-pat00020
Figure 112015021898211-pat00020

에 의해 산출될 수 있다.Lt; / RTI >

여기서, 상기 d 는 표적 신호의 평균 신호대 잡음비를 나타내고 FMCW 레이더로부터 받은 수신 신호 세기에 비례하는 값을 가진다.Here, d represents the average signal-to-noise ratio of the target signal and is proportional to the received signal strength received from the FMCW radar.

본 발명은 표적별 신호 세기를 다르게 적용하는 것이므로, 상기 [수학식 10] 을 상기 [수학식 9]에 적용하게 되면 표적별 상이한 신호대 잡음비(λi)가 상기 데이터 연관 확률(βj,t)에 적용되어, 표적으로부터의 수신 신호 세기 데이터가 반영된 데이터 연관 확률(βj,t)을 획득할 수 있다. The invention, because of applying a different target-specific signal strength, wherein if the [Equation 10] applied to the above [Equation 9] a different target-specific signal-to-noise ratio (λ i) is the data associated probability (β j, t) To obtain the data association probability (? J, t ) in which the received signal strength data from the target is reflected.

여기서, 상기 각 표적별 신호대 잡음비(λi)를 도출하기 위하여 특정 시간동안 누적된 수신 신호 세기의 평균값(dt)은 다음 식,Here, an average value (d t ) of the received signal strength accumulated for a specific time to derive the signal-to-noise ratio (λ i )

Figure 112015021898211-pat00021
Figure 112015021898211-pat00021

에 의해 산출될 수 있다.Lt; / RTI >

여기서, 상기 L 은 신호 세기 평균 연산을 위한 누적 시간을 의미한다.Here, L denotes an accumulated time for calculating the signal intensity average.

또한, 단순 결합 확률 데이터 연관기법에서 연관확률을 구하는 [수학식 9]뿐만이 아니라, 결합확률 데이터 연관기법에서 연관확률을 구하는 [수학식 8] 에 표현된 평균신호대잡음비 d도 표적별 서로 다른 신호세기를 반영한 [수학식 11] 적용이 가능하다. The average signal-to-noise ratio d expressed in Equation (8) for obtaining the association probability in the joint probability data association technique as well as Equation (9) for obtaining the association probability in the simple joint probability data association technique may be different from each other (11) can be applied.

상기 항적 갱신단계(S360)에서는, 상기 표적의 데이터 연관 확률을 추적 필터(4100)에 입력하여 상기 생성된 표적 항적(200)을 갱신한다. 더욱 상세하게는, 상기 항적 갱신단계(S360)에서는, 상기 결합 확률 데이터 연관 기법 또는 단순 결합 확률 데이터 연관 기법에 의해 생성된 상기 데이터 연관 확률을 추적 필터(4100)에 입력하여 k+1 번째 시간에서 t 번째 표적에 대한 상태 벡터 연산치(x(k+1)) 및 상기 상태 벡터의 오차공분산 행렬(P(k+1))을 구함으로써 상기 생성된 표적의 항적(200)을 갱신할 수 있다.In the navigation updating step S360, the data association probability of the target is input to the tracking filter 4100 to update the generated target wake 200. [ More specifically, in the navigation update step (S360), the data association probability generated by the joint probability data association technique or the simple association probability data association technique is input to the tracking filter 4100, the generated trajectory 200 of the generated target can be updated by obtaining the state vector operation value x (k + 1) for the t-th target and the error covariance matrix P (k + 1) of the state vector .

이때 상기 k+1 번째 시간에서 t 번째 표적에 대한 상태 벡터 연산치(x(k+1)) 및 상기 상태 벡터의 오차공분산 행렬(P(k+1))은 다음 식,The state vector calculation value x (k + 1) for the t-th target in the (k + 1) th time and the error covariance matrix P (k + 1)

Figure 112016019039665-pat00041
Figure 112016019039665-pat00041

Figure 112016019039665-pat00042
Figure 112016019039665-pat00042

에 의해 산출될 수 있다.Lt; / RTI >

여기서, 상기 S(k+1)는 이노베이션의 오차공분산 행렬이고, 상기 S-1(k+1)는 상기 S(k+1)의 역함수이며, 상기 W(k+1)은 추적 필터의 필터게인으로 본 발명에서는 칼만 필터의 필터게인이고, 상기 R(k+1)은 측정잡음의 오차공분산 행렬이며, 상기 H(k+1)은 측정치와 상태 변수의 관계를 나타내는 상태변환 행렬이고, 상기 H'(k+1) 은 H(k+1)의 역행렬이며, 상기 vt(k)는 이노베이션의 가중합이고, 상기 β0 는 데이터 연관 확률이며, 상기 I 는 단위행렬이다.Here, the S (k + 1) is the error covariance matrix of the innovation, is the inverse of the S -1 (k + 1) is the S (k + 1), wherein W (k + 1) is a filter of the tracking filter (K + 1) is an error covariance matrix of the measurement noise, H (k + 1) is a state transformation matrix indicating a relationship between a measured value and a state variable, H (k + 1) is an inverse matrix of H (k + 1), v t (k) is a weighted sum of Innovation, β 0 is a data association probability and I is a unit matrix.

상기 [수학식 12] 를 참조하면, 상기 항적 갱신단계(S360)에서, 상기 표적의 항적 갱신은 상기 표적의 데이터 연관 확률을 상기 추적 필터(4100)의 필터 게인 가중치(W(K+1))로 입력하여 표적의 상태 벡터 연산치(x(k+1)) 및 오차공분산 행렬(P(k+1))을 연산함으로써, 상기 생성된 표적 항적(200)을 갱신할 수 있다.In the navigation update step S360, the target navigation update of the target is performed based on the filter gain weight W (K + 1) of the tracking filter 4100 by referring to Equation (12) The generated target wakefull 200 can be updated by calculating the target state vector calculation value x (k + 1) and the error covariance matrix P (k + 1).

상기 항적 갱신단계(S360)에서 그렇게 생성된 추적필터 시스템 상수들을 기반으로 항적이 갱신된다. 항적 제거 및 관리 단계(S370)에서 갱신된 항적 결과에 대한 항적 유지 및/또는 제거 과정을 거치게 된다. 물론, 항적 제거 및 관리 단계(S370)이전에 새 트랙을 생성하는 단계(S340)가 구성될 수 있다.The navigation is updated based on the tracking filter system constants so generated in the navigation update step (S360). In the trajectory elimination and management step (S370), trajectories are maintained and / or eliminated for the updated trajectory results. Of course, a step S340 of creating a new track may be configured before the track removal and management step S370.

상술한 바와 같이, 본 발명의 기본적인 기술적 사상의 범주 내에서, 당업계의 통상의 지식을 가진 자에게 있어서는 다른 많은 변형이 가능하고, 본 발명의 권리범위는 첨부한 특허청구범위에 기초하여 해석되어야 할 것이다.It will be apparent to those skilled in the art that many other modifications and variations are possible in light of the above teachings and the scope of the present invention should be construed on the basis of the appended claims something to do.

1 : 본 발명에 의한 수신 신호 세기를 이용한 다중 표적 추적 장치
10 : 표적
100 : 표적의 신호
110 : 표적의 위치 데이터 111 : t 번째 표적까지의 거리
200 : 표적 항적
210 : 표적이 이동한 이동 항적
211 : 표적이 이동한 항적에서 표적의 최종 항적 위치
220 : 표적이 이동할 항적
300 : 표적 데이터
310 : 표적까지의 거리 데이터 320 : 표적의 속도 데이터
330 : 표적으로부터의 수신 신호 세기 데이터
331 : j 번째 추출된 표적 신호의 신호 세기
332 : 노이즈 신호의 신호 세기
340 : 소정의 유효 거리 내에 존재하는 표적 신호 위치
500 : 유효 행렬(Ω)
510, 520, 530, 540, 550, 560, 570, 580 : 유효 행렬(500)로부터 생성 가능한 발생가능 행렬(θ)
700 : 데이터 연관 확률(βj,t)
1000 : 표적 탐지부 1100 : 레이더 센서
2000 : 게이트닝부
3000 : 데이터 연관부
4000 : 항적 관리부
4100 : 추적 필터
1: Multi-target tracking device using received signal strength according to the present invention
10: Target
100: signal of the target
110: Position data of the target 111: Distance to the t-th target
200: Target wake
210: Movement trajectory to which the target moved
211: the final trajectory position of the target in the trajectory to which the target moved
220: Wand to move the target
300: target data
310: distance data to the target 320: velocity data of the target
330: received signal strength data from the target
331: Signal strength of the j-th extracted target signal
332: Signal strength of noise signal
340: Target signal position existing within a predetermined effective distance
500: effective matrix (Ω)
510, 520, 530, 540, 550, 560, 570, 580: generating matrix (?) That can be generated from the effective matrix (500)
700: Data association probability (? J, t )
1000: target detection unit 1100: radar sensor
2000: Gating section
3000:
4000:
4100: Trace filter

Claims (10)

다중 표적 추적 방법에 있어서,
레이더 센서가 탐지한 표적 신호로부터 표적의 위치 데이터를 추출하는 표적 위치 추출단계;
추적 필터를 이용하여 상기 추출된 표적의 위치 데이터로부터 상기 표적의 항적 및 상기 표적에 대한 상태 벡터 예측치를 생성하는 항적 생성단계;
상기 추출된 표적의 위치 데이터로부터 표적의 거리, 표적의 속도 및 표적으로부터의 수신 신호 세기 데이터를 포함하는 표적 데이터를 생성하는 표적 데이터 생성단계;
상기 표적 데이터가 생성된 표적 신호 중 상기 수신 신호 세기 데이터를 이용하여 생성된 표적의 최종 항적 위치로부터 소정의 유효 거리 이내에 존재하는 표적 신호를 추출하는 게이트닝 단계;
상기 표적으로부터의 수신 신호 세기 데이터를 이용하여 상기 생성된 표적의 최종 항적 위치와 추출된 표적 신호의 위치 사이의 데이터 연관 확률을 생성하는 데이터 연관단계; 및
상기 표적의 데이터 연관 확률을 상기 추적 필터에 입력하여 상기 생성된 표적의 항적을 갱신하는 항적 갱신단계;
를 포함하는 것을 특징으로 하는 수신 신호 세기를 이용한 다중 표적 추적 방법.
In a multiple target tracking method,
A target position extracting step of extracting target position data from the target signal detected by the radar sensor;
Generating a trajectory of the target and a state vector prediction for the target from the extracted positional data of the target using a tracking filter;
A target data generation step of generating target data including the distance of the target, the speed of the target, and the received signal strength data from the target, from the extracted positional data of the target;
A gating step of extracting a target signal existing within a predetermined effective distance from a final navigation position of a target generated using the received signal strength data of the target signals from which the target data is generated;
A data associating step of generating a data associating probability between the final navigation position of the generated target and the position of the extracted target signal using the received signal strength data from the target; And
A navigation update step of inputting the data association probability of the target into the tracking filter and updating the generated map of the generated target;
Wherein the multi-target tracking method is based on received signal strength.
제 1 항에 있어서,
상기 표적 데이터 생성단계에서,
상기 표적으로부터의 수신 신호 세기 데이터는 표적의 종류, 표적의 자세에 따라 달라지는 것을 특징으로 하는 수신 신호 세기를 이용한 다중 표적 추적 방법.
The method according to claim 1,
In the target data generation step,
Wherein the received signal strength data from the target is dependent on the type of target and the attitude of the target.
제 1 항에 있어서,
상기 데이터 연관단계는,
상기 생성된 표적의 최종 항적 위치와 상기 추출된 표적 신호의 위치로부터 연산된 상태 벡터 측정치(xj(k+1))와 상기 생성된 표적에 대한 상태 벡터 예측치(xt(k+1)) 차이인 이노베이션 값(vj,t(k+1))을 연산하고 상기 이노베이션 값(vj,t(k+1))으로부터 유효 행렬(Ω)을 생성하는 유효 행렬 생성단계; 및
상기 유효 행렬(Ω)로부터 발생 가능 행렬(θ)을 도출하고 상기 발생 가능 행렬(θ)을 이용하여 상기 연산된 상태 벡터 측정치(xj(k+1))가 상기 표적에 연관될 데이터 연관 확률(βj,t)을 생성하는 연관 확률 생성단계;
를 포함하는 것을 특징으로 하는 수신 신호 세기를 이용한 다중 표적 추적 방법.
The method according to claim 1,
Wherein the data associating step comprises:
A state vector prediction value x t (k + 1) for the generated target and a state vector measurement value x j (k + 1) calculated from the position of the extracted target signal of the generated target, difference of the effective value Innovation matrix generation step of computing a (v j, t (k + 1)) and to generate an effective matrix (Ω) from the innovation value (v j, t (k + 1)); And
Deriving a generator matrix from the effective matrix Ω and using the generator matrix to calculate the calculated state vector measurements x j (k + 1) (? j, t );
Wherein the multi-target tracking method is based on received signal strength.
제 1 항에 있어서,
상기 항적 갱신단계에서,
상기 표적의 항적 갱신은 상기 표적의 데이터 연관 확률을 상기 추적 필터의 필터 게인 가중치로 입력하여 갱신하고자 하는 표적의 상태 벡터 연산치(xt(k+1)) 및 오차공분산 행렬(Pt(k+1))를 연산함으로써, 상기 생성된 표적의 항적을 갱신하는 것을 특징으로 하는 수신 신호 세기를 이용한 다중 표적 추적 방법.
The method according to claim 1,
In the navigation update step,
Wake update of the target is the filter gain state of the target desired to be updated by entering a weight vector computed value (x t (k + 1) ) and the error covariance matrix (P t (k of the tracking filter the data associated probability of the target +1)), thereby updating the generated trajectory of the target.
제 3 항에 있어서,
상기 연관 확률 생성단계에서는,
상기 도출된 발생 가능 행렬(θ)을 결합 확률 데이터 연관 기법(JPDA : Joint Probability Data-Association)에 적용하여, 상기 연산된 상태 벡터 측정치(xj(k+1))가 상기 표적에 연관될 데이터 연관 확률(βj,t)을 생성하는 것을 특징으로 하는 수신 신호 세기를 이용한 다중 표적 추적 방법.
The method of claim 3,
In the association probability generation step,
Applying the derived generated matrix? To a joint probability data association (JPDA) to determine whether the computed state vector measure x j (k + 1) And generating an association probability (? J, t ) based on the received signal strength.
제 3 항에 있어서,
상기 연관 확률 생성단계에서는,
상기 도출된 발생 가능 행렬(θ)을 단순 결합 확률 데이터 연관 기법(CJPDA : Cheap Joint Probability Data-Association)에 적용하여, 상기 연산된 상태 벡터 측정치(xj(k+1))가 상기 표적에 연관될 데이터 연관 확률(βj,t)을 생성하는 것을 특징으로 하는 수신 신호 세기를 이용한 다중 표적 추적 방법.
The method of claim 3,
In the association probability generation step,
Applying the derived generated matrix? To a simple Joint Probability Data-Association (CJPDA) so that the computed state vector measure x j (k + 1) (? J, t ) to be transmitted to the mobile station based on the received signal strength.
제 4 항에 있어서,
상기 표적의 상태 벡터 연산치(xt(k+1)) 및 오차공분산 행렬(Pt(k+1))은 수학식
Figure 112016019039665-pat00043
및 수학식
Figure 112016019039665-pat00044

(여기서, 상기 S(k+1) 는 이노베이션의 오차공분산 행렬이고, 상기 S-1(k+1)는 상기 S(k+1)의 역함수이며, 상기 W(k+1) 은 추적 필터의 필터게인으로 본 발명에서는 칼만 필터의 필터게인이고, 상기 R(k+1)은 측정잡음의 오차공분산 행렬이며, 상기 H(k+1)은 측정치와 상태 변수의 관계를 나타내는 상태변환 행렬이고, 상기 H'(k+1) 은 H(k+1) 의 역행렬이며, 상기 vt(k) 는 이노베이션의 가중합이고, 상기 β0 는 어떤 측정치와 항적이 연관되지 않는 경우의 데이터 연관 확률이며, 상기 I 는 시스템 행렬임.)에 의해 산출되는 것을 특징으로 하는 수신 신호 세기를 이용한 다중 표적 추적 방법.
5. The method of claim 4,
The state vector computed value (x t (k + 1)) and the error covariance matrix (P t (k + 1)
Figure 112016019039665-pat00043
And Equation
Figure 112016019039665-pat00044

(Here, the S (k + 1) is the error covariance matrix of the innovation, the S -1 (k + 1) is of an inverse of the S (k + 1), wherein W (k + 1) is a tracking filter In the present invention, H (k + 1) is a state transformation matrix indicating a relationship between a measured value and a state variable, and R (k + 1) is an error covariance matrix of the measured noise. (K + 1) is an inverse matrix of H (k + 1), v t (k) is a weighted sum of innovation, and β 0 is a data association probability , Where I is a system matrix). ≪ / RTI >
제 5 항에 있어서,
상기 결합 확률 데이터 연관 기법(JPDA)을 적용한 데이터 연관 확률(βj,t)은 수학식
Figure 112016019039665-pat00026
및 수학식
Figure 112016019039665-pat00027
및 수학식
Figure 112016019039665-pat00028
(여기서, 상기 c 는 정규화상수(normalized coefficent)이고, 상기 V 는 유효 영역 구간으로 πr|S(k)|1/2로 정의되며, r은 게이트 영역을 결정하는 게이트닝 임계치 상수를 의미하며, 상기 φ 는 현 상태에서 신호 개수로 여기서는 m-1 이고, 상기 N[vj,t(k+1);0;S(k+1)] 는 vj,t 의 평균이 0 이고 분산이 S(k+1) 인 가우시안 분포이며, 상기 τj(θ), δt 는 모든 발생가능행렬(θ)에서의 t 번째 표적(트랙)이 연관된 경우와 j 번째 측정치가 연관된 경우에만 1 값을 갖는 변수이고,
Figure 112016019039665-pat00029
는 t 번째 표적(트랙)의 검출 확률로, 사전에 정의된 시스템 변수이며, 상기 p1 η(at j) 는 j 번째 측정치의 신호 세기이고, 상기 at j 는 PDF(Probability density function)이며, 상기 d 는 표적 신호의 평균 신호대 잡음비를 나타내고 FMCW 레이더로부터 받은 수신 신호 세기에 비례하는 값임.)에 의해 산출되는 것을 특징으로 하는 수신 신호 세기를 이용한 다중 표적 추적 방법.
6. The method of claim 5,
The data association probability (β j, t ) applying the joint probability data association technique (JPDA)
Figure 112016019039665-pat00026
And Equation
Figure 112016019039665-pat00027
And Equation
Figure 112016019039665-pat00028
Where c is a normalized coefficient, V is defined as a valid region period, πr S S (k) 1/2 1/2 , r is a gate threshold value for determining a gate region, wherein φ is here, m-1 in number of signals in the current state, the N [v j, t (k + 1); 0; S (k + 1)] is a v j, zero mean of t and variance S (k) is a Gaussian distribution having a value of (k + 1), and τ j (θ) and δ t have a value of 1 only when the t th target (track) in all possible matrices Variable,
Figure 112016019039665-pat00029
Is a t detection probability of the first targets (track), and the system variables are defined in advance, wherein p is 1 η (a t j) is the signal strength of the j-th measurement, wherein a t j is a PDF (Probability density function), and , D is a mean signal-to-noise ratio of the target signal, and is a value proportional to a received signal strength received from the FMCW radar).
제 6 항에 있어서,
상기 단순 결합 확률 데이터 연관 기법(CJPDA)을 적용한 데이터 연관 확률(βj,t)은 수학식
Figure 112015021898211-pat00030
(여기서, 상기 λi 는 신호 대 잡음비로서
Figure 112015021898211-pat00031
이고, 상기 p1 η(aj t) 는 j 번째 추출된 표적 신호의 수신 신호 세기이며, 상기 p0 η(at j) 는 노이즈 신호의 수신 신호 세기임.)에 의해 산출되는 것을 특징으로 하는 수신 신호 세기를 이용한 다중 표적 추적 방법.
The method according to claim 6,
The data association probability (β j, t ) applying the simple joint probability data association technique (CJPDA)
Figure 112015021898211-pat00030
(Where, i is the signal-to-noise ratio
Figure 112015021898211-pat00031
And wherein p 1 η (a j t) is the received signal strength of the j-th extracted the target signal, wherein p 0 η (a t j) is being calculated by the received signal strength Im) of the noise signal Multi - target tracking method using received signal strength.
제 8 항 또는 제 9 항에 있어서,
상기 수신 신호 세기(p1 η(aj t))는 수학식
Figure 112015021898211-pat00032
및 수학식
Figure 112015021898211-pat00033
(여기서, 상기 dt 는 상기 각 표적별 신호대 잡음비(λi)를 도출하기 위하여 특정 시간 동안 누적된 수신 신호 세기의 평균값이고, 상기 L 은 신호 세기 평균 연산을 위한 누적 시간임.)에 의해 산출되는 것을 특징으로 하는 수신 신호 세기를 이용한 다중 표적 추적 방법.






10. The method according to claim 8 or 9,
The received signal strength p 1 ? (A j t )
Figure 112015021898211-pat00032
And Equation
Figure 112015021898211-pat00033
(Where d t is an average value of the received signal strengths accumulated during a specific time in order to derive a signal-to-noise ratio (? I ) for each target, and L is an accumulated time for calculating a signal intensity average) Wherein the multi-target tracking method is based on received signal strength.






KR1020150030943A 2015-03-05 2015-03-05 Multiple target tracking method using received signal strengths KR101628154B1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
KR1020150030943A KR101628154B1 (en) 2015-03-05 2015-03-05 Multiple target tracking method using received signal strengths

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020150030943A KR101628154B1 (en) 2015-03-05 2015-03-05 Multiple target tracking method using received signal strengths

Publications (1)

Publication Number Publication Date
KR101628154B1 true KR101628154B1 (en) 2016-06-08

Family

ID=56194007

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020150030943A KR101628154B1 (en) 2015-03-05 2015-03-05 Multiple target tracking method using received signal strengths

Country Status (1)

Country Link
KR (1) KR101628154B1 (en)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101886856B1 (en) * 2017-09-26 2018-08-09 한화시스템(주) System and method for data combining based on result of non-rigid object tracking on multi-sensor seeker
CN108445480A (en) * 2018-02-02 2018-08-24 重庆邮电大学 Mobile platform based on laser radar adaptively extends Target Tracking System and method
KR20200045356A (en) * 2018-10-22 2020-05-04 국방과학연구소 Method and apparatus for tracking location of target
CN112014814A (en) * 2020-08-06 2020-12-01 西安电子科技大学 Moving target detection and tracking method for screening DP-TBD based on three-frame accumulation speed
CN112101443A (en) * 2020-09-09 2020-12-18 中国航空工业集团公司雷华电子技术研究所 Measurement processing-based small-group track starting method in multi-group target scene
CN112230216A (en) * 2020-10-10 2021-01-15 南京航空航天大学 Vehicle-mounted millimeter wave radar multi-target detection method for cloud control intelligent chassis
CN112327290A (en) * 2020-10-22 2021-02-05 北京理工大学 Low-altitude flight small target tracking method based on multi-dimensional feature fusion JPDA
KR102230132B1 (en) 2020-05-08 2021-03-19 한화시스템 주식회사 Target tracking method and apparatus
CN117109588A (en) * 2023-08-25 2023-11-24 中国船舶集团有限公司第七零七研究所九江分部 Intelligent navigation-oriented multi-source detection multi-target information fusion method
CN117214881A (en) * 2023-07-21 2023-12-12 哈尔滨工程大学 Multi-target tracking method based on Transformer network in complex scene

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20020117580A1 (en) * 2001-02-01 2002-08-29 United Defense, L.P. 2-D projectile trajectory corrector
KR101045348B1 (en) * 2010-11-30 2011-06-30 엘아이지넥스원 주식회사 Target tracking method using doppler velocity in rectangular coordinate system
KR101152987B1 (en) * 2009-06-09 2012-06-04 아틀라스 엘렉트로닉 게엠베하 Direction finding method and direction finding system for detecting and tracking temporally successive bearing angles
KR101280348B1 (en) * 2012-01-18 2013-07-01 국방과학연구소 Multiple target tracking method
KR101402206B1 (en) 2014-04-10 2014-05-30 국방과학연구소 Multiple target tracking method with kinematics and feature information of targets

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20020117580A1 (en) * 2001-02-01 2002-08-29 United Defense, L.P. 2-D projectile trajectory corrector
KR101152987B1 (en) * 2009-06-09 2012-06-04 아틀라스 엘렉트로닉 게엠베하 Direction finding method and direction finding system for detecting and tracking temporally successive bearing angles
KR101045348B1 (en) * 2010-11-30 2011-06-30 엘아이지넥스원 주식회사 Target tracking method using doppler velocity in rectangular coordinate system
KR101280348B1 (en) * 2012-01-18 2013-07-01 국방과학연구소 Multiple target tracking method
KR101402206B1 (en) 2014-04-10 2014-05-30 국방과학연구소 Multiple target tracking method with kinematics and feature information of targets

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
1. Donald Lerro and Yaakov Bar-Shalom, "Automated Tracking with Target Amplitude Information," 2875-2880

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101886856B1 (en) * 2017-09-26 2018-08-09 한화시스템(주) System and method for data combining based on result of non-rigid object tracking on multi-sensor seeker
CN108445480A (en) * 2018-02-02 2018-08-24 重庆邮电大学 Mobile platform based on laser radar adaptively extends Target Tracking System and method
CN108445480B (en) * 2018-02-02 2022-05-03 重庆邮电大学 Mobile platform self-adaptive extended target tracking system and method based on laser radar
KR20200045356A (en) * 2018-10-22 2020-05-04 국방과학연구소 Method and apparatus for tracking location of target
KR102146049B1 (en) * 2018-10-22 2020-08-19 국방과학연구소 Method and apparatus for tracking location of target
KR102230132B1 (en) 2020-05-08 2021-03-19 한화시스템 주식회사 Target tracking method and apparatus
CN112014814A (en) * 2020-08-06 2020-12-01 西安电子科技大学 Moving target detection and tracking method for screening DP-TBD based on three-frame accumulation speed
CN112101443A (en) * 2020-09-09 2020-12-18 中国航空工业集团公司雷华电子技术研究所 Measurement processing-based small-group track starting method in multi-group target scene
CN112101443B (en) * 2020-09-09 2023-10-20 中国航空工业集团公司雷华电子技术研究所 Small group track starting method based on measurement processing under multi-group target scene
CN112230216A (en) * 2020-10-10 2021-01-15 南京航空航天大学 Vehicle-mounted millimeter wave radar multi-target detection method for cloud control intelligent chassis
CN112327290A (en) * 2020-10-22 2021-02-05 北京理工大学 Low-altitude flight small target tracking method based on multi-dimensional feature fusion JPDA
CN112327290B (en) * 2020-10-22 2024-02-27 北京理工大学 Low-altitude flying small target tracking method based on multidimensional feature fusion JPDA
CN117214881A (en) * 2023-07-21 2023-12-12 哈尔滨工程大学 Multi-target tracking method based on Transformer network in complex scene
CN117214881B (en) * 2023-07-21 2024-04-30 哈尔滨工程大学 Multi-target tracking method based on Transformer network in complex scene
CN117109588A (en) * 2023-08-25 2023-11-24 中国船舶集团有限公司第七零七研究所九江分部 Intelligent navigation-oriented multi-source detection multi-target information fusion method
CN117109588B (en) * 2023-08-25 2024-04-30 中国船舶集团有限公司第七零七研究所九江分部 Intelligent navigation-oriented multi-source detection multi-target information fusion method

Similar Documents

Publication Publication Date Title
KR101628154B1 (en) Multiple target tracking method using received signal strengths
US10175348B2 (en) Use of range-rate measurements in a fusion tracking system via projections
KR101402206B1 (en) Multiple target tracking method with kinematics and feature information of targets
CN105372659A (en) Road traffic monitoring multi-target detection tracking method and tracking system
CN103885057A (en) Self-adaptation variable-sliding-window multi-target tracking method
CN107576959B (en) High repetition frequency radar target tracking method before detection based on area mapping deblurring
CN111899568A (en) Bridge anti-collision early warning system, method and device and storage medium
CN108089183A (en) A kind of detecting and tracking integral method for asynchronous multi-static radar system
CN101509972A (en) Wideband radar detecting method for correcting correlation matrix based on high resolution target distance image
CN103727931A (en) Improved logic-based track initiation method
CN104199022A (en) Target modal estimation based near-space hypersonic velocity target tracking method
CN110673130B (en) Moving target track tracking method based on track association
US20200011971A1 (en) Method and apparatus for measuring distance by means of radar
CN107255805A (en) The Forecasting Methodology of radar target RCS based on weighted least-squares
CN111562570A (en) Vehicle sensing method for automatic driving based on millimeter wave radar
CN111830501B (en) HRRP history feature assisted signal fuzzy data association method and system
KR101627419B1 (en) Method for estmating location of mobile node and apparatus thereof
CN109239704A (en) A kind of adaptively sampled method based on Sequential filter interactive multi-model
JP2019138862A (en) Radar signal processor
JP4166651B2 (en) Target tracking device
CN113280821B (en) Underwater multi-target tracking method based on slope constraint and backtracking search
CN108828584A (en) Tracking before multiple frequency target detection based on track folding factor ambiguity solution
Li et al. A track initiation method for the underwater target tracking environment
CN114861725A (en) Post-processing method, device, equipment and medium for perception and tracking of target
CN113866754A (en) Moving target track correlation method based on Gaussian distribution wave gate

Legal Events

Date Code Title Description
E701 Decision to grant or registration of patent right
GRNT Written decision to grant
FPAY Annual fee payment

Payment date: 20190604

Year of fee payment: 4