使用CCS软件查看PID曲线

server/2024/10/18 7:52:46/

  在刚开始学习PID的时候,都需要借助PID曲线来理解比例,积分,微分这三个参数的具体作用。但是这些曲线生成一般都需要借助上位机软件或者在网页上才能实现。如果是在单片机上调试程序的话,想要看曲线,一般就是通过串口将数据发送出去,然后使用串口波形显示软化将曲线打印出来。
如果使用的是CCS软件来开发的话,那么就可以直接使用CCS软件自带的图形显示工具来打印曲线
  比如要测试的代码如下:

typedef struct
{float P;float I;float D;float limit;
}PID;typedef struct
{float Current_Error;			//当前误差float Last_Error;				//上一次误差float Previous_Error;			//上上次误差
}ERROR;
// 位置式PID
float PID_Realize(ERROR *err, PID *pid, float feedback, float reference)
{float iError, Realize;    						// 当前误差 实际输出iError = reference - feedback;					// 计算当前误差  期望值-当前值err->Current_Error += pid->I * iError;    		// 误差积分err->Current_Error = err->Current_Error > pid->limit ? pid->limit : err->Current_Error;	//积分限幅err->Current_Error = err->Current_Error < -pid->limit ? -pid->limit : err->Current_Error;Realize = pid->P * iError       				//比例P+ err->Current_Error    						//积分I+ pid->D * (iError - err->Last_Error);  //微分Derr->Last_Error = iError;		  				// 更新上次误差return Realize;    								// 返回实际值
}// 增量式PID
float PID_Increase(ERROR *err, PID *pid, float feedback, float reference)
{float iError, Increase;    						// 当前误差  最后得出的实际增量iError = reference - feedback;    				// 计算当前误差Increase = pid->P * (iError - err->Last_Error)  //比例P+ pid->I * iError      							//积分I+ pid->D * (iError - 2 * err->Last_Error + err->Previous_Error);    //微分Derr->Previous_Error = err->Last_Error;    		// 更新前次误差err->Last_Error = iError;		  				// 更新上次误差return Increase;								// 返回增量
}

  这是两种很典型的PID算法,一个位置式,一个增量式。下面就在CCS软件中测试这两个函数,这里使用的是dsp28335芯片。首先定义PID和ERROR结构体变量。

PID pid = { 0.1, 0.01, 0.001, 100.0 };
ERROR err = { 0.0, 0.0, 0.0 };

  接下来就可以直接在main函数中调用上面两种pid的算法了。

float ref = 100;		//参考值
float sample1=0,sample2=0;	//采样值
void main()
{InitSysCtrl();InitPieCtrl();IER = 0x0000;IFR = 0x0000;InitPieVectTable();LED_Init();while(1){sample1 += PID_Realize(&err,&pid,sample1,ref);sample2 += PID_Increase(&err,&pid,sample2,ref);}
}

  将需要设定的参考值ref设置为100,然后两个采样值sample1sample2都设置为0,将PID计算出来的结果作为第二次的输入,然后在调用PID计算,这样一直循环。

  代码编写完成之后,直接debug开始调试。

进入调试界面之后,在变量sample1上单击鼠标右键,选择添加到观察窗口
在这里插入图片描述
弹出的窗口上点OK
在这里插入图片描述
此时窗口右边就会出现一个变量显示区域
在这里插入图片描述
同样在变量sample2上也单击右键,添加到观察窗口。
在这里插入图片描述
  添加完成之后,将黄色箭头的图标选中,这个是实时刷新功能,这样在代码全速运行的时候,就可实时看到变量值在变化。接下来全速运行代码。
在这里插入图片描述
  此时可以看到这两个变量值都停在在100作用不变化了,没有看到变换的过程。下面将这两个值添加到图形显示中。在观测窗口变量值上单击鼠标右键,选择Graph。

在这里插入图片描述
此时就会新增一个图形显示窗口
在这里插入图片描述
在另一个变量上右键,选择Graph。
在这里插入图片描述
  此时又增加了一个显示窗口,但是这两个窗口重叠在了一起,看起来不方便,在窗口标题处按住鼠标左键,拖动窗口,将两个图形显示调节到合适的位置。
在这里插入图片描述

  然后再这两个图形显示窗口上,将实时刷新的按钮选中。
在这里插入图片描述
接下里全速运行程序。
在这里插入图片描述
  这时候PID调节的曲线就显示出来了,但是这个曲线PID调节的曲线看起来不太一样,难道是程序有问题吗?当然不是程序的问题,是这个图形显示的刷新率太低了,而我们的程序执行的太快了。下面在代码中添加延时。

while(1){sample1 += PID_Realize(&err,&pid,sample1,ref);sample2 += PID_Increase(&err,&pid,sample2,ref);DELAY_US(1000*200);}

  重新编译运行代码
在这里插入图片描述
在这里插入图片描述
这时候就可看到右侧的变量和曲线都在缓慢的变化,可以看到PID的调节过程了。
在这里插入图片描述
  为了方便修改PID的参数,将结构体变量pid也添加到观察窗口中。在代码中选中pid然后右键添加到观察窗口。
在这里插入图片描述
  此时右边的观察窗口中就多了一个pid的结构体变量,将PID前面的加号点开,这时候就可以看到结构体中每个成员的值了。
在这里插入图片描述
  这时候在变量窗口中就可以直接修改这些参数的值了。鼠标双击要修改的值,选中值之后,直接输入要改的值,然后按回车键。
在这里插入图片描述
这里将P的值修改为0.2在这里插入图片描述
修改完成之后,底下的曲线没有发生变化。要想重新看pid的调节过程,这里将两个采样值都改为0 。
在这里插入图片描述
这时候PID函数就会从0开始重新计算,底下的曲线也会重新绘制。
在这里插入图片描述
  这样只需要在变量观察窗口直接修改PID的参数,就能实时看到PID曲线的效果了。不用在程序里面修改一次参数然后又重新仿真运行。

  还有一种方法,不需要加延时也能实时查看波形。下面就介绍一下这个方法。
定义两个数组来存放采样值,然后增加一个变量来记录数组的下标。在主函数中将采样的两个值分别存放到数组中,当数组个数超过200个之后,将下标清零。此时代码中将延时函数去掉了。没有延时,相当于代码按照正常的计算速度运行。

int cnt=0;
float value1[200];
float value2[200];while(1){sample1 += PID_Realize(&err,&pid,sample1,ref);sample2 += PID_Increase(&err,&pid,sample2,ref);//DELAY_US(1000*200);if(cnt<200){value1[cnt] = sample1;value2[cnt] = sample2;cnt++;}elsecnt =0;}

仿真运行代码,进入仿真界面后再工具栏选择 Tools—Grapg—Single Time
在这里插入图片描述
这时候会弹出一个对话框,在对话框中填入需要观测的数组信息。
在这里插入图片描述
将数组相关信息填写进去,然后点OK按钮,此时就会新增一个图形显示窗口
在这里插入图片描述
按照同样的方法,在工具栏选择 Tools—Grapg—Single Time,将value2数组也添加到图形显示中,然后选中实时刷新按钮。
在这里插入图片描述
然后在 cnt=0这一行打个断点,用来观察一次计算完成之后的波形。全速运行代码。
在这里插入图片描述
  程序在断点处停止后,左边两个图形显示的是正常的PID调节曲线,右边两个是阶跃信号。这是因为取消延时后,图形显示的刷新率太低,不能实时捕获到每个数据,所以直接显示变量的曲线就看不出来细节了,而数组显示比较完整,是因为代码中将过程中的所有数据都先存储到了数组中,然后根据数组中的数据来绘制曲线
  继续全速运行,这时候曲线就会变成一条直线,这是因为PID调节完成之后,数组中的值被更新了。
在这里插入图片描述

如果还想继续观察调节过程,在程序中断之后,将sample1sample2的值改为0,修改完之后,这时变量的曲线值也会从100变为0.
在这里插入图片描述
  全速运行代码
在这里插入图片描述
这时候左边两个数组的曲线就会被重新绘制出来。

  如果要修改PI参数的话,每次程序在断点处停止之后,修改PI参数,然后将sample1sample2的值改为0,接着全速运行,就可观测到新的PI参数生成的曲线了。
在这里插入图片描述
  这样在学习PID的时候就不用借助上位机软件,使用CCS本身的图形显示工具,就能看到不同参数下PID调节的差别了。
在这里插入图片描述
完整代码如下

#include "DSP2833x_Device.h"     // DSP2833x Headerfile Include File
#include "DSP2833x_Examples.h"   // DSP2833x Examples Include File
#include "leds.h"
typedef struct
{float P;float I;float D;float limit;
} PID;typedef struct
{float Current_Error;			//当前误差float Last_Error;				//上一次误差float Previous_Error;			//上上次误差
} ERROR;// 位置式PID控制
float PID_Realize(ERROR *err, PID *pid, float feedback, float reference)
{float iError, Realize;    						// 当前误差 实际输出iError = reference - feedback;					// 计算当前误差  期望值-当前值err->Current_Error += pid->I * iError;    		// 误差积分err->Current_Error = err->Current_Error > pid->limit ? pid->limit : err->Current_Error;    //积分限幅err->Current_Error = err->Current_Error < -pid->limit ? -pid->limit : err->Current_Error;Realize = pid->P * iError       				//比例P+ err->Current_Error    						//积分I+ pid->D * (iError - err->Last_Error);    //微分Derr->Last_Error = iError;		  				// 更新上次误差return Realize;    								// 返回实际值
}// 增量式PID
float PID_Increase(ERROR *err, PID *pid, float feedback, float reference)
{float iError, Increase;    						// 当前误差  最后得出的实际增量iError = reference - feedback;    				// 计算当前误差Increase = pid->P * (iError - err->Last_Error)    //比例P+ pid->I * iError      							//积分I+ pid->D * (iError - 2 * err->Last_Error + err->Previous_Error);    //微分Derr->Previous_Error = err->Last_Error;    		// 更新前次误差err->Last_Error = iError;		  				// 更新上次误差return Increase;								// 返回增量
}PID pid = { 0.1, 0.01, 0.001, 100.0 };
ERROR err = { 0.0, 0.0, 0.0 };float ref = 100;		//参考值
float sample1 = 0, sample2 = 0;    //采样值int cnt = 0;
float value1[200];
float value2[200];
void main()
{InitSysCtrl();InitPieCtrl();IER = 0x0000;IFR = 0x0000;InitPieVectTable();while (1){sample1 += PID_Realize(&err, &pid, sample1, ref);sample2 += PID_Increase(&err, &pid, sample2, ref);//DELAY_US(1000*200);if (cnt < 200){value1[cnt] = sample1;value2[cnt] = sample2;cnt++;}elsecnt = 0;}
}

http://www.ppmy.cn/server/13247.html

相关文章

Ubuntu系统开机长

Ubuntu系统开机长 1. 检查开机自启动软件的所占时间2. 将耗时最高的禁止开机自启动 1. 检查开机自启动软件的所占时间 systemd-analyze blame2. 将耗时最高的禁止开机自启动 sudo systemctl disable networking.service这个耗时是有阈值的&#xff0c;一般大于15s的算&#x…

算法课程笔记——如何前缀和

要静态数组 注意取模的概念 我想的就是把字母变成数字1&#xff0c;-1&#xff0c;直接遍历一次计数为0就行了

美国共和党正式提出国防量子加速法案

内容来源&#xff1a;量子前哨&#xff08;ID&#xff1a;Qforepost&#xff09; 文丨王珩 浪味仙 排版丨沛贤 深度好文&#xff1a;2500字丨15分钟阅读 摘要&#xff1a;国会议员Elise Stefanik和参议员Marsha Blackburn提出了《国防量子加速法案》。立法者表示&#xff0c;…

面试题6 -- 腐烂的苹果(多源DFS)

腐烂的苹果&#xff08;多源DFS&#xff09; 文章目录 腐烂的苹果&#xff08;多源DFS&#xff09;解题思路代码拟写及相关注意点总结 题目链接&#xff1a; 腐烂的苹果_多源DFS 描述 给定一个 nm 的网格&#xff0c;其中每个单元格中可能有三种值中的一个 0 , 1 , 2。 其中 0…

人民法院对外委托司法鉴定的完善

司法鉴定是指司法鉴定人在诉讼活动中运用科学技术或者专门知识对诉讼涉及的专门性问题进行鉴别和判断并提供鉴定意见的活动。人民法院在案件审理过程中&#xff0c;对于涉及的专门性问题&#xff0c;需要司法鉴定机构予以专业辅助&#xff0c;民事诉讼法将司法鉴定意见列为诉讼…

kafka---topic详解

一、分区与高可用 在Kafka中,事件(events 事件即消息)是以topic的形式进行组织的;同时topic是分区(partitioned)的,这意味着一个topic分布在Kafka broker上的多个“存储桶”(buckets)上。这种数据的分布式放置对于可伸缩性非常重要,因为它允许客户端应用程序同时从多个…

常见大厂面试题(SQL)01

知乎问答最大连续回答问题天数大于等于3天的用户及其对应等级 1.描述 现有某乎问答创作者信息表author_tb如下(其中author_id表示创作者编号、author_level表示创作者级别&#xff0c;共1-6六个级别、sex表示创作者性别)&#xff1a; author_id author_level sex 101 …

2024最新大厂C++面试真题合集,大厂面试百日冲刺 bay3

春招-Bilibili - 后端 Bilibili - 研发工程师-推荐搜索方向 一面 DPDK了解吗 DPDK是一套开源高性能网络数据平面库&#xff0c;用以加速网络包的处理速度&#xff0c;实现高性能的包交换。通过绕过传统的操作系统网络堆栈&#xff0c;直接在用户空间处理网络包&#xff0c;D…