Method for preserving a pointing direction of an antenna in the axis of a satellite, characterized in that it executes at least the following steps: for i varying from 1 to n, n being a number of given measurements, measuring a value of a tracking signal SNR i received from the satellite for a given time, expressing the value of the tracking signal SNR i as a function of parameters K 1 and K 2 related to the pattern of the antenna, of a time period j, and of a value of the signal-to-noise ratio SNR 0 obtained without voluntary off-pointing of the antenna, SNR i = K 1. ε z − dτ ∑ j = 1 i − 1 ε zJ. 2 + K 2. ε y − dτ ∑ j = 1 i − 1 ε yJ. 2 + SNR 0 with: ε z and ε y being values of deviations to the satellite, i varying from 1 to n, being the index of an SNR measurement in a given time window, j being the index of the time period preceding that of an index i, all the preceding measurements are summed, ε̇ zj and ε̇ yj being the angular velocities for at least two axes y and z that are different from the pointing axis of the antenna pointing towards the satellite, dτ being a sampling period, defining a cost function from the value SNR i and from its estimated value SNR ^ i SNR ^ i = K 1. ε z ^ − dτ ∑ j = 1 i − 1 ε zJ. 2 + K 2. ε y ^ − dτ ∑ j = 1 i − 1 ε yJ. 2 + SNR 0 Err i = SNR ^ i − SNR J = ∑ j = 1 n p i. Err i 2 p i is a weighting factor chosen by the operator, finding the optimal parameters ε y ^ and ε z ^ in the least-squares sense, i.e. those minimizing the function J in each time interval, using the result of the comparison of the estimated values to the reference value to preserve the pointing direction. for i varying from 1 to n, n being a number of given measurements, measuring a value of a tracking signal SNR i received from the satellite for a given time, expressing the value of the tracking signal SNR i as a function of parameters K 1 and K 2 related to the pattern of the antenna, of a time period j, and of a value of the signal-to-noise ratio SNR 0 obtained without voluntary off-pointing of the antenna, SNR i = K 1. ε z − dτ ∑ j = 1 i − 1 ε zJ. 2 + K 2. ε y − dτ ∑ j = 1 i − 1 ε yJ. 2 + SNR 0 with: ε z and ε y being values of deviations to the satellite, i varying from 1 to n, being the index of an SNR measurement in a given time window, j being the index of the time period preceding that of an index i, all the preceding measurements are summed, ε̇ zj and ε̇ yj being the angular velocities for at least two axes y and z that are different from the pointing axis of the antenna pointing towards the satellite, dτ being a sampling period, ε z and ε y being values of deviations to the satellite, i varying from 1 to n, being the index of an SNR measurement in a given time window, j being the index of the time period preceding that of an index i, all the preceding measurements are summed, ε̇ zj and ε̇ yj being the angular velocities for at least two axes y and z that are different from the pointing axis of the antenna pointing towards the satellite, dτ being a sampling period, defining a cost function from the value SNR i and from its estimated value SNR ^ i SNR ^ i = K 1. ε z ^ − dτ ∑ j = 1 i − 1 ε zJ. 2 + K 2. ε y ^ − dτ ∑ j = 1 i − 1 ε yJ. 2 + SNR 0 Err i = SNR ^ i − SNR J = ∑ j = 1 n p i. Err i 2 p i is a weighting factor chosen by the operator, finding the optimal parameters ε y ^ and ε z ^ in the least-squares sense, i.e. those minimizing the function J in each time interval, using the result of the comparison of the estimated values to the reference value to preserve the pointing direction.