CN101408871B - 数字互相关器的高效实现方法和设备 - Google Patents
数字互相关器的高效实现方法和设备 Download PDFInfo
- Publication number
- CN101408871B CN101408871B CN2007100598388A CN200710059838A CN101408871B CN 101408871 B CN101408871 B CN 101408871B CN 2007100598388 A CN2007100598388 A CN 2007100598388A CN 200710059838 A CN200710059838 A CN 200710059838A CN 101408871 B CN101408871 B CN 101408871B
- Authority
- CN
- China
- Prior art keywords
- data
- counter
- sequence
- random access
- access memory
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Synchronisation In Digital Transmission Systems (AREA)
Abstract
本发明以探求信号处理中普遍存在的接收序列与已知参考序列的相似性为背景,提出了一种数字互相关运算的高效硬件实现方案。其主要思想是使接收序列的第一个数与参考序列的所有数的共轭算乘法,其结果存入随机存取存储器(RAM),第二个数如法炮制其结果与RAM中相应位置的已存数据相加再存入其中,N(序列长度)次同样的操作后RAM中的最终数据即为互相关运算的结果,不同级别处理速度的操作由不同频率的参考时钟控制,该设备可用现场可编程门阵列(FPGA)构建。整个系统用极少的乘法器和很小的存储空间实现原本硬件实现十分庞杂的互相关运算,多时钟以及随加随存的设计极大的减少了硬件资源消耗,降低成本,同时提高了运算速度,适用于通信等信号处理领域中复数序列的互相关运算。
Description
技术领域
本发明涉及一种数字互相关器的高效实现方法和设备,可以测量未知数字信息与已知数字信息的相关程度,在模式识别、密码分析学以及无线通信等数字信号处理领域都有广泛应用。
背景技术
在信号处理领域中,互相关(有时也称为“互协方差”)是用来表示两个信号之间相似性的一个度量,通常通过与已知信号比较用于寻找未知信号中的特性。它是两个信号之间相对于时间的一个函数,有时也称为滑动点积。相关运算的公式如下:
实际中所测信号一般为有限长序列,所以计算公式可改写为:
其中,N为序列长度,共有2N-1个相关值。
互相关器常应用于通信系统的同步环节。正交频分复用(Orthogonal FrequencyDivision Multiplexing,OFDM)通过相互正交的子载波并行传输数据,是一种高效的数据传输方式,被誉为未来移动通信技术的核心。与传统通信手段相比,采用OFDM技术的无线数据通信系统具有抗频率选择性衰落、频带利用率高等特点,能够满足人们对同时传送数据、语音和图象的要求,所以受到人们越来越广泛的关注。近几年来,OFDM在无线通信领域得到了深入研究和广泛应用,如数字音频广播(DAB),数字视频广播(DVB)和高清晰度电视(HDTV)等。采用OFDM方案的协议有IEEE 802.11a和欧洲电信标准化协会的HIPERLAN/2等。
OFDM系统对定时和频率偏移敏感,特别是实际应用中可能与FDMA、TDMA和CDMA等多址方式结合使用时,时域和频率同步显得尤为重要。所以精确的同步措施对正确解调OFDM信号非常重要。在OFDM系统同步中,利用PN序列算整数频偏时,需要算接收序列与本地PN序列的互相关。其中,PN序列由巴克码组成,长度为N。故将巴克码扩展到,则要算序列的互相关。如果信号在时域信道传播过程中没有任何整数频偏,那么接收序列在频域中与本地PN序列副本算互相关后最大相关值应出现在中间m=0时。实际中如果最大相关值位置有了偏移,那么即可根据偏移量估算出整数频偏。
其他通信系统如CDMA、GSM等也离不开互相关的运算。在直接序列扩频通信系统中,发射机对发送数据进行扩频,不同用户拥有不同的相互正交的扩频码。在接收端的相关器中,接收序列与本地参考扩频序列算互相关,从而确定用户,并去掉扩频码,恢复原始数据。
OFDM系统和CDMA系统很多采用MQAM调制方式,其资源利用率高,QAM调制出的数据为复数。传输数据为复数的通信系统将成为一个趋势。复数的互相关器也是不可或缺的。如果系统调制方式采用MPSK等,即如果传输的是实数,那么要算实数序列的互相关。此数字互相关器同样适用,也可通过去掉复数乘法和虚部的相关器件,很容易的简化成实数互相关器。
模式识别主要是对语音波形、地震波、心电图、脑电图、图片、照片、文字、符号、生物的传感器等对象进行测量的具体模式进行分类和辨识,是信息科学和人工智能的重要组成部分。基于模板匹配的人脸识别方法大多是用归一化互相关,直接计算两幅图像之间的匹配程度。
发明内容
本发明目的是适应信号处理领域对互相关器的广泛需求,提出了一种数字互相关运算的高效实现方案。该系统由多参考时钟控制,用RAM、乘法器、加法器等一些简单而少量的资源实现原本硬件实现十分庞杂的相关算法,其各部分可用现场可编程门阵列(FPGA)构建,极大的减少了硬件资源消耗,降低成本,同时提高运算速度,具有广泛的适用性和实用价值。
本发明的技术方案:
接收序列的第一个数与参考序列的所有数的共轭进行复数的乘法运算,其结果分实部和虚部分别存入两个RAM的第一个地址至第N(序列长度)个地址;第二个数如法炮制其结果与RAM中第二个地址至第N+1个地址的已存数据相加再将结果存入其中;依此类推,第k个数与参考序列共轭相乘后,乘积结果与RAM中第k个地址至第k+N-1个地址的已存数据相加再存入其中。这样,接收完N个数后,RAM中的最终数据即为互相关运算的结果。系统中不同级别处理速度的操作由不同频率的参考时钟控制,接收数据的控制时钟频率相对最低,控制参考数据产生的时钟频率高于接收数据的时钟,存取乘积结果的RAM的控制时钟最高。
本发明的有益效果:
本发明设计了一种高效的数字互相关器,用极少的乘法器和很小的存储空间等实现原本硬件实现十分庞杂的互相关运算,多时钟以及随加随存的设计极大的减少了硬件资源消耗,降低成本,同时提高了运算速度,适用于通信等信号处理领域。
附图说明
图1是数字互相关器的工作原理示意图
图2是数字互相关器的结构框图
图3是复数乘法部分的结构框图
图4是实数互相关器的结构框图
具体实施方式
以下结合附图和通过实施例对本发明的具体实施方式作进一步说明:
本发明介绍了一种数字互相关器的高效实现方法,其特征在于:本方法包括以下步骤:
a.每接收一个数据的同时,地址发生器控制只读存储器,从中依次读出参考序列的N个数据,并依次与该接收数据进行复数的乘法运算,则可得到一个长度为N的乘积序列,接收N个数据,就得到N个这样的乘积序列,第k个接收数据与参考序列相乘得到第k个乘积序列,复数的乘法运算由四个无符号数乘法器实现,N为所要进行互相关运算的接收序列与参考序列的长度,k为接收数据的序号,k大于等于1且小于等于N;
b.步骤a中所述第k个乘积序列的N个数据逐个进入有符号数加法器级联随机存取存储器的结构,与随机存取存储器中第k至第k+N-1个地址的共N个数据对应相加,相加结果再存入第k至第k+N-1个地址中,从而实现乘积序列的随加存储功能,当接收完所有N个数据,随机存取存储器中第1至第2N-1个地址中存放的即为所求的2N-1个相关值的实部和虚部;
c,采用多个参考时钟来控制,接收数据的控制时钟CLK1的频率值等于接收数据的速率值,控制参考序列产生的时钟CLK2的频率大于等于倍的CLK1的频率,随机存取存储器的控制时钟CLK3的频率大于等于倍的CLK2的频率,其中1/m为CLK1的占空比,1/n为CLK2的占空比。
步骤a中所述的接收数据的实部与虚部并行传输。
步骤a中所述的只读存储器的数量为2,分别为ROM1和ROM2,其容量都为M×N,M为数据宽度,N为步骤a所述的接收序列与参考序列的长度,将参考序列的实部依次储存在ROM1中,虚部的共轭依次存储在ROM2中,ROM1和ROM2的参考时钟为步骤c所述的CLK2。
步骤a中所述的地址发生器由一计数器构成,该计数器长度为步骤a中所述的N,该计数器由步骤c中所述CLK2的上升沿控制。
步骤a中所述的复数的乘法运算又分为以下步骤:
d.构建四个无符号数乘法器,分别为乘法器1,乘法器2,乘法器3,乘法器4;
e.提取出接收数据实部虚部的符号以及参考数据实部虚部的符号,分别为符号1,符号2,符号3,符号4;
f.将去除符号后的接收数据实部和参考数据实部送入乘法器1相乘,接收数据虚部和参考数据虚部送入乘法器2相乘,接收数据实部和参考数据虚部送入乘法器3相乘,接收数据虚部和参考数据实部送入乘法器4相乘;
g.将符号分别加入四个相乘结果,其中相乘结果1的符号由符号1、3决定,相乘结果2的符号由符号2、4决定,相乘结果3的符号由符号1、4决定,相乘结果4的符号由符号2、3决定;
h.构建一个有符号数加法器和一个有符号数减法器;
i.将相乘结果1和2送入减法器相减得到复数乘法结果的实部,相乘结果3和4送入加法器相加得到复数乘法结果的虚部。
步骤b中所述的随机存取存储器的数量为2,分别为RAM1和RAM2,其容量均为M×(2N-1),M为所述数据宽度,N为步骤a中所述的接收序列与参考序列的长度,RAM1和RAM2的参考时钟为步骤c中所述的CLK3。
步骤b中所述的有符号数加法器级联随机存取存储器的结构的数量为2,包括两个有符号数加法器,分别为加法器1和加法器2,以及所述RAM1和RAM2,RAM1的输出端接回加法器1的一个输入端,加法器1另一个输入端接乘积结果的实部,RAM2的输出端接回加法器2的一个输入端,加法器2另一个输入端接乘积结果的虚部,从而实现乘积序列实部和虚部的累加并存储的功能。
步骤b中需要构建一个随机存取存储器的公共控制器,用于控制地址发生和写使能,该公共控制器的地址发生功能由两个计数器实现,分别为计数器1和计数器2,计数器1和计数器2的长度为步骤a中所述的N,计数器1由步骤c中所述CLK1的上升沿控制,其数值表示的是乘积序列与随机存取存储器中的数据相加再存入的起始地址,计数器2由步骤c中所述CLK2的上升沿控制,其数值表示的是乘积序列中数据的序号,随机存取存储器的地址即为计数器1的数值与计数器2的数值的和,随机存取存储器的写使能由一个长度为3的计数器3控制,其由步骤c中所述CLK3的上升沿控制,当计数器3记第一个数时,写使能禁止,此时读出随机存取存储器相应地址的数据送到所述有符号数加法器中与乘积数据相加;当计数器3记第二个数时,写使能有效,此时将相加结果再写入随机存取存储器的同一地址;当计数器3记第三个数时,写使能禁止。
步骤c所述的多时钟控制,其特征在于:高速时钟所控制的计数器的复位由低速时钟的高电平触发,其计数开始由低速时钟的低电平使能。
实施例
本发明应用于32QAM正交频分复用系统的同步中,系统使用16个子载波,每帧由1个空符号、1个PN序列和62个OFDM符号组成,数据接收频率为24KHz。同步过程包括帧起始位估计、符号定时、小数频偏估计与补偿、整数频偏估计与补偿等。其中利用PN序列来估计整数频偏,在发送端的频域中加入PN序列,经过IFFT后变换到时域进行传输,经过信道产生频偏,接收端将其再做FFT变换到频域,将变换后得到的接收PN序列与已知PN序列本地副本算互相关,看频域中PN序列循环移位的位数,根据数字信号处理中的循环移位定理即可判定时域中的整数频偏。这里,PN序列本地副本选用巴克码的复数域变形形式,如果没有整数频偏,算得的最大互相关值位置应该在中间,如果有了整数频偏,即可根据最大互相关值出现的位置来判断频域中移位的位数。
PN本地副本为[+1,+i,-1,-i,+1,-i,+1,-i,+1,-i,-1,+i,+1,0,0,0],序列长度为16。A/D选用5倍过采样,采样频率为120KHz,则PN本地副本序列长度为80,即80个数与80个数算互相关。接收数据量化后宽度为5bits,本地PN副本数据宽度为2bits。接收CLK1为采样频率120KHz,占空比用1/4,其实现方法是先将占空比为1/2的时钟二分频,把分频后的和分频前的时钟进行与操作,即得占空比为1/4的CLK1。CLK2为15MHz,占空比同为1/4。CLK3为60MHz。
构建两个2bits×128words的ROM,用前80个位置存放80个PN本地参考序列的实部和虚部。其地址发生器中计数器长度为80。
将接收的数据逐个与80个本地PN数据相乘。注意,此处复数乘法中的无符号数乘法器可以简化,因为其中一个乘数,即本地序列的幅值非0即1,所以可用一个二选一数字选择器代替。
Matlab仿真可得互相关值最大数据宽度为9bits,所以构建两个容量为9bits×256words的RAM1和RAM2,用前159个位置存放159个相关值。其控制器中两个计数器的长度分别为80和3。
按照具体实施方式所述构建该数字相关器,当第k个(1≤k≤80)接收数据与已知参考序列相乘得到第k个乘积序列时,其数据实部和虚部逐个进入加法器级联RAM的结构,与RAM中第k个至第k+79个地址的数据对应相加,相加结果再存入第k个至第k+79个地址中。当接收完所有80个数据(k由1依次取到80),RAM1和RAM2中第1个地址至第159个地址中存放的数据即为所求的159个相关值的实部和虚部。
Claims (8)
1.一种数字互相关器的高效实现方法,其特征在于:本方法包括以下步骤:
a.每接收一个数据的同时,地址发生器控制只读存储器,从中依次读出参考序列的N个数据,并依次与该接收数据进行复数的乘法运算,则可得到一个长度为N的乘积序列,接收N个数据,就得到N个这样的乘积序列,第k个接收数据与参考序列相乘得到第k个乘积序列,复数的乘法运算由四个无符号数乘法器实现,N为所要进行互相关运算的接收序列与参考序列的长度,k为接收数据的序号,k大于等于1且小于等于N;
b.步骤a中所述第k个乘积序列的N个数据逐个进入有符号数加法器级联随机存取存储器的结构,与随机存取存储器中第k至第k+N-1个地址的共N个数据对应相加,相加结果再存入第k至第k+N-1个地址中,从而实现乘积序列的随加存储功能,当接收完所有N个数据,随机存取存储器中第1至第2N-1个地址中存放的即为所求的2N-1个相关值的实部和虚部;
2.根据权利要求1所述的一种数字互相关器的高效实现方法,其特征在于:步骤a中所述的接收数据的实部与虚部并行传输。
3.根据权利要求1所述的一种数字互相关器的高效实现方法,其特征在于:步骤a中所述的只读存储器的数量为2,分别为ROM1和ROM2,其容量都为M×N,M为数据宽度,N为步骤a所述的接收序列与参考序列的长度,将参考序列的实部依次储存在ROM1中,虚部的共轭依次存储在ROM2中,ROM1和ROM2的参考时钟为步骤c所述的CLK2。
4.根据权利要求1所述的一种数字互相关器的高效实现方法,其特征在于:步骤a中所述的地址发生器由一计数器构成,该计数器长度为步骤a中所述的N,该计数器由步骤c中所述CLK2的上升沿控制。
5.根据权利要求3所述的一种数字互相关器的高效实现方法,其特征在于:步骤b中所述的随机存取存储器的数量为2,分别为RAM1和RAM2,其容量均为M×(2N-1),M为所述数据宽度,N为步骤a中所述的接收序列与参考序列的长度,RAM1和RAM2的参考时钟为步骤c中所述的CLK3。
6.根据权利要求5所述的一种数字互相关器的高效实现方法,其特征在于:步骤b中所述的有符号数加法器级联随机存取存储器的结构的数量为2,包括两个有符号数加法器,分别为加法器1和加法器2,以及所述RAM1和RAM2,RAM1的输出端接回加法器1的一个输入端,加法器1另一个输入端接乘积结果的实部,RAM2的输出端接回加法器2的一个输入端,加法器2另一个输入端接乘积结果的虚部,从而实现乘积序列实部和虚部的累加并存储的功能。
7.根据权利要求6所述的一种数字互相关器的高效实现方法,其特征在于:步骤b中需要构建一个随机存取存储器的公共控制器,用于控制地址发生和写使能,该公共控制器的地址发生功能由两个计数器实现,分别为计数器1和计数器2,计数器1和计数器2的长度为步骤a中所述的N,计数器1由步骤c中所述CLK1的上升沿控制,其数值表示的是乘积序列与随机存取存储器中的数据相加再存入的起始地址,计数器2由步骤c中所述CLK2的上升沿控制,其数值表示的是乘积序列中数据的序号,随机存取存储器的地址即为计数器1的数值与计数器2的数值的和,随机存取存储器的写使能由一个长度为3的计数器3控制,其由步骤c中所述CLK3的上升沿控制,当计数器3记第一个数时,写使能禁止,此时读出随机存取存储器相应地址的数据送到所述有符号数加法器中与乘积数据相加;当计数器3记第二个数时,写使能有效,此时将相加结果再写入随机存取存储器的同一地址;当计数器3记第三个数时,写使能禁止。
8.根据权利要求1所述的一种数字互相关器的高效实现方法,其特征在于:高速时钟所控制的计数器的复位由低速时钟的高电平触发,其计数开始由低速时钟的低电平使能。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2007100598388A CN101408871B (zh) | 2007-10-12 | 2007-10-12 | 数字互相关器的高效实现方法和设备 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2007100598388A CN101408871B (zh) | 2007-10-12 | 2007-10-12 | 数字互相关器的高效实现方法和设备 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101408871A CN101408871A (zh) | 2009-04-15 |
CN101408871B true CN101408871B (zh) | 2010-07-14 |
Family
ID=40571887
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2007100598388A Expired - Fee Related CN101408871B (zh) | 2007-10-12 | 2007-10-12 | 数字互相关器的高效实现方法和设备 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101408871B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105808206B (zh) * | 2016-03-04 | 2019-01-08 | 广州海格通信集团股份有限公司 | 基于ram实现乘法运算的方法及其系统 |
CN110018456B (zh) * | 2019-04-11 | 2020-11-06 | 上海卫星工程研究所 | 基于Matlab的星载合成孔径雷达数据快速帧头检测方法 |
CN112558921B (zh) * | 2020-12-23 | 2021-08-31 | 广东高云半导体科技股份有限公司 | 一种实现互相关运算的装置及方法 |
CN112597432A (zh) * | 2020-12-28 | 2021-04-02 | 华力智芯(成都)集成电路有限公司 | 基于fft算法的复序列互相关在fpga上的加速实现方法及系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1564505A (zh) * | 2004-04-02 | 2005-01-12 | 清华大学 | 正交频分复用系统中帧同步的动态捕获装置 |
CN1855904A (zh) * | 2005-04-29 | 2006-11-01 | 中国科学院上海微系统与信息技术研究所 | 一种基于ofdm系统的符号定时方法 |
-
2007
- 2007-10-12 CN CN2007100598388A patent/CN101408871B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1564505A (zh) * | 2004-04-02 | 2005-01-12 | 清华大学 | 正交频分复用系统中帧同步的动态捕获装置 |
CN1855904A (zh) * | 2005-04-29 | 2006-11-01 | 中国科学院上海微系统与信息技术研究所 | 一种基于ofdm系统的符号定时方法 |
Also Published As
Publication number | Publication date |
---|---|
CN101408871A (zh) | 2009-04-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101536336B (zh) | 用于快速小区搜索的方法和装置 | |
CN101083645B (zh) | 一种低复杂度ofdm快速同步的方法 | |
CN100477569C (zh) | 使用长码和短码的多载波cdma方案的方法、发送设备和接收设备 | |
CN102185685A (zh) | 移动终端发射随机接入信号方法及系统 | |
CN101068232B (zh) | 得到信道时域响应方法及装置、ofdm符号精同步方法及装置 | |
CN101408871B (zh) | 数字互相关器的高效实现方法和设备 | |
CN101001235A (zh) | 一种弱能量并行pn序列时间同步与频率同步方法 | |
Abdzadeh-Ziabari et al. | An improved timing estimation method for OFDM systems | |
CN101188592A (zh) | 用于多载波系统的同步方法及系统 | |
CN1327642C (zh) | 消除正交频分复用信号时频偏差影响的帧同步电路和方法 | |
CN100591059C (zh) | 接收多载波传输的方法、系统和接收机 | |
CN106160969A (zh) | 一种lte下行同步数据发射配置与接收方法 | |
US7130361B1 (en) | Telecommunication device with analog fourier transformation unit | |
CN101267244B (zh) | 定时跟踪方法 | |
CN101022438A (zh) | 兼容dab的数字广播接收机载波同步方法及系统 | |
CN104486286A (zh) | 一种连续子载波ofdma系统的上行帧同步方法 | |
CN102215200B (zh) | 一种整数频偏估计的方法和装置 | |
CN101729466B (zh) | 对有尖峰干扰的空白频谱进行检测和使用的装置及方法 | |
CN114374589A (zh) | 基于fpga的5g-nr系统pss互相关算法 | |
CN100579101C (zh) | 一种ofdm符号同步方法及装置 | |
CN103095627A (zh) | 一种正交频分复用技术系统同步方法和电子设备 | |
US6950475B1 (en) | OFDM receiver clock synchronization system | |
Kadhum et al. | Digital chunk processing with orthogonal GFDM doubles wireless channel capacity | |
JP5680734B1 (ja) | 帯域分散伝送システム、帯域分散伝送方法、送信装置および受信装置 | |
TWI847620B (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 | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20100714 Termination date: 20131012 |