×

Method and apparatus for finding aircraft position by integrating accelerations less time averages

  • US 5,590,044 A
  • Filed: 02/02/1995
  • Issued: 12/31/1996
  • Est. Priority Date: 02/02/1994
  • Status: Expired due to Fees
First Claim
Patent Images

1. A method for determining the position of an aircraft for radar mapping, comprising:

  • providing a processor in the airplane;

    providing a master navigation system coupled to the processor and including an aircraft speedometer for measuring a forward x-speed vXUAX (t), a lateral horizontal position sensor for measuring a lateral y-position pYUAX (t), and an altimeter for measuring a vertical z-position pZUAX (t);

    providing a slave system including an inertial measurement unit coupled to the processor and further including an x-accelerometer for generating a forward-direction horizontal acceleration signal ax (t), a y-accelerometer for generating a lateral horizontal acceleration signal ay (t), and a z-accelerometer for generating a vertical acceleration signal ax (t);

    (a) inputting the forward-direction acceleration signal ax (t) from the inertial system to the processor;

    time-averaging over a time interval T the acceleration signal ax (t) to obtain a first x-adjustment ax (t);

    subtracting the first x-adjustment ax (t) from the acceleration signal ax (t) to obtain an adjusted acceleration signal a'"'"'x (t);

    integrating the adjusted acceleration signal a'"'"'x (t) over the interval T to obtain a velocity vx "(t) according to ##EQU10## time-averaging the velocity v"x (t) over the interval T to obtain a vx -adjustment vx (t);

    subtracting the vx -adjustment vx (t) from the velocity v"x (t) to obtain a forward airplane speed vXIMU (t) in the desired flight direction according to
    
    
    space="preserve" listing-type="equation">ν

    .sub.XIMU (t)=ν

    .sub.X (t)-ν

    ".sub.X (t) (5)estimating by method of least squares, from a difference between the speed vXIMU (t) and the speed measurement of the speedometer, an average acceleration aX0 and an average speed vX0 ;

    creating a correction term vXFIT (t) from a sum of the average speed vX0 and a time integration of aX0 according to ##EQU11## correcting the speed vXIMU (t) in the desired flight direction by adding the correction term vXFIT (t) thereto to obtain a corrected flight speed vXKORR (t) according to
    
    
    space="preserve" listing-type="equation">ν

    .sub.XKORR (t)=ν

    .sub.XIMU (t)+ν

    .sub.XFIT (t)(b) inputting the lateral-direction acceleration signal ay (t) from the inertial system to the processor;

    time-averaging over a time interval T the acceleration signal ay (t) to obtain a first y-adjustment ay (t);

    subtracting the first y-adjustment ay (t) from the acceleration signal ay (t) to obtain an adjusted acceleration signal a'"'"'y (t);

    integrating the adjusted acceleration signal a'"'"'y (t) over the interval T to obtain a velocity vy "(t) according to ##EQU12## time-averaging the velocity v"y (t) over the interval T to obtain a vy -adjustment vy (t);

    subtracting the first y-adjustment vy (t) from the velocity v"y (t) to obtain an adjusted velocity signal v'"'"'y (t);

    integrating the adjusted velocity signal v'"'"'y (t) over the interval T to obtain a position py '"'"'(t) according to ##EQU13## time-averaging the position p'"'"'y (t) over the interval T to obtain a py -adjustment py (t);

    subtracting the py -adjustment py (t) from the velocity p'"'"'y (t) to obtain a forward airplane speed pYIMU (t) in the desired flight direction according to
    
    
    space="preserve" listing-type="equation">p.sub.YIMU (t)=p'"'"'.sub.Y (t)-p'"'"'.sub.Y (t) (9)estimating by method of least squares, from a difference between the speed pYIMU (t) and the speed measurement of the speedometer, an average acceleration aY0 and an average speed vY0 ;

    creating a correction term pYFIT (t) from a sum of the average speed vY0 and a time integration of aY0 according to ##EQU14## correcting the speed pYIMU (t) in the desired flight direction by adding the correction term pYFIT (t) thereto to obtain a corrected flight speed pYKORR (t) according to
    
    
    space="preserve" listing-type="equation">p.sub.YKORR (t)=p.sub.YIMU (t)+p.sub.YFIT (t)correcting the motion of the aircraft along the desired flight path according to the value of pYKORR (t);

    (c) inputting the lateral-direction acceleration signal az (t) from the inertial system to the processor;

    time-averaging over a time interval T the acceleration signal az (t) to obtain a first z-adjustment az (t);

    subtracting the first z-adjustment az (t) from the acceleration signal az (t) to obtain an adjusted acceleration signal a'"'"'z (t);

    integrating the adjusted acceleration signal a'"'"'z (t) over the interval T to obtain a velocity vz "(t) according to ##EQU15## time-averaging the velocity v"z (t) over the interval T to obtain a vz -adjustment vz (t);

    subtracting the first z-adjustment vz (t) from the velocity vz (t) to obtain an adjusted velocity signal v'"'"'z (t);

    integrating the adjusted velocity signal v'"'"'z (t) over the interval T to obtain a position pz '"'"'(t) according to ##EQU16## time-averaging the position p'"'"'z (t) over the interval T to obtain a pz -adjustment pz (t);

    subtracting the pz -adjustment pz (t) from the velocity p"z (t) to obtain a forward airplane speed pZIMU (t) in the desired flight direction according to
    
    
    space="preserve" listing-type="equation">p.sub.ZIMU (t)=p'"'"'.sub.z (t)-p'"'"'.sub.z (t) (15)estimating by method of least squares, from a difference between the speed pZIMU (t) and the speed measurement of the speedometer, an average acceleration aZ0 and an average speed vZ0 ;

    creating a correction term pZFIT (t) from a sum of the average speed vZ0 and a time integration of aZ0 according to ##EQU17## correcting the speed pZIMU (t) in the desired flight direction by adding the correction term pZFIT (t) thereto to obtain a corrected flight speed pZKORR (t) according to
    
    
    space="preserve" listing-type="equation">p.sub.ZKORR (t)=p.sub.ZIMU (t)+p.sub.ZFIT (t) (29)andcompensating the radar mapping according to the values of vXKORR (t), pYKORR (t), and pZKORR (t).

View all claims
  • 1 Assignment
Timeline View
Assignment View
    ×
    ×