MATLAB与STK联合仿真：轨道数据导出Excel全流程指南

简介：本资源是一套面向航天工程学习者与遥感数据分析初学者的MATLAB-STK协同开发实践方案，聚焦卫星轨道数据自动化提取、处理与Excel标准化导出这一典型跨平台需求。资源共38个文件，主体为30个MATLAB函数（.m），涵盖STK连接控制、卫星六要素与RV矢量解析、轨道数据滤波平滑、多维矩阵转表等核心功能；另含6个备份文件（.zbak）及1个说明文档（.txt），整体压缩包仅36KB，轻量易部署。已有53人学习下载，适合具备基础MATLAB编程能力、正开展轨道力学仿真或任务分析的学生与科研人员。读者可直接复用STKConnect、GetSatelliteData、ConvertToExcel等模块化函数，快速构建从STK场景读取→坐标/速度/时间参数解析→单位归一化→滑动窗口滤波→Excel多Sheet结构化写入的完整流水线，并通过main_script实现日志记录与参数化配置，显著提升轨道数据后处理效率与可重复性。 我经常在项目群看到类似的问题——“STK仿真结果怎么拿给同事看？同事电脑没装STK；轨道数据怎么弄到Excel里？”。其实这个问题的本质，是三个工具各自擅长的事情不一样：STK负责精确建模仿真轨道，MATLAB负责数值分析和算法设计，Excel负责最终交付和汇报。三者之间的数据通道一旦打通，很多看似繁琐的卫星轨道数据处理工作，就能沉淀成一套可复用的自动化流程。

这篇内容我会从工具分工讲起，完整走一遍“MATLAB控制STK生成轨道数据，再把数据导出为Excel”的链路，包括环境准备、接口连接、数据提取、坐标转换、Excel导出这些关键环节，最后附上我自己踩过的几个高频坑。适合刚接触STK与MATLAB联合仿真的学生，也适合需要在项目中交付轨道数据的工程师参考。

1. 为什么 MATLAB 和 STK 需要连起来用（不是选一个就行的）

1.1 各自的强项，正好互补

STK（Systems Tool Kit）在航天领域几乎是轨道仿真和任务分析的标配。它的传播模型很成熟，J2、JGM-3、HPOP这些高阶模型都是经过长年校验的，可见性分析、敏感器覆盖、地面站访问计算也直接内置。你只要在界面里拖几个对象，设置好初始轨道和传播器，就能得到一份具有真实物理含义的轨道预报。

但STK的强项同时也是它的短板：它是一个仿真软件，不是一个数值计算平台。当你需要做批量处理、设计定制算法、做参数扫描，或者跟自己的控制模型耦合时，纯靠STK的界面操作会非常笨重，纯脚本也不是大多数人的主场。

MATLAB则正好相反。矩阵运算、循环批处理、算法验证、绘图可视化都极其顺手。而且绝大多数做轨道动力学、制导导航与控制的人，从学生阶段起就在用MATLAB，几乎零门槛。用MATLAB做轨道外推、姿态控制、载荷数据分析，或者做多圈轨道数据的趋势拟合，都远比比在STK内部点界面高效得多。

1.2 为什么 Excel 也卷了进来

至于Excel，原因很实际：项目评审、阶段报告、数据分析、给不装STK的同事看数据，最通用的交付格式就是表格。一个 .xlsx 文件谁都能打开，不需要license，也不需要培训。尤其是当你做轨道精度比较、覆盖统计、载荷数据分析这类需要给人复核的工作时，Excel的筛选、透视表、图表还能帮上忙。

所以这条链路的典型工作流就是：

• STK根据任务需求构建场景、创建卫星、设置传播模型、生成轨道和访问数据；
• MATLAB通过接口读取这些数据，做自定义计算、绘图、统计、优化；
• MATLAB将最终结果整理成Excel报告，交给下游或者归档。

1.3 打通数据链路的三种常见方式

在讲具体实现之前，先把方案格局说清楚，免得后面看代码一头雾水。

第一种是STK自带的MATLAB集成包，本质上是把STK的Connect命令封装成了 stkInit、stkOpen、stkConnect 等函数。这是最主流、也最稳的方式，跨版本兼容性相对好，支持批处理。

第二种是COM接口，在Windows上用 actxserver 直接创建STK的COM对象，可以控制STK的底层对象模型。灵活度高，但要求你对STK的对象层次比较熟，而且只能在Windows上用。

第三种是文件交换，比如STK导出 .e、.sa、.v 等格式文件，MATLAB再读取。这个方式简单直接，但不适合自动化，适合一次性交换。

本文的核心是第一种，因为它在自动化程度、可维护性、可读性上最平衡，也是我做项目时用得最多的方案。

2. 环境准备：版本匹配、安装路径和接口函数，这三件事最耽误时间

2.1 版本匹配是第一个坑

STK从11之后版本号变动不算大，但MATLAB接口的兼容性是真有讲究。老的STK 10/11配的MATLAB接口，在较新版本的MATLAB上可能加载不成功。我实际遇到过在MATLAB R2023a上调用 stkInit 直接报 “Undefined function 'stkInit'”，后来发现是接口路径没加对，并不是函数本身不存在。

如果你的STK安装目录下有类似 Matlab 或 MATLAB Integration 的文件夹，里面通常包含 stkInit、stkConnect 等函数的源码和编译好的mex文件。你需要把它加入MATLAB搜索路径：

addpath('C:\Program Files\AGI\STK 12\MATLAB');

不同的STK版本路径不同，有的是C:\Program Files\AGI\STK 11\Matlab。最稳妥的方法是打开STK安装目录，直接搜索stkInit.m这个文件在哪，就把它的父目录加进去。

需要提醒一点：如果你的STK是64位，那MATLAB也必须是64位，否则mex文件加载不了。这个错误往往不是“mex文件无法加载”这种直白提示，而是各种莫名其妙的崩溃或者连接失败。查这个问题的时候，先看一眼两侧的位数，能省很多时间。

2.2 启动 STK 之前先确认端口

MATLAB与STK的通信基于TCP/IP，默认的Connect端口是5001。如果你的机器上有安全软件把端口占了，或者以前启动过多个STK实例，都可能导致连接失败。

有个小技巧，在连接之前，可以在MATLAB里执行：

system('netstat -ano | findstr :5001');

看看5001端口有没有被占用。如果被占用，多半是有个残留的STK进程没退出，去任务管理器里把STK.exe结束掉再重来。这个排查步骤虽然简单，但能解决绝大多数“连接半天没反应”的问题。

2.3 接口工具箱是否可用

在较新版本的MATLAB中，如果你的STK版本较老，可能没有官方接口函数。这时可以退而求其次，直接用tcpclient写一个简单的Connect TCP客户端。但说实话，除非你有特殊需求，否则我不建议这么做，因为STK Connect的响应协议虽然不难，但调试起来烦人，远不如直接用官方封装好的 stkConnect 函数。

如果你已经加好了接口路径，可以用which确认一下：

which stkConnect

如果返回路径正确，说明环境就绪。如果返回空白，说明路径没加对，回头再检查安装目录。

另外一个实用技巧：STK自带的“Command Browser”可以在界面里查看和测试Connect命令。你在里面直接搜SetState、GetState，能看到当前版本的准确参数格式，还能在GUI里试跑一遍。这个工具对语法核对特别有用，我后面写的代码如果和你手里的STK版本提示参数不对，请一定打开Command Browser复制对应格式。

3. 用 Connect 命令搭起 MATLAB 与 STK 的通信桥（附可直接跑的代码）

3.1 建立连接的最小可执行代码

先看最核心的三行：

stkInit; % 初始化 STK-MATLAB 接口 remMachine = stkDefaultHost; % 获取默认主机地址 conid = stkOpen(remMachine); % 打开连接，返回连接句柄

执行之后，STK如果没有运行，stkOpen通常会自己启动一个STK实例。注意STK界面会弹出来，这很正常。如果不想看到界面，可以通过STK Engine的模式运行，但桌面版通常都是可视化运行，这对调试反而更友好。

连接建立之后，推荐在脚本开头就把场景清理干净，避免上一次运行残留的对象干扰：

stkClose(conid); % 关闭连接 stkInit; % 重置接口 conid = stkOpen(stkDefaultHost); stkConnect(conid, 'New', 'Scenario/TempScenario');

如果你不想要残留场景，可以直接新建一个场景，STK会覆盖掉之前的场景。

3.2 创建场景、卫星并设置轨道根数

场景创建好了，下一步是添加卫星对象并设置轨道。下面这段代码创建了一个500km高度、倾角53°的低轨卫星：

scenPath = 'Scenario/TempScenario'; stkConnect(conid, 'New', 'Scenario/TempScenario'); % 设置场景时间 stkConnect(conid, 'SetTimePeriod', scenPath, '1 Jun 2024 00:00:00.000', ... '2 Jun 2024 00:00:00.000'); stkConnect(conid, 'Animate', scenPath, 'Time', '1 Jun 2024 00:00:00.000');

创建卫星：

stkConnect(conid, 'New', 'Satellite/LEO_Sat'); % 设置经典轨道根数 stkConnect(conid, 'SetState', 'Satellite/LEO_Sat', 'Classical', 'J2Perturbation', 'J2000', ... '1 Jun 2024 00:00:00.000', ... % 历元 '1 Jun 2024 00:00:00.000', ... % 轨道开始时间 '94.4', ... % 轨道周期（分钟） '6878.0', ... % 半长轴（km） '0.001', ... % 离心率 '53.0', ... % 倾角（deg） '45.0', ... % 升交点赤经（deg） '20.0', ... % 近地点幅角（deg） '0.0'); % 真近点角（deg）

这里有个非常关键的点：SetState的第三个参数是传播模型，我写的是J2Perturbation，表示考虑J2摄动的简单模型。如果你需要更精确的仿真，可以用HPOP，但HPOP对初始条件的要求更多，这里先用J2演示。

关于轨道周期：实际周期需要根据半长轴计算。按 T = 2π√(a³/μ) 计算，a=6878 km，μ=398600.4418 km³/s²，T≈94.4分钟。你可以让STK自己算，也可以在SetState中给一个近似值，它通常会根据其他轨道参数重新计算。

不同STK版本中SetState的参数顺序可能略有差异。我按STK 12.x的常见写法来写，如果你的版本提示参数个数不对，打开STK的Command Browser，搜索SetState，复制当前版本的格式即可。

设置完成之后，一定要执行传播：

stkConnect(conid, 'Propagate', 'Satellite/LEO_Sat', 'Interval', ... '1 Jun 2024 00:00:00.000', '2 Jun 2024 00:00:00.000');

不传播是新手最容易犯的错误——你只设了初始状态，没有做时间积分，STK里卫星就只有初始时刻的一个点，后面要报数据全是空。

3.3 关闭连接的正确姿势

运行结束，记得关闭连接。直接stkClose(conid)即可。如果脚本报错中断了，也可以用stkClose来清理。

提示：STK连接不像文件句柄那样会自动回收，脚本崩了之后残留的连接会占用端口，下一次连接就可能失败。所以建议在try ... catch ... end结构里确保 finally 块执行stkClose(conid)。

4. 轨道数据提取的细节：坐标系、时间基准和 GetState 返回值的解析

4.1 先理解坐标系，否则数据错得你不知道错在哪

轨道数据提取最容易出错的不是命令写错，而是坐标系选错。STK中常见的坐标系有：

• J2000（ECI惯性系）：以J2000历元的平赤道和春分点定义，适合轨道力学计算。
• 地固系（Fixed / ECF）：随地球自转，适合表达地面轨迹、星下点。
• TEME：TLE使用的近似惯性系。
• LLA（经纬高）：大地坐标系，基于WGS84椭球。

如果你要把轨道数据用来做星下点覆盖分析，那必须用地固系或LLA；如果你要做轨道力学、控制或姿态仿真，通常用J2000。两个坐标系之间涉及地球自转角GMST的转换，这不是简单矩阵加减，所以在读取数据之前就要想清楚目的。

4.2 用 GetState 提取位置速度

最常用的取数命令是GetState。下面这行代码返回某时刻卫星在J2000坐标系下的位置和速度：

state = stkConnect(conid, 'GetState', 'Satellite/LEO_Sat', 'Cartesian', 'J2000', ... '1 Jun 2024 01:00:00.000');

返回值是一个字符串，类似：

"6678.123 -1523.456 2356.789 -3.456 5.678 1.234"

这六个数字依次是X、Y、Z（km）和Vx、Vy、Vz（km/s）。注意：单位不是米，是千米。这个细节容易漏，导出Excel后如果自己又做了单位换算，很容易差1000倍。

在MATLAB里解析：

values = textscan(state, '%f', 'Delimiter', ' '); posVel = values{1}'; % 1x6 向量 pos = posVel(1:3); % 位置 km vel = posVel(4:6); % 速度 km/s

4.3 批量提取多个时间点的数据

手动一条条GetState显然不是工程做法。更常见的做法是循环遍历时间数组：

startTime = datetime('2024-06-01 00:00:00', 'Format', 'yyyy-MM-dd HH:mm:ss'); stopTime = datetime('2024-06-02 00:00:00', 'Format', 'yyyy-MM-dd HH:mm:ss'); step = minutes(1); tVec = startTime:step:stopTime; n = length(tVec); states = zeros(n, 6); timeStrings = strings(n, 1); for i = 1:n tStr = char(tVec(i)); % 转换为 STK 支持的格式：01 Jun 2024 00:00:00.000 stkTimeStr = sprintf('%02d %s %04d %02d:%02d:%02d.000', ... day(tVec(i)), month(tVec(i), 'name'), year(tVec(i)), ... hour(tVec(i)), minute(tVec(i)), second(tVec(i))); timeStrings(i) = stkTimeStr; state = stkConnect(conid, 'GetState', 'Satellite/LEO_Sat', ... 'Cartesian', 'J2000', stkTimeStr); vals = textscan(state, '%f', 'Delimiter', ' '); states(i, :) = vals{1}'; end

stkTimeStr的格式要特别注意，STK默认识别1 Jun 2024 00:00:00.000这种带月份缩写和时间毫秒的字符串，纯数字格式不一定能识别，虽然某些版本支持YYYY-MM-DD，但为了兼容性，建议用这种标准UTCG格式。

4.4 用 Report

本文还有配套的精品资源，点击获取