Reduced state estimation with biased measurements
원문보기
IPC분류정보
국가/구분
United States(US) Patent
등록
국제특허분류(IPC7판)
G01S-013/66
G01S-013/00
출원번호
US-0347974
(2006-02-06)
등록번호
US-7277047
(2007-10-02)
발명자
/ 주소
Mookerjee,Purusottam
Reifler,Frank J.
출원인 / 주소
Lockheed Martin Corporation
대리인 / 주소
Duane Morris, LLP
인용정보
피인용 횟수 :
9인용 특허 :
6
초록▼
This invention relates to sate estimation after processing measurements with unknown biases that may vary arbitrarily in time within known physical bounds. These biased measurements are obtained from systems characterized by state variables and by multidimensional parameters, for which the latter ar
This invention relates to sate estimation after processing measurements with unknown biases that may vary arbitrarily in time within known physical bounds. These biased measurements are obtained from systems characterized by state variables and by multidimensional parameters, for which the latter are also known and may vary arbitrarily in time within known physical bounds. The measurements are processed by a filter using a mean square optimization criterion that accounts for random and biased measurement errors, as well as parameters excursions, to produce estimates of the true states of the system. The estimates are applied to one of (a) making a decision, (b) operating a control system, and (c) controlling a process.
대표청구항▼
What is claimed is: 1. A method for recursively estimating the state of a system having multidimensional parameters λ in addition to state variables x(k) at time ik for k=0,1,2, . . . , which parameters λ are unknown, arbitrarily time-varying, but bounded, and driven by the input function
What is claimed is: 1. A method for recursively estimating the state of a system having multidimensional parameters λ in addition to state variables x(k) at time ik for k=0,1,2, . . . , which parameters λ are unknown, arbitrarily time-varying, but bounded, and driven by the input function u(x(k),λ), which may be nonlinear, and expressed by the state equation description="In-line Formulae" end="lead"x(k+1)=Φx(k)+Γu( x(k),λ) (38)description="In-line Formulae" end="tail" where Φ,Γ are system matrices dependent on the discrete time interval T=tk+1-tk, said method comprising the following steps: measuring aspects of the state of the system to produce initial measurements expressed by the measurement equation description="In-line Formulae" end="lead"z(k)=Hx(k)+Jb+n(k) (39)description="In-line Formulae" end="tail" for 1≦k≦k0, where, if no measurements are used in the initialization of the filter, k0=0, where b is an unknown arbitrarily time-varying, but bounded, measurement bias vector with covariance B, whose components correspond to the different sensors, and where the sensor selector matrix J selects the appropriate components of sensor bias, and where n(k) is the measurement noise with covariance N and measurement matrix H at time lk; initializing state estimates {circumflex over (x)}(k0|k0) and the matrices M(k0|k0), D(k0|k0), E(k0|k0) using a priori information and the initial measurements, where vector {circumflex over (x)}(j|k) is defined as the estimate of the state of the system at time tj for j=0,1,2, . . after processing k measurements z(i) for 1≦i≦k; matrix M(j|k) is defined as the covariance of the state estimation errors at time lj for j=0,1,2, . . . due only to the random errors in the k measurements r(f) for 1≦i≦k and a priori initial information that is independent of the parameter uncertainty and measurement bias uncertainty; matrix D(j|k) is defined as the matrix of bias coefficients, which linearly relates state estimation errors to the parameter errors, at time 1j for j=0,1,2, . . . after processing k measurements x(t) for 1≦i≦k; matrix E(j|k) is defined as the matrix of bias coefficients, which linearly relates state estimation errors to the sensor measurement bias, at time tj for j=0,1,2, . . . after processing k measurements z(l) for 1≦i≦k; measuring aspects z(k) of the state of the system expressed by the measurement equation description="In-line Formulae" end="lead"z(k)=Hx(k)+Jb+n(k) (40)description="In-line Formulae" end="tail" where b is an unknown arbitrarily time-varying, but bounded, measurement bias vector with covariance B, whose components correspond to the different sensors, and where the sensor selector matrix J selects the appropriate components of sensor bias, and where n(k) is the measurement noise with covariance N and measurement matrix H at time ik for k≦k0; determining the system transition matrices Φ and Γ using the update interval T=tk+1-tk when a new measurement x(k+1) arrives at time tk+1; determining the mean value λ of unknown but bounded parameters λ, and the input vector u({circumflex over (x)}(k|k), λ); determining F,G using generating a parameter matrix Λ, representing physical bounds on those parameters that are not state variables of the system; extrapolating said state estimates {circumflex over (x)}(k|k) and matrices M(k|k), D(k|k), E(k|k), S(k|k) to {circumflex over (x)}(k+1|k), M(k+1|k), D(k+1|k), E(k+1|k) and S(k+1|k) as description="In-line Formulae" end="lead"{circumflex over (x)}(k+1|k)=Φ{circumflex over (x)})k|k)+Γu({circumflex over (x)}(k|k), λ) (43)description="In-line Formulae" end="tail" description="In-line Formulae" end="lead"M(k+1|k)=FM(k|k)F' (44)description="In-line Formulae" end="tail" description="In-line Formulae" end="lead"D(k+1|k)=FD(k|k)+G (45)description="In-line Formulae" end="tail" description="In-line Formulae" end="lead"E(k+1|k)=FE(k|k) (46)description="In-line Formulae" end="tail" description="In-line Formulae" end="lead"S(k+1|k)=M(k+1|k)+ D(k+1|k(ΛD)k+1|k)'E (k+1|k)BE(k+1|k)' (47)description="In-line Formulae" end="tail" determining covariance of the residual Q as description="In-line Formulae" end="lead"V=HE(k+1|k)+J (48)description="In-line Formulae" end="tail" description="In-line Formulae" end="lead"Q=H[M(k+1|k)+D(k+1|k )ΛD)k+1|k)']H'+VBV+N (49)description="In-line Formulae" end="tail" determining the filter gain matrix K as description="In-line Formulae" end="lead"A=S(k+1|k)H'+E(k+1|k )BJ' (50)description="In-line Formulae" end="tail" description="In-line Formulae" end="lead"K=AQ-1 (51)description="In-line Formulae" end="tail" determining the matrix L as description="In-line Formulae" end="lead"L=1-KH (52)description="In-line Formulae" end="tail" where I is the identity matrix; updating the state estimate {circumflex over (x)}(k+1|k) as description="In-line Formulae" end="lead"{circumflex over (x)}(k+1|k+1)={circumflex over (x)})k+1|k)+K[z)k+1)-H{circumflex over (x)})k+1|k) (53)description="In-line Formulae" end="tail" updating the matrices M(k+1|k), D(k+1|k), E(k+1|k) as description="In-line Formulae" end="lead"M(k+1|k+1)=LM(k+1|k) L'+KNK' (54)description="In-line Formulae" end="tail" description="In-line Formulae" end="lead"D(k+1|k+1)=LD(k+1|k) (55)description="In-line Formulae" end="tail" description="In-line Formulae" end="lead"E(k+1|k+1)=LE)k+1|k)-KJ (56)description="In-line Formulae" end="tail" respectively, and generating the total mean square error S(k+1|k+1) as description="In-line Formulae" end="lead"S(k+1|k+1)=M)k+1|k +1)+D(k+1|k+1)ΛD(k+1|k+1) '+E(k+1|k+1)BE(k+1|k+1)' (57).description="In-line Formulae" end="tail" 2. A method for estimating the state of a system comprising the steps of: observing a system having state variables and also having unknown, multidimensional, arbitrarily time-varying parameters, but which are subject to known bounded values; measuring certain aspects of the state of the system in the presence of sensor measurement biases and random errors to produce initial measurements; initializing state estimates and matrices using a priori information and the initial measurements; using the update interval in determining the system transition matrices and the mean value of unknown but bounded parameters and the input vector; applying the measurements to an estimating filter that explicitly uses a mean square optimization criterion that separately accounts for measurement biases and errors and said bounding values, to produce estimates of the true state of the system; and applying said estimates to one of (a) make a decision (b) operate a control system, and (c) control a process.
연구과제 타임라인
LOADING...
LOADING...
LOADING...
LOADING...
LOADING...
이 특허에 인용된 특허 (6)
Blair William D. (Fredericksburg VA) Watson Gregory A. (Fredericksburg VA) Rice Theodore R. (Fredericksburg VA), Interacting multiple bias model filter system for tracking maneuvering targets.
Zemany Paul D. (Amherst NH) Smith Kirby A. (Derry NH) Smith John A. (Bedford MA), State-deviation-estimation circuit employing a phase-locked-loop phase reference.
Millet, Nicolas; Allam, Sébastien; Klein, Mathieu; Malherbe, Thierry, Multi-target data processing for multi-receiver passive radars in an SFN or MFN mode.
Mookerjee, Purusottam; Reifler, Frank J., Track fusion by optimal reduced state estimation in multi-sensor environment with limited-bandwidth communication path.
※ AI-Helper는 부적절한 답변을 할 수 있습니다.