matlab调用stk生成walker星座，并批量导出tle

gj-boy

已于 2024-02-05 14:24:06 修改

阅读量1.1k

点赞数 8

文章标签： matlab 开发语言

于 2024-02-05 09:38:58 首次发布

本文链接：https://blog.csdn.net/ghboyone/article/details/136036327

版权

uiap = actxserver('STK11.application');
root = uiap.Personality2;
root.NewScenario('TLE');
sc = root.CurrentScenario;

%建立种子卫星，并设置轨道参数。500km、倾角50°的圆轨道，其他参数选择默认值，并添加一个Sensor（Sensor形状为简单圆锥型，半张角为77°）
sat = sc.Children.New(18,'mysat');
kep = sat.Propagator.InitialState.Representation.ConvertTo('eOrbitStateClassical');
kep.SizeShapeType = 'eSizeShapeAltitude';
kep.SizeShape.ApogeeAltitude = 500;
kep.SizeShape.PerigeeAltitude = 500;
kep.Orientation.Inclination = 50;
sat.Propagator.InitialState.Representation.Assign(kep);
sat.Propagator.Propagate;

%添加Sensor，并进行设置
% sen = sat.Children.New('eSensor','mysen');
% sen.CommonTasks.SetPatternSimpleConic(77,1);%这里加一条跟星座创建无关的设置，半张角选用77°
% sen.VO.ProjectionType = 'eProjectionEarthIntersections';%保留了Sensor与地球的截线。

%生成Walker星座了。这里是利用root.ExecuteCommand来是实现。具体语句如下：
% 生成一个Delta类型、4轨道面、每个轨道面3颗星、相位因子2、同轨道面颜色相同的星座。
root.ExecuteCommand('Walker */Satellite/mysat Type Delta NumPlanes 4 NumSatsPerPlane 3 InterPlanePhaseIncrement 2 ColorByPlane Yes');

%两行轨道根数，需要卫星编号，这里假设初始编号60000
ssc = 60000;
%创建星历存储文件
fid = fopen('TLEData.tle','wt');

%获取所有卫星的路径
satpathcollection = root.ExecuteCommand('ShowNames * Class Satellite');
satpathcollection.Item(0);
%Item(0)中包含了所有卫星路径，但是是用一个字符串包含所有路径，
%中间用空格隔开，下面的语句是先将字符串分割，然后去除空格元素，
%则satPaths的每个元素都只包含一个卫星路径
satPaths = regexp(satpathcollection.Item(0),' ','split');
satPaths(cellfun(@isempty,satPaths)) = [];

%生成两行轨道根数星历，并输出至文件中
for i = 1:length(satPaths) % 从1开始排除种子卫星
    sattemp = root.GetObjectFromPath(satPaths{i});
    start = sc.StartTime;
    cmd1 = ['GenerateTLE ',satPaths{i},' Point "',start,'" ', sprintf('%05.0f',ssc) , ' 20 0.01 SGP4 ', ' ',sattemp.InstanceName];
    root.ExecuteCommand(cmd1);
    satDP = sattemp.DataProviders.Item('TLE Summary Data').Exec();
    TLEData = satDP.DataSets.GetDataSetByName('TLE').GetValues;

    fprintf(fid, '%s\n%s\n', TLEData{1,1}, TLEData{2,1});
    ssc = ssc +1 ;
end
fclose(fid);

参考了两篇csdn的文章，具体链接为：

MATLAB与STK互联10：卫星对象操作（2）—卫星轨道参数设置（方法1，通过轨道生成器设置实现）_kep.sizeshapetype = 'esizeshapealtitude'; kep.loca-CSDN博客 MATLAB与STK互联42：创建星座（Walker星座）-CSDN博客