CN111769845B - 一种加权叠加信道化方法 - Google Patents

一种加权叠加信道化方法 Download PDF

Info

Publication number
CN111769845B
CN111769845B CN202010615694.5A CN202010615694A CN111769845B CN 111769845 B CN111769845 B CN 111769845B CN 202010615694 A CN202010615694 A CN 202010615694A CN 111769845 B CN111769845 B CN 111769845B
Authority
CN
China
Prior art keywords
filter
channelization
data
group
obtaining
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.)
Active
Application number
CN202010615694.5A
Other languages
English (en)
Other versions
CN111769845A (zh
Inventor
阎波
陈昳霏
郑行杰
魏磊
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
University of Electronic Science and Technology of China
Original Assignee
University of Electronic Science and Technology of China
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN202010615694.5A priority Critical patent/CN111769845B/zh
Publication of CN111769845A publication Critical patent/CN111769845A/zh
Application granted granted Critical
Publication of CN111769845B publication Critical patent/CN111769845B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B1/00Details of transmission systems, not covered by a single one of groups H04B3/00 - H04B13/00; Details of transmission systems not characterised by the medium used for transmission
    • H04B1/06Receivers
    • H04B1/16Circuits
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B17/00Monitoring; Testing
    • H04B17/30Monitoring; Testing of propagation channels
    • H04B17/382Monitoring; Testing of propagation channels for resource allocation, admission control or handover

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Digital Transmission Methods That Use Modulated Carrier Waves (AREA)

Abstract

本发明公开了一种加权叠加信道化方法,其包括以下步骤:S1、设置信道化参数;S2、获取滤波器的通带参数和阻带参数;S3、获取滤波器的阶数和系数;S4、将输入的中频信号进行复数混频,得到复基带信号;S5、对复基带信号进行分段锁存;S6、将分段锁存的数据和滤波器系数对应加权;S7、将所有加权结果编号后进行分组,并对每个组内数据进行求和,得到每组的和值;S8、对每组的和值进行相位校正,得到相位校正后的数据;S9、将相位校正后的数据进行快速傅里叶变换,得到各输出子信道的输出结果,完成加权叠加信道化。本发明可以有效的降低系统的运算量,并可以提高系统信号处理的实时性。

Description

一种加权叠加信道化方法
技术领域
本发明涉及通信领域,具体涉及一种加权叠加信道化方法。
背景技术
在软件无线电系统中,高速、宽带的数字信号需要经过信道分离和数字变频技术后变换为基带信号。但此时所有的基带信号都处于严重的过采样状态。过采样状态虽然有利于提高信号的信噪比,但其过高的输出数据率要求后续的信号处理要有很高的处理速度。因此数字信号处理前端还需要在保证信号不失真的情况下有效的降低信号的数据率。
在数字中频接收机结构中,在AD采样后需要接一个数字下变频模块将信号混频至基带并降低数据率,再进行后续的数字信号处理与分析。此模块具有两个功能,混频和抽取滤波。当需求的输出子信道数很多时的时候,普通的数字下变频方案对系统的数据处理速度要求高且资源占用多。
发明内容
针对现有技术中的上述不足,本发明提供的一种加权叠加信道化方法解决了现有下变频方案资源消耗大的问题。
为了达到上述发明目的,本发明采用的技术方案为:
提供一种加权叠加信道化方法,其包括以下步骤:
S1、设置信道化参数,包括输入采样率、输出子信道个数和抽取倍数;
S2、根据输入采样率和输出子信道个数获取滤波器的通带参数;根据输入采样率和抽取倍数获取滤波器的阻带参数;
S3、获取滤波器的阶数和系数;
S4、将输入的中频信号进行复数混频,得到复基带信号;
S5、根据抽取倍数和滤波器阶数对复基带信号进行分段锁存;
S6、将分段锁存的数据和滤波器系数对应加权,得到数量与滤波器阶数相等的加权结果;
S7、将所有加权结果编号后进行分组,并对每个组内数据进行求和,得到每组的和值;
S8、对每组的和值进行相位校正,得到相位校正后的数据;
S9、将相位校正后的数据进行快速傅里叶变换,得到各输出子信道的输出结果,完成加权叠加信道化。
进一步地,步骤S2中根据输入采样率和输出子信道个数获取滤波器的通带参数的具体方法为:
根据公式:
Figure BDA0002563669930000021
获取滤波器的通带参数Wp;其中Fs为输入采样率,K为输出子信道个数。
进一步地,步骤S2中根据输入采样率和抽取倍数获取滤波器的阻带参数的具体方法为:
根据公式:
Figure BDA0002563669930000022
获取滤波器的阻带参数Ws;其中Fs为输入采样率,D为抽取倍数。
进一步地,步骤S3的具体方法为:
获取滤波器的阶数L,并基于通带参数和阻带参数采用最小二乘法得到阶数为L时滤波器的系数h(n)。
进一步地,步骤S4的具体方法为:
根据公式:
X=Xi×e-j2πft
对输入的中频信号Xi进行复数混频,得到复基带信号X;其中j为虚数单位,π为圆周率,f为中频信号Xi的频率,t为时间,e为常数。
进一步地,步骤S5的具体方法为:
对复基带信号每隔D个点开始锁存一次,每次锁存的点位为L,完成对复基带信号进行分段锁存;其中D为抽取倍数,L为滤波器阶数。
进一步地,步骤S6的具体方法为:
对每次锁存的L个数据和滤波器系数h(n)一一对应相乘,得到L个加权结果;其中L为滤波器阶数。
进一步地,步骤S7中将所有加权结果编号后进行分组的具体方法为:
将所有加权结果从1开始依次编号,将所有加权结果分为K组,使第k组内包含第k+qK个加权结果;其中K为输出子信道个数,k=0,1,2,...,K,q为整数且
Figure BDA0002563669930000031
L为滤波器阶数。
进一步地,步骤S8的具体方法为:
根据公式:
Figure BDA0002563669930000032
对第k组和值
Figure BDA0002563669930000033
进行相位校正,得到第k组和值
Figure BDA0002563669930000034
相位校正后的数据
Figure BDA0002563669930000035
其中j为虚数单位,π为圆周率,m为数据点序号。
本发明的有益效果为:
1、本发明把运算环节安排在抽取之后,可以有效的降低系统的运算量,并可以提高系统信号处理的实时性。在需求多个输出子信道的情况下,由于部分资源的复用,可以更明显的降低资源消耗。
2、本发明的输出子信道个数、输入信号的采样率、抽取倍数和滤波器阶数均可以在一定范围内任意选择,更加灵活并且适用性强。
3、本发明在最初的混频时得到复基带信号,可以节省一半的子信道数量,并且无需额外进行复指数调制。
4、本发明在进行相位矫正时利用数据移位替代混频,可以减少运算量。
附图说明
图1为本发明的流程示意图;
图2为滤波器幅频响应图;
图3为实施例中输入信号的频谱图;
图4为实施例中所有输出子信道信号的频谱图。
具体实施方式
下面对本发明的具体实施方式进行描述,以便于本技术领域的技术人员理解本发明,但应该清楚,本发明不限于具体实施方式的范围,对本技术领域的普通技术人员来讲,只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内,这些变化是显而易见的,一切利用本发明构思的发明创造均在保护之列。
如图1所示,该加权叠加信道化方法包括以下步骤:
S1、设置信道化参数,包括输入采样率、输出子信道个数和抽取倍数;
S2、根据输入采样率和输出子信道个数获取滤波器的通带参数;根据输入采样率和抽取倍数获取滤波器的阻带参数;
S3、获取滤波器的阶数和系数;
S4、将输入的中频信号进行复数混频,得到复基带信号;
S5、根据抽取倍数和滤波器阶数对复基带信号进行分段锁存;
S6、将分段锁存的数据和滤波器系数对应加权,得到数量与滤波器阶数相等的加权结果;
S7、将所有加权结果编号后进行分组,并对每个组内数据进行求和,得到每组的和值;
S8、对每组的和值进行相位校正,得到相位校正后的数据;
S9、将相位校正后的数据进行快速傅里叶变换,得到各输出子信道的输出结果,完成加权叠加信道化。
步骤S2中根据输入采样率和输出子信道个数获取滤波器的通带参数的具体方法为:根据公式:
Figure BDA0002563669930000051
获取滤波器的通带参数Wp;其中Fs为输入采样率,K为输出子信道个数。
步骤S2中根据输入采样率和抽取倍数获取滤波器的阻带参数的具体方法为:根据公式:
Figure BDA0002563669930000052
获取滤波器的阻带参数Ws;其中Fs为输入采样率,D为抽取倍数。
步骤S3的具体方法为:获取滤波器的阶数L,并基于通带参数和阻带参数采用最小二乘法得到阶数为L时滤波器的系数h(n),在此过程中,生成的滤波器系数的幅频响应、相位响应如图2所示,图2中上下图的横坐标均为归一化的频率,上图的纵坐标为幅度值,单位dB;下图的纵坐标为相位值。
步骤S4的具体方法为:根据公式:
X=Xi×e-j2πft
对输入的中频信号Xi进行复数混频,得到复基带信号X;其中j为虚数单位,π为圆周率,f为中频信号Xi的频率,t为时间,e为常数。混频时采用复数混频已经完成实信号至基带复信号的转换,无需在信道化输出时额外进行复指数调制。并且与实信号相比,复信号没有冗余的负频谱部分,在信道化时不会产生无用的子信道。
步骤S5的具体方法为:对复基带信号每隔D个点开始锁存一次,每次锁存的点位为L,完成对复基带信号进行分段锁存;其中D为抽取倍数,L为滤波器阶数。
步骤S6的具体方法为:对每次锁存的L个数据和滤波器系数h(n)一一对应相乘,得到L个加权结果;其中L为滤波器阶数。
步骤S7中将所有加权结果编号后进行分组的具体方法为:将所有加权结果从1开始依次编号,将所有加权结果分为K组,使第k组内包含第k+qK个加权结果,即第k组包含第k、k+K、k+2K、…、k+L-K个加权结果;其中K为输出子信道个数,k=0,1,2,...,K,q为整数且
Figure BDA0002563669930000061
L为滤波器阶数。由于分组后求得的和值个数与子信道个数相同,因此后续进行的FFT计算的长度也与子信道个数相同。FFT后得到的结果即为各子信道的数据,不存在无效数据占用运算资源。
步骤S8的具体方法为:根据公式:
Figure BDA0002563669930000062
对第k组和值
Figure BDA0002563669930000063
进行相位校正,得到第k组和值
Figure BDA0002563669930000064
相位校正后的数据
Figure BDA0002563669930000065
其中j为虚数单位,π为圆周率,m为输入数据点序号,k=0,1,2,...,K,K为输出子信道个数。
在具体实施过程中,为了简化对每组的和值进行相位校正的运算,将公式
Figure BDA0002563669930000066
带入FFT运算式可得:
Figure BDA0002563669930000071
其中
Figure BDA0002563669930000072
为FFT输出数据,n为输出数据序号,N为FFT点数,WN为旋转因子,mod为取模运算。上述计算方式采用对数据进行距离为(mD)modK的循环移位的方式进行相位矫正,由于混频需要乘法和指数运算,而移位只是单纯的数据搬移无需运算,因此采用移位替代混频的方式对于每个数据点可以减少一次指数计算和一次乘法运算。在算法实现时,对于一个子信道数为8的信道化结构,采用这样的方式可以减少10%的DSP资源和50%的BRAM资源。
在本发明的一个实施例中,输入信号如图3所示,设定输入采样率为6.4M,输入信号为双频信号,信号频率为1.83M和1.77M。个子信道输出信号的频谱如图4所示,设定输出子信道的个数为8,抽取倍数为5。根据计算可得,信号频率为1.83M的信号应在第2子信道输出,输出频率的信号频率应为230KHz。信号频率为1.77M的信号同样应在第2子信道输出,输出频率的信号频率应为170KHz,与图4中显示的结果相符。
综上所述,本发明把运算环节安排在抽取之后,可以有效的降低系统的运算量,并可以提高系统信号处理的实时性。在需求多个输出子信道的情况下,由于部分资源的复用,可以更明显的降低资源消耗。

Claims (8)

1.一种加权叠加信道化方法,其特征在于,包括以下步骤:
S1、设置信道化参数,包括输入采样率、输出子信道个数和抽取倍数;
S2、根据输入采样率和输出子信道个数获取滤波器的通带参数;根据输入采样率和抽取倍数获取滤波器的阻带参数;
S3、获取滤波器的阶数和系数;
S4、将输入的中频信号进行复数混频,得到复基带信号;
S5、根据抽取倍数和滤波器阶数对复基带信号进行分段锁存;
S6、将分段锁存的数据和滤波器系数对应加权,得到数量与滤波器阶数相等的加权结果;
S7、将所有加权结果编号后进行分组,并对每个组内数据进行求和,得到每组的和值;
S8、对每组的和值进行相位校正,得到相位校正后的数据;
S9、将相位校正后的数据进行快速傅里叶变换,得到各输出子信道的输出结果,完成加权叠加信道化;
所述步骤S4的具体方法为:
根据公式:
X=Xi×e-j2πft
对输入的中频信号Xi进行复数混频,得到复基带信号X;其中j为虚数单位,π为圆周率,f为中频信号Xi的频率,t为时间,e为常数。
2.根据权利要求1所述的加权叠加信道化方法,其特征在于,所述步骤S2中根据输入采样率和输出子信道个数获取滤波器的通带参数的具体方法为:
根据公式:
Figure FDA0003420353810000021
获取滤波器的通带参数Wp;其中Fs为输入采样率,K为输出子信道个数。
3.根据权利要求1所述的加权叠加信道化方法,其特征在于,所述步骤S2中根据输入采样率和抽取倍数获取滤波器的阻带参数的具体方法为:
根据公式:
Figure FDA0003420353810000022
获取滤波器的阻带参数Ws;其中Fs为输入采样率,D为抽取倍数。
4.根据权利要求1所述的加权叠加信道化方法,其特征在于,所述步骤S3的具体方法为:
获取滤波器的阶数L,并基于通带参数和阻带参数采用最小二乘法得到阶数为L时滤波器的系数h(n)。
5.根据权利要求1所述的加权叠加信道化方法,其特征在于,所述步骤S5的具体方法为:
对复基带信号每隔D个点开始锁存一次,每次锁存的点位为L,完成对复基带信号进行分段锁存;其中D为抽取倍数,L为滤波器阶数。
6.根据权利要求1所述的加权叠加信道化方法,其特征在于,所述步骤S6的具体方法为:
对每次锁存的L个数据和滤波器系数h(n)一一对应相乘,得到L个加权结果;其中L为滤波器阶数。
7.根据权利要求1所述的加权叠加信道化方法,其特征在于,所述步骤S7中将所有加权结果编号后进行分组的具体方法为:
将所有加权结果从1开始依次编号,将所有加权结果分为K组,使第k组内包含第k+qK个加权结果;其中K为输出子信道个数,k=0,1,2,...,K,q为整数且
Figure FDA0003420353810000031
L为滤波器阶数。
8.根据权利要求1所述的加权叠加信道化方法,其特征在于,所述步骤S8的具体方法为:
根据公式:
Figure FDA0003420353810000032
对第k组和值
Figure FDA0003420353810000033
进行相位校正,得到第k组和值
Figure FDA0003420353810000034
相位校正后的数据
Figure FDA0003420353810000035
其中j为虚数单位,π为圆周率,m为数据点序号。
CN202010615694.5A 2020-06-30 2020-06-30 一种加权叠加信道化方法 Active CN111769845B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010615694.5A CN111769845B (zh) 2020-06-30 2020-06-30 一种加权叠加信道化方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010615694.5A CN111769845B (zh) 2020-06-30 2020-06-30 一种加权叠加信道化方法

Publications (2)

Publication Number Publication Date
CN111769845A CN111769845A (zh) 2020-10-13
CN111769845B true CN111769845B (zh) 2022-03-22

Family

ID=72724218

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010615694.5A Active CN111769845B (zh) 2020-06-30 2020-06-30 一种加权叠加信道化方法

Country Status (1)

Country Link
CN (1) CN111769845B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113546247B (zh) * 2021-07-09 2022-06-21 浙江大学 一种基于uhf rfid的输液过程监控装置及方法
CN117251717B (zh) * 2023-11-17 2024-02-09 成都立思方信息技术有限公司 多个不同信号同步信道化提取方法、装置、设备及介质

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1492613A (zh) * 2002-10-26 2004-04-28 乐金电子(中国)研究开发中心有限公 异步码分多址通信系统的符号映射及信道化装置
CN104735006A (zh) * 2015-03-12 2015-06-24 电子科技大学 一种基于加权叠加的直接中频信道化方法
CN107241107A (zh) * 2017-06-28 2017-10-10 电子科技大学 一种数字信道化滤波器组实现方法
CN108259024A (zh) * 2018-02-26 2018-07-06 中国人民解放军国防科技大学 一种过d采样多相dft滤波器组设计方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101667818B (zh) * 2009-09-03 2012-07-04 中国人民解放军信息工程大学 对实信号进行信道化处理的方法及装置
CN107634768B (zh) * 2017-09-11 2020-02-14 哈尔滨工程大学 一种基于fpga的mwc压缩采样宽带数字接收机pdw形成方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1492613A (zh) * 2002-10-26 2004-04-28 乐金电子(中国)研究开发中心有限公 异步码分多址通信系统的符号映射及信道化装置
CN104735006A (zh) * 2015-03-12 2015-06-24 电子科技大学 一种基于加权叠加的直接中频信道化方法
CN107241107A (zh) * 2017-06-28 2017-10-10 电子科技大学 一种数字信道化滤波器组实现方法
CN108259024A (zh) * 2018-02-26 2018-07-06 中国人民解放军国防科技大学 一种过d采样多相dft滤波器组设计方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
加权叠加结构滤波器组动态信道化结构研究;马宏 等;《遥测遥控》;20190331;第40卷(第2期);全文 *
无线接收机信道数字化的研究与FPGA实现;李超;《中国优秀硕士学位论文全文数据库 信息科技辑》;20180915(第09期);第3.3.2节,3.5.1-3.5.2节,4.5.1节,4.4节,表4-6 *

Also Published As

Publication number Publication date
CN111769845A (zh) 2020-10-13

Similar Documents

Publication Publication Date Title
CN111769845B (zh) 一种加权叠加信道化方法
CN103259754B (zh) 一种用于电力线载波通信的数字前端系统及其实现方法
US11424974B2 (en) Method and system for multi-carrier time division multiplexing modulation/demodulation
CN101378263A (zh) 基于数字中频的多载波数字接收机及多载波数字接收方法
CN101232488A (zh) 能降低正交频分复用系统峰均功率比的通信方法
CN110210101A (zh) 一种基于cem frm的动态非均匀窄过渡带滤波器组及设计方法
CN112532351A (zh) 一种加权分数傅里叶变换频域两分量信号的交织传输方法
CN112714090B (zh) 一种加权分数傅里叶变换扩展混合载波传输方法
CN107070836A (zh) 一种5g系统中基于frm技术的fbmc收发系统的设计方法
CN105553909B (zh) 一种迭代限幅滤波降低ofdm信号峰均比的方法
US8406320B2 (en) Digital clipping method for a transmitter of an orthogonal frequency division multiple access system
CN113872898A (zh) 一种基于有界成分分析的数字域自干扰抑制的方法和系统
CN107210985B (zh) 一种接收机及信号处理的方法
CN115776313B (zh) 弱化宽带短波数字基带信号载波相移的处理方法及装置
CN102223223B (zh) 多载波信号的接收方法及装置
CN109039968B (zh) 一种基于窄带物联网终端的数据接收方法及装置
CN201114162Y (zh) 一种基于数字中频技术的多载波数字接收机系统
CN101917376B (zh) 一种多载波数字接收机中数字下变频系统的两级变频方法
CN107204952A (zh) 一种滤波ofdm系统的子带滤波处理方法
Keerthana et al. FPGA implementation of FBMC baseband modular for 5G wireless communication
CN109067689A (zh) 一种f-ofdm通信系统中part的限幅方法
CN102281224B (zh) Ofdm接收机及其实现方法
CN201114161Y (zh) 多载波信号处理的数字下变频系统
CN100502267C (zh) Ofdm系统中信道估计的方法
CN109450599A (zh) 一种多项分数域Fourier再生变换的兼容扫描方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant