这个案例比较简单,生成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版本使用,请大家注意。
|
请发表评论