Prior Art It is worldwide known a method for satellite identification, based on radiating of unique satellite identification (ID) code that is received from the terrestrial equipment of the satellite antenna, and after decoding it is determine the radiating satellite.

A disadvantage of this familiar method is a long time for identification and the possibility for mistakes appearing in the ID codes, which are due to changes in the environment conditions, fluctuations of the signals in the electronic devices of the antenna and etc. Furthermore, the method requires specific hardware devices for broadcasting, receiving and decoding of the ID signals. These specific hardware modules leads to complication of the receiving equipment that implements the method and decreases the reliability of the telecommunications.

From the patent RU No. 2 14 675 is known a method for defining the position of a broadcasting source in case of multi-beam radiation. The method consists of : - acceptance of the signals from an N-element phased-array antenna, - formation of scanning signals in a given angle sector, -determination of the weight coefficients and multiplication with the respective scanning signals, - memorizing of the received signals, and selecting from them of the angle values with the maximum strength of the received signal.

This familiar method is intended for implement of a mobile radio communication and does not provide identification of the broadcasting source.

Summary of the invention The purpose of the present invention is to create a method for identification of set of satellites located on GEO orbit, that ensuring a correct identification of two or more satellites by simplifying of the transmitting and receiving part of the telecommunication system and respectively increasing the reliability of the telecommunications.

The purpose is achieved with creation of a method for identification of set of satellites, situated on GEO orbit, based on the received signals radiated from the satellites, and characterized by the following step : - determination of the antenna beams signal level - detection of the local maximums according to the antenna structure, defined by the number of beams and their coordinates in the antenna plane ; - approximation of the GEO orbit using the strength of the received signals wherein, the satellites identification is carried out by local maximums located on or near to the approximated GEO curve, and according to the information of the geographic position of each satellite and antenna orientation.

In one variant of the execution of the method for satellites identification are used probability coefficients for the projection of each satellite in a given area of the antenna plane.

. It is expedient, the information of the. geographic position of each satellite to be corrected by Azimuth, Elevation and Roll angles according to the antenna installation inaccuracy.

In addition, it is suitable after determination of the signal level of all the beams to introduce a correction, reflecting to the signal/noise ratio. In this case it is expedient, the signals corresponding to the strength of the received signals, to be amplified by square or bi-square law.

The advantages of the'method for identification of set of satellites, located on GEO orbit, consists of correct identification of each satellite from the set without the need using complex and expensive decoding modules in the receiving satellite equipments. By implementation of the method according to the invention, the

reliability of the telecommunications is increased ; since by decreasing the number of elements comprising the telecommunication system, failure probability of given elements is reduced.

Brief description of the drawings The method for identification of set of satellites, located on GEO orbit, is illustrated by the enclosed drawings :.

Figure 1-flowchart representing the algorithm for identification of set of satellites by the method according to the invention; Figure 2-chart representing the concept"the, first adjacent set of a given beam" Figure 3-flowchart representing the means of calculation of the satellite statistical probability projections in given area of the antenna plane ; Figure 4-chart representing the satellites probability projection in given area of the antenna plane ; Figure 5-chart representing all the received signals by the antenna's beams; Figure 6-flowchart representing local maximums detecting algorithm ; Figure 7-flowchart representing the approximation of the GEO orbit ; Figure 8-chart representing the GEO orbit, approximated with a first order curve (straight line) ; Figure 9-chart representing the GEO orbit, approximated with a curve of second order ; Figure 10-flowchart representing the algorithm for detecting of the local maximums, situated in the range of the approximated GEO orbit ; Figure 11-flowchart representing the satellites identification algorithm ; Figure 12a, 12b and 12c-flowcharts representing the algorithm for satellites identification in critical circumstances ; Figure 13a and 13B-flowcharts representing the block"Finding of missing satellites"in case of critical circumstances ; Figure 14-chart representing the result of satellites identification method.

Examples for execution of the invention The method for identification of set of satellites, located on GEO orbit, according to the invention is illustrated by the algorithm flowchart on Figure 1. Input data, necessary for satellites identification, are variables and constants. The variables carry information for the strength of the received signal from each beam of the antenna and generally shown as RSSI-Received Signal Strength Indicator. The constants or hardware dependent data are: -satellites location-information about the satellites longitudes ; - antenna structure-information about the number of the beams and their coordinates in the antenna plane ; -beams anrangement (LINK_LIST)-it is an array containing indexes of beams from the first adjacent set. Figure 2, illustrates the concept"first adjacent set"of the given beam. On Figure 2, the beams in this first adjacent set are marked in black color ; - additional information about the arrangement of the beams (ENDADJSET)-it is an array filled with count of the beams from the first adjacent set; - projections statistical probability of the set of satellites in the antenna plane (PROB_PROJ)-defines the probability of a given beam to be directed towards a given. satellite. This probability for the desired set of satellites is calculated in advance 'or-a given geographical region, which the antenna will operate in. Determination of the statistical probability is illustrated on Figure 3. Variables as antenna's longitude and latitude, and constants as satellites location, antenna's structure, Earth radius (6378.16 km) and GEO orbit altitude (for example 35786.43 km) are inputs for the block"The antenna orientation". This block calculates satellite's Azimuth and Elevation angles in the antenna plate. The purpose is orientation of boresight (beam oriented in direction at angle 90° with respect to the antenna plane) of the antenna (for multi-beams antenna) or orientation of the antenna (in single beam antenna). to a given satellite. After these operations, the block calculates the probability of a given satellite to fall into a given beam of the antenna. Block"PROB-PROJ"generates

weight coefficients, which are proportional of the probability a given beam to be directed to a given satellite. Figure 4 shows the distribution of probabilitties in the antenna plane, and the areas with lighter spice reflect greater projection probability of a given satellite.

The input data (Figure 1) is processed in the block"Input RSSI processing". For suppressing of the noise and amplification of the useful signal level, it is suitable to normalize the signal. For this purpose the RSSI values are subjected to square or bi- square law amplification.

Further, the normalized values are used for finding of all the local maximums. This action is performed by consecutively processing of all the beams (Figure 6). It checks whether the RSSI values of a given beam are bigger than RSSI values of the adjacent beams. If that condition is true, this beam presents a local maximum. If the RSSI values of a given beam are not bigger than beams of the adjacent set, the next beam is analyzed. This operation is repeated until the list of the beams with the local maximums is filled.

Next step after finding of all local maximums is approximation of the GEO orbit.

First order polynomial approximation of least-square method is used. The choice of method for approximation depends on maximal permissible error for a given case.

For the approximation, the normalized RSSI values are used as weight coefficients of sought function for decreasing the influence of the noise signals, respectively for more precise approximation. Coefficients of the approximated line and the error of approximation are calculated according to preliminary data that contain the maximal admissible error and RSSI values of all the beams. After this operation the algorithm finds the beam that have maximal calculated error and check whether this error is greater than the maximal permission error. If this condition is true, the beam is excluded from the approximate calculation otherwise this beam is included in the approximate calculation. The GEO orbit, approximated with first order curve (straight line) is shown on Figure 8. It is possible to perform aproximation with second order curve, which case is shown on Figure 9.

After approximation of the GEO orbit, the algorithm continues with detecting of local maximums, situated on or in defined range of the approximated GEO orbit (Figure 10). For this purpose are used already calculated coefficients of the approximated GEO orbit and the list of the local maximums. It checks consecutively all the local maximums from the list. The local maximums that are situated at a distance from the approximated GEO curve smaller than the maximal admissible distance are selected. Using this method, the algorithm composes list of local maximums, situated on the approximated GEO orbit.

Besing on the already determinated list of the local maximums, the algorithm realizes satellite's identification that is shown on Figure 11. The local maximums are sorted by their X axis coordinates in the antenna plane. When the number of the local maximums is bigger than 2, the algorithm estimates consecutively the distance between the current set of local maximums and if it is equal of the distance between the desired group of satellites then correlative sum for this current group is calculated. After this procedure, the algorithm selects the group of local maximum that have the biggest calculated correlative sum (maximal correlative sum). This group contains the beams that point to the desired group of satellites.

If described above algorithm identifies count of satellites smaller than desired in the algorithm, begins the algorithm shown in Figure 12a,-12b and 12c. This algorithm like described in Figure 11 aims to. define a group of beams that point to satellites group. Insufficiency of input data like local maximums is complemented with data of type probability. This data consists of beams selected according to the projection probabilities of the given satellites (PROB_PROJ) and situated in defined range of the approximated GEO orbit, called PP valid beams. In case of absence of satellites signals, artificial (non-existing) beams-PP_valid beams are selected, which point to the missing satellites. The algorithm evaluates consecutively the distances between currents local maximums and if it is necessary, the group of beams is complemented by Pu valid beams (Figure 13a and 13b). If the calculated distance is equal of the distance between the corresponding satellites, the correlative sum for

a given group of beams is calculated. After estimation of the correlative sums of all the groups of beams, the algorithm selects the group that has the biggest calculated correlative sum (maximal correlative sum). This group contains beams that point to the desired satellite's group and beams that point to the probable position of the missing satellites.

Determination of the called PP valid beams is performed by the algorithm shown of Figure 13a and 13b (Block"Finding of missing satellites"). Satellites will not received signal for, the algorithm searchs fictitious (non-existing) beams that pointing them. The minimum count of satellites, which should received signal from is two. The fictitious beams must to be located on or in defined range of the approximated GEO orbit and to have a probability projection value (PROB_PROJ) for the given satellites.

On the base of the local maximums (beams that point to satellites, which are received signals from) and Pu valid beams (beams that point to missing satellites) is defined group of beams, which correspond to given conditions. And they are that the distance between Pu valid beams and the other beams (the local maximums) that point to satellites should be in given range, defined from the minimal and maximal permissible distance between each two satellites in the antenna plane.

Figure 14 shows the results of the identification algorithm. In this case, the count of desired satellites is three.