网刊加载中。。。

使用Chrome浏览器效果最佳,继续浏览,你可能不会看到最佳的展示效果,

确定继续浏览么?

复制成功,请在其他浏览器进行阅读

双频测速应答机灵敏度遥测参数异常下跳原因分析  PDF

  • 郭忠来
  • 李红梅
  • 李婷
  • 马向斌
  • 李晶
中国酒泉卫星发射中心,兰州,732750

中图分类号: TN8

发布日期:2019-10-11

DOI:10.16356/j.1005-2615.2019.增刊.013

  • 全文
  • 图表
  • 参考文献
  • 作者
  • 出版信息
目录contents

摘要

针对某航天发射箭载测控合作目标双频测速应答机灵敏度遥测参数异常下跳问题,深入研究了连续波无线电外测系统的高精度测速原理,根据实际飞行弹道计算出地面连续波跟踪测量雷达系统4个测点相对于双频测速应答机接收天线的方向图,得出在灵敏度遥测参数异常下跳期间,由于迭代制导,火箭姿态变化剧烈,使得各测点的赤道面张角α变化范围较大,大大超出了双频测速应答机接收天线的增益要求,出现接收天线接收的信号变弱从而导致灵敏度在此时段异常下跳的结论。本文可为后续航天发射飞行中出现类似问题时的故障快速定位提供分析思路。

在历次CZ‑2D火箭飞行任务中,无论发射方位角如何变化,双频测速应答机灵敏度遥测参数总是随着测控距离的增加缓慢下降,但在某次发射飞行中,火箭双频测速应答机灵敏度遥测参数却在火箭飞行328.1~332.2 s期间出现了异常下跳。

连续波无线电外测系统的功能是由地面连续波跟踪测量雷达系统的发射机产生无线电信号,通过天线发向火箭的双频测速应答机,双频测速应答机接收后经过运算处理,输出相参下行信号,通过发射天线发送到地面连续波雷达,地面雷达经过接收处理,计算出火箭的飞行速度和飞行弹道,为火箭飞行监控和实时安控提供数据源,以及为光学、遥测等地面测量设备提供引导弹道数[

1]。因此,对于双频测速应答机灵敏度遥测参数异常变化的分析对于火箭飞行安全具有重要意义。

1 双频测速应答机的工作原理

双频测速应答机与地面站测控设备协同工作,对火箭进行连续跟踪测量,获取高精度测速数据。

图1为双频测速应答机其中一个通道的工作原理框图。当没有上行信号(频点1或频点2)时,载波环路(频点1或频点2)处于失锁状态,应答和信标切换电路将该路(频点1或频点2)置于信标状态,信标信号由恒温晶振产生参考信号,经过倍频链、滤波、放大输出非相参下行信号,此时,接收机的一本振和二本振处于扫描状态。当接收机接收到上行信号时,上行信号经过低噪声微波放大、下混频、前置放大、变频、中频放大后,进入载波同步环路,当环路锁定时,锁定指示输出高电平,将信号置于应答跟踪状态,同步后相参信号经过倍频、滤波、放大后,输出相参下行信号通过发射机发送到地面连续波雷[2‑3]

图1 双频测速应答机单通道原理框图

Fig.1 Single channel principle block diagram of dual‑frequency speed transponder

当双频测速应答机处于跟踪上行信号的斜锁相工作状态时,对于F1频率情况,接收频率为Fr1=x735 MHz,发送频率Ft1=x318.75 MHz,第一收发频率之间相差一个相参的中频信号416.25 MHz,收发频率由系统总体选定。考虑到方案的简单与避免组合干扰的影响,第一混频器的本振由主环VCXO频率倍频115次提供,第二混频器的本振由主环VCXO频率倍频8次提供,则前置中放频率设定为Fr1‑115FVCXO,中频放大频率设定为FVCXO,主环完成对接收信号Fr1的频率跟踪,故环路FVCXO频率与接收频率Fr1是相参的。

同样对于F2频率情况,接收频率为Fr2=x697.9 MHz,发送频率Ft2=x259.6 MHz,第一收发频率之间差一个相参的中频信号438.3 MHz,收发频率由系统总体选定。考虑到方案的简单与避免组合干扰的影响,第一混频器的本振由主环VCXO频率倍频108次提供,第二混频器的本振由主环VCXO频率倍频8次提供,则前置中放频率设定为Fr2‑108FVCXO,中频放大频率设定为FVCXO,主环完成对接收信号Fr2的频率跟踪,故环路FVCXO频率与接收频率Fr2是相参的。

双频测速应答机在结构上由5个独立的单元组成,即接收锁相组合1单元、接收锁相组合2单元、本振信标组合1单元、本振信标组合2单元、电源变换功放组合单元。设备加电后,由28VDC电源经过二次电源产生发射机、接收机、载波同步环路、应答和切换电路所需的各种电源。

2 连续波雷达测量原理

地面连续波跟踪测量雷达系[

4]发射连续波射频信号对空中目标进行跟踪测量,其特点是容易实现测速和信道的综合利用,测量精度较高。连续波测速方法最终归结为多普勒频移的测量,再由多普勒频移换算出径向速度。

2.1 双程多普勒频移测速

当箭上的双频测速应答机某一通道处于应答状态时,地面连续波跟踪测量雷达系统高稳定度地发射频率为f0的正弦波信号,双频测速应答机锁定f0并向地面转发以实现上下频率相干。

双频测速应答机接收到的信号时延为

τd=Rc (1)

则接收信号相位为

φR=2πf0(t-τ)=2πf0(t-Rc) (2)

f0恒定时,双频测速应答机接收频率为

f1=12πdφRdt=f0(1-R˙c) (3)

地面连续波跟踪测量雷达系统接收频率为

f2=f1(1-R˙c)=f0(1-R˙c)2 (4)

故双程多普勒频率为

fd=f3-f0=-2f0R˙c (5)

式中: c为光速;R˙为径向距离变化率,即径向速度Vd

2.2 多普勒频率的测量

目标回波的多普勒频移与其径向速度成正比,因此只要准确地测出目标的多普勒频移,经换算就可以确定目标运动的径向速度和方[

5],即

vd=-c2f0fd (6)

要完成多普勒频率的测量,首先必须从包含噪声的接收载波中提取多普勒频率。接收载波经过快速傅里叶变换(FFT),可以完成多普勒频率的粗略捕获。随着技术的发展,应用DSP技术的数字锁相环技[

6]愈加成熟,它可以在过滤多普勒频率信息过程中,提取数字化的多普勒频率信息,完成多普勒频率的精确测量。

2.3 高精度测速处[7‑8]

目标在发射系的位置X,Y,Z需通过光测交会或脉冲雷达定位进行计算,使用当前所有状态可用的主站和副站的测速数据,一般选择主站和3个副站中任选2个副站的测速数据,即可完成发射系的速度VxVyVz的测量。

(1)各测站到目标连线的方向余弦计算

limini=1RiX-XoiY-YoiZ-Zoi (7)
Ri=(X-Xoi)2+(Y-Yoi)2+(Z-Zoi)2 (8)

式中:XoiYoiZoi表示各站站址,其中i=1表示当前所用主站发射天线,i =2,3,4,…时表示当前状态可用的副站收天线站址。

(2)矩阵计算

A=l1+l2m1+m2n1+n2l1+l3m1+m3n1+n3.........l1+lim1+min1+ni (9)

(3)速度计算

X˙,Y˙,Z˙T=(ATA)-1ATS1,S2,,SiT (10)

式中:S1S2,…,Si表示当前状态可用的主站和副站的测速数据。

新一代连续波测速设备LD‑3704采用70 MHz中频数字锁相环处理方案,充分利用DSP和FPGA的技术优势,实现数字测速单元及数字角误差解调单元的有机结合,降低设备量,改善系统性能,提高设备可靠性。

测速信息流程如图2所示。图中发射分系统通过频综Ⅰ以时统频标为基准产生上行信号。应答机接收到发射机辐射的上行信号后,相参转发下行信号。接收分系统接收的下行信号与一,二本振混频产生70 MHz中频信号。数字基带分系统进行载波跟踪和窄带滤波,得到载波的相位增量或平均多普勒。最后进行记录存盘,以便事后处理;同时通过数传、通信链路送往指挥中心,完成速度测量任务。

图2 测速信息流程图

Fig.2 Flow chart of speed measurement information

3 双频测速应答机灵敏度遥测参数异常下跳原因分析

3.1 问题现象

图3所示,双频测速应答机灵敏度1遥测参数曲线平稳缓慢的下降是正常的,而在2018年某任务中,灵敏度在328.1~332.2 s由4.38 V迅速下跳至3.26 V,后又在366.6 s缓慢上升到4.11 V,与以往任务规律不相符。

图3 双频测速应答机灵敏度遥测参数曲线

Fig.3 Sensitivity telemetry parameter curve of dual‑frequency speed transponder

3.2 原因分析

3.2.1 双频测速应答机故障

虽然该火箭双频测速应答机的灵敏度遥测参数在328.1~ 366.6 s期间迅速下降,但是功率遥测和主环锁定遥测两个参数正常,如图4所示。

图4 双频测速应答机通道1遥测参数曲线

Fig.4 Telemetry parameter curve of channel 1 in dual‑frequency speed transponder

双频测速应答机灵敏度遥测参数反映的是双频测速应答机接收地面连续波雷达系统发射的连续波射频信号强弱的物理量。随着火箭飞行距离的增大,地面连续波雷达系统离火箭越来越远,双频测速应答机接收到的连续波射频信号越来越弱,灵敏度遥测参数从而越来越小。主环锁定遥测参数表征的是当接收信号进入载波同步环路后,当环路锁定时,锁定指示输出高电平,将信号置于应答跟踪状态,同步后相参信号经过倍频、滤波、放大后,输出相参下行信号通过发射机发送到地面连续波雷达。功率遥测参数表征的是双频测速应答机的工作电压。

由此可以断定双频测速应答机工作正常,不是引起灵敏度遥测参数下跳的原因。

3.2.2 地面连续波跟踪测量雷达系统故障

火箭飞行过程中,地面连续波跟踪测量雷达系统在跟踪至327~370 s期间发现设备接收电平有异常下降但恢复的情况,但没影响跟踪和数据质量。各设备通道一接收电平、多普勒频率三阶差分曲线分别如图5,6,7,8所示。

图5 主站接收电平

Fig.5 Main station receiving level

图6 副站1接收电平

Fig.6 Secondary 1 receiving level

图7 副站2接收电平

Fig.7 Secondary 2 receiving level

图8 副站3接收电平

Fig.8 Secondary 3 receiving level

从事后数据分析可知:各站跟踪情况(接收电平和测速随机误差)一致良好;各站接收电平平稳变化,在327~370 s出现“V”型下降期间,接收电平仍满足设备正常跟踪条件。

由此可见,地面连续波跟踪测量雷达系统不是引起灵敏度遥测参数下跳的原因。

3.2.3 无线传输信道故障

为了分析地面连续波跟踪测量雷达系统和箭上双频测速应答机在灵敏度遥测参数异常下跳期间无线传输信道是否畅通,需要根据实际飞行弹道和地面测站的布局计算双频测速应答机收发天线的方向性图,看实际计算的结果能否满足双频测速应答机收发天线的增益要求。

(1)天线方向性图定义

在箭体坐标系中天线增益的方向性以两个参考角αβ描述,如图9所示。所谓天线方向性图是对天线辐射强度特性(增益)的一种描[

9]

图9 天线方向性与火箭箭体坐标系关联示意图

Fig.9 Schematic diagram of the correlation between the antenna directivity and the coordinate systemof rocket body

图9中,α俗称赤道面张角,以火箭的第Ⅲ象限线(箭体坐标系的y轴)为基准,沿飞行方向观察,顺时针方向为角度增加方向,变化范围为0°~360°;β俗称子午面张角,以火箭的轴线(箭体坐标系的x轴)为基准,从头部向尾部旋转为角度增加方向,变化范围为0°~180°,箭首为0°,箭尾方向为180°[

10]αm为天线在火箭横截面上的安装角,接收天线安装角为187.5°,发射天线安装角为167.5°。

Δα=α-αm (11)

在箭体坐标系中天线增益的分布为

G=f(α,β) (12)

(2)天线方向性图计算方法

地面高精度测速系统4个测量点的箭体坐标为

xjtyjtzjt=FjtTxf-xyf-yzf-z (13)

式中:xjtyjtzjt表示测量点的箭体坐标系坐标;xyz表示箭体坐标系原点的发射坐标系坐标,选用理论弹道参数; xfyfzf表示测量点的发射坐标系坐标,其值为

xfyfzf=F0TXe-X0Ye-Y0Ze-Z0 (14)

式中X0Y0Z0表示射点的地心坐标,计算公式为

X0Y0Z0=(N0+h0)cosB0cosL0(N0+h0)cosB0sinL0[N01-e2+h0]sinB0 (15)
N0=a1-e2sin2B0 (16)

式中:L0B0h0分别表示发射点大地经度、纬度和大地高;ae2分别为地球参考椭球长半径和第一偏心率平方;XeYeZe为测量点的地心坐标,计算方法同上,将测量点的大地经度、纬度和大地高LeBehe替换L0B0h0即可;F0为发射坐标系至地心坐标系的旋转变换矩阵,计算方法为

F0=R(90-L0)R(-B0)R(90+A0) (17)

式中:A0为发射大地方位角;Fjt为发射坐标系至箭体坐标系的旋转变换矩阵,计算方法为

Fjt=Rx(γ)Ry(ψ)Rz(ϕ) (18)

式中ΦΨγ分别为箭体对发射坐标系的俯仰角、偏航角和滚动角。

天线方向性图用箭体球坐标表述,箭体直角坐标至箭体球坐标的转换方法为

R=xjt2+yjt2+zjt2 (19)
a=arctanzjtyjt+0       yjt>0,zjt0π       yjt02π    yjt>0,zjt<0 (20)
β=arctanyjt2+zjt2xjt+0    xjt>0π    xjt0 (21)

(3)天线方向性图计算结果

根据上述计算方法,编程计算出4个测量点相对于双频测速应答机的收发天线方向图,赤道面张角α随时间变化曲线如图10所示,子午面张角β随时间变化曲线如图11所示。

图10 赤道面张角曲线

Fig.10 Equatorial plane angle‑opening curve

图11 子午面张角曲线

Fig.11 Meridian plane angle‑opening curve

图10图11可以看出,在328.1~366.6 s期间,主站和2个副站的α变化范围是189.8°~6.3°,对应的Δα为2.3°~-181.2°;1个副站的α变化范围是199.1°~348.7°,对应的Δα为11.6°~161.2°;β变化范围是173.3°~177.7°(箭尾方向)。

(4)原因分析

经分析,火箭在327.3 s迭代制导开始,在329.3 s控制系统箭机发出-1.55°的偏航程序角,同时发出+9°的俯仰程序角(3°/s的速率,持续3 s),如图12所示。

图12 控制系统程序角

Fig.12 Program angle of the control system

由于火箭姿态变化剧烈,使得双频测速应答机收发天线方向图中的赤道面张角α变化范围较大,查阅《双频连续波应答机收发天线增益方向图测试报告》,当Δα在-80°~90°的范围内变化时,接收天线(工作频点为x735 MHz)在箭尾方向的增益满足要求,接收天线能够收到地面雷达发送的信号,当Δα大大超出了接收天线的增益要求时,就会使得接收天线接收的信号变弱从而导致接收灵敏度在此时段异常下跳。

4 结  论

针对某次火箭发射飞行中双频测速应答机灵敏度遥测参数异常下跳问题,通过对箭上双频测速应答机和地面连续波跟踪测量雷达系统的工作机理分析,以及根据实际飞行弹道计算地面各测点相对于双频测速应答机接收天线的方向图,得出在灵敏度遥测参数异常下跳期间,由于火箭姿态变化剧烈,各测点的赤道面张角α变化范围较大,大大超出了双频测速应答机接收天线的增益要求,使得接收天线接收的信号变弱从而导致接收灵敏度在此时段异常下跳。

另外层层递进、测发测控联合的分析方法可以类推解决其他遥测参数异常变化的问题,增强基地测发和测控两大系统的知识融合,有较高的实际应用价值,可在今后火箭发射飞行类似故障分析和定位中推广应用。

参考文献

1

于志坚. 航天测控系统工程[M]. 北京: 国防大学出版社, 2008.

YU Zhijian. Space TT&C system engineering[M]. Beijing: National Defense University Press, 2008.

2

李孟祥. 运载火箭外测与安全系统[M]. 北京: 国防工业出版社, 2004.

LI Mengxiang. Launch vehicle outboard and safety system[M]. Beijing: National Defense University Press, 2004.

3

成求青. 导弹测控系统总体设计原理与方法[M]. 北京:清华大学出版社, 2014.

CHENG Qiuqing. Principle and method of overall design for missile measurement and control system[M]. Beijing: Tsinghua University Press, 2014.

4

张明友汪学刚. 雷达系统[M]. 北京: 电子工业出版社, 2005.

ZHANG Mingyou, WANG Xuegang. Radar system[M]. Beijing: Electronic Industry Press, 2005.

5

张琳, 盛卫星, 马晓峰. 基于微多普勒效应的雷达目标识别算法[J]. 现代雷达200712): 3538.

ZHANG Lin, SHENG Weixing, MA Xiaofeng. Radar target recognition algorithm based on micro‑Doppler effect[J]. Modern Radar, 2007(12): 3538.

6

方立军, 马骏. 基于数字锁相的雷达频率合成器的研究[J]. 现代电子, 20013): 2125.

FANG Lijun, MA Jun. Research on radar frequency synthesizer based on digital phase lock[J]. Modern Electronics, 2001(3): 2125.

7

江衍煊. 浅析提高雷达精度的数字设计[J]. 科技创新与应用, 2012(29):25.

JIANG Yanxuan. Analysis of digital design for improving radar accuracy[J]. Science and Technology Innovation and Application, 2012(29):25.

8

郭军海, 吴正容, 黄学德. 多测速雷达弹道测量体制研究[J]. 飞行器测控学报, 2002, 21(3) : 511.

GUO Junhai, WU Zhengrong, HUANG Xuede. Research on trajectory measurement system of multi‑velocity radar[J]. Journal of Spacecraft TT & C Technology, 2002, 21(3): 511.

9

陈志杰李永祯戴幻尧. 相控阵天线方向图的建模与实时仿真方法[J]. 计算机仿真, 201128(3) : 3135.

CHEN Zhijie, LI Yongzhen, DAI Huanyao. Modeling and real‑time simulation of phased array antenna pattern[J]. Computer Simulation, 2011, 28(3): 3135.

10

刘蕴才.无线电遥测遥控[M]. 北京国防工业出版社2001.

LIU Yuncai. Radio telemetry and remote control[M]. Beijing: National Defence Industrial Press, 2001.

您是第位访问者
网站版权 © 南京航空航天大学学报
技术支持:北京勤云科技发展有限公司
请使用 Firefox、Chrome、IE10、IE11、360极速模式、搜狗极速模式、QQ极速模式等浏览器,其他浏览器不建议使用!