简介：这份文档面向自动化、机电类专业学生及从事温度控制开发的工程人员，围绕恒温水箱被控对象，给出以单片机为控制器的模糊PID恒温控制系统设计方案，重点解决温度控制滞后、大惯性条件下快速性与准确性难以兼顾的问题，温度范围覆盖40℃至90℃，误差要求控制在0.5℃以内。资源包仅含1个doc文件，大小约698KB，内容以毕业设计论文形式组织，包含中英文摘要、目录、系统总体方案、AT89S52单片机与DS18B20测温电路、温度显示及可控硅加热控制等硬件设计，并讨论模糊控制规则、PID参数调整与两者结合的实现思路。目前已有297人学习，适合需要参考完整课题结构、撰写论文或搭建恒温控制实验原型的读者，可从中获取需求分析、硬件选型、算法设计到系统测试的成文框架与关键设计要点。
1. 从 40℃ 到 90℃ 锁死在 0.5℃：一个模糊 PID 恒温箱的完整拆解
温度控制里有个反直觉的地方：你把 PID 三个参数调到曲线漂漂亮亮，换个加热丝、换台水箱、甚至室温从 20℃ 变成 30℃，那条曲线立刻变形。这不是调参水平问题，是定参数 PID 的固有短板——它假定被控对象模型不变，而恒温箱恰恰是单向性、大惯性、时变的对象。这套基于模糊 PID 恒温算法的毕业设计，就是冲着这个痛点做的：以恒温水箱为对象，51 内核单片机做控制器，DS18B20 测温，双向可控硅调功，要求在 40℃~90℃ 任意设定点上把误差压在 0.5℃ 以内，并且参数能随工况自己动。它适合正在做温控类课设、想搞明白模糊 PID 到底怎么落到 8 位机上的人，也适合手上有台老式恒温设备想改成自适应控制的人。
2. 系统总体架构：一次闭环里五个环节各干什么活
这套系统的骨架不复杂，但每个环节的边界划得很清楚，拆开看是一条完整的前向通道加一条反馈通道。理解这条链路，后面的代码才有地方放。
2.1 闭环链路的五段划分
从能量流看，主回路是「单片机 → 光电隔离 → 驱动 → 双向可控硅 → 加热丝 → 水箱」；从信号流看，反馈回路是「水箱 → DS18B20 → 单总线 → 单片机」。两条回路在水箱这个物理对象上闭合，构成典型的温度闭环。
前向通道的末端是可控硅，它串在 220V 市电和加热丝之间，靠导通角或者导通时间来调节加在加热丝上的平均功率。这里用的是通断式调功——单片机的 P1.7 口输出高低电平，经光电耦合器隔离后送到可控硅控制端，高电平期间可控硅导通、加热丝通电，低电平期间关断。也就是说单片机控制的不是「加热多热」，而是「一个时间窗口内加热多久」，这在热惯性大的水箱上完全够用，因为水箱本身就是一个天然的低通滤波器，把断续的通断平滑成了连续温升。
反馈通道用 DS18B20，单总线器件，一根 DQ 线既传命令又回数据，省掉了 A/D 转换这一整块电路。这一点对 51 单片机很关键：AT89S52 片内没有 ADC，如果选模拟传感器（比如热电偶加运放），就必须外挂 ADC0809 之类的芯片，布线和时序都会复杂一大截。DS18B20 直接出数字量，12 位分辨率对应 0.0625℃，理论精度已经远超 0.5℃ 的指标要求，余量是够的。
人机接口这部分容易被忽略，但它决定了这套东西能不能脱离电脑独立跑。显示用三只共阳数码管，通过两片 74LS164 串行移位寄存器驱动，只占用单片机 P3.0、P3.1 两根线；按键四个，接在 P1 口低四位，分别是启动、百位加、十位加、个位加；报警用一个有源蜂鸣器，P1.6 输出低电平经反相器驱动。停止加热时数码管显示设定温度，启动加热后切换显示当前实测温度，这个显示策略在调试阶段特别有用，一眼就能看出系统是「没加热」还是「加热了但没升温」。
提示：显示切换逻辑一定要做成「启动后显示实测值」，否则现场调试时你根本分不清是传感器坏了还是根本没启动。
2.2 器件选型对照与各自承担的职责
选型这件事，毕设项目和工业产品考虑的角度不完全一样。工业上优先看可靠性和供货，毕设更看「能不能用最少的外围把功能闭环跑通」。这套方案在两端之间取了个平衡点，下面这张表是我拆下来之后整理的职责划分，改项目时照着改就行。
| 器件 | 承担职责 | 关键参数 | 替换时的注意点 |
|---|---|---|---|
| AT89S52 | 主控、算法运算、时序调度 | 8K Flash、256B RAM、32 个 I/O、3 个 16 位定时器 | 换 STC89C52 可直接代换，注意 ISP 下载脚冲突 |
| DS18B20 | 温度采集，数字量输出 | -55℃~125℃，9~12 位可配，12 位时 0.0625℃/LSB | 换 DS18B20 防水封装时线缆电容会拖慢时序 |
| 双向可控硅 | 加热功率调节 | 耐压 ≥600V，电流按加热丝额定值留 2 倍余量 | 感性负载要加 RC 吸收，纯阻性可省 |
| 光电耦合器 | 强电弱电隔离 | 常用双向可控硅型输出 | 别用普通三极管型去驱可控硅，触发极性不对 |
| 74LS164 ×2 | 串行转并行驱动数码管 | 8 位串入并出，带锁存时钟 | 换 74HC164 时注意驱动电流下降 |
| LM7805 | 5V 稳压 | 输入需 ≥7V，压差 2~3V | 输入滤波电容不能省，否则复位会抖 |
| 有源蜂鸣器 | 超温报警 | 5V 直驱，反相器推动 | 无源蜂鸣器要配方波，别买错 |
这张表里最值得说的是可控硅和光耦的搭配。单片机的 I/O 口输出电流只有几十微安到几毫安量级，直接驱可控硅门极是绝对带不动的，必须过一级驱动；而可控硅这侧接的是 220V 市电，一旦击穿就会把高压灌回单片机，烧掉的不只是一个芯片。所以光电耦合器在这里不是「可选优化」，是安全底线，两个方向的作用都得有。
2.3 采样节拍推导：6MHz 晶振怎么凑出 15 秒控制周期
温度是大惯性对象，采样周期定得太短没意义——水箱热时间常数通常在几分钟量级，你 10ms 采一次，相邻两次读数几乎一样，纯属浪费 CPU 和总线。这套系统把采样控制周期定在 15 秒，是个比较务实的值，既能及时响应，又不会让可控硅的通断过于频繁。
问题在于，6MHz 晶振下 51 单片机的机器周期是 2μs，16 位定时器最大计数值 65536，单次最长定时只有 131ms 左右，离 15 秒差了两个数量级。标准解法是「硬件定时 + 软件计数」两级串起来：定时器 0 设成 50ms 中断一次，中断服务程序里对一个全局变量累加，累加到 300 次正好 15 秒。
/* 6MHz 晶振，12T 模式，机器周期 = 2us */ /* 50ms = 25000 个机器周期，初值 = 65536 - 25000 = 40536 = 0x9E58 */ #define T0_H 0x9E #define T0_L 0x58 #define TICK_15S 300 /* 50ms × 300 = 15s */ volatile unsigned int tick_cnt = 0; volatile bit sample_flag = 0; /* 15s 到点标志，主循环轮询 */ void timer0_isr(void) interrupt 1 { TH0 = T0_H; /* 必须软件重装，方式 1 不自动重载 */ TL0 = T0_L; if (++tick_cnt >= TICK_15S) { tick_cnt = 0; sample_flag = 1; /* 只置标志，中断里不做好事 */ } } void timer0_init(void) { TMOD &= 0xF0; /* 清 T0 控制位，保留 T1 配置 */ TMOD |= 0x01; /* T0 方式 1：16 位定时器 */ TH0 = T0_H; TL0 = T0_L; ET0 = 1; /* 开 T0 中断 */ EA = 1; /* 开总中断 */ TR0 = 1; /* 启动计数 */ }这里有两个细节值得说透。第一，方式 1 是「不自动重装」的，中断里必须手动写 TH0/TL0，而且要在中断服务程序的最前面写，否则中断响应本身的那几个机器周期会被算进下一次定时，误差一点点累积，跑几个小时就偏出好几秒。如果你对周期精度有要求，改用方式 2（8 位自动重装）配预分频，或者干脆换带自动重载的定时器 2。第二，中断服务程序里只置标志、不做运算，这是 51 上的铁律——中断里做的事越多，主循环被切得越碎，串口、按键这些对时序敏感的部分就越容易出错。
主程序的骨架就很直白了：初始化 → 等 15 秒标志 → 采一次温度 → 算一次模糊 PID → 更新输出占空比 → 刷新显示 → 回来继续等。
void main(void) { unsigned char set_temp = 60; /* 默认设定 60℃ */ system_init(); /* 端口、定时器、DS18B20 全部初始化 */ display_set_temp(set_temp); /* 上电先显示设定值 */ while (1) { key_scan(&set_temp); /* 按键随时响应，不受采样周期约束 */ if (sample_flag) { int fb = ds18b20_read_temp_x16(); /* 单位 1/16℃ */ if (fb != TEMP_ERR) { /* 读数有效才进算法 */ int duty = fuzzy_pid(set_temp * 16, fb); triac_set_duty(duty); /* 写入全局占空比，由过零中断消费 */ display_temp(fb); check_alarm(fb, set_temp); } sample_flag = 0; } } }注意fuzzy_pid的入参单位是「1/16℃」，也就是直接把 DS18B20 的原始码当定点数用，全程不碰浮点。8 位机上做一个 float 乘法要几百个机器周期，15 秒一次虽然不心疼，但一旦你把采样周期改短、或者加了串口通信，浮点就会变成瓶颈。用定点是这套代码里最要紧的一个习惯。
3. 感知层与执行层：单总线时序和可控硅调功怎么落地
算法再漂亮，底层时序错一位，读回来的就是垃圾数据。这一章把 DS18B20 和可控硅这两块的实现摊开讲。
3.1 DS18B20 单总线时序与 0.0625℃ 的换算
单总线之所以叫单总线，是因为它把所有通信都压在一根线上，靠电平持续的长短来区分 0 和 1。这带来一个直接后果：时序精度全靠软件延时撑着，中断一打断，时序就崩。所以读 DS18B20 的那一段代码，必须保证中间不被高优先级中断切走，实在不行就临时关一下总中断。
基本的三段操作是复位、写字节、读字节。复位是靠主机拉低总线 480μs 以上，然后释放，器件会在 15~60μs 内把总线拉低作为应答。写一位时，拉低 2μs 后立刻送数据位，写 1 就释放，写 0 就继续拉低到 60μs。读一位时，主机先拉低 2μs 发起读时隙，然后在 15μs 内采样总线电平。
sbit DQ = P3^3; /* 单总线数据脚 */ /* 12T 模式，6MHz，一次循环约 2us，具体次数按编译器优化级别实测校准 */ static void delay_us(unsigned char n) { while (n--) { _nop_(); _nop_(); } } bit ds18b20_reset(void) { bit ack; DQ = 0; delay_us(240); delay_us(240); /* 拉低 > 480us 产生复位脉冲 */ DQ = 1; delay_us(30); /* 等 15~60us 的应答窗口 */ ack = DQ; /* 0 表示有器件应答 */ delay_us(210); /* 补足 480us 释放期 */ return ack; } void ds18b20_write_byte(unsigned char dat) { unsigned char i; for (i = 0; i < 8; i++) { DQ = 0; delay_us(1); /* 拉低 > 1us 起写时隙 */ DQ = dat & 0x01; /* 写 1 释放，写 0 继续保持低 */ delay_us(30); /* 整个时隙 >= 60us */ DQ = 1; dat >>= 1; } } unsigned char ds18b20_read_byte(void) { unsigned char i, dat = 0; for (i = 0; i < 8; i++) { dat >>= 1; DQ = 0; delay_us(1); /* 拉低 > 1us 起读时隙 */ DQ = 1; delay_us(5); /* 15us 内采样有效 */ if (DQ) dat |= 0x80; delay_us(30); /* 补足时隙 */ } return dat; }延时函数的循环次数千万别照抄，它跟你用的编译器优化等级、晶振频率、甚至_nop_()被展开的方式都有关系。正确的做法是打开示波器，量一下复位脉冲的实际宽度，落在 480~960μs 之间才算合格，短了器件不应答，长了也不是不行但会拖慢整体节奏。
温度读取的流程是固定的五步：复位 → 跳过 ROM（0xCC）→ 启动转换（0x44）→ 等转换完成 → 复位 → 跳过 ROM → 读暂存器（0xBE）→ 连读两字节。
#define TEMP_ERR (-1000) /* 无效读数标志 */ int ds18b20_read_temp_x16(void) { unsigned char tl, th; int raw; if (!ds18b20_reset()) return TEMP_ERR; ds18b20_write_byte(0xCC); /* 总线上只有一个器件，跳过 ROM 匹配 */ ds18b20_write_byte(0x44); /* 启动温度转换 */ delay_ms(750); /* 12 位分辨率最坏 750ms */ if (!ds18b20_reset()) return TEMP_ERR; ds18b20_write_byte(0xCC); ds18b20_write_byte(0xBE); /* 读暂存器 */ tl = ds18b20_read_byte(); /* 低字节在前 */ th = ds18b20_read_byte(); /* 高字节在后 */ raw = ((int)th << 8) | tl; /* 二进制补码，单位 1/16℃ */ if (raw == 0x0550) return TEMP_ERR; /* 85℃ 上电默认值，视为无效 */ return raw; }最后这个0x0550的判断非常关键。DS18B20 的暂存器上电默认值是 85℃，如果你在转换还没完成的时候就去读，读回来的就是这 550。很多人调试时看到温度恒定显示 85，第一反应是传感器坏了，其实十有八九是延时不 够或者时序被打断。把它当成无效值直接丢掉，比放进算法里强。
3.2 可控硅触发与时间比例调功
可控硅调功有两种做法：调导通角（相位控制）和调通断时间（周波控制）。相位控制能做到连续调压，但会在电网里注入大量谐波，对 51 这样的慢速控制器来说，过零后精确延时触发也很难做准。所以这套系统选的是周波控制——在一个固定时间窗口内（这里取 1 秒，正好 50 个工频周期、100 个半波），决定让多少个半波导通，剩下的一律关断。
1 秒的窗口相对于水箱的热时间常数（分钟级）非常短，加热丝得到的平均功率几乎是连续的，而可控硅每次都是在过零点附近开关，谐波和浪涌都很小。
实现上需要一个过零检测电路，把市电的每个过零点变成一个脉冲，接到单片机的外部中断 0。中断里数半波，数到设定值就停止触发。
sbit TRIAC = P1^7; /* 可控硅触发输出 */ volatile unsigned char hw_cnt = 0; /* 本窗口内已过的半波数 */ volatile unsigned char hw_duty = 0; /* 本窗口应导通的半波数，0~100 */ void ext0_isr(void) interrupt 0 /* 过零检测，每个半波触发一次 */ { if (hw_cnt < hw_duty) { TRIAC = 1; /* 触发导通 */ delay_us(100); /* 脉冲宽度要大于擎住电流建立时间 */ TRIAC = 0; } if (++hw_cnt >= 100) { hw_cnt = 0; /* 1s 窗口结束，重新计 */ hw_duty = pid_duty; /* 载入算法给出的新占空比 */ } } void triac_set_duty(unsigned char d) { if (d > 100) d = 100; pid_duty = d; /* 只写不读，避免和中断打架 */ }这段代码有两个容易翻车的地方。第一，hw_duty是中断和主循环共享的变量，主循环写、中断读，8 位机上一个字节的写是原子的，所以这里不需要关中断保护，但如果你把占空比扩成 16 位，就必须要EA = 0; ... EA = 1;包一下，否则可能读到高字节是新的、低字节是旧的。第二，触发脉冲的宽度不能太窄，可控硅一旦导通，只要电流超过擎住电流就会自锁，但如果脉冲太短、负载电流又没建立起来，就会触发失败，表现出来就是「占空比给了 100 但水箱不升温」。
3.3 显示与按键：74LS164 串行驱动和四键设置
三只数码管如果直接用 I/O 口驱动，至少要占 11 根线（8 段 + 3 位选）。这套方案用两片 74LS164 级联，把串行数据转成并行段码，只占 P3.0（数据）和 P3.1（时钟）两根线，代价是刷新时要把三位的数据依次推出去。
sbit SDATA = P3^0; /* 74LS164 数据输入 */ sbit SCLK = P3^1; /* 74LS164 时钟 */ static const unsigned char seg_code[10] = { 0xC0, 0xF9, 0xA4, 0xB0, 0x99, 0x92, 0x82, 0xF8, 0x80, 0x90 }; /* 共阳数码管 0~9 段码 */ void display_show(unsigned char val) { unsigned char h, t, u, i; h = val / 100; t = val % 100 / 10; u = val % 10; /* 先推百位，最后推个位，个位会落在最靠近输出的那片 164 上 */ for (i = 0; i < 8; i++) { SDATA = (seg_code[h] >> (7 - i)) & 0x01; SCLK = 0; SCLK = 1; /* 上升沿移入 */ } for (i = 0; i < 8; i++) { SDATA = (seg_code[t] >> (7 - i)) & 0x01; SCLK = 0; SCLK = 1; } for (i = 0; i < 8; i++) { SDATA = (seg_code[u] >> (7 - i)) & 0x01; SCLK = 0; SCLK = 1; } }段码是共阳的，所以 0 亮 1 灭，别拿共阴的表直接套，套上去是一片黑。另外 74LS164 没有输出锁存，移位过程中数码管是能看到中间态的，快速连续刷的时候会有一点点鬼影，正常现象，如果在意就在输出端加一片 74LS273 做锁存。
按键这边，四个键接 P1 口低四位，程序上没什么花活，就是扫描加消抖。但设置逻辑要写清楚：百位键只管百位，按一下加一，到 9 回 0；十位个位同理。因为限定范围是 40~90℃，所以百位实际只有 4~9 这一段有意义，越界的组合直接拒掉。
void key_scan(unsigned char *set_temp) { static unsigned char last_key = 0x0F; unsigned char k = P1 & 0x0F; if (k == last_key) return; /* 无变化直接返回，天然消抖 */ last_key = k; if (k == 0x0F) return; /* 全高说明是松开，不是按下 */ if (k == 0x0E) { /* 启动键 */ start_flag = 1; } else if (k == 0x0D) { *set_temp = bump(*set_temp, 100, 40, 90); } else if (k == 0x0B) { *set_temp = bump(*set_temp, 10, 40, 90); } else if (k == 0x07) { *set_temp = bump(*set_temp, 1, 40, 90); } }这种「记录上次状态、变化才处理」的写法比「延时 10ms 再读一次」省事得多，也不会卡住主循环。缺点是对快速连按不敏感，不过温度设定本来就不需要连按，够用。
4. 模糊 PID 核心：从 E、EC 到 ΔKp、ΔKi、ΔKd 的完整链路
这一章是整套系统的重点。模糊 PID 听着玄，拆开就是三件事：把精确量模糊化、查表推理、再解模糊回精确量。难点不在数学，在于怎么在只有 256 字节 RAM 的 51 上把它跑得又快又稳。
4.1 模糊化：论域、量化因子和隶属函数
模糊控制器的两个输入是误差 E 和误差变化率 EC。误差就是设定温度减实测温度，单位用 1/16℃ 的原始码；误差变化率是本次误差减上次误差，单位同样是 1/16℃，因为是按固定 15 秒周期采的，所以这个差值本身就代表了「每 15 秒变化了多少」。
两个输入都要映射到一个统一的整数论域上。这套系统取的论域是 {-3, -2, -1, 0, 1, 2, 3}，对应七个模糊语言值：NB（负大）、NM（负中）、NS（负小）、ZO（零）、PS（正小）、PM（正中）、PB（正大）。七个值是最常用的档位，再多了规则表会膨胀到 81 条甚至 121 条，51 的 Flash 装得下但查表开销上去了；再少了控制精度不够，0.5℃ 的指标就守不住。
量化因子决定了物理量到论域的缩放比例。误差范围按 ±10℃ 设计，10℃ 对应论域边界 3，所以 Ke = 3/10 = 0.3。误差变化率的范围小得多，水箱在 15 秒内的温度变化撑死 2℃，所以 Kec = 3/2 = 1.5。这两个值不是拍脑袋定的，它们直接决定了控制器的「灵敏度」：Ke 太小，小误差全被压成 ZO，系统反应迟钝；Ke 太大，稍微一点偏差就顶到 PB，输出剧烈跳变。
/* 把 1/16℃ 单位的物理量量化到 -3~3，k_num/k_den 是量化因子 */ static signed char fuzzy_quant(int raw, unsigned char k_num, unsigned char k_den) { long t; t = (long)raw * 10L / 16L; /* 先换成 0.1℃ 单位，保留一位小数 */ t = t * k_num / k_den; /* 再乘量化因子 */ if (t > 3) t = 3; if (t < -3) t = -3; return (signed char)t; }这里全部用整数运算，long在 C51 里是 32 位，一次乘除大概几十个机器周期，15 秒跑一次完全不心疼。中间那步「换成 0.1℃ 单位」是为了保留精度——如果直接拿 1/16℃ 的原始码去除，量化误差会吃掉小偏差段的分辨率，而小偏差段恰恰是稳态精度最吃紧的地方。
隶属函数这块，工程上几乎清一色用三角形。原因很简单：三角形只需要三个点就能描述，存储省、计算简单，而且相邻隶属函数交点处的隶属度是 0.5，符合「一个量同时属于两个语言值」的直觉。高斯型虽然光滑，但要算指数，在 51 上纯属自找麻烦。
实际实现时还有一个更省事的做法：既然论域只有 7 个整数点，而三角形隶属函数在这 7 个点上的取值是固定的，干脆把隶属度表离线算好存进code区，运行时直接查，连乘法都省了。常见做法是把 7×7=49 个点的隶属度打平成一张unsigned char数组，索引就是e*7+ec。
4.2 模糊规则表：49 条规则和三张决策表
规则的形式是「IF E is A AND EC is B THEN ΔKp is C」。对 Kp、Ki、Kd 三个参数各要一张 7×7 的表，总共 147 条规则。手写这么多条不现实，实际都是按工程经验总结出规律再填表。
Kp 的整定规律最好记：误差大时 Kp 要大，加快响应；误差中等时 Kp 要降，防止超调；误差很小时 Kp 又要适当加大，保证稳态精度。Kd 的作用是抑制超调，所以在误差变化率大的地方要加大 Kd，但误差极小的稳态段反而要减小 Kd，否则微分项会把传感器噪声放大成输出抖动。Ki 主要用于消除静差，一般在误差中等偏小时才投入，误差大时积分饱和会把系统拖垮。
/* ΔKp 模糊规则表：行 = E(NB..PB)，列 = EC(NB..PB) */ static const signed char rule_dkp[7][7] = { /* EC: NB NM NS ZO PS PM PB */ /*E NB*/ { PB, PB, PM, PM, PS, ZO, ZO }, /*E NM*/ { PB, PB, PM, PS, PS, ZO, NS }, /*E NS*/ { PM, PM, PM, PS, ZO, NS, NS }, /*E ZO*/ { PM, PM, PS, ZO, NS, NM, NM }, /*E PS*/ { PS, PS, ZO, NS, NS, NM, NM }, /*E PM*/ { PS, ZO, NS, NM, NM, NM, NB }, /*E PB*/ { ZO, ZO, NM, NM, NM, NB, NB } }; /* ΔKd 模糊规则表：误差变化率大的区域加大微分作用 */ static const signed char rule_dkd[7][7] = { /* EC: NB NM NS ZO PS PM PB */ /*E NB*/ { PS, NS, NB, NB, NB, NM, PS }, /*E NM*/ { PS, NS, NB, NM, NM, NS, ZO }, /*E NS*/ { ZO, NS, NM, NM, NS, NS, ZO }, /*E ZO*/ { ZO, NS, NS, NS, NS, NS, ZO }, /*E PS*/ { ZO, ZO, ZO, ZO, ZO, ZO, ZO }, /*E PM*/ { PB, NS, PS, PS, PS, PS, PB }, /*E PB*/ { PB, PM, PM, PM, PS, PS, PB } };Ki 那张表结构一样，中心区（E 和 EC 都在 ZO 附近）取正值消除静差，外围取零或者很小的值防止积分饱和，这里就不重复贴了。
表怎么存很讲究。signed char一行 7 个字节，一张表 49 字节，三张 147 字节。AT89S52 有 8K Flash，放这点数据毫无压力，但一定要加code关键字让它待在 Flash 里。如果不加，Keil 会默认分配到片内 RAM，256 字节的 RAM 一下子被吃掉一半多，栈空间立刻告急，程序跑着跑着就飞了——这是 51 上最经典的翻车方式之一。
4.3 去模糊与输出：查表法替代重心计算
标准流程里，去模糊要用重心法把所有激活规则的隶属度加权平均，计算量不小。但在本系统这种只有 7 个论域点的情况下，有个更实用的简化：既然 E 和 EC 都已经量化成整数，那查出来的规则结果直接就是 -3~3 的整数值，把它乘上比例因子就得到实际的参数增量。
这个做法在教科书上叫「单点模糊化 + 加权平均去模糊」的特例，本质上忽略了相邻模糊集的交叠。代价是输出在论域边界处会有台阶，但温度控制对这种小台阶完全不敏感，水箱的热惯性会把它们全都抹平。
/* 基准 PID 参数，Q8 定点：实际值 = 变量 / 256 */ static int kp_base = 512; /* 2.0 */ static int ki_base = 26; /* 0.1 */ static int kd_base = 1280; /* 5.0 */ /* 比例因子：规则表输出 1 个单位，对应基准值的多少比例 */ #define KU_KP_Q8 64 /* 1/4 基准 */ #define KU_KI_Q8 13 /* 1/20 基准 */ #define KU_KD_Q8 128 /* 1/2 基准 */ int fuzzy_pid(int set_x16, int fb_x16) { static int e1 = 0, e2 = 0; /* 前两次误差 */ static int u = 0; /* 累计占空比 0~100 */ int e, ec, du; signed char fe, fce, fkp, fki, fkd; int kp, ki, kd; long t; e = set_x16 - fb_x16; /* 单位 1/16℃ */ ec = e - e1; fe = fuzzy_quant(e, 3, 10); /* Ke = 0.3 */ fce = fuzzy_quant(ec, 15, 10); /* Kec = 1.5 */ fkp = rule_dkp[fe + 3][fce + 3]; fki = rule_dki[fe + 3][fce + 3]; fkd = rule_dkd[fe + 3][fce + 3]; kp = kp_base + (fkp * KU_KP_Q8); ki = ki_base + (fki * KU_KI_Q8); kd = kd_base + (fkd * KU_KD_Q8); if (kp < 64) kp = 64; /* 兜底，防止参数被规则拉到零 */ if (ki < 0) ki = 0; if (kd < 0) kd = 0; /* 增量式 PID：Δu = Kp*(e-e1) + Ki*e + Kd*(e-2*e1+e2) */ t = (long)kp * (e - e1) + (long)ki * e + (long)kd * (e - 2L * e1 + e2); du = (int)(t >> 8); /* 还原 Q8，单位仍是 1/16℃ */ du = du / 16; /* 换算成占空比百分点 */ u += du; if (u > 100) u = 100; /* 输出限幅，抗积分饱和第一道闸 */ if (u < 0) u = 0; e2 = e1; e1 = e; return u; }几点说明。第一，rule_dki表这里没贴，结构跟另外两张一样，实际项目里补上就行。第二，du/16那个换算系数是「1℃ 偏差对应 1% 占空比」的经验起点，实际必须现场标定——水箱容积、加热丝功率、环境温度都会影响这个比例，大了会振荡，小了升不上去。第三，输出限幅u在 0~100 之间，这是最朴素也最有效的抗积分饱和手段，比什么复杂的 back-calculation 都实在。
还有一点容易被忽略：e1、e2用的是「误差」而不是「实测值」。这在设定温度不变时没区别，但一旦用户按了按键改设定值，用误差做历史量会产生一个巨大的阶跃，输出瞬间打满。稳妥做法是在设定值变化时清空e1、e2和u，让系统重新起步。
4.4 定时调度：中断嵌套里谁先谁后
这套系统的中断源不少：定时器 0 管采样节拍，外部中断 0 管过零检测，串口（如果加了）管调试输出。51 只有两级中断优先级，必须排清楚。
过零检测是硬实时的，晚了几毫秒这一个半波就浪费了，而且过零点就那么一瞬间，所以它必须是高优先级。定时器 0 的 50ms 中断晚一点无所谓，设成低优先级即可。这样即便 15 秒到点的那一刻正好撞上过零中断，也只是推迟几十微秒执行，对温控毫无影响。
void interrupt_init(void) { IT0 = 1; /* 外部中断 0 下降沿触发 */ EX0 = 1; /* 开外部中断 0 */ PX0 = 1; /* 过零检测设为高优先级 */ ET0 = 1; /* 定时器 0 中断 */ PT0 = 0; /* 采样节拍设为低优先级 */ EA = 1; }优先级配好之后，还有一个隐性的坑：过零中断里调用了delay_us做脉冲展宽，这段时间 CPU 是被占住的。如果此时定时器 0 也到点，它只能等着。50ms 的中断延迟几十微秒完全没问题，但如果你把过零中断的延时改长，或者在过零中断里加了别的事情，累积起来就会让 15 秒的采样周期偏长。测量方法很简单：在定时器 0 中断里翻转一个空闲 I/O 口，用示波器看方波周期，是不是准确的 50ms 一目了然。
5. 常见问题与排查：恒温箱跑不起来多半是这几个坑
这套东西的电路和代码都不算复杂，但真动手做的时候，翻车点相当集中。下面这几条是我见过频率最高的，按「现象 → 原因 → 解决」拆开写。
5.1 温度恒定显示 85℃ 或者数值乱跳
现象：上电后数码管一直显示 85，或者读数在几个值之间毫无规律地跳。
原因：85℃ 是 DS18B20 暂存器的上电默认值，读到它就说明温度转换根本没完成。要么是delay_ms(750)不够（12 位分辨率最坏情况就是 750ms），要么是这 750ms 里被中断打断导致总线时序乱了。数值乱跳则多半是 DQ 线上少了 4.7kΩ 上拉电阻，开漏输出拉不高，读回来的位全是随机的。
解决：先量复位脉冲宽度，确认在 480~960μs；再把读温度的整个流程用EA=0包起来，避免中断插入；最后检查 DQ 上的上拉电阻，长线缆（超过 1 米）要适当减小阻值或者改用外部供电模式。程序里加一道raw == 0x0550的判据，把无效值挡在算法外面，比事后调试省事得多。
5.2 占空比给到 100%，水箱还是不升温
现象：算法输出正常，数码管显示占空比拉满，但水箱温度纹丝不动。
原因：可控硅没真正导通。常见的有三种：光耦选错型号，用了普通三极管输出型去驱可控硅，触发极性对不上；触发脉冲太窄，可控硅还没来得及建立擎住电流脉冲就撤了；加热丝串在回路里但可控硅的 A1/A2 接反，或者门极电阻取值过大导致触发电流不足。
解决：先用万用表量加热丝两端有没有 220V，没有就往上查可控硅；有电压但不热，说明加热丝本身断了。确认是可控硅问题的话，把触发脉冲展宽到 100μs 以上，光耦换成双向可控硅输出型，门极串一个 300Ω 左右的限流电阻。还有一个隐蔽的点：如果负载是感性或者带变压器，必须加 RC 吸收网络，否则关断瞬间的 dv/dt 会让可控硅误导通，表现出来就是「关了还在加热」。
5.3 温度冲过设定值一大截才回落
现象：设定 60℃，实际冲到 68℃ 才开始降，然后在 55~65 之间来回荡。
原因：积分饱和。启动瞬间误差很大，积分项一路累加，等温度快到设定值时积分已经攒了一大堆，输出退不下来。水箱这种大惯性对象尤其明显，热量灌进去要等好几分钟才能体现在温度上，控制器的「刹车」永远慢半拍。
解决：三道闸一起上。第一道是输出限幅，把累计输出卡在 0~100；第二道是积分分离，误差绝对值大于某个阈值（比如 5℃）时直接把 Ki 置零，只用 PD 快速逼近；第三道是设定值变化时清空积分累计和历史误差。这三条加完，超调量通常能从 8℃ 压到 2℃ 以内。另外，模糊规则表里 E 为 PB、EC 为 PB 那一格如果给的是负的 ΔKp，也能帮上忙——误差大且还在快速扩大时，反而要减小比例增益，防止输出打满。
5.4 数码管有鬼影、按键偶尔失灵
现象：数码管上不该亮的段微微发亮，或者按键要按好几下才有反应。
原因：鬼影是 74LS164 没有输出锁存，移位过程中中间态被数码管看到了；按键失灵多半是扫描周期太长——如果按键扫描放在 15 秒的采样分支里，那按一下得等 15 秒才有反应。
解决：鬼影可以加一片 74LS273 做输出锁存，或者干脆提高移位时钟频率让中间态短到看不见。按键必须放在主循环的最外层，不受sample_flag约束，这样任何时候按都有响应。另外，P1 口低四位做按键输入时要先写 1，让内部上拉生效，否则读回来永远是 0。
5.5 程序烧进去跑几分钟就复位
现象：上电正常，跑一阵子突然重启，或者显示乱码。
原因：十有八九是 RAM 溢出。51 的片内 RAM 只有 256 字节，其中高 128 字节还和 SFR 重叠，实际可用更少。如果规则表忘了加code关键字，147 字节的表直接把 RAM 吃掉大半，加上栈和局部变量，一调用深一点的函数就踩栈。另一个可能是看门狗没喂——AT89S52 片内有看门狗，如果使能了但主循环里没定时清零，就会周期性复位。
解决：所有常量表一律加code，编译完看一眼 Keil 的.m51文件，确认 DATA 段占用别超过 100 字节。局部变量尽量少用大的数组，能用全局就用全局。看门狗要么不用，要么就在主循环里老老实实喂。
6. 二次开发：把算法层剥出来，换平台只改三处
这套代码要拿去做二次开发，最值得动的地方不是硬件，是分层。原始代码把采样、算法、输出全揉在一个main里，改一个参数要翻三个文件。我的做法是把模糊 PID 抽成一个纯函数——输入设定值和实测值，输出占空比，中间不碰任何全局变量和硬件寄存器。抽干净之后，你会发现移植成本低得吓人。
| 改造方向 | 改哪里 | 怎么验证改对了 |
|---|---|---|
| 换主控到 STM32 | 只替换delay_us、定时器和 GPIO 三层封装 | 用逻辑分析仪对比 DS18B20 时序波形 |
| 加自整定 | 在fuzzy_pid外面套一层，稳态时微调kp_base | 阶跃响应超调量应逐步收敛 |
| 换执行机构为固态继电器 | 改triac_set_duty，把 0~100 映射成 PWM 占空比 | 用示波器看 PWM 频率和占空比线性度 |
| 增加串口上位机 | 主循环里加一帧上报，格式S:60.0 F:59.8 U:45 | 上位机曲线和数码管显示一致 |
| 多段温区曲线 | 把set_temp改成查表，按时间推进索引 | 每段切换点温度无突跳 |
移植的时候有个技巧：把fuzzy_pid里的long运算全部改成int32_t，把signed char改成int8_t，代码一行不用动就能在 ARM 上编译通过。真正需要重写的只有底层三个函数——微秒延时、毫秒延时、以及那个读温度的时序函数。DS18B20 的时序在 STM32 上反而更好做，因为主频高、延时精度好，用定时器加 DMA 都能实现，不会再出现「延时函数被优化掉」这种玄学问题。
还有一处值得单独拎出来说：fuzzy_quant里的量化因子目前是硬编码的3/10和15/10。这两个数在整定阶段是一定要反复调的，建议改成两个全局变量，通过串口在线下发，边看曲线边改，比每次改完重新编译烧录快十倍。同理，kp_base、ki_base、kd_base也应该可在线修改，配合一个简单的「保存到 EEPROM」功能，调好的参数断电不丢。
我自己踩过最深的一个坑是：早期版本为了图省事，把rule_dkp表直接用二维数组下标[fe][fce]访问，fe和fce是signed char类型，范围 -3~3。程序在误差为正的时候一切正常，一旦温度超过设定值变成负误差，下标就跑到数组前面去了，读出来的是栈上的随机数据，表现是「降温比升温猛得多」，查了整整两天才定位到。从那以后我每次写查表代码，都把下标偏移量单拎出来加一句断言，或者干脆在量化函数返回前就把范围钳死，绝不让负下标有机会流到数组访问那一层。
这个恒温箱项目本身不大，硬件成本也就百来块，但它把「传感器时序、执行器驱动、模糊推理、定点运算、中断调度」这条链路完整走了一遍，拿来练手比自己臆想要扎实得多。希望帮到你。
本文还有配套的精品资源，点击获取
