五章共 70 余条编号公式、9 张 DH 参数表、3 套逆解公式——它们散在 40 余个小节里。本页把它们按「表示 → 传递 → 反解 → 微分 → 力 → 驱动与测量」的顺序压成一页,每条后面给出正文出处,供做题与复习时单点查阅。
怎么用这一页
- 做题卡壳时:按下面的分节标题定位,不用翻正文;括号里的链接跳回推导与上下文。
- 考前自测时:遮住右栏,只看「要算什么」,默写公式,再对答案。
- 符号约定:cϑ=cosϑ、sϑ=sinϑ、c12=cos(ϑ1+ϑ2);上下标读法「上标是参考系、下标是被描述对象」(Rji = 系 j 相对系 i),见符号表。
| 要算什么 | 公式 | 出处 |
| 正交性 | RTR=I,即 R−1=RT | 定义 2.2.2 |
| 行列式 | detR=+1(=−1 是镜像,不是旋转) | 定义 2.2.2 |
| 独立参数 | 9 量 −6 约束 =3 | 定义 2.2.1 |
| 列的含义 | R 的第 j 列 = 系 j 的 xj 轴在系 i 中的方向余弦 | 定义 2.2.1 |
三重身份(同一矩阵、三种读法):姿态描述 / 坐标变换 / 旋转算子——定理 2.2.B。
Rx(ϑ)=1000cϑsϑ0−sϑcϑ,Ry(ϑ)=cϑ0−sϑ010sϑ0cϑ,Rz(ϑ)=cϑsϑ0−sϑcϑ0001
| 情形 | 做法 | 出处 |
| 绕当前轴(动系)连续转 | 右乘 R=R1R2R3 | 原理 2.3.A |
| 绕固定轴(定系)连续转 | 左乘 R=R3R2R1 | 原理 2.3.B |
| 旋转不可交换 | R1R2=R2R1 | 例 2.3 |
口诀:右乘当前轴、左乘固定轴。ZYZ / RPY 的对偶、DH 的 Ai 连乘、角和轴的五步序列,全是这一条的换装登场。
R(ϕ)=Rz(φ)Ry′(ϑ)Rz′′(ψ)=cφcϑcψ−sφsψsφcϑcψ+cφsψ−sϑcψ−cφcϑsψ−sφcψ−sφcϑsψ+cφcψsϑsψcφsϑsφsϑcϑ
逆解(sϑ=0,原书式 2.19;两支等价):
φ=Atan2(r23,r13),ϑ=Atan2(r132+r232, r33),ψ=Atan2(r32,−r31)
第二支:三个宗量整体变号(ϑ 取 (−π,0))。奇点 sϑ=0:只能定 φ±ψ,见 定理 2.4.B。
R(ϕ)=Rz(φ)Ry(ϑ)Rx(ψ)=cφcϑsφcϑ−sϑcφsϑsψ−sφcψsφsϑsψ+cφcψcϑsψcφsϑcψ+sφsψsφsϑcψ−cφsψcϑcψ
逆解(ϑ∈(−π/2,π/2)):
φ=Atan2(r21,r11),ϑ=Atan2(−r31,r112+r212),ψ=Atan2(r32,r33)
奇点 cϑ=0(俯仰 ±90°)。角名:roll 绕纵轴 x、pitch 绕 y、yaw 绕竖轴 z。
为什么反解一律用 Atan2:arccos 的值域只有 [0,π](不唯一),且导数 (arccosx)′=−1/1−x2 在 x→±1 处无界(误差被无限放大);(arctanx)′=1/(1+x2)≤1 处处有界。见 原理 2.4.A。
R(ϑ,r)=rx2(1−cϑ)+cϑrxry(1−cϑ)+rzsϑrxrz(1−cϑ)−rysϑrxry(1−cϑ)−rzsϑry2(1−cϑ)+cϑryrz(1−cϑ)+rxsϑrxrz(1−cϑ)+rysϑryrz(1−cϑ)−rxsϑrz2(1−cϑ)+cϑ
| 要算什么 | 公式 | 出处 |
| 合成路径 | Rz(α)Ry(β)Rz(ϑ)Ry(−β)Rz(−α)(全左乘) | 定义 2.5.1 |
| 紧凑形式 | R=cϑI+(1−cϑ)rrT+sϑS(r) | 证明 |
| 由 R 求角(取迹) | cosϑ=2r11+r22+r33−1 | 定义 2.5.1 |
| 由 R 求轴(取反对称部分) | r=2sϑ1[r32−r23, r13−r31, r21−r12]T | 定义 2.5.1 |
| 指数映射 | R(ϑ,r)=exp(ϑS(r))=I+sinϑS+(1−cosϑ)S2(幂次三循环 S2=rrT−I、S3=−S) | 定理 2.5.D |
| 对数映射 | logR=ϑS(r)=2sϑϑ(R−RT) | 定理 2.5.E |
| 不唯一 | R(−ϑ,−r)=R(ϑ,r) | 定理 2.5.A |
| 失效 | sϑ=0:ϑ=0 轴任意;ϑ=π 反对称部分为零(对角线 ri2=(rii+1)/2 另解) | 定理 2.5.B |
η=cos2ϑ,ϵ=sin2ϑr,η2+∣ϵ∣2=1
R(η,ϵ)=2(η2+ϵx2)−12(ϵxϵy+ηϵz)2(ϵxϵz−ηϵy)2(ϵxϵy−ηϵz)2(η2+ϵy2)−12(ϵyϵz+ηϵx)2(ϵxϵz+ηϵy)2(ϵyϵz−ηϵx)2(η2+ϵz2)−1
η=21r11+r22+r33+1,ϵ=21sgn(r32−r23)r11−r22−r33+1sgn(r13−r31)r22−r33−r11+1sgn(r21−r12)r33−r11−r22+1
逆 Q−1={η,−ϵ};乘积 Q1⊗Q2={η1η2−ϵ1Tϵ2, η1ϵ2+η2ϵ1+ϵ1×ϵ2}。见 定义 2.6.1、原理 2.6.A。
p~=[p1],A10=[R100To101],A01=(A10)−1=[R10T0T−R10To101]
连乘:p~0=A10A21⋯Ann−1p~n(原理 2.7.A)。
两条纪律:A−1=AT(正交性只在旋转块内成立);方向向量(自由向量)的变换只用左上块 R,不吃平移。
Aii−1(qi)=cϑisϑi00−sϑicαicϑicαisαi0sϑisαi−cϑisαicαi0aicϑiaisϑidi1
两步分解:先「沿 zi−1 平移 di + 绕 zi−1 转 ϑi」,再「沿 xi 平移 ai + 绕 xi 转 αi」,右乘合成(Aii−1:从建系几何直接读出两步分解)。
第 4 列就是三参数的「肉身」:(aicϑi, aisϑi, di);旋转块 = 绕 z 转 ϑ ⊕ 绕 x 转 α。
记忆主线 · 看守模型(先把「系到底挂在哪」说透:看守模型):系 i−1 看守关节 i——zi−1 就是关节 i 的轴线;杆 i 固连系 i(原点在关节 i+1 轴上)。建系四步实例见 例解。
| 规则 | 服务的参数 | 为什么 |
| zi 沿关节 i+1 的轴线 | ϑi,di(变量) | 关节运动就是「绕/沿自己的轴」,在 A 里恰表现为绕 zi−1 转 / 沿 zi−1 移 |
| 原点 oi 取在公垂线交点、xi 沿公垂线 | ai(常数) | 两根 z 轴的相对位置由公垂线一锤定音 |
| 同一对 z 轴的方向关系 | αi(常数) | 方向差只剩绕 xi 的一个转角 |
| yi 右手系补齐 | — | 只保证 R∈SO(3) |
总原则:能自由选取的地方「让参数尽量为零」(总原则:让参数尽量为零)——2.9 各表里成片的 0 就是这条的兑现。
3 常 1 变:ai,αi 恒为常数;di,ϑi 中恰有一个是关节变量 qi(转动关节 qi=ϑi、移动关节 qi=di)——定义 2.8.3。
- 关节轴编号,定 z0,…,zn−1;2. 选系 0;3. 对 i=1..n−1 定 oi;
- 沿公垂线取 xi;5. 右手系补 yi;6. 选系 n;7. 列 4n 参数表;
- 逐个算 Aii−1;9. 连乘得 Tn0(q);10. 串常值变换 Tb0,Ten 得 Teb。
课堂口诀:编号、建系(能零则零)、列参数表、代入、连乘。
xe=k(q),Teb(q)=Tb0Tn0(q)Ten,Tn0=A10A21⋯Ann−1
| 概念 | 内容 | 出处 |
| 可达 ⊇ 灵活 | 可达 = 至少一个方向能到;灵活 = 任意方向都能到 | 定义 2.10.2 |
| 精度 ≠ 重复精度 | 精度 = 指令位姿与真实位姿的差;重复精度 = 多次回同一点的散布 | 定义 2.10.2 |
| 冗余判据 | n 关节数、m 操作空间维数、r 任务实际用到的维数;n>r 即冗余 | 定义 2.10.3 |
Δx=ΦΔξ,Δξ=(Φ∗TΦ∗)−1Φ∗TΔx∗,ξ′=ξ+Δξ
Φ 为 m×4n 标定矩阵,Δξ 为 4n 维参数偏差。量测条件 tm≥4n(宜远大于,防 Φ 病态);因 k 非线性,须迭代。见 标定三步法。
四难点(定理 2.12.A):非线性(解方程而非算函数)、多解、可能无解 / 无穷多解、无通用解析法。
闭合解门槛(定理 2.12.B):球形腕(三轴交一点)或三轴平行 ⇒ 可解耦。
解耦五步(解耦五步)——这是全书最实用的一套流程:
- 由末端位姿反推腕心:
pw=pe−d6ae
(d6 = 腕心到末端原点的「外移长度」,沿接近向量 ae 反推。)
- 解臂:由 pw 求 (q1,q2,q3)(pw 只依赖前三个关节);
- 算 R30=R10R21R32;
- 解腕:由 (R30)TRe=R63 反解 (q4,q5,q6);
- 整理全部解支。
腕 = ZYZ 反解:球形腕的姿态部分与 ZYZ 欧拉角反解同构(A4 的公式直接可用)。
数解支的口诀:位置解给「臂」的支数(平面臂 2 支 = 肘上 / 肘下,拟人臂 4 支),方向解给「腕」的支数(球腕 2 支 = 翻转),总数 = 臂 × 腕。
[p˙eωe]=J(q)q˙,J=[JPJO]∈R6×n
R˙=S(ω)R,S(ω)=0ωz−ωy−ωz0ωxωy−ωx0,S(ω)p≡ω×p
推导链:RRT=I 求导 ⇒ R˙RT 反对称 ⇒ 用常模向量比对 ω×r 认出 S=S(ω)(定理 3.1.A)。共轭关系 RS(ω)RT=S(Rω)(定理 3.1.B,已在 2.5 作为工具证明)。李代数读法:ω=(R˙RT)∨(基系读数),RTR˙ 给刚体系读数(李代数视角:R˙RT 与 RTR˙ 各读出谁的角速度)——「角速度 = so(3) 的元素」(2.5.3)。
J=[J1⋯Jn],Ji=[JPiJOi],{关节 i 转动:关节 i 移动:JPi=zi−1×(pe−pi−1),JOi=zi−1JPi=zi−1,JOi=0
见 定理 3.1.C。所需量全部从正运动学读取:
| 量 | 从哪读 |
| zi−1 | Ti−10 左上 3×3 块的第三列(基系表示,不是 [0,0,1]T!) |
| pi−1 | Ti−10 第四列前三行 |
| pe | Tn0 第四列前三行 |
物理含义:其余关节全部锁死时,单个关节的单位速度对末端速度的贡献。
Ju=[Ru00Ru]J
detJ(q)=0
| 结构 | 判据 | 出处 |
| 两连杆平面臂 | detJ=a1a2s2=0⟺ϑ2=0,π | 定义 3.3.1 |
| 带球腕机械手(解耦) | detJ=detJ11⋅detJ22(J11 只含臂、J22 只含腕) | 定理 3.3.A |
| 腕奇点 | detJ22∝s5 ⇒ ϑ5=0,π | 定理 3.3.A |
| 拟人臂肘 / 肩 | detJ11=−a2a3s3(a2c2+a3c23) ⇒ s3=0(肘)或 a2c2+a3c23=0(肩腕点落轴) | 定理 3.3.A |
n>r ⟹ N(J)={0},P=NNT (秩 n−r),q˙=q˙∗+Pq˙0
Pq˙0 在零空间里,不惊动末端(JPq˙0=0)——用来做内部运动 / 再配置(避障、躲关节极限、远离奇异)。见 定义 3.4.1、定理 3.4.B。
| J 的形状 | 工具 | 公式 | 条件 |
| 方阵满秩 | J−1 | q˙=J−1ve,配数值积分 qk+1=qk+J−1veΔt | n=r、无奇异 |
| 长方(冗余) | J† | q˙=J†ve+(I−J†J)q˙0,J†=JT(JJT)−1 | r<n、行满秩 |
| 退化(奇异邻域) | J∗ | J∗=W−1JT(JW−1JT+K−1)−1;W=I,K=k2I 时 J∗=JT(JJT+k2I)−1 | 容许末端误差 δ |
见 定义 3.5.1、定义 3.5.2。
口诀:方阵求逆、矮胖求伪逆、退化加阻尼。
K 的取法:奇点邻域取大(压住发散),远离奇点取小甚至趋 0(退回伪逆)。
q˙0 的三种取法(都取 k0(∂g/∂q)T,最大化 g):
| 目标 g | 作用 |
| w=det(JJT)(可操作度) | 远离奇异 |
| 关节角到中点的距离 | 躲关节极限 |
| 臂上点到障碍物的最小距离 | 避障 |
ωe=T(φe)φ˙e,T(φ)=001−sφcφ0sϑcφsϑsφcϑ (ZYZ),JA=TA−1J,TA=[I300T]
T 的三列 = 欧拉角三根转轴在基系中的方向(未必正交)。线速度部分 J 与 JA 相同,姿态部分差一个换基。sϑ=0 时 T 奇异 ⇒ JA 失效(与 ZYZ 表示奇点同源)。见 定义 3.6.1、定理 3.6.A。
τ=JT(q)Fe ,Fe=[feμe]∈R6
由虚功原理推出:q˙Tτ=veTFe=(Jq˙)TFe=q˙TJTFe。速度变换的转置就是力变换(对偶),长方阵也成立,无需满秩。见 定义 3.8.1。
两坐标系间的速度 / 力变换:
[p˙10ω10]=[R100−R10S(p~)R10][v1ω1],[f10μ10]=[R10R10S(p~)0R10][f1μ1]
右边(力)恰是左边(速度)的转置——平移块 S(p~) 在速度式里带负号、在力式里移位且不带负号,正是对偶的体现。
三次(4 个系数 ↔ 4 个约束):q(t)=a3t3+a2t2+a1t+a0,由 q(0)=qi, q˙(0)=q˙i, q(tf)=qf, q˙(tf)=q˙f 定系数。
规范化形式(初终速度为零时与数据无关,qi=0,qf=π,tf=1 时 a2=3π,a3=−2π):
q(t)=qi+(qf−qi)[3η2−2η3],η=tft
五次(6 个系数 ↔ 位置/速度/加速度 × 两端):
q(t)=qi+(qf−qi)[10η3−15η4+6η5]
关键点:能量最优解只能是三次;选五次就放弃该指标(定理 4.2.A、定义 4.2.2)。
速度梯形:起段常加速度 q¨c、中段巡航 q˙c、末段常减速度;位置 = 直线段 + 两端抛物线。
tm=2tf,qm=2qf+qi,q¨ctc=tm−tcqm−qc,qc=qi+21q¨ctc2
⟹q¨ctc2−q¨ctftc+(qf−qi)=0
tc=2tf−21q¨ctf2q¨c−4(qf−qi),∣q¨c∣≥tf24∣qf−qi∣
指定巡航速度(替代指定加速度):
tf∣qf−qi∣<∣q˙c∣≤tf2∣qf−qi∣,tc=q˙cqi−qf+q˙ctf,q¨c=qi−qf+q˙ctfq˙c2
分段位置((4.8)):
q(t)=⎩⎨⎧qi+21q¨ct2qi+q¨ctc(t−tc/2)qf−21q¨c(tf−t)20≤t≤tctc<t≤tf−tctf−tc<t≤tf
极限:下界取等 ⇒ 三角形速度曲线(无巡航段);代价:能量比最优高 12.5%(定义 4.2.3、定理 4.2.C)。
① 指定路径点速度:Πk(tk)=qk, Πk(tk+1)=qk+1, Π˙k(tk)=q˙k, Π˙k(tk+1)=q˙k+1(每段独立四方程)。
② 按规则计算速度((4.12),vk=tk−tk−1qk−qk−1):
q˙1=0,q˙k={02vk+vk+1sgnvk=sgnvk+1sgnvk=sgnvk+1,q˙N=0
③ 三次样条(加速度连续)——方程计数与降维:
4N−2 个方程4(N−2)+6 vs 未知数4(N−1) ⇒ 差 2 ⇒ 两个虚拟点 ⇒ 4(N+1)=4(N+1)
Π¨k(t)=ΔtkΠ¨k(tk+1)(tk+1−t)+ΔtkΠ¨k(tk)(t−tk)⟹AΠ=b
A 为三对角、非奇异,且只依赖 Δtk——对所有关节相同(定理 4.2.H)。
④ 抛物线 + 线性混合(经由点):
qk−1,k=Δtkqk−qk−1,q¨k=Δtk′qk,k+1−qk−1,k,bk0=qk+8(q˙k,k+1−q˙k−1,k)Δtk′
轨迹不经过路径点(削角量 = 上式第二项),故称经由点(定义 4.2.7)。
参数化与 Frenet 标架:
p=f(s),t=dsdp,n=∥d2p/ds2∥d2p/ds2,b=t×n
直线(pl→pf):p(s)=pl+∥pf−pl∥s(pf−pl),ds2d2p=0(法向无法定义)。
圆弧(轴 r、轴上点 d、圆上点 pl,δ=pl−d,退化条件 ∣δTr∣<∥δ∥):
c=d+(δTr)r,p(s)=c+Rp′(s),p′(s)=[ρcos(s/ρ), ρsin(s/ρ), 0]T
位置轨迹 = 几何 × 时间:
p˙e=s˙t,p¨e=向心:改方向ρs˙2n+切向:改大小s¨t
多段路径与时间超前:pe=p0+∑j∥pj−pj−1∥sj(pj−pj−1);Δtj=Δtj−1+δtj(定理 4.3.D)。
旋转矩阵不能逐元素插值(破坏正交性)⇒ 两条正路:
欧拉角:φe=φi+∥φf−φi∥s(φf−φi),再经 ωe=T(φe)φ˙e 换回基系。
角和轴(Rif=RiTRf):
θf=arccos(2r11+r22+r33−1),r′=2sinθf1r32−r23r13−r31r21−r12,ω′=θ˙r′
两次降维:位置落在 s(t)、姿态落在 θ(t)——都是一维标量,直接复用关节空间的方法。
能源 → 功率放大器 → 伺服发动机 → 传动装置 → 关节(定义 5.1.1)。功率恒为「力 × 速度」,由此可辨功率种类(功率总是「力 × 速度」)。
Pc (控制信号)⇒Pa (供给发动机)⇒Pm (发动机输出)⇒Pu (关节所需)
三处损失:Pda 放大器、Pdm 发动机、Pdt 传动装置(定义 5.1.2)。
关节要低速度、高转矩;伺服发动机最优工况是高速度、低转矩(定理 5.1.B)⇒ 必须加传动装置。三重作用:定量(降速增矩)/ 定性(转动→平动)/ 减重优化有效载荷(定义 5.1.3)。无间隙刚性传动 ⇒ 输入输出纯比例(定理 5.1.F)。
Va=(Ra+sLa)Ia+Vg,Vg=kvΩm
Cm=(sIm+Fm)Ωm+Cl,Cm=ktIa
kt=kv(标准国际单位下数值相等)(定义 5.2.1)。放大器 VcVa=1+sTvGv,Tv 在 10−5–10−4 s 量级。
M(s)=s(1+sTm)km
| km | Tm | 干扰增益 |
| 速度控制(ki=0) | 1/kv | RaIm/(kvkt) | Ra/(kvkt) |
| 力矩控制(Kki≫Ra) | kt/(kiFm) | Im/Fm | 1/Fm |
速度控制抗干扰更好 ⇒ 独立关节控制不加电流反馈;保护改用死区非线性的电流反馈(定理 5.2.A、定理 5.2.B)。
Q=Qm+Ql+Qc,Ql=klp,Qc=γVsp,Qm=kqΩm,Cm=(sIm+Fm)Ωm+Cl,Cm=kpp
伺服阀 Vcx=1+sTcGc(毫秒级);分布器线性化 p=kxx−kqQ(定义 5.2.3)。旋转式 kc=γV 为常数,直线式随工作点变。
kr=ωωm=rmr (典型 几十∼几百)
Cm=Ieqω˙m+Feqωm+krCl,Ieq=Im+kr2I,Feq=Fm+kr2F
惯量与摩擦按 1/kr2、负载力矩与非线性按 1/kr(定理 5.2.C)。摆的例子:Cm=Ieqω˙m+Feqωm+krmglsinkrϑm——kr 同时缩小幅值与自变量,故大减速比使非线性趋向线性化(例 5.1)。
位置环需要 积分(稳态消除干扰)+ 比例(稳定与阻尼)——纯比例给不出阻尼瞬态,原因是前向通路有两个原点极点(定理 5.2.D)。加测速反馈后:干扰点附近两个回路 ⇒ 瞬态抗扰更强、带宽更宽、振荡更少(方法 5.2.E)。推广到齿轮耦合:用 Ieq,Feq 代替 Im,Fm,干扰按 1/kr(两条推广)。
n 轨绝对式:每圈 2n 个角位置,电机侧 2π/2n;经减速比 kr 后关节侧分辨率提高 kr 倍(例 5.3.A)。12 轨 ⇒ 360∘/4096≈0.0879∘。格雷码:相邻码位只变一位,避免过渡错误的累积(定理 5.3.B)。
[fccμcc]=[RscS(rcsc)RscORsc][fsμs],S(r)=0rz−ry−rz0rxry−rx0
六维腕力传感器的 8 个张力经校准矩阵重构成 6 个分量(定义 5.4.3、定理 5.4.A)。
d=2vt,v≈331.4+0.6T (m/s, T 为摄氏温度)
相机投影:[xIyIλ]T=Ωp~c,Ω=fax000faY0X0Y01;像素 XI=xI/λ, YI=yI/λ。整体合成相机校准矩阵 Ξ=ΩΠTbc(定义 5.4.5、定理 5.4.B)。
若只背这些条,背下面这十七条(老师明确点过考试权重的在最前):
- 绿框逐列构造 Ji(转动 / 移动两式)——「肯定是要考的」;
- Aii−1 完整矩阵与「第 4 列 = 三参数肉身」;
- 右乘当前轴、左乘固定轴;
- τ=JTFe——「主要是记住这个公式就可以」;
- 解耦腕心 pw=pe−d6ae;
- ZYZ 正解矩阵与逆解三式(Atan2);
- detJ=detJ11detJ22(球腕解耦);
- R˙=S(ω)R 与 S(ω) 的三分量;
- J†=JT(JJT)−1 与零空间项 (I−J†J)q˙0;
- DH 四参数「3 常 1 变」与「能零则零」;
- 三次 / 五次的规范化形式 3η2−2η3、10η3−15η4+6η5(η=t/tf)——「约束数 = 系数个数」;
- 梯形速度曲线:tc 二次方程 (4.5)、加速度下界 ∣q¨c∣≥4∣qf−qi∣/tf2、巡航速度约束 (4.9)——「重点介绍混合梯形」;
- 两次降维:位置 p˙e=s˙t(圆弧再加向心项 s˙2/ρn)、姿态(角和轴)ω′=θ˙r′;
- 减速比两个折算因子:Ieq=Im+I/kr2、Feq=Fm+F/kr2,负载力矩 Cl/kr——「1/kr2 管惯量、1/kr 管力矩」;
- 电气驱动四方程与统一传递函数 M(s)=km/[s(1+sTm)],两种发电机的 km,Tm(速度控制 1/kv、RaIm/(kvkt);力矩控制 kt/(kiFm)、Im/Fm);
- 传动装置三重作用(定量 / 定性 / 减重解耦)与三种装置的适用场合,以及「同步带快、正时链重」;
- 传感器的三个公式:编码器 2π/2n 与减速后的分辨率、d=vt/2(声纳,含 v≈331.4+0.6T)、相机 XI=xI/λ。
配套的常见错误清单见易错点总表,术语查不到时用术语索引。