[특허]Method and apparatus for detecting anomalous measurements in a satellite navigation receiver

Topcon GPS, LLC

검색어에 아래의 연산자를 사용하시면 더 정확한 검색결과를 얻을 수 있습니다.
검색연산자

검색도움말
검색연산자	기능	검색시 예
()	우선순위가 가장 높은 연산자	예1) (나노 (기계 \| machine))
공백	두 개의 검색어(식)을 모두 포함하고 있는 문서 검색	예1) (나노 기계) 예2) 나노 장영실
\|	두 개의 검색어(식) 중 하나 이상 포함하고 있는 문서 검색	예1) (줄기세포 \| 면역) 예2) 줄기세포 \| 장영실
!	NOT 이후에 있는 검색어가 포함된 문서는 제외	예1) (황금 !백금) 예2) !image
*	검색어의 *란에 0개 이상의 임의의 문자가 포함된 문서 검색	예) semi*
""	따옴표 내의 구문과 완전히 일치하는 문서만 검색	예) "Transform and Quantization"

과학기술 지식인프라 ScienceON

연합인증

연합인증 가입 기관의 연구자들은 소속기관의 인증정보(ID와 암호)를 이용해 다른 대학, 연구기관, 서비스 공급자의 다양한 온라인 자원과 연구 데이터를 이용할 수 있습니다.

이는 여행자가 자국에서 발행 받은 여권으로 세계 각국을 자유롭게 여행할 수 있는 것과 같습니다.

연합인증으로 이용이 가능한 서비스는 NTIS, DataON, Edison 등이 있습니다.

한번의 인증절차만으로 연합인증 가입 서비스에 추가 로그인 없이 이용이 가능합니다.

다만, 연합인증을 위해서는 최초 1회만 인증 절차가 필요합니다. (회원이 아닐 경우 회원 가입이 필요합니다.)

연합인증 절차는 다음과 같습니다.

최초이용시에는
ScienceON에 로그인 → 연합인증 서비스 접속 → 로그인 (본인 확인 또는 회원가입) → 서비스 이용

그 이후에는
ScienceON 로그인 → 연합인증 서비스 접속 → 서비스 이용

연합인증을 활용하시면 KISTI가 제공하는 다양한 서비스를 편리하게 이용하실 수 있습니다.

특허 상세정보

MyON담기

내보내기 페이스북 공유 트위터 공유 카카오톡 공유

Method and apparatus for detecting anomalous measurements in a satellite navigation receiver

특허상세정보
국가/구분	United States(US) Patent 등록
국제특허분류(IPC7판)	H04B-007/185
미국특허분류(USC)	342/357.03; 701/215
출원번호	US-0798536 (2004-03-11)
발명자 / 주소	Zhodzishsky, Mark I. Veitsel, Victor A. Narizhny, Andrey V.
출원인 / 주소	Topcon GPS, LLC
인용정보	피인용 횟수 : 18 인용 특허 : 6

초록
▼

A method and apparatus for detecting and correcting anomalous measurements in a satellite navigation receiver is disclosed. Anomalous measurements are detected using the redundancy of observed satellite signals, or by analyzing the relationship between phase measurements at two frequencies when using dual frequency receivers. Upon determination that an anomalous measurement exists, the particular channel on which the anomalous measurement has occurred is determined. In addition, the extent of the anomalous measurement is estimated to produce an estimated...

대표
청구항
▼

1. A method for detecting anomalous phase measurements in a satellite differential navigation system in which a first satellite receiver and a second satellite receiver each compute phase measurements on a plurality of satellite channels, said method comprising the steps, performed for each of a plurality of iterations, of:a) generating a combined phase difference vector; b) generating a phase mismatch vector representing the difference between said combined phase difference vector and an estimated combined phase difference vector; c) generating an averaged estimate vector representing the averaged estimate of said phase mismatch vector over said channels; d) generating a residual vector representing the difference between said phase mismatch vector and said averaged estimate vector; e) generating a vector of controlling signals by linear transformation of said residual vector and said averaged estimate vector; f) generating said estimated combined phase difference vector by successively storing components of said vector of controlling signals for each of said channels; and g) detecting anomalous phase measurements by analyzing said residual vector. 2. The method of claim 1 wherein said combined phase difference vector Φc comprises phase differences (Φcl . . . Φcj . . . ΦcN) for a j-th satellite channel, and wherein step a) of generating a combined phase difference vector further comprises the step of:calculating said phase difference vector Φc as Φc=ΦB?ΦR where vector ΦB comprises the full phases measured by one of said satellite receivers for each j-th satellite channel (ΦBl . . . ΦBj . . . ΦBN) and vector ΦR comprises the fall phases measured by the other satellite receiver for each j-th satellite channel (ΦRl . . . ΦRj . . . ΦRN). 3. The method of claim 1 wherein said first and second satellite receivers are dual frequency (f1 and f2) receivers and wherein said combined phase difference vector Φc comprises phase differences (Φc1 . . . Φcj . . . ΦcN) for a j-th satellite channel, and wherein step a) of generating a combined phase difference vector further comprises the step of:calculating said phase differences Φcj of said phase difference vector Φc as Φc1=(ΦB2j?ΦR2j)?(ΦB1j?ΦR1j)(f2/f1) where ΦB2j and ΦB1j represent the full phase measured by one of the satellite receivers at each j-th satellite channel at frequencies f2 and f1 respectively, andwhere ΦR2j and ΦR1j represent the full phase measured by the other of the satellite receivers at each j-th satellite channel at frequencies f2 and f1 respectively.4. The method of claim 1 wherein said step c) of generating an averaged estimate vector ΔΨ further comprises the step of:calculating said average estimate vector ΔΨ as where ΔΦcj represents the components of said phase mismatch vector for a j-th satellite channel, wj represents a weight coefficient for each j-th satellite channel, and N represents the number of satellite channels.5. The method of claim 1 wherein step c) of generating an averaged estimate vector ΔΨ further comprises the step of:calculating said average estimate vector ΔΨ as ΔΨ=H·G·ΔΦc where H is a matrix of directional cosines, ΔΦc represents said phase mismatch vector, and G is a matrix defined by G=(HTR?1H)?1HTR?1 where R is a covariance matrix of phase mismatches. 6. The method of claim 1 wherein said step e) of generating a vector of controlling signals further comprises the step of:generating said vector of controlling signals in combined loop filters of each of said plurality of channels. 7. The method of claim 6 where each of said combined loop filters generates a control signal (Ucj) for a j-th channel according to the following: where Vprj(i) is a predicted value computed on the basis of a priori data, and α, β, γ are coefficients of the loop filter.8. The method of claim 7 wherein said coefficients are time dependent and form a Kalman-type filter.9. The method of claim 7 wherein said coefficients are constants and define a third-order channel loop.10. The method of claim 1 wherein said step g) of detecting anomalous phase measurements farther comprises the steps of:h) comparing residuals of said residual vector with a threshold, and i) determining that a particular channel has an anomalous phase measurement if a residual associated with said particular channel exceeds said threshold. 11. The method of claim 10 wherein said step h) of comparing further comprises the step of comparing the absolute value of said residuals to said threshold, wherein said threshold is in the range of 0.2 to 0.25 cycles.12. The method of claim 10 wherein said step g) of detecting anomalous phase measurements further comprises the steps of:j) eliminating said particular channel from subsequent iterations of step a); and repeating steps b) through j) until no residual exceeds said threshold or until a threshold number of said channels remain. 13. The method of claim 1 wherein said step g) of detecting anomalous phase measurements further comprises the steps of:h) generating a weighted sum of residual squares for said channels; i) comparing said weighted sum of residual squares to a threshold; and j) determining that an anomalous phase measurement exists if said weighted sum of residual squares exceeds said threshold. 14. The method of claim 13 wherein said step g) of detecting anomalous phase measurements further comprises the step of:k) determining that a particular channel has an anomalous phase measurement if a residual associated with said particular channel is the maximum residual in said residual vector. 15. The method of claim 14 wherein said step g) of detecting anomalous phase measurements further comprises the steps of:l) eliminating said particular channel from subsequent iterations of step a); and repeating steps b) through l) until said weighted sum of residual squares does not exceed said threshold or until a threshold number of said channels remain. 16. The method of claim 13 wherein said step g) of detecting anomalous phase measurements further comprises the steps of:k) eliminating x channels from subsequent iterations of step a), varying the channels eliminated for each iteration of step k); and repeating steps b) through k) until either said weighted sum of residual squares does not exceed said threshold or until x reaches a threshold, incrementing x by one after all combinations of channels for a given x have been eliminated during an iteration of step k). 17. The method of claim 1 wherein said steps are performed in one of said satellite receivers.18. The method of claim 17 wherein said one satellite receiver receives phase measurements from the other satellite receiver.19. The method of claim 1 further comprising the step of:h) using information about said anomalous phase measurements detected in step g) in a navigation task. 20. The method of claim 19 further comprising the step of:i) generating a cycle slip estimate for a channel determined to have an anomalous phase measurement and using said cycle slip estimate in said navigation task. 21. An apparatus for detecting anomalous phase measurements in a satellite differential navigation system in which a first satellite receiver and a second satellite receiver each compute phase measurements on a plurality of satellite channels, said apparatus comprising:a combined phase difference generator for generating a combined phase difference vector; means for generating a phase mismatch vector representing the difference between said combined phase difference vector and an estimated combined phase difference vector; an integrated discriminator for generating an averaged estimate vector representing the averaged estimate of said phase mismatch vector over said channels; means for generating a residual vector representing the difference between said phase mismatch vector and said averaged estimate vector; at least one joint loop filter for generating a vector of controlling signals by linear transformation of said residual vector and said averaged estimate vector; an accumulator for generating said estimated combined phase difference vector by successively storing components of said vector of controlling signals for each of said channels; and a residuals analyzer for detecting anomalous phase measurements by analyzing said residual vector. 22. The apparatus of claim 21 wherein said combined phase difference vector Φc comprises phase differences (Φcl . . . Φcj . . . ΦcN) for a j-th satellite channel, and wherein said phase difference generator generates said combined phase difference vector Φc as Φc=ΦB?ΦR where vector ΦB comprises the full phases measured by one of said satellite receivers for each j-th satellite channel (ΦBl . . . ΦBj . . . ΦBN) and vector ΦR comprises the full phases measured by the other satellite receiver for each j-th satellite channel (ΦRl . . . ΦRj . . . ΦRN).23. The apparatus of claim 21 wherein said first and second satellite receivers are dual frequency (f1 and f2) receivers and wherein said combined phase difference vector Φc comprises phase differences (Φcl . . . Φcj . . . ΦcN) for a j-th satellite channel, and wherein said phase difference generator generates said phase differences Φcj of said phase difference vector Φc asΦcj=(ΦB2j?ΦR2j)?(ΦB1j?ΦR1j)(f2/f2) where ΦB2j and ΦB1j represent the full phase measured by one of the satellite receivers at each j-th satellite channel at frequencies f2 and f1 respectively, andwhere ΦR2j and ΦR1j represent the full phase measured by the other of the satellite receivers at each j-th satellite channel at frequencies f2 and f1 respectively.24. The apparatus of claim 21 wherein said integrated discriminator generates said averaged estimate vector ΔΨ as where ΔΦcj represents the components of said phase mismatch vector for a j-th satellite channel, wj represents a weight coefficient for each j-th satellite channel, and N represents the number of satellite channels.25. The apparatus of claim 21 wherein said integrated discriminator generates said averaged estimate vector ΔΨ as ΔΨ=H·G·ΔΦc where H is a matrix of directional cosines, ΔΦc represents said phase mismatch vector, and G is a matrix defined by G=(HTR?1H)?1HTR?1 where R is a covariance matrix of phase mismatches.26. The apparatus of claim 21 wherein said joint loop filter generates said vector of controlling signals in combined loop filters of each of said plurality of channels.27. The apparatus of claim 26 wherein each of said at least one combined loop filters generates a control signal (Ucj) for a j-th channel according to the following: where Vprj(i) is a predicted value computed on the basis of a priori data, and α, β, γ are coefficients of the loop filter.28. The apparatus of claim 27 wherein said coefficients are time dependent and form a Kalman-type filter.29. The apparatus of claim 27 wherein said coefficients are constants and define a third-order channel loop.30. The apparatus of claim 21 wherein said residuals analyzer detects anomalous phase measurements by comparing residuals of said residual vector with a threshold and determining that a particular channel has an anomalous phase measurement if a residual associated with said particular channel exceeds said threshold.31. The apparatus of claim 21 wherein said residuals analyzer detects anomalous phase measurements by generating a weighted sum of residual squares for said channels, comparing said weighted sum of residual squares to a threshold, and determining that an anomalous phase measurement exists if said weighted sum of residual squares exceeds said threshold.32. The apparatus of claim 31 wherein said residuals analyzer further detects anomalous phase measurements by determining that a particular channel has an anomalous phase measurement if a residual associated with said particular channel is the maximum residual in said residual vector.33. A method for detecting anomalous phase measurements in a satellite differential navigation system in which a first satellite receiver and a second satellite receiver each compute phase measurements on a plurality of satellite channels, said method comprising the steps, performed for each of a plurality of iterations, of:a) generating a combined phase difference vector; b) generating an increment vector representing the difference between said combined phase difference vector and a combined phase difference vector of a preceding iteration, c) generating an averaged increment estimate vector representing the averaged estimate of said increment vector; d) generating an incremental residual vector representing the difference between said increment vector and said averaged increment estimate vector; e) generating an integrated residual vector from said incremental residual vector; and f) detecting anomalous phase measurements by analyzing said integrated residual vector. 34. The method of claim 33 wherein step a) of generating a combined phase difference vector Φc(i) further comprises the step of:calculating said phase difference vector Φc(i) as Φc(i)=ΦB(i)?ΦR(i)?ΦBR(i) where vector ΦB(i) comprises the full phases measured by one of said satellite receivers for each i-th iteration and vector ΦR(i) comprises the full phases measured by the other satellite receiver for each i-th iteration and ΦBR(i) is a prediction of satellite movement. 35. The method of claim 33 wherein said step c) of generating an averaged increment estimate vector further comprises the step of:calculating said average estimate vector Δ{circumflex over (Φ)} as Δ{circumflex over (Φ)}(i)=H(i)·G(i)·ΔΦc(i), where Δ{circumflex over (Φ)}(i) represents the components of said average estimate vector for an i-th iteration, H is a matrix of directional cosines, G is a matrix defined by G=(HTR?1H)?1HTR?1 where R is a covariance matrix of phase mismatches, and ΔΦc(i)=Φc(i)?Φc(i?1) where Φc(i) is a vector of combined phase difference at each i-th iteration.36. The method of claim 33 wherein said step e) of generating an integrated residual vector further comprises the step of generating integrated residuals δ(i) asδ(i)=A·δ(i?1)+Δδ(i) for each i-th iteration, where A is a coefficient in the range 0.995 . . . 0.999.37. The method of claim 33 wherein said step f) of detecting anomalous phase measurements further comprises the steps of:g) generating a weighted sum of integrated residual squares for said channels; h) comparing said weighted sum of integrated residual squares to a threshold; and i) determining that an anomalous phase measurement exists if said weighted sum of integrated residual squares exceeds said threshold. 38. The method of claim 37 wherein said step f) of detecting anomalous phase measurements further comprises the steps of:j) eliminating x channels from subsequent iterations of step a), varying the channels eliminated for each iteration of step j); and repeating steps b) through j) until either said weighted sum of integrated residual squares does not exceed said threshold or until x reaches a threshold, incrementing x by one after all combinations of channels for a given x have been eliminated during an iteration of step j). 39. The method of claim 33 wherein said steps are performed in one of said satellite receivers.40. The method of claim 39 wherein said one satellite receiver receives phase measurements from the other satellite receiver.41. The method of claim 33 further comprising the step of:g) using information about said anomalous phase measurements detected in step f) in a navigation task. 42. The method of claim 41 further comprising the step of:h) generating a cycle slip estimate for a channel determined to have an anomalous phase measurement and using said cycle slip estimate in said navigation task. 43. The method of claim 33 further comprising the steps of:g) receiving channel indicator alarms from a channel indicator, said channel indicator alarms marking channels with anomalous phase measurements; h) eliminating channels marked with channel indicator alarms from subsequent iterations of step a); and i) repeating steps a) through f). 44. The method of claim 43 wherein said step f) of detecting anomalous phase measurements further comprises the steps of:j) generating a weighted sum of integrated residual squares for said channels; k) comparing said weighted sum of integrated residual squares to a threshold; and l) determining that an anomalous phase measurement exists if said weighted sum of integrated residual squares exceeds said threshold. 45. An apparatus for detecting anomalous phase measurements in a satellite differential navigation system in which a first satellite receiver and a second satellite receiver each compute phase measurements on a plurality of satellite channels, said apparatus comprising:a) a combined phase difference generator for generating a combined phase difference vector; b) means for generating an increment vector representing the difference between said combined phase difference vector and a combined phase difference vector of a preceding measurement, c) an integrated converter for generating an averaged increment estimate vector representing the averaged estimate of said increment vector; d) means for generating an incremental residual vector representing the difference between said increment vector and said averaged increment estimate vector, e) a digital filter for generating an integrated residual vector from said incremental residual vector; and f) a residuals analyzer for detecting anomalous phase measurements by analyzing said integrated residual vector. 46. The apparatus of claim 45 wherein said combined phase difference generator generates said combined phase difference vector Φc(i) as Φc(i)=ΦB(i)?ΦR(i)?ΦBR(i) where vector ΦB(i) comprises the full phases measured by one of said satellite receivers for each i-th iteration and vector ΦR(i) comprises the full phases measured by the other satellite receiver for each i-th iteration and ΦBR(i) is a prediction of satellite movement.47. The apparatus of claim 45 wherein said integrated converter generates said averaged increment estimate vector Δ{circumflex over (Φ)} asΔ{circumflex over (Φ)}(i)=H(i)·G(i)·ΔΦc(i), where Δ{circumflex over (Φ)}(i) represents the components of said average estimate vector for an i-th iteration, H is a matrix of directional cosines, G is a matrix defined by G=(HTR?1H)?1HTR?1 where R is a covariance matrix of phase mismatches, and ΔΦc(i)=Φc(i)?Φc(i?1) where Φc(i) is a vector of combined phase difference at each i-th iteration.48. The apparatus of claim 45 wherein said residuals analyzer detects anomalous phase measurements by generating a weighted sum of integrated residual squares for said channels, comparing said weighted sum of integrated residual squares to a threshold, and determining that an anomalous phase measurement exists if said weighted sum of integrated residual squares exceeds said threshold.49. The apparatus of claim 45 wherein said digital filter generates integrated residuals δ(i) asδ(i)=A·δ(i?1)+Δδ(i) for each of an i-th iteration, where A is a coefficient in the range 0.995 . . . 0.999.50. A method for detecting anomalous phase measurements in a satellite differential navigation system in which a first satellite receiver and a second satellite receiver each compute phase measurements on a plurality of satellite channels, said method comprising the steps, performed for each of a plurality of iterations, of:a) generating a combined phase difference vector; b) generating an increment vector representing the difference between said combined phase difference vector and a combined phase difference vector of an initial measurement at an initial time, c) generating a corrected increment vector using a cycle slip correction estimate generated in a preceding iteration; d) analyzing channel indicator alarms; e) generating an averaged increment estimate vector representing the averaged estimate of said corrected increment vector using channels not associated with a channel indicator alarm; f) generating a residual vector representing the difference between the corrected increment vector and the averaged increment estimate vector; and g) generating said cycle slip correction estimate using said residual vector. 51. The method of claim 50 wherein said step a) of generating a combined phase difference vector further comprises the step of:including a prediction of phase difference in said difference. 52. The method of claim 50 wherein said initial time is reset periodically.53. The method of claim 50 further comprising the step of:h) using said cycle slip correction estimate in a navigation task. 54. An apparatus for detecting anomalous phase measurements in a satellite differential navigation system in which a first satellite receiver and a second satellite receiver each compute phase measurements on a plurality of satellite channels, said apparatus comprising:a) a combined phase difference generator for generating a combined phase difference vector; b) means for generating an increment vector representing the difference between said combined phase difference vector and a combined phase difference vector of an initial measurement at an initial time; c) a cycle slip correction unit for generating a cycle slip correction estimate; d) means for generating a corrected increment vector using a cycle slip correction estimate generated in a preceding measurement; d) a channel indicator analyzer for analyzing channel indicator alarms; e) an integrated converter for generating an averaged increment estimate vector representing the averaged estimate of said corrected increment vector using channels not associated with a channel indicator alarm; f) means for generating a residual vector representing the difference between the corrected increment vector and the averaged increment estimate vector; and g) a correction unit for generating said cycle slip correction estimates using said residual vector. 55. The apparatus of claim 54 wherein said combined phase difference generator further generates said combined phase difference vector by including a prediction of phase difference in said difference.56. The apparatus of claim 54 wherein said initial time is reset periodically.

인용문헌 (6)

이 특허를 인용한 특허 (18)

IPC 상위 출원인

내보내기

상세정보
내보내기 구분	파일저장 인쇄 메일전송
구성항목	관리번호, 국가코드, 자료구분, 상태, 출원번호, 출원일자, 공개번호, 공개일자, 등록번호, 등록일자, 발명명칭(한글), 발명명칭(영문), 출원인(한글), 출원인(영문), 출원인코드, 대표IPC 관리번호, 국가코드, 자료구분, 상태, 출원번호, 출원일자, 공개번호, 공개일자, 공고번호, 공고일자, 등록번호, 등록일자, 발명명칭(한글), 발명명칭(영문), 출원인(한글), 출원인(영문), 출원인코드, 대표출원인, 출원인국적, 출원인주소, 발명자, 발명자E, 발명자코드, 발명자주소, 발명자 우편번호, 발명자국적, 대표IPC, IPC코드, 요약, 미국특허분류, 대리인주소, 대리인코드, 대리인(한글), 대리인(영문), 국제공개일자, 국제공개번호, 국제출원일자, 국제출원번호, 우선권, 우선권주장일, 우선권국가, 우선권출원번호, 원출원일자, 원출원번호, 지정국, Citing Patents, Cited Patents
저장형식	Text(ASCII format) Excel format PIAS분석(.xls)
메일정보	받는사람 (필수) @ 보내는사람 (선택) @ 제목 내용 KISTI 검색결과 이메일 서비스
안내	> 총 건의 자료가 검색되었습니다. > 다운받으실 자료의 인덱스를 입력하세요. (1-10,000) > 검색결과의 순서대로 최대 10,000건 까지 다운로드가 가능합니다. > 데이타가 많을 경우 속도가 느려질 수 있습니다.(최대 2~3분 소요) > 다운로드 파일은 UTF-8 형태로 저장됩니다. 파일의 내용이 제대로 보이지 않을실 때는 웹브라우저 상단의 보기 -> 인코딩 -> 자동선택 여부를 확인하십시오 > ~ > Text(ASCII format) Excel format

과학기술 지식인프라 ScienceON

ScienceON Chatbot

연합인증

특허 상세정보

Method and apparatus for detecting anomalous measurements in a satellite navigation receiver

초록
▼

대표
청구항
▼

인용문헌 (6)

이 특허를 인용한 특허 (18)

IPC 상위 출원인

과제정보

상세정보조회

원문조회

IPC 상세정보(version : 201201)

연구과제 성과물(0)

참여연구원의 다른 문헌(0)

처리중입니다.

잠시만 기다려주세요.

ScienceON Chatbot

연합인증

특허 상세정보

Method and apparatus for detecting anomalous measurements in a satellite navigation receiver

초록 ▼

대표청구항 ▼

인용문헌 (6)

이 특허를 인용한 특허 (18)

IPC 상위 출원인

과제정보

상세정보조회

원문조회

IPC 상세정보(version : 201201)

연구과제 성과물(0)

참여연구원의 다른 문헌(0)

처리중입니다.

잠시만 기다려주세요.

초록
▼

대표
청구항
▼