首页 > 专利查询

matlab绘制卫星轨道,[转载]卫星的

matlab绘制卫星轨道,[转载]卫星的两⾏星历及转化为轨道六

根数的MATLAB程序

function [oe,epoch,yr,M,E,satname] = TLE2oe(fname);褐煤干燥

% fname is a filename string for a file containing

% a two-line element set (TLE)

% oe is a 1/6 matrix containing the orbital elements

% [a e i Om om nu]

% yr is the two-digit year

% M is the mean anomaly at epoch

% E is the eccentric anomaly at epoch

悠悠球轴承

% satname is the satellite name

% Calls Newton iteration function file EofMe.m

function

v[oe,epoch,yr,M,E,satname] = );

% Open the file

up and scan in the elements

fid =

fopen(fname, 'r');

A = fscanf(fid,'%13c%*s',1);

B =

fscanf(fid,'%d%6d%*c%5d%*3c%2d%f%f%5d%*c%*d%5d%*c%*d%d%5d',[1,10]);

C = fscanf(fid,'%d%6d%f%f%f%f%f%f',[1,8]);

fclose(fid);

satname=A;

% The value of

mu is for the earth

mu = 3.986004415e5;

% Calculate

2-digit year (Oh no!, look out for Y2K bug!)

yr =

B(1,4);

低通滤波器应用实例% Calculate

epoch in julian days车门铰链

epoch = B(1,5);

%ndot = B(1,6);

% n2dot = B(1,7);网络广告监测系统

% Assign

variables to the orbital elements

i =

C(1,3)*pi/180; % inclination

Om =

C(1,4)*pi/180; % Right Ascension of the Ascending Node e =

C(1,5)/1e7; % Eccentricity

om =

C(1,6)*pi/180; % Argument of periapsis

M =

C(1,7)*pi/180; % Mean anomaly

n = C(1,8)*2*pi/(24*3600); % Mean

motion

% Calculate the

semi-major axis

a = (mu/n^2)^(1/3);

% Calculate the

eccentric anomaly using mean anomaly

E = EofMe(M,e,1e-10);

% Calculate

true anomaly from eccentric anomaly

cosnu = (e-cos(E)) / (e*cos(E)-1);

sinnu = ((a*sqrt(1-e*e)) / (a*(1-e*cos(E))))*sin(E);

nu = atan2(sinnu,cosnu);

if (nu<0), nu=nu+2*pi; end

% Return the

orbital elements in a 1x6 matrix

oe = [a e i Om om nu];

本文发布于:2024-09-23 11:19:12，感谢您对本站的认可！

本文链接：https://www.17tex.com/tex/4/175839.html

上一篇：中国大地测量的发展趋势与研究

下一篇：400公里小时与350公里小时高速铁路主要技术参数对比分析

标签：轨道应用网络广告程序铰链实例系统

留言与评论（共有 0 条评论）