简介吉林大学通信工程学院自动化专业大三的《计算机控制系统》课件面向自动化、电气与控制工程专业学生及从业者系统讲解计算机控制系统的核心组成传感器、数据采集、控制器、执行机构、反馈回路与人机交互、常见类型开环/闭环、PLC、DCS、现场总线与关键知识点经典/现代控制理论、数字信号处理、实时操作系统、工业网络通信内容由浅入深既可用于课堂同步复习也可作为备考与入职培训的参考资料。资源为zip压缩包共116个文件以85个PPT课件为主体辅以10个DOC试题与答案、9个MATLAB的m脚本、6个fig图形和6个gif演示动画整体大小仅14.54MB便于快速下载与离线学习。目前已有185人学习下载除课堂讲义外资源还包含试题参考答案、仿真源码和图形动画可帮助读者对照理解闭环反馈、PID控制、数字信号处理等抽象概念并通过动手实践验证控制算法、梳理考试重点。1. 计算机控制系统这门课的资源到底该先看哪块计算机控制系统的资源包信息密度比想象中高得多离散化建模、数字PID、仿真调试、实物联调每一块单独拿出来都能写几十页。下载之后先别急着翻代码我的经验是先想清楚一个问题——你要控制的那个连续对象计算机到底是以什么节奏、按什么模型在看见它。这套资源适合正在做计算机控制课程设计的学生、第一次把连续对象搬到单片机上的入门工程师以及想把手里的理论补上工程细节的从业者。把采样、离散化、PID和联调这条线捋直后面的代码每个字都看得懂。2. 离散化建模采样周期、Z变换与差分方程怎么落2.1 Z变换不是玄学它只回答离散后长什么样计算机控制里大部分翻车根源都在同一个地方拿连续系统的思路去写离散的代码。被控对象的原始模型是连续传递函数G(s)而计算机每个采样周期才看一次对象中间这段时间对象在干什么它不知道。Z变换要解决的问题就是把这个连续模型翻译成计算机真正能执行的形式。这里有个关键器件叫零阶保持器ZOH它假设计算机输出的控制量在两次更新之间保持不变。加了ZOH之后离散模型和连续模型之间的关系是G(z) (1 - z⁻¹) · Z{G(s) / s}这个式子手工推起来非常痛苦实际项目里没人这么干。我一般直接调用现成的离散化工具常用仿真软件里有c2d函数Python生态里有control库几行代码就能拿到离散化结果。但理解ZOH的含义仍然重要——它决定了采样周期一变整个离散模型都要重新算不能只改代码里的一个延时参数。2.2 采样周期怎么定香农定理只是下限工程约束才是真门槛采样周期的选择是计算机控制系统设计里最容易被低估的一步。香农定理说采样频率要大于信号最高频率的两倍但这只是理论下限不是工程推荐值。工程上真正约束采样周期的是三件事上升时间、执行机构响应、以及量化计算。上升时间约束要求采样周期满足Ts ≤ tr/10其中tr是被控对象开环阶跃响应的上升时间这样闭环控制才有足够的纠偏密度。执行机构约束指的是控制量更新频率得和功率器件能响应的频率匹配更新太频繁执行器根本没反应。量化约束更隐蔽——采样周期太小控制器每次只能输出很小的控制增量量化误差直接把控制量吃掉系统表现为低速段跳动。我一般按被控对象的时间常数初选采样周期然后在这个表的基础上微调| 控制对象 | 典型时间常数 | 采样周期初选 | | 温度 | 分钟级 | 110 s | | 液位 | 秒级 | 0.52 s | | 电机转速 | 几十几百 ms | 110 ms | | 电流环 | ms 级 | 0.11 ms |初选之后把采样周期代入离散化工具再看一次开环Bode图如果相位裕度比连续模型掉得太多说明采样周期偏大需要减小。这一步是验证采样周期最快的方式。提示采样周期一旦修改前面算好的离散模型系数全部要重算不要只改定时器而不重新离散化。2.3 从传递函数到差分方程三步拿到能写进单片机的形式拿到采样周期后离散化的完整流程是这样的# 连续被控对象G(s) 1 / (s^2 2*s 3) import control as ct G ct.tf([1], [1, 2, 3]) # 连续传递函数分母系数按降幂排列 Ts 0.01 # 采样周期由上一节的原则选定 Gd ct.c2d(G, Ts, methodzoh) # 用零阶保持器离散化 print(Gd) # 直接打印出离散传递函数这段代码里ct.tf的第一个参数是分子系数第二个是分母系数。c2d第三个参数methodzoh指定零阶保持器这是和物理实现最匹配的离散化方式。Ts就是上一步选定的采样周期。打印出的离散传递函数有分子分母两行系数它们直接对应差分方程y(k) b₀·u(k) b₁·u(k-1) b₂·u(k-2) - a₁·y(k-1) - a₂·y(k-2)把系数写进单片机代码时我习惯先核对两件事第一分母首项必须是1如果不是说明工具做了归一化你得搞清楚它除以了什么第二系数符号不能抄错尤其分母项经常带负号。这个细节错过一次仿真和实物对不上能查一上午。2.4 开环验证离散模型和连续模型必须对得上离散化之后别急着进闭环。先在仿真里把连续模型和离散模型的阶跃响应叠在一起看曲线基本重合说明采样周期选得没问题如果离散模型有明显滞后、发散或者超调变形那就是Ts太大回去重选。这一步五分钟就能结束但能省掉后面所有联调阶段的玄学问题。开环验证通过才算真正拿到了计算机能执行的被控对象模型。3. 数字PID落地位置式、增量式与抗饱和三段代码3.1 位置式PID结构体定义、公式与输出限幅数字PID是所有计算机控制系统课程设计的核心无论题目是温度、液位还是电机转速最终都要落到这一段代码。位置式PID是最直观的写法它直接计算当前控制量u(k) Kp·e(k) Ki·Σe(i)·dt Kd·(e(k) - e(k-1))/dt对应代码typedef struct { float Kp, Ki, Kd; // 三个增益整定阶段会频繁修改 float setpoint; // 目标值 float integral; // 积分累计 float prev_error; // 上一次误差 float out_min, out_max; // 输出限幅 } PID_t; float pid_position_update(PID_t *pid, float measurement, float dt) { float error pid-setpoint - measurement; pid-integral error * dt; // 积分项要乘采样时间才能量纲正确 float derivative (error - pid-prev_error) / dt; float output pid-Kp * error pid-Ki * pid-integral pid-Kd * derivative; if (output pid-out_max) output pid-out_max; if (output pid-out_min) output pid-out_min; pid-prev_error error; return output; }注意两个细节。积分项里error乘以dt是为了把积分从累加误差个数变成误差对时间的积分这样Ki的量纲才是1/s和连续PID公式对得上。输出限幅放在最后但只限幅输出而不处理积分后面会引出积分饱和问题。dt必须是真实采样间隔不能把定时器的名义周期硬填进去实际有抖动时用上一次实测间隔更稳。注意PID结构体里的积分项是长期累加量长时间运行要留意浮点精度必要时定期把积分拉回零点附近。3.2 增量式PID为什么工程现场更常用增量式PID不直接给控制量而是给控制量的增量Δu(k) Kp·(e(k) - e(k-1)) Ki·e(k) Kd·(e(k) - 2·e(k-1) e(k-2))typedef struct { float Kp, Ki, Kd; float setpoint; float e_prev, e_prev2; // 最近两次的历史误差 float output; // 当前控制量 float out_min, out_max; } IncPID_t; float pid_incremental_update(IncPID_t *pid, float measurement, float dt) { float error pid-setpoint - measurement; float delta pid-Kp * (error - pid-e_prev) pid-Ki * error * dt pid-Kd * (error - 2.0f * pid-e_prev pid-e_prev2) / dt; pid-output delta; if (pid-output pid-out_max) pid-output pid-out_max; if (pid-output pid-out_min) pid-output pid-out_min; pid-e_prev2 pid-e_prev; pid-e_prev error; return pid-output; }增量式在带执行机构保持功能的系统里优势明显比如步进电机、带自锁的阀门断电或通信中断时执行器停在最后位置不会猛冲。位置式输出直接给控制量一旦工程现场临时改限幅或者重启可能把执行器突然推到极限。增量式的代价是误差被积分进output里长期运行要确认浮点精度够也可以定期把output强制拉回零点附近重新积累。3.3 积分饱和与微分滤波工程上必加的两个补丁教科书上的PID写到3.1就该停了但实物联调时积分饱和和微分噪声两个问题几乎必然出现。积分饱和的典型场景输出已经限幅误差仍然存在积分还在涨等误差反向积分要很久才能退回来系统表现为反应迟钝、过冲明显。// 抗积分饱和输出未限幅时才允许积分累加 float output_pre pid-Kp * error pid-Ki * pid-integral pid-Kd * derivative; if (output_pre pid-out_min output_pre pid-out_max) { pid-integral error * dt; } // 微分项接一阶低通滤波压掉高频噪声 float alpha dt / (dt Tf); float deriv_raw (error - pid-prev_error) / dt; pid-deriv_filtered alpha * deriv_raw (1.0f - alpha) * pid-deriv_filtered;这里的条件式积分是最常用的抗饱和策略输出饱和时冻结积分等误差反向再放开。微分滤波的时间常数Tf一般取采样周期的1/10到1/5太小没效果太大会把有用的微分作用也滤掉。这两个补丁教科书里常常略过但少任何一个实物调试都会让你怀疑人生。3.4 整定顺序从P走到PID每一步都要留记录参数整定我固定按三步走。第一步只留比例项Ki和Kd全部置零从小到大加Kp直到系统出现等幅振荡记下这个临界Kp第二步取临界值的0.5到0.6倍作为工作Kp然后加积分从很小的Ki开始看到稳态误差缓慢消失为止第三步才碰微分用于压超调。每次改参数都要记录波形变化不要凭感觉连调三个参数那是典型越急越乱的操作方式。4. 调试排查采样、量化与控制量五个高频坑这一章是血泪经验的集中区。计算机控制系统的仿真波形再漂亮实物联调都会碰到下面几个问题基本每个实验室都会遇到一遍。4.1 输出剧烈抖动执行器高频振动现象PID一投入闭环控制量毛刺特别多电机或阀门跟着高频颤抖声音发闷发振。原因微分项把测量噪声放大了。示波器上看着干净的信号在ADC采样里其实带着几毫伏波动微分项对这类波动极其敏感增益一高控制量全是毛刺。解决先在微分项上加一阶低通滤波滤波时间常数从采样周期的1/10起步再检查反馈通道要不要加硬件RC滤波。注意RC滤波的时间常数要远小于采样周期否则采到的信号本身被平均了整个系统会变迟钝。4.2 稳态误差消不掉加大Ki反而震荡现象误差长期存在一个固定值加大Ki之后误差没消失系统反而开始周期震荡。原因输出已经限幅误差还在继续积分这是典型的积分饱和。更隐蔽的情况是执行机构存在死区输出小于某个阈值时执行器根本不动积分项一直累加等累加跨过死区输出猛然跳变形成极限环。解决先加条件式抗积分饱和再针对死区做补偿比如PWM占空比小于最小可动值时直接跳到最小可动值或者做线性化映射。死区阈值必须实测直接抄说明书标称值往往不准。4.3 系统低频爬行动作一卡一卡现象给定值是小幅阶跃时系统很久才动一下动一下又停下像爬行一样。原因采样周期和PWM或执行机构驱动频率匹配不当。采样周期太短每次控制增量小于执行机构分辨率或者PWM分辨率太低控制量变化被量化成几档中间档缺失。解决适当增大采样周期让每次输出变化跨过执行机构分辨率PWM位数从8位提升到10位或12位。这个问题跟PID参数无关是量化和执行机构分辨率的匹配问题改Kp和Ki都没用。4.4 上电瞬间执行器猛冲一下现象系统一上电PID还没正常工作执行器先猛冲一步然后又回来。原因初始化时历史误差全是零而上电瞬间测量值和目标值差距很大第一次调用PID相当于输入一个大阶跃微分项瞬间爆表。如果没置位历史误差输出会异常大。解决初始化时把prev_error和积分项设成当前时刻的实测值差不要从零开始更稳妥的做法是上电先开环输出一段时间确认反馈正常后再切入闭环。这个软起步逻辑是实物系统里必写的保护。4.5 仿真曲线漂亮实物完全对不上现象同一套参数仿真里超调5%实物超调30%调节时间差好几倍。原因仿真模型里的参数来自文档铭牌不是实测值。电机时间常数、阀门死区、传感器响应延迟实际值和标称值通常差20%到50%。解决联调前先做一次开环阶跃辨识记录输入和响应用辨识出来的结果更新模型参数。手工拟合精度不够就做最小二乘。这是一次投入长期受益的步骤做完之后参数整定会顺很多。5. 仿真到实物信号调理、执行机构与联调顺序5.1 输出通道占空比映射、死区补偿与PWM频率仿真里控制量是-1到1的小数实物里必须变成执行器真正吃得下的电平量。常见输出链是控制量 → 归一化映射到PWM占空比或D/A电压 → 功率放大 → 执行机构。这条链有两个细节最容易出事映射方向和控制量正负号的处理。我一般先单独写一段归一化映射代码跑通再进PID// 控制量 u 范围 -1~1输出 PWM 占空比 0~100 // pwm_min 是实测得到的最小启动占空比 float pwm_deadband_map(float u, float pwm_min) { float duty fabsf(u) * 100.0f; // 取绝对值换算成占空比 if (duty pwm_min) return 0.0f; // 低于启动阈值直接关断 return duty; }pwm_min必须是实测值不是说明书标称值。比如某型号电机标注启动电压1V但实物到30%占空比才真正转得动。方向由H桥的方向引脚控制不要在占空比里搞正负两套解释。PWM频率同样有讲究频率太低电机啸叫明显动作一卡一卡频率太高开关损耗大驱动器发热。常见做法按执行机构电气时间常数选电机驱动典型值在1020kHz区间慢对象如阀门和加热器可以降到几百赫兹。还有一个容易被忽略的坑是电源共地。PWM输出和ADC采样如果共用一根地线功率级开关瞬间会把几十毫伏噪声灌进采样电路轻则读数漂移重则系统震荡。常见做法是驱动侧单独供电采样电路的地单点接回主控地开关噪声立刻下降一个数量级。这个不是玄学实测过很多次。5.2 反馈通道传感器调理、ADC量程与过采样平均反馈通道的问题比输出通道更隐蔽因为输出错了执行器会动反馈错了系统只是不对难定位。传感器信号要先调理到ADC量程范围内比如ADC满量程3.3V传感器输出01V就该放大三倍再用否则有效分辨率只有满量程的三分之一。ADC过采样取平均是最常用的抗噪手段// 连续采样 8 次取平均等效增加约 1.5 位有效分辨率 uint16_t adc_read_filtered(uint8_t channel) { uint32_t sum 0; for (uint8_t i 0; i 8; i) { sum adc_single_read(channel); // 逐次读取并累加 } return (uint16_t)(sum 3); // 右移三位等于除以 8 }这个做法的前提是信号里确实混有噪声纯直流信号平均多少次都没用。平均次数一般取4到16次再多就明显拖慢采样周期。采样周期是第2章定死的约束不能因为软件滤波随意加长。如果硬件已经加了RC低通过采样次数可以降到4次。5.3 联调顺序先开环后闭环方向必须第一个确认实物联调最忌讳一上来就全套PID闭环我固定按五个阶段推进。第一阶段固定占空比开环输出确认执行机构方向、速度或变化率符合预期第二阶段直接读反馈手动拨动执行器确认传感器读数方向、量程和噪声水平第三阶段只加比例项给小幅阶跃确认响应方向正确第四阶段逐步加积分消稳态误差第五阶段才考虑微分项压超调。开环阶段最容易抓到的错误就是方向反了。方向反了不是改代码就能绕过去的必须在闭环之前抓出来。曾有同学没做开环方向测试直接闭环Kp一加就正反馈发散输出直接顶到限幅差点烧了驱动器。开环三秒钟就能发现的问题闭环之后要多花一小时定位。每走一个阶段都把时间戳、测量值、输出值通过串口记录成CSV。没有记录后面出了问题全凭感觉猜这是最消耗时间的调试方式。串口打印必须带时间戳否则数据点对不上号等于白打。6. 验证技巧十分钟跑完一次阶跃响应测试6.1 测试准备与数据记录不管前面调得多顺最后一定要跑一次完整的阶跃响应测试这是判断系统达没达标的唯一客观手段。做法系统稳定在某个工作点后把目标值改变一个适中幅度一般取满量程的10%到20%记录测量值和输出值直到系统重新稳定。记录频率至少是采样周期的2到3倍否则曲线算不出超调。数据存成CSV后用脚本算指标比人眼估靠谱得多import numpy as np data np.loadtxt(step_log.csv, delimiter,, skiprows1) t, y data[:, 0], data[:, 1] y0 np.mean(y[:20]) # 阶跃前稳态 yss np.mean(y[-50:]) # 阶跃后稳态 ymax np.max(y) overshoot (ymax - yss) / (yss - y0) * 100 # 超调百分比 idx np.where(np.abs(y - yss) 0.02 * (yss - y0))[0] settle t[idx[0]] if idx.size 0 else float(nan) # 2%调节时间 print(f超调量{overshoot:.1f}% 调节时间{settle:.2f}s)6.2 读曲线回答三个问题超调量大于预期优先减小Kp或加微分压超调调节时间太长看是积分太弱还是执行机构饱和稳态误差不为零查反馈通道的偏置和死区。这三个问题对应三个参数方向一次只改一个改完重新跑一次阶跃测试这才算一轮完整的闭环迭代。这套验收习惯是我在一次实际项目的教训里养成的。当时调温控系统调了两天凭感觉觉得差不多了结果换了台设备所有问题全冒出来超调30%稳态误差两度根本交不了差。从那以后我每次做计算机控制系统联调都强制自己先跑一遍阶跃测试数据落地、指标量化再谈改参数。这套流程和上面的代码希望能在你的课程设计或项目里帮你少走几个坑。本文还有配套的精品资源点击获取