专利名称：一种基于盲信号处理的北斗导航电文提取装置和方法
技术领域：
本发明涉及一种北斗卫星导航系统领域的技术，更具体的说，涉及一种能够在未知北斗卫星导航系统伪随机码情况下，提取北斗导航电文的装置和方法。
背景技术：
全球导航卫星定位系统(GNSS)在过去几十年间在各个领域已经得到了广泛的应用。目前GNSS包括了美国的GPS、俄罗斯的GL0NASS、中国的北斗(Compass)以及欧盟的Galileo系统。作为中国自主研发的卫星导航系统，北斗系统将在市场中获得越来越重要的地位。使用北斗系统进行定位、授时等应用的一个必要步骤是提取北斗导航电文技术，因此如何有效提取北斗导航电文是当前北斗系统实用化必须解决的重要问题。目前传统的北斗导航电文提取方案从GPS等其他GNSS系统电文提取方案发展而来，通过捕获和跟踪北斗信号得到北斗导航电文。为了提取北斗导航电文，北斗接收机需要捕获来自不同卫星的北斗信号，解调出北斗信号的导航电文。不同卫星的北斗信号具有不同的起始时间的伪随机码和不同的多普勒频移。因此，为了搜索某个卫星信号，北斗接收机通常进行二维搜索，在每个可能的多普勒频移上对每个起始时间不同的伪随机码进行搜索。接收机捕获到北斗卫星信号后将伪随机码的起始时间和多普勒频移输入至跟踪�？椋倌？槟芄皇凳备僭夭ê臀彼婊氡浠佣袢【范任彼婊胂嘁坪投嗥绽掌狄埔越獾鞒霰倍沸藕胖邪牡己降缥�。尽管传统方案能够有效提取北斗导航电文，但是北斗系统与其他GNSS系统的重要区别使得传统卫星导航技术并不是针对北斗系统最优。对于北斗GEO卫星，因为其轨道是相对地球静止的，其发射的北斗信号的载波多普勒只由接收机运动引起，因此在已知接收机速度的情况下，北斗信号的载波多普勒基本可以估算出来。这样传统基于捕获跟踪的方案的二维搜索将造成资源的不必要浪费。同时传统接跟踪�？橹性夭ɑ仿泛臀彼婊牖仿肥腔诜蠢〉幕仿罚杓聘丛樱仿凡问舾校诟丛有藕呕肪诚戮Ｊ�。
综上所述，急需一个有效方案来提高北斗导航电文提取效率，简化设计，使其获得更广泛的应用。该方案能克服传统方案三维搜索中资源浪费的缺点，能避开反馈环路设计中的复杂性和不稳定性，并且能够在未知北斗卫星导航系统伪随机码序列情况提取北斗导航电文。
发明内容
针对上述现有技术中存在的技术问题，本发明提供一种基于盲信号处理的北斗导航电文的提取装置和方法。该装置和方法不需要事先知道北斗卫星导航系统的伪随机码，并能绕过捕获和跟踪等传统步骤，直接利用输入数据提取北斗导航电文，算法复杂度低，非常适用于软件接收机。为达到上述发明目的，本发明所采用的技术方案如下一种基于盲信号处理的北斗导航电文提取装置，该装置包括依次连接的中频采样�？�、数据分段�？椤⑿讲罟兰颇？楹推嬉熘捣纸饽？椋鍪莘侄文？榈氖涑龆擞胨銎嬉熘捣纸饽？榈氖涑龆司拥揭恍蜃蛹扑隳？楹鸵坏己降缥幕指茨？樯希遥鲂蜃蛹扑隳？榈氖涑龆肆拥剿龅己降缥幕指茨？樯�。一种基于盲信号处理的北斗导航电文提取方法，采用上述的基于盲信号处理的北斗导航电文提取装置，包括步骤如下I)、所述中频采样�？榻邮帐淙肽Ｄ庵衅敌藕牛媚Ｊ蛔爸枚允淙肽Ｄ庵衅敌藕沤心Ｊ徊裳玫绞种衅敌藕牛�2)、所述数据分段�？榻种衅敌藕乓员倍肺彼婊胄蛄兄芷谖侄沃芷诮蟹侄危玫椒侄魏蟮氖种衅敌藕牛侄魏蟮氖种衅敌藕殴钩梢桓鲋衅凳菥卣螅�3)、所述协方差估计�？槎灾衅凳菥卣蠼行讲罟兰疲玫叫讲罹卣螅�4)、所述奇异值分解�？槎孕讲罹卣蠼衅嬉熘捣纸猓玫阶畲笃嬉熘刀杂Φ姆纸饩卣笙蛄浚孟蛄坑糜谛蜃蛹扑隳？楹偷己降缥幕指茨？椋�5)、所述旋转因子计算�？槔弥衅凳菥卣蠛头纸饩卣笙蛄考扑愠鱿嗔诹蕉畏侄沃衅敌藕诺难≡褚蜃樱�6)、所述导航电文恢复�？槔盟鲂蜃樱衅凳菥卣笠约胺纸饩卣笙蛄拷械己降缥幕指吹玫奖倍返己降缥摹Ｋ鲋衅挡裳？榈牟裳士傻鳎裳菸豢砜傻�。所述奇异值分解�？榻邮招讲罹卣笞魑淙耄銎嬉熘捣纸猓玫狡嬉熘敌蛄校闷嬉熘敌蛄薪薪敌蚺判颍玫阶畲笃嬉熘担⑻崛∽畲笃嬉熘刀杂Φ姆纸饩卣笙蛄俊Ｋ鲂蜃蛹扑隳？榧扑阒衅凳菥卣蠛头纸饩卣蟮南蛄康哪诨⒆鞴橐换缘玫降南蛄磕诨蛄薪锌焖俑盗⒁侗浠唬玫叫蜃�。所述导航电文提取�？槔眯蜃�、数据矩阵和分解矩阵向量得到带有一定相角导航电文，通过分析带有一定相角导航电文的相角和模，得到只有符号�：鹊牡己降缥模ü己降缥牡闹〗峁固匦韵玫缥牡姆拍：�。与现有技术相比，本发明具有如下有益效果1、本发明利用盲信号处理技术，在未知北斗卫星导航系统伪随机码序列情况下仍能提取北斗卫星导航电文；2、本发明避免传统接收机方案中占用很多空间和时间资源的捕获步骤，不需要已知每颗可见卫星信号的伪随机码起始时间和多普勒频移就能对北斗信号进行处理；3、本发明绕开传统接收机方案中设计复杂的载波环路和伪随机码环路的跟踪步骤，实施复杂度降低，尤其适用于软件接收机设计。
图1是本发明所述的基于盲信号处理的北斗导航电文提取装置的系统框图；图2是本发明所述的基于盲信号处理的北斗导航电文提取方法的工作流程图；图3是本发明装置中奇异值分解模块的工作流程图；图4是本发明装置中旋转因子计算�？榈墓ぷ髁鞒掏迹煌�5是本发明装置中导航电文恢复�？榈墓ぷ髁鞒掏迹煌�6是本发明装置中要处理的北斗信号结构图。
具体实施例方式下面结合附图和具体实施例对本发明的技术方案作详细说明，本实施例在以本发明技术方案为前提下进行实施，给出了详细的实施方式和具体的操作过程，但本发明的保护范围不限于下述的实施例。如图1所示，本发明所公开的一种基于盲信号处理的北斗导航电文提取装置，包括依次连接的中频采样�？�102、数据分段�？�104、协方差估计模块106和奇异值分解�？�108，所述数据分段�？�104的输出端与所述奇异值分解�？�108的输出端均连接到一旋转因子计算�？�110和一导航电文恢复�？�112上，且，所述旋转因子计算�？�110的输出端连接到所述导航电文恢复�？�112上。以下是本发明所公开的采用上述装置来对基于盲信号处理的北斗导航电文提取方法的详细解释和说明。图2所示，中频采样模 块102接收输入模拟中频信号y(t)，利用模数转换装置对输入模拟中频信号y(t)进行模数转换采样，得到数字中频信号y(t JKTs);数据分段�？�104将数字中频信号y(tJKTs)以北斗伪随机码序列周期为分段周期进行分段，得到分段后的数字中频信号，分段后的数字中频信号构成一个中频数据矩阵yn ;协方差估计�？�106对中频数据矩阵7 进行协方差估计，得到协方差矩阵R ;奇异值分解模块108对协方差矩阵R进行奇异值分解，得到最大奇异值对应的分解矩阵向量V，该向量用于旋转因子计算�？楹偷己降缥幕指茨？椋恍蜃蛹扑隳？�110利用中频数据矩阵^和分解矩阵向量V计算出两段分段中频信号的旋转因子;导航电文恢复�？�112利用所述旋转因子，中频数据矩阵yn以及分解矩阵向量V进行导航电文恢复得到北斗导航电文s (t)。异值分解�？�108接收协方差矩阵R作为输入，做奇异值分解，得到奇异值序列；将该奇异值序列进行降序排序，得到最大奇异值，并提取最大奇异值对应的分解矩阵向量
Vo旋转因子计算�？�110计算中频数据矩阵yn和分解矩阵向量V的内积，并作归一化对得到的向量内积序列进行快速傅立叶变换，得到旋转因子。导航电文恢复�？�112利用旋转因子、中频数据矩阵7 和分解矩阵向量V得到带有一定相角导航电文，通过分析带有一定相角导航电文的相角和模，得到只有符号�：鹊牡己降缥模ü己降缥牡闹〗峁固匦韵玫缥牡姆拍：龋玫奖倍返己降缥膕 (t)。以下是本发明一具体实施例根据北斗的信号规范，卫星信号结构如图5所示。对北斗信号通过正交接收并转换成中频信号后，中频信号(IQ两路信号组成的复数)表达式为x (t) =eJ(Mt+$)c (t) s (t)，其中Φ是初始相位，ω是多普勒角频率，c (t)是从某个未知初始码相位开始的取值{±1}的扩频码序列，码片速率2. 046Mcps,扩频码周期是lms, s(t)是取值{±1}的导航电文。实际接收的卫星信号信噪比极低，噪声功率是信号的1(Γ500倍，即I (t) =x (t) +u (t)，其中 u (t)是噪声。中频采样�？�102接收的输入模拟中频信号(本实施例中频7.098MHz)，利用模数转换装置以2. 046MHz采样率对输入模拟中频信号采样，采样200ms，得到200X2046个数字中频信号{y (0)，y (Ts)，y (2TS)，y (3TS)，L}，并存储在存储介质，其中Ts=10_6/2. 046是采样间隔。数据分段�？楦堇┢德�1ms的周期性将数字中频信号以1ms周期进行分段，得到分段后的数字中频信号，其中分段后的某段1ms中频数据可表示为yn=[y (2046nTs) y (2046nTs+Ts) L y (2046nTs+2045Ts) ]τ,分段后的数字中频信号构成一个200行2046列的中频数据矩阵yn，我们将7 表示成信号与噪声叠加yn=xn+un，un是对噪声采样构成的向量，满足
权利要求
1.一种基于盲信号处理的北斗导航电文提取装置，其特征在于，该装置包括依次连接的中频采样�？�、数据分段�？�、协方差估计�？楹推嬉熘捣纸饽？椋鍪莘侄文？榈氖涑龆擞胨銎嬉熘捣纸饽？榈氖涑龆司拥揭恍蜃蛹扑隳？楹鸵坏己降缥幕指茨？樯希遥鲂蜃蛹扑隳？榈氖涑龆肆拥剿龅己降缥幕指茨？樯�。
2.一种基于盲信号处理的北斗导航电文提取方法，其特征在于，采用如权利要求1所述的基于盲信号处理的北斗导航电文提取装置，包括步骤如下 1)、所述中频采样�？榻邮帐淙肽Ｄ庵衅敌藕牛媚Ｊ蛔爸枚允淙肽Ｄ庵衅敌藕沤心Ｊ徊裳玫绞种衅敌藕牛� 2)、所述数据分段�？榻种衅敌藕乓员倍肺彼婊胄蛄兄芷谖侄沃芷诮蟹侄危玫椒侄魏蟮氖种衅敌藕牛侄魏蟮氖种衅敌藕殴钩梢桓鲋衅凳菥卣螅� 3)、所述协方差估计�？槎灾衅凳菥卣蠼行讲罟兰疲玫叫讲罹卣螅� 4)、所述奇异值分解模块对协方差矩阵进行奇异值分解，得到最大奇异值对应的分解矩阵向量，该向量用于旋转因子计算�？楹偷己降缥幕指茨？椋� 5)、所述旋转因子计算�？槔弥衅凳菥卣蠛头纸饩卣笙蛄考扑愠鱿嗔诹蕉畏侄沃衅敌藕诺难≡褚蜃樱� 6)、所述导航电文恢复�？槔盟鲂蜃樱衅凳菥卣笠约胺纸饩卣笙蛄拷械己降缥幕指吹玫奖倍返己降缥�。
3.根据权利要求2所述的基于盲信号处理的北斗导航电文提取方法，其特征在于，所述中频采样�？榈牟裳士傻鳎裳菸豢砜傻�。
4.根据权利要求2所述的基于盲信号处理的北斗导航电文提取方法，其特征在于，所述奇异值分解�？榻邮招讲罹卣笞魑淙耄銎嬉熘捣纸猓玫狡嬉熘敌蛄校闷嬉熘敌蛄薪薪敌蚺判颍玫阶畲笃嬉熘担⑻崛∽畲笃嬉熘刀杂Φ姆纸饩卣笙蛄俊�
5.根据权利要求2所述的基于盲信号处理的北斗导航电文提取方法，其特征在于，所述旋转因子计算�？榧扑阒衅凳菥卣蠛头纸饩卣蟮南蛄康哪诨⒆鞴橐换缘玫降南蛄磕诨蛄薪锌焖俑盗⒁侗浠唬玫叫蜃�。
6.根据权利要求2所述的基于盲信号处理的北斗导航电文提取方法，其特征在于，所述导航电文提取�？槔眯蜃�、数据矩阵和分解矩阵向量得到带有一定相角导航电文，通过分析带有一定相角导航电文的相角和模，得到只有符号�：鹊牡己降缥模ü己降缥牡闹〗峁固匦韵玫缥牡姆拍：�。
全文摘要
本发明公开一种基于盲信号处理的北斗导航电文提取装置和方法，包括对输入的中频模拟信号进行采样得到数字中频信号；基于伪随机码周期对数字中频信号进行分段得到中频数据矩阵；对分段数字中频信号进行协方差估计得到协方差矩阵；对协方差矩阵进行奇异值分解得到最大奇异值以及对应的分解矩阵向量；基于分解矩阵向量和所述中频数据矩阵得到旋转因子；利用所述旋转因子，中频数据矩阵以及分解矩阵向量进行导航电文恢复得到北斗导航电文。本发明基于盲信号处理，相比较常规基于捕获、跟踪的电文提取方法和装置，算法复杂度低，不需得到北斗卫星信号中的初始码相位，甚至不需事先知道北斗伪随机码序列。
文档编号G01S19/35GK103064096SQ20121029956
公开日2013年4月24日 申请日期2012年8月21日 优先权日2012年8月21日
发明者许良备, 钱久超, 应忍冬, 刘佩林, 郁文贤 申请人:上海交通大学, 上海美迪索科电子科技有限公司
专业数控铣床产品供应商
20年生产经验，实力工厂，精工品质，质量有保障
