控制系统仿真：从水箱水位到倒立摆的探索-洪萨配资

MATLAB仿真simulink水箱水位控制器设置pid 源码＋报告倒立摆控制器设置源码＋报告

在控制系统的学习与实践中，MATLAB 的 Simulink 是一款强大的工具，它能帮助我们直观地搭建模型并进行仿真分析。今天咱们就来聊聊水箱水位控制器和倒立摆控制器的设置，还会附上源码和相关报告要点。

水箱水位控制器之 PID 设置

理论基础

PID 控制器，即比例（Proportional）、积分（Integral）、微分（Derivative）控制器，通过对误差的比例、积分和微分运算，来调整控制量，使系统输出尽可能接近设定值。比例环节能快速响应误差，积分环节消除稳态误差，微分环节预测误差变化趋势，提前进行调整。

Simulink 模型搭建

搭建基本模型：在 Simulink 库中找到所需模块，比如水箱模型（可以用一个简单的蓄水池模块代表），输入为水流速度，输出为水位高度。再添加一个 PID Controller 模块用于控制。
连接模块：将输入信号连接到 PID 控制器的输入端，PID 控制器的输出连接到水箱的水流输入端口，水箱的水位输出连接到反馈端口，与设定水位进行比较，得出误差信号反馈给 PID 控制器。

代码实现（MATLAB 脚本形式）

% 定义PID参数 Kp = 1; Ki = 0.1; Kd = 0.01; % 初始化变量 time = 0:0.01:10; % 仿真时间范围 setpoint = ones(size(time)) * 5; % 设定水位为5 process_variable = zeros(size(time)); error = zeros(size(time)); integral = 0; derivative = 0; previous_error = 0; for i = 2:length(time) error(i) = setpoint(i) - process_variable(i - 1); integral = integral + error(i) * 0.01; derivative = (error(i) - previous_error) / 0.01; control_signal = Kp * error(i) + Ki * integral + Kd * derivative; % 简单模拟水箱动态，这里假设水位变化与控制信号成正比 process_variable(i) = process_variable(i - 1) + control_signal * 0.01; previous_error = error(i); end figure; plot(time, setpoint, 'r', 'LineWidth', 1.5); hold on; plot(time, process_variable, 'b', 'LineWidth', 1.5); legend('设定水位', '实际水位'); xlabel('时间(s)'); ylabel('水位(m)'); title('水箱水位PID控制仿真');

代码分析

首先定义了Kp（比例系数）、Ki（积分系数）、Kd（微分系数），这些值会显著影响控制效果。
初始化了时间范围time，设定水位setpoint以及用于存储过程变量process_variable、误差error等的数组。
在for循环中，计算每个时间步的误差，积分项通过不断累加误差与时间步长的乘积得到，微分项根据当前误差与上一步误差的差值除以时间步长算出。
然后根据 PID 公式计算控制信号control_signal，并简单模拟水箱水位随控制信号的变化。
最后通过绘图函数直观展示设定水位和实际水位的变化情况。

报告要点

模型描述：详细阐述水箱模型结构，包括输入输出关系，以及 PID 控制器在系统中的作用。
参数调整：记录不同 PID 参数下的仿真结果，分析参数变化对水位控制效果（如响应速度、稳态误差等）的影响。
结论：总结通过 PID 控制实现水箱水位稳定的可行性，以及可能存在的改进方向。

倒立摆控制器设置

理论基础

倒立摆是一个典型的不稳定系统，控制目标是通过对摆杆底部施加力或扭矩，使摆杆保持垂直状态。常用的控制方法有 LQR（线性二次型调节器）、PID 等。这里我们以 PID 控制为例。

Simulink 模型搭建

模型组件：在 Simulink 中构建倒立摆模型，通常包括摆杆模型（可通过刚体动力学模块实现），一个力或扭矩输入模块，以及角度输出反馈模块。
连接与反馈：将角度反馈信号与期望角度（通常为 0 表示垂直）比较得出误差，输入到 PID 控制器，PID 控制器输出的控制信号连接到摆杆底部的力或扭矩输入端口。

代码实现（MATLAB 脚本形式）

% 倒立摆参数 m = 0.1; % 摆杆质量 L = 0.5; % 摆杆长度 g = 9.8; % 重力加速度 % 定义PID参数 Kp = 10; Ki = 0.1; Kd = 1; % 初始化变量 time = 0:0.01:10; theta_setpoint = zeros(size(time)); % 期望角度为0 theta = zeros(size(time)); omega = zeros(size(time)); % 角速度 error = zeros(size(time)); integral = 0; derivative = 0; previous_error = 0; for i = 2:length(time) error(i) = theta_setpoint(i) - theta(i - 1); integral = integral + error(i) * 0.01; derivative = (error(i) - previous_error) / 0.01; control_torque = Kp * error(i) + Ki * integral + Kd * derivative; % 简单的倒立摆动力学模型（Euler法近似） alpha = (control_torque - m * g * L * sin(theta(i - 1))) / (m * L ^ 2); omega(i) = omega(i - 1) + alpha * 0.01; theta(i) = theta(i - 1) + omega(i) * 0.01; previous_error = error(i); end figure; plot(time, theta_setpoint, 'r', 'LineWidth', 1.5); hold on; plot(time, theta * 180 / pi, 'b', 'LineWidth', 1.5); legend('期望角度', '实际角度(°)'); xlabel('时间(s)'); ylabel('角度(°)'); title('倒立摆PID控制仿真');