# Patent

## Title
Communications systems including adaptive antenna systems and methods for inter-system and intra-system interference reduction

## Abstract

Methods of forming a downlink beam in an adaptive antenna system of a communications system that may reduce inter-system and/or intra-system interference include receiving vectors of signals including signals transmitted by user terminals of the communications system and signals transmitted by transceivers of an independent communications system, obtaining spatial information for the user terminals of the communications system and the transceivers of the independent communications system, generating complex transmitting weights that form spatial nulls directed at the transceivers of the independent communications system, and transmitting a downlink communications signal using the complex transmitting weights. Related systems are also disclosed.

## Background
RELATED APPLICATIONS
This application claims the benefit of and priority to U.S. Provisional Patent Application No. 60/692,404 entitled “Systems And Methods Of Interference Reduction To Satellite Terminals Using Ancillary Terrestrial Components With Smart Antennas,” filed Jun. 21, 2005, U.S. Provisional Patent Application No. 60/707,418 entitled “Systems and Methods for Reducing Inter- and Intra-System Interference by Adaptive Beam Forming Antennas,” filed Aug. 11, 2005, and U.S. Provisional Patent Application No. 60/787,052 entitled “Adaptive Antenna Systems And Methods For Inter-System Interference Reduction,” filed Mar. 29, 2006, the disclosures of which are hereby incorporated herein by reference as if set forth in its entirety.


FIELD OF THE INVENTION
This invention relates to wireless communications systems and methods, and more particularly to antenna systems and methods for terrestrial and/or satellite wireless communications systems.
BACKGROUND
Satellite radiotelephone communications systems and methods are widely used for radiotelephone communications. Satellite radiotelephone communications systems and methods generally employ at least one space-based component, such as one or more satellites, that is/are configured to wirelessly communicate with a plurality of satellite radiotelephones.
A satellite radiotelephone communications system or method may utilize a single satellite antenna pattern (beam or cell) covering an entire service region served by the system. Alternatively or in combination with the above, in cellular satellite radiotelephone communications systems and methods, multiple satellite antenna patterns (beams or cells) are provided, each of which can serve a substantially distinct service region in an overall service region, to collectively provide service to the overall service region. Thus, a cellular architecture similar to that used in conventional terrestrial cellular radiotelephone systems and methods can be implemented in cellular satellite-based systems and methods. The satellite typically communicates with radiotelephones over a bidirectional communications pathway, with radiotelephone communications signals being communicated from the satellite to the radiotelephone over a downlink or forward link (also referred to as forward service link), and from the radiotelephone to the satellite over an uplink or return link (also referred to as return service link). In some cases, such as, for example, in broadcasting, the satellite may communicate information to one or more radioterminals unidirectionally.
The overall design and operation of cellular satellite radiotelephone systems and methods are well known to those having skill in the art, and need not be described further herein. Moreover, as used herein, the term “radiotelephone” includes cellular and/or satellite radiotelephones with or without a multi-line display; Personal Communications System (PCS) terminals that may combine a radiotelephone with data processing, facsimile and/or data communications capabilities; Personal Digital Assistants (PDA) that can include a radio frequency transceiver and/or a pager, Internet/Intranet access, Web browser, organizer, calendar and/or a global positioning system (GPS) receiver; and/or conventional laptop and/or palmtop computers or other appliances, which include a radio frequency transceiver. A radiotelephone also may be referred to herein as a “radioterminal”, a “mobile terminal” or simply as a “terminal”. As used herein, the term(s) “radioterminal,” “radiotelephone,” “mobile terminal” and/or “terminal” also include(s) any other radiating user device/equipment/source that may have time-varying or fixed geographic coordinates and/or may be portable, transportable, installed in a vehicle (aeronautical, maritime, or land-based) and/or situated and/or configured to operate locally and/or in a distributed fashion over one or more terrestrial and/or extra-terrestrial location(s). Furthermore, as used herein, the term “space-based component” or “space-based system” includes one or more satellites at any orbit (geostationary, substantially geostationary, medium earth orbit, low earth orbit, etc.) and/or one or more other objects and/or platforms (e.g., airplanes, balloons, unmanned vehicles, space crafts, missiles, etc.) that has/have a trajectory above the earth at any altitude.
Terrestrial networks can enhance cellular satellite radiotelephone system availability, efficiency and/or economic viability by terrestrially using/reusing at least some of the frequencies that are allocated to cellular satellite radiotelephone systems. In particular, it is known that it may be difficult for cellular satellite radiotelephone systems to reliably serve densely populated areas, because satellite signals may be blocked by high-rise structures and/or may not penetrate into buildings. As a result, satellite spectrum may be underutilized or unutilized in such areas. The terrestrial use/reuse of at least some of the satellite system frequencies can reduce or eliminate this potential problem.
Moreover, the capacity of an overall hybrid system, comprising space-based (i.e., satellite) and terrestrial communications capability, may be increased by the introduction of terrestrial frequency use/reuse of frequencies authorized for use by the space-based component, since terrestrial frequency use/reuse may be much denser than that of a satellite-only system. In fact, capacity may be enhanced where it may be mostly needed, i.e., in densely populated urban/industrial/commercial areas. As a result, the overall system may become more economically viable, as it may be able to serve more effectively and reliably a larger subscriber base.
One example of terrestrial reuse of satellite frequencies is described in U.S. Pat. No. 5,937,332 to co-inventor Karabinis entitled Satellite Telecommunications Repeaters and Retransmission Methods, the disclosure of which is hereby incorporated herein by reference in its entirety as if set forth fully herein. As described therein, satellite telecommunications repeaters are provided which receive, amplify, and locally retransmit the downlink/uplink signal received from a satellite/radioterminal thereby increasing an effective downlink/uplink margin in the vicinity of the satellite telecommunications repeater and allowing an increase in the penetration of uplink and downlink signals into buildings, foliage, transportation vehicles, and other objects which can reduce link margin. Both portable and non-portable repeaters are provided. See the abstract of U.S. Pat. No. 5,937,332.
Satellite radiotelephones for a satellite radiotelephone system or method having a terrestrial communications capability by terrestrially using/reusing at least some of the same satellite frequency band and using substantially the same air interface for both terrestrial and satellite communications may be cost effective and/or aesthetically appealing. Conventional dual band/dual mode radiotelephone alternatives, such as the well known Thuraya, Iridium and/or Globalstar dual mode satellite/terrestrial radiotelephones, duplicate some components (as a result of the different frequency bands and/or air interface protocols between satellite and terrestrial communications), which leads to increased cost, size and/or weight of the radiotelephone. See U.S. Pat. No. 6,052,560 to co-inventor Karabinis, entitled Satellite System Utilizing a Plurality of Air Interface Standards and Method Employing Same.
Satellite radioterminal communications systems and methods that may employ terrestrial use and/or reuse of satellite frequencies by an Ancillary Terrestrial Network (ATN) comprising at least one Ancillary Terrestrial Component (ATC) are also described in U.S. Pat. No. 6,684,057 to Karabinis, entitled Systems and Methods for Terrestrial Reuse of Cellular Satellite Frequency Spectrum; U.S. Pat. No. 6,785,543 to Karabinis, entitled Filters for Combined Radiotelephone/GPS Terminals; U.S. Pat. No. 6,856,787 to Karabinis, entitled Wireless Communications Systems and Methods Using Satellite-Linked Remote Terminal Interface Subsystems; U.S. Pat. No. 6,859,652 to Karabinis et al., entitled Integrated or Autonomous System and Method of Satellite-Terrestrial Frequency Reuse Using Signal Attenuation and/or Blockage, Dynamic Assignment of Frequencies and/or Hysteresis; and U.S. Pat. No. 6,879,829 to Dutta et al., entitled Systems and Methods for Handover Between Space Based and Terrestrial Radioterminal Communications, and For Monitoring Terrestrially Reused Satellite Frequencies At a Radioterminal to Reduce Potential Interference, and in U.S. Pat. Nos. 6,892,068, 6,937,857, 6,999,720 and 7,006,789; and Published U.S. Patent Application Nos. US 2003/0054761 to Karabinis, entitled Spatial Guardbands for Terrestrial Reuse of Satellite Frequencies; US 2003/0054814 to Karabinis et al., entitled Systems and Methods for Monitoring Terrestrially Reused Satellite Frequencies to Reduce Potential Interference; US 2003/0073436 to Karabinis et al., entitled Additional Systems and Methods for Monitoring Terrestrially Reused Satellite Frequencies to Reduce Potential Interference; US 2003/0054762 to Karabinis, entitled Multi-Band/Multi-Mode Satellite Radiotelephone Communications Systems and Methods; US 2002/0041575 to Karabinis et al., entitled Coordinated Satellite-Terrestrial Frequency Reuse; US 2003/0068978 to Karabinis et al., entitled Space-Based Network Architectures for Satellite Radiotelephone Systems; US 2003/0153308 to Karabinis, entitled Staggered Sectorization for Terrestrial Reuse of Satellite Frequencies; and US 2003/0054815 to Karabinis, entitled Methods and Systems for Modifying Satellite Antenna Cell Patterns In Response to Terrestrial Reuse of Satellite Frequencies, and in Published U.S. Patent Application Nos. 2004/0121727, 2004/0142660, 2004/0192395, 2004/0192200, 2004/0192293, 2004/0203742, 2004/0240525, 2005/0026606, 2005/0037749, 2005/0041619, 2005/0064813, 2005/0079816, 2005/0090256, 2005/0118948, 2005/0136836, 2005/0164700, 2005/0164701, 2005/0170834, 2005/0181786, 2005/0201449, 2005/0208890, 2005/0221757, 2005/0227618, 2005/0239457, 2005/0239403, 2005/0239404, 2005/0239399, 2005/0245192, 2005/0260984, 2005/0260947, 2005/0265273, 2005/00272369, 2005/0282542, 2005/0288011, 2006/0040613, 2006/040657 and 2006/0040659; all of which are assigned to the assignee of the present invention, the disclosures of all of which are hereby incorporated herein by reference in their entirety as if set forth fully herein.
Many of the above-cited patents and publications can reduce or eliminate intra-system interference that is caused by terrestrial use and/or reuse of satellite frequencies by an ancillary terrestrial network. However, inter-system interference also can be caused by terrestrial use and/or reuse of satellite frequencies by an ancillary terrestrial network and/or radioterminals communicating therewith, due to, for example, the potential production of out-of-band and/or out-of-channel emissions by a radioterminal that is transmitting over a satellite frequency band.
SUMMARY
Methods of forming a downlink beam in an adaptive antenna system of a communications system according to some embodiments of the invention include receiving vectors of signals including signals transmitted by user terminals of the communications system and signals transmitted by transceivers of an independent communications system, obtaining spatial information for the user terminals of the communications system and the transceivers of the independent communications system, generating complex transmitting weights that form spatial nulls directed at the transceivers of the independent communications system, and transmitting a downlink communications signal using the complex transmitting weights.
Receiving vectors of signals may include receiving vectors of signals using an antenna array including a number of antenna feed elements, and a respective one of the complex transmitting weights is associated with a corresponding one of the antenna feed elements.
Obtaining spatial information for the user terminals of the communications system and the transceivers of the independent communications system may include obtaining the spatial information using the received vectors of signals.
Obtaining spatial information for the user terminals of the communications system may include performing pseudo uplink beamforming, and obtaining spatial information for the transceivers of the independent communications system may include using blind spatial channel estimation. Obtaining spatial information may further include generating estimates of an array covariance matrix using the signals transmitted by the user terminals of the communications system and the signals transmitted by the transceivers of an independent communications system.
Generating estimates of the array covariance matrix may include obtaining estimates of spatial signatures of the user terminals of the communications system using a preamble or pilot signal reference.
Generating estimates of the array covariance matrix may include estimating a covariance matrix for the transceivers of the independent communications system, estimating a covariance matrix for the user terminals of the communications system, and deriving a composite covariance matrix by combining the covariance matrix for the transceivers of the independent communications system with the covariance matrix for the user terminals of the communications system.
The methods may further include performing frequency transformation of the covariance matrix for the transceivers of the independent communications system to transform a frequency of the covariance matrix for the transceivers of the independent communications system to a frequency of the communications system.
Generating complex transmitting weights may include generating a set of pseudo uplink beamforming weights by combining the estimates of the spatial signatures of user terminals of the communications system with the composite covariance matrix, and transmitting a downlink communications signal using the complex transmitting weights may include transmitting a downlink communications signal using the pseudo uplink beamforming weights transmitting weights.
Obtaining spatial information for the user terminals of the communications system and/or the transceivers of the independent communications system may include obtaining spatial signatures using an angle of arrival estimation algorithm.
Obtaining spatial information for the transceivers of the independent communications system may further include estimating angle of arrival information of the signals from the transceivers of the independent communications system using a maximum likelihood angle of arrival estimation algorithm, and generating spatial signatures of the transceivers of the independent communications system using the angle of arrival information.
The communications system may be configured to communicate with user terminals using a time division duplex communications protocol.
Obtaining spatial information for the transceivers of the independent communications system may include generating estimates of a covariance matrix using the signals transmitted by the transceivers of an independent communications system and decomposing the covariance matrix to obtain the principal eigenvector of the covariance matrix or the K eigenvectors corresponding to the K largest eigenvalues, where K represents the number of transceivers of the independent communications system.
Receiving signals transmitted by transceivers of an independent communications system may include receiving pilot or carrier signals transmitted by the transceivers of the independent communications system. In some embodiments, receiving signals transmitted by transceivers of an independent communications system may include receiving beacon signals transmitted by the transceivers of the independent communications system.
The methods may further include forming an uplink antenna pattern using the complex transmitting weights, receiving the beacon signals using the uplink antenna pattern, and measuring strength of the beacon signals received using the uplink antenna pattern, the strength of the beacon signals received using the uplink antenna pattern provides an error signal associated with the complex transmitting weights. The methods may further include adjusting the complex transmitting weights to reduce the strength of the beacon signals.
The adaptive antenna system may include a vertical uniform linear array of radiating elements, and generating complex transmitting weights that form spatial nulls directed at the transceivers of the independent communications system may include generating complex transmitting weights that form spatial nulls directed at angles of arrival associated with the transceivers of the independent communications system.
The adaptive antenna system may include a phased array antenna, and generating complex transmitting weights that form spatial nulls directed at the transceivers of the independent communications system may include generating complex transmitting weights that form spatial nulls directed at azimuth/elevation coordinates associated with the transceivers of the independent communications system.
The methods may include receiving the downlink communications signal at a user terminal of the communication system. The communications system comprises a terrestrial communications system.
Methods of generating a downlink transmission signal from a base transceiver system to a user terminal associated with the base transceiver system according to further embodiments of the invention include obtaining location information for a transceiver of an independent communications system not associated with the base transceiver system, generating complex transmitting weights for an adaptive antenna associated with the base transceiver system that form a spatial null directed at the transceiver of the independent communications system, and transmitting a downlink communications signal over the adaptive antenna system to the user terminal using the complex transmitting weights.
The adaptive antenna may include a plurality of vertical uniform linear array antennas, and generating complex transmitting weights may include generating complex transmitting weights that form a spatial null directed at an elevation associated with the transceiver of the independent communications system.
Obtaining location information for the transceiver of the independent communications system may include receiving a beacon signal from the transceiver of the independent communications system. The beacon signal may include an unmodulated carrier signal, a control channel signal and/or a traffic channel signal.
Generating the complex transmitting weights may include forming an uplink antenna pattern using the complex transmitting weights, receiving the beacon signal using the uplink antenna pattern, and measuring strength of the beacon signal received using the uplink antenna pattern, the strength of the beacon signal received using the uplink antenna pattern provides an error signal associated with the complex transmitting weights. The methods may further include adjusting the complex transmitting weights to reduce the strength of the beacon signals.
Transmitting the downlink communications signal may include transmitting a signal that is Left Hand Circularly Polarized (LHCP), Right Hand Circularly Polarized (RHCP), or Linearly Polarized (LP).
The methods may include receiving the downlink communications signal at a user terminal of the communication system. The communications system comprises a terrestrial communications system.
Some embodiments of the invention provide a base transceiver system of a communications system, including an antenna array configured to receive vectors of signals including signals transmitted by user terminals of the communications system and signals transmitted by transceivers of an independent communications system, and an electronics system coupled to the antenna array. The electronics system is configured to obtain spatial information for the user terminals of the communications system and the transceivers of the independent communications system, to generate complex transmitting weights that form spatial nulls directed at the transceivers of the independent communications system, and to transmit a downlink communications signal over the antenna array using the complex transmitting weights. The antenna array may include an array of antenna feed elements, and respective ones of the complex transmitting weights are associated with corresponding ones of the antenna feed elements. The electronics system may be configured to obtain spatial information for the user terminals of the communications system and the transceivers of the independent communications system using the received vectors of signals.
The electronics system may be configured to obtain spatial information for the user terminals of the communications system by performing uplink beamforming, and the electronics system may be configured to obtain spatial information for the transceivers of the independent communications system using blind spatial channel estimation.
The electronics system may be configured to generate estimates of an array covariance matrix using the signals transmitted by the user terminals of the communications system and the signals transmitted by the transceivers of an independent communications system.
The electronics system may be configured to obtain estimates of spatial signatures of the user terminals of the communications system using a preamble or pilot signal reference.
The electronics system may be configured to generate estimates of the array covariance matrix by estimating a covariance matrix for the transceivers of the independent communications system, estimating a covariance matrix for the user terminals of the communications system, and deriving a composite covariance matrix by combining the covariance matrix for the transceivers of the independent communications system with the covariance matrix for the user terminals of the communications system.
The electronics system may be further configured to transform a frequency of the covariance matrix for the transceivers of the independent communications system to a frequency of the communications system.
The electronics system may be configured to generate a set of pseudo uplink beamforming weights by combining the estimates of the spatial signatures of user terminals of the communications system with the composite covariance matrix, and the electronics system may be further configured to transmit the downlink communications signal using the pseudo uplink beamforming weights transmitting weights.
The electronics system may be configured to obtain spatial signatures for the user terminals of the communications system and/or the transceivers of the independent communications system using an angle of arrival estimation algorithm.
The electronics system may be configured to estimate angle of arrival information of the signals from the transceivers of the independent communications system using a maximum likelihood angle of arrival estimation algorithm, and to generate spatial signatures of the transceivers of the independent communications system using the angle of arrival information. The communications system may be configured to communicate with user terminals using a time division duplex communications protocol.
The electronics system may be configured to generate estimates of a covariance matrix using the signals transmitted by the transceivers of an independent communications system and to decompose the covariance matrix to obtain the principal eigenvector of the covariance matrix or the K eigenvectors corresponding to the K largest eigenvalues, where K represents the number of transceivers of the independent communications system.
The antenna array may be configured to receive pilot or carrier signals transmitted by the transceivers of the independent communications system.
The antenna array may be configured to receive beacon signals transmitted by the transceivers of the independent communications system. The electronics system may be configured to form an uplink antenna pattern using the complex transmitting weights, to receive the beacon signals using the uplink antenna pattern, and to measure strength of the beacon signals received using the uplink antenna pattern, the strength of the beacon signals received using the uplink antenna pattern provides an error signal associated with the complex transmitting weights. The electronics system may be configured to adjust the complex transmitting weights to reduce the strength of the beacon signals.
The antenna array may include a vertical uniform linear array of radiating elements, and the electronics system may be configured to generate complex transmitting weights that form spatial nulls directed at angles of arrival associated with the transceivers of the independent communications system.
The antenna array may include a phased array antenna, and the electronics system may be configured to generate complex transmitting weights that form spatial nulls directed at azimuth/elevation coordinates associated with the transceivers of the independent communications system.
A base transceiver system of a communication system according to some embodiments of the invention includes an antenna array configured to receive vectors of signals including signals transmitted by user terminals of the communications system and signals transmitted by a transceiver of an independent communications system, and an electronics system coupled to the antenna array. The electronics system is configured to obtain location information for a transceiver of the independent communications system, to generate complex transmitting weights for the antenna array that form a spatial null directed at the transceiver of the independent communications system, and to transmit a downlink communications signal over the adaptive antenna system to the user terminal using the complex transmitting weights.
The adaptive antenna may include a plurality of vertical uniform linear array antennas, and the electronics system may be configured to generate complex transmitting weights that form a spatial null directed at an elevation associated with the transceiver of the independent communications system.
The antenna array may be configured to receiving a beacon signal from the transceiver of the independent communications system. The beacon signal may include an unmodulated carrier signal, a control channel signal and/or a traffic channel signal.
The electronics system may be configured to form an uplink antenna pattern using the complex transmitting weights, receive the beacon signal using the uplink antenna pattern, and to measure strength of the beacon signal received using the uplink antenna pattern, the strength of the beacon signal received using the uplink antenna pattern provides an error signal associated with the complex transmitting weights. The electronics system may be configured to adjust the complex transmitting weights to reduce the strength of the beacon signals.
The antenna array may be configured to transmit the downlink communications signal may include transmitting a signal that is Left Hand Circularly Polarized (LHCP), Right Hand Circularly Polarized (RHCP), or Linearly Polarized (LP).
A communications system according to some embodiments of the invention includes user terminals configured to generate communication signals in a first frequency band, and a base transceiver system configured as described above including an antenna array configured to receive vectors of signals including signals transmitted by the user terminals in the first frequency band and signals transmitted by transceivers of an independent communications system in a second frequency band and an electronics system coupled to the antenna array.



BRIEF DESCRIPTION OF THE DRAWINGS
The accompanying drawings, which are included to provide a further understanding of the invention and are incorporated in and constitute a part of this application, illustrate certain embodiment(s) of the invention. In the drawings:
FIG. 1 is a schematic diagram of a cellular satellite communications system and methods according to embodiments of the invention;
FIG. 2 is a schematic block diagram of an adaptive antenna array system at a base transceiver system (BTS) including 6 Vertical Uniform Linear Arrays;
FIG. 3 is a schematic block diagram of an architecture of a BTS with the adaptive antenna array system;
FIG. 4 is a flowchart of an adaptive downlink beamforming algorithm according to some embodiments of the invention;
FIG. 5 is a flowchart of an adaptive downlink beamforming algorithm according to further embodiments of the invention;
FIG. 6 is a plot of a downlink beam pattern according to some embodiments of the invention when a terrestrial communications system (TCS) user is located at −5° elevation and an ICS satellite is located at 5° elevation (−20 dB constraint for the satellite);
FIG. 7 is a plot of a downlink beam pattern according to some embodiments of the invention when a TCS user is located at −5° elevation and ICS satellites are located at 5 & 18° elevations respectively (−20 dB constraint for the satellites);
FIG. 8 is a plot of a downlink beam pattern according to some embodiments of the invention when a TCS user is located at −5° elevation and ICS satellites are located at 7, 17 & 20° elevations respectively (−20 dB constraint for the satellites);
FIG. 9 is a plot of a downlink beam pattern according to some embodiments of the invention when a TCS user is located at −5° elevation and ICS satellites are located at 7, 17 & 20° elevations respectively (−28 dB constraint for the satellites);
FIG. 10 is a plot of a downlink beam pattern according to some embodiments of the invention when a TCS user is located at −5° elevation and ICS satellites are located at 7, 17 & 20° elevations respectively (−34 dB constraint for the satellites);
FIG. 11 is a plot of average gains of formed beams according to some embodiments of the invention toward a TCS user and ICS satellites when the TCS user is located at −5° elevation and satellites are uniformly randomly distributed in [0°, 30°] elevation; and
FIG. 12 is a plot of average gains of formed beam according to some embodiments of the invention toward TDD user and ICS satellites when AOA of TDD user varies from −1 to −5° elevation and satellites are uniformly randomly distributed in [0°, 10°] elevation (Algorithm I).



DETAILED DESCRIPTION OF EMBODIMENTS OF THE INVENTION
Embodiments of the present invention now will be described more fully hereinafter with reference to the accompanying drawings, in which embodiments of the invention are shown. This invention may, however, be embodied in many different forms and should not be construed as limited to the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the scope of the invention to those skilled in the art. Like numbers refer to like elements throughout.
It will be understood that, although the terms first, second, etc. may be used herein to describe various elements, these elements should not be limited by these terms. These terms are only used to distinguish one element from another. For example, a first element could be termed a second element, and, similarly, a second element could be termed a first element, without departing from the scope of the present invention. As used herein, the term “and/or” includes any and all combinations of one or more of the associated listed items.
The terminology used herein is for the purpose of describing particular embodiments only and is not intended to be limiting of the invention. As used herein, the singular forms “a”, “an” and “the” are intended to include the plural forms as well, unless the context clearly indicates otherwise. It will be further understood that the terms “comprises” “comprising,” “includes” and/or “including” when used herein, specify the presence of stated features, integers, steps, operations, elements, and/or components, but do not preclude the presence or addition of one or more other features, integers, steps, operations, elements, components, and/or groups thereof.
Unless otherwise defined, all terms (including technical and scientific terms) used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs. It will be further understood that terms used herein should be interpreted as having a meaning that is consistent with their meaning in the context of this specification and the relevant art and will not be interpreted in an idealized or overly formal sense unless expressly so defined herein.
As will be appreciated by one of skill in the art, the present invention may be embodied as a method, data processing system, and/or computer program product. Accordingly, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects all generally referred to herein as a “circuit” or “module.” Furthermore, the present invention may take the form of a computer program product on a computer usable storage medium having computer usable program code embodied in the medium. Any suitable computer readable medium may be utilized including hard disks, CD ROMs, optical storage devices, a transmission media such as those supporting the Internet or an intranet, or magnetic storage devices.
The present invention is described below with reference to flowchart illustrations and/or block diagrams of methods, systems and computer program products according to embodiments of the invention. It will be understood that each block of the flowchart illustrations and/or block diagrams, and combinations of blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions/acts specified in the flowchart and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer readable memory produce an article of manufacture including instruction means which implement the function/act specified in the flowchart and/or block diagram block or blocks.
The computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions/acts specified in the flowchart and/or block diagram block or blocks.
It is to be understood that the functions/acts noted in the blocks may occur out of the order noted in the operational illustrations. For example, two blocks shown in succession may in fact be executed substantially concurrently or the blocks may sometimes be executed in the reverse order, depending upon the functionality/acts involved. Although some of the diagrams include arrows on communication paths to show a primary direction of communication, it is to be understood that communication may occur in the opposite direction to the depicted arrows.
Some embodiments of the present invention may arise from realization that a transceiver of an Independent Communications System (ICS), such as, for example, a satellite or user terminal of a satellite communications systems, may operate in an adjacent (or proximate) frequency band to a frequency band used by one or more Terrestrial Cellular Systems (TCS). For example, a user device of an ICS may be physically proximate to (e.g., within a cell of) the TCS, and/or a satellite of an ICS may be located at an azimuth/elevation relative to a Base Transceiver System (BTS) of the TCS such that emissions of the BTS may be received by the satellite. As such, the ICS transceiver may be subjected to overload/intermodulation interference from a transmitter of the TCS. A desire to reduce or avoid interference, such as, for example, inter-modulation and/or overload interference, from a TCS transmitter (e.g., base transceiver system) to an ICS transceiver, may be used to impose an additional constraint to adaptive beam forming that may be used by the TCS transmitter.
In some embodiments of the present invention, systems and methods are provided that may reduce or minimize inter-system interference from a Terrestrial Cellular System (TCS) to an Independent Communications System (ICS) that may be operating in a frequency band that is proximate to a frequency band used by the TCS. The technique also may reduce or minimize the intra-system interference. A Base Transceiver System (BS) of the TCS may be configured with an adaptive antenna array to thereby receive vectors of ICS and TCS signals, respectively. Estimates of array covariance matrices may be derived from the received vectors of an ICS signal(s) and TCS signal(s), respectively. An estimate of an ICS covariance matrix (or spatial signature) may be transformed to the TCS frequency band. The frequency transformed ICS covariance matrix or spatial signature, along with the TCS covariance matrix or a spatial signature, may be exploited to form a composite covariance matrix or spatial signature matrix that contains spatial information associated with both TCS user terminals and ICS transceivers. With an a priori known preamble and/or pilot signal for each of the TCS user terminals and the composite TCS/ICS covariance or spatial signature matrix from uplink processing of received signals, an adaptive beam-former may determine optimum transmitting weights that may be used to form a beam that substantially reduces or minimizes intra-system interference among the TCS user terminals while substantially forming a “null” toward the ICS. By using the composite covariance matrix, the adaptive beamformer introduces spatial information relative to ICS transceivers, including for example, ICS user terminals and/or ICS satellites, as an additional constraint on top of intra-system multi-user spatial constraints when forming a beam on a downlink.
FIG. 1 is a schematic diagram of cellular satellite communications systems and methods according to embodiments of the invention. As shown in FIG. 1, these cellular satellite communications systems and methods 100 include a Space-Based Component (SBC) 110, such as one or more geostationary or non-geostationary orbiting satellites. The space-based component 110 may be configured to selectively use/reuse geographically a set of frequencies and to transmit wireless communications signals via a communication link 114 to a plurality of wireless terminals, only one of which is illustrated in FIG. 1 (terminal 120a), in a satellite footprint including one or more satellite cells 130-130′, over one or more satellite forward service link (downlink) frequencies fD. The space-based component 110 may also be configured to receive wireless communications from a plurality of wireless terminals, such as wireless terminal 120a in the satellite cell 130, over one or more satellite return service link (uplink) frequencies fU.
A cellular satellite communications system 100 may include a terrestrial communications system (TCS) including at least one Base Transceiver System (BTS) 140, which may include an antenna 140a and an electronics system 140b. The BTS 140 is configured to receive wireless communications signals from, for example, at least one wireless terminal 120b via a communication link 116. In some embodiments, the TCS is an ancillary terrestrial network (ATN) of the cellular satellite communications system 100 that uses/reuses satellite frequencies terrestrially. A cellular satellite communications system may provide users with increased mobility/functionality and/or may provide increased frequency reuse. However, it will be appreciated that embodiments of the invention can be used to decrease inter-system interference and/or intra-system interference in systems that are terrestrial only, space-based only, and/or include combinations of terrestrial and space-based components.
As noted above, a user terminal 120a of the cellular satellite communications system 100 may communicate with the SBC 110 using an uplink frequency fU within the satellite frequency band. In order to provide increased frequency utilization in a cellular satellite communications system 100, a user terminal 120b of the cellular satellite communications system 100 may communicate with the terrestrial BTS 140 using an uplink frequency f′U within the satellite frequency band. The frequency f′U may be the same as an uplink or downlink frequency used for communicating with the space-based component (SBC) 110 in the satellite cell 130 in which the wireless terminal 120b is located and/or in an adjacent or remotely-located satellite cell 130. Thus, as illustrated in FIG. 1, the wireless terminal 120a may be communicating with the space-based component 110 using a frequency in the satellite frequency band while the wireless terminal 120b may be communicating with the ancillary terrestrial component 140, also using a frequency in the satellite frequency band. The space-based component 110 may also undesirably receive a component of the wireless communications from the wireless terminal 120b and/or the BTS 140 in the satellite cell 130 as interference. In addition, the space based component 110 may receive, as interference, a component of wireless communications from a wireless terminal and/or BTS (not shown) located in a different satellite cell over a satellite frequency that may be the same as (and/or overlapping with) fU and/or f′U.
This interference may generally be strongest when the transmitted signal uses the same carrier frequency as the cell in question (e.g., f′U=fU), because, in that case, the same return link frequency would be used for space-based component and ancillary terrestrial component communications and, if used over the same satellite cell, no substantial spatial discrimination between satellite cells would appear to exist to reduce a level of interference. Even with spatial separation, however, interference may impair the signal from the first wireless terminal 120a. 
Still referring to FIG. 1, embodiments of satellite communications systems/methods 100 can include at least one satellite gateway 160 that can include an antenna 160a and an electronics system 160b. The satellite gateway 160 may be connected to other networks 162, including terrestrial and/or other wired and/or wireless communications networks such as, for example, a public switched telephone network and/or the Internet. The satellite gateway 160 communicates with the space-based component 110 over a satellite feeder link 112. The satellite gateway 160 may also be configured to communicate with ancillary terrestrial components 140 in the ancillary terrestrial network, generally over a terrestrial link 142.
The satellite gateway 160 may communicate with the BTS 140 via a communications link 142, which may be used to coordinate communications of the BTS 140 and the SBC 110 in order to reduce interference between the ancillary terrestrial network and the SBC 110.
Still referring to FIG. 1, a satellite 210 of an independent communications system (ICS) 200 may also communicate via a communications link 212 with an ICS user terminal 220 of the ICS 200 using an uplink and/or downlink frequency that is the same as, adjacent to or overlapping with a frequency used by the BTS 140. The satellite 210 may be part of a separate satellite communications system (i.e. one that is not under common control and/or coordination with the cellular satellite communications system 100). Accordingly, the use of satellite frequencies may not be coordinated between the external satellite 210 and the BTS 140. For example, the external satellite 210 may be a satellite in the existing Inmarsat satellite communications system, which includes eleven geostationary communications satellites, and the ICS user terminal 220 may include, for example, a fixed or mobile satellite terminal.
For a BTS operating in time division duplex (TDD) mode, it may be desirable to require that an antenna of the BTS that transmits using a frequency used by an ICS should provide at least 20 dB of “aggregate” discrimination relative to its bore-sight EIRP (Effective Isotropically Radiated Power) at and above the horizontal plane of the antenna. In fact, such requirements have been implemented by regulation with respect to ICS uplink frequencies used by a BTS. That is, under certain regulatory provisions, the power radiated by the BTS antenna at and above the horizontal plane of the antenna in the direction of an ICS satellite must be at least 20 dB lower than the EIRP of the maximum signal transmitted by the antenna. EIRP refers to the amount of power that would have to be emitted by an isotropic antenna (i.e., one that evenly distributes power in all directions) to produce the peak power density observed in the direction of maximum antenna gain. The aggregate discrimination may come from true antenna gain discrimination and/or signal blockage due to terrain/buildings.
Furthermore, it may be desirable to reduce the EIRP in a direction of an ICS user terminal 220 in order to reduce or avoid interference, such as intermodulation and/or overload interference at the ICS user terminal 220 from a transmitter (such as a BTS) in the TCS.
Some embodiments of the invention may improve true antenna discrimination (as opposed to discrimination provided by terrain/building blockage) by using an adaptive antenna system (AAS) at a BTS. By exploiting the spatial properties of the communications signals of an ICS satellite communications system and BTS end user signals, an adaptive antenna system according to some embodiments of the invention may adaptively form spatial “null(s)” directed toward ICS satellite(s) and/or ICS user terminals, while forming a beam toward a desired end user of the TCS. The term “null” or “spatial null” is used herein in a general sense to refer to directions/areas of reduced signal strength response, and may correspond to, but does not necessarily indicate, a zero or minimum signal strength response at a particular direction. The performance of an adaptive antenna system according to some embodiments of the invention has been evaluated through simulations under various scenarios. The simulation results show that providing an adaptive antenna system at a BTS may provide more than 20 dB of true antenna discrimination irrespective of any signal blockage due to buildings or terrain.
Some embodiments of the invention may arise from realization that for a BTS utilizing a time division duplex (TDD) multiple access scheme, the BTS may transmit a downlink signal using the same frequency band as the one used by the external satellite uplink, such as an ICS uplink. Although only sidelobes of the ICS transceivers may pick up BTS interference, the aggregate interference from a large number of BTSs may still cause an interference problem at the ICS transceiver, e.g. on the uplink of the external satellite. To alleviate this problem, it may be desirable for any BTS antenna transmitting on an uplink frequency of the external satellite to provide at least 20 dB of “aggregate” discrimination (relative to its bore-sight EIRP) at the horizontal plane and above the horizontal plane, from which the external satellite(s) may pick up BTS interference. While possible terrain/building blockage may provide some cushion for the required antenna discrimination, it is difficult to anticipate the effectiveness of such blockage. An antenna array system for a BTS according to some embodiments of the invention may reduce and/or minimize inter-system interference from the BTS to ICS transceivers that may be operating in the same frequency band as the one used by the terrestrial BTS, regardless of the presence or absence of terrain/building blockage.
Adaptive antenna array processing is a technique that may be used to exploit the spatial dimension by adaptively forming a maximum, near maximum, or relatively high power beam (antenna pattern) toward a desired user, and nulls toward other unintended users or systems, such as ICS transceivers. Many algorithms have appeared in the literature for beamforming and estimating signal parameters based on measurements of outputs of an antenna sensor array. Spatial filtering techniques represent early approaches to performing space-time processing of data at the output of an antenna array. Conventional beamformer and Capon's methods are examples of spatial filtering techniques. A conventional beam-former is a spatial filter matched to an impinging signal. Capon's beam-former attempts to reduce or minimize the power received from the noise and any signal coming from directions other than the direction of interest, while maintaining a fixed gain in the direction of interest. This method may be useful in situations where information about the statistical properties of data is not available.
Adaptive beam-forming algorithms can be generally classified into two categories: 1) training or decision-directed adaptive algorithms, and 2) blind adaptive algorithms. Depending on a performance criterion, different beamforming techniques lead to different adaptive algorithms. Typically, beam-forming techniques based on Minimum Mean Squared Error (MMSE) or Least Squares (LS) criteria lead to the Least Mean Squares (LMS) and the Recursive Least Squares (RLS) adaptive algorithms. The MMSE or LS algorithms require a reference signal for the desired user, which could be based either on a training sequence or on decisions about received symbols. The beamforming techniques based on Linear Constrained Minimum Variance (LCMV) and Maximum Signal-to-Noise Ratio (MSNR) criteria do not require a reference signal to determine optimal weight(s). See, e.g., B. D. Van Veen and K. M. Buckley, “Beamforming: A Versatile Approach to Spatial Filtering,” IEEE Acoustics, Speech, and Signal Processing Magazine, April 1988. However, both techniques need to have knowledge of the Angle of Arrival (AOA) of the desired signal. High-resolution AOA estimation algorithms are typically based on model-based parametric approaches. These algorithms assume some knowledge about the statistical properties of the data (such as second order ergodicity).
Another model-based approach is the Maximum Likelihood (ML) method. The ML method estimates are those values of the unknown parameters that maximize a likelihood function. Two different ML approaches may be used, namely, deterministic and stochastic, depending on the model assumption regarding the signal waveform. The unknown, deterministic signal model (which is the case for mobile communications) leads to the deterministic ML method, while the Gaussian random signal model results in the stochastic ML method. Although the ML approaches may provide optimal solutions to the antenna array problem, these techniques are often exceedingly complex. For the case of spatially white noise, subspace-based techniques may provide high resolution and/or computationally efficient algorithms. These algorithms may include the Multiple Signal Classification (MUSIC), Estimation of Signal Parameters via Rotational Invariance Techniques (ESPRIT) as well as root-MUSIC for a uniform linear array. The subspace-based techniques also lead to the subspace-based ML approximation approaches, such as efficient Method Of Direction Estimation (MODE). Some of these high resolution AOA estimation techniques, and/or other beamforming techniques, may be used in embodiments of the invention.
In general, the receive (uplink) beamforming and transmit (downlink) beamforming may be quite different for a BTS, depending on the type of system under consideration. Much of beamforming work has been done in the uplink because the knowledge of the wireless channel condition can be easily and robustly estimated. A spatial processing algorithm for a transmitter (such as, for example, a base station transmitter) may be more challenging, since downlink transmission occurs without explicit knowledge of the wireless channel condition(s).
For a Time Division Duplex (TDD) system, as long as the time duplex distance (that is, the temporal separation between transmit and receive time slots) is smaller than the channel coherence time, the uplink channel state typically provides a valid estimate of the downlink channel state, and knowledge of the uplink channel conditions may be used to substantially optimize the downlink transmission. In a Frequency Division Duplex (FDD) system, because the frequency duplex distance (that is, the frequency separation between the uplink channel and the downlink channel) is usually much larger than a channel coherence bandwidth, uplink information may not be useable for downlink beamforming. Moreover, the separation between uplink and downlink carrier frequencies may lead to different antenna spacing in terms of wavelengths for both links even though the same physical antenna array may be used. Consequently, steering vectors in the uplink and downlink may not be the same. Open-loop or closed-loop calibration may be required to derive downlink parameters from measurements of uplink vector channel information in such cases.
According to some embodiments of the invention, systems and methods are provided that may reduce or minimize inter-system interference from a terrestrial BTS (which may or may not be part of an ancillary terrestrial network of a terrestrial/satellite communications system) to an independent communications system (such as Inmarsat satellites) that may be operating in a frequency band that is the same as or proximate to the frequency band used by the terrestrial BTS. Techniques according to some embodiments of the invention may also reduce and/or minimize intra-system interference. A terrestrial BTS may include an adaptive antenna array configured to receive vectors of BTS end user signals and ICS signals, respectively. A vector of signals is a group of related signals, such as signals received over multiple elements of a phased-array antenna. For example, assuming a M-length vector rM,k denotes the received signal vector for the kth signal by an antenna array including M feed elements, the total of K signals received over the M feed element paths may be collectively denoted as vector




  
    r
    M
  
  =
  
    
      ∑
      
        k
        =
        1
      
      K
    
    ⁢
    
      
        r
        
          M
          ,
          k
        
      
      .
    
  




For purposes of this discussion, BTS end user terminals will be referred to herein as “ICS user terminals” or “TDD user terminals.” For purposes of discussion and analysis, the independent communication system may be assumed to be the Inmarsat satellite communications system. However, it will be recognized that embodiments of the invention may be employed to reduce and/or minimize inter-system interference with any separate communications system.

Estimates of array covariance matrices may be derived from the received vectors of a TDD signal and ICS signals, respectively. The frequency transformed ICS covariance matrix and/or spatial signature/AOA, along with the TDD covariance matrix and/or spatial signature/AOA, may be exploited to form a composite covariance matrix and/or a spatial signature matrix that contains spatial information associated with both ICS user terminals and ICS transceiver(s). With an a priori known preamble/pilot for each of the TCS users and the composite TDD/ICS covariance and/or spatial signature matrix obtained from uplink processing of received signals, an adaptive beam-former may determine optimum transmitting weights to form a beam that substantially reduces and/or minimizes inter-system interference by forming “spatial nulls” toward the ICS transceivers. By using the composite covariance matrix, the adaptive beamformer introduces spatial information of ICS transceivers as an additional constraint on top of the desired TDD spatial constraint when forming a beam on a downlink.
Still referring to FIG. 1, in situations where the satellite 210 is associated with a first system operator and the BTS 140 and TCS terminals 120a, 120b are associated with an unrelated second system operator, the frequencies used by the ICS transceivers to receive and/or transmit communications information may be substantially different than the frequencies used by the BTS 140 and/or the TCS terminals 120a, 120b. As such, receivers associated with the BTS 140, operatively configured to receive and process communications signals transmitted by the TCS terminals 120a, 120b, may not substantially detect transmissions of the ICS transceivers to thereby reduce or minimize a radiation pattern of the antenna system of the BTS 140 (smart antenna) in a direction associated with the ICS transceivers to thereby reduce an overload and/or intermodulation interference at the ICS transceivers. Furthermore, even if a receiver of the BTS 140 is configured to detect a communications transmission of the ICS transceivers, an ICS transceiver, prior to deciding to transmit a communications signal, may experience overload and/or intermodulation interference which may prevent the ICS transceiver from transmitting conventional control channel and/or traffic channel signals. Accordingly, in some embodiments of the invention, an ICS transceiver, such as a satellite 210 and/or an ICS user terminal 220, is configured to transmit a signal (e.g. a “beacon”) independent of any overload and/or intermodulation interference condition at the ICS transceiver. The beacon signal may be used by the BTS 140 to obtain angle of arrival (“AOA”) and/or channel estimates for the ICS transceiver.
For example, according to some embodiments of the invention, an ICS transceiver such as an ICS user terminal 220 is configured to transmit a beacon signal as soon as the ICS user terminal 220 is activated (turned-on) by a user. In other embodiments, a special switch/button may be provided at the ICS user terminal 220 which, when activated by a user, initiates a transmission of a beacon signal by the ICS user terminal 220. The beacon signal may be a substantially periodic and/or non-periodic signal and it may be radiated by the ICS user terminal 220 at a relatively low power level (e.g., at about 1 watt or less). The bandwidth of the beacon signal may be as narrow as a single sinusoid or wider than a single sinusoid. The beacon signal may represent an un-modulated carrier or a modulated carrier and the frequency of the carrier may be in any frequency band, such as, for example, an L-band, S-band, ISM band, or any other band. The signal may be radiated by the ICS user terminal 220 using a substantially omni-directional antenna and/or any other antenna configuration. The beacon signal need not correspond to a conventional control channel or traffic channel signal. In some embodiments, the beacon signal is radiated by a device that is proximate to the ICS user terminal 220 but may not be attached to the ICS user terminal 220. In some embodiments, the device may be attached or carried by the user of the ICS user terminal 220 and may be activated by the user of the ICS user terminal 220 to transmit the signal when the user desires to communicate, using the ICS user terminal 220, from a location that is served by a TCS and/or that may be susceptible to interference caused by a TCS.
Referring again to FIG. 1, the dotted arrow 225 that is pointed towards the BTS 140 from the ICS user terminal 220 represents the beacon signal that is radiated by the ICS user terminal 220 (or by the device that is proximate to the ICS user terminal 220). Similarly, the dotted arrow 215 that is pointed towards the BTS 140 from the ICS satellite 210 represents the beacon signal that is radiated by the ICS satellite 210. The BTS 140 is configured to detect a measure of the beacon signal(s) 215, 225 and use it as an error quantity for each BTS 140 receiver in forming respective antenna patterns associated with the plurality of active TCS terminals 120A, 120B, thereby allowing each antenna pattern to generate a low gain characteristic (a “null”) in a direction associated with the received signal from an ICS transceiver. Since, in accordance with a TDD air interface protocol, a transmit antenna pattern is substantially the same as a corresponding receive antenna pattern (at least over a TDD interval/frame), having configured the BTS 140 to form receive antenna patterns with substantially low gain (nulls) in directions associated with the ICS transceivers, the corresponding transmit antenna patterns of the BTS 140 may also provide low gain (nulls) in the directions associated with the ICS transceivers, thereby reducing an overload and/or intermodulation interference at the ICS transceivers. In some embodiments, the BTS 140 is configured to detect a measure of the beacon signal(s) transmitted by the ICS transceivers and use it as an error quantity for at least one BTS 140 receiver in forming at least one antenna pattern associated with at least one active TCS terminal thereby allowing the at least one antenna pattern to generate a low gain characteristic (a “null”) in a direction associated with the received beacon signal(s) 215, 225 from the ICS transceivers.
It will be understood that some embodiments of the present invention that are described herein include ICS transceivers that transmit a beacon signal and allow the BTS 140 to be configured to receive this signal. However, in other embodiments, the BTS 140 may be configured to receive a standard control or traffic channel signal from the ICS transceiver(s) and/or may be configured to determine a location of an ICS transceiver independent of a signal that is transmitted from the ICS transceiver. For example, GPS and/or other techniques may be used by the BTS 140 and/or another infrastructure element, to provide a measure of the location of the ICS transceivers to the BTS 140.
It also will be understood that embodiments of the invention are described herein relative to TDD air interface protocols because in a TDD air interface protocol, a transmit antenna pattern generally is substantially the same as a corresponding receive antenna pattern (at least over a TDD interval/frame). However, embodiments of the invention may also be employed with non-TDD air interfaces, such as frequency division duplex (FDD) air interfaces, wherein the differences between the propagation characteristics and/or other system/component characteristics associated with the differing transmit and receive frequencies and/or with different components associated therewith may be ignored or compensated for, for example by using modeling techniques and/or other measures of these differences.
Furthermore, it will be understood by those having skill in the art that embodiments of the present invention describe the use of “smart” antennas and/or adaptive beam forming to reduce or minimize overload and/or intermodulation interference to an ICS transceiver that is communicating in a TCS service area. In other embodiments of the present invention, however, an BTS 140 that is configured to communicate with one or more TCS terminals by using and/or reusing satellite frequencies, may be provided with a smart antenna to increase or maximize a radiation that is directed to the TCS terminal(s), while minimizing or reducing transmissions that are not directed to the TCS terminal(s) and which may be wasted and/or cause interference to other devices (including other TCS terminals, cellular and/or PCS terminals). Smart antennas also may be used to increase or maximize a gain for signals that are received from TCS terminal(s) while reducing a BTS 140 receiver gain for signals that are not associated with TCS terminals and/or are associated with co-frequency/co-channel TCS terminals.
Finally, it will be understood that a smart antenna or an adaptive antenna may include a linear (one-dimensional) array of antenna elements or a multi-dimensional (e.g., two-dimensional) array of antenna elements and that a smart antenna may, in some embodiments, be configured to radiate signals that are polarized substantially in one polarization sense such as, for example, Left Hand Circularly Polarized (LHCP), Right Hand Circularly Polarized (RHCP), Linearly Polarized (LP) and may be configured to receive signals over at least two polarizations that are substantially in spatial quadrature with each other.
1. System Configuration and Model
FIG. 2A is a schematic illustration of a BTS 140 serving a cell 130 of a terrestrial communications system (TCS) 100 which may or may not be an ancillary terrestrial network of a terrestrial/satellite communications system. The cell 130 is divided into three sectors 130A, 130B, 130C. However, the cell 130 may be divided into more or less than three sectors. Each sector is served by a pair of or more antennas 180. The BTS 140 may communicate on a downlink signal 118 with a TCS user terminal 120 located in a sector 130C of the cell 130. The downlink signal 118 is drawn in FIG. 2A to illustrate a main lobe 118A of the signal 118 directed at the TCS user terminal 120 and sidelobes 118B directed in other directions. Thus, as shown in FIG. 2A, while a beam formed to the TCS user terminal 120 may have a maximum gain in the direction of the TCS user terminal 120, the beam may have sidelobes aimed in other directions.
Also illustrated in FIG. 2A are a pair of satellites 210A, 210B from an independent communications system (ICS) 200, which may be Inmarsat satellites, for example. The ICS satellites 210A, 210B communicate with an ICS user terminal 220. While only two ICS satellites 210A, 210B and a single ICS user terminal 220 are illustrated in FIG. 2A, it will be understood that more or less than two ICS satellites 210A, 210B may be visible from a particular terrestrial location at a given time and/or that more or less than one ICS user terminal 220 may be in or near the cell 130 served by the BTS 140. For purposes of the present discussion, the ICS satellites 210A, 21B and/or the ICS user terminal 220 may be referred to collectively as “ICS transceivers” or individually as an “ICS transceiver.”
The ICS satellites 210A, 2101B, which may operate in FDD mode, may be on or above the horizontal plane of the antennas 180, and may be receiving at carrier frequency frx and transmitting at carrier frequency ftx, as shown in FIG. 2A. Likewise, the ICS user terminal 220 may be receiving signals from one or more of the ICS satellites at carrier frequency ftx and may be transmitting signals at carrier frequency frx. The terrestrial cellular system (TCS) 100 including the BTS 140 is assumed to be a time division duplex (TDD) system that operates at carrier frequency fc for both down-link (DL) and uplink (UL) communications. In general, the carrier frequency fc used by the BTS 140 may not be equal to the transmitting carrier frequency ftx for the ICS system 200. However, the carrier frequency fc may be equal or approximately equal to the receiving carrier frequency frx for the ICS 200. Without loss of generality, it may be assumed that frx is equal to fc. In order to reduce or minimize interference to the ICS satellites 210A, 210B and/or the ICS user terminal 220 from transmissions by the BTS 140, the BTS 140 may obtain knowledge of where the ICS satellites 210A, 210B and the ICS user terminal 220 are located. The BTS 140 may be configured with an antenna array 170 and a detector 175 that is configured to detect the transmit carrier frequency ftx of the ICS satellites 210A, 210B and the transmit carrier frequency frx (or fc) and, responsive to the detection, to estimate the AOAs and/or spatial signatures of the ICS satellites 210A, 210B and/or the ICS user terminal 220 relative to the orientation of the BTS antenna array 170, even though the BTS 140 may not have any knowledge about the ICS signal other than its transmitting carrier frequency and/or its beacon signal.
Referring to FIGS. 2A and 2B, the antenna array 170 at the BTS 140 may include 3 sectors of linear arrays. The linear array in each sector may include a number of antennas 180, each of which is a Vertical Uniform Linear Array (VULA) of radiating elements 182. However, other structures, including non-sector circular array with a greater or lesser number of VULAs 180 are also possible. In addition, an antenna in the antenna array 170 may include a phased-array antenna including a two dimensional array of antenna elements capable of forming a beam having both azimuth and elevation components.
FIG. 2A illustrates a BTS having 3 sectors of VULAs 180. Each VULA 180 may include M antenna elements 182 that have adjacent element spacing d. The gain and phase of the antenna elements 182 (with the first element as a reference), are given by a complex-valued steering vector a(θ) as a function of AOA in elevation. For a VULA 180 tuned to the transmit carrier frequency ftx, a steering vector may be given by:



  
    
      
        a
        ICS
        UL
      
      ⁡
      
        (
        θ
        )
      
    
    =
    
      
        [
        
          1
          ⁢
          
            
          
          ⁢
          
            e
            jΩ
          
          ⁢
          
            
          
          ⁢
          ...
          ⁢
          
            
          
          ⁢
          
            e
            
              
                j
                ⁡
                
                  (
                  
                    M
                    -
                    1
                  
                  )
                
              
              ⁢
              Ω
            
          
        
        ]
      
      T
    
  




  
    
      where
      ⁢
      
        
      
      ⁢
      Ω
    
    =
    
      
        
          2
          ⁢
          π
          ⁢
          
            
          
          ⁢
          
            f
            tx
          
        
        c
      
      ⁢
      d
      ⁢
      
        
      
      ⁢
      
        sin
        ⁡
        
          (
          θ
          )
        
      
      ⁢
      
        
      
      ⁢
      is
      ⁢
      
        
      
      ⁢
      a
      ⁢
      
        
      
      ⁢
      spatial
      ⁢
      
        
      
      ⁢
      
        frequency
        .
      
    
  



If there are K ICS transceivers each transmitting a signal, the array output vector can be expressed as a linear combination of the K incident waveforms and noise as below:




  
    
      
        
          
            y
            ICS
          
          ⁡
          
            (
            t
            )
          
        
        =
        
          
            
              ∑
              
                k
                =
                0
              
              
                K
                -
                1
              
            
            ⁢
            
              
                
                  a
                  ICS
                  UL
                
                ⁡
                
                  (
                  
                    θ
                    k
                  
                  )
                
              
              ⁢
              
                
                  s
                  k
                
                ⁡
                
                  (
                  t
                  )
                
              
            
          
          +
          
            n
            ⁡
            
              (
              t
              )
            
          
        
      
    
  
  
    
      
        =
        
          
            
              
                A
                ICS
              
              ⁢
              
                s
                ⁡
                
                  (
                  t
                  )
                
              
            
            +
            
              n
              ⁡
              
                (
                t
                )
              
            
          
          ∈
          
            C
            
              M
              ⁢
              
                
              
              ⁢
              x
              ⁢
              
                
              
              ⁢
              1
            
          
        
      
    
  




where

AICS=[aICSUL(θ1)aICSUL(θ2) . . . aICSUL(θK)]

is the array manifold, and

s(t)=[s1(t)s2(t) . . . sK(t)]T 

is a vector of signal waveforms modeled as deterministic unknowns.

For an array tuned to L TCS user terminals 120 in the terrestrial system, a received array vector may be given by:



  
    
      
        
          
            
              y
              TDD
            
            ⁡
            
              (
              t
              )
            
          
          =
          
            
              
                ∑
                
                  p
                  =
                  0
                
                
                  L
                  -
                  1
                
              
              ⁢
              
                
                  
                    a
                    TDD
                    UL
                  
                  ⁡
                  
                    (
                    
                      θ
                      p
                    
                    )
                  
                
                ⁢
                
                  
                    
                      s
                      _
                    
                    p
                  
                  ⁡
                  
                    (
                    t
                    )
                  
                
              
            
            +
            
              
                n
                _
              
              ⁡
              
                (
                t
                )
              
            
          
        
      
    
    
      
        
          =
          
            
              
                
                  A
                  TDD
                
                ⁢
                
                  
                    s
                    _
                  
                  ⁡
                  
                    (
                    t
                    )
                  
                
              
              +
              
                
                  n
                  _
                
                ⁡
                
                  (
                  t
                  )
                
              
            
            ∈
            
              C
              
                M
                ⁢
                
                  
                
                ⁢
                x
                ⁢
                
                  
                
                ⁢
                1
              
            
          
        
      
    
  




  
    
      where
      ⁢
      
        
      
      ⁢
      
        A
        TDD
      
    
    =
    
      
        [
        
          
            
              a
              TDD
              UL
            
            ⁡
            
              (
              
                θ
                1
              
              )
            
          
          ⁢
          
            
          
          ⁢
          
            
              a
              TDD
              UL
            
            ⁡
            
              (
              
                θ
                2
              
              )
            
          
          ⁢
          
            
          
          ⁢
          ⋯
          ⁢
          
            
          
          ⁢
          
            
              a
              TDD
              UL
            
            ⁡
            
              (
              
                θ
                L
              
              )
            
          
        
        ]
      
      ⁢
      
        
      
      ⁢
      and
    
  




  
    
      
        s
        _
      
      ⁡
      
        (
        t
        )
      
    
    =
    
      
        
          [
          
            
              
                
                  s
                  _
                
                1
              
              ⁡
              
                (
                t
                )
              
            
            ⁢
            
              
            
            ⁢
            
              
                
                  s
                  _
                
                2
              
              ⁡
              
                (
                t
                )
              
            
            ⁢
            
              
            
            ⁢
            ⋯
            ⁢
            
              
            
            ⁢
            
              
                
                  s
                  _
                
                L
              
              ⁡
              
                (
                t
                )
              
            
          
          ]
        
        T
      
      ⁢
      
        
      
      ⁢
      with
      ⁢
      
        :
      
    
  




  
    
      
        a
        TDD
        UL
      
      ⁡
      
        (
        θ
        )
      
    
    =
    
      
        [
        
          1
          ⁢
          
            
          
          ⁢
          
            e
            
              j
              ⁢
              
                
                  2
                  ⁢
                  π
                  ⁢
                  
                    
                  
                  ⁢
                  
                    f
                    c
                  
                
                c
              
              ⁢
              d
              ⁢
              
                
              
              ⁢
              
                sin
                ⁡
                
                  (
                  θ
                  )
                
              
            
          
          ⁢
          
            
          
          ⁢
          ⋯
          ⁢
          
            
          
          ⁢
          
            e
            
              
                j
                ⁡
                
                  (
                  
                    M
                    -
                    1
                  
                  )
                
              
              ⁢
              
                
                  2
                  ⁢
                  π
                  ⁢
                  
                    
                  
                  ⁢
                  
                    f
                    c
                  
                
                c
              
              ⁢
              d
              ⁢
              
                
              
              ⁢
              
                sin
                ⁡
                
                  (
                  θ
                  )
                
              
            
          
        
        ]
      
      T
    
  



[Note that for purposes of the present discussion, the subscript “TDD” refers to the TCS user terminals, since the exemplary embodiments involve a BTS operating in TDD mode. However, other communications protocols could be employed.]

In order to perform the AAS operations, the BTS 140 may be configured to:
a) estimate the spatial signature and/or AOA of the ICS satellites 210A, 210B and the ICS user terminal(s) 220;
b) calibrate an antenna array vector channel response difference between receiving and transmitting due to the use of different carrier frequencies for the uplink and downlink by the ICS satellites 210A, 210B and the ICS user terminal(s) 220;
c) form a composite array covariance matrix or composite spatial signature matrix; and
d) form a downlink beam based on the composite covariance matrix or spatial signature matrix that includes constraint information from both the TCS user terminals and the ICS transceivers.
2. Description of Exemplary Embodiments
In some embodiments of the invention, a BTS 140 uses an antenna array 170 and detector 175 tuned to the ICS transmitting carrier frequencies ftx to estimate the spatial signatures for the ICS satellites 210A, 210B and/or the tuned to the carrier frequencies frx to estimate the spatial signatures of the ICS user terminals 220. The spatial signature estimates are used as partial inputs for adaptive downlink (DL) beamforming by the BTS 140. The adaptive DL beam-forming uses the spatial signature estimates for the ICS transceivers along with the information obtained from the adaptive uplink (UL) beamforming for a desired TDD terminal to form a beam that has a peak substantially pointed to the desired mobile terminal and a “spatial null” against the ICS satellite(s) and/or other TCS user terminals. Adaptive UL beam-forming is described in the literature and need not be discussed further herein.
Since the ICS transceivers belong to an independent communications system, the BTS 140 may estimate the spatial signatures of the ICS transceivers by blind spatial channel estimation techniques, because the ICS signals may not be known by the BTS 140. To use the received ICS vector channel estimates of the ICS downlink channel (which has a carrier frequency at ftx) or the ICS uplink channel (which has a carrier frequency at frx) in the BTS 140 downlink transmit beamforming process, a calibration/transformation from uplink to downlink may be performed, because the ICS downlink may operate at a different carrier frequency from the uplink. The frequency duplex distance between the uplink and the downlink signals used by the ICS satellites 210A, 210B may result in different steering vectors between uplink and downlink. A linear transformation through rotation of subspaces that provides an estimate of the downlink steering vectors from the uplink vectors may be used. A downlink calibration table of an antenna array is provided as a matrix that contains the downlink steering vector for a set of azimuth directions {θ1, θ2, . . . , θL}

ADL=[aICSDL(θ1), aICSDL(θ2), . . . , aICSDL(θL)]

where





  
    
      a
      ICS
      DL
    
    ⁡
    
      (
      θ
      )
    
  
  =
  
    
      [
      
        1
        ⁢
        
          
        
        ⁢
        
          e
          
            j
            ⁢
            
              
                2
                ⁢
                π
                ⁢
                
                  
                
                ⁢
                
                  f
                  c
                
              
              c
            
            ⁢
            d
            ⁢
            
              
            
            ⁢
            
              sin
              ⁡
              
                (
                θ
                )
              
            
          
        
        ⁢
        
          
        
        ⁢
        ⋯
        ⁢
        
          
        
        ⁢
        
          e
          
            
              j
              ⁡
              
                (
                
                  M
                  -
                  1
                
                )
              
            
            ⁢
            
              
                2
                ⁢
                π
                ⁢
                
                  
                
                ⁢
                
                  f
                  c
                
              
              c
            
            ⁢
            d
            ⁢
            
              
            
            ⁢
            
              sin
              ⁡
              
                (
                θ
                )
              
            
          
        
      
      ]
    
    T
  




Similarly for the uplink calibration table:

AUL=[aICSUL(θ1), aICSUL(θ2), . . . , aICSUL(θL)]

A transformation problem is to find the best transformation matrix T that minimizes the difference between the array of downlink steering vectors and the transformed array of uplink steering vectors, namely, a transformation matrix T that satisfies the following criterion:

min∥ADL−TAUL∥2 subject to THT=IM 

With the known calibration table matrices for both downlink and uplink, it can be shown that the optimal orthogonal matrix TεCM×M may be given by:

T=UVH 

where U and V are the left and right singular vectors of ADL(AUL)H by calculating the Singular Value Decomposition (SVD) of ADL(AUL)H, i.e.,

ADL(AUL)H=UΣVH 

Applying the transformation matrix T to the estimate of a spatial steering vector ĉICSUL(θk), from the uplink, we obtain the steering vector estimate of the satellites for the downlink as follows:

ĉICSDL(θk)=TĉICSUL(θk), k=1, . . . , K 

The steering vector estimates for the ICS user terminals 220 may or may not need to be transformed.

For a TDD-based TCS 100, the same or a similar transformation technique may be used to compensate a mismatch between multiple receiver chains and multiple transmit chains. An UL calibration table and a corresponding DL calibration may be obtained from the measurements of receive and transmit signals through a calibration unit.
FIG. 3 illustrates, in block diagram form, a base station architecture and associated methods for a BTS 140 according to some embodiments of the invention. A TDD BTS 140 with an adaptive antenna system may adaptively form a DL beam that may increase or maximize the signal-to-interference plus noise ratio for the desired TCS user while forming a “null” toward the ICS transceiver by using the UL received vectors. In the uplink, the TDD BTS 140 may receive both the TCS signals and ICS signals though a single antenna array. Signals received at the antenna array pass through a transmit/receive duplexer (block 310) to an RF receive front end processor (block 315). The received signals are down-converted and sampled by an A/D and downconverter (block 320). The sampled digital signals are then filtered by both a TCS channel digital filter (block 325) along a first processing path 322 and an ICS channel digital filter (block 310) along a second processing path 324.
The TCS signal path 322 and the ICS signal path 324 form respective array vector outputs. From the TCS array vector output and the ICS array vector output, respectively, the TCS array covariance matrix (block 330) and ICS array covariance matrix (block 345) are estimated. The TCS covariance matrix is input to an adaptive UL beamformer (block 335), which generates a beamformed baseband received TCS signal.
Depending on the form of the embodiment, the spatial signatures may be estimated. For the ICS path 324, the covariance matrix may be processed for frequency transformation/calibration (block 350) before it is used with the counterpart from the TCS path 322 for the DL adaptive beamforming by an adaptive DL beamformer (block 355). The weights output by the adaptive beamformer (block 355) are compensated by a weights compensator (block 365) based on the covariance matrix (block 330) that has been processed by a calibration unit (block 360). The calibration unit (block 360) estimates any differential gain and phase introduced by receiver and transmitter hardware mismatch among the multi-antenna paths. The differential gain and phase shifts among the multi-antenna paths are compensated through the weight compensator (block 365).
The compensated weights are applied to the baseband transmit signal (block 370), and the resulting signal is converted to analog and upconverted (block 375). The analog signal is then amplified by a front end transmit processor (block 380) and transmitted through the duplexer (block 310).
FIG. 4 shows embodiments of an adaptive downlink beamformer using MMSE and/or LCMV based algorithms. Uplink adaptive beamforming may be accomplished optimally with MMSE and/or LCMV based algorithms without any concern about the ICS signals, since the ICS downlink transmission does not cause interference to the base station. However, if the weights obtained from the uplink beamforming were directly applied to downlink beamforming, there would generally be no guarantee of providing a null in a direction associated with an ICS satellite in the downlink. To create a null in the direction of an ICS satellite in the downlink beamforming, we can create a pseudo uplink beamforming that is constrained by both TCS user terminals and ICS transceivers at the TDD carrier frequency fc, and then apply the pseudo uplink beamforming weights for the downlink beamforming. Specifically, the following procedures, which are referred to as Algorithm I, may be performed:

    Step 1: Estimate the uplink covariance matrix for ICS transceivers {circumflex over (R)}ICS that is given by:






  
    
      
        R
        ^
      
      ICS
    
    ⁡
    
      (
      
        
          f
          tx
        
        ,
        θ
      
      )
    
  
  =
  
    
      1
      N
    
    ⁢
    
      
        ∑
        
          t
          =
          1
        
        N
      
      ⁢
      
        
          
            y
            ICS
          
          ⁡
          
            (
            t
            )
          
        
        ⁢
        
          
            y
            ICS
            H
          
          ⁡
          
            (
            t
            )
          
        
      
    
  




    Step 2: Perform frequency transformation for the ICS uplink covariance matrix from the ICS carrier frequency ftx to the TCS carrier frequency fc:

{circumflex over (R)}ICS(fc,θ)=[Tâ(ftx,θ)][Tâ(ftx,θ)]H=T{circumflex over (R)}ICS(ftx,θ)TH, k=1, . . . , K 

    Step 3: Derive the composite uplink covariance matrix by adding the frequency calibrated ICS uplink covariance matrix to the TCS uplink covariance matrix to get the composite uplink covariance matrix, which yields:





  
    
      
        R
        ^
      
      ⁡
      
        (
        
          
            f
            c
          
          ,
          θ
        
        )
      
    
    =
    
      
        
          
            R
            ^
          
          TDD
        
        ⁡
        
          (
          
            
              f
              c
            
            ,
            θ
          
          )
        
      
      +
      
        
          
            R
            ^
          
          ICS
        
        ⁡
        
          (
          
            
              f
              c
            
            ,
            θ
          
          )
        
      
    
  




  
    
      where
      ⁢
      
        
      
      ⁢
      
        
          
            R
            ^
          
          TDD
        
        ⁡
        
          (
          
            
              f
              c
            
            ,
            θ
          
          )
        
      
    
    =
    
      
        1
        N
      
      ⁢
      
        
          ∑
          
            t
            =
            1
          
          N
        
        ⁢
        
          
            
              y
              TDD
            
            ⁡
            
              (
              t
              )
            
          
          ⁢
          
            
              y
              TDD
              H
            
            ⁡
            
              (
              t
              )
            
          
        
      
    
  



    Step 4: Obtain the estimate of the desired user's spatial signature using preamble or pilot reference rk as:






  
    
      a
      ^
    
    ⁡
    
      (
      
        
          f
          c
        
        ,
        
          θ
          k
        
      
      )
    
  
  =
  
    
      E
      ⁡
      
        [
        
          y
          ⁢
          
            
          
          ⁢
          
            r
            k
            *
          
        
        ]
      
    
    =
    
      
        1
        N
      
      ⁢
      
        
          ∑
          
            t
            =
            1
          
          N
        
        ⁢
        
          
            y
            ⁡
            
              (
              t
              )
            
          
          ⁢
          
            
              r
              k
              *
            
            ⁡
            
              (
              t
              )
            
          
        
      
    
  




    Step 5: Determine the weights of pseudo uplink beamforming by:






  
    
      w
      ^
    
    k
  
  =
  
    
      
        
          
            R
            ^
          
          
            -
            1
          
        
        ⁡
        
          (
          
            
              f
              c
            
            ,
            θ
          
          )
        
      
      ⁢
      
        
          a
          ^
        
        ⁡
        
          (
          
            
              f
              c
            
            ,
            
              θ
              k
            
          
          )
        
      
      ⁢
      
        
      
      ⁢
      or
      ⁢
      
        
      
      ⁢
      
        
          w
          ^
        
        k
      
    
    =
    
      
        
          
            
              R
              ^
            
            
              -
              1
            
          
          ⁡
          
            (
            
              
                f
                c
              
              ,
              θ
            
            )
          
        
        ⁢
        
          
            a
            ^
          
          ⁡
          
            (
            
              
                f
                c
              
              ,
              
                θ
                k
              
            
            )
          
        
      
      
        
          
            
              a
              ^
            
            H
          
          ⁡
          
            (
            
              
                f
                c
              
              ,
              
                θ
                k
              
            
            )
          
        
        ⁢
        
          
            
              R
              ^
            
            
              -
              1
            
          
          ⁡
          
            (
            
              
                f
                c
              
              ,
              θ
            
            )
          
        
        ⁢
        
          
            a
            ^
          
          ⁡
          
            (
            
              
                f
                c
              
              ,
              
                θ
                k
              
            
            )
          
        
      
    
  




    Step 6: Apply the weights of pseudo uplink beamforming obtained in Step 4 to the downlink transmit beamforming:

xk=ŵk*sk 



Since the adaptively formed beam is constrained by both intra-system and inter-system constraints, the adaptively formed beam may provide an increased gain toward the desired user while minimizing or reducing the inter-system interference and intra-system interference to the extent allowed by the number of available degrees of freedom.
FIG. 5 illustrates embodiments of an adaptive downlink beam-former using an algorithm according to further embodiments of the invention, referred to herein as Algorithm II. The array covariance matrix is defined by:

R=E{y(t)yH(t)}=ARsAH+σ2I 

where Rs=E{s(t)sH(t)} denotes the source covariance matrix, and σ2 is the thermal noise variance. The eigenvectors of the covariance matrix R belong to either of the two orthogonal subspaces—the principal eigen subspace called the signal subspace and the non-principal eigen subspace called the noise subspace. The steering vectors related to AOAs lie in the signal subspace and hence are orthogonal to the noise subspace. Therefore, it is possible to estimate signal steering vectors by performing eigen-decomposition of the array covariance matrix. The estimate of the covariance matrix {circumflex over (R)} is given by:





  
    R
    ^
  
  =
  
    
      
        1
        N
      
      ⁢
      
        
          ∑
          
            t
            =
            1
          
          N
        
        ⁢
        
          
            y
            ⁡
            
              (
              t
              )
            
          
          ⁢
          
            
              y
              H
            
            ⁡
            
              (
              t
              )
            
          
        
      
    
    =
    
      
        
          
            E
            ^
          
          s
        
        ⁢
        
          
            Λ
            ^
          
          s
        
        ⁢
        
          
            E
            ^
          
          s
          H
        
      
      +
      
        
          
            E
            ^
          
          n
        
        ⁢
        
          
            Λ
            ^
          
          n
        
        ⁢
        
          
            E
            ^
          
          n
          H
        
      
    
  




The eigen-decomposition of {circumflex over (R)} results in the representation:




  
    
      R
      ^
    
    =
    
      
        
          
            
              E
              ^
            
            s
          
          ⁢
          
            
              Λ
              ^
            
            s
          
          ⁢
          
            
              E
              ^
            
            s
            H
          
        
        +
        
          
            
              E
              ^
            
            n
          
          ⁢
          
            
              Λ
              ^
            
            n
          
          ⁢
          
            
              E
              ^
            
            n
            H
          
        
      
      =
      
        
          ∑
          
            i
            =
            1
          
          M
        
        ⁢
        
          
            
              λ
              ^
            
            i
          
          ⁢
          
            
              e
              ^
            
            i
          
          ⁢
          
            
              e
              ^
            
            i
            H
          
        
      
    
  
  ,




where {circumflex over (λ)}1, > . . . >{circumflex over (λ)}K>{circumflex over (λ)}K+1= . . . ={circumflex over (λ)}M=σ2. The matrix Ês=[ê1, . . . , êK] contains the K eigenvectors corresponding to the largest eigenvalues. The range space of Ês is the signal subspace. Its orthogonal complement is the noise subspace and is spanned by the columns of Ên=[êK+1, . . . , êM].

For K ICS satellites with covariance matrix {circumflex over (R)}ICS, the estimated spatial matched filter (eigen-filter) bank is given by estimates of signal subspace, i.e:

[ĉICSUL(θ1), . . . , ĉICSUL(θK)]=Ês=[ê1, . . . , êK]εCM×K 

By using the transformation matrix T in the estimation of the spatial steering vector ĉICSUL(θk), from the uplink, the downlink steering vector estimates of the ICS satellites can be obtained by ĉICSDL(θk)=TĉICSUL(θk), k=1, . . . , K. Alternatively, a noise subspace or signal subspace based high-resolution AOA estimation algorithms, such as MUSIC, and MODE, may be used to estimate θk. With the estimation of θk, the downlink steering vector estimates of the ICS satellites can be obtained by





  
    
      
        
          c
          ^
        
        ICS
        DL
      
      ⁡
      
        (
        
          θ
          k
        
        )
      
    
    =
    
      
        [
        
          
            
              1
            
            
              
                e
                
                  j
                  ⁢
                  
                    
                      2
                      ⁢
                      π
                      ⁢
                      
                        
                      
                      ⁢
                      
                        f
                        c
                      
                    
                    c
                  
                  ⁢
                  d
                  ⁢
                  
                    
                  
                  ⁢
                  
                    sin
                    ⁡
                    
                      (
                      
                        θ
                        k
                      
                      )
                    
                  
                
              
            
            
              ⋯
            
            
              
                e
                
                  
                    j
                    ⁡
                    
                      (
                      
                        M
                        -
                        1
                      
                      )
                    
                  
                  ⁢
                  
                    
                      2
                      ⁢
                      π
                      ⁢
                      
                        
                      
                      ⁢
                      
                        f
                        c
                      
                    
                    c
                  
                  ⁢
                  d
                  ⁢
                  
                    
                  
                  ⁢
                  
                    sin
                    ⁡
                    
                      (
                      
                        θ
                        k
                      
                      )
                    
                  
                
              
            
          
        
        ]
      
      T
    
  
  ,
  
    

  
  ⁢
  
    k
    =
    1
  
  ,
  ⋯
  ⁢
  
    
  
  ,
  K




The downlink beamforming may be based on the following input information: 1) estimate of the steering vector for the ICS satellites at downlink ĉICSDL(θk); 2) estimate of spatial signatures (steering vectors) for the desired TCS user ĥDL(φ1).
Given that these estimates are available for the downlink beamforming, in the case of single TCS user and K ICS satellites, to form a beam that has unity response in the direction of the desired TCS user, and null response in the directions of the ICS satellites, the desired weight vector is subject to the following constraints:

wHĥDL(φ1)=1

wHĉICSDL(θk)=g, k=1, . . . , K 

where g is the constraint beam gain toward the ICS transceivers.

By forming a matrix B whose columns are the steering vectors associated with directions of desired TCS user and ICS transceivers:

B≡[ĥDL(φ1), ĉICSDL(θ1), . . . , ĉICSDL(θK)]εCM×(K+1) 

We can rewrite the above constraint equations as:

wHB=i1T 

where i1=[1, g, . . . , g]T.

To provide an antenna array that has enough degrees of freedom to form K nulls, one may require K≦M−1. A general form of solution for the estimate of the weighting vector may be given by:

wH=i1T(BHB)−1BH εC1×M 

where the term (BHB)−1BH is the pseudo inverse of matrix B.

In summary, downlink beam-forming based on Algorithm II may include the following steps:

    Step 1: Estimate the uplink covariance matrix {circumflex over (R)}ICS for the ICS transceivers.
    Step 2: Perform eigen-decomposition of {circumflex over (R)}ICS.
    Step 3: Obtain the estimates of the spatial signatures for the ICS transceivers from the principal eigenvector or the K eigenvectors corresponding to K largest eigenvalues if K ICS satellites are involved.
    Step 4: Perform frequency transformation/calibration by using the transformation matrix to obtain the ICS spatial signature at the TDD frequency band. Alternatively, estimate the AOAs using MUSIC or MODE, and then construct the ICS spatial signature at the TDD frequency band with the AOA information.
    Step 5: Estimate the spatial signature for the TCS user terminals by performing correlation between the received TDD array vectors and the preamble/pilot sequence if preamble/pilot is available.
    Step 6: Form a matrix B whose columns are filled with the estimated spatial signature vectors associated with directions of the TCS user terminals and the ICS transceivers.
    Step 7: Obtain the pseudo-inverse of matrix B, and obtain weights.
    Step 8: Apply the weighting vector obtained in Step 7 to the TCS user for the downlink transmit beamforming.


3. Performance Simulations
Simulation models consider the TCS operation in TDD mode where the BTS may be configured with an adaptive antenna array to thereby receive vectors of terrestrial TDD end user signals and ICS satellite signals, respectively. The ICS satellite(s) may operate in FDD mode where the downlink carrier frequency ftx may not overlap the TCS carrier frequency fc. However the ICS uplink carrier frequency frx may be co-channeled with the TCS carrier frequency fc. The ICS signals may impinge on the BTS antenna array 170 from low elevation angles at the horizontal plane and above the horizontal plane. The ICS signals are assumed to be modulated using QPSK modulation. The vertical uniform linear antenna array (VULA) 180 includes M antenna elements 182 with a half-wave length spacing between adjacent elements. The simulation models assume the following key system parameters













TCS mode
TDD


BTS carrier frequency fc
1634.5 MHz


ICS downlink frequency ftx
1549.85 MHz


ICS uplink frequency frx
1634.5 MHz


ICS signal modulation
QPSK


Number of antenna elements of VULA
12


Spacing between two adjacent elements
Half-wavelength (9.18 cm)


BTS received ICS SNR
5 dB


AAA beam gain constraint g toward ICS
−20 dB, −28 dB, −34 dB


(for Algorithm II)


Adaptive beamforming algorithms
1) Algorithm II with AOA



estimation



2) Algorithm I









The BTS AAS receiver is 175 configured to operate in dual bands in parallel. One band is for the ATC carrier frequency fc. The other band is for the ICS downlink frequency ftx. To ensure that the BTS 140 receives a reasonable level of ICS signals, the BTS adaptive antenna array 170 is configured to receive the ICS signals at a local ICS downlink frequency ftx. The ICS angle-of-arrival (AOA) estimates are obtained using the MUSIC algorithm with 1000 samples of received vector signals at a time.
In a first example, one ICS satellite is located at an elevation angle of 5 degrees while the desired TCS user terminal is located at elevation angle of −5 degrees. By processing the received vectors of ICS signals and the TDD end user signals as described in the foregoing section, the BTS adaptive antenna array may adaptively form a maximum, or near maximum beam toward the desired TCS user and a beam that has a radiated power toward the ICS satellite that is −20 dB down from the maximum direction. FIG. 6A shows the simulated downlink transmission beam pattern of the BTS in a polar plot, and FIG. 6B shows the simulated downlink transmission beam pattern in a linear plot.
The next example considers two ICS satellites at elevation angle of 5 degrees and 8 degrees and the desired TCS user terminal at elevation angle of −5 degrees. As shown in FIGS. 7A and 7B, the adaptively formed downlink beam pattern has a gain of at least −20 dB toward the satellites and a near maximum gain toward the desired user. For most of the unconstrained region above the horizon, the sidelobes are also well under −20 dB. This is because the constraint gain for the satellites is set to a limited value. The deeper the constraint gain is, the higher the unconstrained sidelobes. Thus there is a trade-off between the depth of nulls and gains of the unconstrained sidelobes.
The third example assumes a scenario where three ICS satellites are located at angles of 7, 17 and 20 degrees respectively with the desired TDD end user at angle of −5 degrees. FIGS. 8A/8B, 9A/9B and 10A/10B show the formed beam patterns when the constraint gain against the satellites is set to −20 dB, −28 dB and −34 dB respectively. The resulting average gains for the satellites are about the same as the constraint values respectively for each of the three cases. It can be seen that while the “nulls” against satellites become deeper as deeper constraints are used, the sidelobes in the unconstrained region become higher as well. However, in this example, the sidelobes are still well below −20 dB.
When the number of visible ICS satellites becomes large (the high resolution AOA estimation algorithm can handle up to 11 satellites simultaneously), the number of constraint points also increases. To satisfy more constraint points, the beam gain toward the desired TCS user may be sacrificed to some degree due to the adaptive nature of the algorithm, depending on various scenarios. Monte-Carlo simulation results show the statistical performance of the adaptive array in terms of the average beam gain toward the desired TCS user and the average beam gain against the ICS satellites. In the simulation, the number of satellites is varied from 1 to 5 with uniform random distributed within angles of 0 to 30 degrees. The results are based on 1000 independent trials with the beam gain constraint for the satellites being set to −20 dB. FIG. 11 shows the average gains of the formed beam toward the desired user and the satellites. It can be seen that the average gain for the desired user has very little degradation when the number of visible satellites is equal to or less than three. A noticeable beam gain reduction for the desired user begins when number of the visible satellites is equal or greater than four. However, the overall beam gain against the satellites also becomes more negative as the number of the visible satellite increases.
The last example considers the case that a TDD user is located at near horizon with elevation angle from −1 degree to −5 degrees and 1 to 3 ICS satellites are at elevation angles with uniformly random distributed in [0, −10] degrees. For this case, where a TDD user is very close to horizon and the angular separation between the TDD user and one or more ICS satellites is likely very small, using Algorithm II may cause severe gain loss for the desired TDD user because the beamformer tries to satisfy the satellite gain constraint which respects to the desired TDD user and consequently has to overly tilt down the beam pointing position. Algorithm I may provide better performance in this case because the algorithm puts the desired TDD user as priority while trying to minimize the interference with no extra specific constraints. FIG. 12 shows the performance results of the Monte Carlo simulation with Algorithm I. It appears that the algorithm provides more robust performance than Algorithm II, especially for the TDD user at a low elevation angle. The algorithm offers less gain degradation for the TDD user while providing decent nulls toward the satellites. The gain degradation is mild for the TDD user at −1 degree, and approaches zero dB as the TDD user gets −5 degrees. Meanwhile the nulls toward satellite become deeper as the TDD user approaches −5 degree. The algorithm performs the best for the one satellite scenario. For the two and three satellites situations, the beamformer provides slightly better performance for the TDD user in the case of two satellites.
An adaptive antenna system for a BTS according to some embodiments of the invention may reduce inter-system inference from terrestrial cellular systems to independent satellite systems (such as ICS transceivers) that may be operating in the same frequency band as the one used by terrestrial cellular systems. The performance of the adaptive antenna system has been analyzed through simulations. The simulation results have showed that using an adaptive antenna system according to some embodiments of the invention with a vertical uniform linear array of 12 antenna elements along may provide more than 20 dB of isolation toward a satellite regardless of any terrain/building blockage.
In the drawings and specification, there have been disclosed typical embodiments of the invention and, although specific terms are employed, they are used in a generic and descriptive sense only and not for purposes of limitation, the scope of the invention being set forth in the following claims.

## Claims

1. A method of forming a downlink beam in an adaptive antenna system of a communications system, comprising:
receiving vectors of signals including signals transmitted by user terminals of the communications system and signals transmitted by transceivers of an independent communications system;
obtaining spatial information for the user terminals of the communications system and the transceivers of the independent communications system;
generating complex transmitting weights that form spatial nulls directed at the transceivers of the independent communications system; and
transmitting a downlink communications signal using the complex transmitting weights,
wherein obtaining spatial information for the user terminals of the communications system and the transceivers of the independent communications system comprises obtaining the spatial information using the received vectors of signals,
wherein obtaining spatial information comprises generating estimates of an array covariance matrix using the signals transmitted by the user terminals of the communications system and the signals transmitted by the transceivers of an independent communications system, and
wherein generating estimates of the array covariance matrix comprises:
estimating a covariance matrix for the transceivers of the independent communications system;
estimating a covariance matrix for the user terminals of the communications system; and
deriving a composite covariance matrix by combining the covariance matrix for the transceivers of the independent communications system with the covariance matrix for the user terminals of the communications system.

2. The method of claim 1, wherein receiving vectors of signals comprises receiving vectors of signals using an antenna comprising an array of antenna feed elements, and wherein a respective one of the complex transmitting weights is associated with a corresponding one of the antenna feed elements.

3. The method of claim 1, wherein obtaining spatial information for the user terminals of the communications system comprises performing pseudo uplink beamforming, and wherein obtaining spatial information for the transceivers of the independent communications system comprises using blind spatial channel estimation.

4. The method of claim 1, wherein generating estimates of the array covariance matrix comprises obtaining estimates of spatial signatures of the user terminals of the communications system using a preamble or pilot signal reference.

5. The method of claim 1, further comprising:
performing frequency transformation of the covariance matrix for the transceivers of the independent communications system to transform a frequency of the covariance matrix for the transceivers of the independent communications system to a frequency of the communications system.

6. The method of claim 1, wherein generating complex transmitting weights comprises generating a set of pseudo uplink beamforming weights by combining the estimates of the spatial signatures of user terminals of the communications system with the composite covariance matrix, and wherein transmitting a downlink communications signal using the complex transmitting weights comprises transmitting a downlink communications signal using the pseudo uplink beamforming weights transmitting weights.

7. The method of claim 1, wherein obtaining spatial information for the user terminals of the communications system and/or the transceivers of the independent communications system comprises obtaining spatial signatures using an angle of arrival estimation algorithm.

8. The method of claim 7, wherein obtaining spatial information for the transceivers of the independent communications system further comprises estimating angle of arrival information of the signals from the transceivers of the independent communications system using a maximum likelihood angle of arrival estimation algorithm, and generating spatial signatures of the transceivers of the independent communications system using the angle of arrival information.

9. The method of claim 1, wherein the communications system is configured to communicate with user terminals using a time division duplex communications protocol.

10. The method of claim 1, wherein obtaining spatial information for the transceivers of the independent communications system comprises generating estimates of a covariance matrix using the signals transmitted by the transceivers of an independent communications system and decomposing the covariance matrix to obtain the principal eigenvector of the covariance matrix or the K eigenvectors corresponding to the K largest eigenvalues, where K represents the number of transceivers of the independent communications system.

11. The method of claim 1, wherein receiving signals transmitted by transceivers of an independent communications system comprises receiving pilot or carrier signals transmitted by the transceivers of the independent communications system.

12. The method of claim 1, wherein receiving signals transmitted by transceivers of an independent communications system comprises receiving beacon signals transmitted by the transceivers of the independent communications system.

13. The method of claim 12, further comprising:
forming an uplink antenna pattern using the complex transmitting weights;
receiving the beacon signals using the uplink antenna pattern; and
measuring strength of the beacon signals received using the uplink antenna pattern, wherein the strength of the beacon signals received using the uplink antenna pattern provides an error signal associated with the complex transmitting weights.

14. The method of claim 13, further comprising adjusting the complex transmitting weights to reduce the strength of the beacon signals.

15. The method of claim 1, wherein the adaptive antenna system comprises a vertical uniform linear array of radiating elements, and wherein generating complex transmitting weights that form spatial nulls directed at the transceivers of the independent communications system comprises generating complex transmitting weights that form spatial nulls directed at angles of arrival associated with the transceivers of the independent communications system.

16. The method of claim 1, wherein the adaptive antenna system comprises a phased array antenna, and wherein generating complex transmitting weights that form spatial nulls directed at the transceivers of the independent communications system comprises generating complex transmitting weights that form spatial nulls directed at azimuth/elevation coordinates associated with the transceivers of the independent communications system.

17. The method of claim 1, further comprising:
receiving the downlink communications signal at a user terminal of the communication system.

18. The method of claim 1, wherein the communications system comprises a terrestrial communications system.

19. A base transceiver system of a communications system, comprising:
an antenna array configured to receive vectors of signals including signals transmitted by user terminals of the communications system and signals transmitted by transceivers of an independent communications system; and
an electronics system coupled to the antenna array and configured to obtain spatial information for the user terminals of the communications system and the transceivers of the independent communications system, to generate complex transmitting weights that form spatial nulls directed at the transceivers of the independent communications system, and to transmit a downlink communications signal over the antenna array using the complex transmitting weights,
wherein the electronics system is configured to obtain spatial information for the user terminals of the communications system and the transceivers of the independent communications system using the received vectors of signals,
wherein the electronics system is configured to generate estimates of an array covariance matrix using the signals transmitted by the user terminals of the communications system and the signals transmitted by the transceivers of an independent communications system, and
wherein the electronics system is configured to generate estimates of the array covariance matrix by estimating a covariance matrix for the transceivers of the independent communications system, estimating a covariance matrix for the user terminals of the communications system, and deriving a composite covariance matrix by combining the covariance matrix for the transceivers of the independent communications system with the covariance matrix for the user terminals of the communications system.

20. The base transceiver system of claim 19, wherein the antenna array includes an array of antenna feed elements, and wherein respective ones of the complex transmitting weights are associated with corresponding ones of the antenna feed elements.

21. The base transceiver system of claim 19, wherein the electronics system is configured to obtain spatial information for the user terminals of the communications system by performing pseudo uplink beamforming, and wherein the electronics system is configured to obtain spatial information for the transceivers of the independent communications system using blind spatial channel estimation.

22. The base transceiver system of claim 19, wherein the electronics system is configured to obtain estimates of spatial signatures of the user terminals of the communications system using a preamble or pilot signal reference.

23. The base transceiver system of claim 19, the electronics system is further configured to transform a frequency of the covariance matrix for the transceivers of the independent communications system to a frequency of the communications system.

24. The base transceiver system of claim 19, wherein the electronics system is configured to generate a set of pseudo uplink beamforming weights by combining the estimates of the spatial signatures of user terminals of the communications system with the composite covariance matrix, and wherein the electronics system is further configured to transmit the downlink communications signal using the pseudo uplink beamforming weights transmitting weights.

25. The base transceiver system of claim 19, wherein the electronics system is configured to obtain spatial signatures for the user terminals of the communications system and/or the transceivers of the independent communications system using an angle of arrival estimation algorithm.

26. The base transceiver system of claim 21 wherein the electronics system is configured to estimate angle of arrival information of the signals from the transceivers of the independent communications system using a maximum likelihood angle of arrival estimation algorithm, and to generate spatial signatures of the transceivers of the independent communications system using the angle of arrival information.

27. The base transceiver system of claim 19, wherein the communications system is configured to communicate with user terminals using a time division duplex communications protocol.

28. The base transceiver system of claim 19, wherein the electronics system is configured to generate estimates of a covariance matrix using the signals transmitted by the transceivers of an independent communications system and to decompose the covariance matrix to obtain the principal eigenvector of the covariance matrix or the K eigenvectors corresponding to the K largest eigenvalues, where K represents the number of transceivers of the independent communications system.

29. The base transceiver system of claim 19, wherein the antenna array is configured to receive pilot or carrier signals transmitted by the transceivers of the independent communications system.

30. The base transceiver system of claim 19, wherein the antenna array is configured to receive beacon signals transmitted by the transceivers of the independent communications system.

31. The base transceiver system of claim 30, wherein the electronics system is configured to form an uplink antenna pattern using the complex transmitting weights, to receive the beacon signals using the uplink antenna pattern, and to measure strength of the beacon signals received using the uplink antenna pattern, wherein the strength of the beacon signals received using the uplink antenna pattern provides an error signal associated with the complex transmitting weights.

32. The base transceiver system of claim 31, wherein the electronics system is configured to adjust the complex transmitting weights to reduce the strength of the beacon signals.

33. The base transceiver system of claim 19, wherein the antenna array comprises a vertical uniform linear array of radiating elements, and wherein the electronics system is configured to generate complex transmitting weights that form spatial nulls directed at angles of arrival associated with the transceivers of the independent communications system.

34. The base transceiver system of claim 19, wherein the antenna array comprises a phased array antenna, and wherein the electronics system is configured to generate complex transmitting weights that form spatial nulls directed at azimuth/elevation coordinates associated with the transceivers of the independent communications system.

35. A communications system, comprising:
user terminals configured to generate communication signals in a first frequency band; and
a base transceiver system comprising:
an antenna array configured to receive vectors of signals including signals transmitted by the user terminals in the first frequency band and signals transmitted by transceivers of an independent communications system in a second frequency band and an electronics system coupled to the antenna array and configured to obtain spatial information for the user terminals of the communications system and the transceivers of the independent communications system, to generate complex transmitting weights that form spatial nulls directed at the transceivers of the independent communications system, and to transmit a downlink communications signal to at least one of the user terminals of the communications system over the antenna array using the complex transmitting weights,
wherein the electronics system is configured to obtain spatial information for the user terminals of the communications system and the transceivers of the independent communications system using the received vectors of signals,
wherein the electronics system is configured to generate estimates of an array covariance matrix using the signals transmitted by the user terminals of the communications system and the signals transmitted by the transceivers of an independent communications system, and
wherein the electronics system is configured to generate estimates of the array covariance matrix by estimating a covariance matrix for the transceivers of the independent communications system, estimating a covariance matrix for the user terminals of the communications system, and deriving a composite covariance matrix by combining the covariance matrix for the transceivers of the independent communications system with the covariance matrix for the user terminals of the communications system.

36. The communications system of claim 35, wherein the antenna array includes an array of antenna feed elements, and wherein respective ones of the complex transmitting weights are associated with corresponding ones of the antenna feed elements.

37. The communications system of claim 35, wherein the electronics system is configured to obtain spatial information for the user terminals of the communications system by performing pseudo uplink beamforming, and wherein the electronics system is configured to obtain spatial information for the transceivers of the independent communications system using blind spatial channel estimation.

38. The communications system of claim 35, wherein the electronics system is configured to obtain estimates of spatial signatures of the user terminals of the communications system using a preamble or pilot signal reference.

39. The communications system of claim 35, wherein the electronics system is further configured to transform a frequency of the covariance matrix for the transceivers of the independent communications system to a frequency of the communications system.

40. The communications system of claim 35, wherein the electronics system is configured to generate a set of pseudo uplink beamforming weights by combining the estimates of the spatial signatures of user terminals of the communications system with the composite covariance matrix, and wherein the electronics system is further configured to transmit the downlink communications signal using the pseudo uplink beamforming weights transmitting weights.

41. The communications system of claim 35, wherein the electronics system is configured to obtain spatial signatures for the user terminals of the communications system and/or the transceivers of the independent communications system using an angle of arrival estimation algorithm.

42. The communications system of claim 41, wherein the electronics system is configured to estimate angle of arrival information of the signals from the transceivers of the independent communications system using a maximum likelihood angle of arrival estimation algorithm, and to generate spatial signatures of the transceivers of the independent communications system using the angle of arrival information.

43. The communications system of claim 35, wherein the communications system is configured to communicate with user terminals using a time division duplex communications protocol.

44. The communications system of claim 35, wherein the electronics system is configured to generate estimates of a covariance matrix using the signals transmitted by the transceivers of an independent communications system and to decompose the covariance matrix to obtain the principal eigenvector of the covariance matrix or the K eigenvectors corresponding to the K largest eigenvalues, where K represents the number of transceivers of the independent communications system.

45. The communications system of claim 35, wherein the antenna array is configured to receive pilot or carrier signals transmitted by the transceivers of the independent communications system.

46. The communications system of claim 35, wherein the antenna array is configured to receive beacon signals transmitted by the transceivers of the independent communications system.

47. The communications system of claim 46, wherein the electronics system is configured to form an uplink antenna pattern using the complex transmitting weights, to receive the beacon signals using the uplink antenna pattern, and to measure strength of the beacon signals received using the uplink antenna pattern, wherein the strength of the beacon signals received using the uplink antenna pattern provides an error signal associated with the complex transmitting weights.

48. The communications system of claim 47, wherein the electronics system is configured to adjust the complex transmitting weights to reduce the strength of the beacon signals.

49. The communications system of claim 35, wherein the antenna array comprises a vertical uniform linear array of radiating elements, and wherein the electronics system is configured to generate complex transmitting weights that form spatial nulls directed at angles of arrival associated with the transceivers of the independent communications system.

50. The communications system of claim 35, wherein the antenna array comprises a phased array antenna, and wherein the electronics system is configured to generate complex transmitting weights that form spatial nulls directed at azimuth/elevation coordinates associated with the transceivers of the independent communications system.