第二节 基础数学预备：向量、矩阵、点乘、叉乘、行列式、特征值等基本概念回顾

气动坐标系转换的核心，本质上是一系列数学运算的工程化封装。无论是从机体坐标系到地面坐标系的姿态解算，还是从传感器坐标系到执行机构坐标系的力与力矩映射，其底层都离不开向量与矩阵的运算规则。本节不追求纯数学的完备性推导，而是从工程应用视角，系统梳理后续章节将反复调用的数学工具：向量的几何意义与代数表达、矩阵作为线性变换的载体、点乘与叉乘的物理含义、行列式对可解性的判定，以及特征值在系统稳定性分析中的角色。这些概念并非孤立存在，它们共同构成了坐标转换的“语法规则”。

2.1 向量：从几何箭头到代数坐标

核心原理

在气动坐标系中，向量是最基本的描述单元。一个向量 v 同时包含大小（模长）和方向两个属性。在三维笛卡尔坐标系下，向量通常表示为列矩阵形式：v = [v_x, v_y, v_z]^T。这里的下标对应坐标轴上的投影分量。例如，来流速度向量可以分解为机体轴系下的三个速度分量 (u, v, w)，分别代表前向、侧向和垂向速度。向量的加法满足平行四边形法则，数乘则改变其大小或反向。在坐标转换中，向量本身不随坐标系旋转而改变其物理意义，但其坐标分量会随参考系变化——这正是我们需要转换矩阵的原因。

判断方法或操作流程

在实际工程中，处理向量时需遵循以下步骤：

1. 明确参考坐标系：向量必须标注其所在的坐标系（如地面系 E、机体系 B、风轴系 W），否则数值无意义。
2. 分量分解：将向量投影到所选坐标系的三个正交轴上，得到有序三元组。
3. 模长与单位向量：计算模长 ||v|| = √(v_x² + v_y² + v_z²)，单位向量 v̂ = v / ||v||，用于表示纯方向。
4. 向量运算：加减法对应分量运算，数乘对应每个分量缩放。

案例拆解或实战场景

以空速管测量为例，空速管安装在机头，测得的是沿机体轴 x_b 方向的气流速度分量 u。但气动分析需要的是相对于风轴的速度向量。假设飞机以 50 m/s 的速度平飞，迎角 α = 5°，侧滑角 β = 2°。则来流速度在机体轴系下的分量为：u = V·cosα·cosβ ≈ 49.7 m/s，v = V·sinβ ≈ 1.7 m/s，w = V·sinα·cosβ ≈ 4.4 m/s。这个向量 V_b = [49.7, 1.7, 4.4]^T 就是后续气动力计算的输入。若直接使用标量速度 50 m/s 而忽略方向，则升力、阻力的计算将产生显著误差。

常见误区或总结复盘

2.2 矩阵：线性变换的代数载体

核心原理

矩阵是向量的“加工厂”。一个 m×n 矩阵 A 乘以一个 n 维向量 v，得到 m 维向量 w = Av。从几何角度看，矩阵乘法实现了对向量的线性变换：旋转、缩放、剪切、投影等。在气动坐标转换中，最核心的矩阵是旋转矩阵 R，它属于特殊正交群 SO(3)，满足 R^TR = I 且 det(R) = +1。这意味着旋转矩阵保持向量长度不变，且不改变手性（右手坐标系仍为右手系）。

判断方法或操作流程

构造和使用旋转矩阵时，需遵循以下流程：

• 确定旋转顺序：气动领域常用 Z-Y-X 欧拉角顺序（偏航-俯仰-滚转），对应矩阵 R = R_x(φ)·R_y(θ)·R_z(ψ)。
• 写出基本旋转矩阵：绕 x 轴旋转 φ 的矩阵为 [[1,0,0],[0,cosφ,-sinφ],[0,sinφ,cosφ]]，绕 y 轴、z 轴类似。
• 矩阵乘法顺序：注意矩阵乘法不满足交换律，R_x·R_y ≠ R_y·R_x。必须严格按照定义的旋转顺序左乘。
• 验证正交性：计算 R^TR 是否等于单位阵，若偏差大于 1e-6，说明矩阵构造有误。

案例拆解或实战场景

考虑一个无人机从地面系到机体系的转换。地面系定义为北-东-地（NED），机体系定义为前-右-下（FRD）。假设无人机当前偏航角 ψ = 30°，俯仰角 θ = 10°，滚转角 φ = 0°。则从地面系到机体系的旋转矩阵为：

R_B2E = R_x(0)·R_y(10°)·R_z(30°)
= [[cos10·cos30, cos10·sin30, -sin10],
   [sin30,        cos30,        0],
   [sin10·cos30,  sin10·sin30,  cos10]]

代入数值后，若地面系下有一个速度向量 V_E = [20, 0, 0]^T m/s（纯北向），则机体系下的分量为 V_B = R_B2E·V_E ≈ [17.1, 10.0, 3.5]^T m/s。这意味着纯北向的风在机体坐标系下表现为前向 17.1 m/s、右侧 10 m/s、向下 3.5 m/s 的分量组合。

常见误区或总结复盘

2.3 点乘与叉乘：投影与法向量的工程意义

核心原理

点乘（内积）a·b = ||a||·||b||·cosθ，结果是一个标量，表示 b 在 a 方向上的投影长度乘以 a 的模长。叉乘（外积）a×b 的结果是一个向量，其方向垂直于 a 和 b 所在的平面，遵循右手定则，模长等于 ||a||·||b||·sinθ。在气动坐标转换中，点乘用于计算力在某个方向上的分量（如升力垂直于来流方向），叉乘用于计算力矩（力臂 × 力）或确定旋转轴。

判断方法或操作流程

• 点乘应用：当需要求一个向量在另一个向量方向上的投影时，使用点乘。例如，求气动力在机体轴 x 方向的分量：F_x = F·i_b，其中 i_{b 是机体轴 x 方向的单位向量。}
• 叉乘应用：当需要求垂直于两个向量的方向，或计算力矩时，使用叉乘。力矩 M = r×F，其中 r 是力作用点相对于参考点的位置向量。
• 叉乘的矩阵形式：叉乘可以写成反对称矩阵与向量的乘法：a×b = [a]_×·b，其中 [a]_× 是由 a 的分量构成的 3×3 反对称矩阵。这在推导旋转矩阵的导数时非常有用。

案例拆解或实战场景

以飞机气动焦点计算为例。气动焦点是气动力增量的作用点。假设机翼在机体坐标系下的位置向量为 r = [2.0, 0, 0.3]^T m（从重心指向气动中心），气动力增量 ΔF = [0, 0, -500]^T N（垂直向下）。则气动力增量对重心产生的俯仰力矩为：

ΔM = r × ΔF = [2.0, 0, 0.3] × [0, 0, -500]
= [0·(-500) - 0.3·0, 0.3·0 - 2.0·(-500), 2.0·0 - 0·0]
= [0, 1000, 0]ᵀ N·m

结果 [0, 1000, 0]^T 表示产生了一个 +1000 N·m 的俯仰力矩（抬头力矩）。这个计算直接使用了叉乘，且避免了三角函数，计算效率高。

常见误区或总结复盘

2.4 行列式与特征值：矩阵的“体检报告”

核心原理

行列式 det(A) 是一个标量，它衡量矩阵 A 对空间体积的缩放比例。对于旋转矩阵，det(R) = 1，表示体积不变且方向不变（右手系）。若 det(A) = 0，则矩阵奇异，不可逆，意味着该变换将空间压缩到更低维度，信息丢失。特征值 λ 和特征向量 v 满足 Av = λv，即经过变换后，特征向量方向不变，仅被缩放 λ 倍。在气动坐标转换中，特征值分解常用于分析旋转矩阵的等效旋转轴和旋转角，也用于判断线性系统的稳定性。

判断方法或操作流程

矩阵性质	行列式含义	特征值含义	工程应用
旋转矩阵 R	det(R) = +1	λ = 1, e±iθ	提取等效旋转轴与角度
缩放矩阵 S	det(S) = 各缩放因子乘积	λ = 各缩放因子	判断各轴增益
投影矩阵 P	det(P) = 0	λ = 0 或 1	识别信息丢失维度
系统矩阵 A	det(A) ≠ 0 保证可控性	实部为负 → 稳定	飞行动力学稳定性分析

案例拆解或实战场景

在飞行动力学中，纵向运动的系统矩阵 A 的特征值决定了短周期和长周期模态的稳定性。假设某飞机纵向状态矩阵的特征值为 λ_1,2 = -0.5 ± 2.5i（短周期模态）和 λ_3,4 = -0.02 ± 0.1i（长周期模态）。所有特征值实部均为负，因此该飞机纵向运动是稳定的。实部值越大，模态衰减越快。短周期模态实部 -0.5 意味着扰动在约 2 秒内衰减到 37%，而长周期模态需要约 50 秒。这个分析直接指导了飞行控制律的设计——需要增加阻尼以加快长周期模态的收敛。

常见误区或总结复盘

2.5 综合应用：从数学工具到坐标转换闭环

核心原理

将上述工具串联起来，就构成了气动坐标转换的完整闭环。首先，用向量描述物理量（速度、力、力矩）；其次，用旋转矩阵实现不同坐标系之间的变换；再次，用点乘和叉乘计算投影与力矩；最后，用行列式和特征值验证变换的正确性与系统稳定性。这个闭环在飞控代码中体现为：传感器数据（向量）→ 姿态矩阵（旋转矩阵）→ 导航解算（点乘/叉乘）→ 稳定性监控（特征值）。

操作流程清单

1. 定义坐标系：明确地面系、机体系、风轴系的定义与原点。
2. 构建旋转矩阵：根据欧拉角或四元数生成 R，并验证 R^TR = I。
3. 向量变换：v_B = R_B2E·v_E，注意平移向量 t 的处理。
4. 计算气动力/力矩：使用点乘求分量，叉乘求力矩。
5. 稳定性校验：计算系统矩阵特征值，判断实部符号。

案例拆解或实战场景

以某型无人机风洞试验数据处理为例。试验中，天平测得机体轴系下的气动力分量 F_x = 120 N，F_y = -15 N，F_z = -800 N。已知试验时迎角 α = 8°，侧滑角 β = 0°。需要将这些力转换到风轴系（x 轴沿来流方向，z 轴在对称面内垂直于来流）。转换矩阵为从机体轴到风轴的旋转：R_W2B = R_y(-α)·R_z(β)。代入 α = 8°，β = 0°，得到：

R_W2B = [[cos8, 0, sin8],
         [0,    1, 0],
         [-sin8, 0, cos8]]

风轴系下的力为 F_W = R_W2B·F_B。计算得：阻力 D = F_x·cos8 + F_z·sin8 ≈ 120·0.990 + (-800)·0.139 ≈ 118.8 - 111.2 = 7.6 N；升力 L = -F_x·sin8 + F_z·cos8 ≈ -120·0.139 + (-800)·0.990 ≈ -16.7 - 792 = -808.7 N（负号表示向下，实际升力向上为 +808.7 N）。这个转换直接使用了旋转矩阵和点乘，且验证了 det(R) = 1。

常见误区或总结复盘