单片机算法c语言程序,51单片机PID的算法实现程序C语言

/*

用整型变量来实现PID算法，由于是用整型数来做的,所以也不是很精确,但是对于很多的使用场合,这个精度也够了，

关于系数和采样电压全部是放大10倍处理的.所以精度不是很高. 但是也不是那么低,大部分的场合都够了. 实在觉得精度不够,

可以再放大10倍或者100倍处理,但是要注意不超出整个数据类型的范围就可以了.本程序包括PID计算和输出两部分.

当偏差>10度全速加热,偏差在10度以内为PID计算输出.   具体的参考代码参见下面：

*/

//================================================================

// pid.H

// Operation about PID algorithm procedure

// C51编译器  Keil 7.08

//================================================================

// All rights reserved.

//================================================================

#include #include typedef   unsigned   char        uint8;

typedef   unsigned   int         uint16;

typedef   unsigned   long int    uint32;

/函数声明/

void     PIDOutput ();

void     PIDOperation ();

/*/

typedef struct PIDValue

{

uint32      Ek_Uint32[3];                  //差值保存，给定和反馈的差值

uint8       EkFlag_Uint8[3];              //符号，1则对应的为负数，0为对应的为正数

uint8       KP_Uint8;

uint8       KI_Uint8;

uint8       KD_Uint8;

uint16      Uk_Uint16;                 //上一时刻的控制电压

uint16      RK_Uint16;                //设定值

uint16      CK_Uint16;               //实际值

}PIDValueStr;

PIDValueStr  PID;

uint8        out ;                 // 加热输出

uint8        count;               // 输出时间单位计数器

/*

PID = Uk + KP*[E(k)-E(k-1)]+KI*E(k)+KD*[E(k)-2E(k-1)+E(k-2)];(增量型PID算式)

函数入口: RK(设定值),CK(实际值),KP,KI,KD

函数出口: U(K)

//PID运算函数

/

void     PIDOperation (void)

{

uint32       Temp[3];                                        //中间临时变量

uint32       PostSum;                                       //正数和

uint32       NegSum;                                       //负数和

Temp[0] = 0;

Temp[1] = 0;

Temp[2] = 0;

PostSum = 0;

NegSum  = 0;

if( PID.RK_Uint16 > PID.RK_Uint16 )                    //设定值大于实际值否？

{

if( PID.RK_Uint16 - PID.RK_Uint16 >10 )            //偏差大于10否？

{

PID.Uk_Uint16 = 100;    //偏差大于10为上限幅值输出(全速加热)

}

else

{

Temp[0] = PID.RK_Uint16 - PID.CK_Uint16;       //偏差<=10,计算E(k)

PID.EkFlag_Uint8[1]=0;                        //E(k)为正数

//数值移位

PID.Ek_Uint32[2] = PID.Ek_Uint32[1];

PID.Ek_Uint32[1] = PID.Ek_Uint32[0];

PID.Ek_Uint32[0] = Temp[0];

//

if( PID.Ek_Uint32[0] >PID.Ek_Uint32[1] )                            //E(k)>E(k-1)否？

{

Temp[0]=PID.Ek_Uint32[0] - PID.Ek_Uint32[1];           //E(k)>E(k-1)

PID.EkFlag_Uint8[0]=0;                                       //E(k)-E(k-1)为正数

}

else

{

Temp[0]=PID.Ek_Uint32[0] - PID.Ek_Uint32[1];        //E(k)

PID.EkFlag_Uint8[0]=1;                                                //E(k)-E(k-1)为负数

}

//

Temp[2]=PID.Ek_Uint32[1]*2 ;                                             // 2E(k-1)

if( (PID.Ek_Uint32[0]+ PID.Ek_Uint32[2])>Temp[2] )            //E(k-2)+E(k)>2E(k-1)否？

{

Temp[2]=(PID.Ek_Uint32[0]+ PID.Ek_Uint32[2])-Temp[2];     //E(k-2)+E(k)>2E(k-1)

PID.EkFlag_Uint8[2]=0;                                           //E(k-2)+E(k)-2E(k-1)为正数

}

else

{

Temp[2]=Temp[2]-(PID.Ek_Uint32[0]+ PID.Ek_Uint32[2]);  //E(k-2)+E(k)<2E(k-1)

PID.EkFlag_Uint8[2]=1;                                      //E(k-2)+E(k)-2E(k-1)为负数

}

//

Temp[0] = (uint32)PID.KP_Uint8 * Temp[0];                        // KP*[E(k)-E(k-1)]

Temp[1] = (uint32)PID.KI_Uint8 * PID.Ek_Uint32[0];              // KI*E(k)

Temp[2] = (uint32)PID.KD_Uint8 * Temp[2];                      // KD*[E(k-2)+E(k)-2E(k-1)]

/*以下部分代码是讲所有的正数项叠加，负数项叠加*/

/KP*[E(k)-E(k-1)]/

if(PID.EkFlag_Uint8[0]==0)

PostSum += Temp[0];                                    //正数和

else

NegSum += Temp[0];                                    //负数和

/* KI*E(k)/

if(PID.EkFlag_Uint8[1]==0)

PostSum += Temp[1];                                 //正数和

else

;                                                 //空操作，E(K)>0

/KD*[E(k-2)+E(k)-2E(k-1)]/

if(PID.EkFlag_Uint8[2]==0)

PostSum += Temp[2];                               //正数和

else

NegSum += Temp[2];                             //负数和

/*U(K)*/

PostSum += (uint32)PID.Uk_Uint16;

if(PostSum > NegSum )                         // 是否控制量为正数

{

Temp[0] = PostSum - NegSum;

if( Temp[0] < 100 )                         //小于上限幅值则为计算值输出

PID.Uk_Uint16 = (uint16)Temp[0];

else

PID.Uk_Uint16 = 100;                     //否则为上限幅值输出

}

else                                     //控制量输出为负数，则输出0(下限幅值输出)

PID.Uk_Uint16 = 0;

}

}

else

{

PID.Uk_Uint16 = 0;

}

}

/*

函数入口: U(K)

函数出口: out(加热输出)

//PID运算植输出函数

/

void     PIDOutput (void)

{

static  int i;

i=PID.Uk_Uint16;

if(i==0)

out=1;

else out=0;

if((count++)==5)//如定时中断为40MS,40MS*5=0.2S(输出时间单位),加热周期20S(100等份)

{              //每20S PID运算一次

count=0;

i--;

}

}

关注eeworld公众号

快捷获取更多信息

关注eeworld服务号

享受更多官方福利

推荐阅读

STC15单片机实验名称：使用NTC电阻测量温度实验内容：        使用NTC电阻测量温度并显示在数码管上        读取DS18B20模块测量温度作为参考实验器材：        STC15W408AS_DIP16 x1        DS18B20  x1        TM1637数码管 x1

发表于 2020-12-17

/*一款电子钟程序，此电子钟是四位电子钟数码管组成的，功能分别有走时、响闹、走时间调整、闹钟时间调整等功能！调整部分分别由三个按键实现，一个功能键，一个＋(加)键，另一个为－(减)键，其它功能键是复用功能！*/硬键介绍：数码管个位.十位.百位.千位.小数点分别接P0.0 P0.1 P0.2 P0.3 P0.4。蜂鸣器接P3.5脚，功能按键接P3.2 加按键接P3.3 减按键接P3.4以下是两张图片！！！#include //头文件#define uchar unsigned char//宏定义#define uint unsigned int//宏定义uchar shi,fen,nshi,nfen,ma

发表于 2020-12-17

51单片机做的温度显示，温度显示在LCD1602液晶屏上。然后按键可以调整温度阈值，温度高于或低于所设温度，蜂鸣器就会响……单片机源程序如下:/***************************************************************************************                              基于单片机的实时温度报警系统设计        &nbsp

发表于 2020-12-17

程序代码 仿真图单片机源程序如下:#include#include/*********************宏定义************************/#define uchar unsigned char#define uint unsigned int/************位定义*****************************/sbit Trig=P3^7; //超声波模块触发位sbit LEDRED=P1^0;sbit LEDGREEN=P1^2;sbit speaker=P1^7;sbit person=P1^6;sbit

发表于 2020-12-09

电子秤是基于STC89C51单片机，以及HX711称重模块，利用单片机控制A/D转换，屏幕数据显示，键盘实时输入，设计方案能够易于各种场合使用。单片机源程序如下：#include #include #include #include "main.h"#include "HX711.h"#include "eeprom52.h"#define uchar unsigned char#define uint  unsigned intunsigned long HX711

发表于 2020-12-09