空间目标轨道分布特性分析实验报告

空间目标轨道分布特性分析实验报告一、实验目的1、了解空间目标轨道分布规律；2、掌握TLE数据格式分析方法；3、掌握空间目标高度分布特性分析方法与过程。二、实验环境Matlab或C语言三、实验原理1、空间目标及其分布空间目标广义是指离地球表面120公里以外空间的所有目标，包括自然天体和人造天体。本研究报告中的空间目标系指环绕在地球周围数万公里内的人造天体，包括卫星、平台和运载，以及上述目标解体后形成的空间碎片。对这些人造目标进行监视属空间目标监视系统的范畴。根据有关研究，环绕地球的空间目标数目大约为35,000,000，其中大小在1～10cm的约110,000个，大于10cm的在8000个以上。目前美国空间目标监视系统可对30cm以上的空间目标进行例行的日常观测，对10cm以上的目标可能观测到，但不能保证例行的日常跟踪。上述空间目标中，到2008年8月24日，被美国空间目标监视系统编目过的空间目标数目为33311个，其中21597个已经陨落，11714个仍在轨。空间目标都有一定大小、形状，运行在一定轨道上，使得每一空间目标都有其独特的轨道特性、几何特性和物理特性。这些特性奠定了对空间目标进行定轨和识别的基础，尤其是在用航天器一般都有特定的外形、稳定的轨道、姿态、温度等特性，是空间目标识别的主要技术支撑。空间目标监视的核心任务是对空间目标进行探测、跟踪和识别。获取空间目标的几何特征、物理特征和运动参数等重要目标信息，进而确定目标威胁度、警戒空间碰撞、提供安全告警信息，是实施防御性空间对抗和进攻性空间对抗的基础。其中在空间目标的识别过程中，空间目标的轨道特性是主要依据，而其几何特性和物理特性则是对其轨道特性的进一步补充。空间目标的分布是指空间分布，空间分布由轨道分布体现，了解空间目标的轨道分布特征对探测系统的设计等都具有非常重要的参考价值。这里主要通过分析轨道根数的分布来完成空间目标分布的统计，主要从空间目标半长轴、偏心率和轨道倾角三个方面进行考虑。2、TLE文件格式分析NASA发布航天器轨道根数的文本文件是TLE格式，TLE文件有多种扩展名，包括*.wr、*.fr、*.gd、*.om、*.sd和*.tce等六种，NORAD的双行轨道根数保存在tce文件中，由两行组成，有效字符为数字0-9，大写字母A-Z、正负号、空格和句点。其格式如下[4]：比特位1381017193234434552546163691NNNNNCNNNNNAAANNNNN.NNNNNNNN+.NNNNNNNN+NNNNN-N+NNNNN-NNNNNNN2NNNNNNNN.NNNNNNN.NNNNNNNNNNNNNN.NNNNNNN.NNNNNN.NNNNNNNNNNNNNN包括的数据项有：行号、卫星编号、卫星密级、卫星国际编号（发射年份、年中序号、发射中产生的目标序列）、历元之年份、历元之天数、平动一阶导数、平动二阶导数、大气阻力项、定规模型的类型、和校验结果、轨道倾角、升交点赤经、偏心率、近地点幅角、平近点角、平均角速度、相对于历元的圈数。具体定义见表1、2。表1双行根数第一行格式定义序号比特位含义示例1.101行号，取值为111.203-07卫星编目号，5位十进制数表示，最多可编目99999个目标234551.308卫星密级分类标识，U表示非密，S表示秘密U1.410-11卫星国际编号之发射年份，2位十进制数表示，如94代表1994年941.512-14卫星国际编号之年中序号，3位十进制数表示，如089表示当年的第89次次发射0891.615-17卫星国际编号之本次发射中产生的目标序列，字符表示，如A表示本次发射中形成的第一个目标A1.719-20历元之年份，2位十进制数表示，如97代表1997年971.821-32历元之天数，表示根数时刻在该年中的天数（DayoftheYear），小数点后保留8位有效数字（精确到1ms）320.909460191.934-43平动一阶导数20.n，单位为圈数/天2.000001401.1045-52平动二阶导数60n，单位为圈数/天3，表示方法为一种调整的指数形式，前6位为0～1的小数部分，后2位为指数部分，如-12345-6为-0.12345*10-600000-01.1154-61大气阻力项之B*0*21ρmSCBD，ρ0为参考点大10191-3气密度)，单位为1/ae（ae为地球半径），表示方法同1.101.1263生成该星历所采用的定轨模型的类型，对外公布的编目中将此值统一设为0，表示由SGP4/SDP4模型生成，NORAD内部为1=SGP，2=SGP4，3=SDP4，4=SGP8，5=SDP801.1365-68一般情况下表示该目标的轨道根数编号(即更新次数)，但在空间监视中心从夏延山主用中心切换到海军备用中心时，该值可能有误2621.1469和校验结果，具体过程是：略去所有空格、加号、字符、小数点，视负号为1，将所有数字相加，除以10后的余数(即和的最后一位)。该校验约可查出90%的错误1表2双行根数第二行格式定义序号比特位含义示例2.101行号，取值为222.203-07卫星编号234552.309-16轨道倾角，单位：度99.00902.418-25轨道升交点赤经，单位：度272.67452.527-33轨道偏心率，小数表示（即1234567为0.1234567）00085462.635-42近地点幅角，单位：度223.16862.744-51平近点角，单位：度136.88162.853-63平均角速度，单位：圈/天14.117117472.964-68相对于历元的圈数（发射后首次过升交点为第一圈），单位：圈148492.1069同1.145其中，需要说明的是：（1）1.2为NORAD根据目标初次发现所制定的编目号；（2）1.4～1.6是遵守国际公约依据最终发射时间的编号，与美国宇航局（NASA）的国家空间科学数据中心（NSSDC）等一致；（3）1.7～1.8的历元采用平太阳时。以下是NORAD公布的第一颗全球星的双行根数，供对比参考。GLBALSTAR01125162U9800800069.18987500-.0000010400000-000000-00342222516252.0028117.12480000635142.9613217.129812.6227205995313通过对空间目标轨道信息的分析，确定将从tce文件中提取的数据项有卫星编号、轨道倾角、升交点赤经、偏心率、近地点幅角、平近点角、平均角速度。四、实验内容1、分析两行根数（TLE）的数据格式；用matlab读取*.TCE文件，本实验主要用的是fopen、importdata和str2num三个函数，先用fopen函数打开文件，然后用importdata函数将原文件中的数据以struct型读取出来，并存入S矩阵，其结构如下：代码如下：S=fopen('stkSatDbAll.tce');S=importdata('stkSatDbAll.tce');采用此种方法读取的数据，可以将原文件中的卫星编号、轨道倾角和轨道偏心率等信息直接以类似矩阵的形式存放，只是数据格式为cell型，暂无法直接进行运算：利用以下代码提取含有六个轨道根数的偶数行：[m1n1]=size(S.textdata);[m2n2]=size(S.data);fori=2:2:m1forj=1:n1A(i/2,j)=S.textdata(i,j);%提取偶数行endend又观察看出，卫星的角速度数据存放在S.data矩阵的第一列，卫星的轨道倾角数据和偏心率数据分别存放在cell型数据S.textdata的第3列和第5列，将其分别提取出来：fori=2:2:m2forj=1:n2B(i/2,j)=S.data(i,j);%角速度endendfori=1:m1/2A1(i)=A(i,3);endA1=A(:,3);%轨道倾角A2=A(:,5);%偏心率然后再利用str2num函数将数据格式转换为可以计算的矩阵形式：fori=1:length(A1);A_qingjiao(i,1)=str2num(A1{i,1});endfori=1:length(A2);A_pianxinlv(i,1)=str2num(A2{i,1})*1e-7;end2、对目前在轨目标半长轴、偏心率和轨道倾角的分布规律进行统计分析；（1）在轨目标轨道倾角分布规律统计首先用直方图进行倾角的频数分布统计，分为90个区间，统计图表如下:05010015005001000150020002500从上图可以看出轨道倾角的大致分布情况，有一部卫星的倾角在0度至30度之间，而倾角在50度附近的卫星数量较少，60至80之间也有较多卫星分布，90至100度倾角之间卫星的分布密度最大，到了110度往上，就很少有卫星存在了。编写代码更细一步统计卫星分布情况，列出表格如下：倾角/度数量倾角/度数量0—1084775-807210—2048780-8220820—3033882-83102430—4015883-9037540—5012390-9849750-6020898-99163960-6556299-100110665-70963100-11075370-72384110-1502972-75822总计10595（2）在轨目标偏心率分布规律统计下图为用直方图进行倾角的频数分布统计，分为20个区间，统计图表如下:00.10.20.30.40.50.60.70.80.90100020003000400050006000700080009000由表格可以看出，卫星轨道偏心率在0——0.1之间分布最多，最为密集，同时也反映出大多数卫星的轨道都是近圆轨道，更详细的分布情况如下所示：偏心率/度数量等于000-0.00113610.001-0.00530330.005-0.0116710.01-0.131110.1-0.32100.3-0.51210.5-0.77050.7-1383等于10总计10595（3）在轨目标半长轴分布规律统计首先用直方图进行倾角的频数分布统计，分为100个区间，统计图表如下:024681012x10401000200030004000500060007000更详细的半长轴分布情况可以从下图的表格看出：半长轴/千米数量0-800078338000-1600065316000-2200024022000以上1869总计105953、分别统计低轨、中轨和高轨目标的分布规律。从上一实验内容中已经求出卫星轨道半长轴的数值，又知地球的平均半径为6378.145千米，故两者相减即可求出各个卫星的轨道高度。已知轨道高度小于等于2000千米的卫星为地轨卫星，轨道高度在2000千米至33000千米之间的为中轨卫星，轨道高度在33000千米以上的为高轨卫星。统计出低中高轨卫星数量分布情况如下：轨道类型数量低轨8041中轨1608高轨946可以看出，大多数卫星为低轨卫星，有8000多个，中轨卫星有1600多个，最少的是高轨卫星，只有946个。五、实验程序clc;clearall;closeall;S=fopen('stkSatDbAll.tce');S=importdata('stkSatDbAll.tce');[m1n1]=size(S.textdata);[m2n2]=size(S.data);fori=2:2:m1forj=1:n1A(i/2,j)=S.textdata(i,j);%提取偶数行endendfori=2:2:m2forj=1:n2B(i/2,j)=S.data(i,j);%角速度endendfori=1:m1/2A1(i)=A(i,3);endA1=A(:,3);%轨道倾角A2=A(:,5);%偏心率u=398600.47;A_banchangzhou=((24*3600./B(:,1)/2/pi).^2*u).^(1/3);%半长轴，单位千米%A_banchangzhou=(u./(B