读取TLE卫星两行根数

该mathlab程序运用于读取.txt格式的卫星两行根数文件，并统计卫星信息

clc;
clear;

fidin = fopen（‘resource.txt‘‘r‘）;
% fidout=fopen（‘tongji.xls‘‘w‘）;    %创建一个表格文件，用于后面的存储数据
% data = textscan（fid ‘%s‘  ‘Delimiter‘ ‘\n‘）;
% data1 = data{:};  
% n_line = size（data1 1）;  %行数
data=‘‘;
N=0;
miu=3.986e5; 
while ~feof（fidin）        %判断是否文件读完
    tline=fgetl（fidin）;  %逐行读入数据
    data=[data tline];  
     N=N+1;  %统计文件的行数           
end
satlitNum=N/3;      %卫星数目
fclose（fidin）;
%+++++++++++++++++提取数据+++++++++++++++++
for i=1:satlitNum 
    bianhao=str2NUM（data（-136+162*i:-131+162*i））;
    X（1i）=bianhao;        %将卫星的编号存到矩阵第一行中

    qingjiao=str2NUM（data（-60+162*i:-53+162*i））;% 09-16轨道倾角，单位:度
    X（2i）=qingjiao;      %将倾角值保存到矩阵第二行中
    
    Omiga=str2NUM（data（-52+162*i:-44+162*i））;
    X（3i）=Omiga;         %将升交点赤经存到矩阵第三行
    
    e=str2double（data（-42+162*i:-36+162*i））;
    e=e*10^（-7）;
    X（4i）=e;         %将偏心率保存到矩阵第三行中
    
    omiga=str2NUM（data（-34+162*i:-27+162*i））;
    X（5i）=omiga;          %将近地点角幅存入矩阵的第五行
    
    M=str2NUM（data（-26+162*i:-18+162*i））;
    X（6i）=M;          %将近地点角幅存入矩阵的第六行
    
    n=str2double（data（-16+162*i:-6+162*i））;   %取出平均角速度
    a=（86400^2*miu/（4*pi^2*n^2））^（1/3）
    X（7i）=a;
end
%+++++++++++++++++画散点图++++++++++++++++
figure（1）;
x=1:satlitNum;
y=X（2x）;
scatter（xy‘.‘‘K‘）;
title（‘卫星的轨道倾角分布情况散点图‘）;  
xlabel（‘卫星排序（第x颗卫星）‘）; ylabel（‘轨道倾角（度）‘）;
figure（2）;
x=1:satlitNum;
y=X（4x）;
scatter（xy‘.‘‘K‘）;
title（‘卫星的轨道偏心率分布情况散点图‘）;  
xlabel（‘卫星排序（第x颗卫星）‘）; ylabel（‘偏心率‘）;
figure（3）;
x=1:satlitNum;
y=X（7x）;
scatter（xy‘.‘‘K‘）;
title（‘卫星的长半轴分布情况散点图‘）;  
xlabel（‘卫星排序（第x颗卫星）‘）; ylabel（‘长半轴（千米）‘）;

%将相关的数据保存起来。
save shujutongji.mat X ;
save shujutongji.txt X -ASCII;
% xlswrite（‘tongji.xls‘X‘sheet1‘）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       2001  2014-01-07 17:23  TLEread.m

     文件      16968  2013-11-05 20:53  resource.txt

----------- ---------  ---------- -----  ----

                18969                    2