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