多径信号对卫星接收机的影响及时域抑制算法研究.pdf

C o m m u n i c a t i o n&I n f o r m a t i o n T e c h n o l o g y N o.4.2 0 2 3通信与信息技术2 0 2 3 年第4 期(总第2 6 4 期)多径信号对卫星接收机的影响及时域抑制算法研究朱奎宝，郭广源，杜云，温紫晴，邓承宾，康浩楠河北科技大学电气工程学院，河北石家庄0 5 0 0 1 8摘 要：针对卫星导航系统多径误差问题，建立了多径效应模型，采用二象限反正切鉴相函数和E M L 鉴相器，在多径环境下分析了载波跟踪环和伪码跟踪环受到的影响。在抑制多径算法方面采用M E D L L 和W R E L A X 算法抑制多径信号并进行仿真，并就两种算法对多径信号抑制的性能进行对比分析。实验表明，不同的多径信号对接收机的影响不同，对接收机解算真实信号造成的影响程度不同。在抑制多径方面，W R E L A X 算法要比M E D L L 算法的抑制多径信号效果更好，估计误差更小。关键词：卫星导航系统；多径信号；接收机；M E D L L,W R E L A X中图分类号：T N 9 2 7+.2 1文献标识码：A文章编号：1 6 7 2-0 1 6 4(2 0 2 3)0 4-0 0 0 7-0 71 引言全球卫星导航系统(G l o b a l N a v i g a t i o n S a t e l l i t e S y s t e m,G N S S)是各国极其重要的基础设施，在国防安全、经济建设和社会发展等方面起着重要的作用。随着社会的不断进步，用户对导航定位的需求不断增加，对其抗干扰和抑制多径效应的性能要求也越来越高。城市建筑、山川河流等产生的峡谷效应，是产生多径效应的重要因素四。为了能够得到全地域高精度导航定位服务信息，需要接收机具有检测干扰信号和有效抗干扰的能力，有针对性地消除干扰所带来的误差。通过仿真手段构造多径信号，与卫星信号叠加，形成卫星多径干扰效应，采用多径估计延迟锁定环(M u l t i p a t hE s t i m a t i n g D e l a y L o c k L o o p,M E D L L)和基于W R E L A X(W e i g h t e d F o u r i e r T r a n s f o r m a n d R E L A X a t i o n)参数估计的时域多径干扰抑制方法，对卫星多径干扰信息进行识别、分离，并对两种算法对多径信号的识别和分离性能进行对比分析。2 多径效应多径效应是指卫星接收机除了接收卫星经过直线传播的电磁波信号外，还收到天线附近的高楼、山川、河流、树木等对卫星信号的反射信号，甚至当直达信号受到阻碍时，接收机只能接收多径信息 2 。多径误差与其他误差不同，在不同时间和不同环境下，表现出不同的特点，即使同样的多径干扰在不同接收机上也会产生不同的影响。在测量定位中，多径效应误差的大小，取决于反射物离接收机天线的距离和反射系数及卫星信号的方向等多因素。对于双天线卫星导航系统，基线两端的多径效应不具备相关性，无法用差分技术来消除其影响。假设卫星的直射波信号s q(t)为：s a(1)=D(t)c(t)c o s(o o.t)(1)式(1)中：D(1)为卫星信号的导航电文；c(1)为C/A 码；w.为载波频率。采用两路信号来构造受多径信息干扰的卫星信号模型，其中一路为直达信号，另外一路为经由四周物体反射后的多径信号。两路信号的载波频率相同，幅值、延时和相位与直达信号不同。则卫星接收机收到的中频信号的模型可表示为：x(t)=a m D(t)c(t-r m)c o s(o。(t-t m)+m)+(2)a?D(t)c(t-t?)c o s(o。(t-x?)+q)式(2)中：m、4 m、T m 为直达信号的幅值、初相位和时延；1、1、T?表示为多径信号的幅值、相位和时延，w.为信号的载波频率。3 多径信号对卫星接收机的影响多径与直达信号相比有明显的时间延迟，码相位和载波相位发生偏差，进而造成定位误差。下文针对两路多径模型对卫星接收机中的伪码跟踪环路和载波跟踪环路造成的影响进行讨论 3。收稿日期：2 0 2 2 年1 1 月1 7 日；修回日期：2 0 2 3 年6 月6 日基金项目：中央军委科学技术委基础加强基金项目(2 0 2 0-J C J Q-J J-2 1 7)7通信与信息技术2 0 2 3 年第4 期(总第2 6 4 期)3.1 多径信号对伪码跟踪环的影响以双径模型(2)式为例，接收机早晚码相关器的间隔时间为d,经过载波环路解调后成为基带信号，再通过积分累加模块与本地码进行相关，得到超前、即时、滞后的相关值为 3 :(da mI E=Rc o s(q m)+T+22,(da?Rc o s(q m+0)T-T?+22.(3)a m1,=R(r)c o s(q m)+2?R(x-x?)c o s(q m+0)2da mc o s(q m)+L=RT-22.dq?Rc o s(q m+0)T-T?-22,(4)(5)d a mQ=R|s i n(q m)+T+22,dR|x-?+s i n(q m+0)22,(6)a mQ=R(r)s i n(q m)+2?R(r-x?)s i n(q m+0)2da mQ=R|T-s i n(q m)+22,da?Rs i n(q m+0)t-T?-22,(7)(8)其中，公式(3)、(4)、(5)中I g、l p、为同相支路超前、即时、滞后的相关值，公式(6)、(7)、(8)中Q g、Q p、Q L 为正交支路超前、即时、滞后的相关值，t=m-T m=和m=4 m-4 m 分别表示多径信号影响下的跟踪环路中的码跟踪误差和载波相位跟踪误差，t?为多径相比于直达信号所产生的码延迟；为多径相比于直达信号所产生的载波相位跟踪误差，R()为C/A 码的自相关函数。多径延迟t?为0.5 c h i p s,多径与直达信号之间的幅度比值a 为0.5,图1(a)和图1(b)分别表示了同相和反相多径对相关峰的影响。在卫星接收机的伪码跟踪环路中，多径信号先是通过相同相多径情况下的相关函数1.5合成信号多径信号-真实信号1.0归一化幅值0.50.0-20码偏移(码片)(a)同相反向多径情况下相关函数21.0-归一化幅值0.50.0-0.5-真实信号多径信高合成信号-202码偏移(码片)(b)反相图1 存在同相和反相多径时的相关函数超前减滞后包络鉴相器鉴相函数-d=1.0 c h i pd=0.8 c h i pd=0.4 c h i p鉴相器输出(chips)0-0.0 10.-1 0.0 0a.0 20.00.2-2+-4-2024输入误差图2 含有多径信号的超前减滞后鉴相曲线关器得到相关值，在伪码鉴相器输入相关值，进而得到码相位误差。以超前减滞后包络鉴相器为例，其表达式为式(9):D=(1 2+0)-(T 2+0 2)(9)图2 所示的是在相关器不同间隔下C/A 码调制的卫星信号的鉴相曲线。在多径环境下会造成伪码相位延迟，鉴相曲8运营一线多径信号对卫星接收机的影响及时域抑制算法研究线产生畸变，曲线过零点时偏离原点。由图2 可以看出，由于受到相关器间隔大小的影响，过原点的偏移程度随着相关器的间隔变化而变化。接下来讨论在不同的多径延迟下，伪码跟踪环路和载波跟踪环路受到的影响 4-5 。假设已经将载波全部剥离，在忽略直达信号载波的影响，采用E M L 鉴相器，其表达式为公式(1 0):D。r,=I E-I(1 0)对于连续的伪码跟踪环路，令D e m,=0,可以解出伪码相跟踪误差和多径参数之间的关系间。a r,c o s(Q)0 T?T1+a c o s(o)a d c o s()T w T?T y2T=(1 1)a c o s(o)d dT-T?+T 22公式(1 1)中，T。为码片的宽度，t 为受多径信号影响的码相跟踪误差，?为多径相对于直达信号的载波相位误差，T 1 为多径延迟。假设只有一路多径信号干扰，在多径与直达信号之间的幅度比值为=1/2,d=T c 的情况下，图3 所示?分别取0.1,0.3,0.6 和0.7 时码相跟踪误差与多径延迟之间的关系。由图3 可以看出，当相关器的间隔d=T?时，多径相比于直达信号的延迟T?1.5 c h i p s 时，会造成严重的码相跟踪误差；当T i 1.5 c h i p s 时，相关性就消失了，此时多径信号对跟踪环路几乎不存在影响。码相跟踪误差与多径延迟的关系0.3-0.2-码相跟踪误差/Tc0.10.0-0.1-0.2-00.1*p i0.3*p i0.6*p i0.7*p i多径延迟/T C2图3 码相跟踪误差与多径延迟的关系图0.2 5-0.2 0 0.1 5码相跟踪误差与多径相位的关系码相跟踪误差/Tc0.1 00.0 5-0.0 0-0.0 5-0.1 0-0.1 5-0.2 0-0.2 500.1*d0.3*d0.6*d0.7*d2多径相位/r a d6图4 码相跟踪误差与多径相位之间的关系图4 所示的幅度比=1/2,d=T c,T?分别取0.1 d,0.3 d,0.6 d 和0.7 d 时码跟踪误差与多路径相位之间的关系。由图4 可以看出，码相跟踪误差与多径相位之间是符合周期变化的。图5 所示的相位1=0.3,d=T e,T?分别取0.1 d,0.3 d,0.6 d 和0.7 d 时，码跟踪误差与多径信号之间的幅度比的关系。由图5 可以看出，当多径相对于直达信号的幅值比变大时，码相跟踪误差也会随之变大。3.2 多径信号对载波跟踪环的影响在载波跟踪的过程中，通常采用的即时支路是1 p 和Q p,对于连续的跟踪D e r r=0,采用二象限反正切鉴相函数可以得到载波跟踪误差与多径信号各参数之间的关系：1R(t-t?)s i n a m(1 2)o =a r c t a nR(t)+R(x-x?)c o s a m根据式(1 2)假设t=0,j/a m=0.5,d=T o,图6 所示为?分别取0.1,0.3,0.6 和0.7 时的载波相位误差与多径信号延迟之间的关系。由图6 可以看出，由于多径信号的延迟不断增加，直达信号和多径信号之间的相关程度在逐渐变小，载波相位误差也在不断变小。根据(1 2)假设t=0,1/a m=0.5,d=T o,图7 所示为T?分别取0.1 d,0.3 d,0.6 d 和0.7 d 时的载波相位误差与多径信号相对相位之间的关系。由图7 可以看出，载波跟踪环路误差与多径信号相位之间是呈现周期变化的。假设t=0,=0.3,d=T c,图8 所示为T?分别取0.1 d,0.3 d,0.6 d,和0.7 d 时的载波相位跟踪误差与多径相对于直达信号之间的幅度比的关系。由图8 可以看出，多径相对于直达信号的幅度比的变大会导致载波相位误差变大。9通信与信息技术2 0 2 3 年第4 期(总第2 6 4 期)码相跟踪误差与多径幅度的关系0.5 5-0.5 0-0.1*d0.3*d0.4 5-0.6*d0.4 00.7*d码相跟踪误差/Te0.3 5-0.3 0-0.2 5*0.2 0-0.1 5-0.1 0-0.0 5-0.0 0-0.0 5-0.00.51.0多径信号幅度图5 码相跟踪误差与多径幅度比的关系图载波相位跟踪误差与多径延迟的关系0.5 50.5 0-0.4 5-载波相位跟踪误差/rad0.4 00.3 50.3 0-0.2 5-0.2 0-0.1 50.1 00.0 5-0.0 0-0.0 5 40.00.1*d-0.3*d0.6*d-0.7*d0.5老经延迟/T s1.0图6 载波相位跟踪误差与多径延迟的关系图载波相位跟踪误差与多径相位的关系0.5*载波相位跟踪误差/rad0.0-0.5-0.1*d0.3*d0.6*d0.7*d026多径相位/r a d图7 载波跟踪误差与多径相位的关系图载波相位跟踪误差与多径幅度的关系0.4-载波相位跟踪误差/rad0.2-a.0-0.1*d 0.3*d-0.6*d-0.7*d0.00.51.0多径幅度图8 载波相位跟踪误差与多径幅度比的关系图4 多径抑制算法分析4.1 M E D L L 算法假设卫星接收机收到的多径信号模型如下：Mx(t)=a m c(t-i m)c o s(w t+m)+e(t)m=0(1 3)构造最大似然函数估计信号x(t),假设估计的信号为s(t),则可通过最大似然函数，使得式(1 3)得到最小值，即式(1 4):Ms(1)=Z a m c(t-r m)c o s(o t+)(1 4)m=0估计信号的多路径数目为M,对于第m 路信号的参数可以通过对L(m,T m,m)求偏导数，并且令偏导数等于0,从而得D L(a m.T m 4 m),a L(a m T m.4 m)到估计参数。即令=0,=0,d a md a mO L(a m.t m(p m)=0,从而得到：0 4 mJ_应t=m aR eR.(t)-a,R(t-t:)ee-i om-(1 5)应R.(r h)-a R(t m-t!)e i n e-mm =R ei=0i m儿(1 6)MR.(r h)-Z a/R(r m-t)e m im =a r gi=0i m(1 7)式中，R(t)为C I A 码的自相关函数，归一化幅值的最大值是1,相对延时是0,作为基准相关函数，R x(T)为实测相关函数。M E D L L 利用循环迭代来获取不同路径的参数估计值，当信号中只存在一路直达信号和一路多径信号?时，其循环迭代过程如下 8 1:(1)当只存在一路多径信号时，则可得到R o(T)=R?(T),将其带入估计模块后，可以得到第一组估计参数a。,。1 01运营一线多径信号对卫星接收机的影响及时域抑制算法研究根据这组参数可以得到n?()的相关函数R o(T)=a o R(t-)e 0)。(2)将实测相关信号R?(T)减去(1)中所得的R o(T),假设得到的差值为第二路信号的n?相关函数，即R()=R,(T)-R o(T),将其代入估计模块后，可以得到第二路多径信号2 估计参数a?,。根据得到的参数可以重构n?()的相关函数R?(T)-a?R(t-分)e 啊。(3)用实测函数减去(2)中所得的R?(T)的结果，假设为?()的相关函数，即R o(T)=R x(T)-R 1(T)代入估计模块，得到估计参数a o,o,根据估计参数构造?()的相关函数R。(t)(4)重复(2)和(3)的操作过程，即可完成两路信号的迭代过程，直到不满足设定的循环条件，跳出循环。此迭代过程可推广到多路径干扰情况，但是在传输过程中，信号随着折射次数增加而变弱，故通常情况下只考虑一路多径信号。4.2 W R E L A X 算法W R E L A X 算法是一种基于最小二乘准则的有效地抑制多径干扰的算法，优点是将多维转换成一维进行优化，可有效地减少运算量9。将u(t)进行载波环路解调后可得式(1 8):Mu(t)=Z a,c(t-m)c o s(4 m)+e(t)(1 8)m=0对解调后的u(1)进行H i l b e r t 变换并进行数据采样后得到：Mx(n)=a,c(n-i m)e m+e(n),n=1,2 Nm=0(1 9)将式(1 9)中的a m 和e j?m 合成为新的变量，用m 表示，则待估计参数为a m,T m。同时对其进行F F T 变换，则有：MX(p)=C(p)Z a m e x p(j o m p)+E(p)B=0(2 0)式(2 0)中p 的范围为(-N/2,N/2-1),X(p)、C(p)、E(p)为x(n)、c(n)、e(n)的F F T 变换。为了得到估计值，可以利用以下非线性最小二乘代价函数最小化求得式(2 1):K/2-1MQ(M=2X(k)-C(k)Z a m e f)m,mm=0,p=-K/2m=0(2 1)假设 a m,0 m M=0,M p 已知或者已被估计出，定义式(2 2):MR m=R-a m R a(o)p=p m(2 2)其中x d i a g R(-K/2),R(-K/2-1),R=R(一N/2),R(N/2-1),最小化上述代价函数可得1 0:1 2|a (a.)(X R)=a r g m a x(2 3)ma”(0 m)(X R m)m =I I x I Il o m=c i m(2 4)其中，止临表示为F-范数，m 为可以通过找周期图1 2.|a (w m)(X*R m)I 的最大值对应的位置得到，而这可通过l a(a m)(X R m)l一维补零F F T 运算实现，a m 为-的复幅度。I X I当只存在一路真实信号和一路多径信号时，可通过如下几个步骤实现：(1)假设只存在卫星直达信号，通过R 再根据(2 3)式和(2 4)式即可估计出o,W o。(2)利用(1)中的a o,O o,根据(2 2)式计算得到R?,再利用R?通过(2 3)式和(2 4)式即可估计出a 1,1。(3)同理，利用(2)中的a 1,0 1,计算得到R g,再利用R o 通过(2 3)式和(2 4)式即可估计出a o,w o。(4)循环上述步骤，直至收敛求得a m,0 m。4.3 M E D L L 与W R E L A X 算法多径信号估计误差对比采用两路信号模型，设真实信号的幅值为1,多径信号的幅值为0.5,多径延时为0.1 c h i p s。在有无噪声的情况下，M E D L L 和W R E L A X 算法得到的伪距跟踪误差如图9 所示。图9(a)可以看出，在无噪声情况下，随着多径延迟的不断增加，W R E L A X 始终比M E D L L 估计误差小。在有噪声情况下，如图9(b)-9(e)所示，信噪比分别为-5 d B、-1 1 d B、-1 6 d B、-2 0 d B 情况下，M E D L L 和W R E L A X 的伪码跟踪误差变化情况，可以看出M E D L L 的伪码跟踪误差比W R E L A X 大，在信噪比为-2 0 d B 情况下，M E D L L 的误差甚至达到了0.2 c h i p s 左右，1 1通信与信息技术2 0 2 3 年第4 期(总第2 6 4 期)无噪声情况下N E D L L 和K R E L A X 误差对比0.5-M E D L L-W R E L A X0.4伪码跟踪误差(chips)0.3-0.2.中a.00.00.51.0多径延迟(c h i p s)(a)无噪声情况下的误差对比1.5噪声情况下M E D L L 和W R E L A X 误差对比0.5-0.4伪码跟踪误差(chips)0.3-0.20.1*M E D L LW R E L A X0.0+0.00,51.0多径强迟(c h i p s)1 s(b)信噪比为-5 d B 时的误差对比0.5-噪声情况下M E D L L 和W R E L A X 误差对比0 4伪码银踪误差(chips)0.30.20.1-0.0-M E D L LW R E L A X0.00.51.0多径延迟(c h i p s)L s噪声情况下M E D L L 和V R E L A X 误差对比0.5-M E D L-W R E L A X0.4-伪码跟踪误差(chips)0.30.2+0.10.0-0.00.51.0多径延迟(c h i p s)1.5(c)信噪比为-1 1 d B 时的误差对比0.5噪声情况下M E D L L 和W R E L A X 误差对比-M E M L-W R E L A X0.4-伪码跟踪误差(chips)0.3-0.20.1-6.00.00.51.01.5多径延迟(c h i p s)(e)信噪比为-2 0 d B 时的误差对比图9 M E D L L 与W R E L A X 算法的估计误差对比5 结论以卫星导航系统为例，分析了多径信号的来源，构造了多路径信号的模型，得到以下结论：(1)利用E M L 鉴相器和二象限反正切鉴相函数分析多径信号对卫星接收机的影响，对于不同的多径信号，对卫星接收机跟踪环造成的影响程度不同。(2)利用M E D L L 和W R E L A X 算法对多径信号的伪码跟踪误差进行仿真分析，在信噪比大于-2 0 d B 时，W R E L A X(d)信噪比为-1 6 d B 时的误差对比算法对多径信号伪码跟踪误差的估计精度优于M E D L L 算法。参考文献 1 J i a Q i o n g q i o n g,W u R e n b i a o,W a n g W e n y i,e-t a l.M u l t i p a t hi n t e r f e r e n c e m i t i g a t i o n i n G N S S v i a W R E L A X J .G P S S o l u t i o n s,2 0 1 7,2 1(2):4 8 7-4 9 8.2 张梦尧，陈熙源.改进的M E D L L 技术在多径误差估计中的应用 J .导航定位学报，2 0 1 9,7(0 3):9 6-1 0 1.3 方玲.G P S 接收机多径抑制算法研究 D .南昌：华东交通大学，2 0 1 4.4 王志远.导航系统中的多径误差抑制算法研究 D .太原：太原理工大学，2 0 1 7.5 周博海，郑建生，陈鲤文，朱玉建.G P S 接收机在高动态下的抗干扰技术实现 J .科学技术与工程，2 0 1 5,1 5(2 1):1 3 6-1 4 1.6 刘雪萌.导航系统中基于参数估计的多径抑制算法研究 D .哈尔滨：哈尔滨工程大学，2 0 2 1.7 张晓婷.G P S 系统多径抑制技术方法研究与优化 D .上海：上海交通大学，2 0 1 6.8 G a o Y,L i u F,L o n g T.A n a l y s i s o f m u l i p a t h p a r a m e t e r e s t i m a t i o n a c c u r a c y i nM E D L L a l g-o r i t h m C .C h i n a S a t e l l i t e N a v i g a t i o n C o n f e r e-n c e(C S N C)2 0 1 3P r o c e e d i n g s.S p r i n g e r B e r l i n H e i d e l b e r g,2 0 1 3:5 9 7-6 0 6.(下转第1 6 页)1 2通信与信息技术2 0 2 3 年第4 期(总第2 6 4 期)1 0 黄国兴，王东，陈光生.医院数据中心虚拟化架构探讨 J .医学信息学杂志，2 0 1 9,4 0(0 1):2 6-2 9.1 1 宋阳，郝泽余，刘坤尧.医院跨平台容灾数据库的设计与实现 J .中国数字医学，2 0 1 9,1 4(0 2):8 8-9 0.1 2 王建英，黄士琴.基于C D P 技术的医院信息系统数据备份与保护 J .数字技术与应用，2 0 1 9,3 7(0 8):1 6 7-1 7 0.作者简介宋阳(1 9 8 9-),男，硕士，工程师，主要研究方向：无线通信、医院信息化建设。T h e r e s e a r c h o f d i s a s t e r r e c o v e r y s y s t e m f o r t h e h o s p i t a l i n f o r m a t i o n s y s t e mS O N G Y a n g,Z O U J i a n,Z H O U X i nT h e A f f i l i a t e d H o s p i t a l o f S o u t h W e s t M e d i c a l U n i v e r s i t y;L u Z h o u,6 4 6 0 0 0,C h i n aA b s t r a c t:I n o r d e r t o e n s u r e h o s p i t a l d a t a s e c u r i t y a n d b u s i n e s s c o n t i n u i t y a n d s o l v e t h e d i s a s t e r r e c o v e r y p r o b l e m o f h o s p i t a l,t h i s p a p e r s t u d i e s t h e d i s a s t e r r e c o v e r y s y s t e m c o n s t r u c t i o n s c h e m e o f h o s p i t a l i n f o r m a t i o n s y s t e m.B y c o l l e c t i n g f a u l t c a s e s o fh o s p i t a l i n f o r m a t i o n s y s t e m,s u m m a r i z i n g t h e f a u l t t y p e s a n d c a u s e s,s t u d y i n g d i s a s t e r r e c o v e r y t e c h n o l o g y a n d m a n a g e m e n t s y s t e m,f o r m i n g a d i s a s t e r r e c o v e r y s y s t e m a n d v e r i f y i n g w h e t h e r t h e e x p e c t e d g o a l i s a c h i e v e d.A c c o r d i n g t o t h i s s y s t e m,t h e d i s a s t e rr e c o v e r y s y s t e m o f t h e i n f o r m a t i o n s y s t e m o f t h e A f f i l i a t e d H o s p i t a l o f S o u t h w e s t M e d i c a l U n i v e r s i t y w a s d e s i g n e d a n d c o m p l e t e d,a n d t h e b u s i n e s s w a s q u i c k l y r e c o v e r e d i n t h e r e a l f a u l t,t h e d a t a w a s n o t l o s t,a n d t h e e x p e c t e d d i s a s t e r r e c o v e r y e f f e c t w a s a c h i e v e d.T h e r e s u l t s s h o w t h a t i t i s n e c e s s a r y t o b u i l d a n i n f o r m a t i o n s y s t e m d i s a s t e r r e c o v e r y s y s t e m i n h o s p i t a l,w h i c h c a n e f f e c t i v e l y d e a lw i t h i n f o r m a t i o n s y s t e m f a i l u r e s.K e y w o r d s:D i s a s t e r r e c o v e r y s y s t e m,I n f o r m a t i o n s y s t e m f a i l u r e s,D a t a s e c u r i t y,S e r v i c e c o n t i n u i t y(上接第1 2 页)1 9 1 吴仁彪，王文益，卢丹，王璐，贾琼琼.卫星导航自适应抗干扰技术 M .北京：科学出版社，2 0 1 5:1 4 7-1 6 1.1 0 贾琼琼，吴仁彪，王文益，卢丹，王璐.基于W R E L A X 的小运算量G N S S多径干扰抑制算法 J .数据采集与处理，2 0 1 8,3 3(0 3):4 6 9-4 7 6.作者简介朱奎宝(1 9 7 2 一),男，博士，教授。主要研究方向：主要从事系统协同导航技术研究。E-m a i l:z h u k u i b a o 1 6 3.c o m郭广源(1 9 9 6 一),男，主要研究方向：G P S 多径信号的抑制研究。E-m a i:8 1 2 2 5 8 0 1 4 a q g.c o mI n f l u e n c e o f m u l t i p a t h s i g n a l o n s a t e l l i t e r e c e i v e r a n d r e s e a r c h o n t i m e d o m a i ns u p p r e s s i o n a l g o r i t h mZ H U K u i b a o,G U O G u a n g y u a n,D U y u n,W E N Z i q i n g,D E N G C h e n g b i n,K A N G H a o n a nS c h o o l o f E l e c t r i c a l E n g i n e e r i n g,H e b e i U n i v e r s i t y o f S c i e n c e a n d T e c h n o l o g y;S h i j i a z h u a n g,0 5 0 0 1 8,C h i n aA b s t r a c t:A i m i n g a t t h e m u l t i p a t h e r r o r p r o b l e m o f s a t e l l i t e n a v i g a t i o n s y s t e m,a m u l t i p a t h e f f e c t m o d e l i s e s t a b l i s h e d,a n d t h ei n f l u e n c e o f c a r r i e r t r a c k i n g l o o p a n d p s e u d o-r a n d o m c o d e t r a c k i n g l o o p i n m u l t i p a t h e n v i r o n m e n t i s a n a l y z e d b y u s i n g t w o q u a d r a n ta r c t a n g e n t p h a s e d i s c r i m i n a t o r a n d E M L p h a s e d i s c r i m i n a t o r.I n t h e a s p e c t o f m u l t i p a t h s u p p r e s s i o n a l g o r i t h m,M E D L L a n dW R E L A X a l g o r i t h m a r e u s e d t o s u p p r e s s m u l t i p a t h s i g n a l s a n d s i m u l a t e d,a n d t h e p e r f o r m a n c e o f t h e t w o a l g o r i t h m s f o r m u l t i p a t hs i g n a l s u p p r e s s i o n i s c o m p a r e d a n d a n a l y z e d.T h e e x p e r i m e n t s h o w s t h a t d i f f e r e n t m u l t i p a t h s i g n a l s h a v e d i f f e r e n t e f f e c t s o n t h er e c e i v e r a n d h a v e d i f f e r e n t d e g r e e s o f i n f l u e n c e o n t h e r e c e i v e r t o s o l v e t h e r e a l s i g n a l.I n t e r m s o f m u l t i p a t h s u p p r e s s i o n,W R E L A Xa l g o r i t h m i s b e t t e r t h a n M E D L L a l g o r i t h m i n s u p p r e s s i n g m u l t i p a t h s i g n a l s,a n d t h e e s t i m a t i o n e r r o r i s s m a l l e r.K e y w o r d s:S a t e l l i t e n a v i g a t i o n s y s t e m,M u l t i p a t h s i g n a l,R e c e i v e r,M E D L L,W R E L A X1 6