本文主要是介绍使用卫星仿真软件STK的一些应用和思考(星地链路、星间链路),希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!
目录
- 任务描述
- 利用STK建模星地协同系统
- 3个GEO高轨卫星
- 240/20/1 Walker-Star Constellation 低轨卫星星座
- 地面站或者地面设备
- 链路建模与数据提取处理
- 星地链路
- 星间链路数据读取的几种方法
- 最麻烦的方法
- 使用Matlab与STK互联接口使用
- 大规模使用Chain
- 总结
任务描述
在一个星地协同的空天地网络中,科研中可能需要建模星地链路以及星间链路,主要受到距离、仰角等参数的影响。
利用STK建模星地协同系统
此处略过Scenario的创建什么的
3个GEO高轨卫星
插入地球同步轨道卫星
240/20/1 Walker-Star Constellation 低轨卫星星座
先插入种子卫星,再根据种子卫星利用Walker Tool生成想要的星座,注意要选择create constellation选项以创建星座对象,后面可能会用到。
地面站或者地面设备
按坐标或者facility插入地面实体即可
链路建模与数据提取处理
星地链路
地面设备将本地信息上传至距离最近的卫星:需得出的结果是地面设备上传本地信息给了哪一颗卫星。
方法是将刚刚创建的星座对象与某一地面设备对象建立chain:注意先后顺序是地面IoT设备在前为发射端,且IoT设备需要加入Transmitter对象。
值得注意的是,STK中可以详细设置Transmitter和Antennna的参数,以获取通信链路的相关参数,这部分暂时省略。
通过Chain->Compute Accesses得到如下可视化结果:
再通过Report & Graph Manager获取所需的数据,可定制想要的数据然后导出为csv/txt文件(如距离、链路时间等)。
星间链路数据读取的几种方法
最麻烦的方法
在每颗卫星对象上插入transmitter/receiver,利用compute access计算所有连接的相关数据,将数据保存后进行后处理用于所需要的场景,这种方法复杂度简直太高了,果断抛弃。
使用Matlab与STK互联接口使用
STK11.6仅仅兼容Matlab2018b以下的版本,但实测2022b、2024a版本都能运行,以下是一个简单demo。
% 创建与STK的连接
clear,clc;
stkApp = actxserver('STK11.Application');% 打开STK场景
stkRoot = stkApp.Personality2;
stkRoot.LoadScenario('D:\STK\Walker-Delta\Walker-Delta.sc');
sc = stkRoot.CurrentScenario;% 获取场景内所有卫星
allChildren = sc.Children;
allSatellites = allChildren.GetElements('eSatellite');
satNum = allSatellites.Count;
% 初始化空数组来存储GEO和LEO卫星
% geoSatellites = {};
% leoSatellites = {};GEO1 = allSatellites.Item('GEO1');
GEO2 = allSatellites.Item('GEO2');
GEO3 = allSatellites.Item('GEO3');% for i = 3:1:satNum-1
% eval(['LEOS', num2str(i), ' = allSatellites.Item(cast(i, ''int32''));'])
% end
count = 1;
for i = 3:1:satNum-1if mod(i-2, 20) == 1 && i ~= 1count = count + 1;endindex = mod(i-3, 20) + 1;eval(['LEOS', num2str((count-1) * 100 + index), '= allSatellites.Item(cast(i, ''int32''));'])
% satellites.(key) = ;
end% 使用 Data Provider 接口获取位置数据
dp101 = LEOS101.DataProviders.Item('Cartesian Position').Group.Item('ICRF').Exec(sc.StartTime, sc.StopTime, 10);
pos101 = dp101.DataSets.ToArray();dp102 = LEOS102.DataProviders.Item('Cartesian Position').Group.Item('ICRF').Exec(sc.StartTime, sc.StopTime, 10);
pos102 = dp102.DataSets.ToArray();% 提取位置坐标
for i = 1:1:length(pos101)x1(i,1) = pos101{i,2};y1(i,1) = pos101{i,3}; z1(i,1) = pos101{i,4};x2(i,1) = pos102{i,2}; y2(i,1) = pos102{i,3}; z2(i,1) = pos102{i,4};
end% 计算距离
distance = sqrt((x2 - x1).^2 + (y2 - y1).^2 + (z2 - z1).^2);% fprintf('The distance between LEOS101 and LEOS102 is %.2f km\n', distance / 1000); % 将距离转换为公里% 断开与STK的连接
stkApp.Quit;
delete(stkApp);
这里读取了一个80/4/1星座命名为LEOS101-LEOS120, LEOS201-LEOS220, LEOS301-LEOS320, LEOS401-LEOS420同轨道相邻的两颗卫星在场景开始和结束时间之内,每隔10s的距离数据,如下所示:
大规模使用Chain
首先添加Constellation对象,将1个轨道的卫星放进去,即将每个轨道建成一个Constellation对象。
第二步,建立Chain,包含的对象是Orbit的Constellation对象和整个星座的Constellation对象Walker-Star。
通过Report & Graph Manager定制所需的数据即可。此时获取了每颗卫星的所有链路信息,但实际上最多4-5条星间链路是能够建立的,我在python中对保存数据进行读取时进行了一些筛选。
总结
本贴只是给出了一些获取数据的方法,但是后处理的代码暂时还未给出,下次再更新吧。
这篇关于使用卫星仿真软件STK的一些应用和思考(星地链路、星间链路)的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!