发明名称 相控阵雷达电控波束稳定方法
摘要 本发明提供一种相控阵雷达电控波束稳定方法,解决载体扰动带来的天线波束指向变化对目标探测与跟踪的影响这一问题。第一步:初始化,相控阵雷达上电,波束稳定模块各参量初始化;第二步:利用上一时刻四元数与当前时间间隔内角速度测量值,通过数值积分法解算天线相对姿态;第三步:根据上一步得到的天线相对姿态以及上一时刻目标跟踪向量,通过坐标变换反推得到保持天线波束稳定所需的波束补偿角;第四步:获取目标测量值;第五步:利用计算得到的波束补偿角、目标方位误差角和俯仰误差角共同实现波束控制,并更新目标跟踪向量、四元数用于算法迭代,从而抵消当前时间间隔内载体扰动对波束的影响,保持波束始终正对目标。
申请公布号 CN102981151B 申请公布日期 2015.03.04
申请号 CN201210488000.1 申请日期 2012.11.26
申请人 北京理工大学 发明人 李阳;肖增利;陆顶洪
分类号 G01S7/40(2006.01)I;G01S13/66(2006.01)I 主分类号 G01S7/40(2006.01)I
代理机构 北京理工大学专利中心 11120 代理人 高燕燕
主权项 一种相控阵雷达电控波束稳定方法,其特征在于,假定坐标系o‑xyz由原坐标系o‑xyz按如下方式得到:先绕y轴旋转角度ψ,接着绕经ψ变换后的z轴旋转角度θ,最后绕经ψ、θ变换后的x轴旋转角度γ,则可得如下关系:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mfenced open='[' close=']'><mtable><mtr><mtd><msup><mi>x</mi><mo>&prime;</mo></msup></mtd></mtr><mtr><mtd><msup><mi>y</mi><mo>&prime;</mo></msup></mtd></mtr><mtr><mtd><msup><mi>z</mi><mo>&prime;</mo></msup></mtd></mtr></mtable></mfenced><mo>=</mo><msub><mi>A</mi><mrow><mi>z</mi><mo>&RightArrow;</mo><msup><mi>z</mi><mo>&prime;</mo></msup></mrow></msub><mfenced open='[' close=']'><mtable><mtr><mtd><mi>x</mi></mtd></mtr><mtr><mtd><mi>y</mi></mtd></mtr><mtr><mtd><mi>z</mi></mtd></mtr></mtable></mfenced><mo>=</mo><mi>K</mi><mrow><mo>(</mo><mi>&gamma;</mi><mo>)</mo></mrow><mi>L</mi><mrow><mo>(</mo><mi>&theta;</mi><mo>)</mo></mrow><mi>M</mi><mrow><mo>(</mo><mi>&psi;</mi><mo>)</mo></mrow><mfenced open='[' close=']'><mtable><mtr><mtd><mi>x</mi></mtd></mtr><mtr><mtd><mi>y</mi></mtd></mtr><mtr><mtd><mi>z</mi></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000568623910000011.GIF" wi="788" he="245" /></maths>其中,坐标变换矩阵A<sub>z→z′</sub>=K(γ)L(θ)M(ψ),K(γ)、L(θ)、M(ψ)分别为三步变换的变换矩阵,它们的具体表达式为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>K</mi><mrow><mo>(</mo><mi>&gamma;</mi><mo>)</mo></mrow><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mi>cos</mi><mi>&gamma;</mi></mtd><mtd><mi>sin</mi><mi>&gamma;</mi></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mo>-</mo><mi>sin</mi><mi>&gamma;</mi></mtd><mtd><mi>cos</mi><mi>&gamma;</mi></mtd></mtr></mtable></mfenced><mo>,</mo><mi>L</mi><mrow><mo>(</mo><mi>&theta;</mi><mo>)</mo></mrow><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mi>cos</mi><mi>&theta;</mi></mtd><mtd><mi>sin</mi><mi>&theta;</mi></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mo>-</mo><mi>sin</mi><mi>&theta;</mi></mtd><mtd><mi>cos</mi><mi>&theta;</mi></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd></mtr></mtable></mfenced><mo>,</mo><mi>M</mi><mrow><mo>(</mo><mi>&psi;</mi><mo>)</mo></mrow><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mi>cos</mi><mi>&psi;</mi></mtd><mtd><mn>0</mn></mtd><mtd><mo>-</mo><mi>sin</mi><mi>&psi;</mi></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mi>sin</mi><mi>&psi;</mi></mtd><mtd><mn>0</mn></mtd><mtd><mi>cos</mi><mi>&psi;</mi></mtd></mtr></mtable></mfenced><mo>;</mo></mrow>]]></math><img file="FDA0000568623910000012.GIF" wi="1838" he="252" /></maths>将三轴陀螺仪安装在相控阵天线背面,信号处理模块进行目标的测量,稳定控制模块进行波束稳定控制,信号处理周期为T<sub>s</sub>,波束稳定模块数据更新周期为T<sub>d</sub>,T<sub>s</sub>与T<sub>d</sub>满足T<sub>s</sub>/T<sub>d</sub>=N,其中N为整数,该方法包括以下步骤:第一步:相控阵雷达上电时,波束稳定模块各参量初始化,天线的姿态(ψ<sub>0</sub>,θ<sub>0</sub>,γ<sub>0</sub>)=(0,0,0)四元数Q<sub>0</sub>=[1 0 0 0],目标跟踪向量<img file="FDA0000568623910000013.GIF" wi="667" he="108" />ψ<sub>0</sub>为初始时刻偏航、θ<sub>0</sub>初始时刻俯仰、γ<sub>0</sub>初始时刻横滚;<img file="FDA0000568623910000014.GIF" wi="191" he="77" />分别为天线波束的初始时刻的方位角与俯仰角;第二步:利用t<sub>n‑1</sub>时刻四元数与t<sub>n</sub>时刻与t<sub>n‑1</sub>时刻时间间隔内的角速度测量值,通过数值积分法解算天线相对姿态;t<sub>n‑1</sub>时刻为t<sub>n</sub>时刻的上一测量时刻;第三步:根据上一步得到的天线相对姿态以及t<sub>n‑1</sub>时刻目标跟踪向量g<sub>n‑1</sub>,通过坐标变换反推得到保持天线波束稳定所需的波束补偿角<img file="FDA0000568623910000015.GIF" wi="548" he="78" />分别为t<sub>n</sub>时刻补偿平台扰动所需的方位角和俯仰角;为了使扰动后波束仍旧指向g<sub>n‑1</sub>,需要满足以下关系式:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>K</mi><mrow><mo>(</mo><msub><mi>&gamma;</mi><mi>n</mi></msub><mo>)</mo></mrow><mi>L</mi><mrow><mo>(</mo><msub><mi>&theta;</mi><mi>n</mi></msub><mo>)</mo></mrow><mi>M</mi><mrow><mo>(</mo><msub><mi>&psi;</mi><mi>n</mi></msub><mo>)</mo></mrow><mo>&times;</mo><msub><mi>g</mi><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>=</mo><mi>M</mi><mrow><mo>(</mo><mo>-</mo><msub><mi>&lambda;</mi><mi>y</mi></msub><mo>)</mo></mrow><mi>L</mi><mrow><mo>(</mo><mo>-</mo><msub><mi>&lambda;</mi><mi>z</mi></msub><mo>)</mo></mrow><mo>&times;</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mn>1</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd></mtr></mtable></mfenced><mo>;</mo></mrow>]]></math><img file="FDA0000568623910000016.GIF" wi="1006" he="233" /></maths>其中,ψ<sub>n</sub>为t<sub>n‑1</sub>时刻偏航、θ<sub>n</sub>为t<sub>n‑1</sub>时刻俯仰、γ<sub>n</sub>为t<sub>n‑1</sub>时刻横滚;λ<sub>y</sub>为天线波束的t<sub>n‑1</sub>时刻的方位角;λ<sub>z</sub>为天线波束的t<sub>n‑1</sub>时刻的俯仰角;令<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>&chi;</mi><mo>=</mo><mi>K</mi><mrow><mo>(</mo><msub><mi>&gamma;</mi><mi>n</mi></msub><mo>)</mo></mrow><mi>L</mi><mrow><mo>(</mo><msub><mi>&theta;</mi><mi>n</mi></msub><mo>)</mo></mrow><mi>M</mi><mrow><mo>(</mo><msub><mi>&psi;</mi><mi>n</mi></msub><mo>)</mo></mrow><mo>&times;</mo><msub><mi>g</mi><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>&chi;</mi><mn>1</mn></msub></mtd></mtr><mtr><mtd><msub><mi>&chi;</mi><mn>2</mn></msub></mtd></mtr><mtr><mtd><msub><mi>&chi;</mi><mn>3</mn></msub></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000568623910000021.GIF" wi="742" he="235" /></maths>则<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>&lambda;</mi><mi>y</mi></msub><mo>=</mo><mi>arctan</mi><mrow><mo>(</mo><mo>-</mo><mfrac><msub><mi>&chi;</mi><mn>3</mn></msub><msub><mi>&chi;</mi><mn>1</mn></msub></mfrac><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000568623910000022.GIF" wi="395" he="164" /></maths>λ<sub>z</sub>=arcsin(‑χ<sub>2</sub>);当前时刻抵消阵面扰动所需的波束补偿角为<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mfenced open='{' close=''><mtable><mtr><mtd><mi>&Delta;</mi><msub><mi>&lambda;</mi><msub><mi>y</mi><mi>n</mi></msub></msub><mo>=</mo><msub><mi>&lambda;</mi><mi>y</mi></msub><mo>-</mo><msub><mi>&lambda;</mi><msub><mi>y</mi><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></msub></msub></mtd></mtr><mtr><mtd><mi>&Delta;</mi><msub><mi>&lambda;</mi><msub><mi>z</mi><mi>n</mi></msub></msub><mo>=</mo><msub><mi>&lambda;</mi><mi>z</mi></msub><mo>-</mo><msub><mi>&lambda;</mi><msub><mi>z</mi><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></msub></msub></mtd></mtr></mtable></mfenced><mo>;</mo></mrow>]]></math><img file="FDA0000568623910000023.GIF" wi="385" he="177" /></maths>第四步:获取目标测量值,若t<sub>n</sub>时刻为信号处理帧信号到来时刻,通过信号处理相关算法和相参积累,得到t<sub>n</sub>时刻目标方位误差角ε<sub>y</sub>、俯仰误差角ε<sub>z</sub>;否则,ε<sub>y</sub>=0、ε<sub>z</sub>=0;第五步:利用计算得到的波束补偿角<img file="FDA0000568623910000024.GIF" wi="281" he="78" />目标方位误差角ε<sub>y</sub>和俯仰误差角ε<sub>z</sub>共同实现波束控制,并更新目标跟踪向量、四元数用于算法迭代,即:若目标测量值不为0,则<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msub><mi>g</mi><mi>n</mi></msub><mo>=</mo><mi>M</mi><mrow><mo>(</mo><mo>-</mo><msub><mi>&psi;</mi><mi>n</mi></msub><mo>)</mo></mrow><mi>L</mi><mrow><mo>(</mo><mo>-</mo><msub><mi>&theta;</mi><mi>n</mi></msub><mo>)</mo></mrow><mo>&CenterDot;</mo><mi>K</mi><mrow><mo>(</mo><mo>-</mo><msub><mi>&gamma;</mi><mi>n</mi></msub><mo>)</mo></mrow><mi>M</mi><mrow><mo>(</mo><mo>-</mo><msub><mi>&lambda;</mi><msub><mi>y</mi><mi>n</mi></msub></msub><mo>)</mo></mrow><mi>L</mi><mrow><mo>(</mo><mo>-</mo><msub><mi>&lambda;</mi><msub><mi>z</mi><mi>n</mi></msub></msub><mo>)</mo></mrow><mo>&CenterDot;</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mn>1</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd></mtr></mtable></mfenced><mo>;</mo></mrow>]]></math><img file="FDA0000568623910000025.GIF" wi="1045" he="236" /></maths>否则,g<sub>n</sub>=g<sub>n‑1</sub>,其中g<sub>n</sub>表示为t<sub>n</sub>时刻目标跟踪向量,从而抵消t<sub>n</sub>时刻与t<sub>n‑1</sub>时刻时间间隔内载体扰动对波束的影响,保持波束始终正对目标,相控阵雷达天线波束控制角度为:<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>&lambda;</mi><msub><mi>y</mi><mi>n</mi></msub></msub><mo>=</mo><msub><mi>&lambda;</mi><msub><mi>y</mi><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></msub></msub><mo>+</mo><mi>&Delta;</mi><msub><mi>&lambda;</mi><msub><mi>y</mi><mi>n</mi></msub></msub><mo>+</mo><msub><mi>&epsiv;</mi><mi>y</mi></msub></mtd></mtr><mtr><mtd><msub><mi>&lambda;</mi><msub><mi>z</mi><mi>n</mi></msub></msub><mo>=</mo><msub><mi>&lambda;</mi><msub><mi>z</mi><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></msub></msub><mo>+</mo><mi>&Delta;</mi><msub><mi>&lambda;</mi><msub><mi>z</mi><mi>n</mi></msub></msub><mo>+</mo><msub><mi>&epsiv;</mi><mi>z</mi></msub></mtd></mtr></mtable></mfenced><mo>;</mo></mrow>]]></math><img file="FDA0000568623910000026.GIF" wi="493" he="194" /></maths><img file="FDA0000568623910000027.GIF" wi="190" he="78" />分别为t<sub>n</sub>时刻天线需要指向的方位角和俯仰角,<img file="FDA0000568623910000028.GIF" wi="252" he="78" />分别为t<sub>n‑1</sub>时刻天线所处的方位角和俯仰角。
地址 100081 北京市海淀区中关村南大街5号