温度控制的PID算法的C语言程序.doc

我的题目是：基于PID算法的温度控制系统
89C51单片机，通过键盘输入预设值，与DS18B20测得的实际值做比较，然后驱动制冷或加热电路。用keil C语言来实现PID的控制。

最佳答案
//PID算法温控C语言2008-08-17 18:58
#include<reg51.h>
#include<intrins.h>
#include<math.h>
#include<string.h>
struct PID {
unsigned int SetPoint; // 设定目标 Desired Value
unsigned int Proportion; // 比例常数 Proportional Const
unsigned int Integral; // 积分常数 Integral Const
unsigned int Derivative; // 微分常数 Derivative Const
unsigned int LastError; // Error[-1]
unsigned int PrevError; // Error[-2]
unsigned int SumError; // Sums of Errors
};
struct PID spid; // PID Control Structure
unsigned int rout; // PID Response (Output)
unsigned int rin; // PID Feedback (Input)
sbit data1=P1^0;
sbit clk=P1^1;
sbit plus=P2^0;
sbit subs=P2^1;
sbit stop=P2^2;
sbit output=P3^4;
sbit DQ=P3^3;
unsigned char flag,flag_1=0;
unsigned char high_time,low_time,count=0;//占空比调节参数
unsigned char set_temper=35;
unsigned char temper;
unsigned char i;
unsigned char j=0;
unsigned int s;
/***********************************************************
延时子程序,延时时间以12M晶振为准,延时时间为30us×time
***********************************************************/
void delay(unsigned char time)
{
unsigned char m,n;
for(n=0;n<time;n++)
for(m=0;m<2;m++){}
}
/***********************************************************
写一位数据子程序
***********************************************************/
void write_bit(unsigned char bitval)
{
EA=0;
DQ=0; /*拉低DQ以开始一个写时序*/
if(bitval==1)
{
_nop_();
DQ=1; /*如要写1，则将总线置高*/
}
delay(5); /*延时90us供DA18B20采样*/
DQ=1; /*释放DQ总线*/
_nop_();
_nop_();
EA=1;
}
/***********************************************************
写一字节数据子程序
***********************************************************/
void write_byte(unsigned char val)
{
unsigned char i;
unsigned char temp;
EA=0; /*关中断*/
TR0=0;
for(i=0;i<8;i++) /*写一字节数据，一次写一位*/
{
temp=val>>i; /*移位操作，将本次要写的位移到最低位*/
temp=temp&1;
write_bit(temp); /*向总线写该位*/
}
delay(7); /*延时120us后*/
// TR0=1;
EA=1; /*开中断*/
}
/***********************************************************
读一位数据子程序
***********************************************************/
unsigned char read_bit()
{
unsigned char i,value_bit;
EA=0;
DQ=0; /*拉低DQ，开始读时序*/
_nop_();
_nop_();
DQ=1; /*释放总线*/
for(i=0;i<2;i++){}
value_bit=DQ;
EA=1;
return(value_bit);
}
/***********************************************************
读一字节数据子程序
***********************************************************/
unsigned char read_byte()
{
unsigned char i,value=0;
EA=0;
for(i=0;i<8;i++)
{
if(read_bit()) /*读一字节数据，一个时序中读一次，并作移位处理*/
value|=0x01<<i;
delay(4); /*延时80us以完成此次都时序，之后再读下一数据*/
}
EA=1;
return(value);
}
/***********************************************************
复位子程序
***********************************************************/
unsigned char reset()
{
unsigned char presence;
EA=0;
DQ=0; /*拉低DQ总线开始复位*/
delay(30); /*保持低电平480us*/
DQ=1; /*释放总线*/
delay(3);
presence=DQ; /*获取应答信号*/
delay(28); /*延时以完成整个时序*/
EA=1;
return(presence); /*返回应答信号，有芯片应答返回0,无芯片则返回1*/
}
/***********************************************************
获取温度子程序
***********************************************************/
void get_temper()
{
unsigned char i,j;
do
{
i=reset(); /*复位*/
}while(i!=0); /*1为无反馈信号*/
i=0xcc; /*发送设备定位命令*/
write_byte(i);
i=0x44; /*发送开始转换命令*/
write_byte(i);
delay(180); /*延时*/
do
{
i=reset(); /*复位*/
}while(i!=0);
i=0xcc; /*设备定位*/
write_byte(i);
i=0xbe; /*读出缓冲区内容*/
write_byte(i);
j=read_byte();
i=read_byte();
i=(i<<4)&0x7f;
s=(unsigned int)(j&0x0f);
s=(s*100)/16;
j=j>>4;
temper=i|j; /*获取的温度放在temper中*/
}
/*====================================================================================================
Initialize PID Structure
=====================================================================================================*/
void PIDInit (struct PID *pp)
{
memset ( pp,0,sizeof(struct PID));
}
/*====================================================================================================
PID计算部分
=====================================================================================================*/
unsigned int PIDCalc( struct PID *pp, unsigned int NextPoint )
{
unsigned int dError,Error;
Error = pp->SetPoint - NextPoint; // 偏差
pp->SumError += Error; // 积分
dError = pp->LastError - pp->PrevError; // 当前微分
pp->PrevError = pp->LastError;
pp->LastError = Error;
return (pp->Proportion * Error//比例
+ pp->Integral * pp->SumError //积分项
+ pp->Derivative * dError); // 微分项
}
/***********************************************************
温度比较处理子程序
***********************************************************/
compare_temper()
{
unsigned char i;
if(set_temper>temper)
{
if(set_temper-temper>1)
{
high_time=100;
low_time=0;
}
else
{
for(i=0;i<10;i++)
{ get_temper();
rin = s; // Read Input
rout = PIDCalc ( &spid,rin ); // Perform PID Interation
}
if (high_time<=100)
high_time=(unsigned char)(rout/800);
else
high_time=100;
low_time= (100-high_time);
}
}
else if(set_temper<=temper)
{
if(temper-set_temper>0)
{
high_time=0;
low_time=100;
}
else
{
for(i=0;i<10;i++)
{ get_temper();
rin = s; // Read Input
rout = PIDCalc ( &spid,rin ); // Perform PID Interation
}
if (high_time<100)
high_time=(unsigned char)(rout/10000);
else
high_time=0;
low_time= (100-high_time);
}
}
// else
// {}
}
/*****************************************************
T0中断服务子程序，用于控制电平的翻转 ,40us*100=4ms周期
******************************************************/
void serve_T0() interrupt 1 using 1
{
if(++count<=(high_time))
output=1;
else if(count<=100)
{
output=0;
}
else
count=0;
TH0=0x2f;
TL0=0xe0;
}
/*****************************************************
串行口中断服务程序，用于上位机通讯
******************************************************/
void serve_sio() interrupt 4 using 2
{
/* EA=0;
RI=0;
i=SBUF;
if(i==2)
{
while(RI==0){}
RI=0;
set_temper=SBUF;
SBUF=0x02;
while(TI==0){}
TI=0;
}
else if(i==3)
{
TI=0;
SBUF=temper;
while(TI==0){}
TI=0;
}
EA=1; */
}
void disp_1(unsigned char disp_num1[6])
{
unsigned char n,a,m;
for(n=0;n<6;n++)
{
// k=disp_num1[n];
for(a=0;a<8;a++)
{
clk=0;
m=(disp_num1[n]&1);
disp_num1[n]=disp_num1[n]>>1;
if(m==1)
data1=1;
else
data1=0;
_nop_();
clk=1;
_nop_();
}
}
}
/*****************************************************
显示子程序
功能：将占空比温度转化为单个字符，显示占空比和测得到的温度
******************************************************/
void display()
{
unsigned char code number[]={0xfc,0x60,0xda,0xf2,0x66,0xb6,0xbe,0xe0,0xfe,0xf6};
unsigned char disp_num[6];
unsigned int k,k1;
k=high_time;
k=k%1000;
k1=k/100;
if(k1==0)
disp_num[0]=0;
else
disp_num[0]=0x60;
k=k%100;
disp_num[1]=number[k/10];
disp_num[2]=number[k%10];
k=temper;
k=k%100;
disp_num[3]=number[k/10];
disp_num[4]=number[k%10]+1;
disp_num[5]=number[s/10];
disp_1(disp_num);
}
/***********************************************************
主程序
***********************************************************/
main()
{
unsigned char z;
unsigned char a,b,flag_2=1,count1=0;
unsigned char phil[]={2,0xce,0x6e,0x60,0x1c,2};
TMOD=0x21;
TH0=0x2f;
TL0=0x40;
SCON=0x50;
PCON=0x00;
TH1=0xfd;
TL1=0xfd;
PS=1;
EA=1;
EX1=0;
ET0=1;
ES=1;
TR0=1;
TR1=1;
high_time=50;
low_time=50;
PIDInit ( &spid ); // Initialize Structure
spid.Proportion = 10; // Set PID Coefficients
spid.Integral = 8;
spid.Derivative =6;
spid.SetPoint = 100; // Set PID Setpoint
while(1)
{
if(plus==0)
{
EA=0;
for(a=0;a<5;a++)
for(b=0;b<102;b++){}
if(plus==0)
{
set_temper++;
flag=0;
}
}
else if(subs==0)
{
for(a=0;a<5;a++)
for(b=0;a<102;b++){}
if(subs==0)
{
set_temper--;
flag=0;
}
}
else if(stop==0)
{
for(a=0;a<5;a++)
for(b=0;b<102;b++){}
if(stop==0)
{
flag=0;
break;
}
EA=1;
}
get_temper();
b=temper;
if(flag_2==1)
a=b;
if((abs(a-b))>5)
temper=a;
else
temper=b;
a=temper;
flag_2=0;
if(++count1>30)
{
display();
count1=0;
}
compare_temper();
}
TR0=0;
z=1;
while(1)
{
EA=0;
if(stop==0)
{
for(a=0;a<5;a++)
for(b=0;b<102;b++){}
if(stop==0)
disp_1(phil);
// break;
}
EA=1;
}
}
//DS18b20 子程序
#include <REG52.H>
sbit DQ=P2^1; //定义端口

typedef unsigned char byte;
typedef unsigned int word;
//延时
void delay(word useconds)
{
for(;useconds>0;useconds--);
}
//复位
byte ow_reset(void)
{
byte presence;
DQ=0; //DQ低电平
delay(29); //480us
DQ=1; //DQ高电平
delay(3); //等待
presence=DQ; //presence信号
delay(25);
return(presence);
} //0允许，1禁止
//从1-wire 总线上读取一个字节
byte read_byte(viod)
{
byte i;
byte value=0;
for (i=8;i>0;i--)
{
value>>=1;
DQ=0;
DQ=1;
delay(1);

if(DQ)value|=0x80;
delay(6);
}
return(value);
}

//向1-wire总线上写一个字节
void write_byte(char val)
{
byte i;
for (i=8;i>0;i--) //一次写一个字节
{
DQ=0;
DQ=val&0x01;
delay(5);
DQ=1;
val=val/2;
}
delay(5);
}
//读取温度

char Read_Temperature(void)
{
union{
byte c[2];
int x;
}temp;

ow_reset();
write_byte(0xcc);
write_byte(0xBE);
temp.c[1]=read_byte();
temp.c[0]=read_byte();
ow_reset();
write_byte(0xCC);
write_byte(0x44);
return temp.x/2;
}

参考资料：你把这两个程序组合就可以了
PID算法

    PID算法是本程序中的核心部分。我们采用PID模糊控制技术,通过Pvar、Ivar、Dvar（比例、积分、微分）三方面的结合调整形成一个模糊控制来解决惯性温度误差问题。其原理如下：
    本系统的温度控制器的电热元件之一是发热丝。发热丝通过电流加热时，内部温度都很高。当容器内温度升高至设定温度时，温度控制器会发出信号停止加热。但这时发热丝的温度会高于设定温度，发热丝还将会对被加热的器件进行加热，即使温度控制器发出信号停止加热，被加热器件的温度还往往继续上升几度，然后才开始下降。当下降到设定温度的下限时，温度控制器又开始发出加热的信号，开始加热，但发热丝要把温度传递到被加热器件需要一定的时候，这就要视发热丝与被加热器件之间的介质情况而定。通常开始重新加热时，温度继续下降几度。所以，传统的定点开关控制温度会有正负误差几度的现象，但这不是温度控制器本身的问题，而是整个热系统的结构性问题，使温度控制器控温产生一种惯性温度误差。
    增量式PID算法的输出量为
               ΔUn = Kp[(en-en-1)+(T/Ti)en+(Td/T)(en-2*en-1+en-2)]
    式中，en、en-1、en-2分别为第n次、n-1次和n-2次的偏差值，Kp、Ti、Td分别为比例系数、积分系数和微分系数，T为采样周期。
    计算机每隔固定时间 T将现场温度与用户设定目标温度的差值带入增量式PID算法公式，由公式输出量决定PWM方波的占空比，后续加热电路根据此PWM方波的占空比决定加热功率。现场温度与目标温度的偏差大则占空比大，加热电路的加热功率大，使温度的实测值与设定值的偏差迅速减少；反之，二者的偏差小则占空比减小，加热电路加热功率减少，直至目标值与实测值相等，达到自动控制的目的。
    PID参数的选择是实验成败的关键，它决定了温度控制的准确度。数字PID调节器参数的整定可以仿照模拟PID调节器参数整定的各种方法，根据工艺对控制性能的要求，决定调节器的参数。各个参数对系统性能的影响如下：
    ①比例系数P对系统性能的影响：比例系数加大，使系统的动作灵敏，速度加快，稳态误差减小；P偏大，振荡次数加多，调节时间加长；P太大时，系统会趋于不稳定；P太小，又会使系统的动作缓慢。P可以选负数，这主要是由执行机构、传感器以及控制对象的特性决定的。如果P的符号选择不当对象测量值就会离控制目标的设定值越来越远，如果出现这样的情况P的符号就一定要取反。
    ②积分控制I对系统性能的影响：积分作用使系统的稳定性下降，I小（积分作用强）会使系统不稳定，但能消除稳态误差，提高系统的控制精度。
③微分控制D对系统性能的影响：微分作用可以改善动态特性，D偏大时，超调量较大，调节时间较短；D偏小时，超调量也较大，调节时间也较长；只有D合适，才能使超调量较小，减短调节时间。

温底控制PID的算法设计及实现

-
PID 简介
    PID（Proportional Integral Derivative）控制是控制工程中技术成熟、应用广泛的一种控制策略，经过长期的工程实践，已形成了一套完整的控制方法和典型的结构。它不仅适用于数学模型已知的控制系统中，而且对于大多数数学模型难以确定的工业过程也可应用，在众多工业过程控制中取得了满意的应用效果。
    PID 工作基理：由于来自外界的各种扰动不断产生，要想达到现场控制对象值保持恒定的目的，控制作用就必须不断的进行。若扰动出现使得现场控制对象值(以下简称被控参数)发生变化，现场检测元件就会将这种变化采集后经变送器送至PID 控制器的输入端，并与其给定值(以下简称SP 值)进行比较得到偏差值(以下简称e 值)，调节器按此偏差并以我们预先设定的整定参数控制规律发出控制信号，去改变调节器的开度，使调节器的开度增加或减少，从而使现场控制对象值发生改变，并趋向于给定值(SP 值)，以达到控制目的 ，如图 1 所示，其实PID 的实质就是对偏差（e 值）进行比例、积分、微分运算，根据运算结果控制执行部件的过程。

图1 模拟PID 控制系统原理图
    PID 控制器的控制规律可以描述为：
  (1)
    比例（P）控制能迅速反应误差，从而减小稳态误差。但是，比例控制不能消除稳态误差。比例放大系数的加大，会引起系统的不稳定。积分（I）控制的作用是：只要系统有误差存在，积分控制器就不断地积累，输出控制量，以消除误差。因而，只要有足够的时间，积分控制将能完全消除误差，使系统误差为零，从而消除稳态误差。积分作用太强会使系统超调加大，甚至使系统出现振荡。微分（D）控制可以减小超调量，克服振荡，使系统的稳定性提高，同时加快系统的动态响应速度，减小调整时间，从而改善系统的动态性能。根据不同的被控对象的控制特性，又可以分为P、PI、PD、PID 等不同的控制模型。
数字PID 的实现
    在连续-时间控制系统（模拟PID 控制系统）中，PID 控制器应用得非常广泛。其设计技术成熟，长期以来形成了典型的结构，参数整定方便，结构更改灵活，能满足一般的控制要求。随着计算机的快速发展，人们将计算机引入到PID 控制领域，也就出现了数字式PID 控制。
    由于计算机基于采样控制理论，计算方法也不能沿袭传统的模拟PID 控制算法（如公式1 所示），所以必须将控制模型离散化，离散化的方法：以T 为采样周期，k 为采样序号，用求和的形式代替积分，用增量的形式（求差）代替微分，这样可以将连续的PID 计算公式离散：
 (2)
式1 就可以离散为：
 (3)
或者：
 (4)
    这样就可以让计算机或者单片机通过采样的方式实现PID 控制，具体的PID 控制又分为位置式PID 控制和增量式PID 控制，公式4 给出了控制量的全部大小，所以称之为全量式或者位置式控制；如果计算机只对相邻的两次作计算，只考虑在前一次基础上，计算机输出量的大小变化，而不是全部输出信息的计算，这种控制叫做增量式PID 控制算法，其实质就是求Δμ的大小，而   Δμk =μk -μk-1 ；所以将式4 做自减变换有：
 (5)
其中

温度控制PID 算法设计
    本设计利用了上面所介绍的位置式PID 算法，将温度传感器采样输入作为当前输入，然后与设定值进行相减得偏差ek，然后再对之进行PID 运算产生输出结果fOut，然后让fOut 控制定时器的时间进而控制加热器。为了方便PID 运算，首先建立一个PID 的结构体数据类型，该数据类型用于保存PID 运算所需要的P、I、D 系数，以及设定值，历史误差的累加和等信息：

typedef struct PID
{
float SetPoint; // 设定目标 Desired Value
float Proportion; // 比例系数 Proportional Const
float Integral; // 积分系数 Integral Const
float Derivative; // 微分系数 Derivative Const
int LastError; // 上次偏差
int SumError; // 历史误差累计值
} PID;
PID stPID; // 定义一个stPID 变量

    下面是PID 运算的算法程序，通过PID 运算返回fOut，fOut 的值决定是否加热，加热时间是多少。
PID 运算的C 实现代码：

float PIDCalc( PID *pp, int NextPoint )
{
int dError,Error;
Error = pp->SetPoint*10 - NextPoint; // 偏差，设定值减去当前采样值
pp->SumError += Error; // 积分，历史偏差累加
dError = Error-pp->LastError; // 当前微分，偏差相减
pp->PrevError = pp->LastError; // 保存
pp->LastError = Error;
+ pp->Integral * pp->SumError // 积分项
- pp->Derivative * dError // 微分项
);
}

    在实际运算时，由于水具有很大的热惯性，而且PID 运算中的I（积分项）具有非常明显的延迟效应所以不能保留，我们必须把积分项去掉，相反D（微分项）则有很强的预见性，能够加快反应速度，抑制超调量，所以积分作用应该适当加强才能达到较佳的控制效果，系统最终选择PD 控制方案，下面C 代码所示为PD 控制的实现过程：

float PIDCalc( PID *pp, int NextPoint )
{
int dError,Error;
Error = pp->SetPoint*10 - NextPoint; // 偏差，设定值减去当前采样值
dError = Error-pp->LastError; // 当前微分，偏差相减
pp->PrevError = pp->LastError; // 保存
pp->LastError = Error;
return (pp->Proportion * Error // 比例项
- pp->Derivative * dError // 微分项
);
}

温度控制实现
    通过温度的PID 运算，产生结果fOut，该参数决定是否加热，加热时间是多长。该程序如下：

stPID.Proportion = 2; //设置PID 比例值
stPID.Integral = 0; //设置PID 积分值
stPID.Derivative = 5; //设置PID 微分值
fOut = PIDCalc ( &stPID,(int)(fT*10) ); //PID 计算
if(fOut<=0)
*P_IOA_Buffer &= 0xff7f; //温度高于设定值，关闭加热器
else
*P_IOA_Buffer |= 0x0080; //温度低于设定值，打开加热器

    加热时间由主函数计算，由TimerB 中断控制。主程序中通过PIDCalc 函数得到fOut 参数，如果该参数大于“0”，则开启加热器。IRQ2_TMB 中断一直处于允许状态，每进入一次IRQ2_TMB 中断，fOut 参数减1，直到fOut = 0，停止加热。如果PIDCalc 计算结果比较大说明离目标温度相差较大，则加热时间比较长，如果计算结果比较小，说明离目标温度相差较小，加热时间相对较短。

基于PID算法和89C52单片机的温度控制系统
作者：张艳艳 安徽电子信息职业技术学院  来源：现代电子技术  发布时间：2009-12-22 17:36:09  [收 藏] [评 论]
基于PID算法和89C52单片机的温度控制系统
0 引 言
    温控技术无论是在工业生产，还是日常生活中都起着非常重要的作用。在冶金、石油、化工、电力和现代农业等行业，温度是极为重要而又普遍的热工参数之一，在普通家庭里热水器、电饭煲、电烤箱等依赖于温控技术的家电设备也是必不可少。可以说温度控制技术无处不在。
    常规的温度控制方法以设定温度为临界点，超出设定允许范围即进行温度调控：低于设定值就加热，反之就停止或降温。这种方法实现简单、成本低，但控制效果不理想，控制温度精度不高、容易引起震荡，达到稳定点的时间也长，因此，只能用在精度要求不高的场合。
    而采用PID算法进行温度控制，它具有控制精度高，能够克服容量滞后的特点，特别适用于负荷变化大、容量滞后较大、控制品质要求又很高的控制系统。
    单片机作为控制系统中必不可少的部分，在各个领域得到了广泛的应用，用单片机进行实时