Abstract:
A system and method for determination of a vehicle or crafts velocity uses on-board time dilation measurements of a moving signal generator and signal generator movement measurements to determine a velocity vector of a craft or vehicle.

Description:
BACKGROUND OF THE INVENTION 
       [0001]    1. Field of the Invention 
         [0002]    This invention relates to a system and method for determination of a vehicle&#39;s or craft&#39;s velocity using time dilation measurements made by onboard instruments. 
         [0003]    2. Description of Related Art 
         [0004]    Reliance upon external references and signals in measuring a mobile craft&#39;s motion creates a dependence on the availability, accuracy and promptness of the external links. Applications that are especially vulnerable to this dependency include, 
         [0005]    Space Probes: 
         [0006]    Transmission delays cannot be avoided. Signals from an earth based control center may be delayed by minutes. Measurements leading to emergency quick reaction corrections of the flight path are limited to those made by on-board instruments; none of these instruments measures the craft&#39;s velocity. 
         [0007]    Submarine Navigation: 
         [0008]    When submerged the submarine is isolated from all external references and signals except long wavelength radio communication, sonar and the earth&#39;s magnetic field. The measured velocity vector of the craft with respect to its surrounding medium is inaccurate, and the conversion of this velocity to the geodetic coordinate frame (the frame attached to the spinning earth) must be corrected for the ocean&#39;s current. Navigation by sonar requires a topological map of the ocean floor. 
         [0009]    Reentry and Flight Test Vehicles: 
         [0010]    When a reentry vehicle or an aircraft loses control of its flight path, the principal source of information regarding its motion is from ground based observation of its position. On board instruments measuring air speed and rate of ascent/descent are unreliable. This lack of on-board, real time measurements jeopardizes the recovery maneuvers and limits the effectiveness of the subsequent diagnostics. 
         [0011]    Velocity can be measured indirectly by an inertial navigation system (INS) containing three orthogonal accelerometers mounted on a stable table defining directions. The measured vector acceleration is time integrated producing the accumulated change in velocity. The initial velocity vector, however, must be determined by external signals or fields. Examples of this dependency are typified in the referenced patents: 
         [0012]    U.S. Pat. No. 7,890,260 by Ring [1]  describes a method to improve location accuracy by combining the velocity change obtained by integration of acceleration and the artificially generated velocity obtained from GPS. 
         [0013]    U.S. Pat. No. 7,248,967 B2 by Hagstedt [2]  describes a quasi-autonomous method in which the velocity change obtained by integration of acceleration is combined with the velocity obtained by a natural external source, the earth&#39;s magnetic field. He employs a magnetic measurement unit to register a field vector of the terrestrial magnetic field. The accuracy of the velocity depends upon the accuracy of the magnetic field map at hand. 
         [0014]    The prior art contains no techniques for measuring vector velocity independent of artificial signals and natural fields. 
       SUMMARY OF THE INVENTION 
       [0015]    Through years of operation, the Global Positioning System (GPS) has established that clocks slow with motion, that time intervals at a moving system are dilated by the velocity of the system relative to the GPS coordinate frame. This frame is three-dimensional, Euclidean, non-rotating in inertial space and anchored at the earth&#39;s barycenter (the earth-moon center of gravity). Time dilation varies with the magnitude of the GPS velocity therefore time dilation can be utilized as a measure of this velocity. 
         [0016]    Sensing velocity by time dilation is simple in concept. The sensor determines the magnitude and direction of its velocity by measuring the time dilation of a stable time interval such as the period of an atomic oscillator. This oscillator is moving at a known velocity vector that adds to the unknown velocity. As this controlled velocity vector rotates the magnitude of the sum velocity varies, reaching a maximum when the controlled and the unknown vector velocities point in the same direction, and a minimum a half turn later. The maximum and minimum of the sum velocity can be identified by the maximum and minimum of the velocity induced time dilation of the oscillator period, or the reciprocal, the minimum and maximum oscillator frequencies. The magnitude of the unknown velocity is a first order function of the maximum and minimum dilations. This measurement requires no artificial or natural external references. The measurement is in real time. 
     
    
     
       BRIEF DESCRIPTION OF THE DRAWINGS 
         [0017]    The method of sensing and measuring velocity is explained more clearly by example embodiments and with reference to the attached drawings. 
           [0018]      FIG. 1  shows the sensor geometry with basic GPS coordinates. 
           [0019]      FIG. 2  shows the arrangement of components referred to in the description of the sensor. 
           [0020]      FIG. 3  presents a flow chart of the signals and signal processing. 
       
    
    
     DETAILED DESCRIPTION OF THE INVENTION 
       [0021]    In Section 6.4 of GPS Theory and Practice by Hofmann-Wellenhof et al [1] , the authors describe the relativistic properties of the Global Positioning System (GPS), clocks must be corrected for both motion-induced and gravity-induced time dilation. The motion-induced dilation factor is drawn from Einstein&#39;s Special Theory of Relativity and the gravity-induced dilation factor is drawn from his General Theory of Relativity. To isolate the motion effects, all clocks are assumed to be at the same gravity potential 
         [0022]    Velocity is relative; the GPS reference frame is Euclidean, three-dimensional, non-rotating in inertial space, and anchored at the earth&#39;s barycenter. All GPS velocities reference this GPS coordinate frame. The GPS velocity creates a motion induced dilation factor 
         [0000]    
       
         
           
             
               γ 
               V 
             
             = 
             
               
                 1 
                 
                   
                     1 
                     - 
                     
                       
                         V 
                         2 
                       
                       / 
                       
                         c 
                         2 
                       
                     
                   
                 
               
               ≈ 
               
                 1 
                 + 
                 
                   
                     1 
                     2 
                   
                    
                   
                     
                       V 
                       2 
                     
                     
                       c 
                       3 
                     
                   
                 
               
             
           
         
       
     
         [0000]    where c is the velocity of light and V is the GPS velocity of the oscillator. Motion dilates the oscillator&#39;s period and contracts the oscillator&#39;s frequency. When stationary, the oscillator has a frequency f 0 , when moving at GPS velocity V the oscillator&#39;s frequency is contracted to f 1 =γ V   −1 f 0 , where 
         [0000]    
       
         
           
             
               γ 
               V 
               
                 - 
                 1 
               
             
             = 
             
               
                 
                   1 
                   - 
                   
                     
                       V 
                       2 
                     
                     
                       c 
                       2 
                     
                   
                 
               
               ≈ 
               
                 1 
                 - 
                 
                   
                     1 
                     2 
                   
                    
                   
                     
                       V 
                       2 
                     
                     
                       c 
                       2 
                     
                   
                 
               
             
           
         
       
     
         [0023]    The sensor measures the magnitude of its GPS vector velocity and the vector direction referencing the sensor axes. Typically these axes are, 1) GPS coordinates, 2) geodetic coordinates. or 3) body fixed coordinates. The angles measured with one reference can be converted to another system by conventional transformations. For example a GPS velocity vector is transformed to a geodetic velocity vector by the subtraction of the earth&#39;s peripheral velocity vector from the GPS vector. 
         [0024]    Whereas the current GPS operates in a region in which the earth&#39;s gravity is dominant, a GPS could be installed at any object sufficiently massive to create an extended gravity field. For example, a sun-centered GPS would have a non-rotating frame anchored at the sun&#39;s barycenter as a reference for its GPS velocities. Within the region in which the sun&#39;s gravity potential is dominant, these velocities determine the motion induced time dilation factors of moving clocks as compared to a stationary clock. 
         [0025]    In summary, the method of sensing velocity described herein measures the magnitude of the GPS velocity vector, and the vector angles referencing the sensor axes. This vector velocity can be converted to other reference frames by conventional coordinate transformations. Velocities referencing the earth-centered GPS coordinate frame dilate time in the region in which the earth&#39;s gravity is dominant. Velocities referencing other masses, such as Jupiter or the sun, dilate time in the regions in which the mass sponsored gravity field is dominant. 
       Sensor Geometry 
       [0026]    The sensor may have one, two or three modules with scan planes perpendicular to each other. In the basic embodiment (see  FIG. 1 ), the XYZ axes of the coordinate frame ( 101 ), defining the direction of the GPS vector velocities at the module, are parallel to the GPS reference frame. Typically the Z-axis is parallel to the earth&#39;s spin axis, and the XY plane is non-rotating in inertial space and parallel to the earth&#39;s equatorial plane. The measured velocity is the GPS vector velocity of the sensor. 
         [0027]    The two module embodiment, also shown in  FIG. 1 , is a typical example, The XY plane is the scan plane of the XY Module and the YZ plane is the scan plane of the YZ Module. The unknown velocity V ( 102 ) projects as V XY  ( 103 ) on the XY plane and as V YZ  ( 105 ) on the YZ plane. These projections create two right triangles. The XY Module triangle has sides V XY  and V Z  ( 104 ) and hypotenuse V, and the YZ Module triangle has sides V YZ  and V X  ( 106 ) and hypotenuse V. Not shown, the ZX Module has sides V ZX  and V Y  and hypotenuse V. 
         [0028]    Angle φ lies in the XY plane with φ=0 at the X-axis. The vector V XY  lies at angle φ V . ( 107 ). Angle θ lies in the YZ plane with θ=0 at the Y-axis. The vector V YZ  lies at angle θ V . ( 108 ). Not shown, angle ψ lies in the ZX-plane with ψ=0 at the Z-axis. The vector V ZX  lies at angle ψ V . 
         [0029]    Projected components V XY , V YZ  and V ZX  of the GPS vector velocity V can be transformed to Cartesian components by, 
         [0000]        V   X   =V   XY  cos α V   =V   ZX  sin ψ V    V   Y   =V   XY  sin α V   =V   YZ  cos β V    V   Z   =V   YZ  sin β V   =V   ZX  cos ψ V  
 
         [0030]    The scanning mechanism of the XY Module moves its oscillator at velocity ν 1  ( 109 ) where ν 1  is in the XY plane and rotates 2π radians each scan revolution. The sum velocity V XY +ν 1  is a maximum when ν 1  points in the same direction as V XY  and a minimum, when ν 1  points opposite to V XY . When V XY +ν 1  is maximum the GPS velocity V+ν 1  of the signal generator is maximum creating maximum dilation of the oscillation period, and minimum oscillation frequency. At minimum frequency, ν 1  points at angle φ V . When V XY +ν 1  is minimum the GPS velocity V+ν 1  of the signal generator is minimum creating minimum dilation of the oscillation period and maximum frequency. At maximum frequency, ν 1  points at angle φ V +π. 
         [0031]    The scanning mechanism of the YZ Module moves its oscillator at velocity ν 2  ( 110 ), where ν 2  is in the YZ plane and rotates 2π radians each scan revolution. The sum velocity V YZ +ν 2  is a maximum when ν 2  points in the same direction as V YZ  and a minimum when ν 2  points opposite to V YZ . When V YZ +ν 2  is maximum, the GPS velocity V+ν 2  of the signal generator is maximum creating maximum time dilation. At this maximum, ν 2  points at angle θ V . When V YZ +ν 2  is minimum the GPS velocity V+ν 2  of the signal generator is minimum creating minimum dilation of the oscillation period and maximum frequency. At maximum frequency, ν 2  points at angle θ V +π. 
         [0032]    The geometry of single module embodiment consists of the XY Module scanning a horizontal XY plane. The geometry of the three module embodiment includes the ZX scanning plane and projected velocity vector V ZX . 
         [0033]    When the sensor coordinate system is other than GPS, coordinate transformations are necessary to properly interpret the sensor measurements. 
         [0034]    Basic GPS Coordinate Frame: 
         [0035]    Typically the XY plane is the earth&#39;s equatorial plane. The Z-axis is the earth&#39;s spin axis. The extension of the X-axis is a fixed point on the surrounding celestial sphere. Axes of the coordinate frame defining GPS velocities at the sensor are parallel to the GPS coordinate frame axes 
         [0036]    Geodetic Coordinates: 
         [0037]    Typically the sea level earth is represented as an ellipsoid of revolution with coordinates longitude, latitude and elevation above sea level. A point at rest on the earth&#39;s surface moves at GPS velocity V PER , where the magnitude of the peripheral velocity V PER  is defined by its latitude and elevation, and its direction is due east. The geodetic horizontal plane is tangent to the ellipsoid at the sensor latitude and longitude, or parallel to this plane at sensor elevation. Typically the horizontal plane is the XY plane and is compass oriented with X pointing north. 
         [0038]    In this embodiment, the sensor measures the magnitude of its GPS velocity vector V and the direction of the velocity vector referencing its horizontal coordinate frame. The difference vector velocity V EARTH =V−V PER  is the geodetic velocity of the sensor. The horizontal component of V EARTH  is the sensor “ground velocity”; the time integral of the velocity is the ground track. The vertical component of V EARTH  is the sensor&#39;s rate of ascent/descent. 
         [0039]    Body Axes: 
         [0040]    This is the simplest installation. If the body is flexible, the two modules are mounted on a common supporting plate. If the body is rigid the modules can be separated. In one embodiment (referencing the axis designation of  FIG. 1 ) X is the roll axis, Y is the pitch axis and Z is the yaw axis. To determine the sensor&#39;s V EARTH , the body-fixed coordinates must be transformed to geodetic coordinates. This transformation requires a three-axis, gyroscopically controlled horizontal table and a directional gyro slaved to north by the earth&#39;s magnetic field or by the earth&#39;s spin. 
         [0041]    Extra-Terrestrial Systems: 
         [0042]    In the region in which a massive body&#39;s gravity is dominant, time is dilated by motion with respect to a non-rotation frame attached to the body. This body could be another planet, the sun, or the black holes at the core of our galaxy. 
       Sensor Oscillator 
       [0043]    At the heart of the sensor system, the oscillator defines the system&#39;s configuration, its feasibility and its performance. The scanning mechanism rotates at Ω radians/second with an arm length L·meters, moving the oscillator in a circle. In the XY Module, for example, the controlled vector velocity of the oscillator is ν 1  where the amplitude is ν 1 =LΩ and its direction is φ=Ωt. 
         [0044]    As shown in  FIG. 1 , the XY Module oscillator GPS velocity V OSC-1  has three orthogonal components: 1) V XY +ν 1  cos(Ωt−φ V ) pointing in the φ V  direction, 2) ν 1  sin(Ωt−φ V ) in the XY-plane and perpendicular to φ V , and 3) V Z  perpendicular to the XY-plane. Thus, 
         [0000]    
       
         
           
             
               
                 
                   
                     V 
                     
                       OSC 
                       - 
                       1 
                     
                     2 
                   
                   = 
                     
                    
                   
                     
                       
                         [ 
                         
                           
                             V 
                             XY 
                           
                           + 
                           
                             
                               v 
                               1 
                             
                              
                             
                               cos 
                                
                               
                                 ( 
                                 
                                   
                                     Ω 
                                      
                                     
                                         
                                     
                                      
                                     t 
                                   
                                   - 
                                   
                                     ϕ 
                                     V 
                                   
                                 
                                 ) 
                               
                             
                           
                         
                         ] 
                       
                       2 
                     
                     + 
                     
                       
                         [ 
                         
                           
                             v 
                             1 
                           
                            
                           
                             sin 
                              
                             
                               ( 
                               
                                 
                                   Ω 
                                    
                                   
                                       
                                   
                                    
                                   t 
                                 
                                 - 
                                 
                                   ϕ 
                                   V 
                                 
                               
                               ) 
                             
                           
                         
                         ] 
                       
                       2 
                     
                     + 
                     
                       
                         [ 
                         
                           V 
                           Z 
                         
                         ] 
                       
                       2 
                     
                   
                 
               
             
             
               
                 
                   = 
                     
                    
                   
                     
                       V 
                       XY 
                       2 
                     
                     + 
                     
                       V 
                       Z 
                       2 
                     
                     + 
                     
                       
                         v 
                         1 
                         2 
                       
                        
                       
                         [ 
                         
                           
                             
                               cos 
                               2 
                             
                              
                             
                               ( 
                               
                                 
                                   Ω 
                                    
                                   
                                       
                                   
                                    
                                   t 
                                 
                                 - 
                                 
                                   ϕ 
                                   V 
                                 
                               
                               ) 
                             
                           
                           + 
                           
                             
                               sin 
                               2 
                             
                              
                             
                               ( 
                               
                                 
                                   Ω 
                                    
                                   
                                       
                                   
                                    
                                   t 
                                 
                                 - 
                                 
                                   ϕ 
                                   V 
                                 
                               
                               ) 
                             
                           
                         
                         ] 
                       
                     
                     + 
                   
                 
               
             
             
               
                 
                     
                    
                   
                     2 
                      
                     
                       V 
                       XV 
                     
                      
                     
                       v 
                       1 
                     
                      
                     
                       cos 
                        
                       
                         ( 
                         
                           
                             Ω 
                              
                             
                                 
                             
                              
                             t 
                           
                           - 
                           
                             ϕ 
                             V 
                           
                         
                         ) 
                       
                     
                   
                 
               
             
             
               
                 
                   = 
                     
                    
                   
                     
                       V 
                       2 
                     
                     + 
                     
                       v 
                       1 
                       2 
                     
                     + 
                     
                       2 
                        
                       
                         V 
                         XV 
                       
                        
                       
                         v 
                         1 
                       
                        
                       
                         cos 
                          
                         
                           ( 
                           
                             
                               Ω 
                                
                               
                                   
                               
                                
                               t 
                             
                             - 
                             
                               ϕ 
                               V 
                             
                           
                           ) 
                         
                       
                     
                   
                 
               
             
           
         
       
     
         [0045]    Rather than focus on the oscillator period, the preference of technical literature is its reciprocal, oscillator frequency. The signal from the moving oscillator is beamed to a vehicle-stationary frequency analyzer. The reciprocal time dilation factor of the oscillator is γ OSC-1   −1 =√{square root over (1−V OSC-1   2 /c 2 )} and the reciprocal time dilation factor of the analyzer is √{square root over (1−V 2 /c 2 )}, so the signal at the analyzer is, 
         [0000]    
       
         
           
             
               
                 
                   
                     f 
                     XY 
                   
                   = 
                     
                    
                   
                     
                       
                         γ 
                         
                           OSC 
                           - 
                           1 
                         
                         
                           - 
                           1 
                         
                       
                       
                         γ 
                         V 
                         
                           - 
                           1 
                         
                       
                     
                      
                     
                       f 
                       0 
                     
                   
                 
               
             
             
               
                 
                   = 
                     
                    
                   
                     
                       f 
                       0 
                     
                      
                     
                       
                         
                           1 
                           - 
                           
                             
                               V 
                               
                                 OSC 
                                 - 
                                 1 
                               
                               2 
                             
                             / 
                             
                                 
                             
                              
                             
                               c 
                               2 
                             
                           
                         
                       
                       
                         
                           1 
                           - 
                           
                             
                               V 
                               2 
                             
                             / 
                             
                               c 
                               2 
                             
                           
                         
                       
                     
                   
                 
               
             
             
               
                 
                   ≈ 
                     
                    
                   
                     
                       f 
                       0 
                     
                      
                     
                       [ 
                       
                         1 
                         - 
                         
                           
                             1 
                             2 
                           
                            
                           
                             
                               v 
                               1 
                               2 
                             
                             
                               c 
                               3 
                             
                           
                         
                         - 
                         
                           
                             
                               
                                 V 
                                 XY 
                               
                                
                               
                                 v 
                                 1 
                               
                             
                             
                               c 
                               2 
                             
                           
                            
                           
                             cos 
                              
                             
                               ( 
                               
                                 
                                   Ω 
                                    
                                   
                                       
                                   
                                    
                                   t 
                                 
                                 - 
                                 
                                   ϕ 
                                   V 
                                 
                               
                               ) 
                             
                           
                         
                       
                       ] 
                     
                   
                 
               
             
           
         
       
     
         [0000]    The controlled velocity ν 1  is less than 10 meters/second, the term 
         [0000]    
       
         
           
             
               1 
               2 
             
              
             
               
                 v 
                 1 
                 2 
               
               
                 c 
                 2 
               
             
           
         
       
     
         [0000]    is less than 10 −15  and can be disregarded. The measured frequency is 
         [0000]    
       
         
           
             Module 
              
             
                 
             
              
             XY 
           
         
       
       
         
           
             
               
                 f 
                 XY 
               
               ≈ 
               
                 
                   f 
                   0 
                 
                 - 
                 
                   
                     f 
                     0 
                   
                    
                   
                     
                       
                         V 
                         XY 
                       
                        
                       
                         v 
                         1 
                       
                     
                     
                       c 
                       2 
                     
                   
                    
                   
                     cos 
                      
                     
                       ( 
                       
                         
                           Ω 
                            
                           
                               
                           
                            
                           t 
                         
                         - 
                         
                           ϕ 
                           V 
                         
                       
                       ) 
                     
                   
                 
               
             
             = 
             
               
                 f 
                 
                   0 
                    
                   
                       
                   
                 
               
               - 
               
                 δ 
                  
                 
                     
                 
                  
                 
                   f 
                   XY 
                 
                  
                 
                   cos 
                    
                   
                     ( 
                     
                       
                         Ω 
                          
                         
                             
                         
                          
                         t 
                       
                       - 
                       
                         ϕ 
                         V 
                       
                     
                     ) 
                   
                 
               
             
           
         
       
     
         [0000]    Controlled motion of the oscillator frequency modulates the oscillator carrier f 0  with an amplitude δf XY . Similar analyses yield, 
         [0000]    
       
         
           
             Module 
              
             
                 
             
              
             YZ 
           
         
       
       
         
           
             
               f 
               YZ 
             
             = 
             
               
                 
                   f 
                   0 
                 
                 - 
                 
                   
                     f 
                     0 
                   
                    
                   
                     
                       
                         V 
                         YZ 
                       
                        
                       
                         v 
                         2 
                       
                     
                     
                       c 
                       2 
                     
                   
                    
                   
                     cos 
                      
                     
                       ( 
                       
                         
                           Ω 
                            
                           
                               
                           
                            
                           t 
                         
                         - 
                         
                           θ 
                           V 
                         
                       
                       ) 
                     
                   
                 
               
               = 
               
                 
                   f 
                   
                     0 
                      
                     
                         
                     
                   
                 
                 - 
                 
                   δ 
                    
                   
                       
                   
                    
                   
                     f 
                     YZ 
                   
                    
                   
                     cos 
                      
                     
                       ( 
                       
                         
                           Ω 
                            
                           
                               
                           
                            
                           t 
                         
                         - 
                         
                           θ 
                           V 
                         
                       
                       ) 
                     
                   
                 
               
             
           
         
       
       
         
           
             Module 
              
             
                 
             
              
             ZX 
           
         
       
       
         
           
             
               f 
               ZX 
             
             = 
             
               
                 
                   f 
                   0 
                 
                 - 
                 
                   
                     f 
                     0 
                   
                    
                   
                     
                       
                         V 
                         ZX 
                       
                        
                       
                         v 
                         3 
                       
                     
                     
                       c 
                       2 
                     
                   
                    
                   
                     cos 
                      
                     
                       ( 
                       
                         
                           Ω 
                            
                           
                               
                           
                            
                           t 
                         
                         - 
                         
                           ψ 
                           V 
                         
                       
                       ) 
                     
                   
                 
               
               = 
               
                 
                   f 
                   
                     0 
                      
                     
                         
                     
                   
                 
                 - 
                 
                   δ 
                    
                   
                       
                   
                    
                   
                     f 
                     ZX 
                   
                    
                   
                     cos 
                      
                     
                       ( 
                       
                         
                           Ω 
                            
                           
                               
                           
                            
                           t 
                         
                         - 
                         
                           ψ 
                           V 
                         
                       
                       ) 
                     
                   
                 
               
             
           
         
       
     
         [0000]    The amplitudes of the frequency modulations are: 
         [0000]    
       
         
           
             
               δ 
                
               
                   
               
                
               
                 f 
                 XY 
               
             
             = 
             
               
                 
                   f 
                   0 
                 
                  
                 
                   
                     
                       V 
                       XY 
                     
                      
                     
                       v 
                       1 
                     
                   
                   
                     c 
                     2 
                   
                 
                  
                 
                   : 
                 
                  
                 
                     
                 
                  
                 δ 
                  
                 
                     
                 
                  
                 
                   f 
                   YZ 
                 
               
               = 
               
                 
                   
                     f 
                     0 
                   
                    
                   
                     
                       
                         V 
                         YZ 
                       
                        
                       
                         v 
                         2 
                       
                     
                     
                       c 
                       2 
                     
                   
                    
                   
                       
                   
                    
                   δ 
                    
                   
                       
                   
                    
                   
                     f 
                     ZX 
                   
                 
                 = 
                 
                   
                     f 
                     0 
                   
                    
                   
                     
                       
                         V 
                         ZX 
                       
                        
                       
                         v 
                         3 
                       
                     
                     
                       c 
                       2 
                     
                   
                 
               
             
           
         
       
     
         [0046]    Considerations in selecting the oscillator include stability and modulation measurability. In Module XY, for example, the modulation amplitude depends upon the projected velocity V XY , the controlled velocity ν 1  and the carrier frequency f 0 . Typically V XY  lies between 0 and 5000 meter/second, and ν 1  is of the order 10 meter/second: So δf XY  varies from 0 to (5×10 −13 )f 0 . The frequency modulation amplitudes for representative ultra-stable oscillators are: 
         [0000]                                                        Wavelength   Frequency   Amplitude δf           Oscillator   (meters)   (Hertz)   (Hertz)                           Rubidium    43 × 10 −3     6.8 × 10 9       0-.003           HeNe   630 × 10 −9     473 × 10 12     0-240           Mercury ion   282 × 10 −9      1.1 × 10 15     0-530                        
The rubidium microwave oscillator offers a maximum amplitude of only 3 milliradians. By contrast, the optical HeNe and mercury oscillators provide maximum amplitudes of greater than 100 Hertz.
 
       Component Arrangement 
       [0047]    The sensor modules are identical, oriented with their scan planes perpendicular to each other. The arrangement of components within a module is shown in  FIG. 2 , An oscillator ( 202 ) generates an optical signal ( 203 ) that is collimated and directed radially to a frequency analyzer ( 204 ) atop a hub ( 205 ). The oscillator is moved by a scan mechanism comprised of an arm ( 201 ) attached to the hub, a hollow shaft ( 206 ) rotating the hub, a motor ( 207 ) rotating the shaft, hub and arm. Angular position of the shaft is measured by a shaft encoder ( 208 ) and encoder electronics ( 210 ). The analyzer output signal is transferred to a computer ( 211 ) by a conduit (not shown) following the centerline of the shaft to a de-rotator ( 209 ) and thence to the computer. Alternatively, the connection is wireless. The encoder electronics connect to the computer. 
         [0048]    In one embodiment referencing the GPS coordinate frame, a gyro stabilized table (not shown) is included to transform body coordinates to GPS coordinates. In another embodiment the stable table supports the modules. 
         [0049]    In another embodiment referencing the geodetic coordinates, the body-fixed sensor system includes a horizontal reference and a directional gyro (not shown). In another embodiment the sensor axes are slaved to the geodetic horizontal plane and geodetic direction. In these embodiments the sensor may consist of one, two or three modules. 
       Signal Processing 
       [0050]    Signal generation and processing are identical in the sensor modules, but their scanning planes are orthogonal. Typically, in the three module embodiment, the XY Module the scan plane is parallel to the GPS XY plane, the YZ Module scan plane is parallel to the GPS YZ plane, and the ZX Module scan plane is parallel to the GPS ZX plane. 
         [0051]      FIG. 3  shows, as an example, XY Module processing. From its circular motion, the oscillator ( 201 ) frequency modulates its optical carrier: 
         [0000]        f   XY   =f   0   −δf   XY  cos(Ω t−φ   V ).
 
         [0000]    This optical signal is demodulated at the analyzer ( 303 ) creating a series of samples of the modulation sinusoidal envelope. This series is sent to the computer ( 305 ). In parallel the shaft encoder and its electronics ( 304 ) measures the direction of the oscillator motion, and enters this data in the computer. 
         [0052]    There are three methods of demodulation;
   1) Direct demodulation of the optical carrier. This presumes a new technology will become available.   2) Reduction of the optical carrier to a radio frequency carrier with conventional radio frequency modulation detection.   3) Reduction of the optical carrier to an audio frequency carrier with conventional audio frequency analysis.
 
The carrier frequency can be reduced by generating a beat frequency between the modulated carrier and a second reference carrier. Typically this beat is produced by an interferometer or a heterodyne mixer.
   
 
         [0056]    The first operation of the XY computer is to tag each frequency sample with the scan angle φ at the sample time. The second step is to make a best fit of a sinusoidal function, typically representing a cycle or longer, to sets of the sampled data. Next, the computer determines from the smoothed sinusoid its amplitude δf XY  and its phase φ V . 
         [0057]    The analysis is simple. The projected GPS vector velocity V XY  has an amplitude: 
         [0000]    
       
         
           
             
               V 
               XY 
             
             = 
             
               
                 
                   c 
                   2 
                 
                 
                   v 
                   1 
                 
               
                
               
                 
                   δ 
                    
                   
                       
                   
                    
                   
                     f 
                     XY 
                   
                 
                 
                   f 
                   0 
                 
               
             
           
         
       
     
         [0000]    pointing in the φ V  direction. By similar analyses, the YZ Module and ZX Module determine the magnitudes: 
         [0000]    
       
         
           
             
               V 
               ZX 
             
             = 
             
               
                 
                   c 
                   2 
                 
                 
                   v 
                   2 
                 
               
                
               
                 
                   δ 
                    
                   
                       
                   
                    
                   
                     f 
                     YZ 
                   
                 
                 
                   f 
                   0 
                 
               
                
               
                   
               
                
               pointing 
                
               
                   
               
                
               at 
                
               
                   
               
                
               
                 θ 
                 V 
               
             
           
         
       
       
         
           
             
               V 
               ZX 
             
             = 
             
               
                 
                   c 
                   2 
                 
                 
                   v 
                   3 
                 
               
                
               
                 
                   δ 
                    
                   
                       
                   
                    
                   
                     f 
                     ZX 
                   
                 
                 
                   f 
                   0 
                 
               
                
               
                   
               
                
               pointing 
                
               
                   
               
                
               at 
                
               
                   
               
                
               
                 ψ 
                 V 
               
             
           
         
       
     
         [0058]    The final operation is to transform the projected velocities into vector components of the GPS, geodetic, or body-fixed coordinate system. 
         [0059]    Although specific embodiments have been described above, any arrangement that achieves the same purpose may be substituted for the specific embodiment. This application is intended to cover any adaptation of this method of sensing velocity. The descriptions of the invention provided are not intended to limit the invention, but rather the claims define the invention. 
       Feasibility 
       [0060]    The described velocity sensor is feasible. There is a worldwide effort to develop optical atomic clocks for satellites. Any of the optical oscillators resulting from this effort would be suitable for the velocity sensor. This is a rapidly evolving technology. As reported by the National Institute for Standards and Technology [2] , a comparison of two mercury ion oscillators shows a remarkable short term fractional frequency instability of 5×10 −16 . With these or similar oscillators, proof of concept of the velocity sensor could be tested today. 
         [0061]    One of the planned space experiments is the measurement of the gravitational redshift comparing a clock aboard the International Space Station and a ground clock [3] . Both are high accuracy (10 −16 ) clocks. The velocity sensor requires a similar accuracy, but the link between it oscillators is several centimeters rather than 400 kilometers through the earth&#39;s atmosphere. The satellite velocity sensor would be amongst the first applications of he optical frequency oscillators.