GNSS系统中一种低复杂度的载噪比估计算法

载噪比(carrier to noise ratio，C/N0)是描述接收机性能的一个很重要的输出参数，并且可以衡量跟踪环路的跟踪状态和各个通道的跟踪策略。载噪比不仅是接收机输出给用户的一个测量值，同时也是接收机信号处理中一个重要的控制量，用于定位解算中的加权最小二乘法、设置信号检测门限以抑制互相关干扰，以及利用载噪比来抑制多路径干扰[1]。信号捕获门限与信噪比(signal to noise ratio, SNR)有着直接关系，信号跟踪环路的锁定检测和接收机性能的预估也都与载噪比的测定相关[2]。信噪比是信号功率与噪声功率之比，而载噪比是信噪比的一种归一化表现形式[3]。输出的载噪比估计值可以反映该通道跟踪信号的强度，通常情况下，信号的载噪比越高意味着信号强度越高，则将会得到更精确的伪距和载波相位测量精度以及更好的定位精度。

接收机的载噪比估计依赖于信号功率估计与噪声功率估计，大部分估计算法都是通过I、Q两路累加值估计得到的[4-5]。比较常见的载噪比估计方法有宽窄带功率比值法(narrow to wide power ratio method, NWPRM)[6-8]、方差求和法(variance summing method, VSM)[9-11]、I支路功率方差比方法[12]和平方信噪比方差方法[13]。其中，应用最为广泛的是宽窄带功率比值法，由于其在估计过程中采用了20 ms的相干积分，因此需要先进行比特同步后再估计载噪比值，并且不适用于北斗系统以2 ms调制一个导航电文的GEO卫星；方差求和法在虽然能够应用于GEO卫星，在一些多信号体制接收机中有所应用，但由于其在估计过程中用到了求方差运算，而方差的样本量至少要20 ms以上的数据估计得到的载噪比才具有一定的可信度，因此需要存满20 ms的数据后才能进行运算，这就占用了存储器资源，在实际DSP实现的时候是一种对存储器资源的浪费。

为了解决上述问题，本文提出了一种载噪比估计算法，该方法利用I支路二阶矩和I、Q两支路四阶矩估计的分布特性，推导出了载噪比估计的具体表达式。本文从估计偏差、实现复杂度、适应性和输出时延四个方面对NWPRM、VSM和本文提出的算法进行了分析比较。

1 信号模型

当接收机跟踪某个信号时，由该信号通道上I支路和Q支路输出的相干积分值IP(n)和QP(n)的表达式如下[14]

 

(1)

 

(2)

式中：D(n)是值为±1的数据比特电平值，τP为即时复制C/A码与接收C/A码之间的相位差异，R(·)代表最大值为1的C/A码自相关函数；且信噪比SNR=A2/2σ2,载噪比C/N0与信噪比SNR的关系为C/N0=SNR/2T;nI和nQ均为均值为零，方差为σ2/2的随机正态噪声。

假设码对齐假设码对齐且码相关过程中卫星电文数据D(t)=1，则上面两式可简化为

 

(3)

 

(4)

上面两式还可以写成如下一个关于载噪比C/N0的函数：

 

(5)

 

(6)

进一步对以上两式单位化：

 

(7)

 

(8)

此时，nI和nQ变为均值为零，方差为1/2的随机正态噪声。

式中：ei(k)为第k位专家给出的第i项指标的重要性序数；ej(k)为第k位专家给出的第j项指标的重要性序数。于是得到13个元素只有0和1的设计矩阵Aij(k)。

另外由式(3)、(4)可知，若φe=0，则有

 

(9)

式中：nI～N(0,σ2/2)，nQ～N(0,σ2/2)。

在推导出I、Q两路每1 ms累加值的概率分布之后，就可以对这些载噪比估计方法进行推导。

2 典型载噪比估计算法

2.1 宽窄带功率比值法

宽窄带功率比值法的基本原理是根据信号加噪声的功率在不同噪声带宽上的差异，得到宽带功率和窄带功率并利用二者的比值来计算载噪比C/N0。

带宽为1/T的宽带功率PWB(k)被定义为

 

(10)

带宽为1/MT的窄带功率PNB(k)则被定义为

方案一和方案二在本质上是一样的，方案二的实质就是方案一中式(31)等式两边分子分母同时调换解出的结果。下面从理论和测试两个方面来验证其一致性。

 

(11)

将窄带功率PNB(k)用相应的宽带功率PWB(k)来进行单位化得到比率值PNW(k)为并对其做平均，得到

 

(12)

由此可求得载噪比C/N0为

 

(13)

2.2 方差求和法

因此，当nI～N(0,σ2/2)，令则有Y～N(0,1),可得Y2～χ2(1)。

式(33)得到的值为一个假值。

 

(14)

则参量Z的期望和方差分别为

 

(15)

 

(16)

因此载噪比估计值C/N0为

 

(17)

3 一种新型载噪比估计算法

3.1 算法的理论推导

该方法是基于I路信号累加值的二阶矩和I、Q两路信号累加值的四阶矩估计，具体的推导如下。

学生顶岗实习管理系统主要面向高职院校学生的顶岗实习管理，需要根据教师、学生的职责进行任务的分配，对学生实习过程中涉及的学生基本信息、周志、实习总结、教师分阶段评价等信息进行详细的记录和梳理，从而有效的提高学生顶岗实习过程管理的质量，做到规范性。

定义I路信号累加值的二阶矩：

证法5 (构造辅助圆法)如图7，以点E为圆心，EB为半径作⊙E，作EF⊥OC，垂足为F，连接OB，则四边形OAEF是矩形．因而AE=OF．由EB′=EB，知⊙E过点B′．由矩形OABC的轴对称性，有∠BAC=∠BOB′．

 

(18)

由概率论知识可知期望和方差的关系如下[15]

E(X2)=E2(X)+D(X)

(19)

由式(9)、(19)可得到MI2的均值为

 

(20)

定义I、Q两路信号累加值的四阶矩：

 

(21)

将式(9)代入式(21)可得

表1详细描述了四个班级(较高中等听力水平和较低中等听力水平)学生在主测1和主测2中的平均分和标准差表现。

 

(22)

进一步展开可得

 

(23)

令则由式(9)、(19)可得

 

(24)

由概率论知识,可知若X～N(0,1)则有其均值为n方差为2n。

方差求和法的基本原理是通过即时支路输出的I、Q的累加值的期望和方差来分别计算有用信号功率和噪声功率,从而得到载噪比估计值。

由χ2分布性质，可知

D(Y2)=2即

(25)

可得到

 

(26)

同理可得

 

(27)

3.2.2 数据测试验证

 

(28)

将式(28)代入式(19)可得

 

(29)

MI2的均值的平方E(MI4)为

相关研究发现有控制排水设施的排水沟排放的氮浓度比无控制排水设施的要低[12]，为此，在农田铺设地下排水管道，并安装控制排水阀，与其集水的渗滤排水明沟及沟道出口处水量、水位控制构筑物形成控制排水系统。实现精确调控排水沟渠中的水位，增加田间入渗量，控制农田径流的排出过程。

 

(30)

1)方案一：

令得到

 

(31)

多Agent仿真方法是团队仿真的通用技术，相关研究和案例较多，如王剑等[6]基于贝叶斯决策网络对团队协作进行了仿真研究，提出双层团队协商模型以实现个体利益需求的均衡；王君华等[7]从合作博弈角度研究了创新研发团队的群体和个体收益，并给出了合作的收益函数。

（3）人民币汇率预期与证券投资跨境资金流入之间的脉冲响应分析。对人民币汇率预期、境内外利差与证券投资跨境资金流入之间建立VAR模型，确定VAR模型滞后阶数为5阶。根据所建立的VAR（5）可以得到脉冲响应函数图（如图3所示）。

 

(32)

 

(33)

接收信号中信号和噪声总功率为

因此，可以得到载噪比C/N0：

 

(34)

2)方案二：

令得到

 

(35)

解式(35)得到

 

(36)

 

(37)

式(36)得到的值为一个假值。

因此，可以得到载噪比C/N0：

EVA考核中资本成本率的设置直接影响EVA的高低，影响指标的准确性。 2016年12月12日国务院国资委印发的《考核办法》(国资委第33号)突出经济增加值考核，引导企业资本投向更加合理，注重企业分类及差异化因素，根据功能不同、资本结构、风险程度不同的企业，提出差异化资本回报的要求。

解式(31)得到

 

 

(38)

3.2 方案一致性验证

本次调查共回收问卷201份，其中有效问卷175份，约有26名受访者表示根本不了解运费险。根据本文所要研究的问题和数据结果，在此主要从以下三方面进行分析。

3.2.1 理论验证

将方案一中的代入式(34)中得到

 

(39)

将方案二中的代入式(38)得到

 

(40)

对比式(39)、(40)可以发现两式是相同的。因此，从理论上来讲两种方案结果一致得到验证。

鸦片战争前后，为了解决严峻的民族危机，湖湘学派的代表人物魏源大举倡导 “变法改革的经世致用”的思想，提出“师夷长技”的观点。魏源认为这里的“致”有着施行、取得的双重涵义，这里的“用”也有着日用、功用的双重涵义；因而“经世致用”就是把经义施行于日用实际，取得经国济民之功用。这无疑包含着知行问题，因为经义属于“知”的范畴，而施行之义的“致”则与“行”相通。为此魏源揭开了“知行之辩”在近代的序幕。

因此M′的方差为

仿真采用真实GPS中频数据，中心频率为1 572.42 MHz，数字基带中频为4.092 MHz，采样频率为16.368 MHz，取1 ms相干积分，对1 000次估计结果取平均。无多普勒频率残差的情况下，对两种方案的载噪比估计进行仿真，得到的仿真结果如表1。

 

表1 两种方案实测数据验证Table 1 The validation of two schemes

  

信号功率-138 dB·m-135 dB·m-134 dB·m-133 dB·m方案一-138.067 9-135.256 6-134.361 8-133.280 3方案二-138.067 9-135.256 6-134.361 8-133.280 3

由表1可知，两种方案得到的载噪比相同。但是方案一中计算载噪比的表达式更为简单，计算量相对来讲更小。因此在后面的测试比较中采用方案一来与其他几种算法进行测试比较性能。并将新方法命名为I支路矩方法(I branch moment method, IBMM)。

2.3.1 防治CINV用药存在问题 DXM的问题集中在给药剂量不合理和疗程不足。应用5-HT3RA存在给药疗程不合理，占总例数的71.17%。104例（93.69%）存在H2RA/PPI无指征用药的情况。见表5。

4 性能分析

几种算法在载噪比估计偏差、算法实现复杂度、输出时延、适应性等方面存在差异。下面将分别从这四个方面进行仿真分析。

4.1 估计偏差比较

4.1.1 仿真数据测试分析

理论仿真采用Matlab生成GPS 19号卫星数据，仿真参数：采样频率为16.368 MHz，系统中频为4.092 MHz。为了更精准的对比分析上述三种算法，本文采用长度相同的数据来进行性能分析，具体仿真方案为：宽窄带法采用20 ms相干积分，对50次估计结果取平均；其他两种算法均是取1 ms相干积分，对1 000次估计结果取平均。无多普勒频率残差的情况下，对信号功率为-138～-120 dB·m的数据三种算法分别进行载噪比估计仿真，得到的载噪比估计偏差如下：

由图1可知，在弱信号条件下，方差求和法的估计偏差为0.5 dB, IBMM偏差为0.6 dB, NWPRM偏差为0.65 dB,说明IBMM和NWPRM弱信号应用时性能较差；在强信号条件下，四种算法偏差均在0.5 dB以内，算法性能相当。总的来说，本文提出的算法性能介于宽窄带和方差求和法之间。

  

图1 仿真数据载噪比估计偏差Fig.1 The C/N0 estimation deviation of simulated data

4.1.2 实测数据测试分析

仿真采用真实GPS中频数据，中心频率为1 572.42 MHz，数字基带中频为4.092 MHz，采样频率为16.368 MHz，其他参数和理论仿真一致。对各算法所能达到的载噪比估计偏差进行仿真，得到的仿真结果如图2。

由图2可知，几种算法的估计偏差相差不大，都在0.5 dB以内。其中，在信号较强的三组数据时，I支路矩估计法的偏差介于宽窄带和矩估计法之间。上述实测结果和理论分析一致，说明本文提出算法的正确性。总的来说，三种算法的估计性能相差不大。

  

图2 实测数据载噪比估计偏差Fig.2 The C/N0 estimation deviation of measured data

4.2 算法实现复杂度比较

4.2.1 运算量比较

对三种算法的运算量进行分析，为了精确的比较三种算法的运算量，采用长度相同的数据来进行性能分析，具体方案为：宽窄带法采用20 ms相干积分，对50次估计结果取平均；其他两种算法均是取1 ms相干积分，对1 000次估计结果取平均。具体结果如表2。

 

表2 三种算法运算量比较Table 2 The calculation comparison of three algorithm

  

运算形式AddMulsqrt10log10NWPRM3 9512 15401VSM4 0013 00721IBMM3 0023 00811

由表2可知，三种算法中，NWPRM和IBMM算法的运算量相当，VSM与IBMM相比不仅多了开方运算而且加法次数接近多了1 000次，因此本文提出的算法运算量较为适中。

4.2.2 占用存储器资源比较

在实际应用到DSP中时，由于VSM在求方差时会用到均值，而均值是在存满1 000次结果得到的，因此VSM需要先将1 000个数据存储后再进行运算。其中，存储一组数据需要的存储器大小为1 000×8 byte，这会对存储器资源造成巨大的浪费。但NWPRM和IBMM在DSP中应用时可以做到实时的运算，存储一组数据需要的存储器大小仅为2×8 byte，并不会占用太大的存储器空间。因此，本文提出的算法对存储器资源消耗较少。

4.3 输出时延比较

宽窄带法计算载噪比时是每20 ms输出一个结果。方差求和法计算载噪比时含有求方差运算，至少是2 ms输出一个结果。通常情况下，要达到比较准确的估计输出，方差求和法需要20 ms的数据来准确估计样本的方差。而本文提出的方法则使用比较灵活，可以1 ms输出一个结果，也可以设置为任意时间长度输出。

4.4 适应性比较

目前GNSS包含了美国的GPS、俄罗斯的GLONASS、欧盟的Galileo、中国的Compass(北斗)。其中北斗的GEO卫星2 ms调制一个导航电文，宽窄带法采用20 ms相干积分，因此它不能直接应用于BD系统。方差求和法虽然能2 ms输出一个结果，但要达到比较准确的载噪比估计值输出时，方差求和法则需要20 ms的数据来准确估计样本的方差。而本文提出的方法可以任意ms输出一个结果，提升了适应性，可广泛用于北斗等GNSS系统。

5 结论

1) 仿真数据和实测数据都表明，本文提出的IBMM算法的估计偏差介于NWPRM和VSM之间，三种算法的估计偏差都维持在0.5 dB以内，总的来说，三种算法的估计性能相差不大。

2) 从实现复杂度方面来看，三种算法中，VSM的运算量最大且实现复杂度较高，存储数据过程中会对存储器资源造成极大的浪费；NWPRM和本文提出的IBMM运算量相当，且都能做到实时的运算，不会占用太大的存储器空间。

3) 输出时延和适应性比较结论：NWPRM和VSM存在20 ms的输出时延，而本文提出的算法则可以做到任意毫秒输出一个结果，而NWPRM不能应用于北斗系统，要能够准确估计GEO卫星的载噪比值时，VSM就不适用了，而本文提出的IBMM算法可广泛应用于北斗等GNSS系统，在这方面具有明显的优势。

参考文献：

[1] 何文涛, 徐建华, 叶甜春. GPS弱信号的自适应载噪比估计算法[J]. 电子技术应用, 2010, 36(6): 111-114.

HE Wentao, XU Jianhua, YE Tianchun. Adaptive C/N estimation method of GPS weak signal[J]. Application of electronic technique, 2010, 36(6): 111-114.

[2] 袁建国, 欧松林. GPS接收机并行捕获算法的研究[J]. 重庆邮电大学学报(自然科学版), 2016, 25(4): 470-474.

YUAN Jianguo, OU Songlin. Study on parallel acquisition algorithm for GPS receiver[J]. Journal of Chongqing university of posts and telecommunications (natural science edition), 2016, 25(4): 470-474.

[3] SHARAWI M S, AKOS D M, ALOI D N. GPS C/N0 estimation in the presence of interference and limited quantization levels[J]. IEEE transactions on aerospace and electronic systems, 2007, 43(1): 227-238.

[4] BRAASCH M S, VAN DIERENDONCK A J. GPS receiver architectures and measurements[J]. Proceedings of the IEEE, 1999, 87(1): 48-64.

[5] BORIO D, GIOIA C, BALDINI G. Asynchronous pseudolite navigation using C/N0 measurements[J]. Journal of navigation, 2016, 69(3): 639-658.

[6] 傅金琳, 赵子阳, 李醒飞. 一种改进的功率比值法载噪比估算算法[J]. 中国惯性技术学报, 2015, 23(5): 642-647.

FU Jinlin, ZHAO Ziyang, LI Xingfei. Improved power ratio method of carrier-to-noise ratio estimation algorithm[J]. Journal of Chinese inertial technology, 2015, 23(5): 642-647.

[7] GROVES P D. GPS signal to noise measurement in weak signal and high interference environments[C]//Proceedings of the 18th International Technical Meeting of the Satellite Division of the Institute of Navigation. Long Beach, CA, 2005: 644-645.

[8] BHUIYAN M Z H, SÖDERHOLM S, THOMBRE S, et al. Performance evaluation of carrier-to-noise density ratio estimation techniques for BeiDou Bl signal[C]//IEEE Ubiquitous Positioning Indoor Navigation and Location Based Service. Corpus Christ, TX, USA, 2015: 19-25.

[9] LI Li, CHEONG J W, WU Jinghui, et al. Improvement to multi-resolution collective detection in GNSS receivers[J]. Journal of navigation, 2014, 67(2): 277-293.

[10] PSIAKI M L, AKOS D M, THOR J. A comparison of “Direct RF Sampling” and “Downconvert & sampling” GNSS receiver architectures[C]//Proceedings of the 16th International Technical Meeting of the Satellite Division of the Institute of Navigation. Portland, OR, 2003: 1942-1952.

[11] 蒋君伟, 段晓辉, 林阳. GPS信号载噪比估计算法的研究[J]. 北京大学学报(自然科学版), 2009, 45(3): 409-414.

JIANG Junwei, DUAN Xiaohui, LIN Yang. The research of estimation method for GPS C/N0[J]. Acta scientiarum naturalium universitatis pekinensis, 2009, 45(3): 409-414.

[12] 巴晓辉, 刘海洋, 郑睿, 等. 一种有效的GNSS接收机载噪比估计方法[J]. 武汉大学学报(信息科学版), 2011, 36(4): 457-460.

BA Xiaohui, LIU Haiyang, ZHENG Rui, et al. An effective carrier-to-noise ratio estimation method for GNSS receiver[J]. Geomatics and information science of Wuhan university, 2011, 36(4): 457-460.

[13] FALLETTI E, PINI M, PRESTI L L. Low complexity carrier-to-noise ratio estimators for GNSS digital receivers[J]. IEEE transactions on aerospace and electronic systems, 2011, 47(1): 420-437.

[14] 谢钢. GPS原理与接收机设计[M]. 北京: 电子工业出版社, 2017: 306-364.

XIE Gang. Principles of GPS and receiver design[M]. Beijing: Publishing House of Electronics Industry, 2017: 306-364.

[15] 盛骤, 谢式千, 潘承毅. 概率论与数理统计[M]. 4版. 北京: 高等教育出版社, 2010: 46-103.

SHENG Zhou, XIE Shiqian, PAN Chengyi. Probability theory and mathematical statistics[M]. 4th ed. Beijing: Higher Education Press, 2010: 46-103.