CN101859292A - 通信装置、信号处理方法、信号处理装置以及移动体 - Google Patents
通信装置、信号处理方法、信号处理装置以及移动体 Download PDFInfo
- Publication number
- CN101859292A CN101859292A CN201010151265A CN201010151265A CN101859292A CN 101859292 A CN101859292 A CN 101859292A CN 201010151265 A CN201010151265 A CN 201010151265A CN 201010151265 A CN201010151265 A CN 201010151265A CN 101859292 A CN101859292 A CN 101859292A
- Authority
- CN
- China
- Prior art keywords
- mentioned
- high speed
- fourier transform
- speed fourier
- execution section
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000012545 processing Methods 0.000 title claims abstract description 23
- 238000003672 processing method Methods 0.000 title claims abstract description 8
- 238000007781 pre-processing Methods 0.000 claims abstract description 117
- 238000002203 pretreatment Methods 0.000 claims abstract description 47
- 238000005070 sampling Methods 0.000 claims abstract description 22
- 238000012360 testing method Methods 0.000 claims description 23
- 230000000737 periodic effect Effects 0.000 claims description 21
- 238000009825 accumulation Methods 0.000 claims description 18
- 238000005520 cutting process Methods 0.000 claims description 8
- 238000001228 spectrum Methods 0.000 claims description 8
- 230000001360 synchronised effect Effects 0.000 description 76
- 238000010586 diagram Methods 0.000 description 38
- 230000015572 biosynthetic process Effects 0.000 description 29
- 238000005755 formation reaction Methods 0.000 description 29
- 230000009471 action Effects 0.000 description 23
- 230000035508 accumulation Effects 0.000 description 12
- 238000006243 chemical reaction Methods 0.000 description 8
- 230000006870 function Effects 0.000 description 6
- 238000000034 method Methods 0.000 description 6
- 230000008569 process Effects 0.000 description 6
- 238000003860 storage Methods 0.000 description 5
- 230000009467 reduction Effects 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 230000008676 import Effects 0.000 description 3
- 238000012423 maintenance Methods 0.000 description 3
- 230000000996 additive effect Effects 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000000284 extract Substances 0.000 description 2
- 230000010363 phase shift Effects 0.000 description 2
- 230000007480 spreading Effects 0.000 description 2
- 238000003892 spreading Methods 0.000 description 2
- 238000010897 surface acoustic wave method Methods 0.000 description 2
- 239000000654 additive Substances 0.000 description 1
- 230000003750 conditioning effect Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000010453 quartz Substances 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N silicon dioxide Inorganic materials O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/24—Acquisition or tracking or demodulation of signals transmitted by the system
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/35—Constructional details or hardware or software details of the signal processing chain
- G01S19/37—Hardware or software details of the signal processing chain
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Signal Processing (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
- Analogue/Digital Conversion (AREA)
- Arrangements For Transmission Of Measured Signals (AREA)
Abstract
本发明提供通信装置、信号处理方法、信号处理装置及移动体。通信装置具备:采样部,其将接收信号采样;前处理执行部,其将利用所述采样部得到的采样值一边实施为了分割进行高速傅立叶变换而进行的前处理一边取得;高速傅立叶变换部,其基于来自所述前处理执行部的输出,分割进行高速傅立叶变换。
Description
技术领域
本发明涉及通信装置、信号处理方法、信号处理装置以及移动体。
背景技术
近年来,能够接收从人造卫星发送的导航消息来计算当前位置的GPS(Global Positioning System)接收机被应用于携带电话或导航系统中等而得到广泛地普及。
具体来说,在由人造卫星发送的导航消息中,包含表示人造卫星的轨道的轨道信息、信号的发送时刻等信息。GPS接收机从4个以上的人造卫星接收上述导航消息,根据包含于导航消息中的轨道信息计算出各人造卫星的位置。而且,GPS接收机可以基于各人造卫星的位置、导航消息的发送时刻与接收时刻的差值,利用联立方程计算当前的三维位置。而且,在计算三维位置时需要从4个以上的人造卫星发送的导航消息是因为,在内置于GPS接收机中的时钟与设于人造卫星中的原子钟间存在误差。
另外,人造卫星利用被称作L1频段C/A码的信号,即码长度为1,023而码片速率为1.023MHz的伪随机(PRN:Pseudo-Random Noise)码将50bps的数据进行频谱扩展,利用频谱扩展后的信号将1,575.42MHz的载波进行BPSK(Binary Phase Shift Keying,相位移相键控)调制,利用该调制后的信号发送上述导航消息。
所以,GPS接收机为了接收来自人造卫星的信号,就需要取得PRN码、载波及数据的同步。接收信号的PRN码的同步是通过使用例如高速傅立叶变换运算接收信号的PRN码与在本装置中产生的PRN码之间的相关性来实现的。
另外,在专利文献1中,记载了通过将用于获得全部频率的1/L(L为整数)的频率成分的高速傅立叶变换进行L次来得到全部频率成分的接收装置。根据该接收装置,可以削减保持高速傅立叶变换结果的存储器的使用量。
专利文献1日本专利第3906913号
但是,在以往的接收装置进行的L次的高速傅立叶变换中,包含很多重复的运算。所以,在以往的接收装置中,存在运算量及运算时间增加的问题。
发明内容
所以,本发明是鉴于上述问题而完成的,本发明的目的在于,提供可以缩短用于高速傅立叶变换的运算时间的、新型且改进的通信装置、信号处理方法、信号处理装置及移动体。
为了解决上述问题,根据本发明的一个观点,提供一种通信装置,具备:采样部,其对接收信号进行采样;前处理执行部,其将利用上述采样部得到的采样值一边实施为了分割进行高速傅立叶变换而进行的前处理一边取得;和高速傅立叶变换部,其基于来自上述前处理执行部的输出,分割进行高速傅立叶变换。
也可以是,上述前处理执行部具备:信号生成部,其生成周期信号;乘法部,其将利用上述采样部得到的采样值和利用上述信号生成部生成的上述周期信号的值相乘;和累积部,其将利用上述乘法部得到的各个乘积值依照利用上述乘法部获得的顺序以规定间隔累积,得到多个累积值,其中,上述高速傅立叶变换部对利用上述累积部得到的上述多个累积值进行高速傅立叶变换。
也可以是,上述通信装置具备多个上述前处理执行部,每个上述前处理执行部所包含的上述信号生成部生成具有不同的周期的周期信号。
也可以是,每个上述前处理执行部的上述累积部以接收信号的规定长度为单位进行各个上述乘积值的累积,上述通信装置还具备与每个上述前处理执行部对应的存储器,在上述存储器中,记录利用对应的前处理执行部得到的上述多个累积值,上述高速傅立叶变换部将依次记录于每个上述存储器中的上述多个累积值进行高速傅立叶变换。
也可以是,上述接收信号被利用扩展码进行了频谱扩展,上述通信装置还具备:第二前处理执行部,其包括生成与第一前处理执行部的上述信号生成部具有相同的周期的周期信号的信号生成部,接收信号被以上述规定长度为单位划分的位置与上述第一前处理执行部中不同;相关点检测部,其基于上述高速傅立叶变换部的高速傅立叶变换结果,检测上述扩展码的相关点;和选择部,其选择基于利用上述第一前处理执行部得到的所述多个累积值的高速傅立叶变换结果由相关点检测部检测出的相关点和基于利用上述第二前处理执行部得到的上述多个累积值的高速傅立叶变换结果由相关点检测部检测出的相关点中的任意一个。
也可以是,上述累积部以接收信号的规定长度为单位进行各个上述乘积值的累积,上述信号生成部以接收信号的上述规定长度为单位变更要生成的周期信号的周期。
另外,为了解决上述问题,根据本发明的其他观点,提供一种信号处理方法,包括:对接收信号进行采样的步骤;将通过采样得到的采样值一边实施为了分割进行高速傅立叶变换而进行的前处理一边记录于存储器中的步骤;和基于来自上述存储器的输出分割进行高速傅立叶变换的步骤。
另外,为了解决上述问题,根据本发明的其他观点,提供一种信号处理装置,具备:采样部,其对接收信号进行采样;前处理执行部,其将利用上述采样部得到的采样值一边实施为了分割进行高速傅立叶变换而进行的前处理一边取得;和高速傅立叶变换部,其基于来自上述前处理执行部的输出分割进行高速傅立叶变换。
另外,为了解决上述问题,根据本发明的其他观点,提供一种移动体,其搭载有如下的通信装置,即,该通信装置具备:采样部,其对接收信号进行采样;前处理执行部,其将利用上述采样部得到的采样值一边实施为了分割进行高速傅立叶变换而进行的前处理一边取得;和高速傅立叶变换部,其基于来自上述前处理执行部的输出分割进行高速傅立叶变换。
如上说明所示,根据本发明的通信装置、信号处理方法、信号处理装置及移动体,可以缩短用于高速傅立叶变换的运算时间。
附图说明
图1是表示本发明的实施方式的人造卫星系统的构成的说明图。
图2是表示导航消息的帧构成的说明图。
图3是表示本发明的实施方式涉及的接收机的硬件构成的说明图。
图4是表示本发明的实施方式涉及的接收机的动作例的流向的流程图。
图5是表示与本发明相关的同步捕捉部的构成的说明图。
图6是表示输入数据为16个时的高速傅立叶变换的运算的流向的说明图。
图7是表示分割进行高速傅立叶变换时的运算的流向的说明图。
图8是表示分割进行高速傅立叶变换时的运算的流向的说明图。
图9是表示分割进行高速傅立叶变换时的运算的流向的说明图。
图10是表示分割进行高速傅立叶变换时的运算的流向的说明图。
图11是定义在高速傅立叶变换的运算过程中得到的中间数据的说明图。
图12是表示第一实施方式涉及的同步捕捉部的构成的说明图。
图13是表示第一实施方式涉及的同步捕捉部的动作的流向的说明图。
图14是表示第二实施方式涉及的同步捕捉部的构成的说明图。
图15是表示第二实施方式涉及的同步捕捉部的动作的流向的说明图。
图16是表示第三实施方式涉及的同步捕捉部的构成的说明图。
图17是表示第三实施方式涉及的同步捕捉部的动作的流向的说明图。
图18是表示第四实施方式涉及的同步捕捉部的构成的说明图。
图19是表示第四实施方式涉及的同步捕捉部的动作的流向的说明图。
图20是表示第五实施方式涉及的同步捕捉部的构成的说明图。
图21是表示第五实施方式涉及的同步捕捉部的动作的流向的说明图。
图中符号:
240、241、242、243、244同步捕捉部
310、310A、310B、310C、310D前处理执行部
312信号生成部
314乘法部
316加法部
318、320、320A、320B、320C、320D、332、348存储器
330、344FFT
336运算控制部
356反FFT
364选择部
具体实施方式
下面,在参照附图的同时,对本发明的优选的实施方式进行详细说明。而且,本说明书及附图中,对于具有实质上相同的功能构成的构成要素,通过使用相同的符号而省略重复说明。
另外,本说明书及附图中,也有将具有实质上相同的功能构成的多个构成要素在相同的符号后附加不同的字母而进行区别的情况。例如,将具有实质上相同的功能构成的多个构成根据需要区分为人造卫星10A、10B及10C。但是,在不需要特意地区别各个具有实质上相同的功能构成的多个构成要素时,仅标以相同的符号。例如,在不需要特意地区别人造卫星10A、10B及10C的情况下,单纯地称作人造卫星10。
另外,依照以下所示的项目顺序说明对应的“用于实施发明的方式”。
1.人造卫星系统的概要
1-1.利用GPS的位置定位的概要
1-2.导航消息的构成
1-3.接收机的构成及动作
2.本发明的实施方式的背景
3.第一实施方式
4.第二实施方式
5.第三实施方式
6.第四实施方式
7.第五实施方式
8.本发明的实施方式的效果
<1.人造卫星系统的概要>
首先,参照图1~图4,对本发明的实施方式涉及的人造卫星系统1进行说明。
[1-1.利用GPS的位置定位的概要]
图1是表示本发明的实施方式涉及的人造卫星系统1的构成的说明图。如图1所示,该人造卫星系统1具备多个人造卫星10A~10D、接收机20。
人造卫星10(GPS卫星)构成定位卫星系统(GNSS:GlobalNavigation Satellite System,全球卫星定位系统),环绕在地球8的上空。图1中,仅图示了4个人造卫星10A~10D,然而例如在6个轨道面上各4个共计24个人造卫星环绕在地球8的上空。
另外,人造卫星10发送包括人造卫星的轨道信息、导航消息的发送时刻等广播星历(Ephemeris)信息的导航消息(关于详情在[1-2.导航消息的构成]中说明。)。在人造卫星10中设有原子钟,发送时刻是依照设于人造卫星10中的原子钟例如以1秒为单位表现的。
而且,人造卫星10利用被称作L1频段C/A码的频谱扩展信号发送上述导航消息。该频谱扩展信号是将50bps的数据利用伪随机序列的扩展码(以下称作PRN码。)进行频谱扩展,利用频谱扩展后的信号将1,575.42MHz的载波进行BPSK调制后的信号。另外,PRN码的码长度为1,023,码片速率为1.023MHz。
存在于地球8上的接收机20可以接收从人造卫星10A~10D发送的导航消息,基于接收到的导航消息计算本装置的当前位置。
更具体来说,接收机20接收从人造卫星10A~10D发送的导航消息,从导航消息中取得广播星历信息。而且,接收机20根据广播星历信息计算出人造卫星10A~10D的位置。另外,接收机20根据广播星历信息中所包含的发送时刻与导航消息的接收时刻的差值,计算出人造卫星10A~10D与接收机20之间的距离。其后,接收机20利用计算出的人造卫星10A~10D的各自的位置、和各个人造卫星10A~10D与接收机20之间的距离,运算以接收机20的当前的三维位置作为未知数的方程式。
在像这样运算接收机20的当前的三维位置时,需要从4个以上的人造卫星10发送的导航消息。这是因为,在内置于接收机20中的时钟(RTC:Real Time Clock)与设置于人造卫星10中的原子钟问存在误差。
另外,人造卫星10以规定周期对广播星历信息进行更新,发送包含更新后的广播星历信息的导航消息。这里,由于人造卫星10总是在移动着,因此从广播星历信息的更新时起随着时间的经过,基于广播星历信息计算出的人造卫星10的位置与实际的人造卫星10的位置之间的误差不断变大。由此,对于导航消息中包含的广播星历信息例如设定2小时左右的有效期限。
以上,参照图1对利用GPS的位置定位进行了概略性说明。而且,图1中将接收机20作为通信装置的一例以圆圈标记表示,然而接收机20也可以是PC(Personal Computer)、家用影像处理装置(DVD录像机、盒式录像机等)、携带电话、PHS(Personal Handyphone System,低功率移动电话系统)、携带用音乐再现装置、携带用影像处理装置、PDA(Personal Digital Assistants,个人数码助理)、家用游戏机、携带用游戏机、家电机器及车载机器等信息处理装置。
另外,接收机20也可以搭载于能够使人或货物等搭乘而移动的移动体中。例如,作为移动体,可以举出摩托车或三轮车等汽车、自行车、公共汽车、电力火车、新干线、有轨电车、飞机、船、小艇等多种多样的交通工具。
[1-2.导航消息的构成]
下面,参照图2,对从人造卫星10发送的导航消息的构成进行说明。
图2是表示导航消息的帧构成的说明图。如图2所示,导航消息的1个帧由5个子帧构成。另外,1个帧长度为30秒,具有1500bits的信息量。另外,在从最初的子帧1起第三个子帧3中,包含用于算出轨道长半径、离心率、平均近点角、升降点经度、近地点变量以及轨道倾斜角等要素的参数、和导航消息的发送时刻toc等广播星历信息。另一方面,在第四个子帧4及第五个子帧5中,包含在所有的人造卫星10中共用的预报星历(Almanac)信息。预报星历信息包含所有的人造卫星10的6个要素、和表示哪个人造卫星10可以使用的信息等。
另外,在各子帧中,在作为固定模式的前同步码之后记载数据。各子帧由10个字(word)构成,长度为6秒,具有300bits的信息量。
另外,各字由30个比特构成,长度为600msec。另外,各比特的长度是作为C/A码(扩展码)的20个周期的20msec。所以,数据的传输速度为50bps。另外,C/A码如图2所示,1个周期为1msec,由PRN码的1023个码片(chip)构成。
[1-3.接收机的构成及动作]
图3是表示本发明的实施方式涉及的接收机20的硬件构成的说明图。如图3所示,接收机20具备:包含天线212、频率变换部220、同步捕捉部240及同步保持部250的接收处理部210;CPU(CentralProcessing Unit)260;RTC(Real Time Clock,实时时钟)264;计时器268;存储器270;XO(石英晶体振荡器、X’tal Oscillator)272;TCXO(Temperature Compensated X’tal Oscillator,温补振荡器)274和倍增/分频器276。
XO272振荡出具有规定的频率(例如32.768kHz左右)的信号D1,将振荡出的信号D1向RTC264供给。TCXO274振荡出具有与XO272不同的频率(例如16.368MHz左右)的信号D2,将振荡出的信号D2向倍增/分频器276或频率合成器228供给。
倍增/分频器276对由TCXO274提供的信号D2基于来自CPU260的指示进行倍增、分频或其两者。而且,倍增/分频器276将进行了倍增、分频或其两者的信号D4向频率变换部220的频率合成器228、CPU260、计时器268、存储器270、同步捕捉部240及同步保持部250提供。
天线212接收从人造卫星10发送的导航消息等无线信号(例如1575.42MHz的载波被扩展的RF信号),将该无线信号变换为电信号D5并向频率变换部220提供。
频率变换部220具备:LNA(Low Noise Amplifier,低噪声放大器)222、BPF(Band Pass Filter,带通滤波器)224、放大器226、频率合成器228、乘法器230、放大器232、LPF(Low Pass Filter,低通滤波器)234和ADC(Analog Digital Converter,模数转换器)236。该频率变换部220如以下所示那样,为了容易进行数字信号处理,把利用天线212接收到的具有1575.42MHz的较高频率的信号D5降频变换为例如具有1.023MHz左右的频率的信号D14。
LNA222将从天线212提供的信号D5放大,向BPF224提供。BPF224由SAW滤波器(Surface Acoustic Wave Filter,声表面波滤波器),在由LNA放大的信号D6的频率成分中只提取特定的频率成分向放大器226提供。放大器226将具有利用BPF224提取的频率成分的信号D7(频率FRF)放大,向乘法器230提供。
频率合成器228利用由TCXO274提供的信号D2,基于来自CPU260的指示D9生成具有频率FLO的信号D10。而且,频率合成器228将所生成的具有频率FLO的信号D10向乘法器230提供。
乘法器230将由放大器226提供的具有频率FRF的信号D8与由频率合成器228提供的具有频率FLO的信号D10相乘。即,乘法器230将高频信号降频变换为IF(Intermediate Frequency,中间频率)信号D11(例如具有1.023MHz左右的频率的中间频率信号)。
放大器232将利用乘法器230降频变换后的IF信号D11放大,向LPF234提供。
LPF234在利用放大器232放大后的IF信号D12的频率成分中提取低频成分,将具有提取出的低频成分的信号D13向ADC236提供。而且,图3中虽然说明了在放大器232与ADC236之间配置LPF234的例子,然而也可以在放大器232与ADC236之间配置BPF。
ADC236(采样部)对由LPF234提供的模拟形式的IF信号D13进行采样,由此变换为数字形式,将变换为数字形式的IF信号D14向同步捕捉部240及同步保持部250提供。
同步捕捉部240根据CPU260的控制,利用由倍增/分频器276提供的信号D4,进行从ADC236提供的IF信号D14的PRN码的同步捕捉,另外,检测出IF信号D14的IF频率。而且,同步捕捉部240将PRN码的相位或IF信号D14的载波频率等向同步保持部250及CPU260提供。
同步保持部250根据CPU260的控制,利用由倍增/分频器276提供的信号D4,进行从ADC236提供的IF信号D14的PRN码与载波的同步保持。更具体来说,同步保持部250将由同步捕捉部240提供的PRN码的相位或IF信号D14的载波频率作为初始值动作。另外,同步保持部250对由ADC236提供的IF信号D14中所包含的导航消息进行解调,向CPU260提供。而且,同步保持部250具备多个信道用电路,每个信道用电路进行从对应的人造卫星10发送的信号的同步保持及导航消息的解调。
CPU260基于从同步保持部250提供的导航消息,计算出各人造卫星10的位置或速度,计算接收机20的位置。另外,CPU260基于导航消息进行RTC264的时间信息的修正,或与控制端子、I/O端子及附加功能端子等连接而进行各种控制。
RTC264利用从XO272提供的具有规定的频率的信号D1计测时间。由CPU260对RTC264计测的时间适当地进行修正。
计时器268利用从倍增/分频器276提供的信号D4进行计时。这样的计时器268在决定基于CPU260的各种控制的开始定时等时被参照。例如,CPU260在基于由同步捕捉部240捕捉到的PRN码的相位决定开始同步保持部250的PRN码发生器的动作的定时时,参照计时器268。
存储器270由RAM(Random Access Memory)或ROM(Read-OnlyMemory)等构成,具有作为CPU260的操作空间、程序的存储部、导航消息的存储部等的功能。在存储器270中,作为进行基于CPU260等的各种处理时的工作区使用RAM。另外,在将所输入的各种数据进行缓冲时、或为了保持利用同步保持部250得到的作为人造卫星10的轨道信息的广播星历及预报星历、在运算过程中生成的中间数据及运算结果数据,也使用RAM。另外,在存储器270中,作为存储各种程序或固定数据等的单元使用ROM。另外,在存储器270中,有时使用非易失性存储器,作为即使在接收机20的电源被切断的期间,也能存储作为人造卫星10的轨道信息的广播星历及预报星历、以及定位结果的位置信息、TCXO12的误差量等的单元。
而且,也可以在由1个芯片构成的集成电路中安装图3所示的接收机20的构成中的除了XO272、TCXO274、天线212、LNA222及BPF224以外的构成。
接下来,参照图4,对上述的接收机20的动作例进行说明。
图4是表示本发明的实施方式涉及的接收机20的动作例的流向的流程图。如图4所示,若起动了接收机20,CPU260就进行初始设定(S42)。接下来,当利用RTC264计数1秒时(S44),则CPU260就对同步保持部250中所包含的各信道用电路分配人造卫星10(S46)。
其后,当利用接收处理部210取得导航消息时(S48),则CPU260就选择实际上保持的至少4个以上的人造卫星10(S50)。而且,CPU260计算出所选择的人造卫星10的当前位置及速度(S52),基于计算出的人造卫星10的当前位置及速度计算接收机20的当前位置及速度(S54)。
接下来,CPU260制成表示计算出的接收机20的当前位置及速度的输出消息(S56),在执行了与该输出消息相应的指令处理后返回到S44的处理(S58)。
<2.本发明的实施方式的背景>
下面,参照图5,对本发明的实施方式的背景进行说明。
图5是表示与本发明相关的同步捕捉部240’的构成的说明图。如图5所示,与本发明相关的同步捕捉部240’具备:存储器32、FFT34、存储器364、读出地址控制部38、PRN码产生部40、FFT42、存储器44、乘法器46、反FFT48以及相关点检测部50。
如图5所示,由ADC236提供的IF信号D14被记录于存储器32中,FFT34对记录于存储器32中的IF信号D14进行高速傅立叶变换。FFT34的高速傅立叶变换结果被记录于存储器36中。记录于存储器36中的高速傅立叶变换结果被基于读出地址控制部38的控制,依次向乘法器46提供。
另一方面,PRN码产生部40产生在接收信号的扩展中被利用的PRN码,FFT42对由PRN码产生部40产生的PRN码进行高速傅立叶变换。FFT42的高速傅立叶变换结果被记录于存储器44中。
乘法器46将由存储器36提供的高速傅立叶变换结果、与由存储器44提供的高速傅立叶变换结果的复共轭相乘。即,乘法器46运算频率域中的接收信号的PRN码与由PRN码产生部40产生的PRN码之间的相关程度。
反FFT48对由乘法器46提供的乘法结果进行反高速傅立叶变换。反FFT48的反高速傅立叶变换结果表示时间域中的接收信号的PRN码与由PRN码产生部40产生的PRN码之间的相关性。相关点检测部50基于该反高速傅立叶变换结果,检测接收信号的PRN码与由PRN码产生部40产生的PRN码同步的相关点。
这里,FFT34也可以通过将用于获得全部频率的1/L(L为整数)的频率成分的高速傅立叶变换进行L次来获得全部频率成分。该情况下,能够通过在进行新的高速傅立叶变换期间由乘法器46乘以记录于存储器36中的上次的高速傅立叶变换结果,来利用新的高速傅立叶变换结果将存储器36更新。所以,能够通过分割进行高速傅立叶变换,来削减存储器36的容量。
但是,在FFT34进行的L次的高速傅立叶变换中,包含很多重复的运算。所以,在FFT34单独地分割进行高速傅立叶变换时,存在运算量及运算时间增加的问题。
所以,以上述情况作为一个着眼点而制造出本发明的实施方式涉及的接收机20。本发明的实施方式涉及的接收机20可以缩短用于高速傅立叶变换的运算时间。以下,对此种接收机20的第一实施方式~第五实施方式进行说明。
<3.第一实施方式>
实际的接收机20例如由ADC236以4.096MHz进行采样,同步捕捉部240将由ADC236得到的数msec的采样值作为输入数据进行高速傅立叶变换。但是,下面为了说明的方便,对利用16个输入数据进行高速傅立叶变换的例子进行说明。
图6是表示输入数据为16个时的高速傅立叶变换的运算的流向的说明图。时间域的输入数据x0~x15被通过图6所示的运算,变换为频率域的输出数据X0~X15。而且,图6中的“W”是“exp(-j2π/16)”。
另一方面,如图7~图10所示,也可以通过将高速傅立叶变换分割为4次进行,来得到输出数据X0~X15。
图7~图10是表示分割进行高速傅立叶变换时的运算的流向的说明图。具体来说,图7表示用于获得输出数据X0、X8、X4及X12的运算,图8表示用于获得输出数据X2、X10、X6及X14的运算。同样地,图9表示用于获得输出数据X1、X9、X5及X13的运算,图10表示用于获得输出数据X3、X11、X7及X15的运算。
这里,如果将运算过程中得到的各中间数据如图11所示地定义,则成为针对用于获得图8所示的输出数据X2、X10、X6及X14的高速傅立叶变换的输入的中间数据x”4~x”7可以如下所示地表现。(W^0=W^16=1、W^8=-1)
数学式1:
x”4=(x’0-x’4)×w^0
=((x0+x8)-(x4+x12))×w^0
=x0×w^0+x4×w^8+x8×w^0+x12×w^8
x”5=(x’1-x’5)×w^2
=((x1+x9)-(x5+x13))×w^2
=x1×w^2+x5×w^10+x9×w^2+x13×w^10
x”6=(x’2-x’6)×w^4
=((x2+x10)-(x6+x14))×w^4
=x2×w^4+x6×w^12+x10×w^4+x14×w^12
x”7=(x’3-x’7)×w^6
=((x3+x11)-(x7+x15))×w^6
=x3×w^6+x7×w^14+x11×w^6+x15×w^14
此种用于获得中间数据x”4~x”7的运算与如下的处理等价,即,将输入数据与相位按每个输入数据变化W^2的信号相乘,将乘积值每隔4个累积。而且,由于乘积值的累积间隔为输入数据数/分割数,因此在输入数据数为16个、分割数为4的本实施方式中,如上所述成为4个。
同样地,成为针对用于获得图7所示的输出数据X0、X8、X4及X12的高速傅立叶变换的输入的中间数据x”0~x”3可以如以下的数学式2所示地表现。
数学式2
x”0=(x’0+x’4)
=(x0+x8)+(x4+x12)
=x0×w^0+x4×w^0+x8×w^0+x12×w^0
x”1=(x’1+x’5)
=(x1+x9)+(x5+x13)
=x1×w^0+x5×w^0+x9×w^0+x13×w^0
x”2=(x’2+x’6)
=(x2+x10)+(x6+x14)
=x2×w^0+x6×w^0+x10×w^0+x14×w^0
x”3=(x’3+x’7)
=(x3+x11)+(x7+x15)
=x3×w^0+x7×w^0+x11×w^0+x15×w^0
此种用于获得中间数据x”0~x”3的运算与如下的处理等价,即,将输入数据与相位不变化的信号(或者相位按每个输入数据变化W^16的信号)相乘,将乘积值每隔4个累积。
另外,成为针对用于获得输出数据X2、X10、X6及X14的高速傅立叶变换的输入的中间数据x”8~x”11可以如以下的数学式3所示地表现。
数学式3
x”8=(x’8+x’12)
=(x0-x8)×w^0+(x4-x12)×w^4
=x0×w^0+x4×w^4+x8×w^8+x12×w^12
x”9=(x’9+x’13)
=(x1-x9)×w^1+(x5-x13)×w^5
=x1×w^1+x5×w^5+x9×w^9+x13×w^13
x”10=(x’10+x’14)
=(x2-x10)×w^2+(x6-x14)×w^6
=x2×w^2+x6×w^6+x10×w^10+x14×w^14
x”11=(x’11+x’15)
=(x3-x11)×w^3+(x7-x15)×w^7
=x3×w^3+x7×w^7+x11×w^11+x15×w^15
此种用于获得中间数据x”8~x”11的运算与如下的处理等价,即,将输入数据与相位按每个输入数据变化W^1的信号相乘,将乘积值每隔4个累积。
另外,成为针对用于获得输出数据X3、X11、X7及X15的高速傅立叶变换的输入的中间数据x”12~x”15可以如以下的数学式4所示地表现。
数学式4
x”12=(x’8-x’12)×w^0
=((x0-x8)×w^0-(x4-x12)×w^4)×w^0
=x0×w^0+x4×w^12+x8×w^8+x12×w^4
x”13=(x’9-x’13)×w^2
=((x1-x9)×w^1-(x5-x13)×w^5)×w^2
=x1×w^3+x5×w^15+x9×w^11+x13×w^7
x”14=(x’10-x’14)×w^4
=((x2-x10)×w^2-(x6-x14)×w^6)×w^4
=x2×w^6+x6×w^2+x10×w^14+x14×w^10
x”15=(x’11-x’15)×w^6
=((x3-x11)×w^3-(x7-x15)×w^7)×w^6
=x3×w^9+x7×w^5+x11×w^1+x15×w^13
此种用于获得中间数据x”12~x”15的运算与如下的处理等价,即,将输入数据与相位按每个输入数据变化W^3的信号相乘,将乘积值每隔4个累积。
如上说明所示,在分割进行高速傅立叶变换的情况下,可以通过规定的前处理获得成为高速傅立叶变换的对象的中间数据x”0~x”15。本实施方式涉及的同步捕捉部240,取得从ADC126作为IF信号D14提供的输入数据的同时,进行规定的前处理,由此就可以缩短用于获得全部频率成分的高速傅立叶变换结果的处理时间。以下,对此种同步捕捉部240的详细的构成及动作进行说明。
[3-1同步捕捉部的构成]
图12是表示第一实施方式涉及的同步捕捉部240的构成的说明图。如图12所示,第一实施方式涉及的同步捕捉部240具备:前处理执行部310、存储器320、FFT330、存储器332、读出地址控制部334、运算控制部336、PRN码产生部340、FFT344(高速傅立叶变换部)、存储器348、乘法器352、反FFT356和相关点检测部360。
另外,前处理执行部310具备信号生成部312、乘法部314、加法部316和存储器318。该前处理执行部310执行为了分割进行高速傅立叶变换而进行的前处理。
信号生成部312生成具有由运算控制部336指示的周期的周期信号。此外,乘法部314将作为输入数据由ADC236提供的IF信号D14与由信号生成部312生成的周期信号相乘。即,乘法部314将输入数据与相位按每个数据(每个采样间隔)变化规定量的信号相乘。
加法部316及存储器318作为将利用乘法部314得到的乘积值依照由乘法部314得到的顺序以规定间隔进行累积的累积部发挥功能。具体来说,加法部316将利用乘法部314得到的乘法值与从存储器318的某个地址中读出的累积值相加,利用相加结果将上述地址的累积值更新,对存储器318的地址1个1个错位地反复进行该处理。
例如,在将存储器318以可以存储4个数据的方式构成,信号生成部312生成相位按每个输入数据变化W^2的周期信号,输入数据数为16的情况下,就可以将数学式1中所示的中间数据x”4~x”7作为累积值获得。
同样地,在将存储器318以可以存储4个数据的方式构成,信号生成部312生成相位不按每个输入数据变化的周期信号,输入数据数为16的情况下,就可以将数学式2中所示的中间数据x”0~x”3作为累积值获得。另外,在将存储器318以可以存储4个数据的方式构成,信号生成部312生成相位按每个输入数据变化W^1的周期信号,输入数据数为16的情况下,就可以将数学式3中所示的中间数据x”8~x”11作为累积值获得。此外,在将存储器318以可以存储4个数据的方式构成,信号生成部312生成相位按每个输入数据变化W^3的周期信号,输入数据数为16的情况下,就可以将数学式4中所示的中间数据x”12~x”15作为累积值获得。
即,通过每当输入数据数达到16时,运算控制部336就使信号生成部312中生成的周期信号的周期变化,将存储器318复位,由此就可以依次得到中间数据x”0~x”3、x”4~x”7、x”8~x”11、x”12~x”15。
在存储器320中,依次记录利用前处理执行部310得到的中间数据x”0~x”3、x”4~x”7、x”8~x”11或x”12~x”15。FFT(Fast FourierTransform)330对记录于存储器320中的中间数据x”0~x”3、x”4~x”7、x”8~x”11或x”12~x”15的任意一个进行高速傅立叶变换。FFT330的高速傅立叶变换结果被记录于存储器332中。记录于存储器332中的高速傅立叶变换结果被基于读出地址控制部334的控制,依次向乘法器352提供。
另一方面,PRN码产生部340产生在接收信号的扩展中所利用的PRN码,FFT344对由PRN码产生部340产生的PRN码进行高速傅立叶变换。FFT344的高速傅立叶变换结果被记录于存储器348中。
乘法器352将由存储器332提供的高速傅立叶变换结果、与由存储器348提供的高速傅立叶变换结果的复共轭相乘。即,乘法器352运算频率域中的接收信号的PRN码与由PRN码产生部40产生的PRN码之间的相关程度。
反FFT356对由乘法器352提供的乘法结果进行反高速傅立叶变换。反FFT356的反高速傅立叶变换结果表示时间域中的接收信号的PRN码与由PRN码产生部340产生的PRN码之间的相关性。相关点检测部360基于该反高速傅立叶变换结果,检测接收信号的PRN码与由PRN码产生部340产生的PRN码同步的相关点。
[3-2同步捕捉部的动作]
接下来,参照图13,对本实施方式涉及的同步捕捉部240的动作的流向进行说明。而且,下面为了说明的方便,将中间数据x”0~x”3称作中间数据A,将中间数据x”4~x”7称作中间数据B,将中间数据x”8~x”11称作中间数据C,将中间数据x”12~x”15称作中间数据D。
图13是表示第一实施方式涉及的同步捕捉部240的动作的流向的说明图。如图13所示,前处理执行部310基于运算控制部336的控制,依次进行用于获得中间数据A、B、C及D的前处理。
另外,结束用于获得各个中间数据A~D的前处理执行部310的各前处理的同时,向存储器320转送通过前处理得到的中间数据。而且,当中间数据的存储器转送结束时,FFT330就对转送到存储器320的中间数据进行高速傅立叶变换。这里,由于在利用FFT330进行中间数据的高速傅立叶变换的期间,前处理执行部310进行用于获得新的中间数据的前处理,因此可以缩短用于获得全部频率的输出数据X0~X15的处理时间。
<4.第二实施方式>
如上说明的那样,第一实施方式涉及的同步捕捉部240具备前处理执行部310及存储器320各一个。与之不同,第二实施方式涉及的同步捕捉部240具备前处理执行部310及存储器320各两个,从而可以进一步缩短处理时间。
[4-1同步捕捉部的构成]
图14是表示第二实施方式涉及的同步捕捉部241的构成的说明图。如图14所示,第二实施方式涉及的同步捕捉部241具备:前处理执行部310A及310B、存储器320A及320B、FFT330、存储器332、读出地址控制部334、运算控制部336、PRN码产生部340、FFT344、存储器348、乘法器352、反FFT356和相关点检测部360。而且,由于PRN码产生部340、FFT344、存储器348、乘法器352、反FFT356及相关点检测部360等可以与第一实施方式的同步捕捉部240实质上相同地构成,因此省略说明。
本实施方式涉及的前处理执行部310A及310B基于运算控制部336的控制,生成具有不同的周期的周期信号。由此,从各个前处理执行部310A及310B中,得到不同的中间数据。
在存储器320A中记录利用前处理执行部310A得到的中间数据,在存储器320B中记录利用前处理执行部310B得到的中间数据。
[4-2同步捕捉部的动作]
图15是表示第二实施方式涉及的同步捕捉部241的动作的流向的说明图。如图15所示那样,前处理执行部310A基于运算控制部336的控制,交替地进行用于获得中间数据A的前处理、和用于获得中间数据C的前处理。另外,前处理执行部310B基于运算控制部336的控制,交替地进行用于获得中间数据B的前处理、和用于获得中间数据D的前处理。
而且,结束基于前处理执行部310A的各前处理的同时,向存储器320A转送利用前处理执行部310A得到的中间数据。进而,当中间数据的存储器转送结束时,FFT330就对转送到存储器320A的中间数据进行高速傅立叶变换。
同样地,结束基于前处理执行部310B的各前处理的同时,向存储器320B转送利用前处理执行部310B得到的中间数据。进而,当中间数据的存储器转送和存储器320A的中间数据的高速傅立叶变换结束时,FFT330就对转送到存储器320B的中间数据进行高速傅立叶变换。
这样,第二实施方式涉及的同步捕捉部241,由于具备多个前处理执行部310及存储器320,因此可以并列地获得不同的中间数据。其结果是,根据第二实施方式,与第一实施方式相比可以进一步缩短处理时间。
而且,虽然在图15中,表示了基于前处理执行部310A的输入数据的取得开始定时(图中的上升边)与基于前处理执行部310B的输入数据的取得开始定时一致的例子,然而两者的定时也可以不同。
<5.第三实施方式>
如上说明的那样,第一实施方式涉及的同步捕捉部240及第二实施方式涉及的同步捕捉部241具备与前处理执行部310相同数目的存储器320。与之不同,第三实施方式涉及的同步捕捉部242中,由于可以使存储器320的数目比前处理执行部310少,因此可以实现电路规模的削减。
[5-1同步捕捉部的构成]
图16是表示第三实施方式涉及的同步捕捉部242的构成的说明图。如图16所示那样,第三实施方式涉及的同步捕捉部242具备:前处理执行部310A~310D、存储器320、FFT330、存储器332、读出地址控制部334、运算控制部336、PRN码产生部340、FFT344、存储器348、乘法器352、反FFT356和相关点检测部360。而且,由于PRN码产生部340、FFT344、存储器348、乘法器352、反FFT356及相关点检测部360等可以与第一实施方式涉及的同步捕捉部240实质上相同地构成,因此省略说明。
每个前处理执行部310A~310D进行用于获得不同的中间数据的前处理。利用前处理执行部310A~310D得到的中间数据的每个依次被向存储器320转送。FFT330对向存储器320转送的中间数据进行高速傅立叶变换。
[5-2同步捕捉部的动作]
图17是表示第三实施方式涉及的同步捕捉部242的动作的流向的说明图。如图17所示那样,前处理执行部310A进行用于获得中间数据A的前处理,前处理执行部310B进行用于获得中间数据B的前处理。同样地,前处理执行部310C进行用于获得中间数据C的前处理,前处理执行部310D进行用于获得中间数据D的前处理。
而且,结束基于前处理执行部310A的前处理的同时,向存储器320转送利用前处理执行部310A得到的中间数据A。进而,当中间数据A的存储器转送结束时,FFT330就对转送到存储器320的中间数据A进行高速傅立叶变换。
其后,利用前处理执行部310B得到的中间数据B被向存储器320转送,FFT330对转送到存储器320的中间数据B进行高速傅立叶变换。接下来,利用前处理执行部310C得到的中间数据C被向存储器320转送,FFT330对转送到存储器320的中间数据C进行高速傅立叶变换。进而,利用前处理执行部310D得到的中间数据D被向存储器320转送,FFT330对转送到存储器320的中间数据D进行高速傅立叶变换。
这样,根据本实施方式,利用不同的前处理执行部310A~310D得到的中间数据A~D被依次向相同的存储器320转送。所以,由于在本实施方式中不需要设置与前处理执行部310相同数目的存储器320,因此可以实现电路规模的削减。
<6.第四实施方式>
另外,通过如以下说明的第四实施方式涉及的同步捕捉部243那样,设置与分割数相同数目的前处理执行部310及存储器320,就可以进一步缩短用于得到全部频率的输出数据X0~X15的处理时间。
[6-1同步捕捉部的构成]
图18是表示第四实施方式涉及的同步捕捉部243的构成的说明图。如图18所示那样,第四实施方式涉及的同步捕捉部243具备:前处理执行部310A~310D、存储器320A~320D、FFT330、存储器332、读出地址控制部334、运算控制部336、PRN码产生部340、FFT344、存储器348、乘法器352、反FFT356和相关点检测部360。而且,由于PRN码产生部340、FFT344、存储器348、乘法器352、反FFT356及相关点检测部360等可以与第一实施方式涉及的同步捕捉部240实质上相同地构成,因此省略说明。
前处理执行部310A~310D的每个进行用于获得不同的中间数据的前处理。利用前处理执行部310A得到的中间数据被向存储器320A转送,利用前处理执行部310B得到的中间数据被向存储器320B转送。同样,利用前处理执行部310C得到的中间数据被向存储器320C转送,利用前处理执行部310D得到的中间数据被向存储器320D转送。
FFT330对向存储器320A~存储器320D转送的中间数据依次进行高速傅立叶变换。
[6-2同步捕捉部的动作]
图19是表示第四实施方式涉及的同步捕捉部243的动作的流向的说明图。如图19所示那样,前处理执行部310A进行用于获得中间数据A的前处理,前处理执行部310B进行用于获得中间数据B的前处理。同样地,前处理执行部310C进行用于获得中间数据C的前处理,前处理执行部310D进行用于获得中间数据D的前处理。
而且,结束基于前处理执行部310A的前处理的同时,向存储器320A转送利用前处理执行部310A得到的中间数据A。另外,结束基于前处理执行部310B的前处理的同时,向存储器320B转送利用前处理执行部310B得到的中间数据B。另外,结束基于前处理执行部310C的前处理的同时,向存储器320C转送利用前处理执行部310C得到的中间数据C。同样地,结束基于前处理执行部310D的前处理的同时,向存储器320D转送利用前处理执行部310D得到的中间数据D。
其后,FFT330对向存储器320A~存储器320D转送的中间数据依次进行高速傅立叶变换。
这样,根据本实施方式,由于可以并列地进行用于获得各中间数据的前处理和各中间数据的存储器转送,因此可以实现用于获得全部频率的输出数据X0~X15的处理时间的缩短。
<7.第五实施方式>
在上述说明中,为了说明的方便,对将16个输入数据用于高速傅立叶变换的例子进行了说明,然而在实际上,是将数msec的输入数据用于高速傅立叶变换。这里,如参照图2所说明的那样,构成导航消息的各比特的长度为20msec。所以,可以设想如下的情况,即,在要利用的输入数据的途中比特值被切换,从而无法正常地检测出PRN码的相关点。第五实施方式涉及的同步捕捉部244可以防止此种情况。
[7-1同步捕捉部的构成]
图20是表示第五实施方式涉及的同步捕捉部244的构成的说明图。如图20所示那样,第五实施方式涉及的同步捕捉部244具备:前处理执行部310A及310B、存储器320、FFT330、存储器332、读出地址控制部334、运算控制部336、PRN码产生部340、FFT344、存储器348、乘法器352、反FFT356、相关点检测部360和选择部364。而且,由于PRN码产生部340、FFT344、存储器348、乘法器352、反FFT356及相关点检测部360等可以与第一实施方式涉及的同步捕捉部240实质上相同地构成,因此省略说明。
前处理执行部310A及310B在不同的定时开始输入数据的取得,进行用于基于具有相同周期的周期信号获得中间数据的前处理。利用前处理执行部310A得到的中间数据及利用前处理执行部310B得到的中间数据被向存储器320转送。FFT330对向存储器320转送的中间数据依次进行高速傅立叶变换。FFT330的高速傅立叶变换结果被向存储器332提供。
乘法器352将由存储器332提供的高速傅立叶变换结果与由存储器348提供的高速傅立叶变换结果的复共轭相乘。反FFT356对由乘法器352提供的乘法结果进行反高速傅立叶变换。反FFT356的反高速傅立叶变换结果表示时间域中的接收信号的PRN码与利用PRN码产生部340产生的PRN码之间的相关性。相关点检测部360基于该反高速傅立叶变换结果,检测接收信号的PRN码与利用PRN码产生部340产生的PRN码同步的相关点。
选择部364比较来自相关点检测部360的基于利用前处理执行部310A得到的中间数据的输出、和基于利用前处理执行部310B得到的中间数据的输出,将良好的一方的输出作为相关点选择。
[7-2同步捕捉部的动作]
图21是表示第五实施方式涉及的同步捕捉部244的动作的流向的说明图。如图21所示那样,前处理执行部310A基于运算控制部336的控制,依次进行用于获得中间数据A、B、C及D的前处理。另外,前处理执行部310B基于运算控制部336的控制,在与前处理执行部310A不同的定时开始输入数据的取得,依次进行用于获得中间数据A’、B’、C’及D’的前处理。而且,中间数据A’是输入数据的取得开始定时与中间数据A不同,但是基于与中间数据A相同的周期信号得到的数据。中间数据B’、C’及D’也相同。
而且,结束基于前处理执行部310A的各前处理的同时,向存储器320转送利用前处理执行部310A得到的中间数据。另外,结束基于前处理执行部310B的各前处理的同时,向存储器320转送利用前处理执行部310B得到的中间数据。
FFT330对向存储器320转送的中间数据依次进行高速傅立叶变换。其后,在相关点检测部360中,进行基于FFT330的中间数据的高速傅立叶变换结果和由存储器348提供的高速傅立叶变换结果的相关点的检测。而且,选择部364选择来自相关点检测部360的基于中间数据A~D的高速傅立叶变换结果的输出和基于中间数据A’~D’的高速傅立叶变换结果的输出中的任意一个。所以,根据本实施方式,可以防止因输入数据的途中的比特值的切换而无法正常地检测出PRN码的相关点的情况。
<8.本发明的实施方式的效果>
(1)根据本发明的实施方式,如图13、图15等所示那样,由于可以同时地进行用于得到中间数据的基于前处理执行部310的前处理、和基于FFT330的高速傅立叶变换,因此可以实现处理时间的缩短。
(2)另外,根据本发明的实施方式,可以削减配置于FFT330的前级的存储器320所需要的容量。以下,说明其理由。
在与本发明相关的图5所示的同步捕捉部240’中,利用ADC236变换为数字形式的数据被原样不动地记录于存储器32中。由此,在采样频率为N(kHz)、比特数为I(Bit)、输入数据的长度为M(ms)的情况下,存储器32的大小就需要N×M×I(Bit)。
但是,输入数据的能量的大半为热噪音,可以看作是遵循高斯分布的随机数。由此,通过如本发明的实施方式那样进行输入数据的累积,就可以削减存储器320所需要的容量。
例如,可以考虑将高速傅立叶变换分割为16次进行的情况。另外,将输入数据的标准偏差设为σ0,将前处理执行部310中的通过16次的加法运算得到的中间数据的标准偏差设为σ1。并且,如果将通过由ADC236进行适当的阈值处理而以I(Bit)表现的最大值Max0表现为3×σ0,则以下的数学式成立。
σ1=√16×σ0=4×σ0(根据方差的可加性)
另外,如果将中间数据的最大值Max1表现为3×σ1,则以下的数学式成立。
Max1=3×σ1=12×σ0=4×Max0
根据上述式子,由于Max1是Max0的4倍,因此只要使用于存储各中间数据的比特数比存储各输入数据时增加2比特即可。所以,在与本发明相关的同步捕捉部240’的存储器32中需要N×M×I(Bit)的容量,而在本发明的实施方式的存储器320中则为N×M÷16×(I+2)(Bit)的容量,由于I为整数,因此可以大幅度削减存储器320的容量。而且,即使在将高速傅立叶变换分割为4次进行的情况下,由于存储器320需要的容量为N×M÷4×(I+1)(Bit),因此也可以削减存储器320的容量。
(3)通常来说,从FFT的性质考虑,向FFT输入的数据数需要为2的幂。所以,在以1.024MHz的倍数对1023chip/ms的接收信号进行采样时,以往只能将8ms的采样值、或16ms的采样值等作为针对FFT的输入数据而利用。
与之不同,本发明的实施方式中,虽然向FFT330输入的中间数据的数需要是2的幂,然而前处理执行部310中的加法次数只要是2以上的整数,就没有特别限制。所以,通过适当地设定前处理执行部310中的加法次数,即高速傅立叶变换的分割数,就可以将任意的整数ms的采样值作为用于高速傅立叶变换的输入数据而利用。
例如,通过将前处理执行部310的加法次数,即,高速傅立叶变换的分割数设定为5,就可以将5ms的采样值作为用于高速傅立叶变换的输入数据而利用。
(4)如图18所示那样,第四实施方式涉及的同步捕捉部243具备多个前处理执行部310和针对FFT330的输入用的存储器320。由此,在与前处理相比高速地进行高速傅立叶变换的情况下,如图19所示那样,前处理执行部310一旦结束了向存储器320的中间数据的转送,就可以开始下面的输入数据的取得和前处理。
(5)另外,根据第五实施方式,多个前处理执行部310在不同的定时开始输入数据的取得,进行用于基于具有相同周期的周期信号取得中间数据的前处理。并且,选择部364比较基于利用各前处理执行部310得到的中间数据的来自相关点检测部360的输出,选择良好的一方的输出。所以,根据第五实施方式,可以防止因输入数据的途中的比特值的切换而无法正常地检测出PRN码的相关点的情况。
[补充]
以上虽然在参照附图的同时对本发明的优选的实施方式进行了详细说明,然而本发明并不限定于这样的例子。显而易见,只要是具有本发明所属的技术领域中的通常的知识的人,就可以在技术方案的范围中记载的技术思想的范畴内,想到各种变更例或修正例,可以理解,它们当然也属于本发明的技术范围。
Claims (9)
1.一种通信装置,其特征在于,具备:
采样部,其对接收信号进行采样;
前处理执行部,其将利用上述采样部得到的采样值一边实施为了分割进行高速傅立叶变换而进行的前处理一边取得;和
高速傅立叶变换部,其基于来自上述前处理执行部的输出,分割进行高速傅立叶变换。
2.根据权利要求1所述的通信装置,其特征在于,
所述前处理执行部具备:
信号生成部,其生成周期信号;
乘法部,其将利用上述采样部得到的采样值和利用上述信号生成部生成的上述周期信号的值相乘;和
累积部,其将利用上述乘法部得到的各个乘积值依照利用上述乘法部获得的顺序以规定间隔累积,得到多个累积值,
其中,上述高速傅立叶变换部对利用上述累积部得到的上述多个累积值进行高速傅立叶变换。
3.根据权利要求2所述的通信装置,其特征在于,
上述通信装置具备多个上述前处理执行部,
每个上述前处理执行部所包含的上述信号生成部生成具有不同的周期的周期信号。
4.根据权利要求3所述的通信装置,其特征在于,
每个上述前处理执行部的上述累积部以接收信号的规定长度为单位进行各个上述乘积值的累积,
上述通信装置还具备与每个上述前处理执行部对应的存储器,
在上述存储器中,记录利用对应的前处理执行部得到的上述多个累积值,
上述高速傅立叶变换部将依次记录于每个上述存储器中的上述多个累积值进行高速傅立叶变换。
5.根据权利要求4所述的通信装置,其特征在于,
上述接收信号被利用扩展码进行了频谱扩展,
上述通信装置还具备:第二前处理执行部,其包括生成与第一前处理执行部的上述信号生成部具有相同的周期的周期信号的信号生成部,接收信号被以上述规定长度为单位划分的位置与上述第一前处理执行部中不同;
相关点检测部,其基于上述高速傅立叶变换部的高速傅立叶变换结果,检测上述扩展码的相关点;和
选择部,其选择基于利用上述第一前处理执行部得到的上述多个累积值的高速傅立叶变换结果由相关点检测部检测出的相关点和基于利用上述第二前处理执行部得到的上述多个累积值的高速傅立叶变换结果由相关点检测部检测出的相关点中的任意一个。
6.根据权利要求2所述的通信装置,其特征在于,
上述累积部以接收信号的规定长度为单位进行各个上述乘积值的累积,
上述信号生成部以接收信号的上述规定长度为单位变更要生成的周期信号的周期。
7.一种信号处理方法,其特征在于,包括:
对接收信号进行采样的步骤;
将通过采样得到的采样值一边实施为了分割进行高速傅立叶变换而进行的前处理一边记录于存储器中的步骤;和
基于来自上述存储器的输出分割进行高速傅立叶变换的步骤。
8.一种信号处理装置,其特征在于,具备:
采样部,其对接收信号进行采样;
前处理执行部,其将利用上述采样部得到的采样值一边实施为了分割进行高速傅立叶变换而进行的前处理一边取得;和
高速傅立叶变换部,其基于来自上述前处理执行部的输出分割进行高速傅立叶变换。
9.一种移动体,其特征在于,
搭载有如下的通信装置,即,该通信装置具备:
采样部,其对接收信号进行采样;
前处理执行部,其将利用上述采样部得到的采样值一边实施为了分割进行高速傅立叶变换而进行的前处理一边取得;和
高速傅立叶变换部,其基于来自上述前处理执行部的输出分割进行高速傅立叶变换。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2009-096195 | 2009-04-10 | ||
JP2009096195A JP5353396B2 (ja) | 2009-04-10 | 2009-04-10 | 通信装置、信号処理方法、信号処理装置、および移動体 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101859292A true CN101859292A (zh) | 2010-10-13 |
CN101859292B CN101859292B (zh) | 2013-01-02 |
Family
ID=42543065
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2010101512653A Expired - Fee Related CN101859292B (zh) | 2009-04-10 | 2010-04-02 | 通信装置、信号处理方法、信号处理装置以及移动体 |
Country Status (6)
Country | Link |
---|---|
US (1) | US8249202B2 (zh) |
EP (1) | EP2239593B1 (zh) |
JP (1) | JP5353396B2 (zh) |
CN (1) | CN101859292B (zh) |
AT (1) | ATE549639T1 (zh) |
TW (1) | TWI425242B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2022077938A1 (zh) * | 2020-10-15 | 2022-04-21 | 华为技术有限公司 | 一种星历数据传输方法、装置、存储介质、芯片及系统 |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20220007420A1 (en) * | 2020-07-01 | 2022-01-06 | Qualcomm Incorporated | Spreading aspects of random access channel procedure |
US20240094410A1 (en) * | 2022-09-20 | 2024-03-21 | Onenav, Inc. | Methods and apparatuses for improved acquisition of modern gnss signals |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2003001698A1 (fr) * | 2001-06-25 | 2003-01-03 | Sony Corporation | Procede de demodulation de signaux d'etalement du spectre et appareil correspondant |
US20030187576A1 (en) * | 2002-02-28 | 2003-10-02 | Jason Sanmiya | GPS receiver and receiving method |
CN101093254A (zh) * | 1995-10-09 | 2007-12-26 | 快速追踪有限公司 | 处理gps信号的gps接收机及方法 |
US20090013021A1 (en) * | 2007-07-06 | 2009-01-08 | Mediatek Inc. | Variable length fft system and method |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
AU2001281180A1 (en) * | 2000-08-08 | 2002-02-18 | Enuvis, Inc. | Methods and apparatus for dynamically determining location using a simplified client device |
JP3906913B2 (ja) | 2001-06-25 | 2007-04-18 | ソニー株式会社 | スペクトラム拡散信号復調方法および装置 |
JP2003316763A (ja) * | 2001-07-04 | 2003-11-07 | Sony Corp | 周波数分析方法および装置並びにスペクトラム拡散復調方法および装置 |
US7639181B2 (en) * | 2005-07-01 | 2009-12-29 | Sirf Technology Holdings, Inc. | Method and device for tracking weak global navigation satellite system (GNSS) signals |
US8031816B2 (en) * | 2006-07-17 | 2011-10-04 | Mediatek Inc. | Method and apparatus for determining boundaries of information elements |
JP2009096195A (ja) | 2007-09-27 | 2009-05-07 | Fujifilm Corp | 感熱転写受像シート |
-
2009
- 2009-04-10 JP JP2009096195A patent/JP5353396B2/ja not_active Expired - Fee Related
-
2010
- 2010-02-23 US US12/710,624 patent/US8249202B2/en not_active Expired - Fee Related
- 2010-03-03 TW TW099106184A patent/TWI425242B/zh not_active IP Right Cessation
- 2010-03-19 EP EP10157016A patent/EP2239593B1/en not_active Not-in-force
- 2010-03-19 AT AT10157016T patent/ATE549639T1/de active
- 2010-04-02 CN CN2010101512653A patent/CN101859292B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101093254A (zh) * | 1995-10-09 | 2007-12-26 | 快速追踪有限公司 | 处理gps信号的gps接收机及方法 |
WO2003001698A1 (fr) * | 2001-06-25 | 2003-01-03 | Sony Corporation | Procede de demodulation de signaux d'etalement du spectre et appareil correspondant |
US20030187576A1 (en) * | 2002-02-28 | 2003-10-02 | Jason Sanmiya | GPS receiver and receiving method |
US20090013021A1 (en) * | 2007-07-06 | 2009-01-08 | Mediatek Inc. | Variable length fft system and method |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2022077938A1 (zh) * | 2020-10-15 | 2022-04-21 | 华为技术有限公司 | 一种星历数据传输方法、装置、存储介质、芯片及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN101859292B (zh) | 2013-01-02 |
ATE549639T1 (de) | 2012-03-15 |
US8249202B2 (en) | 2012-08-21 |
JP2010251850A (ja) | 2010-11-04 |
EP2239593B1 (en) | 2012-03-14 |
US20100260239A1 (en) | 2010-10-14 |
TWI425242B (zh) | 2014-02-01 |
TW201126195A (en) | 2011-08-01 |
EP2239593A1 (en) | 2010-10-13 |
JP5353396B2 (ja) | 2013-11-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1457435B (zh) | Gps接收机装置及其数据存储装置 | |
CN101356449B (zh) | 为接收自多种卫星系统的卫星信号执行信号相关的方法和装置 | |
CN101084632B (zh) | 用于确定码相位的方法、设备、系统和模块 | |
CN102378921B (zh) | 用于gps软件接收器的方法和装置 | |
US8384592B1 (en) | FFT based acquisition techniques for satellite based navigation systems | |
CN101359044B (zh) | 计算并消除gnss强信号对弱信号干扰的方法 | |
JP2004527763A (ja) | 信号相関計算方法と装置 | |
JP2005500731A (ja) | スペクトラム拡散レシーバのアーキテクチャおよびその方法 | |
CN101329396B (zh) | 相干累积增进方法及系统、定位方法及系统、电子设备 | |
JP2003255040A (ja) | Gps受信機および受信方法 | |
US8362953B2 (en) | Sequential chip correlation array | |
CN101957450A (zh) | 接收器、信号处理方法和程序 | |
CN102928854A (zh) | 一种基于匹配滤波器的gps捕获单元设计方法 | |
CN101169478B (zh) | 定位装置、电子设备 | |
US20100172448A1 (en) | Method and receiver apparatus for determining a correlation value considering frequency correction | |
CN101859292B (zh) | 通信装置、信号处理方法、信号处理装置以及移动体 | |
JP4848146B2 (ja) | 測位信号を送信するための装置、その装置を備える測位システムおよび測位信号を送信するシステム | |
KR101015890B1 (ko) | 위성항법시스템 수신기의 신호획득 방법 및 장치 | |
CN102565816B (zh) | 一种载波平滑伪距分组平滑方法和装置 | |
CN102193096A (zh) | 信号捕获方法及信号捕获装置 | |
CN102081161A (zh) | 卫星信号捕获方法及卫星信号接收装置 | |
US7498981B2 (en) | Method and apparatus for real-time digital processing of satellite positional signals for fast acquisition and low SNR tracking | |
CN102162853B (zh) | 接收信号累计方法和接收装置 | |
JP2008151611A (ja) | 測位回路、電子機器、測位方法及びプログラム | |
CN103069300B (zh) | 接收器设备和接收方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20130102 Termination date: 20200402 |
|
CF01 | Termination of patent right due to non-payment of annual fee |