发明名称 一种基于SAR卫星姿态控制实现滑动聚束模式的方法
摘要 本发明公开了一种基于SAR卫星姿态控制实现滑动聚束模式的方法,属于信号处理技术领域。所述方法通过对卫星姿态进行调控,实现卫星天线指向角连续变化,不存在角度量化问题,通过建立较完善的空间关系几何模型,可更精确地实现滑动聚束模式,提高SAR方位向分辨率。而且,天线结构简单,可以降低卫星成本,减小卫星体积和重量,具有广阔的应用前景。本发明提出的方法相比于目前通过相控阵实现滑动聚束模式的方法,天线指向连续变化,准确度更高;可以大大简化卫星天线结构,降低卫星成本和减轻卫星重量,具有十分重要的应用价值。
申请公布号 CN103076607B 申请公布日期 2014.07.30
申请号 CN201310001124.7 申请日期 2013.01.04
申请人 北京航空航天大学 发明人 陈杰;邹德意;王鹏波;朱燕青
分类号 G01S13/90(2006.01)I;G01S7/02(2006.01)I 主分类号 G01S13/90(2006.01)I
代理机构 北京永创新实专利事务所 11121 代理人 姜荣丽
主权项 一种基于SAR卫星姿态控制实现滑动聚束模式的方法,其特征在于,包括以下几个步骤: 步骤一:计算卫星在转动地心坐标系E<sub>g</sub>中的坐标(x<sub>s</sub>,y<sub>s</sub>,z<sub>s</sub>); 发射第m个脉冲时,平均近心角M为: <img file="FDA0000512254140000011.GIF" wi="1284" he="138" />进而计算偏心角E为: <img file="FDA0000512254140000012.GIF" wi="1549" he="318" />计算真近心角θ为: <img file="FDA0000512254140000013.GIF" wi="1223" he="166" />计算矢 径r为: <img file="FDA0000512254140000015.GIF" wi="1115" he="130" />由上可得卫星在轨道平面坐标系E<sub>v</sub>中的坐标(x<sub>vs</sub>,y<sub>vs</sub>,z<sub>vs</sub>)为: <img file="FDA0000512254140000016.GIF" wi="1127" he="232" />通过坐标系转换,转换至转动地心坐标系E<sub>g</sub>,得到卫星在转动地心坐标系E<sub>g</sub>中的坐标(x<sub>s</sub>,y<sub>s</sub>,z<sub>s</sub>)为: <img file="FDA0000512254140000017.GIF" wi="1157" he="232" />从中取T<sub>0</sub>时刻卫星在转动地心坐标系E<sub>g</sub>的坐标为<img file="FDA0000512254140000018.GIF" wi="273" he="84" />步骤二:计算T<sub>0</sub>时刻天线坐标系E<sub>a</sub>中单位矢量(0,1,0)在转动地心坐标系E<sub>g</sub>中的坐标(x<sub>1</sub>,y<sub>1</sub>,z<sub>1</sub>); 通过坐标系转换,转换至卫星星体坐标系E<sub>e</sub>,转换至卫星平台坐标系E<sub>r</sub>,转换至轨道平面坐标系E<sub>v</sub>,转换至不转动的地心轨道坐标系E<sub>o</sub>,再转换至转动地心坐标系E<sub>g</sub>,得到天线坐标系E<sub>a</sub>中单位矢量(0,1,0)在转动地心坐标系E<sub>g</sub>中的坐标(x<sub>1</sub>,y<sub>1</sub>,z<sub>1</sub>)为: <img file="FDA0000512254140000021.GIF" wi="1301" he="231" />步骤三:计算T<sub>0</sub>时刻天线坐标系E<sub>a</sub>中单位矢量(0,1,0)到地面的距离R<sub>0</sub>; <img file="FDA0000512254140000022.GIF" wi="1185" he="78" /><img file="FDA0000512254140000023.GIF" wi="1273" he="98" /><img file="FDA0000512254140000024.GIF" wi="1221" he="85" />其中,r<sub>1</sub>,r<sub>2</sub>,r<sub>3</sub>分别为中间变量; 天线坐标系E<sub>a</sub>中单位矢量(0,1,0)到地面的距离R<sub>0</sub>: <img file="FDA0000512254140000025.GIF" wi="1151" he="167" />步骤四:计算T<sub>0</sub>时刻旋转点C到地面指向点的斜距ΔR<sub>0</sub>; 由滑动聚束方位向分辨率计算公式<img file="FDA0000512254140000026.GIF" wi="436" he="136" />得:<img file="FDA0000512254140000027.GIF" wi="1084" he="140" />步骤五:计算旋转点C在转动地心坐标系E<sub>g</sub>中的坐标<img file="FDA00005122541400000211.GIF" wi="278" he="83" /><img file="FDA0000512254140000028.GIF" wi="1172" he="267" />步骤六:计算全滑聚过程中卫星与旋转点C在转动地心坐标系E<sub>g</sub>中的相对矢量(R<sub>x</sub>,R<sub>y</sub>,R<sub>z</sub>); <img file="FDA0000512254140000029.GIF" wi="1076" he="264" />步骤七:计算全滑聚过程中天线指向在转动地心坐标系E<sub>g</sub>中的坐标(x<sub>1</sub>′,y<sub>1</sub>′,z<sub>1</sub>′); <img file="FDA00005122541400000210.GIF" wi="1217" he="550" />步骤八:计算全滑聚过程中天线指向在天线坐标系E<sub>a</sub>中的坐标(x<sub>1</sub>",y<sub>1</sub>",z<sub>1</sub>"); <img file="FDA0000512254140000031.GIF" wi="1323" he="246" />步骤九:计算天线坐标系E<sub>a</sub>下单位矢量(0,1,0),经姿态控制后的坐标(x<sub>ZT</sub>,y<sub>ZT</sub>,z<sub>ZT</sub>); 先控制俯仰角再控制偏航角得到姿态控制后的坐标(x<sub>ZT</sub>,y<sub>ZT</sub>,z<sub>ZT</sub>)为: <img file="FDA0000512254140000032.GIF" wi="1457" he="261" />化简得: <img file="FDA0000512254140000033.GIF" wi="1289" he="284" />步骤十:计算全滑聚过程中卫星姿态控制角:俯仰角θ<sub>p</sub>和偏航角θ<sub>y</sub>; <img file="FDA0000512254140000034.GIF" wi="1301" he="284" />解方程(20)得: <img file="FDA0000512254140000035.GIF" wi="1220" he="497" />并根据<img file="FDA0000512254140000036.GIF" wi="62" he="66" />得到<img file="FDA0000512254140000037.GIF" wi="64" he="70" />如下,<img file="FDA0000512254140000038.GIF" wi="541" he="167" />其中,<img file="FDA0000512254140000039.GIF" wi="398" he="326" />步骤十一:计算全滑聚过程中卫星姿态控制规律,及各轴角速度控制曲线; 以<img file="FDA00005122541400000310.GIF" wi="112" he="80" />分别表示顺序转轴矢量<img file="FDA00005122541400000311.GIF" wi="314" he="89" /><img file="FDA00005122541400000312.GIF" wi="318" he="89" /><img file="FDA00005122541400000313.GIF" wi="326" he="89" />则有:<img file="FDA0000512254140000041.GIF" wi="1448" he="649" />由此得姿态控制规律,x,y,z三轴的旋转角速度ω<sub>x</sub>,ω<sub>y</sub>,ω<sub>z</sub>分别为: <img file="FDA0000512254140000042.GIF" wi="1182" he="270" />经过以上十一个步骤,完成了对滑动聚束模式姿态控制规律的计算; 上述步骤中,字母参数的定义如下: T<sub>0</sub>为卫星正侧视波束中心指向观测目标区域中心的时刻;A<sub>num</sub>为满足观测带方位向长度要求,发射脉冲数;f<sub>PRF</sub>为脉冲重复频率PRF;θ<sub>L</sub>为雷达天线视角;ρ<sub>a</sub>为方位向分辨率;K<sub>a</sub>为方位向波束展宽因子;L<sub>s</sub>为等效天线长度;<img file="FDA0000512254140000043.GIF" wi="82" he="68" />为卫星在轨运动平均角速度;<img file="FDA0000512254140000044.GIF" wi="109" he="68" />为地球平均自转角速度;a为地球长半轴长度;b为地球短半轴长度;M为平均近心角;E为偏心角;θ为真近心角;e为偏心率;r为矢径;P为半正焦距;H<sub>G</sub>为春分点的格林威治时角;Ω为轨道升交点赤经;i为轨道倾角;ω为轨道近地点幅角;γ为卫星的航迹角;θ<sub>r</sub>为卫星绕x轴的滚转角;θ<sub>y</sub>为卫星绕y轴的偏航角;θ<sub>p</sub>为卫星绕z轴的俯仰角;θ<sub>y_c</sub>为偏航控制后的卫星偏航角;ω<sub>32</sub>为等效旋转角速度矢量;R<sub>2</sub>为卫星偏航时旋转矩阵;R<sub>3</sub>为卫星俯仰时旋转矩阵;用c表示余弦cos,用s表示正弦sin;<img file="FDA0000512254140000045.GIF" wi="462" he="250" />为不转动地心坐标系E<sub>o</sub>到转动地心坐标系E<sub>g</sub>的转换矩阵;<img file="FDA0000512254140000046.GIF" wi="1049" he="232" />为轨道平面坐标系E<sub>v</sub>到不转动地心坐标系E<sub>o</sub>的转换矩阵;<img file="FDA0000512254140000047.GIF" wi="564" he="258" />为卫星平台坐标系E<sub>r</sub>到轨道平面坐标系E<sub>v</sub>的转换矩阵;<img file="FDA0000512254140000051.GIF" wi="1093" he="285" />为卫星星体坐标系E<sub>e</sub>到卫星平台坐标系E<sub>r</sub>的转换矩阵;<img file="FDA0000512254140000052.GIF" wi="435" he="269" />为天线坐标系E<sub>a</sub>到卫星星体坐标系E<sub>e</sub>的转换矩阵;相反转换通过矩阵求逆实现。 
地址 100191 北京市海淀区学院路37号