一种用于安检成像的VCO扫频非线性校正方法与流程

本发明涉及一种用于安检成像的vco扫频非线性校正方法。

背景技术：

压控振荡器(voltagecontrolledoscillator，vco)是一种振荡频率随控制电压变化的可调信号源，广泛应用于雷达收发机、降频变换器、雷达、自动测试系统和数字无线通信系统中。作为安检成像的频率源的核心元件，vco的输出频率线性度决定了成像结果的距离分辨率和副瓣电平，然而受电调元件变容二极管固有调谐非线性和振荡器结构中高稳与宽带的矛盾的影响，vco的射频输出在扫频过程中产生扫频非线性和功率起伏，其自身的调频线性度一般只能达到2％～5％，严重影响实际测距精度和距离分辨力。为了提高安检仪系统的距离分辨能力，必须对vco进行扫频非线性校正。

目前vco线性校正技术主要分为3类：一是电抗补偿线性校正法，二是闭环线性校正法，三是开环线性校正法。电抗补偿线性校正的主要缺点是调试工作较繁琐，经线性校正后的vco电调特性的波纹起伏较大，线性校正精度不易得到可靠的保证；闭环线性校正法的缺点是实现方案复杂、成本较大；开环线性校正的优点是其实现方案简单，成本小，在对vco电调线性度有一定要求但又不是很高的场合特别适用。

技术实现要素：

本发明解决的技术问题是：克服现有技术的不足，提供了一种用于安检成像的vco扫频非线性校正方法，该方法在系统发射线性调频连续波脉冲的同时采集vco输出的m倍分频信号，采用最小二乘法获得n阶多项式拟合系数，并根据多项式拟合结果调整下一脉冲的vco的输入电压，从而迭代的完成线性校正。由于该方法在系统发射线性调频连续波的同时采集校准信号，不需要增加单独的校准时间，实时性能好。该方法校准的是vco输出信号的相位信息，可以长时间保持输出脉冲的初始相位的稳定，有效抑制了vco温度漂移问题，适合用于长时间合成孔径成像中。

本发明的技术解决方案是：一种用于安检成像的vco扫频非线性校正方法，包括如下步骤：

第一步，设置vco调整电压随时间t变化的初始电压波形；

第二步，将vco调整电压进行m倍分频，通过ad采集分频信号，并通过相位解缠绕得到分频信号的相位pi(t)；其中，i的初值为0；

第三步，求解第i+1次电压调整信号所需的差值δvi(t)；

第四步，生成新的vco电压控制信号vi+1(t)。

所述的设置vco调整电压随时间t变化的初始电压波形v0(t)为

式中，vc为待输出线性调频连续波信号中心频率对应的vco调整电压，vd为待输出线性调频连续波信号最高频率和最低频率对应的vco调整电压之差，t为时间变量，t为待输出线性调频连续波信号的脉冲时长。

所述的相位解缠绕得到分频信号的相位pi(t)为

式中，为单位纯虚数，p0为对相位信号pi(t)进行taylor展开后的第0项系数，vi(t)为迭代修正第i次后的vco调整电压随时间t变化的电压波形和，i为正整数且初值为1，f(vi(t))为迭代修正第i次后的vco调整电压对应的电调特性函数f(vi(t))。

所述的电调特性函数f(vi(t))的计算方法为

式中，表示对电调特性函数taylor展开后第n阶多项式系数an的估计值，vc为待输出线性调频连续波信号中心频率对应的vco调整电压，n为taylor展开后多项式的阶数，为对相位信号pi(t)进行taylor展开后的第0项系数p0的估计值，通过最小化相位误差

计算得到。

所述的第i+1次电压调整信号所需的差值δvi(t)为

所述的生成新的vco电压控制信号vi+1(t)的计算方法为

vi+1(t)＝vi(t)+αδv(t),0<α≤1

其中，p(t)为vco输出m分频信号相位对应的理想二次相位函数。

本发明与现有技术相比的优点在于：

(1)本发明通过对正常工作的vco输出信号的分频信号进行采集分析，计算得到下次线性调频连续波脉冲的电压调整信号。整个过程不需要额外的信号采集时间，具有连续校正的优点，适合应用于雷达占空比比较高的场合；

(2)本发明通过对vco输出信号的相位进行实时校正，相比只对频率进行校正的方法，具有温度漂移抑制好的优点，适合应用于长时间合成孔径成像中。

附图说明

图1为一种用于安检成像的vco扫频非线性校正方法硬件图；

图2为一种用于安检成像的vco扫频非线性校正方法流程图；

图3为基于多项式拟合的vco线性校正效果图，a为vco输入为锯齿波时vco输出频率误差，b为第1次迭代校正后的频率误差，多项式阶数为5，c为第5次迭代校正后的频率误差，d为校正前后雷达脉冲压缩效果。

具体实施方式

本发明针对现有技术的不足，提出一种基于多项式拟合的vco开环线性校正方法，该方法可以实时的完成线性校正，有效抑制了vco温度漂移问题。

如图1所示为本发明的硬件原理框图，图2为一种用于安检成像的vco扫频非线性校正方法流程图，其通过以下步骤实现：

第一步，通过fpga+da设置vco调整电压随时间t变化的初始电压波形v0(t)，0表示该电压为未经过校正的初始电压。设置vco调整电压为锯齿波信号，信号可以表示为：

式中，vc为待输出线性调频连续波信号中心频率对应的vco调整电压，vd为待输出线性调频连续波信号最高频率和最低频率对应的vco调整电压之差，t为时间变量，t为待输出线性调频连续波信号的脉冲时长。由于vco电调特性为非线性函数，在理想锯齿波信号的控制下，vco的振荡频率在一个脉冲时长t中不是按照线性规律变换的。所以vco的线性度校正的目标是通过调节vco的控制信号，使vco的振荡频率在一个脉冲时长t中按线性规律变换。

第二步，将vco输出的信号进行m倍分频，通过ad采集分频信号，并通过相位解缠绕得到分频信号的相位pi(t)。通常情况下vco输出信号的频率比较高，不适合通过ad直接采集。为了既能采集vco输出信号的相位信息，又能降低采样率，本发明中将vco输出信号进行m倍分频，经过带通滤波器滤除谐波分量后，通过ad采集得到vco的输出分频信号s(t)＝exp{jpi(t)}，并通过相位解缠绕方法得到分频信号的相位pi(t)。pi(t)体现了迭代修正第i次后vco的输出信号m倍分频后信号的相位随时间t的变化，与迭代修正第i次后的vco调整电压随时间t变化的电压波形vi(t)和vco的电调特性函数f(vi(t))相关：

式中，为单位纯虚数。p0为对相位信号pi(t)进行taylor展开后的第0项系数。

第三步，通过最小二乘法求解电调特性函数f(vi(t))的n阶多项式系数，得到电调特性函数f(vi(t))。vco电调特性函数f(vi(t))为缓变的单调增函数，其高分频较小，对成像雷达系统的脉冲压缩影响也小，可以忽略。因此对vco电调特性函数用有限n阶多项式进行近似，即电调特性函数可以近似为：

式中表示对电调特性函数taylor展开后第n阶多项式系数an的估计。式中，vc为待输出线性调频连续波信号中心频率对应的vco调整电压。

本发明中通过最小二乘法求解估计参数和通过最小化相位误差：

求多项式系数的估计参数从而得到电调特性函数f(vi(t))，式中表示对相位信号pi(t)进行taylor展开后的第0项系数p0的估计。

第四步，求解第i+1次电压调整信号vi+1(t)。设δvi(t)＝vi+1(t)-vi(t)为第i+1次电压调整信号vi+1(t)与第i次电压调整信号vi(t)的差值。假设第i+1次电压调整后vco输出的信号为理想的线性调频信号，vco输出m分频信号的相位为理想的二次相位函数p(t)，可以通过解方程求解得到电压调整信号δvi(t)：

第五步，生成新的vco电压控制信号vi+1(t)。在δvi(t)求解的过程中，存在测量误差、taylor截断误差、多项式系数估计误差，为了增强该方法的鲁棒性，增加调整步长系数α。vco电压调整信号vi+1(t)通过下列公式求解：

vi+1(t)＝vi(t)+αδv(t),0<α≤1(6)

式中α的取值，决定了迭代速度和稳定度，α越大迭代速度越快，α越小稳定度越好。

第六步，通过fpga+da设置vco调整电压。将生成好的vco调整电压存储至fpga中，等待下一次触发到来时，通过da输出至vco，并跳转至第二步进行下一轮迭代。

图3为该方法的安检仪系统中线性调频连续波生成模块中vco扫频非线性校正的效果。线性调频连续波生成模块中vco输出频率为3～3.75ghz，vco输出信号经过耦合输出，并进行20倍分频后，通过ad采集校正。如图所示，图a为vco输入为锯齿波时vco输出频率误差，图b为第1次迭代校正后的频率误差，多项式阶数为5，图c为第5次迭代校正后的频率误差，可以看到通过校正，线性调频连续波脉冲信号的线性度得到了很大的提高。

本发明说明书中未作详细描述的内容属本领域技术人员的公知技术。