空间目标轨道分布特性分析实验报告一、实验目的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的双行轨道根数保存在tee文件中，由两行组成，有效字符为数字0-9，大写字母A-Z、正负号、空格和句点。

其格式如下⑷：比特位1 3 8 10 17 19 32 34 43 45 52 54 61 63 691 NNNNNC NNNNNAAA NNNNN.NNNNNNNN +.NNNNNNNN +NNNNN-N +NNNNN-N N NNNNN2 NNNNN NNN.NNNN NNN.NNNN NNNNNNN NNN.NNNN NNN.NNNNNN.NNNNNNNNNNNNNN 包括的数据项有：行号、卫星编号、卫星密级、卫星国际编号(发射年份、年中序号、发射中产生的目标序列)、历元之年份、历元之天数、平动一阶导数、平动二阶导数、大气阻力项、定规模型的类型、和校验结果、轨道倾角、升交点赤经、偏心率、近地点幅角、平近点角、平均角速度、相对于历元的圈数。

具体定义见表1、2。

表双行根数第二行格式定义（1）1.2为NORAD根据目标初次发现所制定的编目号；（2）1.4〜1.6是遵守国际公约依据最终发射时间的编号，与美国宇航局（NASA的国家空间科学数据中心（NSSDC等一致；（3）1.7〜1.8的历元采用平太阳时。

以下是NORAD公布的第一颗全球星的双行根数，供对比参考。

GLBALSTAR 011 25162U 98008 00069.18987500 -.00000104 00000-0 00000-0 0 34222 25162 52.0028 117.1248 0000635 142.9613 217.1298 12.62272059 95313通过对空间目标轨道信息的分析，确定将从tee文件中提取的数据项有卫星编号、轨道倾角、升交点赤经、偏心率、近地点幅角、平近点角、平均角速度。

四、实验内容1、分析两行根数(TLE)的数据格式；用matlab读取*.TCE文件，本实验主要用的是fopen、importdata和str2num 三个函数, 先用fopen函数打开文件，然后用importdata函数将原文件中的数据以struct型读取出来，并存入S矩阵，其结构如下：F^l S <1x1 struct»Held *Value5data<21150x2 doublestextdata<211^0x7 cell>代码如下：S=fope n( 'stkSatDbAll.tce');S=importdata('stkSatDbAll.tce');采用此种方法读取的数据，可以将原文件中的卫星编号、轨道倾角和轨道偏心率等信息直接以类似矩阵的形式存放，只是数据格式为cell型，暂无法直接进行运算：回S.textdata <21190x7 cell>1234567 1100005U58002B07133.209...*.0000010500000-0-1110S 3220000534.2375173.40671B52O49104.05S4277.1SB43100011U59OO1A07139.835...-.0000004400000-0-26374-4420001132.8792302,54521433772356,1867 2.84165100012U59001B07140.985...-.0000003800000-0-44544-4620001232,904597.98021681770231.0666112.77947100016U58002A0713&765....0000020600000-025323-3820001634.2793264.019720334S9261.727074.9179q100020U59007A07140.613....0000028900000-0105^2-3in7祝2%门礒 1 巧口gn ansi TTNfM Afi利用以下代码提取含有六个轨道根数的偶数行:[m1 n1]=size(S.textdata);[m2 n2]=size(S.data); for i=2:2:m1for j=1: n1A(i/2,j)=S.textdata(i,j);% 提取偶数行endend又观察看出，卫星的角速度数据存放在S.data矩阵的第一列，卫星的轨道倾角数据和偏心率数据分别存放在cell型数据S.textdata的第3列和第5列，将其分别提取出来：for i=2:2:m2for j=1: n2B(i/2,j)=S.data(i,j); % 角速度endendfor i=1:m1/2A1(i)=A(i,3);endA仁A(:,3); %轨道倾角A2=A(:,5); % 偏心率然后再利用str2num函数将数据格式转换为可以计算的矩阵形式：for i=1:le ngth(A1);A_qi ngjiao(i,1)=str2 num(A1{i,1});endfor i=1:le ngth(A2);A_pia nxi nlv(i,1)=str2 num(A2{i,1})*1e-7;end2、对目前在轨目标半长轴、偏心率和轨道倾角的分布规律进行统计分析；(1)在轨目标轨道倾角分布规律统计首先用直方图进行倾角的频数分布统计，分为90个区间，统计图表如下:倾角/度 数量 0 —10 847 10—20 487 20—30 338 30—40 158 40—50 123 50-60 208 60-65 562 65-70 963 70-72 384 72-75822倾角/度 数量 75-80 72 80-82 208 82-83 1024 83-90 375 90-98 497 98-99 1639 99-100 1106 100-110 753 110-15029总计10595(2)在轨目标偏心率分布规律统计图为用直方图进行倾角的频数分布统计，分为20个区间，统计图表如下从上图可以看出轨道倾角的大致分布情况，有一部卫星的倾角在0度至30度之间，而倾角在50度附近的卫星数量较少， 60至80之间也有较多卫星分布， 90至100度倾角之间 卫星的分布密度最大，到了 110度往上，就很少有卫星存在了。

由表格可以看出，卫星轨道偏心率在0——0.1之间分布最多，最为密集，同时也反映出大多数卫星的轨道都是近圆轨道，更详细的分布情况如下所示：(3 )在轨目标半长轴分布规律统计首先用直方图进行倾角的频数分布统计，分为100个区间，统计图表如下7000更详细的半长轴分布情况可以从下图的表格看出:半长轴/千米 数量 0-8000 7833 8000-16000 653 16000-22000 240 22000以上1869 总计105953、分别统计低轨、中轨和高轨目标的分布规律。

从上一实验内容中已经求出卫星轨道半长轴的数值，又知地球的平均半径为6378.145千米，故两者相减即可求出各个卫星的轨道高度。

已知轨道高度小于等于 2000千米的卫星为地轨卫星，轨道高度在2000千米至33000千米之间的为中轨卫星，轨道高度在33000千轨道类型 数量低轨 8041 中轨1608高轨946可以看出，大多数卫星为低轨卫星，有 多个，中轨卫星有 多个，最少的是高轨卫星，只有946个。

60005000 •400030002000 -1000 -2 4 10 12 4x 10五、实验程序clc;clearall;close all;S=fope n( 'stkSatDbAll.tce');S=importdata('stkSatDbAll.tce');[ml n1]=size(S.textdata);[m2 n2]=size(S.data);fori=2:2:m1for j=1: n1A(i/2,j)=S.textdata(i,j);% 提取偶数行endendfori=2:2:m2for j=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)A2*u)A(1 /3); %半长轴，单位千米% A_banchangzhou=(u./(B(:,1)*2*pi/ 24/3600).A2).A1 /3;fori=1:le ngth(A1);A_qi ngjiao(i,1)=str2 nu m(A1{i,1}); %转换为矩阵形式endfori=1:le ngth(A2);A_pia nxi nlv(i,1)=str2 num(A2{i,1})*1e-7;end%%%统计轨道倾角figure;hist(A_qi ngjiao,90);max(A_ qin gjiao)min(A_qin gjiao)e1=fi nd((A_qi ngjiao>=0)&( A_qi ngjiao<10));le ngth(e1);e2=fi nd((A_qi ngjiao>=10)&( A_qi ngjiao<20));le ngth(e2);e3=fi nd((A_qi ngjiao>=20)&( A_qi ngjiao<30));le ngth(e3);e4=fi nd((A_qi ngjiao>=30)&( A_qi ngjiao<40));le ngth(e4);e5=fi nd((A_qi ngjiao>=40)&( A_qi ngjiao<50));le ngth(e5);e6=fi nd((A_qi ngjiao>=50)&( A_qi ngjiao<60));le ngth(e6);e7=fi nd((A_qi ngjiao>=60)&( A_qi ngjiao<65));le ngth(e7);e8=fi nd((A_qi ngjiao>=65)&( A_qi ngjiao<70));le ngth(e8);e9=fi nd((A_qi ngjiao>=70)&( A_qi ngjiao<72));le ngth(e9);e10=fi nd((A_qi ngjiao>=72)&( A_qi ngjiao<75));le ngth(e10);e11=fin d((A_qi ngjiao>=75)&( A_qi ngjiao<80));le ngth(e11);e12=fi nd((A_qi ngjiao>=80)&( A_qi ngjiao<82));le ngth(e12);e13=fi nd((A_qi ngjiao>=82)&( A_qi ngjiao<83));le ngth(e13);e14=fi nd((A_qi ngjiao>=83)&( A_qi ngjiao<90));le ngth(e14);e15=fi nd((A_qi ngjiao>=90)&( A_qi ngjiao<98));le ngth(e15);e16=fi nd((A_qi ngjiao>=98)&( A_qi ngjiao<99));le ngth(e16);e17=fi nd((A_qi ngjiao>=99)&( A_qi ngjiao<100));le ngth(e17);e18=fi nd((A_qi ngjiao>=100)&( A_qi ngjiao<110));le ngth(e18);e19=fi nd((A_qi ngjiao>=110)&( A_qi ngjiao<150));le ngth(e19);%%%%%俭计偏心率figure;hist(A_pia nxi nlv,20);max(A_pia nxin lv)mi n( A_pia nxi nlv)fO=fi nd(A_pia nxinl v==0);le ngth(fO);f1=fi nd((A_pia nxi nl v>=0)&( A_pia nxi nlv<0.001));le ngth(f1);f2=fi nd((A_pia nxi nlv>=0.001)&( A_pia nxi nlv<0.005));le ngth(f2);f3=fi nd((A_pia nxi nlv>=0.005)&( A_pia nxinl v<0.01));le ngth(f3);f4=fi nd((A_pia nxin lv>=0.01)&( A_pia nxi nlv<0.1));le ngth(f4);f5=fi nd((A_pia nxi nl v>=0.1)&( A_pia nxi nlv<0.3));le ngth(f5);f6=fi nd((A_pia nxi nl v>=0.3)&( A_pia nxi nlv<0.5));le ngth(f6);f7=fi nd((A_pia nxi nl v>=0.5)&( A_pia nxi nlv<0.7));le ngth(f7);f8=fi nd((A_pia nxi nlv>=0.7)&( A_pia nxi nlv<1));le ngth(f8);f9=fi nd(A_pia nxinl v==1);le ngth(f9);%%%%统计轨道半长轴figure;hist(A_ba ncha ngzhou,100);max(A_ba ncha ngzhou)min( A_ba ncha ngzhou)g1=fi nd((A_ba nchan gzhou>=0)&( A_ba nchan gzhou<8000));le ngth(g1);g2=fi nd((A_ba ncha ngzhou>=8000)&( A_ba ncha ngzhou<16000));le ngth(g2); g3=fi nd((A_ba nchan gzhou>=16000)&( A_ba ncha ngzhou<22000));le ngth(g3); g4=fi nd(A_ba nchan gzhou>=22000);le ngth(g4);%%%%嗽道高度A_gaodu=A_ba ncha ngzhou-6378.145;h1=fi nd(A_gaodu<2000);digui=le ngth(h1)h2=fi nd((A_gaodu>=2000)&( A_gaodu<=33000));zho nggui=le ngth(h2)h3=fi nd(A_gaodu>33000);gaogui=le ngth(h3)。