一种大型设备多点微位移测量方法与流程

文档序号:11914437阅读:来源:国知局

技术特征:

1.一种大型设备多点微位移测量方法,其特征在于,具体方法如下:

S1、在被测物表面上安装角发射器,在远端布置雷达,雷达发射多子带信号,其中子带内为线性调频信号,子带间为步进频率信号;

S2、对接收到的基频回波做脉冲压缩处理;

S3、对回波信号做匹配滤波、频域拼接及IDFT运算等处理进行相干合成;

S4、获得高分辨率距离像;

S5、对前后两次距离像做相位干涉测量,测得每个角反射器在雷达视线方向上的位移量。

2.根据权利要求1所述的一种大型设备多点微位移测量方法,其特征在于,所述的测量过程用公式表述为,

第一步将t1和t2时刻两次分别测得的信号S(R,t1)和S(R,t2)进行复相关运算,第二步通过公式(8)将第i个散射点载波干涉相位计算出来,第三步根据公式(11)将散射点的位移计算出来。

其中:

在t1和t2时刻,能够通过计算得出散射点i的干涉相位

<mrow> <msubsup> <mi>&Delta;&theta;</mi> <mn>12</mn> <mi>i</mi> </msubsup> <mo>=</mo> <mi>a</mi> <mi>n</mi> <mi>g</mi> <mrow> <mo>(</mo> <msup> <mi>S</mi> <mo>*</mo> </msup> <mo>(</mo> <msub> <mi>R</mi> <mi>i</mi> </msub> <mo>,</mo> <msub> <mi>t</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mo>&CenterDot;</mo> <mi>S</mi> <mrow> <mo>(</mo> <msubsup> <mi>R</mi> <mi>i</mi> <mo>&prime;</mo> </msubsup> <mo>,</mo> <msub> <mi>t</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> <mo>)</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>8</mn> <mo>)</mo> </mrow> </mrow>

第i个散射点在雷达距离像的位移量就可以表示为:

<mrow> <msubsup> <mi>&Delta;R</mi> <mn>12</mn> <mi>i</mi> </msubsup> <mo>=</mo> <mo>-</mo> <mfrac> <mi>&gamma;</mi> <mrow> <mn>4</mn> <mi>&pi;</mi> </mrow> </mfrac> <msubsup> <mi>&Delta;&theta;</mi> <mn>12</mn> <mi>i</mi> </msubsup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>11</mn> <mo>)</mo> </mrow> <mo>.</mo> </mrow>

3.根据权利要求1所述的一种大型设备多点微位移测量方法,其特征在于,所述的相位干涉测量,即通过信号往返于雷达和角反射器之间的相位差来测量角反射器的位移量。

4.根据权利要求1所述的一种大型设备多点微位移测量方法,其特征在于,公式具体过程如下:

s(t,k)=exp(-j4πf0(k)R/c)Sa(Tp0K(t-2R/c)) (1)

其中f0(k)为发射信号子脉冲的中心频率,R为目标点到雷达的距离,c为电磁波的传播速度,K为调频斜率,Tp0为子脉冲的长度,公式(1)中不考虑幅度的影响;

对公式(1)做傅里叶变换,同样的忽略幅度,得到各个子脉冲的脉冲压缩信号的频谱表达式为:

<mrow> <mi>S</mi> <mrow> <mo>(</mo> <mi>f</mi> <mo>,</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>G</mi> <msub> <mi>B</mi> <mn>0</mn> </msub> </msub> <mrow> <mo>(</mo> <mi>f</mi> <mo>)</mo> </mrow> <mi>exp</mi> <mrow> <mo>(</mo> <mo>-</mo> <mi>j</mi> <mn>4</mn> <msub> <mi>&pi;f</mi> <mn>0</mn> </msub> <mo>(</mo> <mi>k</mi> <mo>)</mo> <mi>R</mi> <mo>/</mo> <mi>c</mi> <mo>)</mo> </mrow> <mi>exp</mi> <mrow> <mo>(</mo> <mi>K</mi> <mn>4</mn> <mi>&pi;</mi> <mi>f</mi> <mi>R</mi> <mo>/</mo> <mi>c</mi> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow>

其中然后再进行频谱的平移,得到各个子脉冲经过脉冲压缩信号的频移量是Δf,

<mrow> <mi>&Delta;</mi> <mi>f</mi> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>f</mi> <mn>0</mn> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>f</mi> <mn>0</mn> </msub> <mo>=</mo> <mrow> <mo>(</mo> <mi>k</mi> <mo>-</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <mo>-</mo> <mfrac> <mi>n</mi> <mn>2</mn> </mfrac> <mo>)</mo> </mrow> <msub> <mi>B</mi> <mn>0</mn> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow>

经过频移以后的频谱表示为:

<mrow> <mtable> <mtr> <mtd> <mrow> <mi>S</mi> <mrow> <mo>(</mo> <mi>f</mi> <mo>,</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>G</mi> <msub> <mi>B</mi> <mn>0</mn> </msub> </msub> <mrow> <mo>(</mo> <mi>f</mi> <mo>-</mo> <mi>&Delta;</mi> <mi>f</mi> <mo>(</mo> <mi>k</mi> <mo>)</mo> <mo>)</mo> </mrow> <mi>exp</mi> <mrow> <mo>(</mo> <mo>-</mo> <mfrac> <mrow> <mi>j</mi> <mn>4</mn> <msub> <mi>&pi;f</mi> <mn>0</mn> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mi>R</mi> </mrow> <mi>c</mi> </mfrac> <mo>)</mo> </mrow> <mi>exp</mi> <mrow> <mo>(</mo> <mo>-</mo> <mfrac> <mrow> <mi>j</mi> <mn>4</mn> <mi>&pi;</mi> <mrow> <mo>(</mo> <mi>f</mi> <mo>-</mo> <mi>&Delta;</mi> <mi>f</mi> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>)</mo> <mi>R</mi> </mrow> <mi>c</mi> </mfrac> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>=</mo> <msub> <mi>G</mi> <msub> <mi>B</mi> <mn>0</mn> </msub> </msub> <mrow> <mo>(</mo> <mi>f</mi> <mo>-</mo> <mi>&Delta;</mi> <mi>f</mi> <mo>(</mo> <mi>k</mi> <mo>)</mo> <mo>)</mo> </mrow> <mi>exp</mi> <mrow> <mo>(</mo> <mo>-</mo> <mfrac> <mrow> <mi>j</mi> <mn>4</mn> <msub> <mi>&pi;f</mi> <mn>0</mn> </msub> <mi>R</mi> </mrow> <mi>c</mi> </mfrac> <mo>)</mo> </mrow> <mi>exp</mi> <mrow> <mo>(</mo> <mo>-</mo> <mfrac> <mrow> <mi>j</mi> <mn>4</mn> <mi>&pi;</mi> <mi>f</mi> <mi>R</mi> </mrow> <mi>c</mi> </mfrac> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow>

当频率步进量与子脉冲带宽相等的时候,即是Δf=B0,把频移之后的频谱进行叠加,得到叠加之后的数学表达式为:

<mrow> <mi>S</mi> <mrow> <mo>(</mo> <mi>f</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>G</mi> <mi>B</mi> </msub> <mrow> <mo>(</mo> <mi>f</mi> <mo>)</mo> </mrow> <mi>exp</mi> <mrow> <mo>(</mo> <mo>-</mo> <mfrac> <mrow> <mi>j</mi> <mn>4</mn> <msub> <mi>&pi;f</mi> <mn>0</mn> </msub> <mi>R</mi> </mrow> <mi>c</mi> </mfrac> <mo>)</mo> </mrow> <mi>exp</mi> <mrow> <mo>(</mo> <mo>-</mo> <mfrac> <mrow> <mi>j</mi> <mn>4</mn> <mi>&pi;</mi> <mi>f</mi> <mi>R</mi> </mrow> <mi>c</mi> </mfrac> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> </mrow>

最后,对叠加得到的信号频谱进行傅里叶逆变换,得到时域上的距离压缩信号;

各子带信号的带宽为B0,因此各个子带信号的距离分辨率可表示为c/2B0,而经过合成之后的信号带宽为B,合成以后的信号距离分辨率就可以表示为c/2B;

当在设备上安装了P个角反射器时,第i个角反射器的位置离雷达天线阵的真实距离表示为Ri,定义一个虚拟位置为R′i,R′i表示角反射器延迟线离雷达天线阵的虚拟距离,即R′i在真实距离Ri上叠加一固定距离;当系统运作时,雷达使用了天线阵发射无线电波信号,照射角反射器,延迟τt时间后,接受回波信号,进行匹配滤波运算处理;以下均用散点表示角反射器;

当只有一个散射点i时,其信号的数学表达式为:

<mrow> <msub> <mi>S</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>R</mi> <mo>,</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mo>|</mo> <msub> <mi>S</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>R</mi> <mo>,</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>|</mo> <msup> <mi>e</mi> <mrow> <msub> <mi>j&theta;</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow>

当存在P个散射点时,复合信号可表示为:

<mrow> <mi>S</mi> <mrow> <mo>(</mo> <mi>R</mi> <mo>,</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <msubsup> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>p</mi> </msubsup> <msub> <mi>S</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>R</mi> <mo>,</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <msubsup> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>p</mi> </msubsup> <mo>|</mo> <msub> <mi>S</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>R</mi> <mo>,</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>|</mo> <msup> <mi>e</mi> <mrow> <msub> <mi>j&theta;</mi> <mi>i</mi> </msub> </mrow> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow>

在t1和t2时刻,通过计算得出散射点i的干涉相位

<mrow> <msubsup> <mi>&Delta;&theta;</mi> <mn>12</mn> <mi>i</mi> </msubsup> <mo>=</mo> <mi>a</mi> <mi>n</mi> <mi>g</mi> <mrow> <mo>(</mo> <msup> <mi>S</mi> <mo>*</mo> </msup> <mo>(</mo> <msub> <mi>R</mi> <mi>i</mi> </msub> <mo>,</mo> <msub> <mi>t</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mo>&CenterDot;</mo> <mi>S</mi> <mrow> <mo>(</mo> <msubsup> <mi>R</mi> <mi>i</mi> <mo>&prime;</mo> </msubsup> <mo>,</mo> <msub> <mi>t</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> <mo>)</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>8</mn> <mo>)</mo> </mrow> </mrow>

<mrow> <msubsup> <mi>&Delta;&theta;</mi> <mn>12</mn> <mi>i</mi> </msubsup> <mo>=</mo> <mi>a</mi> <mi>n</mi> <mi>g</mi> <mrow> <mo>(</mo> <msubsup> <mi>S</mi> <mi>i</mi> <mo>*</mo> </msubsup> <mo>(</mo> <msub> <mi>R</mi> <mi>i</mi> </msub> <mo>,</mo> <msub> <mi>t</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mo>&CenterDot;</mo> <msub> <mi>S</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <msubsup> <mi>R</mi> <mi>i</mi> <mo>&prime;</mo> </msubsup> <mo>,</mo> <msub> <mi>t</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> <mo>)</mo> <mo>+</mo> <msubsup> <mi>&Sigma;</mi> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mi>k</mi> <mo>&NotEqual;</mo> <mi>i</mi> </mrow> <mi>p</mi> </msubsup> <msubsup> <mi>S</mi> <mi>k</mi> <mo>*</mo> </msubsup> <mrow> <mo>(</mo> <msub> <mi>R</mi> <mi>i</mi> </msub> <mo>,</mo> <msub> <mi>t</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mo>&CenterDot;</mo> <msub> <mi>S</mi> <mi>k</mi> </msub> <mrow> <mo>(</mo> <msubsup> <mi>R</mi> <mi>i</mi> <mo>&prime;</mo> </msubsup> <mo>,</mo> <msub> <mi>t</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> <mo>)</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>9</mn> <mo>)</mo> </mrow> </mrow>

其中,ang()表示求相位运算,为方便计算,当

但k≠1,在这种情况下就可以得到第i个散射点的干涉相位:

<mrow> <msubsup> <mi>&Delta;&theta;</mi> <mn>12</mn> <mi>i</mi> </msubsup> <mo>=</mo> <mi>ang</mi> <mrow> <mo>(</mo> <msubsup> <mi>S</mi> <mi>i</mi> <mo>*</mo> </msubsup> <mrow> <mo>(</mo> <msub> <mi>R</mi> <mi>i</mi> </msub> <mo>,</mo> <msub> <mi>t</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mo>&CenterDot;</mo> <msub> <mi>S</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <msubsup> <mi>R</mi> <mi>i</mi> <mo>&prime;</mo> </msubsup> <mo>,</mo> <msub> <mi>t</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> <mo>)</mo> </mrow> <mo>&cong;</mo> <mrow> <mo>(</mo> <msub> <mi>&theta;</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>t</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>&theta;</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>t</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mo>)</mo> </mrow> <mo>&cong;</mo> <mo>-</mo> <mfrac> <mrow> <mn>4</mn> <mi>&pi;</mi> </mrow> <mi>&gamma;</mi> </mfrac> <msubsup> <mi>&Delta;R</mi> <mn>12</mn> <mi>i</mi> </msubsup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>10</mn> <mo>)</mo> </mrow> </mrow>

其中,在公式(10)中,γ代表信号的波长;将得出第i个散射点在雷达距离像的位移量表示为:

<mrow> <msubsup> <mi>&Delta;R</mi> <mn>12</mn> <mi>i</mi> </msubsup> <mo>=</mo> <mo>-</mo> <mfrac> <mi>&gamma;</mi> <mrow> <mn>4</mn> <mi>&pi;</mi> </mrow> </mfrac> <msubsup> <mi>&Delta;&theta;</mi> <mn>12</mn> <mi>i</mi> </msubsup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>11</mn> <mo>)</mo> </mrow> <mo>.</mo> </mrow>

5.根据权利要求1所述的一种大型设备多点微位移测量方法,其特征在于,

所述的微位移的测量精度与载波波长、载波干涉相位测量精度相关,即当微位移测量精度提高,载波相位测量精度越高,信号的波长越短。

6.一种大型设备多点微位移测量装置,其特征在于,包括雷达和角发射器,其中,角发射器安装在被测大型设备表面;

还包括脉冲压缩模块,用于对接收到的基频回波做脉冲压缩处理;匹配滤波模块,用于对回波信号做匹配滤波处理;频域拼接模块,用于对于对回波信号做频域拼接处理;IDFT运算模块,用于对于对回波信号做IDFT运算处理;距离像形成模块,用于高分辨率一维距离像的生成;相位干涉测量模块,用于对前后两次距离像做相位干涉测量,测得每个角反射器在雷达视线方向上的位移量。

当前第2页1 2 3 
网友询问留言 已有0条留言
  • 还没有人留言评论。精彩留言会获得点赞!
1