資源簡介
該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
- 上一篇:離散數學.pdf
- 下一篇:iexm
ltls.exe
評論
共有 條評論