汽车二自由度动力学仿真与状态空间分析源码

资 源 简 介

本项目旨在构建一个经典的汽车二自由度（Lateral and Yaw dynamics）线性动力学模型，专为车辆工程及控制理论初学者设计，用于分析汽车在转向过程中的瞬态响应和稳态特性。模型基于“自行车模型”简化理论，假设车辆纵向速度恒定，忽略悬架侧倾和俯仰运动，仅考虑并在数学上描述车辆的侧向运动（沿Y轴）和横摆运动（绕Z轴）。功能包括：1. 参数化建模，允许用户定义车辆质量、绕Z轴转动惯量、前后轴距、前后轮胎侧偏刚度等物理参数；2. 状态空间方程构建，将二阶微分方程转化为标准状态空间形式 (dx/dt = Ax + Bu)，便于MATLAB求解；3. 仿真求解模块，利用MATLAB的ode45求解器或Simulink模块搭建仿真框图，计算在阶跃输入、正弦输入等不同前轮转角工况下的系统响应；4. 数据可视化，自动绘制横摆角速度、质心侧偏角随时间变化的响应曲线，以及侧向加速度图表。该项目能够帮助学习者深入理解车辆操纵稳定性、欠转向与过转向特性，并为后续学习四自由度或七自由度模型打下坚实基础。

详 情 说 明

汽车二自由度动力学仿真模型系统

项目简介

本项目构建了一个经典的汽车二自由度（Degrees of Freedom, DOF）线性动力学仿真模型。该系统专为车辆工程及控制理论的学习者设计，基于“自行车模型”简化理论，忽略了悬架侧倾、俯仰及纵向速度变化，聚焦于分析汽车在恒定车速下的侧向运动（Lateral Dynamics）和横摆运动（Yaw Dynamics）。通过数学建模与数值求解，本系统能够直观展示车辆对转向输入的瞬态响应与稳态特性，是理解车辆操纵稳定性（Handling Stability）的基础工具。

功能特性

• 参数化车辆建模：支持自定义车辆物理参数，包括整车质量、绕Z轴转动惯量、前后轴距分布以及前后轮胎的侧偏刚度。
• 状态空间方程自动构建：根据输入的物理参数，自动推导并生成线性系统的状态空间矩阵（系统矩阵 A 和 输入矩阵 B）。
• 系统稳定性判定：通过计算系统矩阵 A 的特征值，自动判断当前车速下车辆系统的稳定性。
• 多种仿真工况：

* 角阶跃输入（Step Input）：模拟驾驶员快速转动方向盘并保持不动的情景，分析瞬态响应时间及稳态值。 * 正弦输入（Sine Input）：模拟蛇形机动或避障测试，分析车辆在连续转向下的动态跟随能力。

• 稳态特性深度分析：基于动力学参数计算稳定性因数（Stability Factor），自动判别车辆属于不足转向、过度转向还是中性转向，并计算相应的特征车速或临界车速。
• 多维度数据可视化：自动生成包含前轮转角输入、横摆角速度响应、质心侧偏角响应及侧向加速度的对比图表。

系统要求

• MATLAB R2016b 或更高版本
• 无需额外工具箱（使用标准 MATLAB 数值计算与绘图函数）

使用方法

1. 确保 MATLAB 工作路径包含本项目的 main.m 文件。
2. 直接运行 main 函数。
3. 程序将自动在控制台输出车辆参数、稳定性检查结果及稳态特性分析数据。
4. 程序运行结束后，将弹出一个包含六个子图的窗口，分别展示阶跃响应和正弦响应下的各项动力学指标。

代码实现逻辑详解

本项目核心逻辑封装在 main.m 文件中，通过主函数与多个子函数协同工作，具体实现流程如下：

1. 初始化与参数定义

程序首先清理工作环境，然后定义了一个名为 vehicle 的结构体。该结构体封装了车辆的核心物理参数（质量、惯量、轴距等）以及仿真工况参数（如纵向车速）。特别注意，代码中采用了轮胎侧偏刚度为负值的符号惯例（恢复力方向与侧偏角相反）。

2. 状态空间矩阵构建

通过 build_state_space 子函数，系统将二自由度微分方程转化为标准的 $dot{x} = Ax + Bu$ 形式：

• 状态变量 $x$：选取质心侧偏角（beta）和横摆角速度（r）。
• 输入变量 $u$：选取前轮转角（delta）。
• 矩阵生成：利用牛顿第二定律对车辆侧向受力和绕Z轴力矩进行平衡推导，计算出矩阵 $A$（系统动态矩阵）和 $B$（输入控制矩阵）的各个元素。

3. 稳定性检查

主函数调用 MATLAB 内置的 eig 函数计算矩阵 $A$ 的特征值。若所有特征值的实部均小于 0，则判定系统在当前车速下是稳定的。

4. 时域仿真求解

程序设定了两种标准工况，利用 ode45（四五阶龙格-库塔法）数值求解器进行解算：

• 阶跃响应工况：定义了一个延迟 1 秒触发、幅值为 0.1 弧度的阶跃函数。
• 正弦响应工况：定义了一个频率为 0.5 Hz、幅值为 0.05 弧度的正弦函数。
• 求解核心：vehicle_dynamics 子函数作为 ODE 求解器的回调函数，实时计算 $Ax + Bu$ 的导数。

5. 输出量后处理

由于状态变量仅包含侧偏角和横摆角速度，calculate_outputs 子函数负责根据求解出的状态序列，进一步计算侧向加速度。该计算并非简单对速度求导，而是基于物理公式 $F_y/m$（即前后轮胎侧向力之和除以质量）得出，更能反映真实的动力学受力情况。同时，该函数还计算了前后轮的滑移角。

6. 稳态特性分析

analyze_steady_state 子函数基于代数方程对车辆固有的操纵特性进行分析：

• 稳定性因数 K：根据质量、轴距分配及轮胎刚度计算 $K$ 值。
• 转向特性判定：

* 若 $K > 0$，判定为不足转向，并计算特征车速。 * 若 $K < 0$，判定为过度转向，并计算临界车速（失稳边界）。

• 稳态增益：计算横摆角速度相对于前轮转角的稳态增益比。

7. 数据可视化

visualize_results 子函数创建一个包含 3行2列 子图的图形窗口：

• 左侧列展示阶跃响应结果：输入转角、车辆状态（横摆角速度和侧偏角双轴显示）、侧向加速度。
• 右侧列展示正弦响应结果：对应的输入、状态及加速度曲线。
• 所有角度数据在绘图时均自动从弧度转换为角度（deg），以便于直观阅读。

关键算法与数学模型

• 二自由度自行车模型：假设车辆左右对称，将双轨车辆简化为单轨模型，仅考虑侧向运动和横摆运动两个自由度。
• 线性轮胎模型：假设轮胎侧偏力与侧偏角成线性关系（$F_y = k cdot alpha$），适用于此时域仿真中的小角度假设范围（< 5度）。
• 状态空间法：利用现代控制理论中的状态空间表达式处理多输入多输出（虽此处为单输入）动态系统，便于后续扩展控制算法（如 LQR 控制等）。