MATLAB与STK互联28：仿真案例3—读取轨道六根数（DataProviders使用示例）

这个案例比较简单，生成20个卫星，然后获取第一个轨道历元的卫星轨道六要素数据。
约束：卫星轨道高度1000~6000km，偏心率0-0.2，轨道倾角0-50°，升交点赤经0-360，近地点辐角、真近点角不做要求。
这里重点看一下DataProviders的用法。
下面上代码：
uiap = actxserver(‘STK11.application’);
root = uiap.Personality2;
root.NewScenario(‘elements’);
sc = root.CurrentScenario;

for j = 1:20
satname = [‘str_’,num2str(j)];
sat = sc.Children.New(18,satname);
%设置卫星轨道动力学类型，这里选用高精度HPOP
sat.SetPropagatorType(‘ePropagatorHPOP’);

%设置卫星轨道参数
kep = sat.Propagator.InitialState. Representation.ConvertTo(‘eOrbitStateClassic’);
kep.SizeShapeType = ‘eSizeShapeSemimajorAxis’;
kep.LocationType = ‘eLocationTrueAnomaly’;
kep.Orientation.AscNodeType = ‘eAscNodeLAN’;

kep.SizeShape.SemimajorAxis = 6378.14 + 1000 + 5000 * rand(1);
kep.SizeShape.Eccentricity = 0.2 * rand(1);
kep.Orientation.Inclination = 50 * rand(1);
kep.Orientation.ArgOfPerigee = 360 * rand(1);
kep.Orientation.AscNode.Value = 360 * rand(1);
kep.Location.Value = 360 * rand(1);
sat.Propagator.InitialState.Representation.Assign(kep);
sat.Propagator.Propagate;
end

outputData = cell(21,7);
outputData{1,1} = {‘卫星名称’};
outputData{1,2} = {‘半长轴’};
outputData{1,3} = {‘偏心率’};
outputData{1,4} = {‘倾角’};
outputData{1,5} = {‘升交点赤经’};
outputData{1,6} = {‘近地点辐角’};
outputData{1,7} = {‘真近点角’};

%获取所有卫星的路径
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)
sat = root.GetObjectFromPath(satPaths{i});%注意是大括号
%将时间改为历元形式
root.UnitPreferences.Item(‘DateFormat’).SetCurrentUnit(‘EpSec’);
%计算多个历元语句为times = {0;1000;2000};
times = {0};
elems = {‘Semi-major Axis’;‘Eccentricity’;‘Inclination’;‘RAAN’;‘Arg of Perigee’;‘True Anomaly’};

satDP = sat.DataProviders.Item(‘Classical Elements’).Group.Item(‘ICRF’).ExecSingleElementsArray(times,elems);

outputData{i+1,1} = sat.InstanceName;
outputData{i+1,2} = cell2mat(satDP.GetArray(int32(0)));
outputData{i+1,3} =cell2mat(satDP.GetArray(int32(1)));
outputData{i+1,4} = cell2mat(satDP.GetArray(int32(2)));
outputData{i+1,5} = cell2mat(satDP.GetArray(int32(3)));
outputData{i+1,6} = cell2mat(satDP.GetArray(int32(4)));
outputData{i+1,7} = cell2mat(satDP.GetArray(int32(5)));
end

输出结果如下图所示：

上面代码如果拷贝执行，请注意中文符号的问题。这个案例，主要是介绍DataProviders的用法，可以计算多个历元的多个元素，使用方法如上。有其他的用法，也请大家分享。以上代码，有部分只能在STK10/STK11版本使用，请大家注意。