CN114826456A - 一种针对多陪集采样中通道时钟偏移的校准方法 - Google Patents
一种针对多陪集采样中通道时钟偏移的校准方法 Download PDFInfo
- Publication number
- CN114826456A CN114826456A CN202210553439.1A CN202210553439A CN114826456A CN 114826456 A CN114826456 A CN 114826456A CN 202210553439 A CN202210553439 A CN 202210553439A CN 114826456 A CN114826456 A CN 114826456A
- Authority
- CN
- China
- Prior art keywords
- sampling
- matrix
- frequency
- signal
- channel
- 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
- 238000005070 sampling Methods 0.000 title claims abstract description 78
- 238000000034 method Methods 0.000 title claims abstract description 14
- 239000011159 matrix material Substances 0.000 claims description 62
- 238000001228 spectrum Methods 0.000 claims description 20
- 239000013598 vector Substances 0.000 claims description 19
- 230000014759 maintenance of location Effects 0.000 claims description 2
- 238000006243 chemical reaction Methods 0.000 abstract description 2
- 238000004891 communication Methods 0.000 description 4
- 230000008447 perception Effects 0.000 description 4
- 230000006870 function Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000011218 segmentation Effects 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 230000017105 transposition Effects 0.000 description 2
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 238000000844 transformation Methods 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B17/00—Monitoring; Testing
- H04B17/30—Monitoring; Testing of propagation channels
- H04B17/382—Monitoring; Testing of propagation channels for resource allocation, admission control or handover
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03M—CODING; DECODING; CODE CONVERSION IN GENERAL
- H03M1/00—Analogue/digital conversion; Digital/analogue conversion
- H03M1/10—Calibration or testing
- H03M1/1009—Calibration
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Theoretical Computer Science (AREA)
- Analogue/Digital Conversion (AREA)
Abstract
本发明针对次奈奎斯特采样方案中普遍存在的ADC时钟失配的问题,在现有的次奈奎斯特多陪集采样方案的基础上,提出一种针对多陪集采样中通道时钟偏移的校准方法。本发明的采样结构由I、Q两个采样支路构成,所有ADC以相同的速率对输入信号进行采样,同时两组ADC的输入信号在采样前分别进行上下变频,保证两个支路采样后得到的信号之间存在特定的频率偏移,基于提出的采样结构本发明建立了表达次奈奎斯特采样信号频谱与奈奎斯特频谱之间关系的方程,最后通过稀疏重构得到输入信号的频谱。本发明解决了现有方法中存在时间偏差的问题,能够准确估计偏移量,并消除偏移量带来的性能衰减。
Description
技术领域
本发明涉及次奈奎斯特宽带频谱感知技术,具体涉及一种多通道采样中时间偏差的校准及频谱感知技术。
技术背景
频谱感知的目的在于区分电磁频谱各个频带的占用情况,以便通过更好的频谱分配方式来充分利用空闲的频段,缓解频谱资源的短缺,提高通信系统的传输速率和质量。同时,由于频谱资源的稀缺、以及对更高速率更高质量通信系统的需求,已经有大量的通信系统使用GHZ以上的频段,甚至于毫米波频段。信号频率的提升给频谱感知的硬件实现带来了极大的挑战。在通信系统中,需要将模拟信号通过模数转换器(ADC,Analog to DigitalConverter)变成离散的数字信号,才能够被大多数电子设备识别处理。由香农采样定理可知,需要大于信号带宽两倍的采样率采样才能够将信号完全恢复,过高的信号带宽也就意味着过高的ADC采样率,这对ADC的性能提出了极高的要求,高速ADC往往功耗巨大、价格昂贵,而且分辨率也较低。同时使用高速ADC采样会得到大量的采样数据,对后续的数据处理提出了很高的要求。一系列次奈奎斯特采样方案能够避免上述问题,它们能够以较低采样率对稀疏宽带信号进行低成本地采样和高效地处理。
现有的次奈奎斯特采样方案都采用多通道的采样结构,如多陪集采样方案,其所有通道的采样率都远低于奈奎斯特采样率,每一个采样通道都包含特定的延时处理单元。采样得到的数据以不同的组合方式将信号的各个频带混叠到基带,通过输入信号和次奈奎斯特采样信号之间的关系建立感知模型,并使用稀疏重构方法估计输入信号的频谱。然而在硬件实现中,时延不准和时钟偏移的问题往往会造成严重的感知和重构误差。
现有的次奈奎斯特采样方案都假定各个通道的时延准确已知,多个通道间完全同步,没有考虑时延不准确和各通道间存在未知时钟偏移的情况,而以上两种原因引起的采样时间偏差在实际中普遍存在,在数十GHz的信号频率下,微小的时钟偏移也会给频谱感知带来性能衰减。
发明内容
本发明所要解决的技术问题是,针对现有的次奈奎斯特采样的采样时间普遍存在的偏差,提供一种不需要外部输入参考信号,就能够对各通道间的时间偏差进行盲校准的方法。
本发明为解决上述技术问题所采用的技术方案是,一种针对多陪集采样中通道时钟偏移的校准方法,包括以下步骤:
1、一种针对多陪集采样中通道时钟偏移的校准方法,其特征在于,包括以下步骤:
步骤1、信号输入先分别与频率为的余弦信号和正弦信号相乘后输入至由I、Q两组采样支路构成的采样结构,每组采样支路中共有p个通道,每条通道连接同一模数转换器ADC通过同一个时钟源控制以相同的频率对输入信号进行采样得到采样信号XA(f)与XB(f),采样频率为并且ADC控制I、Q两组采样支路间保持频率偏移T为奈奎斯特采样率下的时间间隔,L为降采样系数,t为时间变量,f为频率变量;
再将两组采样信号XA(f)与XB(f)分别分成L个段得到L个小段组成的向量xA(f)与xB(f);L为降采样系数;
最后利用对角矩阵D(f)和矩阵A得到两组采样信号XA(f)与XB(f)的矩阵形式yA(f)与yB(f):
其中,
diag表示对角矩阵,γi+τi为第i个通道的时间偏移量,γi为已知的延迟,τi为未知的采样时间偏差,φi为第i个通道的中间变量,通道序号i=1,2,…,p,
步骤2、计算利用yA(f)和yB(f)的关系先计算相关矩阵QAA、QAB:
设置相关矩阵QAA、QAB与对角矩阵Λ的关系:QAA=QABΛ;相关矩阵QAA、QAB以及对角矩阵Λ的各元素分别表示为(QAA)mn=amn、(QAB)mn=bmn、(Λ)mm=λm,()mn表示矩阵第m行第n列的元素,()mm表示矩阵第m行第m列的对角元素,m与n的取值范围为1,2,…,p;
对于第i个通道,其相关矩阵QAA、QAB以及对角矩阵Λ的各元素的关系可表示为:
步骤3、根据估计的时间偏差估计值得到偏差校准矩阵:
S为向量xA(f)中非零元素下标的集合,即支撑集
本发明的有益效果是,不需要外部输入参考信号,能够对各通道间的时间偏差进行盲校准,尽可能降低时间偏差对频谱感知性能的影响。
附图说明
图1为本发明的采样结构示意图;
图2为本发明方法与多陪集采样重构的归一化均方误差随平均时间偏移量的比较图;
图3为本发明方法的时间偏移量估计值的归一化均方误差随样本点数的变化图。
具体实施方式
下面对本发明具体实施方法作进一步推导和描述:
步骤1、使用如图1所示的采样结构,由I、Q两组采样支路构成,所有模数转换器ADC以相同的速率对输入信号进行采样,同时I、Q两组ADC的输入信号在采样前分别通过与正余弦函数相乘进行上下变频,保证两个支路采样后得到的信号之间存在特定的频率偏移使用同一个时钟源CLK就可以保证采样速率相同;
奈奎斯特采样率下的时间间隔为T,降采样系数为L,频率范围为的输入信号x(t)分为两组采样通道,分别与频率为的余弦信号和正弦信号相乘后输入采样频率为的ADC完成两种正交变换的XI(f)和XQ(f)得到两组观测数据xI[n]和xQ[n],每组数据中共有p个通道,即以及
正交变换得到的XI(f)和XQ(f)为:
其中,X(f)为傅里叶变换,f表示频率变量,T为奈奎斯特采样率下的时间间隔,L为降采样系数。
在不考虑噪声的情况下,XA(f)与XB(f)的关系有:
XA,i(f)为第i个通道采样后的信号频谱XA(f),n为采样点序号变量,γi+τi为第i个通道的时间偏移量,γi为已知的延迟,τi为未知的采样时间偏差。信号频谱XA,i(f)有
m为运算变量。
p个通道的信号频谱XA(f)通过对角矩阵D(f)和矩阵A写成矩阵形式yA(f)为:
同理可得,多个通道的信号频谱XB(f)的矩阵形式yB(f)为:
其中
yA(f)=[XA,1(f),XA,2(f),……,XA,p(f)]T
yB(f)=[XB,1(f),XB,2(f),……,XB,p(f)]T
xA(f)为将XA(f)分段后各得到L个小段组成的向量,xB(f)为将XB(f)分段后各得到L个小段组成的向量:
diag表示对角矩阵,T表示矩阵转置;
矩阵A为p×L的矩阵,A的第i行第k列元素为:
即
其中
步骤2、利用yA(f)和yB(f)的关系构造一组方程,并估计各通道时间偏差。将A代入yA(f)=D(f)AxA(f)和yB(f)=D(f)AxB(f),得
其中
有
定义相关矩阵QAA、QAB
H表示矩阵共轭转置;
有
即相关矩阵QAA、QAB以及对角矩阵Λ的关系为:
QAA=QABΛ
设相关矩阵QAA、QAB以及对角矩阵Λ的各元素(QAA)mn=amn、(QAB)mn=bmn、(Λ)mm=λm,三者的关系可表示为:
包含了时间偏差的误差函数为:
将时间偏差的估计转化为一个最优化问题:
则
令
可得
得到时间偏差估计值
步骤3、根据估计的时间偏差,得到校准后的感知模型
步骤4、为了重构得到XA(f),定义相关矩阵Q为:
H为矩阵共轭转置;
显然有,
由于Q、Z均为正定的赫米特矩阵,记Z和Q的秩分别为m、n。故可以分解得到
Z=PPH
Q=VVH
其中矩阵P的大小是L×m,矩阵V的大小是p×n,P和V的列向量分别是Z和Q的非零特征值对应的特征向量。S也就是P的非零行下标的集合。
另外,存在酉矩阵R,使
那么有
本发明的效果通过仿真图2和图3进行说明:
Claims (1)
1.一种针对多陪集采样中通道时钟偏移的校准方法,其特征在于,包括以下步骤:
步骤1、信号输入先分别与频率为的余弦信号和正弦信号相乘后输入至由I、Q两组采样支路构成的采样结构,每组采样支路中共有p个通道,每条通道连接同一模数转换器ADC通过同一个时钟源控制以相同的频率对输入信号进行采样得到采样信号XA(f)与XB(f),采样频率为并且ADC控制I、Q两组采样支路间保持频率偏移T为奈奎斯特采样率下的时间间隔,L为降采样系数,t为时间变量,f为频率变量;
再将两组采样信号XA(f)与XB(f)分别分成L个段得到L个小段组成的向量xA(f)与xB(f);
最后利用对角矩阵D(f)和矩阵A得到两组采样信号XA(f)与XB(f)的矩阵形式yA(f)与yB(f):
diag表示对角矩阵,γi+τi为第i个通道的时间偏移量,γi为已知的延迟,τi为未知的采样时间偏差,φi为第i个通道的中间变量,通道序号i=1,2,…,p,
步骤2、计算利用yA(f)和yB(f)的关系先计算相关矩阵QAA、QAB:
设置相关矩阵QAA、QAB与对角矩阵Λ的关系:QAA=QABΛ;估计对角矩阵Λ的各元素值从而得到各通道时间偏差的估计值:
步骤3、根据估计的时间偏差估计值得到偏差校准矩阵:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210553439.1A CN114826456B (zh) | 2022-05-20 | 2022-05-20 | 一种针对多陪集采样中通道时钟偏移的校准方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210553439.1A CN114826456B (zh) | 2022-05-20 | 2022-05-20 | 一种针对多陪集采样中通道时钟偏移的校准方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114826456A true CN114826456A (zh) | 2022-07-29 |
CN114826456B CN114826456B (zh) | 2022-12-13 |
Family
ID=82517468
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210553439.1A Active CN114826456B (zh) | 2022-05-20 | 2022-05-20 | 一种针对多陪集采样中通道时钟偏移的校准方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114826456B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117294337A (zh) * | 2023-10-30 | 2023-12-26 | 西安电子科技大学 | 一种csi载波频率偏移估计方法及装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060034381A1 (en) * | 2004-08-12 | 2006-02-16 | Nokia Corporation | Method and apparatus using coordinate interleaving to increase diversity in a MIMO system |
CN110365437A (zh) * | 2019-07-30 | 2019-10-22 | 电子科技大学 | 基于次奈奎斯特采样的快速功率谱估计方法 |
US20200259534A1 (en) * | 2019-02-07 | 2020-08-13 | California Institute Of Technology | Systems and Methods for Communicating by Modulating Data on Zeros in the Presence of Channel Impairments |
CN111901058A (zh) * | 2020-07-31 | 2020-11-06 | 电子科技大学 | 基于亚奈奎斯特采样的多天线辅助宽带频谱感知方法 |
-
2022
- 2022-05-20 CN CN202210553439.1A patent/CN114826456B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060034381A1 (en) * | 2004-08-12 | 2006-02-16 | Nokia Corporation | Method and apparatus using coordinate interleaving to increase diversity in a MIMO system |
US20200259534A1 (en) * | 2019-02-07 | 2020-08-13 | California Institute Of Technology | Systems and Methods for Communicating by Modulating Data on Zeros in the Presence of Channel Impairments |
CN110365437A (zh) * | 2019-07-30 | 2019-10-22 | 电子科技大学 | 基于次奈奎斯特采样的快速功率谱估计方法 |
CN111901058A (zh) * | 2020-07-31 | 2020-11-06 | 电子科技大学 | 基于亚奈奎斯特采样的多天线辅助宽带频谱感知方法 |
Non-Patent Citations (2)
Title |
---|
张瑜等: "基于多陪集采样的射电天文信号功率谱密度估计", 《云南大学学报(自然科学版)》 * |
林伟毅: "次奈奎斯特采样在超声波成像中的应用", 《电子设计工程》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117294337A (zh) * | 2023-10-30 | 2023-12-26 | 西安电子科技大学 | 一种csi载波频率偏移估计方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN114826456B (zh) | 2022-12-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
AU2003253088B2 (en) | Improvements relating to time-interleaved samplers | |
CN104779989B (zh) | 一种宽带阵列校正滤波器系数计算方法 | |
US20090033529A1 (en) | System and method of multi-channel signal calibration | |
US10873376B2 (en) | OAM multiplexing communication system, OAM multiplexing transmission device, OAM multiplexing receiving device, and OAM multiplexing communication method | |
US12057852B2 (en) | Sampling clock phase mismatch error estimation method and apparatus, and storage medium | |
CN114826456B (zh) | 一种针对多陪集采样中通道时钟偏移的校准方法 | |
CN110515033A (zh) | 一种基于Toeplitz矩阵恢复的欠通道测向系统及方法 | |
CN108881074A (zh) | 一种低精度混合架构下宽带毫米波信道估计方法 | |
CN112733476B (zh) | 基于ADCs的大规模MIMO接收阵列DOA估计方法 | |
CN108183710B (zh) | 基于fpga的时间交织模数转换器的校正方法及系统 | |
CN103780294B (zh) | 一种最大信噪比输出相控阵天线加权方法 | |
CN111314010A (zh) | 一种多通道延时可调的fx相关器及其实现方法 | |
CN111537947A (zh) | 单射频通道空间谱估计测向系统与方法 | |
CN110034759B (zh) | 前馈式全数字tiadc系统的采样时间误差校准模块及其方法 | |
CN114189293B (zh) | 一种宽带接收阵列天线通道幅相校准方法及系统 | |
CN107359877B (zh) | 超宽带信号的时间交织采样adc全数字盲补偿方法 | |
CN115913231B (zh) | 一种tiadc的采样时间误差数字估计方法 | |
CN103746699B (zh) | 基于旋转矩阵误差估计的交替采样系统信号重构方法 | |
CN114268346B (zh) | 非高斯噪声下电力线载波通信压缩感知信道估计方法 | |
CN108847909A (zh) | 一种基于压缩感知的宽带块稀疏频谱恢复方法 | |
CN111464184A (zh) | 一种基于压缩感知的时间交织adc | |
CN114553228A (zh) | 时间交织采样方法、装置、设备及存储介质 | |
CN113346902A (zh) | 基于tiadc复合输出的全数字校准结构及其校准方法 | |
CN115327484B (zh) | 基于单通道非均匀时间调制阵列的波达方向估计方法 | |
CN103780261A (zh) | 一种基于旋转矩阵的并行交替采样系统误差估计方法 |
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 |