MATLAB Aerospace Toolbox飞行仿真与卫星轨道实践

简介：这份 MATLAB Aerospace Toolbox R2023b 官方用户手册面向航空航天工程、飞控与导航仿真方向的 MATLAB 使用者，包括高校学生、科研人员和工程技术人员，用于解决工具箱函数查询、模块调用与官方示例对照的问题。手册由 MathWorks 发布，系统介绍该工具箱在航空航天数据分析、飞行与导航模拟、系统优化等场景中的函数与工具用法，并附函数语法说明、典型示例和历次版本修订记录，适合离线检索，可作为常备工具书。资源包为单一 PDF 文件，共 1 个文件，约 15.14MB，内容完整、目录层级清晰，便于按主题定位章节。目前已有 687 人学习，说明其在同类资料中具备一定参考价值。读者可据此快速熟悉工具箱的接口规范与开发流程，减少查阅在线文档的时间，并为飞行数据处理、导航系统设计等实践提供可直接借鉴的示例思路。

1. 手里有一堆 .dcm、CSV 和 TLE，却不知道从哪一页翻起

前阵子帮人收拾一个飞行数据处理的小摊子：桌面上躺着 DATCOM 导出的.dcm、一段带时间戳的六自由度轨迹 CSV、几颗卫星的 TLE，还有一份导师留下的旧脚本。打开 MATLAB R2023b 的 Aerospace Toolbox 用户手册，目录从坐标系一路排到卫星场景和固定翼分析，看着全，真到用的时候又不知道先看哪节。这份 R2023b 手册的价值就在这里——它不是教程，是把「坐标转换、单位、气动系数、动画、轨道」这些环节按主题摆好的参考书。适合两类人：已经会用 MATLAB、但没系统用过这个工具箱的工程师；以及课程设计或项目里要快速出可视化结果的开发者。下面按我实际拆的顺序走一遍，从坐标系到 DATCOM，再到动画、卫星场景，最后把踩过的坑摊开。

2. 坐标系与单位：Body、Wind、NED、ECI、ECEF 的转换链怎么搭

2.1 五套坐标系各自管什么

手册第 2 章开头就把坐标系分了三组：建模用的机体系（Body）和风轴系（Wind），导航用的 NED、ECI、ECEF，以及显示用的那套。很多人第一次翻会跳过这节，直接去找函数，结果在后面拼旋转矩阵时来回翻车。我的建议是先把下面这张表对着手册 2-4 到 2-8 页过一遍，再动手写代码。

这里最容易混的是大地纬度（geodetic）和地心纬度（geocentric）。手册在 2-6 页专门用一小节讲两者的差别，做远距离导航或卫星地面站计算时，用错纬度会让地面轨迹偏出几公里。另一个常被忽略的是「被动变换」的概念：同一个矢量在不同坐标系下描述，旋转矩阵的方向取决于你是转矢量还是转坐标系，angle2dcm和dcm2angle的方向约定要跟手册对齐。

2.2 用 angle2dcm 搭姿态角到机体系的转换链

姿态角到方向余弦矩阵是整条链的起点。手册里给的顺序约定是 Z-Y-X，也就是偏航-俯仰-滚转，对应angle2dcm的第四个参数'ZYX'。这段代码我一般会先跑一遍，确认矩阵正交再往下接。

% 姿态角（弧度）转机体系到 NED 的方向余弦矩阵 yaw = deg2rad(30); % 偏航角 pitch = deg2rad(5); % 俯仰角 roll = deg2rad(-2); % 滚转角 dcm_b2n = angle2dcm(yaw, pitch, roll, 'ZYX'); % 机体系到风轴系，输入迎角与侧滑角 alpha = deg2rad(4); beta = deg2rad(1); dcm_b2w = dcmbody2wind(alpha, beta); % 机体系下的速度矢量转到风轴系 v_b = [230; 0; 0]; % 机头方向 230 m/s v_w = dcm_b2w * v_b; % 正交性检查：旋转矩阵乘自身转置应接近单位阵 err = norm(dcm_b2n' * dcm_b2n - eye(3)); assert(err < 1e-10, '旋转矩阵不正交，检查角度顺序');

angle2dcm的三个角度输入单位是弧度，顺序参数写错会让滚转和偏航对调，飞机在动画里会侧着飞。dcmbody2wind的输入顺序是迎角在前、侧滑角在后，跟有些教材反过来，这一点手册 2-5 页写得很清楚。最后那行正交性检查看着多余，但它能挡住大部分「矩阵拼错但结果看着差不多」的情况。真到排错时，先看err，比盯着一堆角度值猜要快。

从 NED 往 ECEF 走，会用到dcm2latlon这类跟经纬度挂钩的函数；再往 ECI 走，手册给的是dcmeci2ecef，输入一个 UTC 时间点，返回 ECI 到 ECEF 的旋转矩阵。这条链的顺序不能乱：先机体系到 NED，再 NED 到 ECEF，最后 ECEF 到 ECI，反过来算轨迹回放时会差一个地球自转角。

2.3 单位与常量：ft 和 m 之间别来回手算

手册 2-10 页的「Aerospace Units」单独列了一节，不是凑页数。DATCOM 输出默认是英尺、节、 slugs 这一套，而 MATLAB 里大部分气动公式按米、米每秒、千克写。我见过最典型的一次翻车，是直接把 DATCOM 的升力系数曲线拿去乘按米算的动压，结果升力量级大了十倍，排查了一下午才发现是参考面积单位没换。

工具箱给了几个换算函数，convlength、convvel、convmass、convforce，输入输出都带单位字符串，比自己乘 0.3048 稳。

% 进模型前统一到 SI alt_ft = 35000; alt_m = convlength(alt_ft, 'ft', 'm'); v_kts = 250; v_mps = convvel(v_kts, 'kts', 'm/s'); m_slug = 0.8; m_kg = convmass(m_slug, 'slug', 'kg'); % 量级自检：巡航高度应在 10 km 量级，空速在 100 m/s 量级 assert(alt_m > 8e3 && alt_m < 15e3, '高度单位可能没换'); assert(v_mps > 80 && v_mps < 300, '空速单位可能没换');

这几个assert不是形式主义，单位混用最麻烦的地方是程序不报错，只是结果偏。把量级检查写在数据入口，比事后对着曲线怀疑人生强。手册在单位这节还列了角度、温度、压力的换算，做发动机或风洞示例时会用到，翻到 5-46 页那几个案例前先看一眼有好处。

2.4 DATCOM 导入：datcomimport 的输入输出与缺失块

DATCOM 那块是手册里被翻得最多的部分之一，2-11 到 2-13 页讲了文件格式、导入后的数据结构，以及「缺失数据」怎么处理。datcomimport返回的是 cell 数组，每个元素对应文件里的一个配置，字段名跟着 DATCOM 输出块走。

% 读取 DATCOM 输出文件，第三个参数控制日志详细程度 data = datcomimport('demo_aircraft.dcm', true, 0); case1 = data{1}; % 查看这个文件里记录了哪些马赫数与迎角 machList = case1.mach; alphaList = case1.alpha; % 取马赫数 0.6 那一档的升力、阻力系数 idx = find(abs(machList - 0.6) < 1e-6, 1); CL = case1.cl(idx, :); CD = case1.cd(idx, :); % 画一条升力曲线，先确认数据不是空的 plot(alphaList, CL, 'o-'); xlabel('迎角 (deg)'); ylabel('C_L');

第二个参数usefil控制是否沿用文件里的配置名，第三个参数是日志开关。排查导入问题时把它改成 1，能看到解析到第几个 CASE 块。字段名和 DATCOM 的输出块一一对应，.cl、.cd、.cY这些是常用项，但不同版本的输出文件字段可能略有差异，导入后先fieldnames(case1)看一眼最稳。

缺失数据是这块的主要坑。手册 2-13 页讲了：当某个马赫数或迎角组合在文件里没有对应块时，字段会是空数组或 NaN，直接绘图会出现断点。处理方式有两种，一是用相邻马赫数线性插值补齐，二是把分析范围收窄到有数据的区间。手册里还提了外推的边界，迎角超出风洞数据范围太远时，插值出来的气动系数没有物理意义，这一点在做大迎角机动分析时尤其要注意。我一般会先画一张「马赫数 × 迎角」的覆盖图，把有数据的格子标出来，再决定分析区间。

3. 飞行可视化：Aero.Animation 与 FlightGear 接口怎么跑通

3.1 Aero.Animation 的对象模型

手册 2-22 到 2-30 页讲动画对象，核心是Aero.Animation。它不是一个函数，而是一个句柄对象，里面挂着若干Bodies，每个 Body 有几何文件、轨迹源、坐标源这几组属性。理解这个结构之后，很多「窗口打开了但飞机不动」的问题就有着落点了。

几何文件一般分两部分：.mat里放顶点和面片数据，.ac3是 AC3D 格式的模型。轨迹源可以是 MAT、CSV 或自定义函数，手册里给了TimeSeriesSourceType的几种取值。坐标源决定机体在场景里的位置和姿态，跟轨迹源可以分开设置，这一点在做「同一段轨迹、不同视角回放」时很有用。

3.2 轨迹数据准备与最小可运行例子

先把轨迹准备成时间列在前的 CSV，第一列是秒，后面依次是位置和姿态。下面这段是最小可运行版本，跑通之后再往里加相机、透明度这些。

% 创建动画对象并设定播放节奏 anim = Aero.Animation; anim.FramesPerSecond = 20; anim.TimeScaling = 2; % 2 倍速播放 % 绑定几何与轨迹文件，路径尽量用绝对路径 anim.createBody('aircraft.mat', 'aircraft.ac3'); anim.Bodies{1}.TimeSeriesSource = 'trajectory.csv'; anim.Bodies{1}.TimeSeriesSourceType = 'CSV'; anim.Bodies{1}.CoordSource = 'trajectory.csv'; anim.Bodies{1}.CoordSourceType = 'CSV'; % 打开窗口并播放 anim.show(); anim.play();

FramesPerSecond控制的是渲染帧率，不是数据采样率。轨迹采样密、帧率低，回放会跳；轨迹采样稀、帧率高，插值出来的姿态可能抖。我的经验是先把FramesPerSecond设在 20 左右，TimeScaling设 1，确认轨迹本身没问题再加速。createBody的两个参数顺序是 MAT 在前、AC3 在后，写反了不会立刻报错，但窗口里可能什么都不显示。

轨迹 CSV 的时间列必须是单调递增的秒值。有些飞控日志导出来是毫秒，直接喂进去，动画会在开头几帧就播完，看着像「飞机不动」，其实是时间轴跑完了。手册 2-24 页讲「Loading Recorded Data for Flight Trajectories」时提了时间序列对象的构造方式，可以先转成timeseries再挂上去，方便检查时间轴。

3.3 FlightGear 接口：配置、运行脚本与回放

FlightGear 那条路比内置动画麻烦，但场景真实感强，做演示时更拿得出手。手册 2-42 到 2-49 页讲了接口配置、运行脚本生成和轨迹回放。核心对象是Aero.FlightGearAnimation，用法跟Aero.Animation类似，多了一步「生成启动脚本」。

% 构造 FlightGear 动画对象 fg = Aero.FlightGearAnimation; fg.FramesPerSecond = 10; fg.TimeScaling = 1; % 绑定几何与轨迹 fg.createBody('aircraft.mat', 'aircraft.ac3'); fg.Bodies{1}.TimeSeriesSource = 'trajectory.csv'; fg.Bodies{1}.TimeSeriesSourceType = 'CSV'; % 生成启动脚本，Windows 下是 .bat，Linux 下是 .sh fg.generateRunScript('run_fg.bat');

generateRunScript会把 FlightGear 需要的命令行参数写进脚本，包括机场、初始位置、网络端口、飞机模型路径。先生成脚本、再手动双击启动，比让 MATLAB 直接拉起 FlightGear 稳，因为你能看到 FlightGear 自己的报错窗口。端口对不上是这条链路最常见的问题：MATLAB 侧发数据，FlightGear 侧没收到，两边都不报错，只是飞机停在跑道上。

手册 2-43 页讲了「Install and Start FlightGear」的步骤，2-44 页讲怎么装额外的地景。如果场景加载失败，先确认 FlightGear 自带的默认机场能不能正常打开，再查脚本里的机场代码和跑道号。另外，FramesPerSecond在 FlightGear 这条路上不建议设太高，10 到 15 之间比较稳，网络传输和渲染都吃得消。

3.4 飞行仪表组件挂到 App Designer

手册 2-54 到 2-58 页讲的是把飞行仪表组件拖进 App Designer。R2023b 里这些组件在 App Designer 的组件库里能直接找到，空速表、高度表、姿态仪都有对应函数，也可以在代码里创建。

% 在 uifigure 里放空速表与高度表 fig = uifigure('Name', 'Flight Instruments'); spd = uiaeroairspeed(fig); spd.Position = [20 20 180 180]; alt = uiaeroaltimeter(fig); alt.Position = [220 20 180 180]; % 数据回放时更新指针 spd.Value = 250; % 单位 knots alt.Value = 35000; % 单位 ft

这两个组件的Value单位是固定的，空速表默认节，高度表默认英尺，从 SI 数据源喂进来之前要换一次。手册在 2-52 页讲「Create Flight Instrument Components」时列了每个组件的可用属性，做布局时先看一眼能省不少调整位置的时间。App Designer 里拖组件的方式和代码创建等价，区别是拖进去之后属性面板能直接改，调试阶段更直观。

4. 卫星场景与轨道要素：TLE/OMM 到可视化场景的落地路径

4.1 轨道六要素与 keplerian2ijk/ijk2keplerian

手册 2-65 到 2-73 页讲卫星场景的关键概念，轨道要素那节把半长轴、偏心率、倾角、升交点赤经、近地点幅角、真近点角这六个量挨个解释了一遍。做轨道可视化之前，先在这两个函数之间来回转一次，确认自己对角度约定的理解没错。

% 轨道六要素转 ECI 位置速度 a = 7000; % 半长轴，km e = 0.001; % 偏心率 inc = deg2rad(51.6); % 轨道倾角 raan = deg2rad(30); % 升交点赤经 argp = deg2rad(45); % 近地点幅角 nu = deg2rad(10); % 真近点角 mu = 398600.4418; % 地球引力常数，km^3/s^2 [r, v] = keplerian2ijk(a, e, inc, raan, argp, nu, mu); % 反算前三个要素，检查转换是否自洽 [a2, e2, i2] = ijk2keplerian(r, v, mu); assert(abs(a2 - a) < 1e-6 && abs(i2 - inc) < 1e-6, '轨道要素转换不自洽');

keplerian2ijk的角度输入是弧度，第七个参数可以省略，省略时用默认地球引力常数。ijk2keplerian返回多个输出，只取前三个做检查就够了。这段代码的价值不在计算本身，而在建立「角度单位是弧度、距离单位是公里」的肌肉记忆。后面加载 TLE 时，TLE 里的倾角、升交点赤经都是度，直接喂进去会得到一条完全不对的轨道，而且不报错。

4.2 satelliteScenario 加载 TLE 与 OMM

手册 2-75 到 2-78 页讲场景的创建和播放。satelliteScenario的三个输入是起始时间、结束时间、采样步长，时间用datetime且显式指定 UTC 时区。

% 场景时间与采样步长 startTime = datetime(2024, 3, 1, 0, 0, 0, 'TimeZone', 'UTC'); stopTime = startTime + hours(6); sc = satelliteScenario(startTime, stopTime, 30); % 从 TLE 文件批量加载卫星 sats = satellite(sc, 'constellation.tle'); % 添加地面站，经纬度用示例值 gs = groundStation(sc, 10.0, 20.0, 'Name', 'StationA'); % 计算第一颗星对地面站的访问窗口 ac = access(sats(1), gs); intervals = ac.Intervals; % 播放三维场景 play(sc);

satellite函数接受 TLE 文件路径，也接受 OMM 文件，手册 2-73 页专门讲了 OMM 格式的字段。TLE 的历元时间很关键：场景起始时间离 TLE 历元太远，轨道预报误差会累积，卫星位置会跳。做演示时把场景时长控制在 TLE 历元前后几天内，位置看着才顺。access计算的是可见性窗口，Intervals里存的是起止时间对，导出到表格之前先确认时区还是 UTC。

4.3 地面站与访问分析

地面站的经纬高三个输入里，高度容易被漏掉。手册在讲groundStation时给了默认值，但做低轨卫星可见性分析时，地面站海拔对仰角门限有影响，尤其在高纬度地区。访问分析的结果除了时间窗口，还能拿到最大仰角、距离这些量，做链路预算时用得上。

手册 2-78 页的「Model, Visualize, and Analyze Satellite Scenario」示例里把这几步串了一遍，照着跑一次比看文字快。场景播放时如果卫星轨迹看着像折线，把采样步长从 60 秒降到 30 秒或 10 秒，视觉上会顺很多，代价是内存占用上升。

4.4 星历与大地水准面附加数据

手册第 3 章单独讲了「Add Ephemeris and Geoid Data for Aerospace Products」。这部分数据不是默认安装的，要单独装附加数据包。装完之后，行星历表、大地水准面、磁场模型这些示例才能跑。

这些示例的共同点是依赖外部数据文件，路径没配好就会报「找不到文件」。我的做法是把数据包装完之后，先跑手册里最短的那个示例，确认路径通了，再改参数做自己的分析。磁场模型和大地水准面这两个示例还涉及版本年份，2020 历元和 1996 模型不要混用，否则等值线画出来跟手册对不上。

5. 避坑与排查：手册目录里没写透的 5 个翻车点

5.1 现象：datcomimport 读进来字段全空或全是 NaN

原因通常有两个：.dcm文件里对应马赫数或迎角的 CASE 块本来就没输出，或者固定列宽格式在解析时列对不齐。先打开文件确认有没有那个 CASE，再把datcomimport的第三个参数改成 1 看解析日志，日志会告诉你读到第几块、哪一块被跳过。缺失块用相邻马赫数线性插值补齐，但外推范围不要超过数据边界太远，手册 2-13 页对这一点有明确提醒。

5.2 现象：动画窗口打开了，机体停在原点不动

九成是轨迹文件的时间列没被正确识别。先检查 CSV 第一列是不是单调递增的秒值，如果是毫秒，整列除以 1000 再试。其次看TimeSeriesSourceType有没有设成'CSV'，设成'MAT'但文件是 CSV，对象不会报错，只是读不到数据。还有一种情况是CoordSource和TimeSeriesSource指向了不同文件，位置和姿态对不上，机体在原地打转。

5.3 现象：FlightGear 启动了，但飞机不动或轨迹乱飞

先看端口。generateRunScript生成的脚本里写了网络端口，MATLAB 侧和 FlightGear 侧必须一致，改了一边忘了另一边是常事。再看机场和跑道号，场景加载失败时 FlightGear 会退回默认位置，飞机看着像在原地。最后查飞机模型的路径，.ac3文件不在 FlightGear 的 aircraft 目录下，界面会显示默认模型，但轨迹数据其实在发。排查顺序建议从「手动启动 FlightGear 能不能正常飞默认机型」开始，这一步通了再查数据链路。

5.4 现象：卫星场景播放时位置突然跳变

时间基准问题。TLE 的历元是 UTC，场景起始时间如果用了本地时区，两者差几个小时，轨道预报就会偏。构造datetime时显式写'TimeZone', 'UTC'，场景时长控制在 TLE 历元前后几天内。还有一种跳变来自采样步长太大，卫星在两步之间走过一段弧线，视觉上像瞬移，把步长降到 10 秒或 30 秒就顺了。

5.5 现象：静稳定性或气动系数量级离谱

单位混用。DATCOM 输出是英尺、节、slug 这一套，模型按米、米每秒、千克写，中间少换一次，系数就差一个量级。另一个容易漏的是参考面积和参考弦长，DATCOM 文件里给的是平方英尺和英尺，模型里如果按平方米和米填，力矩系数会对不上。进模型前统一用convlength、convvel、convmass过一遍，再写两个assert卡量级，比事后对着曲线猜原因省时间。

6. 进阶：用函数和对象两条路做固定翼静稳定性与线性分析

手册第 2 章末尾给了两条固定翼分析的路径，一条是函数工作流，一条是对象工作流，对应的章节分别是「Work with Fixed-Wing Aircraft Using Functions」和「Analyze Fixed-Wing Aircraft with Objects」。两条路的差别不在算得准不准，而在复用方式。

对象工作流的骨架大致是这样，部件数据从手册示例的 MAT 文件里取：

% 载入示例几何与气动参数 S = load(' <p> <a href="https://download.csdn.net/download/fengtianxiaoyao/88382725" style="color:#ec7500;font-size:14px;"> 本文还有配套的精品资源，点击获取 </a> <img alt="menu-r.4af5f7ec.gif" src="https://csdnimg.cn/release/wenkucmsfe/public/img/menu-r.4af5f7ec.gif" style="width:16px;margin-left:4px;vertical-align:text-bottom;cursor:text;"> </p>