【STM32】单级与串级PID控制的C语言实现

ops/2024/11/13 9:54:17/

【STM32】单级与串级PID的C语言实现

  • 前言
  • PID理论
    • 什么是PID
    • PID计算过程
    • PID计算公式
    • Pout、Iout、Dout的作用
    • 单级PID与串级PID
  • PID应用
    • 单级PID
    • 串级PID

前言

笔者最近在学习PID控制器,本文基于Blog做以总结。CSDN上已有大量PID理论知识的优秀文章,因此本文将略写理论部分,重点放在应用。

PID理论

什么是PID

PID(比例-积分-微分)控制器是一种广泛用于自动控制系统的闭环反馈控制器,能够实现系统的稳定、快速和精确控制。

PID计算过程

单级PID
在这里插入图片描述
串级PID
串级PID分为外环和内环,是两个单级PID嵌套构成的

在这里插入图片描述

PID计算公式

连续型
在这里插入图片描述
离散型
在这里插入图片描述

Pout、Iout、Dout的作用

  • Pout:主要负责控制,使反馈量快速接近目标值,但可能引起振荡
  • Iout:消除稳态误差,但会增加超调
  • Dout:提供阻尼,抑制振荡和超调,但可能降低响应速度

单级PID与串级PID

我们来思考一下:为什么自动控制系统中要引入PID?能否直接控制目标值?先来看个常见的PID应用例子

例子一:水阀控制水量
在这个例子中,目标值是水槽的水量(也即水面高度,量纲为长度,单位为m)。而可直接控制的变量是水阀的开度(即水流速度,量纲为速度,单位为m/s)。
注意,控制量相较于目标值为低一阶物理量,故我们无法直接控制高阶物理量(水面高度),只能控制低阶物理量(水流速度)。
(此时因为控制量与目标值只差一阶,所以单级PID即可。假设我们不能直接控制水流速度,只能直接控制水流的加速度,那么就需要串级PID控制了)
因此引入PID控制器,为了更好的通过控制低阶物理量完成对高阶物理量的控制
ps:水面高度和水体积呈线性关系,水阀开度和水流速也是线性相关,因此通过合理调节PID参数,系统可以自适应这种内含的线性关系

例子二:循迹小车

在这个例子中,目标是控制小车位置,目标值是小车目标位置,而可直接控制的变量是电机电流(可以理解为加速度或力)。由于电流(加速度)与位置之间相差两阶,因此需要使用串级PID控制

具体做法:

  • 外环PID:根据目标位置和当前实际位置的误差,输出一个目标速度给内环PID
  • 内环PID:根据目标速度和当前实际速度的误差,输出控制量(如PWM占空比)来给执行器

这种串级PID控制方法,通过外环控制速度、调节位置,内环控制PWM、调节速度,从而实现了精准的循迹控制

PID应用

单级PID

pid.h

#ifndef __PID_H__
#define __PID_H__/* USER CODE BEGIN Includes *//* USER CODE END Includes *///首先定义PID结构体用于存放一个PID的数据
typedef struct
{float kp, ki, kd; //三个系数float error, lastError; //误差、上次误差float integral, maxIntegral; //积分、积分限幅float output, maxOutput; //输出、输出限幅
}PID;void PID_Init(PID *pid, float p, float i, float d, float maxI, float maxOut);
void PID_Calc(PID *pid, float reference, float feedback);#endif

pid.c

#include "pid.h"PID speed_pid_speed = {0}; //定义并初始化一个PID结构体变量,单级速度环//初始化pid参数的函数
void PID_Init(PID *pid, float p, float i, float d, float maxI, float maxOut)
{pid->kp = p;pid->ki = i;pid->kd = d;pid->maxIntegral = maxI;pid->maxOutput = maxOut;
}//进行一次PID控制器的计算,更新控制量:pid->output
//参数为(pid结构体,目标值,反馈值),计算结果放在pid结构体的output成员中
void PID_Calc(PID *pid, float reference, float feedback)
{//更新数据pid->lastError = pid->error; //将旧error存起来pid->error = reference - feedback; //计算新error//计算比例float pout = pid->error * pid->kp;//计算积分pid->integral += pid->error * pid->ki;//计算微分float dout = (pid->error - pid->lastError) * pid->kd;//积分限幅if(pid->integral > pid->maxIntegral) pid->integral = pid->maxIntegral;else if(pid->integral < -pid->maxIntegral) pid->integral = -pid->maxIntegral;//计算输出pid->output = pout+dout + pid->integral;//输出限幅if(pid->output > pid->maxOutput) pid->output =  pid->maxOutput;else if(pid->output < -pid->maxOutput) pid->output = -pid->maxOutput;
}

my_main.c

void setup(void)
{PID_Init(&speed_pid,1,0,0,0,1000);	//举例
}void loop(void)
{//更新被控对象反馈值和目标值(反正要不断更新、获取PID_Calc(&speed_pid,targetValue,feedbackValue);设置执行器输出大小(speed.output);//在控制电机转速的单级速度环pid中,控制执行器的操作为set_pwm(speed.output);HAL_Delay(10);	//控制采样周期,防止频繁计算导致系统抖动与资源消耗过高(根据不同情况合理配置即可
}

串级PID

pid.h

#ifndef __PID_H__
#define __PID_H__/* USER CODE BEGIN Includes *//* USER CODE END Includes *///首先定义单级PID结构体
typedef struct
{float kp, ki, kd; //三个系数float error, lastError; //误差、上次误差float integral, maxIntegral; //积分、积分限幅float output, maxOutput; //输出、输出限幅
}PID;//再定义串级PID的结构体,嵌套两个单级PID
typedef struct
{PID inner; //内环PID outer; //外环float output; //串级输出,等于inner.output
}CascadePID;void PID_Init(PID *pid, float p, float i, float d, float maxI, float maxOut);
void PID_Calc(PID *pid, float reference, float feedback);
void PID_CascadeCalc(CascadePID *pid, float outerRef, float outerFdb, float innerFdb);#endif

pid.c

//初始化pid参数的函数
void PID_Init(PID *pid, float p, float i, float d, float maxI, float maxOut)
{pid->kp = p;pid->ki = i;pid->kd = d;pid->maxIntegral = maxI;pid->maxOutput = maxOut;
}//进行一次PID控制器的计算,更新控制量:pid->output
//参数为(pid结构体,目标值,反馈值),计算结果放在pid结构体的output成员中
void PID_Calc(PID *pid, float reference, float feedback)
{//更新数据pid->lastError = pid->error; //将旧error存起来pid->error = reference - feedback; //计算新error//计算比例float pout = pid->error * pid->kp;//计算积分pid->integral += pid->error * pid->ki;//计算微分float dout = (pid->error - pid->lastError) * pid->kd;//积分限幅if(pid->integral > pid->maxIntegral) pid->integral = pid->maxIntegral;else if(pid->integral < -pid->maxIntegral) pid->integral = -pid->maxIntegral;//计算输出pid->output = pout+dout + pid->integral;//输出限幅if(pid->output > pid->maxOutput) pid->output =  pid->maxOutput;else if(pid->output < -pid->maxOutput) pid->output = -pid->maxOutput;
}//串级PID的计算函数(调用两次单级PID计算,先算外环再算内环,串级cascade_pid->output被更新
//参数(串级PID结构体,外环目标值,外环反馈值,内环反馈值)
void PID_CascadeCalc(CascadePID *pid, float outerRef, float outerFdb, float innerFdb)
{PID_Calc(&pid->outer, outerRef, outerFdb); //先计算外环PID_Calc(&pid->inner, pid->outer.output, innerFdb); //再计算内环pid->output = pid->inner.output; //控制量 = 内环输出 = 串级PID的输出(所以
}

my_main.c
内环 PID 控制频率应高于外环,以保证内环能快速响应动态变化。建议避免内外环在同一周期更新,最好通过定时中断分别计算单级和串级 PID 控制,确保系统稳定性和响应速度

void setup(void)
{CascadePID cascade_pid = {0}; //创建串级PID结构体变量PID_Init(&cascade_pid.inner, 10, 0, 0, 0, 1000); //初始化内环参数PID_Init(&cascade_pid.outer, 5, 0, 5, 0, 100); //初始化外环参数	
}void loop(void)
{//更新、获取外环目标值、外环反馈值、内环反馈值PID_CascadeCalc(&cascade_pid, outerTarget, outerFeedback, innerFeedback); //进行串级PID计算设定执行机构输出大小(cascade_pid.output);	//在循迹控制中,控制执行输出为set_pwm(cascade_pid.output)HAL_Delay(10);	//控制采样周期,防止频繁计算导致系统抖动与资源消耗过高(根据不同情况合理配置即可    
}

http://www.ppmy.cn/ops/111703.html

相关文章

【Kubernetes】(K8S)彻底卸载详细教程

以下全部操作都是使用root用户进行&#xff08;非root用户可以使用sudo&#xff09;&#xff0c;并且全部命令都需要在Kubernetes集群的所有节点分别执行&#xff1a; 第一步、停止K8S 所有节点执行&#xff1a; 1 2 3 systemctl stop kubelet systemctl stop etcd systemct…

[Web安全 网络安全]-XSS跨脚本攻击

文章目录&#xff1a; 一&#xff1a;前言 1.定义 2.漏洞出现的原因 3.鉴别可能存在XSS漏洞的地方 4.攻击原理 5.危害 6.防御 7.环境 7.1 靶场 7.2 自动扫描工具 7.3 手工测试工具 8.payload是什么 二&#xff1a;常用的标签语法 三&#xff1a;XSS的分类 反射…

ubuntu系统服务器离线安装python包

一、根据工程需要本地下载所需python包 1. 下载环境requirements.txt pip freeze > requirements.txt2. 根据requirements.txt下载python包 注意&#xff1a;查看服务器属于x_86架构还是arm架构、cpython还是pypy 2.1 确定服务器架构&#xff08;终端输入&#xff09; …

SpringBoot用kafka.listener监听接受Kafka消息

1.创建kafka监听配置并进行注册 import org.apache.kafka.clients.consumer.ConsumerConfig; import org.apache.kafka.common.serialization.StringDeserializer; import org.springframework.beans.factory.annotation.Value; import org.springframework.context.annotation…

连续时间,离散频率 傅里叶

时域周期——不是把一个信号周期化&#xff0c;而是周期信号取一个周期是x(t),对其周期化不会发生时域的重叠。故当接收到信号&#xff0c;在DFT时&#xff0c;以整个接收到的时间信号为周期进行延拓 推导公式时思路&#xff1a;时域卷积周期冲击&#xff0c;用傅里叶变换推导出…

react 基础语法

前置知识 类的回顾 通过class关键字定义一个类 类名首字母大写 class类有constructor构造器 new 一个类得到一个实例 类还有方法&#xff0c;该方法也会在其原型上 static静态数据&#xff0c;访问静态属性通过 类名.id getter和setter getter&#xff1a;定义一个属性&…

Android 11(API 级别 30)及以上版本中,将Bitmap保存到设备上

调用 saveBitmapToMediaStore(getContentResolver(),bitmap,“图片名”,mimeType); 参数解析&#xff1a; Bitmap myBitmap ...; // 这里应该是你获取或创建Bitmap的代码 private String mimeType "image/jpeg"; // 或者"image/png"&#xff0c;取决于…

CUDA及GPU学习资源汇总

CUDA C Programming Guide 的中文翻译版GPU中的SM和warp的关系推荐几个不错的CUDA入门教程