Method and system for determining a position of a communication satellite utilizing twoway ranging

0Associated
Cases 
0Associated
Defendants 
0Accused
Products 
124Forward
Citations 
0
Petitions 
1
Assignment
First Claim
1. A method for determining a position of an orbiting transceiver in a communications network including at least a first and second transceiver at a first and second known location, respectively, on Earth, the first and second transceivers adapted to transmit and receive communications signals to and from the orbiting transceiver, the method comprising:
 determining a first and second range measurement between each of the first and second transceivers and the orbiting transceiver, respectively;
determining a first and second range rate corresponding to a time rate of change of the first and second range measurements, respectively;
determining a circle of intersection representative of the set of possible positions for the orbiting transceiver based on the first range measurement and the first range rate, the circle of intersection having a specific orientation in space, a specific radius and a center in a specific, three dimensional position relative to the known position of the first transceiver;
determining an angular position of the orbiting transceiver along the circle of intersection based on the known position of the second transceiver and the second range measurement; and
determining the position of the orbiting transceiver based on the circle of intersection and the angular position.
1 Assignment
0 Petitions
Accused Products
Abstract
A method and system for determining a position of an orbiting transceiver in a communications network includes at least a first and second transceiver at a first and second known location, respectively, on Earth. The first and second transceivers are adapted to transmit and receive communications signals to and from the orbiting transceiver. A processor coupled to one of the first and second transceivers determines a first and second range measurement between each of the first and second transceivers and the orbiting transceiver, respectively, as well as corresponding first and second range rates representative of a time rate of change of the first and second range measurements. The processor then determines a circle of intersection representative of the set of possible positions for the orbiting transceiver based on the first range measurement and the first range rate wherein the circle of intersection includes a specific orientation in space, a specific radius and a center in a specific, threedimensional position relative to the known position of the first transceiver. An angular position of the orbiting transceiver is then determined along the circle of intersection based on the known position of the second transceiver and the second range measurement. Finally, the position of the orbiting transceiver is determined based on the circle of intersection and the angular position.
124 Citations
System and method for vital signs alerting privileged recipients  
Patent #
US 10,194,293 B2
Filed 01/06/2018

Current Assignee
Johnson Michael J., William J. Johnson

Sponsoring Entity
Johnson Michael J., William J. Johnson

Location blocking service from a web advertiser  
Patent #
US 10,217,137 B2
Filed 12/11/2017

Current Assignee
Google LLC

Sponsoring Entity
Google LLC

Location blocking service from a web advertiser  
Patent #
US 9,852,450 B2
Filed 10/10/2016

Current Assignee
Google LLC

Sponsoring Entity
Google LLC

System and method for vector processing on behalf of image aperture aim  
Patent #
US 9,584,993 B2
Filed 08/01/2016

Current Assignee
Billjco LLC

Sponsoring Entity
Billjco LLC

System and method for sound wave triggered content  
Patent #
US 10,111,034 B2
Filed 07/24/2016

Current Assignee
Billjco LLC

Sponsoring Entity
Billjco LLC

System and method for confirming data processing system target(s)  
Patent #
US 9,445,238 B2
Filed 11/13/2015

Current Assignee
Billjco LLC

Sponsoring Entity
Jason M. Johnson, William J. Johnson

System and method for permission to access mobile location information  
Patent #
US 9,648,454 B2
Filed 07/06/2015

Current Assignee
Google LLC

Sponsoring Entity
ATT Intellectual Property I LP

System and method for service access via hopped wireless mobile device(s)  
Patent #
US 9,456,303 B2
Filed 06/28/2015

Current Assignee
Billjco LLC

Sponsoring Entity
William J. Johnson

Global positioning system (GPS) and doppler augmentation (GDAUG) and space location inertial navigation geopositioning system (SPACELINGS)  
Patent #
US 9,766,339 B2
Filed 04/28/2015

Current Assignee
Raytheon Company

Sponsoring Entity
Raytheon Company

System and method for targeting data processing system(s) with data  
Patent #
US 9,204,275 B2
Filed 11/24/2014

Current Assignee
Billjco LLC

Sponsoring Entity
Jason M. Johnson, William J. Johnson

Location blocking service from a wireless service provider  
Patent #
US 9,763,091 B2
Filed 10/27/2014

Current Assignee
Google LLC

Sponsoring Entity
Google LLC

System and method for remote control of appliances utilizing mobile locationbased applications  
Patent #
US 9,537,997 B2
Filed 08/29/2014

Current Assignee
ATT Intellectual Property I LP

Sponsoring Entity
ATT Intellectual Property I LP

System and method for remote control of appliances utilizing mobile locationbased applications  
Patent #
US 9,584,647 B2
Filed 08/29/2014

Current Assignee
Google LLC

Sponsoring Entity
ATT Intellectual Property I LP

Surveying wireless device users by location  
Patent #
US 9,460,445 B2
Filed 08/08/2014

Current Assignee
ATT Intellectual Property I LP

Sponsoring Entity
ATT Intellectual Property I LP

Surveying wireless device users by location  
Patent #
US 9,501,780 B2
Filed 08/08/2014

Current Assignee
Google LLC

Sponsoring Entity
ATT Intellectual Property I LP

System and method for location based exchange vicinity interest specification  
Patent #
US 9,113,295 B2
Filed 06/17/2014

Current Assignee
Billjco LLC

Sponsoring Entity
William J. Johnson

System and method for targeting data processing system(s) with data  
Patent #
US 8,942,693 B2
Filed 03/18/2014

Current Assignee
Billjco LLC

Sponsoring Entity
William J. Johnson

Anonymous location service for wireless networks  
Patent #
US 9,571,958 B2
Filed 12/13/2013

Current Assignee
Google LLC

Sponsoring Entity
ATT Intellectual Property I LP

Location based exchange operating system  
Patent #
US 8,942,732 B2
Filed 11/22/2013

Current Assignee
William J. Johnson

Sponsoring Entity
William J. Johnson

System and method for location based exchanges of data facilitating distributed location applications  
Patent #
US 8,942,733 B2
Filed 11/22/2013

Current Assignee
Billjco LLC

Sponsoring Entity
William J. Johnson

Serverless synchronized processing across a plurality of interoperating data processing systems  
Patent #
US 9,055,406 B2
Filed 11/22/2013

Current Assignee
Billjco LLC

Sponsoring Entity
William J. Johnson

Location based exchange permissions  
Patent #
US 9,088,868 B2
Filed 11/22/2013

Current Assignee
Billjco LLC

Sponsoring Entity
William J. Johnson

System and method for application search results by locational conditions  
Patent #
US 9,088,869 B2
Filed 11/22/2013

Current Assignee
Billjco LLC

Sponsoring Entity
William J. Johnson

System and method for servicefree location based applications  
Patent #
US 9,100,792 B2
Filed 11/22/2013

Current Assignee
Billjco LLC

Sponsoring Entity
William J. Johnson

System and method for determining mobile users of interest  
Patent #
US 9,253,597 B2
Filed 11/22/2013

Current Assignee
Billjco LLC

Sponsoring Entity
William J. Johnson

System and method for location based exchanges of data facilitating distributed locational applications  
Patent #
US 9,392,408 B2
Filed 11/22/2013

Current Assignee
Billjco LLC

Sponsoring Entity
William J. Johnson

System and method for situational proximity observation alerting privileged recipients  
Patent #
US 9,894,489 B2
Filed 09/30/2013

Current Assignee
Johnson Michael J., William J. Johnson

Sponsoring Entity
Johnson Michael J., William J. Johnson

System and method for location based exchanges of data facilitating distributed locational applications  
Patent #
US 8,750,823 B2
Filed 09/23/2013

Current Assignee
Billjco LLC

Sponsoring Entity
William J. Johnson

System and method for location based exchanges of data facilitating distributed locational applications  
Patent #
US 8,761,804 B2
Filed 09/23/2013

Current Assignee
Billjco LLC

Sponsoring Entity
William J. Johnson

System and method for timely whereabouts determination by a mobile data processing system  
Patent #
US 8,886,226 B2
Filed 09/23/2013

Current Assignee
Billjco LLC

Sponsoring Entity
William J. Johnson

Providing promotion information to a device based on location  
Patent #
US 8,990,287 B2
Filed 09/23/2013

Current Assignee
Wayport Inc.

Sponsoring Entity
Wayport Inc.

Surveying wireless device users by location  
Patent #
US 8,805,414 B2
Filed 09/16/2013

Current Assignee
Google LLC

Sponsoring Entity
ATT Intellectual Property I LP

Methods and systems of aggregating information of geographic context regions of social networks based on geographical locations via a network  
Patent #
US 9,477,991 B2
Filed 08/27/2013

Current Assignee
Snap Trends Inc.

Sponsoring Entity
Snap Trends Inc.

System and method for location based exchange vicinity interest specification  
Patent #
US 8,718,598 B2
Filed 08/21/2013

Current Assignee
Billjco LLC

Sponsoring Entity
William J. Johnson

System and method for mobile device usability by locational conditions  
Patent #
US 8,897,741 B2
Filed 08/21/2013

Current Assignee
Billjco LLC

Sponsoring Entity
William J. Johnson

System and method for sudden proximal user interface  
Patent #
US 8,897,742 B2
Filed 08/21/2013

Current Assignee
Billjco LLC

Sponsoring Entity
William J. Johnson

System and method for presenting application data by data processing system(s) in a vicinity  
Patent #
US 8,923,806 B2
Filed 08/21/2013

Current Assignee
Billjco LLC

Sponsoring Entity
William J. Johnson

System and method for application context location based configuration suggestions  
Patent #
US 9,014,658 B2
Filed 08/21/2013

Current Assignee
Billjco LLC

Sponsoring Entity
William J. Johnson

System and method for location based inventory management  
Patent #
US 9,078,095 B2
Filed 08/21/2013

Current Assignee
Billjco LLC

Sponsoring Entity
William J. Johnson

Location blocking service from a wireless service provider  
Patent #
US 8,874,140 B2
Filed 08/12/2013

Current Assignee
Google LLC

Sponsoring Entity
ATT Intellectual Property I LP

Identity blocking service from a wireless service provider  
Patent #
US 8,718,618 B2
Filed 07/22/2013

Current Assignee
ATT Intellectual Property I LP

Sponsoring Entity
ATT Intellectual Property I LP

Identity blocking service from a wireless service provider  
Patent #
US 8,718,618 W0
Filed 07/22/2013

Current Assignee
ATT Intellectual Property I LP

Sponsoring Entity
ATT Intellectual Property I LP

Identity blocking service from a wireless service provider  
Patent #
US 8,755,777 B2
Filed 07/22/2013

Current Assignee
Google LLC

Sponsoring Entity
ATT Intellectual Property I LP

System and method for automated content distribution objects  
Patent #
US 8,887,177 B2
Filed 07/08/2013

Current Assignee
Billjco LLC

Sponsoring Entity
Jason M. Johnson, William J. Johnson

Third party location query for wireless networks  
Patent #
US 9,078,094 B2
Filed 07/08/2013

Current Assignee
Google LLC

Sponsoring Entity
ATT Intellectual Property I LP

Providing an advertisement based on a geographic location of a wireless access point  
Patent #
US 8,892,736 B2
Filed 05/30/2013

Current Assignee
Wayport Inc.

Sponsoring Entity
Wayport Inc.

Providing information to a computing device based on known location and user information  
Patent #
US 8,929,915 B2
Filed 03/06/2013

Current Assignee
Wayport Inc.

Sponsoring Entity
Wayport Inc.

Location blocking service from a web advertiser  
Patent #
US 9,466,076 B2
Filed 03/05/2013

Current Assignee
Google LLC

Sponsoring Entity
ATT Intellectual Property I LP

Anonymous location service for wireless networks  
Patent #
US 8,645,505 B2
Filed 02/27/2013

Current Assignee
Google LLC

Sponsoring Entity
ATT Intellectual Property I LP

Methods and systems of aggregating information of social networks based on geographical locations via a network  
Patent #
US 8,843,515 B2
Filed 02/05/2013

Current Assignee
Snap Trends Inc.

Sponsoring Entity
Snap Trends Inc.

Methods and systems of advertising based on aggregated information of social networks within geographical locations via a network  
Patent #
US 9,626,446 B2
Filed 11/14/2012

Current Assignee
Snap Trends Inc.

Sponsoring Entity
Snap Trends Inc.

Receiving location based advertisements on a wireless communication device  
Patent #
US 8,583,723 B2
Filed 10/08/2012

Current Assignee
Wayport Inc.

Sponsoring Entity
Wayport Inc.

System and method for remote control of appliances utilizing mobile locationbased applications  
Patent #
US 8,825,035 B2
Filed 07/23/2012

Current Assignee
Google LLC

Sponsoring Entity
ATT Intellectual Property I LP

Method and apparatus for geographicbased communications service  
Patent #
US 8,631,128 B2
Filed 07/17/2012

Current Assignee
Wayport Inc.

Sponsoring Entity
Wayport Inc.

Ground location inertial navigation geopositioning system (groundlings)  
Patent #
US 9,075,126 B2
Filed 06/28/2012

Current Assignee
Raytheon Company

Sponsoring Entity
Raytheon Company

System and method for using location information to execute an action  
Patent #
US 9,020,489 B2
Filed 06/25/2012

Current Assignee
Google LLC

Sponsoring Entity
ATT Intellectual Property I LP

Position and Elevation Acquisition for Orbit Determination  
Patent #
US 20130293413A1
Filed 05/03/2012

Current Assignee
Raytheon Company

Sponsoring Entity
Raytheon Company

Position and elevation acquisition for orbit determination  
Patent #
US 9,365,303 B2
Filed 05/03/2012

Current Assignee
Raytheon Company

Sponsoring Entity
Raytheon Company

Location blocking service from a wireless service provider  
Patent #
US 8,509,813 B2
Filed 05/02/2012

Current Assignee
ATT Intellectual Property I LP

Sponsoring Entity
ATT Intellectual Property I LP

Identity blocking service from a wireless service provider  
Patent #
US 8,494,501 B2
Filed 04/02/2012

Current Assignee
ATT Intellectual Property I LP

Sponsoring Entity
ATT Intellectual Property I LP

Method and apparatus for geographicbased communications service  
Patent #
US 8,606,851 B2
Filed 12/06/2011

Current Assignee
Wayport Inc.

Sponsoring Entity
Wayport Inc.

Third party location query for wireless networks  
Patent #
US 8,483,724 B2
Filed 11/04/2011

Current Assignee
Google LLC

Sponsoring Entity
ATT Intellectual Property I LP

Anonymous location service for wireless networks  
Patent #
US 8,402,117 B2
Filed 09/27/2011

Current Assignee
ATT Intellectual Property I LP

Sponsoring Entity
ATT Intellectual Property I LP

Surveying wireless device users by location  
Patent #
US 8,538,456 B2
Filed 08/01/2011

Current Assignee
ATT Intellectual Property I LP

Sponsoring Entity
ATT Intellectual Property I LP

Identity blocking service from a wireless service provider  
Patent #
US 8,165,568 B2
Filed 02/22/2011

Current Assignee
ATT Intellectual Property I LP

Sponsoring Entity
ATT Intellectual Property I LP

System and method for automated content presentation objects  
Patent #
US 8,566,839 B2
Filed 05/14/2010

Current Assignee
Billjco LLC

Sponsoring Entity
Jason M. Johnson, William J. Johnson

Location Blocking Service for Wireless Networks  
Patent #
US 20100151815A1
Filed 12/29/2009

Current Assignee
ATT Intellectual Property I LP

Sponsoring Entity
ATT Intellectual Property I LP

Location blocking service for wireless networks  
Patent #
US 7,873,369 B2
Filed 12/29/2009

Current Assignee
ATT Intellectual Property I LP

Sponsoring Entity
ATT Intellectual Property I LP

System and method for location based exchanges of data facilitating distributed location applications  
Patent #
US 8,634,796 B2
Filed 11/13/2009

Current Assignee
Billjco LLC

Sponsoring Entity
William J. Johnson

Third Party Location Query for Wireless Networks  
Patent #
US 20100035636A1
Filed 10/16/2009

Current Assignee
Google LLC

Sponsoring Entity
Google LLC

Third party location query for wireless networks  
Patent #
US 7,844,284 B2
Filed 10/16/2009

Current Assignee
Google LLC

Sponsoring Entity
ATT Intellectual Property I LP

System and method for location based exchanges of data facilitating distributed locational applications  
Patent #
US 8,639,267 B2
Filed 10/03/2008

Current Assignee
Billjco LLC

Sponsoring Entity
William J. Johnson

SYSTEM AND METHOD FOR USING LOCATION INFORMATION TO EXECUTE AN ACTION  
Patent #
US 20080299957A1
Filed 08/08/2008

Current Assignee
Google LLC

Sponsoring Entity
Google LLC

System and method for using location information to execute an action  
Patent #
US 8,639,235 B2
Filed 08/08/2008

Current Assignee
Google LLC

Sponsoring Entity
ATT Intellectual Property I LP

System and method for using location information to execute an action  
Patent #
US 8,260,239 B2
Filed 06/03/2008

Current Assignee
ATT Intellectual Property I LP

Sponsoring Entity
ATT Intellectual Property I LP

System and method for location based exchanges of data facilitating distributed locational applications  
Patent #
US 8,600,341 B2
Filed 03/14/2008

Current Assignee
Billjco LLC

Sponsoring Entity
William J. Johnson

Locationbased security rules  
Patent #
US 8,644,506 B2
Filed 12/18/2007

Current Assignee
Google LLC

Sponsoring Entity
ATT Intellectual Property I LP

METHOD AND APPARATUS FOR GEOGRAPHICBASED COMMUNICATIONS SERVICE  
Patent #
US 20080049696A1
Filed 10/29/2007

Current Assignee
Wayport Inc.

Sponsoring Entity
Wayport Inc.

Method and apparatus for geographicbased communications service  
Patent #
US 8,199,733 B2
Filed 10/29/2007

Current Assignee
Wayport Inc.

Sponsoring Entity
Wayport Inc.

Method and apparatus for geographicbased communications service  
Patent #
US 8,095,647 B2
Filed 10/25/2007

Current Assignee
Wayport Inc.

Sponsoring Entity
Wayport Inc.

Method and apparatus for geographicbased communications service  
Patent #
US 8,250,204 B2
Filed 10/25/2007

Current Assignee
Wayport Inc.

Sponsoring Entity
Wayport Inc.

Method and apparatus for geographicbased communications service  
Patent #
US 8,509,246 B2
Filed 10/25/2007

Current Assignee
ATT Intellectual Property I LP

Sponsoring Entity
ATT Intellectual Property I LP

Surveying wireless device users by location  
Patent #
US 20070142060A1
Filed 12/12/2006

Current Assignee
ATT Intellectual Property I LP

Sponsoring Entity
ATT Intellectual Property I LP

Surveying wireless device users by location  
Patent #
US 8,010,126 B2
Filed 12/12/2006

Current Assignee
ATT Intellectual Property I LP

Sponsoring Entity
ATT Intellectual Property I LP

Third party location query for wireless networks  
Patent #
US 20070105565A1
Filed 10/30/2006

Current Assignee
Google LLC

Sponsoring Entity
Google LLC

Third party location query for wireless networks  
Patent #
US 7,636,575 B2
Filed 10/30/2006

Current Assignee
Google LLC

Sponsoring Entity
ATT Intellectual Property I LP

System and method for using location information to execute an action  
Patent #
US 20070042789A1
Filed 09/28/2006

Current Assignee
ATT Intellectual Property I LP

Sponsoring Entity
ATT Intellectual Property I LP

System and method for using location information to execute an action  
Patent #
US 7,941,130 B2
Filed 09/28/2006

Current Assignee
ATT Intellectual Property I LP

Sponsoring Entity
ATT Intellectual Property I LP

System and method for using location information to execute an action  
Patent #
US 20070010260A1
Filed 09/15/2006

Current Assignee
ATT Intellectual Property I LP

Sponsoring Entity
ATT Intellectual Property I LP

System and method for using location information to execute an action  
Patent #
US 7,412,234 B2
Filed 09/15/2006

Current Assignee
ATT Intellectual Property I LP

Sponsoring Entity
ATT Intellectual Property I LP

Method and System for Locating Objects  
Patent #
US 20090153305A1
Filed 07/07/2006

Current Assignee
Ambrosetti Antonio

Sponsoring Entity
Ambrosetti Antonio

Location blocking service from a wireless service provider  
Patent #
US 20060189327A1
Filed 04/19/2006

Current Assignee
ATT Intellectual Property I LP

Sponsoring Entity
ATT Intellectual Property I LP

Location blocking service from a wireless service provider  
Patent #
US 7,509,133 B2
Filed 04/19/2006

Current Assignee
ATT Intellectual Property I LP

Sponsoring Entity
ATT Intellectual Property I LP

Distributed network communication system to provide wireless access to a computing device at a reduced rate  
Patent #
US 8,588,130 B2
Filed 04/12/2006

Current Assignee
Wayport Inc.

Sponsoring Entity
Wayport Inc.

Dynamically modifying the display of a computing device to provide advertisements  
Patent #
US 7,840,689 B2
Filed 04/11/2006

Current Assignee
Wayport Inc.

Sponsoring Entity
Wayport Inc.

Providing advertisements to a computing device based on a predetermined criterion of a wireless access point  
Patent #
US 8,478,887 B2
Filed 03/27/2006

Current Assignee
Wayport Inc.

Sponsoring Entity
Wayport Inc.

Anonymous location service for wireless networks  
Patent #
US 20060195570A1
Filed 03/15/2006

Current Assignee
ATT Intellectual Property I LP

Sponsoring Entity
ATT Intellectual Property I LP

Anonymous location service for wireless networks  
Patent #
US 7,418,503 B2
Filed 03/15/2006

Current Assignee
ATT Intellectual Property I LP

Sponsoring Entity
ATT Intellectual Property I LP

System and method for using location information to execute an action  
Patent #
US 7,245,925 B2
Filed 12/15/2005

Current Assignee
Google LLC

Sponsoring Entity
ATT Intellectual Property I LP

System and method for using location information to execute an action  
Patent #
US 20060089134A1
Filed 12/09/2005

Current Assignee
ATT Intellectual Property I LP

Sponsoring Entity
ATT Intellectual Property I LP

System and method for using location information to execute an action  
Patent #
US 20060099966A1
Filed 12/09/2005

Current Assignee
ATT Intellectual Property I LP

Sponsoring Entity
ATT Intellectual Property I LP

System and method for using location information to execute an action  
Patent #
US 7,383,052 B2
Filed 12/09/2005

Current Assignee
ATT Intellectual Property I LP

Sponsoring Entity
ATT Intellectual Property I LP

System and method for using location information to execute an action  
Patent #
US 7,593,712 B2
Filed 12/09/2005

Current Assignee
ATT Intellectual Property I LP

Sponsoring Entity
ATT Intellectual Property I LP

Identity blocking service from a wireless service provider  
Patent #
US 20060030335A1
Filed 10/17/2005

Current Assignee
ATT Intellectual Property I LP

Sponsoring Entity
ATT Intellectual Property I LP

Identity blocking service from a wireless service provider  
Patent #
US 7,603,110 B2
Filed 10/17/2005

Current Assignee
ATT Intellectual Property I LP

Sponsoring Entity
ATT Intellectual Property I LP

Locationbased security rules  
Patent #
US 20050272445A1
Filed 07/21/2005

Current Assignee
Google LLC

Sponsoring Entity
Google LLC

Locationbased security rules  
Patent #
US 7,428,411 B2
Filed 07/21/2005

Current Assignee
Google LLC

Sponsoring Entity
ATT Intellectual Property I LP

Providing information to a computing device based on known location and user information  
Patent #
US 8,417,763 B2
Filed 02/11/2004

Current Assignee
Wayport Inc.

Sponsoring Entity
Wayport Inc.

Location blocking service for wireless networks  
Patent #
US 7,664,509 B2
Filed 11/12/2003

Current Assignee
Google LLC

Sponsoring Entity
ATT Intellectual Property I LP

Ranging system for determining ranging information of a spacecraft  
Patent #
US 20040140930A1
Filed 09/29/2003

Current Assignee
SES ASTRA S.A.

Sponsoring Entity
SES ASTRA S.A.

Ranging system for determining ranging information of a spacecraft  
Patent #
US 7,512,505 B2
Filed 09/29/2003

Current Assignee
SES ASTRA S.A.

Sponsoring Entity
SES ASTRA S.A.

Stratospheric platform based surface vehicle tracking and mobile data network  
Patent #
US 20020107034A1
Filed 12/14/2001

Current Assignee
Hughes Electronics Corporation

Sponsoring Entity
Hughes Electronics Corporation

Stratospheric platform based surface vehicle tracking and mobile data network  
Patent #
US 6,941,107 B2
Filed 12/14/2001

Current Assignee
Hughes Electronics Corporation

Sponsoring Entity
Hughes Electronics Corporation

Method and system for determining a position of an object using twoway ranging in a polystatic satellite configuration  
Patent #
US 6,563,457 B2
Filed 11/14/2001

Current Assignee
Hughes Electronics Corporation

Sponsoring Entity
Hughes Electronics Corporation

Position location of multiple transponding platforms and users using twoway ranging as a calibration reference for GPS  
Patent #
US 6,785,553 B2
Filed 04/24/2001

Current Assignee
Hughes Electronics Corporation

Sponsoring Entity
Hughes Electronics Corporation

System and method for using location information to execute an action  
Patent #
US 7,116,977 B1
Filed 12/19/2000

Current Assignee
Google LLC

Sponsoring Entity
Bellsouth Intellectual Property Corporation

Communication system for providing information on position of communication party  
Patent #
US 7,065,348 B1
Filed 09/28/2000

Current Assignee
Toshiba Corporation

Sponsoring Entity
Toshiba Corporation

Concurrent communications between a user terminal and multiple stratospheric transponder platforms  
Patent #
US 6,941,138 B1
Filed 09/05/2000

Current Assignee
Hughes Electronics Corporation

Sponsoring Entity
Hughes Electronics Corporation

Synchronization method and apparatus for multiplatform communication system  
Patent #
US 6,757,546 B1
Filed 08/21/2000

Current Assignee
Hughes Electronics Corporation

Sponsoring Entity
Hughes Electronics Corporation

Method and system for delivery of a calling party's location  
Patent #
US 7,796,998 B1
Filed 08/01/2000

Current Assignee
ATT Intellectual Property I LP

Sponsoring Entity
ATT Intellectual Property I LP

User positioning technique for multiplatform communication system  
Patent #
US 6,920,309 B1
Filed 06/06/2000

Current Assignee
Hughes Electronics Corporation

Sponsoring Entity
Hughes Electronics Corporation

Resource allocation method for multiplatform communication system  
Patent #
US 7,215,954 B1
Filed 06/06/2000

Current Assignee
Hughes Electronics Corporation

Sponsoring Entity
Hughes Electronics Corporation

Multinode wireless communication system with multiple transponding platforms  
Patent #
US 7,089,000 B1
Filed 05/30/2000

Current Assignee
Hughes Electronics Corporation

Sponsoring Entity
Hughes Electronics Corporation

Multinode pointtopoint satellite communication system employing multiple geo satellites  
Patent #
US 6,990,314 B1
Filed 05/22/2000

Current Assignee
Hughes Electronics Corporation

Sponsoring Entity
Hughes Electronics Corporation

No References
22 Claims
 1. A method for determining a position of an orbiting transceiver in a communications network including at least a first and second transceiver at a first and second known location, respectively, on Earth, the first and second transceivers adapted to transmit and receive communications signals to and from the orbiting transceiver, the method comprising:
 determining a first and second range measurement between each of the first and second transceivers and the orbiting transceiver, respectively;
determining a first and second range rate corresponding to a time rate of change of the first and second range measurements, respectively;
determining a circle of intersection representative of the set of possible positions for the orbiting transceiver based on the first range measurement and the first range rate, the circle of intersection having a specific orientation in space, a specific radius and a center in a specific, three dimensional position relative to the known position of the first transceiver;
determining an angular position of the orbiting transceiver along the circle of intersection based on the known position of the second transceiver and the second range measurement; and
determining the position of the orbiting transceiver based on the circle of intersection and the angular position.  View Dependent Claims (2, 3, 4, 5, 6, 7, 8, 9, 10, 11)
 determining a first and second range measurement between each of the first and second transceivers and the orbiting transceiver, respectively;
 12. A system for determining a position of an orbiting transceiver in a communications network, the system comprising:
 at least a first and second transceiver at a first and second known location, respectively, on Earth, the first and second transceivers adapted to transmit and receive communications signals to and from the orbiting transceiver; and
a processor coupled to one of the first and second transceivers operative to determine a first and second range measurement between each of the first and second transceivers and the orbiting transceiver, respectively, determine a first and second range rate corresponding to a time rate of change of the first and second range measurements, respectively, determine a circle of intersection representative of the set of possible positions for the orbiting transceiver based on the first range measurement and the first range rate, wherein the circle of intersection includes a specific orientation in space, a specific radius, and a center in a specific, three dimensional position relative to the known position of the first transceiver, determine an angular position of the orbiting transceiver along the circle of intersection based on the known position of the second transceiver and the second range measurement, and determine the position of the orbiting transceiver based on the circle of intersection and the angular position.  View Dependent Claims (13, 14, 15, 16, 17, 18, 19, 20, 21, 22)
 at least a first and second transceiver at a first and second known location, respectively, on Earth, the first and second transceivers adapted to transmit and receive communications signals to and from the orbiting transceiver; and
1 Specification
This application is related to U.S. patent application Ser. No. 08/803,935, filed Feb. 21, 1997, entitled "Method And System For Determining A Position Of A Transceiver Unit Utilizing TwoWay Ranging in a Polystatic Satellite Configuration".
TECHNICAL FIELDThis invention relates to methods and systems for determining a position of a communication satellite utilizing twoway ranging through multiple transceivers.
BACKGROUND ARTCurrent Automatic Dependent Surveillance (ADS) technology, such as Global Positioning System (GPS), Wide Area Augmentation System (WAAS) or GLONASS, provides positioning information utilizing satellite transmissions. For example, the GPS, developed and deployed by the U.S. Department of Defense, consists of 24 satellites orbiting the earth twice a day at an altitude of 12,000 miles, as well as five ground stations to monitor and manage the satellite constellation. Using atomic clocks and location data, GPS satellites transmit continuous time and position information 24 hours a day to a GPS receiver, which listens to four or more satellites at once to determine a user's position. By measuring the time interval between the transmission and the reception of a satellite signal, the GPS receiver calculates the distance between the user and each satellite, and then uses the distance measurements of at least four satellites to arrive at a position.
Such systems, however, utilize oneway ranging in which an accurate, synchronized clock is required at each station. Any synchronization error or error regarding the location of one of the satellites results in an error in the determined position of the target vehicle.
Any type of satellite based navigation system must determine the position and velocity of its ranging satellites to extremely high accuracy. Position and velocity determination are vital parts of the stationkeeping function for any satellite, even one which is not part of a navigation system.
One known system requires the satellite to have simultaneous contact with all of its ground stations in order to obtain position information. This, however, may not always be possible for a satellite in a nongeostationary orbit which moves relative to the Earth and, therefore, periodically loses access to any given point on the Earth's surface.
DISCLOSURE OF THE INVENTIONIt is thus a general object of the present invention to provide a method and system for determining the position of a communication satellite utilizing twoway ranging.
In carrying out the above object and other objects, features, and advantages of the present invention, a method is provided for determining a position of an orbiting transceiver in a communications network including at least a first and second transceiver at a first and second known location, respectively, on Earth. The first and second transceivers are adapted to transmit and receive communications signals to and from the orbiting transceiver. The method includes determining a first and second range measurement between each of the first and second transceivers and the orbiting transceiver, respectively. The method further includes determining a first and second range rate corresponding to a time rate of change of the first and second range measurements, respectively. The method also includes determining a circle of intersection representative of the set of possible positions for the orbiting transceiver based on the first range and the first range rate wherein the circle of intersection includes a specific orientation in space, a specific radius and a center in a specific, three dimensional position relative to the known position of the first transceiver. Still further, the method includes determining an angular position of the orbiting transceiver along the circle of intersection based on the known position of the second transceiver and the second range measurement. Finally, the method includes determining the position of the orbiting transceiver based on the circle of intersection and the angular position.
In further carrying out the above object and other objects, features, and advantages of the present invention, a system is also provided for carrying out the steps of the above described method. The system includes at least a first and second transceiver at a first and second known location, respectively, on Earth. The first and second transceivers are adapted to transmit and receive communications signals to and from the orbiting transceiver. A processor coupled to one of the first and second transceivers is operative to determine a first and second range measurement between each of the first and second transceivers and the orbiting transceiver, respectively. The processor is further operative to determine a first and second range rate corresponding to a time rate of change of the first and second range measurements, respectively. Still further, the processor is operative to determine a circle of intersection representative of the set of possible positions for the orbiting transceiver based on the first range measurement and the first range rate wherein the circle of intersection includes a specific orientation in space, a specific radius, and a center in a specific, threedimensional position relative to the known position of the first transceiver. The processor further determines an angular position of the orbiting transceiver along the circle of intersection based on the known position of the second transceiver and the second range measurement. Finally, the processor determines the position of the orbiting transceiver based on the circle of intersection and the angular position.
The above object and other objects, features and advantages of the present invention are readily apparent from the following detailed description of the best mode for carrying out the invention when taken in connection with the accompanying drawings.
BRIEF DESCRIPTION OF THE DRAWINGSFIG. 1 is a diagram illustrating a communication system employing the method and system of the present invention;
FIG. 2 is a diagram illustrating the geometrical plane containing two ranging nodes of the present invention;
FIG. 3 is a diagram illustrating an enlarged view of the plane of FIG. 2;
FIG. 4 is a diagram illustrating a circle of intersection extending above and below the XY plane of FIG. 2; and
FIG. 5 is a diagram similar to FIG. 3, but adding x and y unit vectors, and a point denoted R.sub.q associated with a later step in the position determination.
BEST MODES FOR CARRYING OUT THE INVENTIONFIG. 1 diagrammatically illustrates a communication system with a typical geometry for practicing the present invention. This system, denoted generally by reference numeral 10, includes a target communication satellite 12 whose position (denoted R.sub.0) is to be determined. The system 10 in FIG. 1 also includes a primary ranging node (RN) 14such as a satellite ground stationat the position denoted R.sub.1, as well as secondary ranging nodes 11 at positions R.sub.2 and R.sub.3. There may be many other secondary RNs at positions R.sub.4, R.sub.5, etc. These additional RNs are not shown in FIG. 1; however, their modes of operation are identical to those described for the RNs at R.sub.2 and R.sub.3. A secondary RN may be any device that includes a transceiver such as, for example, automobiles, mobile telephones, aircraft, or the like, as long as they are stationary during operation of the present invention. As shown in FIG. 1, the distances between the satellite 12 at R.sub.0 and all the RNs at R.sub.1, R.sub.2, R.sub.3, etc. are denoted a.sub.i =.vertline.R.sub.i R.sub.0.vertline., where i>0.
For example, a parked commercial airliner can act as a secondary RN. When a commercial airliner is parked at an airport gate, its pilot enters the precise position of the aircraft into the onboard navigation system. This information allows the inertial navigation system to correct for drift and accumulated position error. As long as the aircraft is parked at the gate, its position is known to a high degree of accuracy. If it is equipped with an appropriate transceiver, then it can act as a secondary RN at position R.sub.2 or R.sub.3. If a large number of aircraft are thus equipped, and a significant percentage of them are parked at a gate at any given time, then the satellite 12 at R.sub.0 can use these parked aircraft at R.sub.2 and R.sub.3 as RNs.
The RNs at R.sub.1, R.sub.2, R.sub.3, etc. communicate with each other by relaying messages through the satellite 12 at position R.sub.0. In the course of operation of the present invention, the primary RN 14 receives two distinct kinds of signals: ranging codes broadcast by the satellite 12, and reply codes which also are relayed through the satellite 12, but originate at secondary RNs 11. A processor 16 located at the primary RN 14 at R.sub.1 distinguishes the two types of signals, and makes the computation required for the determination of the position of the satellite 12.
The position of the satellite 12 is determined according to the following sequence of steps. First, at a precisely measured time, and on a precisely measured carrier frequency, the primary RN 14 at R.sub.1 transmits a ranging signal to the satellite 12 at R.sub.0, which broadcasts it to all of the RNs, including the primary RN 14 at R.sub.1, and secondary RNs 11 at R.sub.2, R.sub.3, etc. The secondary RNs 11 then transmit their distinctive reply code signals back to the satellite 12 for transmission back to primary RN 14. Immediately upon receiving the ranging and reply code signals, the processor 16 at the primary ranging node 14 (R.sub.1) measures the carrier frequency and arrival time of the ranging signal. It then computes the difference between the ranging code's arrival time and the time of its original transmission from the primary ranging node 14. This time difference is a propagation time, denoted t.sub.p1, and is associated with the distance of the transmission of ranging signals from the primary RN 14 at R.sub.1, to the satellite 12 at R.sub.0, and back to R.sub.1. This distance, denoted D.sub.1, is equal to 2a.sub.1. The processor 16 also measures the difference in carrier frequency between the received reply code signal and that of the originally transmitted ranging code. From that frequency difference the processor computes the Doppler frequency shift, as described below.
Before continuing with the description of the processor's computations, we describe the measurements made on signals passing through the secondary RNs 11. As mentioned above, immediately upon receiving a ranging code, the secondary RNs 11 at positions R.sub.2, R.sub.3, R.sub.4. etc. each transmit a reply ranging signal to the primary RN 14 at R.sub.1 via the satellite 12 at R.sub.0. These reply signals are designed such that they uniquely identify the particular ranging code to which they are replying, and the secondary RN 11 from which they are transmitted. The processor 16 at R.sub.1 can then determine the position R.sub.i of the secondary RN 11 which initiated the particular reply code, where R.sub.i is an element of the set of all known secondary RN positions {R.sub.2, R.sub.3, . . . }. The processor 16 also measures the arrival time of each reply signal, and computes the difference between that time and the time of original transmission from the primary ranging node 14 of the ranging code which elicited the reply signal. That difference is denoted as a propagation time, t.sub.pi. Thus, the measured times for round trip communication with secondary RNs at R.sub.2, R.sub.3, etc. are denoted t.sub.p2, t.sub.p3, etc. respectively.
Referring again to FIG. 1, we see that the paths associated with propagation times t.sub.pi for i>1 are of a different shape than the path associated with t.sub.p1 and D.sub.1. For i>1, D.sub.i is the path length from the primary ranging node 14 at R.sub.1, to the satellite 12 at R.sub.0, to the secondary ranging node 11 at R.sub.i, then back to satellite 12 at R.sub.0, and finally, back to primary ranging node 14 at R.sub.1. To a first approximation, D.sub.i =2(a.sub.i +a.sub.1) for i>1. This equality is only approximate because it neglects the fact that the satellite 12 moves during the time of the signal's trip from R.sub.0 to R.sub.i and back to R.sub.0. In reality, the satellite 12 will be in a different place when relaying the initial ranging signal from R.sub.1 to R.sub.i than when relaying the reply signal from R.sub.i back to R.sub.1. We can correct for this effect, but in the interest of clarity, we neglect it in the present description of the system, but not in the final implementation.
We now outline a computational procedure which transforms a set of frequency shifts and propagation time measurements first into a set of communication path lengths (D.sub.1, D.sub.2, D.sub.3, etc.), then into a set of ranges (a.sub.1, a.sub.2, a.sub.3, etc.), and finally, into a three dimensional, vector position of the satellite 12. The measurements of propagation time include time for an electromagnetic wave to travel the prescribed path, and also delay times associated with the electronic detection and retransmission of messages. These intervals are all measured with a single clock 18 at R.sub.1 ; therefore, secondary RNs at R.sub.2, R.sub.3, etc. do not need clocks, and the system depicted in FIG. 1 does not have any clock bias errors such as those in The Global Positioning System (GPS).
Signal Delay Times and Distances
For the path corresponding to each measured propagation time t.sub.p1, t.sub.p2, t.sub.p3, etc., the processor 16 at R.sub.1 computes an estimate of the total time delay, t.sub.di, for all transceivers as well as tropospheric and ionospheric propagation along that path. The corresponding delay times t.sub.d1, t.sub.d2, t.sub.d3, etc. are then used in the computation of the corresponding distances between the transceivers.
The ionosphere and atmosphere are two elements in the signal path that slow the signal to less than the vacuum speed of light, which is defined as c=299,792,458 meters/second. The delay time, t.sub.di, is defined as the propagation time of a signal over the i.sup.th path minus the time required for an electromagnetic wave to traverse a path of the same length in a perfect vacuum. Given that D.sub.i is the length of the i.sup.th communication path, D.sub.i =c(t.sub.pi t.sub.di). Delay time is the sum of many components, due to the ionosphere, atmosphere, transmitter, receiver, etc., and is never negative.
The atmospheric and ionospheric signal delays have been thoroughly studied, and mathematical models have been constructed which can predict these components of delay time (for example: Black, H. D., "An Easily Implemented Algorithm for the Tropospheric Range Correction," JOURNAL OF GEOPHYSICAL RESEARCH, vol. 83, no. B4, p. 1825, (1978); and for the ionospheric effect: Klobuchar, J. A., "Design and Characteristics of the GPS Ionospheric Time Delay Algorithm for Single Frequency Users," IEEE PLANS '186: Position Location and Navigation Symposium, Las Vegas, Nev., Nov. 4, 1986, p. 280).
The ionospheric delay is harder to predict with high accuracy than that of the troposphere because the ionosphere undergoes more radical variations over time. If a particular application demands the highest possible accuracy, then the system 10 will rely upon additional measurements rather than a fixed mathematical model in order to determine the ionospheric propagation delay. Because that delay is inversely proportional to the square of the radio carrier frequency, and because it varies with location and time, additional measurements of roundtrip signal propagation time through different ground reference points and/or at different carrier frequencies will provide the continuous measurements of the ionosphere which are needed for the highest accuracy range measurements. The availability of a large number of secondary RNs provides extra information which can be used to determine the characteristics of the ionosphere, the troposphere, and other sources of signal delay in the communication path. The present invention can compensate for the time varying tropospheric and ionospheric signal propagation delays more accurately than a system which contains a very small number of RNs or turnaround ranging stations.
Other sources of delay involving the system's antennas, transmitters, and receivers (including multipath, etc.) are similarly estimated. All of the estimated time delays associated with the i.sup.th communication path are estimated, and added together to produce the total estimated time delay, t.sub.di. All such delay times t.sub.d1, t.sub.d2, t.sub.d3, etc. are then used in the computation of the corresponding communication path lengths D.sub.1, D.sub.2, D.sub.3, etc., via the formula
D.sub.i =c(t.sub.pi t.sub.di).
Based on D.sub.1, D.sub.2, D.sub.3, etc., the system 10 computes the ranges a.sub.1, a.sub.2, a.sub.3, etc., via formulas such as: a.sub.1 =D.sub.1 /2, and a.sub.i =(D.sub.i D.sub.1)/2. As mentioned previously, the final implementation involves more complicated equations than these, due to the motion of the satellite 12 during the time interval between its broadcast of the ranging signal and its relay of the reply signal. The correction is omitted here to simplify the presentation. The processor 16 then uses those ranges and the known positions of the RNs 14, 11, (a.sub.i and R.sub.i for i>0) to compute the satellite's position R.sub.0 by means of one of the two procedures described below.
We derive two algorithms for determination of the (three dimensional) vector position of the satellite 12; a "Three Sphere Algorithm" which uses range measurements involving at least three ground transceivers, and a "Range and Range Rate Algorithm" which uses range measurements from only two ground transceivers. "Range Rate" refers to the time rate of change of range (da.sub.i /dt), and is derived from a combination of the time sequence of range measurements and measurements of the Doppler frequency shift of the radio signal, which is described briefly below. The "Range and Range Rate Algorithm" applies only to satellites which are not in geostationary orbit; whereas, the "Three Range Algorithm" applies to geostationary or nongeostationary satellites.
Doppler Frequency Shift
If a transmitter at position R.sub.1 sends a radio signal at carrier frequency f.sub.1 to a receiver at position R.sub.0 which is approaching or receding, then the receiver will detect the signal at a different frequency f.sub.0 due to the Doppler effect. The size of the Doppler frequency shift depends on the radial velocity, which we denote as v.sub.r, and which equals the time rate of change of the distance between the transmitter and receiver. If the velocity of the receiver relative to the transmitter is much smaller in magnitude than the speed of light, then (approximately),
f.sub.0 /f.sub.1 =1v.sub.r /c.
If another signal is transmitted at frequency f.sub.0 from R.sub.0 back to R.sub.1, then the Doppler shift is compounded, so that the frequency detected at R.sub.1 is
f.sub.11 f.sub.0 (1v.sub.r /c)=f.sub.1 (1v.sub.r /c).sup.2,
which is approximately
f.sub.11 =f.sub.1 (12v.sub.r /c.
The Doppler frequency shift (f.sub.d) is the difference between the frequency f.sub.11 received at point 1 after the roundtrip, and that of the signal initially transmitted from point 1, which is f.sub.1 ; therefore,
f.sub.d =f.sub.11 f.sub.1 =2f.sub.1 v.sub.r /c.
The Doppler shift measurement therefore yields a measurement of the radial velocity,
v.sub.r =cf.sub.d /(2f.sub.1),
and this information contributes to the orbit determination, as described below.
Geometry
In this invention, we present two algorithms which may be used to convert a set of scalar measurements into the vector position of the satellite 12. The "Three Sphere" positioning algorithm employs round trip signal propagation between a transceiver at an unknown location (which we denote R.sub.0) and a set of three or more RNs. The "Range and Range Rate" positioning algorithm requires round trip communication between the transceiver at R.sub.0 and only two RNs, but also requires that R.sub.0 be moving relative to the RNs. In both of these approaches, the round trip pattern of communication implies that all round trip propagation time intervals can be measured from a single clock located at the primary RN. There is no need for highly accurate or synchronized clocks at any other location in the system.
In contrast, a typical GPS receiver requires measurements of signals from four different transmitters in order to solve for the receiver's position. Due to the one way nature of the communication, the GPS transmitters must have highly accurate and precisely synchronized clocks, and each GPS receiver must solve for the precise local time while determining its three dimensional position.
An additional advantage of the present invention is that these two algorithms are direct, closed form solutions to the nonlinear system of equations which define the satellite's position; whereas, GPS uses an iterative solution to a linearized approximation of an intrinsically nonlinear problem. The direct solution is more computationally efficient, faster, and more accurate.
"Three Sphere" Algorithm
The three sphere algorithm transforms the ranges (a.sub.1, a.sub.2, a.sub.3, etc.) computed above and the known positions of three or more RNs (R.sub.1, R.sub.2, R.sub.3, etc.) into the three dimensional, vector position of the satellite 12 relative to the Earth. We use the following notation to aid in the concise expression of the algorithm.
Vector Notation
In all of the following equations, upper case letters denote three component vectors, and lower case letters denote scalar (nonvector) quantities.
A.times.B=cross product (vector product) of vectors A and B, such that if C=A.times.B then the three components of C are:
c.sub.x =a.sub.y b.sub.z a.sub.z b.sub.y
c.sub.y =a.sub.z b.sub.x a.sub.x b.sub.z
c.sub.z =a.sub.x b.sub.y a.sub.y b.sub.x
Note that the resulting C is perpendicular to A and to B.
A.B=a.sub.x b.sub.x +a.sub.y b.sub.y +a.sub.z b.sub.z =dot product (also known as the scalar product or inner product) of vectors A and B.
A.sup.2 =A.A=square of the length of vector A, which equals the sum of the squares of the three components (x, y, z) of the vector (as per the Pythagorean theorem in three dimensions).
.vertline.A.vertline.=(A.A).sup.1/2 =length of vector A.
R.sub.j =vector which points from the center of the Earth to the position of the j.sup.th RN, where j=[1, 2, 3, . . . ].
R.sub.0 =vector which points from the center of the Earth to the position of the satellite to be located.
R.sub.jk =R.sub.j R.sub.k =vector from point k to point j.
a.sub.j =.vertline.R.sub.0j.vertline.=distance measurement from the j.sup.th RN to the satellite.
Intersection of Two Spheres
The satellite's position, R.sub.0, will be determined relative to the RN positions (R.sub.1, R.sub.2, . . . ). The following geometrical and mathematical derivation of R.sub.0 is performed not directly in Earth centered coordinates, but rather, in a coordinate system defined by the positions of the RNs. This choice of coordinate system separates the vector R.sub.01 (which points from the known R.sub.1 to the unknown R.sub.0) into three mutually perpendicular components (denoted R.sub.p1, R.sub.qp, R.sub.0q). In this carefully selected coordinate system, each of these vector components can be computed directly (noniteratively) from the measurements (a.sub.1, a.sub.2, a.sub.3). The sum of these three vectors and the position vector of the primary RN (R.sub.1) is the desired location of the satellite (R.sub.0) Symbolically:
R.sub.0 =R.sub.0q +R.sub.qp +R.sub.p1 +R.sub.1.
The following discussion describes the geometrical process of locating the intermediate points R.sub.p and R.sub.q in space, and then deriving algebraic formulas for those positions.
FIG. 2 shows a plane containing the primary RN 14 at R.sub.1 and the secondary RN at R.sub.2. It also shows the vector R.sub.21 =R.sub.2 R.sub.1, which points from R.sub.1 to R.sub.2. The measurement al limits the possible positions of R.sub.0 to the points on a sphere of radius a.sub.1 which is centered at R.sub.1. The intersection of that sphere with the plane of FIG. 2 is a circle of radius a.sub.1. That circle is labeled .vertline.R.sub.01.vertline.=a.sub.1. Although R.sub.0 might lie on this circle, it is far more likely to lie above or below the plane depicted in FIG. 2.
Similarly, the measurement a.sub.2 limits the possible positions of R.sub.0 to the points on a sphere of radius a.sub.2 which is centered at R.sub.2. The intersection of that sphere with the plane of FIG. 2 is the circle labeled .vertline.R.sub.02.vertline.=a.sub.2. Since R.sub.0 must lie on both of these spheres, the only possible solution points for R.sub.0 in the plane of FIG. 2 are the points labeled R.sub.b and R.sub.c, where the two circles intersect; however, there are other possible solution points above and below this plane.
In order to visualize the location of all possible points which lie on both of the aforementioned spheres, one may slowly flip over the picture in FIG. 2, doing so in a way that does not move vector R.sub.21. During this slow rotation of the picture, one sees that the circles in the plane trace out spheres in three dimensions, and that R.sub.b and R.sub.c trace out a circle of possible positions for R.sub.0. That circle lies in another plane which is perpendicular to the vector R.sub.21. In other words, the circle defined by the intersection of the two spheres sticks out at a right angle from the plane depicted in FIG. 2.
Having considered the geometrical construction of the circle of intersection, we now develop equations which quantify the position of the structures discussed above. FIG. 3 is an enlargement of the central region of FIG. 2, with additional details added. The vector from R.sub.1 to R.sub.b is denoted .vertline.R.sub.b1.vertline.=a.sub.1, where a.sub.1 is its length. Similarly, the vector from R.sub.2 to R.sub.b is denoted .vertline.R.sub.b2.vertline.=a.sub.2, where a.sub.2 is its length.
The vector R.sub.bc which points from R.sub.c to R.sub.b is perpendicular to the vector R.sub.21. The point of intersection of R.sub.bc with R.sub.21 is denoted R.sub.p. This definition of R.sub.p allows us to break the triangle with vertices at R.sub.1, R.sub.2, R.sub.b into one right triangle with vertices at R.sub.1, R.sub.p, R.sub.b and another right triangle with vertices at R.sub.2, R.sub.p, R.sub.b. The lengths of the sides of right triangles are related by the Pythagorean theorem: x.sup.2 +y.sup.2 =r.sup.2, where x and y represent the lengths of the two perpendicular sides, and r represents the length of the third side. We now use this theorem to determine the position of R.sub.p relative to R.sub.1 and R.sub.2.
a.sub.1.sup.2 =R.sub.b1.sup.2 =R.sub.bp.sup.2 +R.sub.p1.sup.2.
a.sub.2.sup.2 =R.sub.b2.sup.2 =R.sub.bp.sup.2 +R.sub.p2.sup.2.
The difference of these two equations is
a.sub.1.sup.2 a.sub.2.sup.2 =R.sub.p1.sup.2 R.sub.p2.sup.2.
We note that R.sub.p2 =R.sub.p1 R.sub.21 ; therefore,
R.sub.p2.sup.2 =R.sub.p1.sup.2 2R.sub.21.R.sub.p1 +R.sub.21.sup.2.
Combining these two equations gives:
a.sub.1.sup.2 a.sub.2.sup.2 =2R.sub.21.R.sub.p1 R.sub.21.sup.2.
or equivalently,
R.sub.21.R.sub.p1 =(R.sub.21.sup.2 +a.sub.1.sup.2 a.sub.2.sup.2)/2.
This means R.sub.21.R.sub.p1 can be either positive or negative, depending on the sign of R.sub.21.sup.2 +a.sub.1.sup.2 a.sub.2.sup.2. For example, if R.sub.1 were inside the sphere denoted .vertline.R.sub.02.vertline.=a.sub.2 and a.sub.1 were much less than a.sub.2, then R.sub.21.R.sub.p1 <0. Whether R.sub.21.R.sub.p1 is positive or negative, R.sub.21.times.R.sub.p1 =0, meaning that the two vectors are always either parallel or antiparallel. In both cases, the vector extending from R.sub.1 to the center of the circle of intersection is defined by the following equation: ##EQU1##
This equation defining R.sub.p1 differs from those which preceded it. The earlier equations were general mathematical relationships and partial clues pointing toward a solution, but because this equation contains only known (measured or previously computed) quantities on the right hand side and a single variable on the left, it represents not just a mathematical relationship between a set of known and unknown quantities, but rather, an explicit procedure (neither implicit, recursive, nor iterative) for calculating the left hand unknown in terms of a set of known quantities. Explicit procedural equations such as this will be enclosed in boxes to distinguish them from other equations, which only describe how they are derived. The set of all these explicit, procedural equations evaluated in the specified order, constitute the procedure which determines the satellite position R.sub.0. Now that we know R.sub.p1, we may compute R.sub.bp.sup.2 from the earlier equation for a.sub.1.sup.2 :
R.sub.bp.sup.2 =a.sub.1.sup.2 R.sub.p1.sup.2.
As mentioned above, R.sub.b is only one possible solution for R.sub.0, but we have not yet included the a.sub.3 data in our calculation. Without that additional information, we could not distinguish between R.sub.b and any other point on the circle of intersection lying above or below the plane of FIG. 3. R.sub.bp.sup.2 is the square of the radius of the circle of intersection, which is also equal to R.sub.0p.sup.2. ##EQU2##
At this step of the calculation, we know three key things about the circle of intersection: first, the position of its center, R.sub.p =R.sub.p1 +R.sub.1 ; second, its radius, .vertline.R.sub.0p.vertline.; and third, its orientation in three dimensional space, R.sub.21.R.sub.0p =0. We still must determine the angular position of R.sub.0 on the circle of intersection in order to obtain the full vector solution: R.sub.0 =R.sub.0p +R.sub.p. Determining that angle requires a third measurement.
Constructing the Coordinate System
In order to define the angular position of R.sub.0, we define a set of coordinate axes relative to which the angle is measured. We denote the line connecting R.sub.1 and R.sub.2 as the xaxis, and define the xdirected unit vector as ##EQU3##
(Dividing any vector by its own length produces a unit vector, that is, a vector with length equal to one.) This x unit vector (U.sub.x) is the "unit normal vector" of the circle of intersection, meaning that it is perpendicular to the plane containing that circle.
The cross product of two parallel (or antiparallel) vectors is zero. If R.sub.3 does not lie on the xaxis, then R.sub.31 is not parallel (or antiparallel) to R.sub.21 and U.sub.x. In that case, R.sub.z1 =U.sub.x.times.R.sub.31 is nonzero and is normal to the plane containing the RNs. If R.sub.3 lies on or very near to the xaxis, then .vertline.R.sub.z1.vertline.is zero or very close to zero, and the position computation becomes inaccurate. This is true regardless of the computational algorithm used, but this algorithm detects this condition explicitly, so that appropriate action may be taken, and the results of the computation will not be misinterpreted.
It is sometimes convenient to express R.sub.z1 in reference to R.sub.p rather than R.sub.1. This can be done because
R.sub.31 =R.sub.3p +R.sub.p1
R.sub.z1 =U.sub.x.times.R.sub.3p +U.sub.x.times.R.sub.p1
U.sub.x is parallel to R.sub.21 and R.sub.p1 ; therefore, U.sub.x.times.R.sub.p1 =0. We define the remaining two orthonormal vectors (U.sub.y and U.sub.z) of this coordinate system as follows: ##EQU4##
There is no need to divide U.sub.y by its own length, as we did for U.sub.x and U.sub.z, because U.sub.x and U.sub.z are perpendicular, and both are of length one; therefore, their cross product will be of length one.
Intersection with Sphere 3
The measurement a.sub.3 defines another sphere upon which R.sub.0 must lie. If point 3 is not on the xaxis, then this sphere will cross the previously defined circle of intersection in two points, one above and one below the xy plane, as shown in FIG. 4. We denote these crossing points R.sub.d and R.sub.e, and note that R.sub.de is parallel to U.sub.z. Denote as R.sub.q the point of intersection between R.sub.de and the xy plane (which is the plane containing U.sub.x and U.sub.y, and also the plane depicted in FIGS. 2 and 3).
We now develop a procedure for locating R.sub.d and R.sub.e. The satellite position R.sub.0 is equal to either R.sub.d or R.sub.e. Since one of these will lie below the surface of the Earth, there is no difficulty in choosing the proper one as the last step of the computation.
As shown in FIG. 4, the vector R.sub.dq is parallel to U.sub.z. FIG. 5, which is similar to FIG. 3, but with point R.sub.q added, shows that R.sub.p1 is parallel to U.sub.x, and R.sub.qp is parallel to U.sub.y. The sum of these three vectors R.sub.d1 =R.sub.dq +R.sub.qp+R.sub.p1 defines the position of point R.sub.d relative to point R.sub.1, and likewise for R.sub.e.
The Pythagorean theorem for the right triangle dpq decomposes R.sub.dp into two components which are aligned so that R.sub.dq is parallel to U.sub.z (see FIG. 4), and R.sub.qp is parallel to U.sub.y (see FIG. 5).
R.sub.dp.sup.2 =R.sub.dq.sup.2 +R.sub.qp.sup.2.
The decomposition of R.sub.d3 is similar; however, there is an added wrinkle. U.sub.z is perpendicular to R.sub.q3, as shown in FIG. 4; therefore, R.sub.q3 has no z component; however, R.sub.q3 may contain nonzero components in both the x and y directions. This is in contrast to R.sub.qp which has only a y component, as shown in the previous paragraph.
a.sub.3.sup.2 =R.sub.d3.sup.2 =R.sub.dq.sup.2 +R.sub.q3.sup.2.
The difference of these two equations is
a.sub.3.sup.2 R.sub.dp.sup.2 =R.sub.q3.sup.2 R.sub.qp.sup.2.
We note that R.sub.q3 =R.sub.qp +R.sub.p3 ; therefore,
R.sub.q3.sup.2 =R.sub.qp.sup.2 +R.sub.p3.sup.2 +2R.sub.qp.R.sub.p3.
Combining these last two equations yields:
a.sub.3.sup.2 R.sub.dp.sup.2 =R.sub.p3.sup.2 +2R.sub.qp.R.sub.p3.
Note that R.sub.d is one of the possible solutions for R.sub.0 ; therefore, R.sub.dp.sup.2 =R.sub.0p.sup.2. This is the radius of the circle of intersection, which we calculated previously. Substituting this into the previous equation yields:
R.sub.qp.R.sub.p3 =(a.sub.3.sup.2 R.sub.0p.sup.2 R.sub.p3.sup.2)/2.
As in a similar equation in the two sphere intersection procedure, this dot product, R.sub.qp.R.sub.p3 can be either positive or negative. The same procedure applies to both: ##EQU5##
Note that U.sub.y.R.sub.p3 =U.sub.y.R.sub.31. This quantity appears in the denominator of the equation above; hence, if it is zero, the algorithm breaks down. This is not an artifact of this particular algebraic derivation, but a fundamental, geometrical constraint. If U.sub.y.R.sub.31 =0 then point 3 lies on the xaxis. In that case we can determine only that R.sub.0 lies on a circle whose unit normal vector is U.sub.x. We can determine the radius of that circle, .vertline.R.sub.0p.vertline., and the position of its center, R.sub.p, but not the angular position of R.sub.0 around that circle.
The third RN must lie far enough from the xaxis to break the cylindrical symmetry, and resolve the angular position of R.sub.d and R.sub.e, one of which is equal to R.sub.0. The accuracy of the position computation degrades as the third RN approaches the xaxis. We have applied this analysis to actual orbital configurations of proposed systems in order to assess the accuracy of the algorithm in the face of measurement errors and suboptimal RN alignments. In almost all cases, U.sub.y.R.sub.31 is large enough, and the algorithm can continue as follows: The position of point q relative to the center of the Earth is ##EQU6##
Knowing R.sub.qp, we can use the Pythagorean theorem for triangle dpq (see FIG. 4) to obtain R.sub.dp.sup.2 =R.sub.dq.sup.2 +R.sub.qp.sup.2, which leads to the two candidate solutions for the satellite position: ##EQU7##
Notice that the two solutions are located symmetrically above and below the plane containing the RNs. Because the RNs are all on the surface of the Earth, one of these solutions will lie either inside, or on the opposite side of the Earth. In virtually all cases, only one of the solutions, R.sub.d or R.sub.e, will be within the field of view of all of the RNs which supplied the data; therefore, the algorithm can reject the spurious solution, and identify the other as R.sub.0.
Intersection of Sphere and Concentric Cone
If only two RNs are in view of a satellite, and provided that the satellite is not stationary with respect to the Earth, then measurements of range and range rate from a single RN can localize the satellite's possible positions to a circle of intersection, just as does the two sphere intersection discussed above. In this case, the circle of intersection corresponds to the intersection of a sphere with a cone whose vertex is located at the center of the sphere. After this circle of intersection is defined by a new procedure, one can apply the same spherecircle intersection procedure as is used with the Three Sphere positioning algorithm to solve for the satellite's position.
In the Sphere and Concentric Cone algorithm, the range measurement again constrains possible solution points to the surface of a sphere. The new facet is that the range rate measurement further limits solution points to a cone whose vertex lies at the center of that same sphere. Because the cone is concentric with the sphere, their curve of intersection is a circle. Such a circle is analogous to a line of constant latitude on a globe of the Earth; however, its orientation is defined by the motion of the satellite, and not fixed to the Earth.
If the satellite were in uniform, linear motion relative to the RN, then the relative velocity vector (V.sub.01) divided by its own length would define the axis of the cone (which is U.sub.x), and the ratio of radial to total relative velocity would determine the cone's angle=arccos(v.sub.r /.vertline.V.sub.01.vertline.); however, the actual situation is a bit more complicated. The satellite of interest in the present invention, and also the RNs on the surface of the Earth are in nearly uniform, circular motion relative to inertial space.
Despite the added complexity in this case, we still obtain a closed form solution. The generalization to elliptical satellite orbits is a trivial extension of this algorithm, and while not discussed in detail, should be considered included in the present invention. Both the satellite and the RNs are moving along trajectories of constant altitude. In the case of the RNs, this is due to the rotation of the surface of the Earth. One can express the velocity in any circular motion as V.sub.j =W.sub.j.times.R.sub.j ; therefore,
V.sub.01 =W.sub.0.times.R.sub.0 W.sub.1.times.R.sub.1,
where W.sub.j is the vector angular velocity. For RNs fixed to the surface of the Earth, W.sub.1 is identical to the Earth's angular velocity. If the satellite is not geostationary, but is moving in an orbit which is approximately circular, then on the time scale of the calculation process, we consider W.sub.0 to be constant.
As described above, the system measures the distance between the primary RN and the satellite (a.sub.1) from the signal propagation time, and the radial velocity (v.sub.r) from the time history of range measurements and from the Doppler shift. These scalar measurements are related to the satellite's vector position and velocity as follows:
V.sub.01.R.sub.01 =v.sub.r a.sub.1.
The range and range rate measurements define V.sub.01.R.sub.01 ; however, since R.sub.0 is unknown, and is involved in the definitions of both R.sub.01 and V.sub.01, we must do some vector algebra to disentangle R.sub.0 from the cross product.
V.sub.01.R.sub.01 =(W.sub.0.times.R.sub.0 (W.sub.0 W.sub.01).times.R.sub.1).R.sub.01 =(W.sub.0.times.R.sub.01 +W.sub.01.times.R.sub.1).R.sub.01.
Note that W.sub.0.times.R.sub.01 is perpendicular to R.sub.01, while W.sub.01.times.R.sub.1 has components both perpendicular and parallel to R.sub.01. The dot product (W.sub.0.times.R.sub.01).R.sub.01 is zero; therefore,
V.sub.01.R.sub.01 =(W.sub.01.times.R.sub.1).R.sub.01.
We assume that the satellite's angular velocity (W.sub.0) is known to high accuracy, and that what is not known at the time of the measurement is its precise position (R.sub.0). The angular velocity of the Earth (W.sub.1) and the position of the primary RN (R.sub.1) are assumed to be known. This implies that the vector W.sub.01.times.R.sub.1 is known, and we use it to define the x unit vector (U.sub.x) in this Range and Range Rate algorithm. It plays the same role here that R.sub.21 does in the Three Sphere algorithm. ##EQU8##
The center of the circle of intersection is again denoted R.sub.p, and in this algorithm, ##EQU9##
As in the Three Sphere algorithm, ##EQU10##
Via these equations, i.e., range plus range rate measurements from a single RN, the position, orientation, and radius of the circle of intersection are provided without using a second RN. Once one has the circle of intersection, either from Two Sphere or Sphere and Concentric Cone method, then one can use the same method to intersect that circle with the final sphere to get the satellite's position. In the Range and Range Rate method, all references to the "third RN" in the final sphere intersection method actually refer to the second RN, since no third RN is required. After the assignments R.sub.3 =R.sub.2 and a.sub.3 =a.sub.2, the Range and Range Rate algorithm proceeds in the same way as the Three Sphere algorithm.
While the best modes for carrying out the invention have been described in detail, those familiar with the art to which this invention relates will recognize various alternative designs and embodiments for practicing the invention as defined by the following claims.