1.一种GNSS接收机的抑制多径的方法，其特征在于，包括以下步骤：S1、利用GNSS接收机接收多颗GNSS卫星信号；

S2、对步骤S1接收的GNSS卫星信号进行射频处理；

S3、对步骤S2处理后的信号进行基带预处理；

S4、采用TK算法对步骤S3预处理后的信号进行处理，得到多径的条数和每条多径对应的时延，作为先验信息；

S5、根据步骤S4得到的先验信息对各条多径的幅度参数进行预估计；

S6、采用迭代方法对各条多径的相关函数进行极大似然估计，得到直达波信号的相关函数；

S7、对抑制多径后的相关函数进行导航定位解算。

2.如权利要求1所述的GNSS接收机的抑制多径的方法，其特征在于，所述步骤S4具体包括以下分步骤：S41、获取离散相关函数的多个点的横坐标、纵坐标，并以此构建相关函数的完整包络；

S42、对离散相关函数的所有点进行TK运算，得到对应的TK曲线；

S43、设定幅度门限，将步骤S42得到的TK曲线上所有点的幅度依次与幅度门限进行比较，选取高于幅度门限的极值点的个数作为多径的条数；

S44、根据步骤S43选取的每一个极值点所对应的时刻，确定该极值点对应的时延的值，得到所有多径的时延值；

S45、保存多径的条数和每条多径对应的时延，作为先验信息。

3.如权利要求1所述的GNSS接收机的抑制多径的方法，其特征在于，所述步骤S5具体包括以下分步骤：S51、比较每条多径对应的时延大小，将时延最小的多径判断为直达波，将其它多径按照时延从小到大分别判断为第2条至最后一条多径；

S52、根据直达波和各条多径的相关函数的重叠程度，对各条多径的幅度参数进行预置。

4.如权利要求1所述的GNSS接收机的抑制多径的方法，其特征在于，所述步骤S6具体包括以下分步骤：S61、找到合成信号相关函数的幅度最高点，对该点的实部、虚部求反正切，得到直达波的相位估计值；

S62、根据直达波的时延值、幅度估计值及步骤S61得到的相位估计值重构直达波的相关函数；

S63、利用合成信号相关函数减去步骤S62重构的直达波的相关函数所得到的差值估计出相位，结合第2条径的时延值、幅度估计值对第2条径的相关函数进行重构；

S64、利用合成信号相关函数减去步骤S63重构的第2条径的相关函数所得到的差值估计出相位，利用该差值中的能量最大点的幅度作为该差值对应的直达波的幅度估计值，根据得到的相位估计值、幅度估计值及直达波的时延值再次重构直达波的相关函数；

S65、根据多径的条数判断是否已经穷尽所有多径；若是，则步骤结束；若否，则执行下一步骤；

S66、利用合成信号相关函数减去重构的直达波的相关函数、已重构径的相关函数，根据所得到的差值估计出相位，结合下一条径的时延值、幅度估计值对下一条径的相关函数进行重构；

S67、利用合成信号相关函数减去步骤S66重构的径的相关函数所得到的差值估计出相位，利用该差值中的能量最大点的幅度作为该差值对应的直达波的幅度估计值，根据得到的相位估计值、幅度估计值及直达波的时延值再次重构直达波的相关函数；

S68、根据多径的条数判断是否已经穷尽所有多径；若是，则步骤结束；若否，则循环重复步骤S66和S67，最终得到直达波信号的相关函数。

5.一种GNSS接收机的抑制多径的系统，其特征在于，包括天线模块，分别与天线模块连接的多个射频处理模块，分别与射频处理模块连接的基带预处理模块、多径预处理模块、多径幅度参数预估计模块、极大似然估计模块和导航定位数据处理模块；

所述天线模块用于接收多颗GNSS卫星信号，并将接收信号分别传输至多个射频处理模块；

所述射频处理模块用于依次对接收信号进行低噪声放大处理、下变频处理和A/D转换处理，并将射频处理后的接收信号传输至基带预处理模块；

所述基带预处理模块用于依次对接收信号进行数字下变频处理和抗干扰处理，并将基带预处理后的接收信号传输至多径预处理模块；

所述多径预处理模块用于采用TK算法对接收信号进行预处理，得到多径的条数和每条多径对应的时延，作为先验信息，并将预处理后的接收信号传输至多径幅度参数预估计模块；

所述多径幅度参数预估计模块用于根据先验信息对各条多径的幅度参数进行预估计，并将幅度估计值传输至极大似然估计模块；

所述极大似然估计模块用于采用迭代方法对各条多径的相关函数进行极大似然估计，得到直达波信号的相关函数，并传输至导航定位数据处理模块；

所述导航定位数据处理模块用于对抑制多径后的相关函数进行导航定位解算。

6.如权利要求5所述的GNSS接收机的抑制多径的系统，其特征在于，所述多径预处理模块采用TK算法对接收信号进行预处理，得到多径的条数和每条多径对应的时延，作为先验信息，具体为：采用平行相关器采样单元获取离散相关函数的多个点的横坐标、纵坐标，并以此构建相关函数的完整包络；

对离散相关函数的所有点进行TK运算，得到对应的TK曲线；

设定幅度门限，将TK曲线上所有点的幅度依次与幅度门限进行比较，选取高于幅度门限的极值点的个数作为多径的条数；

根据选取的每一个极值点所对应的时刻，确定该极值点对应的时延的值，得到所有多径的时延值；

保存多径的条数和每条多径对应的时延，作为先验信息。

7.如权利要求5所述的GNSS接收机的抑制多径的系统，其特征在于，所述多径幅度参数预估计模块根据先验信息对各条多径的幅度参数进行预估计，具体为：比较每条多径对应的时延大小，将时延最小的多径判断为直达波，将其它多径按照时延从小到大分别判断为第2条至最后一条多径；

根据直达波和各条多径的相关函数的重叠程度，对各条多径的幅度参数进行预置。

8.如权利要求5所述的GNSS接收机的抑制多径的系统，其特征在于，所述极大似然估计模块采用迭代方法对各条多径的相关函数进行极大似然估计，得到直达波信号的相关函数，具体为：找到合成信号相关函数的幅度最高点，对该点的实部、虚部求反正切，得到直达波的相位估计值；

根据直达波的时延值、幅度估计值及相位估计值重构直达波的相关函数；

利用合成信号相关函数减去重构的直达波的相关函数所得到的差值估计出相位，结合第2条径的时延值、幅度估计值对第2条径的相关函数进行重构；

利用合成信号相关函数减去重构的第2条径的相关函数所得到的差值估计出相位，利用该差值中的的能量最大点的幅度作为该差值对应的直达波的幅度估计值，根据得到的相位估计值、幅度估计值及直达波的时延值再次重构直达波的相关函数；

根据多径的条数判断是否已经穷尽所有多径；若是，则步骤结束；若否，则执行下一步骤；

利用合成信号相关函数减去重构的直达波的相关函数、已重构径的相关函数，根据所得到的差值估计出相位，结合下一条径的时延值、幅度估计值对下一条径的相关函数进行重构；

利用合成信号相关函数减去重构的径的相关函数所得到的差值估计出相位，利用该差值中的的能量最大点的幅度作为该差值对应的直达波的幅度估计值，根据得到的相位估计值、幅度估计值及直达波的时延值再次重构直达波的相关函数；

根据多径的条数判断是否已经穷尽所有多径；若是，则步骤结束；若否，则对下一条多径进行估计，最终得到直达波信号的相关函数。