导航卫星卷积码序列 MAP 译码算法分析与改进
胡 毅1 宋茂忠1,2
孟 斌2
1 南京航空航天大学电子信息工程学院 ,江 苏 南京 ,210016 2 中国航天科技集团公司卫星应用研究院 ,北 京 ,100086
摘 要 :在 对导航卫星卷积码序列应用 MAP 算法译码时 ,计 算机会产生数值稳定性问题 ,输 出错误的译码误 码率 ,进 而对定位结果的有效性和定位精度产生重要影响。 针对该问题进行了具体分析 ,并 对算法进行了 改 进 ,给 出了一种基于最大状态度量的归一化 MAP 译码 算 法 ,在 保 持 原 算 法 最 优 性 能 的 同 时 很 好 地 解 决 了 数 值稳定性问题。通过仿真进一步验证了这种改进算法的有效性。 关键词 :导 航卫星卷积码 ;MAP 译码算法 ;数 值稳定性 ;最 大状态度量 中图法分类号 :P228〃41
文献标志码 :A
卷积码是现代导航卫星广泛使用的一种差错
控制编码[1-4
],如在 GPS卫星现 代化中,L2C 信号 导航电文就使用限制长度为 7、码 率为 1/2 的 (2,
1,6
)前向纠错卷积 码[3-5
]。 其常规译码算法采用 的是 维 特 比 (Viterbi
)译 码,而 相 对 于 维 特 比 译 码,最 大 后 验 概 率 MAP(maximum aposteriori
probability
)译码则是一种性能最 优译码[6-7]。 但 另一 方 面,应 用 中 发 现,在对 较 长卷积码序列如
L2C 导航电文数据块(600bit)应 用 MAP 算法译
码时,会产生数值稳定性问题 ,进而对定位结果的
有效性和定位精度产生重要影响。对此给出了一 种改进的归 一 化 MAP 译 码 算 法,很 好 地 解 决 了 该问题。
图 1 (2,1,3)
卷 积码编码器 (2,1,3)ConvolutionalCodeEncoder Fig〃
1 [7]: uk 的 MAP 算法估计^uk 为
∑ αk-1(p)γk(p,q)βk(q)燄 熿 L(uk)=log
(p,q)∈S1 () 1 ∑ αk-1(p)γk(p,q)βk(q) 燀(p,q)∈S
0
^uk =
燅 1,L(uk)≥0
{0,L(uk)<0
(2
) MA
P 译码算法及其不足 1 式中,p、q 分别为k-1 和 k 时 刻 的 编 码 器 状 态, αk-1(p)、β
k(q)和γk(p,q)分 别为前向和后向状态 度量及状态转移支路度量 ,其满足递推关系:
αk(q)=
为便于分析问题 ,给 出 1 位输入 2 位输出,3
位移位寄存器的(2,1,3)前向纠错卷积码 ,如图 ∑α
k-1
(p)γk(p,q)
(3) 1 所 示。 图 中 uk 为 k 时 刻 输 入 信 息 位,uk ∈
(p,q)∈S
{
0,1},ck=(ck (1) (2)
,ck )为 编码器对应的输出码字, k(
p)= γ
k+1
(p,q)β
k+1(q) (4) (5
) ∑ (p,q
) β c
(j) ∈ {0,1},j=1,2。 假设编码器 输 出 序 列 为c= (c1,c2,…,cN ), ∈S k γk(p,q)= Pr {uk }p{yk|ck }
其中1≤k≤N。上面的递推初始条件一般为 : 对应的译码器接收序列为 y=(y1,y2,…,yN ),yk
1,p =0 1,q =0 =(
y
(1) (2) )。令编码器状态集为S,根据输入uk k ,yk {
{
0,q ≠0 0( ) ,βN () () 6 α p = q = 0,p ≠ 0 的值可以 把 S 划 分 成S0 和S1 两
个 子 集,分 别 与 uk=0和uk=1 相对应。 这样k 时刻输入信息位
应用中发现,在对较长卷积码序列进行译 码
时,经 过多次状态转移后 ,αk (p)、β
k (q)的 值就会 收稿日期 :2013-01-15
项目来源 :中国航天科技集团公司卫星应用研究院创新基金资助项目(20121512)
。 第一作者 :胡 毅 ,博 士生 ,现主要从事导航卫星信号处理方面的研究 。E-mail:hy
starcom@163〃com
武 汉 大 学 学 报 · 信 息 科 学 版
2014年 3月
2
68 变得越来越小,最后产生 溢出而变为 0,从 而导致 式(1
)无法给出正确的判决结果 ,进而给出错误的 译码误码率。
表单大师升级容量右边的数值会有所不同。
算法的改进及其应用性能分析
3 算法数值不稳定性分析
2 算法的改进
通过归一化,在不改变判决过程的情况下 ,可
使译码过程数值稳定 ,避免溢出的发生[
7
]。 3〃1 由于βk(q)的分析过 程与αk (p)类 似,因 此这
里只对αk(p)
的情况进行分析 。 对于式 (3),令k 时的状态转移为 (p,q)k,同 时令γk=γk(p,q)(1≤k≤
N),则有: 归一 化 系 数 可 通 过
q
α′k(q)= p
β
′k(p) ∑ ∑ =1来获得[7]
,也 可 由α ( )和β (
)的 最 大 k-1 p k q 来 确 定。 另 外,还 可 选 用 中 值 数 α = ex p
mid αk(pk)=γ1γ2γ3 … [
γ4γ5…γk-1γk ] ∑ ∑ [8]
g max)/2]
来 得 到 。 考 虑 到 计 算 的 [(logα min+loα (p,q) ( ,) 4
∈S
pqk
∈
S 简单性及减少舍入误差的影响 ,文中取状态度量
的最大值来得到归一化系数。令
(7
) 忽略前3项有限值,重点考虑后面的求和项 。令
αk-1 = max{αk-1(p)},βk = max{k(q)} β …
[
γ4γ5…γk-1γk ] (8) Sk-3 =
∑
∑ p∈S q∈S 对于任一k(1≤k≤N)
下的所有可能状态,总 能 保证αk-1>0,βk>0。取 (p,q) ( ,) 4
∈S
pqk
∈
S 并取γt=max{γ4,γ5,…,γk },其 中 4≤t≤k
,同 时 注意到(8)式每个求和项都有 2支,分别与0、1输 入相对应,这样可得:
k-1(p) k-1( )/ k-1,βk(q)=βk(q)/βk α′ α ′ =α p 则根据式(1
)有: k-
3
(9) Sk-3 ≤ (2γt )
k-1 k() α p q (α ) γk(p,q)β (β
) 熿
燄 ∑ 假设信道为无衰落零均值理想高斯白噪声信
道,即yt=ct+nt。
由 于 先 验 概 率 未 知,可 假 设 各 输入信息位等先验概率[8]
,从而有:
(p,q)∈S1java和golang优缺点
-1 k ( ) L uk =log
= k-1 β
k() α p q (
)γk(p,q)( )
∑ 燀(p,q)∈S0
燅 β
k k-1 γt =exp(- ‖yt -ct‖2
/N0 )
(10
) 熿 ∑ α′k-1(p)γk(p,q)β′k(q)燄 式中,N0 为噪声功率谱密度。将nt 代入上式得: 2
(p,q)∈S1 (14
) log
( ) ( ,) () ∑
α′k-1 p γk pqβ′k q 燅 n(1) 烄 熿(t / 槡N0/2) 燀 燄烌 γt =exp - 1 (p,q)∈S (11) 0
烅 烍
2 2 由于上面的处理过程是恒等变换 ,因 此不会
改变算法的最优性能。 据此可得改进后的 MAP 译码算法的主要步骤为:
n
(2) + (t / 槡N0/2)燅烎 烆 燀 由于n(1) (2) t 、nt 相互独立且 ~N(0, n
(i)/ (t 槡N0 2) / 1)分布,i=1,2,这样有: 初始化α0(p)和βN (q),并取k=
1; 步骤1 步骤 2 (1
) 2
(2
) 2
=
+
]~χ
2
(2) 2
χ
t
( )及 式 (3) ( )计 算 出 时 的γk ( , )、 10 ~ 4 k p q αk-1
(12
) 将式(12)代入式(11),当 2γt<1,即 2exp(-χt (p)和βk
();2/ 分别计算出 p 状态下的最大状态度
步骤3 2)<1时,
有: 量αk-1和q 状态下的最大状态度量β2
(13) >2ln2}≈0〃51,即 >2ln2 χt
将 得 到 的 αk-1 和 βk 代 入 (14)式 计 由χ (2)性质可知,Pr {χ
步骤4 2 2 算出α′k-1(p)、β′k(q)及 L(uk ),然 后 再 按 照 式 (2) 进行译码判决;
步骤 5 令 k=k+1,重 复 上 面 的 步 骤 2 到
4,直到k=N。 3〃2 改进算法导航电文译码性能分析 为考察算法的改进效果 ,对所给 (2,1,3)卷 积码进行了算法改进前后的 MonteCarlo仿真。 仿真条 件:采 用 等 先 验 概 率 ,信 号 采 用 0→ -1、1→1映射 双极性信号,两 种 MAP 算法中的
支路转移度量均按式(10
)来计算。仿真按照接收 (13)式会以很 大概率得到满足。 反之,由 式 (13) 可得2γt<1,且该结果不随k 的变化而变化,这样 由式(9)和式(7)可 知,当输入码序列的长度越来 越大时,会 使αk (pk
)越 来 越 小,最 终 导 致 计 算 机 产生溢出而 变 为 0。 同 时,由 式 (10)~ 式 (13)
可 知,2γt<1 成立的条件只与噪声 ‖nt‖2 有关,
而 与信 号 无 关,这 样 当 信 号 有 衰 落 时,上 述 结 论 不
变。由此可知,在用 MAP 算法译码时,
随 着码长 的增加,计算机就会产生数值稳 定性问题。 对于
其他卷积码也可以得到类似的结论 ,只 是式 (13
)
第39卷第3期胡毅等:导航卫星卷积码序列MAP译码算法分析与改进269码序列信号y有无衰落两种情况进行。为检验
算法的正确性,同时给出维特比软译码算法仿真
结果作为对比。
Eb
(N0)dB
号。在每个点,仿真200次。由于导航电
文多普勒频移非常小,因此仿真中将其取为0。
图3(a)是在具有直达信号的莱斯平坦衰落
信道中MAP改进算法译码结果,这是在室外或
具有较少遮挡物下的导航接收机典型工作场景。
从图中可以看出,随着直达信号与多径信号功率
比值K的增加,衰落对译码性能造成的影响得到
明显改善。
图3(b)是在瑞利平坦衰落信道即没有直达
信号,并在不同码速率时改进的MAP算法译码
结果,这是在室内或具有严重遮挡下的导航接收
机的工作场景。从图中可以看出,在信号有较大
衰落的情况下,改进的MAP译码算法也能给出正
常的算法结果,但这时的误码率比无衰落时有明显
的增加,并且这种性能变差随着码速率R的增加
而恶化,这与接收机的实际工作特性相符。
信号无衰落时的算法性能仿真
3〃2〃1
这是一种理想情形,此时接收码序列可表示
为yk=ck+nk(k=1,2,…,N)。不同码序列长度
的译码仿真结果如图2所示,图中MAP(R)表示
改进的归一化MAP译码算法。仿真中编码器输java视频教程免费下载
入信号采用的是从GPS软件接收机输出中提取
的C/A码单子帧导航电文即300bit数据信号,
和从中任意截取的200个数据信号。在每个
Eb
(N0)dB点,(a)仿真1000次,(b)仿真600次。
学python还是c++好从图中可以看出,当码序列长度N较小时,
MAP(R)和MAP译码算法性能相同,如图2(a)
所示。但是,当输入码序列超过一定长度时,
MAP算法由于溢出,输出了错误的译码误码率,
在仿真图中表现为一条误码率水平起伏或基本不
变的水平线,如图2(b)所示。而此时MAP(R)算
法由于采用了状态度量归一化,改变了数值变化
范围,因而仍能给出正确的算法结果。
信号有衰落时的算法性能仿真
3〃2〃2
此时接收码序列可表示为yk=μkck+nk(k=
1,2,…,N),其中μk是k时刻服从一定分布的信
号衰落复包络乘性系数,它直接反映了信号衰落
的程度,值越小,则衰落越大[9]。(a)对于莱斯图2无衰落时不同码序列长度的算法性能仿真
Fig〃2AlgorithmPerformanceSimulationswith
DifferentCodeLengthandwithoutFading
(Rice)平坦衰落信道,
μk可表示为
2
(ξ
k
+槡2K)+δk2
,其中K为莱斯因数,常以dB
槡
2(K+1)
为单位,ξk、δk相互独立且服从N(0,1)分布。
[
9]
(b)对于具有多普勒频移的瑞利(Rayleigh)平坦
衰落信道,μk可表示其中ξk、δk可
下式近似得到[9]:
M
2
cos(ωdkcosαn+φn)
槡Mn∑=1
ξk=
(15)图3信号有衰落时的算法性能仿真
M
2
cos(ωdkcosαn+φn)Fig〃3AlgorithmPerformanceSimulations
withSignalFading
槡Mn∑=1
ξk=
式中,ωd为最大多普勒频移,αn
=2πfd=
2πn-π+θn,n
=1,2,…,M,θ和φn服从
n
、φ
n结语
4M4
[-π,π)上的均匀分布且相互统计独立,M为所
用正弦波数目。
两种衰落信道下的译码仿真结果如图3所
示。仿真中输入信号采用的是从软件接收机获取
的C/A码一帧完整导航电文即1500bit数据信
本文针对导航卫星卷积码序列MAP译码算
法的不足进行了具体分析,并给出了一种改进的
归一化MAP译码算法,同时对改进后的算法进
行了不同码长和不同信道条件下的译码性能仿
武汉大学学报·信息科学版2014年3月270
真。结果表明,改进后的MAP译码算法对于不同情况下的卷积码序列都能保证数值稳定性,给出正确的算法结果,因而更适合于导航卫星卷积码MAP译码的实际应用。
参考文献
GPS-200F[S]〃GPSJointProgramOffice,2011
[5]PenaAG,BoucheretML,MacabiauC,etal〃AnImprovedMethodtoDecodeGPSL2C/L5Naviga-tionMessage:CombinationoftheInnerandtheOuterChannelCodes[C]〃The2009InternationalTechnicalMeetingoftheInstituteofNavigation,Anaheim,2009
[6]BerkanD,SinanG〃DetectorRandomizationandStochasticSignalingforMinimumProbabilityofEr-r
orReceivers[J]〃IEEETransactionsonCommuni-ca
tions,2012,60(4):923-928
[1]ChenNan〃PerformanceEvaluationoftheStructureofGNSS[J]〃GeomaticsandInformationScienceo
fWuhanUniversity,2008,33(5):512-515(陈
南〃卫星导航系统导航电文结构的性能评估
[J]〃武汉大学学报·信息科学版,2008,33
(5):512-515)
TurunenS〃CanCodeRedundancybeUsedtoIm-proveGNSSReceiverAcquisitionSensitivity[C]〃2006IEEE/IONPosition,LocationandNavigationSymposium,SanDiego,2006ChenJinping,WangMengli,QianShuguang〃A-nalysisofModernizationGNSSNavigationMes-
sage’sDesigning[J]〃JournalofElectronics&In-
formationTechnology,2011,33(1):211-217(陈
金平,王梦丽,钱曙光〃现代化GNSS导航电文设计
分析[J]〃电子与信息学报,2011,33(1):211-217)
GlobalPositingSystemDirectorateSystemsEngi-neering&IntegrationInterfaceSpecification,IS-[7]RobertsonP,VillebrunE,HoecherP〃ACompari-sonofOptimalandSub-optimalMAPDecodingAl-gorithmsOperationintheLogDomain[C]〃Interna-tionalConferenceonCommunications1995,Seattle,
1995
[2]
[8]StevenSP,AdrianSB〃ASimplificationofthe
ModifiedBahlDecodingAlgorithmforSystematicC
onvolutionalCodes[C]〃InternationalSymposiumonI
nformationTheory&itsApplications(ISITA),Sydn
ey,1994
[3]
[9]NikolayK〃MobileRadioChannelsModelingin
MATLAB[J]〃RadioEngineering,2003,12(7):1
2-16
[4]
AnalysisandImprovementoftheMAPAlgorithmforNavigation
SatelliteConvolutionalCodesDecoding
HUYi1SONGMaozhong1,2MENGBin2
CollegeofElectronicInformationEngineering,NanjingUniversityofAeronautics&Astronautics,Nanjing210016,China1
TheAcademyofSatelliteApplication,ChinaAerospaceScienceandTechnologyCorporation,Beijing100086,China2position和location的区别
Abstract:Computersmayencountertheproblemofnumericalinstabilityandoutputawrongbiterror
authorizeattribute validatetokenrate(BER)whendecodingnavigationsatelliteconvolutionalcodeswiththemaximumaposterioriprobability(MAP)algorithm〃Subsequently,thiswillaffectthevalidityandtheaccuracyoftheposi-tioningresults〃Toaddressthisproblem,firstadetailedanalysiswasexecuted〃Basedonthisanaly-sis,animprovedMAPalgorithmnormalizedwiththemaximum
statemetricisproposed,thusthequestionofnumericalinstabilitycanbesolvedwhilestillmaintainingtheoriginalalgorithm’soptimalperformance〃Simulationsfurthervalidatetheeffectivenessofthisalgorithm’simprovement〃Keywords:navigationsatelliteconvolutionalcodes;MAPdecodingalgorithm;numericalstability;maximumstatemetric
Firstauthor:HUYi,PhDcandidate,specializesinnavigationsatellitesignalprocessing〃E-mail:hystarcom@163〃com
Foundationsupport:TheOpenResearchFundoftheAcademyofSatelliteApplication,No〃20121512〃
版权声明:本站内容均来自互联网,仅供演示用,请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系QQ:729038198,我们将在24小时内删除。
发表评论