1. A method of multipath error suppression in a satellite navigation receiver, comprising steps of navigation satellites signals searching, receiving and processing complex signals from each satellite, tracking found signals with a following-up open loop and determining coordinates, receiver velocity and exact time based on measurements of direct and reflected signal delay and Doppler frequency, the method comprising: on the basis of navigation parameters including coordinates, receiver velocity and exact time, and also information on satellite movement obtained on the basis of measured time and received ephemerides, calculating delay and Doppler frequency predictions for orientation of a two-dimensional accumulated power grid according to delay and Doppler frequency, which grid formed on the basis of received and processed complex signals; wherein, in the course of grid elements forming, a summation of complex signal is performed in a coherent accumulation interval synchronously with data bits of the given satellite, and the coherent accumulation time is matched with a bit length; further on the basis of processing the accumulated power grid including its centering, calculating corrections for delay and frequency predictions, and calculating delay and Doppler frequency measurements as the sum of predictions and corresponding corrections, then extrapolating obtained measurements for the preset moment of time and calculating pseudorange and pseudovelocity measurements, and determining a receiver coordinates and velocity on the basis of the measurements obtained for all satellites; wherein calculation of delay and Doppler frequency corrections comprises steps of: calculating single-path signal corrections using the centered accumulated power grid; determining, whether a multipath is presented; and in case of positive result of this determination, performing improvement of corrections using weights; wherein, performing improvement of corrections using weights is carried out by forming a likelihood function on the basis of a centered accumulated powers vector, and correction is performed according to the found likelihood function global maximum.