一种基于正频率频带复信号的MWC系统传感矩阵系数校准方法与流程

一种基于正频率频带复信号的mwc系统传感矩阵系数校准方法
技术领域
1.本发明涉及一种基于正频率频带复信号的调制宽带转换器(modulated wideband converter，mwc)传感矩阵系数校准方法，属于高速模拟信息转换技术领域。


背景技术：

2.mwc是一种基于压缩感知理论的模拟信息转换器系统。mwc由多个通道并行的随机解调结构组成，每一通道包含混频器、模拟低通滤波器以及adc等器件，用周期性变化的随机在高低电平间跳变的高频信号，即伪随机序列进行混频，然后低通滤波，最后低速采样，将高频的模拟信号进行压缩并用相对低速的模数转换器(analog to digital converter,adc)进行采样。所以mwc可以认为是一种采样结构，该结构可以在输入信号频谱稀疏且未知的情况下，能够以远低于信号奈奎斯特频率的采样频率完成对输入信号的采样，并且利用采样信息和重构算法可以完美重构出原始输入信号波形。信号的重构实质上就是利用重构算法求解压缩感知问题y＝ax,其中y是mwc系统输出的采样序列构成的矩阵，a是理论传感矩阵代表着mwc系统输入和输出之间的关系，x就代表我们需要求解的原始输入信号。
3.然而，由于mwc系统硬件带有的例如器件非线性、延迟和噪声等各种非理想因素的影响，使输入输出之间的关系变为y＝cx，实际传感矩阵c中的系数相对于a中理论值产生了很大的偏差，如果直接用理论传感矩阵a进行信号重构，会导致信号重构失败甚至无法重构信号。因此我们需要利用校准方法校准出实际传感矩阵中的系数，只有利用校准之后的传感矩阵c来进行信号重构才能准确重构原始信号。
4.目前提出的校准方法采用的测试信号类型有：单个正弦或余弦信号的频点信号、多个正弦或余弦叠加的多频点信号、单个频带信号。以上的测试信号都是实数信号，它们的频谱正负频率轴上都有频谱内容并且关于零频率轴对称，利用这样的测试信号去校准mwc系统传感矩阵时会浪费低通滤波后mwc划分子带内一半的空间，即这一半的空间对于传感矩阵系数校准的贡献为零。本发明目的在于设计一种正频率频带测试信号，以充分利用子带内的空间进行mwc传感矩阵系数的校准。


技术实现要素：

5.本发明的目的是提供一种基于正频率频带测试信号的mwc系统传感矩阵系数校准方法，该种测试信号可以充分利用子带内的空间来进行传感矩阵系数的校准，一个测试信号校准出的传感矩阵系数的数量是原先双边频率单频带测试信号的两倍。
6.本发明是采用以下技术方案实现的：
7.由于mwc系统将信号的奈奎斯特频率范围划分为若干个宽度相等的子带，可以设计一种只在正频率子带内存在频谱内容的频带复信号作为测试信号，而且每个测试信号的频带分别位于每两个相邻子带的不同侧，每个频带的宽度设置为子带宽度的37.5％(该占比可变，但需小于50％)。然后将设计好的测试信号分实部和虚部分别加入mwc系统，将每个
通道输出的两组采样序列分别合成一组复数序列，将得到的复数采样序列作快速傅里叶变换(fast fourier transform，fft)得到其fft序列然后与滤波器补偿序列对应点相乘，截取每个频带位置对应的补偿后fft序列与测试信号的对应频带位置的fft序列进行对应点相除并取均值，并根据传感矩阵系数的共轭关系将求得的均值取共轭即可校准得到每个通道的两组传感矩阵系数。此外，还需要注意的是在校准之前应先测得实际电路中低通滤波器的频域响应并计算出滤波器的补偿序列，然后将实际低通滤波器在数字域补偿成理想低通滤波器。
8.具体步骤如下：
9.步骤一：采集或设置mwc系统的各项参数：奈奎斯特频率设置为f
nyq
，物理通道数为m，理想低通滤波器的通带范围设置为[
‑
f
s
/2，f
s
/2]，则adc的采样率和采样点数分别设置为f
s
和n，伪随机序列的周期设置为t
p
，其在每个周期内高低电平变化的次数设置为m，其中m设置为奇数且令(m 1)/4为整数，此外还要满足f
s
＝f
p
＝1/t
p
,其中f
p
为mwc系统划分子带的宽度。由以上参数可知，mwc输入信号的频谱在范围内被均分为m份子带，且每个子带的宽度为f
p
。先给每个子带添加序号，从负频率到正频率每个子带的序号令为
‑
m0,
…
,
‑
1,0,1,
…
,m0，其中m0＝(m
‑
1)/2，进而确定需要设计的测试信号的个数为(m 1)/4。
[0010]
步骤二：用频点扫描法测量低通滤波器的双边频率响应并按fft结果的排列形式(正频率频点值在左，负频率频点值在右)将其离散化为n个点的序列h[k]，k的取值为0到n
‑
1。根据b[k]＝i[k]/h[k]的关系求出滤波器补偿序列b[k],其中i[k]为理想低通滤波器的fft序列。
[0011]
步骤三：设计长度为mn的测试信号fft序列f
r
[n],其中下标r表示测试信号的序号，r的取值为1,2,
…
,r，其中r＝(m 1)/4为测试信号的总个数，n的取值范围为0到mn
‑
1。设计的频带宽度为b0(b0＝0.375f
p
)，则其在fft序列f
r
[n]中所等效的点数为一半子带宽度等效点数为间隔等效点数为以上计算符表示向下取整。则有fft序列f
r
[n]的构成规则为：首先r＝1时，令fft序列f1[n]的第n
g
个点到第n
g
n
b
‑
1个点的值都为1，再令fft序列f1[n]的第个点到第个点的值都为1，最后令fft序列f1[n]其他位置点的值都为0；然后当r≥2时，将第r
‑
1次的fft序列f
r
‑1[n]循环右移2n个点得到fft序列f
r
[n]。
[0012]
步骤四：将按以上规则构造的共r个测试信号的fft序列进行傅里叶逆变换得到测试信号的时域序列x
r
[n]，对求得的时域序列分别取实部和取虚部得到实部序列α
r
[n]和虚部序列β
r
[n]。
[0013]
步骤五：通过信号发生器将实部序列α
r
[n]生成的测试信号输入mwc系统，采集到每个通道的输出序列a
r,i
[k]，其中k的取值范围为0到n
‑
1，i表示每个通道的序号其取值为1到m；然后再通过信号发生器将虚部序列β
r
[n]生成的测试信号输入mwc系统，采集到每个通道的输出序列b
r,i
[k]其中k的取值范围为0到n
‑
1，i代表每个通道的序号取值为1到m。按照y
r,i
[k]＝a
r,i
[k] b
r,i
[k]
·
j的关系构造序列y
r,i
[k]，其中j为虚数单位。
[0014]
步骤六：对序列y
r,i
[k]进行快速傅里叶变换得到其fft序列y
r,i
[k]，再将y
r,i
[k]与补偿滤波器序列b[k]进行对应点相乘得到补偿后的fft序列y'
r,i
[k]。
[0015]
步骤七：根据mwc系统的频谱搬移关系，截取序列y'
r,i
[k]的第n
g
个点到第n
g
n
b
‑
1个点的值构成长度为n
b
的序列f
r,i,outl
，再截取序列y'
r,i
[k]的第个点到第个点的值构成长度为n
b
的序列f
r,i,outh
。
[0016]
步骤八：由于测试信号的频带幅值都为1，则可根据以下简化公式计算出mwc系统传感矩阵系数：
[0017][0018][0019]
其中下标l1＝2r
‑
2，运算符“sum(
·
)”表示求向量的所有元素之和，上标“*”表示复数的共轭。
[0020][0021][0022]
其中下标l2＝2r
‑
1，运算符“sum(
·
)”表示求向量的所有元素之和，上标“*”表示复数的共轭。求出所有系数后就完成了mwc系统传感矩阵c的校准。
[0023]
本发明的有益效果为：本发明利用正频率频带复信号作为mwc系统的校准测试信号，校准所得的传感矩阵系数为子带部分范围内的平均值，校准后的传感矩阵系数更具有可靠性，同时这种校准方法充分利用子带内的空间，减少了测试信号的数量，提高了校准的效率。
附图说明
[0024]
图1为mwc系统结构图；
[0025]
图2为mwc频谱划分和测试信号频谱示意图；
[0026]
图3(a)为原始信号时域波形；
[0027]
图3(b)为未采用本发明方法重构信号时域波形；
[0028]
图3(c)为采用本发明方法校准后重构信号时域波形；
[0029]
图4(a)为原始信号频谱；
[0030]
图4(b)为未采用本发明方法重构信号频谱；
[0031]
图4(c)为采用本发明方法校准后重构信号频谱。
具体实施方式
[0032]
下面将结合图与一个实例详细说明本发明的具体实施方式。
[0033]
如图1所示为调制宽带转换器(mwc)系统结构示意图。以通道数m＝25为例，每个通道由一个混频器、一个模拟低通滤波器和一个低速模数转换器(adc)构成。信号x(t)经过分
束器后同时进入每个通道，与伪随机序列p
i
(t)混频后再经过滤波采样输出采样序列y
i
[k]，设置信号的奈奎斯特频率f
nyq
＝1ghz，伪随机序列周期t
p
＝7.5
×
10
‑8s,每个周期内高低电平变化的次数m＝75，即划分子带的数量为75，每个子带的宽度f
p
＝13.33mhz，模拟低通滤波器通带范围设置为[
‑
6.67mhz，6.67mhz]，adc采样率也设置为13.33mhz，采样序列长度设置为n＝325，则可知输入信号长度为m*n＝24375。
[0034]
如图2所示为频谱划分和测试信号频谱示意图，测试信号频谱中存在两个频带，每个频带的宽度b0＝5mhz，对应到fft序列中的长度为n
b
＝121，频带所在的fft序列值设置为1。
[0035]
根据步骤一，首先采集并设置mwc系统的各项参数，在本次示例中通道数m＝25，奈奎斯特频率f
nyq
＝1ghz，伪随机序列周期t
p
＝7.5
×
10
‑8s,每个周期内高低电平变化的次数m＝75，即划分子带的数量为75，每个子带的宽度f
p
＝13.33mhz，给每个子带添加序号，从负频率到正频率每个子带的序号令为
‑
m0,
…
,
‑
1,0,1,
…
,m0，其中m0＝37，进而确定需要的测试信号的个数r＝19。低通滤波器通带范围设置为[
‑
6.67mhz，6.67mhz]，adc采样率也设置为13.33mhz，采样序列长度设置为n＝325，则可知输入信号长度为m*n＝24375。
[0036]
根据步骤二，用频点扫描法测量低通滤波器的双边频率响应并按fft结果的排列形式(正频率频点值在左，负频率频点值在右)将其离散化为325个点的序列h[k]，k的取值为0到n
‑
1。根据b[k]＝i[k]/h[k]的关系求出滤波器补偿序列b[k],其中i[k]为截止频率范围内幅值为1其余位置幅值为0的理想低通滤波器的fft序列。
[0037]
根据步骤三，设计长度为24375的测试信号fft序列f
r
[n],其中下标r表示测试信号的序号，r的取值为1,2,
…
,r，其中r＝(75 1)/4＝19为测试信号的总个数，n的取值范围为0到24374。设计的频带宽度为b0＝0.375f
p
＝5mhz,则其在fft序列f
r
[n]中所等效的点数为一半子带宽度等效点数为间隔等效点数为以上计算符表示向下取整。则有fft序列f
r
[n]的构成规则为：首先r＝1时，令fft序列f1[n]的第20个点到第140个点的值都为1，再令fft序列f1[n]的第182个点到第302个点的值都为1，最后令fft序列f1[n]其他位置点的值都为0；然后当r≥2时，将第r
‑
1次的fft序列f
r
‑1[n]循环右移650个点得到fft序列f
r
[n]。
[0038]
根据步骤四，将按以上规则构造的19个测试信号的fft序列进行傅里叶逆变换得到测试信号的时域序列x
r
[n]，对求得的时域序列分别取实部和取虚部得到实部序列α
r
[n]和虚部序列β
r
[n]。
[0039]
根据步骤五，通过信号发生器将实部序列α
r
[n]生成的测试信号输入mwc系统，采集到每个通道的输出序列a
r,i
[k]，其中k的取值范围为0到324，i表示每个通道的序号其取值为1到25；然后再通过信号发生器将虚部序列β
r
[n]生成的测试信号输入mwc系统，采集到每个通道的输出序列b
r,i
[k]其中k的取值范围为0到324，i代表每个通道的序号取值为1到25。按照y
r,i
[k]＝a
r,i
[k] b
r,i
[k]
·
j的关系构造序列y
r,i
[k]，其中j为虚数单位。
[0040]
根据步骤六，对序列y
r,i
[k]进行fft变换得到其fft序列y
r,i
[k]，再将y
r,i
[k]与补偿滤波器序列b[k]进行对应点相乘得到补偿后的fft序列y'
r,i
[k]。
[0041]
根据步骤七，根据mwc系统的频谱搬移关系，截取序列y'
r,i
[k]的第20个点到第140个点的值构成长度为121的序列f
r,i,outl
，再截取序列y'
r,i
[k]的第182个点到第302个点的值
构成长度为121的序列f
r,i,outh
。
[0042]
根据步骤八，由于测试信号的频带幅值都为1，则可根据以下简化公式计算出mwc系统传感矩阵系数：
[0043][0044][0045]
其中下标l1＝2r
‑
2，运算符“sum(
·
)”表示求向量的所有元素之和，上标“*”表示复数的共轭。
[0046][0047][0048]
其中下标l2＝2r
‑
1，运算符“sum(
·
)”表示求向量的所有元素之和，上标“*”表示复数的共轭。求出所有系数后就完成了mwc系统传感矩阵c的校准。
[0049]
校准完成后需要对校准矩阵进行信号重构测试，评价指标为信号重构信噪比，为重构的信号，x为输入的原信号，表示向量二范数的平方。原始信号设为载频中心在100mhz、200mhz、300mhz的频带信号，每个频带宽度为15mhz，输入信号的信噪比为60db。如图3(a)所示，为原始信号时域波形，图3(b)利用未校准矩阵重构信号时域波形，图3(c)为利用校准后矩阵重构信号时域波形。图4(a)为原始信号的频域响应，图4(b)为利用未校准矩阵重构信号的频域响应，图4(c)为利用校准后矩阵重构信号的频域响应。结果利用校准矩阵重构信号信噪比为56db。