发明名称 一种双基雷达固定目标时域定位方法
摘要 本发明公开了一种双基雷达固定目标时域定位方法:步骤1,建立双基定位系统模型;步骤2,推导单散射点目标的斜距历程;步骤3,推导解调后单散射点目标回波信号模型,并将其离散化;步骤4,构造三维信号匹配矩阵;步骤5,通过回波信号与信号匹配矩阵的运算来估计目标位置,实现目标定位。本发明采用了双基雷达模式,生存能力较强;提出的定位系统的定位精度不依赖于信号带宽,可以采用窄带信号甚至是单频连续波信号,对发射机要求不高,降低了整个系统的成本;定位算法运算复杂度较低,鲁棒性高,有较高的定位精度,可以同时测得观察区域内的多个目标位置。
申请公布号 CN105137417A 申请公布日期 2015.12.09
申请号 CN201510503917.8 申请日期 2015.08.17
申请人 周口师范学院 发明人 刘玉春;朱海;程全;姚遥;张利红;程琳;李晋;豆桂平
分类号 G01S13/06(2006.01)I 主分类号 G01S13/06(2006.01)I
代理机构 郑州中原专利事务所有限公司 41109 代理人 霍彦伟
主权项 一种双基雷达固定目标时域定位方法,其特征在于包括以下步骤:步骤1,建立双基定位系统模型:在二维空间内,目标固定尺寸远小于目标至接收机的距离,并且小于雷达波束宽度,可以看作点目标;接收机和发射机分别绕着各自的转动中心对应做半径为r<sub>1</sub>、r<sub>2</sub>的匀速圆周运动,其角速度分别为ω<sub>1</sub>和ω<sub>2</sub>,并且有ω<sub>1</sub>≠ω<sub>2</sub>;以接收机的转动中心为原点,以接收机转动中心到发射机转动中心的方向为X轴正向在二维空间内建立笛卡尔坐标系;步骤2,推导单散射点目标的斜距历程R得到:R(t)≈R<sub>r0</sub>+R<sub>t0</sub>‑r<sub>1</sub>cos(ω<sub>1</sub>t+α<sub>1</sub>‑θ<sub>1</sub>)‑r<sub>2</sub>cos(ω<sub>2</sub>t+α<sub>2</sub>‑θ<sub>2</sub>)其中,R<sub>t0</sub>为发射机旋转中心到目标的距离,R<sub>r0</sub>为接收机转动中心到目标的距离,θ<sub>1</sub>为目标与接收机转动中心的连线与X轴正向的夹角,θ<sub>2</sub>为目标与发射机转动中心的连线与X轴正向的夹角,θ(t)=θ<sub>0</sub>+ωt,θ<sub>0</sub>为初始角度,ω为接收机旋转角速度,α<sub>1</sub>和α<sub>2</sub>分别是接收机和发射机在零时刻的初始方位角;步骤3,推导解调后单散射点目标回波信号模型,得到:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>s</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mi>&sigma;</mi><mi>exp</mi><mrow><mo>{</mo><mrow><mo>-</mo><mi>j</mi><mfrac><mrow><mn>2</mn><mi>&pi;</mi></mrow><mi>&lambda;</mi></mfrac><mrow><mo>(</mo><mrow><msub><mi>R</mi><mrow><mi>t</mi><mn>0</mn></mrow></msub><mo>+</mo><msub><mi>R</mi><mrow><mi>r</mi><mn>0</mn></mrow></msub></mrow><mo>)</mo></mrow></mrow><mo>}</mo></mrow><mi>exp</mi><mrow><mo>{</mo><mrow><mi>j</mi><mfrac><mrow><mn>2</mn><mi>&pi;</mi></mrow><mi>&lambda;</mi></mfrac><msub><mi>r</mi><mn>1</mn></msub><mi>cos</mi><mrow><mo>(</mo><mrow><msub><mi>&omega;</mi><mn>1</mn></msub><mi>t</mi><mo>+</mo><msub><mi>&alpha;</mi><mn>1</mn></msub><mo>-</mo><msub><mi>&theta;</mi><mn>1</mn></msub></mrow><mo>)</mo></mrow></mrow><mo>}</mo></mrow><mi>exp</mi><mrow><mo>{</mo><mrow><mi>j</mi><mfrac><mrow><mn>2</mn><mi>&pi;</mi></mrow><mi>&lambda;</mi></mfrac><msub><mi>r</mi><mn>2</mn></msub><mi>cos</mi><mrow><mo>(</mo><mrow><msub><mi>&omega;</mi><mn>2</mn></msub><mi>t</mi><mo>+</mo><msub><mi>&alpha;</mi><mn>2</mn></msub><mo>-</mo><msub><mi>&theta;</mi><mn>2</mn></msub></mrow><mo>)</mo></mrow></mrow><mo>}</mo></mrow></mrow>]]></math><img file="FDA0000782956780000011.GIF" wi="1922" he="150" /></maths>其中,σ为目标散射强度,λ为信号波长;假设在定位过程中共采样K次,则回波信号离散化形式为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>s</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mi>&sigma;</mi><mi>exp</mi><mo>{</mo><mo>-</mo><mi>j</mi><mfrac><mrow><mn>2</mn><mi>&pi;</mi></mrow><mi>&lambda;</mi></mfrac><mrow><mo>(</mo><msub><mi>R</mi><mrow><mi>t</mi><mn>0</mn></mrow></msub><mo>+</mo><msub><mi>R</mi><mrow><mi>r</mi><mn>0</mn></mrow></msub><mo>)</mo></mrow><mo>}</mo><mi>exp</mi><mo>{</mo><mi>j</mi><mfrac><mrow><mn>2</mn><mi>&pi;</mi></mrow><mi>&lambda;</mi></mfrac><msub><mi>r</mi><mn>1</mn></msub><mi>c</mi><mi>o</mi><mi>s</mi><mrow><mo>(</mo><msub><mi>k&Delta;&omega;</mi><mn>1</mn></msub><mo>+</mo><msub><mi>&alpha;</mi><mn>1</mn></msub><mo>-</mo><msub><mi>&theta;</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>}</mo><mi>exp</mi><mo>{</mo><mi>j</mi><mfrac><mrow><mn>2</mn><mi>&pi;</mi></mrow><mi>&lambda;</mi></mfrac><msub><mi>r</mi><mn>2</mn></msub><mi>c</mi><mi>o</mi><mi>s</mi><mrow><mo>(</mo><msub><mi>k&Delta;&omega;</mi><mn>2</mn></msub><mo>+</mo><msub><mi>&alpha;</mi><mn>2</mn></msub><mo>-</mo><msub><mi>&theta;</mi><mn>2</mn></msub><mo>)</mo></mrow><mo>}</mo></mrow>]]></math><img file="FDA0000782956780000012.GIF" wi="2031" he="151" /></maths>其中,k=1,2,…,K;步骤4,构造三维信号匹配矩阵:假设用β<sub>1</sub>搜索θ<sub>1</sub>,用β<sub>2</sub>搜索θ<sub>2</sub>,在β<sub>1</sub>方向将(0,2π]分成M份,则搜索步长为Δβ<sub>1</sub>=2π/M,在方向β<sub>2</sub>将(0,2π]分成N份,则搜索步长为Δβ<sub>2</sub>=2π/N,根据离散化的回波信号形式可得三维信号匹配矩阵为:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>g</mi><mrow><mo>(</mo><mi>m</mi><mo>,</mo><mi>n</mi><mo>,</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mi>exp</mi><mo>{</mo><mo>-</mo><mi>j</mi><mfrac><mrow><mn>2</mn><mi>&pi;</mi></mrow><mi>&lambda;</mi></mfrac><msub><mi>r</mi><mn>1</mn></msub><mi>c</mi><mi>o</mi><mi>s</mi><mrow><mo>(</mo><msub><mi>k&Delta;&omega;</mi><mn>1</mn></msub><mo>+</mo><msub><mi>&alpha;</mi><mn>1</mn></msub><mo>-</mo><msub><mi>m&Delta;&beta;</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>}</mo><mi>exp</mi><mo>{</mo><mo>-</mo><mi>j</mi><mfrac><mrow><mn>2</mn><mi>&pi;</mi></mrow><mi>&lambda;</mi></mfrac><msub><mi>r</mi><mn>2</mn></msub><mi>c</mi><mi>o</mi><mi>s</mi><mrow><mo>(</mo><msub><mi>k&Delta;&omega;</mi><mn>2</mn></msub><mo>+</mo><msub><mi>&alpha;</mi><mn>2</mn></msub><mo>-</mo><msub><mi>n&Delta;&beta;</mi><mn>2</mn></msub><mo>)</mo></mrow><mo>}</mo></mrow>]]></math><img file="FDA0000782956780000013.GIF" wi="1810" he="148" /></maths>其中,m=1,2,…,M;n=1,2,…,N;k=1,2,…,K;步骤5,将回波信号s(k)中的第k个元素与三维信号匹配矩阵g(m,n,k)中的第k页中所有元素相乘来估计方位角对θ<sub>1</sub>和θ<sub>2</sub>,得到对方位角对θ<sub>1</sub>和θ<sub>2</sub>的估计分别为β<sub>1</sub>和β<sub>2</sub>,进而得到对目标位置的估计:<img file="FDA0000782956780000021.GIF" wi="696" he="162" />实现目标定位,其中,L为接收机转动中心到发射机转动中心之间的距离。
地址 466000 河南省周口市川汇区文昌东路周口师范学院
您可能感兴趣的专利