3. Simulink入门：界面、模块库、信号线、子系统封装与仿真参数设置
在完成MATLAB基础语法与脚本编程的学习后，我们正式进入Simulink环境。Simulink是MATLAB家族中基于模型设计（Model-Based Design）的核心平台，它采用图形化建模方式，通过模块（Block）与信号线（Signal Line）构建动态系统模型。对于机械臂仿真而言，Simulink的价值在于：它能够将机械臂的动力学方程、运动学解算、轨迹规划算法与电机驱动模型以模块化的方式组织起来，形成可复用、可调试、可自动生成代码的完整仿真链路。本章将围绕Simulink的界面布局、模块库体系、信号线操作、子系统封装技术以及仿真参数配置五个核心维度展开，为后续机械臂运动学与轨迹仿真搭建坚实的操作基础。
3.1 Simulink界面布局与建模环境认知
Simulink的界面设计遵循"所见即所得"的交互理念，其核心工作区域由三大部分构成：模块库浏览器（Library Browser）、模型画布（Model Canvas）以及仿真控制工具栏。模块库浏览器位于界面左侧，以树状结构组织所有可用模块，涵盖连续系统、离散系统、数学运算、信号路由、逻辑判断等基础类别，同时集成Simscape、Robotics System Toolbox等专业工具箱模块。模型画布是用户搭建模型的区域，支持鼠标拖拽、缩放、框选等操作，画布底部显示模型状态栏，可查看当前求解器类型、仿真时间与模型版本信息。
对于机械臂仿真，我们需要重点关注以下模块库分类：Simulink/Continuous（积分器、传递函数、状态空间）、Simulink/Math Operations（矩阵运算、三角函数、代数约束）、Simscape/Multibody（三维刚体、关节、约束、传感器与执行器）以及Robotics System Toolbox（机器人模型导入、逆运动学求解器、轨迹规划函数）。这些模块库共同构成了从数学建模到物理建模的完整工具链。
| 界面区域 | 核心功能 | 机械臂仿真中的典型用途 |
|---|---|---|
| 模块库浏览器 | 按类别浏览、搜索、拖拽模块 | 查找关节模块、刚体模块、轨迹生成器 |
| 模型画布 | 搭建模型、连接信号线、编辑模块参数 | 构建机械臂多体动力学模型框图 |
| 仿真工具栏 | 启动/暂停/停止仿真、设置仿真时间 | 运行轨迹跟踪仿真、实时调参 |
| 模型资源管理器 | 管理模型变量、回调函数、模型属性 | 定义机械臂DH参数、质量属性变量 |
| 诊断查看器 | 显示错误、警告、信息提示 | 排查代数环、信号维度不匹配等问题 |
3.2 模块库体系与信号线操作
Simulink模块库采用层级化组织方式，顶层为"Simulink"基础库，其下按功能划分为20余个子库。在机械臂仿真中，最常用的模块组合包括：信号源模块（Step、Sine Wave、Signal Builder、From Workspace）、连续系统模块（Integrator、Derivative、Transfer Fcn、State-Space）、数学运算模块（Gain、Sum、Product、Matrix Multiply、Trigonometric Function）以及信号路由模块（Mux、Demux、Goto、From、Data Store Memory）。对于多体动力学建模，Simscape/Multibody库提供了Revolution Joint、Prismatic Joint、Rigid Transform、Solid等物理建模模块，这些模块与Simulink标准模块通过Simulink-PS Converter和PS-Simulink Converter进行物理信号与数值信号的转换。
信号线是Simulink模型中数据流动的载体，其操作技巧直接影响建模效率。基本操作包括：连线（从输出端口拖拽至输入端口）、分支（在已有信号线上右键拖拽引出分支）、信号标签（双击信号线输入名称，便于识别）、信号维度显示（在Display菜单中开启Port Data Types和Signal Dimensions）。对于机械臂模型，信号线通常承载关节角度向量（6×1）、角速度向量（6×1）或力矩向量（6×1），建议使用Mux模块将多路标量信号合并为向量信号，或使用Bus对象组织结构化数据。
% 在MATLAB工作区定义机械臂参数
% 六轴机械臂DH参数（标准DH表）
a = [0, 0.432, 0.020, 0, 0, 0]; % 连杆长度 (m)
d = [0.150, 0, 0, 0.433, 0, 0.056]; % 连杆偏距 (m)
alpha = [pi/2, 0, pi/2, -pi/2, pi/2, 0]; % 连杆扭角 (rad)
% 在Simulink中通过From Workspace模块加载轨迹数据
t = 0:0.01:10;
q1 = 0.5*sin(0.5*t); % 关节1正弦轨迹
q2 = 0.3*cos(0.3*t); % 关节2余弦轨迹
traj_data = [t', q1', q2']; % 时间戳 + 关节位置3.3 子系统封装技术
子系统（Subsystem）是Simulink实现模型层次化、模块化设计的核心机制。当机械臂模型包含多个关节的动力学方程、复杂的逆运动学解算或轨迹规划算法时，直接平铺所有模块会导致模型混乱、难以调试。子系统封装（Subsystem Mask）则进一步提升了子系统的可用性，它允许我们为子系统创建自定义参数对话框、图标和文档，使其行为类似于一个独立的模块。
创建子系统有两种常用方法：方法一，在模型画布中框选需要封装的模块组，右键选择"Create Subsystem from Selection"（或快捷键Ctrl+G），Simulink自动将选中模块封装为一个子系统，并自动创建输入输出端口；方法二，从模块库中拖入"Subsystem"模块，双击进入子系统内部，手动搭建内部逻辑。对于机械臂仿真，建议采用方法一，因为其保留了原有的信号连接关系，减少了重新连线的错误风险。
封装（Mask）操作流程如下：右键点击子系统，选择"Mask > Create Mask"，打开封装编辑器。在Icon选项卡中，可以绘制自定义图标（如机械臂简图）；在Parameters选项卡中，定义用户可编辑的参数（如连杆长度、质量），并关联到子系统内部的模块参数；在Initialization选项卡中，编写参数初始化代码；在Documentation选项卡中，填写子系统的描述信息。封装后的子系统在机械臂仿真中具有重要价值：例如，将单关节的"电机-减速器-连杆"模型封装为一个带质量参数和减速比参数的子系统，在搭建六轴机械臂时只需拖入6个该子系统并设置不同参数即可。
背景与约束：某六轴协作机械臂，每个关节包含直流电机、谐波减速器（减速比100:1）和连杆。需要在Simulink中搭建完整的动力学仿真模型，要求模型可配置不同关节参数。
初始现象：直接平铺所有模块导致模型包含超过200个模块，信号线交叉严重，调试困难。
关键参数：电机电枢电阻R=2.6Ω，电感L=0.18mH，转矩常数Kt=0.12Nm/A，反电动势常数Ke=0.12V/(rad/s)，减速比N=100，连杆质量m=3.5kg，质心位置r=0.15m。
操作步骤：
1. 搭建单关节模型：电源→PWM调制→电机电枢回路→电磁转矩→机械转动方程→减速器→连杆负载；
2. 框选所有模块，Ctrl+G创建子系统；
3. 右键子系统→Mask→Create Mask，添加参数：R、L、Kt、Ke、N、m、r；
4. 在Initialization中编写：
J_load = m*r^2; B_eff = B_motor + B_load/N^2;5. 将子系统内部对应的Gain模块参数改为
1/(R*J_eff)等表达式；6. 保存封装，在模型画布中复制5份，分别设置不同参数。
结果对比：封装后模型顶层仅包含6个子系统+轨迹生成器+示波器，信号线清晰，参数修改只需双击子系统在对话框中输入数值。
风险与代价：封装后内部信号无法直接观测，需在子系统内部添加To Workspace模块或使用Simulink数据探查器；参数关联错误时排查难度增加，建议在Initialization中添加参数范围检查代码。
3.4 仿真参数设置与求解器选择
仿真参数设置是决定机械臂仿真精度与效率的关键环节。在Simulink工具栏点击"Model Settings"（或快捷键Ctrl+E）打开配置对话框，核心设置项包括：仿真时间（Start Time与Stop Time）、求解器类型（变步长/固定步长）、求解器算法（ode45、ode15s、ode23t等）、步长参数（最大步长、初始步长）、误差容限（相对误差、误差）以及代数环诊断设置。
对于机械臂动力学仿真，求解器选择需要特别谨慎。机械臂动力学方程通常表现为强耦合、非线性的二阶常微分方程组，且可能包含刚性（Stiff）特性——例如当关节减速比很大时，电机高速旋转与连杆低速运动形成时间常数差异悬殊的刚性系统。此时，非刚性求解器（如ode45）可能计算效率极低甚至无法收敛，应改用刚性求解器（如ode15s、ode23t）。下表给出了常用求解器的适用场景对比：
| 求解器 | 类型 | 适用场景 | 机械臂仿真建议 |
|---|---|---|---|
| ode45 | 变步长，非刚性 | 大多数连续系统，精度适中 | 仅用于忽略电机动态的简化运动学仿真 |
| ode15s | 变步长，刚性 | 含快速/慢速动态的系统 | 推荐用于含减速器、电机电枢动力学的完整模型 |
| ode23t | 变步长，适度刚性 | 中等刚性，要求较低数值阻尼 | 适用于关节柔性建模 |
| ode4（RK4） | 固定步长 | 实时仿真、硬件在环 | 需配合信号采样时间设置，步长取1ms~10ms |
在机械臂轨迹仿真中，还需要注意代数环问题。当模型中存在信号直接通过无延迟模块形成闭环时（例如关节角速度同时用于计算动力学和反馈控制），Simulink会在每个仿真步长内迭代求解代数方程，降低仿真速度甚至导致仿真失败。解决方法包括：在反馈回路中插入Memory模块或Unit Delay模块打破代数环，或者使用Simulink的代数环诊断功能定位问题。
0.01s，观察仿真结果后再逐步调整。同时，误差容限不宜设置过小（如1e-9），否则会显著增加计算时间，对于机械臂仿真，1e-5到1e-6通常已足够。
3.5 机械臂仿真模型搭建实战流程
综合以上知识，我们给出一个标准的机械臂仿真模型搭建流程，该流程可作为后续章节中运动学与轨迹仿真实验的操作模板：
- 参数准备：在MATLAB工作区定义机械臂DH参数、质量参数、电机参数，保存为结构体或脚本文件。
- 模型框架搭建：新建Simulink模型，从模块库拖入"From Workspace"（加载轨迹数据）、"Demux"（分解关节向量）、"Subsystem"（运动学解算）、"Subsystem"（动力学模型）、"Scope"（结果可视化）。
- 子系统内部实现：双击进入运动学子系统，使用Trigonometric Function、Matrix Multiply等模块实现正运动学与逆运动学解算；进入动力学子系统，使用Integrator、Gain、Sum等模块实现关节动力学方程。
- 信号连接与检查：连接所有信号线，使用信号标签标注各端口含义；点击"Display > Signals & Ports > Port Data Types"检查信号维度是否匹配。
- 仿真参数配置：设置仿真时间（如10s），选择求解器（推荐ode15s），设置最大步长0.01s，相对误差1e-5。
- 运行与调试：点击"Run"按钮，观察Scope输出；若出现错误，查看诊断查看器中的提示信息，常见错误包括信号维度不匹配、代数环、参数未定义等。
- 结果后处理：使用To Workspace模块将仿真数据导出至MATLAB工作区，利用plot函数绘制关节角度跟踪曲线、力矩曲线，计算跟踪误差。