PID控制实验报告
pid控制实验报告
pid控制实验报告实验报告：PID控制一、实验目的通过本实验，我们的目的是深入了解PID（比例、积分、微分）控制算法，理解其在实际控制中的应用，掌握PID参数的调整方法。
二、实验原理PID控制是依据被控对象的误差（偏差）与时间的积分、微分关系来确定控制器输出的控制方式。
具体来说，PID控制器输出的控制量=Kp*（当前误差+上次误差*dt+所有误差的积分），其中Kp、Ki和Kd分别为比例系数、积分系数和微分系数。
它通过对偏差的补偿，使得被控对象能够在振荡绕过设定值、稳定达到设定值的过程中快速、准确定位设定值。
三、实验设备本实验采用的设备为PID控制器、液晶显示屏、电压控制电机和传感器。
四、实验步骤1. 首先，我们需要将系统设为手动调节状态，关闭控制器。
2. 然后，我们将传感器和记录仪建立起连接。
3. 将系统调整为自动控制状态，让控制器自行计算控制量、作出相应控制。
4. 调整PID控制器的Kp系数，以调整控制精度。
5. 调整PID控制器的Ki系数，以调整控制的灵敏度。
6. 调整PID控制器的Kd系数，以调整控制器的稳定性。
7. 最终完成调整后，我们可以用振荡器数据展示出来实验结果。
五、实验结果在完成调整后，我们得出的控制器输出的控制量稳定在理论值附近，在控制精度与控制的灵敏度达到较好平衡的情况下，控制器的稳定性得到了保证。
实验结果具有较好指导意义。
六、结论本实验通过掌握PID控制算法的实际应用方法，以及对参数的合理设置为基础，完成了对PID控制器各参数调整技巧的掌握，极大地丰富了实验基础技能。
同时，实验结果为之后的实际应用提供了参考，有着极其重要的现实意义。
《自动控制原理》自动控制PID实验报告
《自动控制原理》自动控制PID实验报告课程名称自动控制原理实验类型：实验项目名称：自动控制PID一、实验目的和要求1、学习并掌握利用MATLAB 编程平台进行控制系统复数域和频率域仿真的方法。
2、通过仿真实验研究并总结PID 控制规律及参数对系统特性影响的规律。
3、实验研究并总结PID 控制规律及参数对系统根轨迹、频率特性影响的规律，并总结系统特定性能指标下根据根轨迹图、频率响应图选择PID 控制规律和参数的规则。
二、实验内容和原理一）任务设计如图所示系统，进行实验及仿真程序，研究在控制器分别采用比例（P）、比例积分（PI）、比例微分（PD）及比例积分微分（PID）控制规律和控制器参数（Kp、Ki、Kd）不同取值时，控制系统根轨迹和阶跃响应的变化，总结pid 控制规律及参数变化对系统性能、系统根轨迹、系统阶跃响应影响的规律。
具体实验容如下：1、比例（P）控制，设计参数Kp 使得系统处于过阻尼、临界阻尼、欠阻尼三种状态，并在根轨迹图上选择三种阻尼情况的Kp 值，同时绘制对应的阶跃响应曲线，确定三种情况下系统性能指标随参数Kp 的变化情况。
总结比例（P）控制的规律。
2、比例积分（PI）控制，设计参数Kp、Ki 使得由控制器引入的开环零点分别处于：1）被控对象两个极点的左侧；2）被控对象两个极点之间；3）被控对象两个极点的右侧（不进入右半平面）。
分别绘制三种情况下的根轨迹图，在根轨迹图上确定主导极点及控制器的相应参数；通过绘制对应的系统阶跃响应曲线，确定三种情况下系统性能指标随参数Kp 和Ki 的变化情况。
总结比例积分（PI）控制的规律。
3、比例微分（PD）控制，设计参数Kp、Kd 使得由控制器引入的开环零点分别处于：1）被控对象两个极点的左侧；2）被控对象两个极点之间；66 3）被控对象两个极点的右侧（不进入右半平面）。
分别绘制三种情况下的根轨迹图，在根轨迹图上确定控制器的相应参数；通过绘制对应的系统阶跃响应曲线，确定三种情况下系统性能指标随参数Kp 和Kd 的变化情况。
pid控制实验报告
pid控制实验报告PID控制实验报告引言PID控制是一种常用的控制算法，广泛应用于工业自动化系统中。
本实验旨在通过实际的PID控制实验，验证PID控制算法的效果和优势，并对PID控制的原理、参数调节方法等进行探讨和分析。
一、实验目的本次实验的目的是通过一个简单的温度控制系统，使用PID控制算法来实现温度的稳定控制。
通过实验，验证PID控制算法的有效性和优越性，掌握PID控制的基本原理和参数调节方法。
二、实验设备和原理本实验所用的设备为一个温度控制系统，包括一个温度传感器、一个加热器和一个控制器。
温度传感器用于实时检测环境温度，加热器用于调节环境温度，控制器用于实现PID控制算法。
PID控制算法是基于误差的反馈控制算法，其主要原理是通过不断地调整控制器的输出信号，使得系统的实际输出与期望输出之间的误差最小化。
PID控制算法由比例控制、积分控制和微分控制三部分组成。
比例控制通过比例系数调整控制器的输出信号与误差的线性关系；积分控制通过积分系数调整控制器的输出信号与误差的积分关系；微分控制通过微分系数调整控制器的输出信号与误差的微分关系。
通过合理调节这三个系数，可以实现对系统的精确控制。
三、实验步骤1. 搭建温度控制系统：将温度传感器、加热器和控制器连接在一起，确保信号传输的正常。
2. 设置期望温度：根据实验要求，设置一个期望的温度作为控制目标。
3. 调节PID参数：根据实验的具体要求和系统的特性，调节PID控制器的比例系数、积分系数和微分系数，使得系统的响应速度和稳定性达到最佳状态。
4. 开始实验：启动温度控制系统，观察实际温度与期望温度的变化情况，记录实验数据。
5. 数据分析：根据实验数据，分析PID控制算法的效果和优势，总结实验结果。
四、实验结果与讨论通过实验，我们得到了一系列的实验数据。
根据这些数据，我们可以进行进一步的分析和讨论。
首先，我们观察到在PID控制下，温度的稳定性得到了显著的提高。
PID自动控制系统参数整定实验报告
T13. PID自动控制系统参数整定（化工仪表与自动化，指导教师：卢红梅）实验一：一阶单容上水箱对象特性测试实验实验二：上水箱液位PID整定实验一、实验目的1）、通过实验熟悉单回路反馈控制系统的组成和工作原理。
2）、分析分别用P、PI和PID调节时的过程图形曲线。
3）、定性地研究P、PI和PID调节器的参数对系统性能的影响。
4）、通过实验熟悉单回路反馈控制系统的组成和工作原理。
5）、分析分别用P、PI和PID调节时的过程图形曲线。
6）、定性地研究P、PI和PID调节器的参数对系统性能的影响。
二、实验设备THKJ100-1型过程控制实验装置配置：上位机软件、计算机、RS232-485转换器1只、串口线1根、实验连接线。
型参数为串联釜数N三、实验原理实验一原理：阶跃响应测试法是系统在开环运行条件下，待系统稳定后，通过控制器或其他操作器，手动改变对象的输入信号（阶跃信号）。
同时，记录对象的输出数据或阶跃响应曲线，然后根据已给定对象模型的结构形式，对实验数据进行处理，确定模型中各参数。
实验二原理：图13.1单回路上水箱液位控制系统图13.1为单回路上水箱液位控制系统，单回路调节系统一般指在一个调节对象上用一个调节器来保持一个参数的恒定，而调节器只接受一个测量信号，其输出也只控制一个执行机构。
本系统所要保持的恒定参数是液位的给定高度，即控制的任务是控制上水箱液位等于给定值所要求的高度。
根据控制框图，这是一个闭环反馈单回路液位控制，采用工业智能仪表控制。
当调节方案确定之后，接下来就是整定调节器的参数，一个单回路系统设计安装就绪之后，控制质量的好坏与控制器参数选择有着很大的关系。
合适的控制参数，可以带来满意的控制效果。
反之，控制器参数选择得不合适，则会使控制质量变坏，达不到预期效果。
因此，当一个单回路系统组成好以后，如何整定好控制器参数是一个很重要的实际问题。
一个控制系统设计好以后，系统的投运和参数整定是十分重要的工作。
PID控制电机实验报告
PID控制电机实验报告【摘要】本实验通过PID控制电机，对系统进行控制，实现系统的速度调节和位置调节。
首先通过对系统的建模和参数辨识，得到了系统的数学模型和参数，并根据模型设计了合适的PID控制器。
然后通过实验验证了设计的控制器的有效性，实现了对电机速度和位置的调节。
实验结果表明，PID控制器对于系统的速度调节和位置调节具有良好的性能，能够实现较好的控制效果。
【关键词】PID控制；电机；速度调节；位置调节一、实验目的1.通过PID控制器实现对电机的速度调节和位置调节；2.验证PID控制器的有效性和性能。
二、实验原理PID控制器是一种经典的控制策略，由比例（P）、积分（I）和微分（D）三个部分组成。
PID控制器的数学表达式为输出信号u(t) = Kp*e(t) + Ki*∫e(t)dt + Kd*de(t)/dt，其中e(t)为控制偏差，Kp、Ki和Kd分别为比例、积分和微分系数。
在电机控制中，可以将电机看作一个被控对象，输入电机的电压u(t)通过电机的转矩转化为输出角速度ω(t)。
通过对电机的数学建模，可以得到电机的传递函数为G(s)=k/(Ts+1)，其中k为系统增益，T为系统时间常数。
根据系统传递函数的性质，可以得到电机系统的速度和位置闭环模型为Kv(s)=1/(Ts+1)和Kp(s)=Ks/(Ts+1)，分别对应于速度和位置的调节。
三、实验装置1.PC机；2.PID控制器板卡；3.直流电机；4.电压放大电路；5.角度传感器。
四、实验步骤1.建立电机的数学模型，并利用实验数据辨识系统的参数；2.根据模型设计PID控制器的参数；3.连接实验装置，将PC机与PID控制器板卡连接，通过板卡控制电机的电压，实现速度和位置调节；4.设置不同的目标速度和目标位置，进行实验并记录实验数据；5.分析实验数据，评价控制器的性能和有效性。
五、实验结果与分析通过实验得到了电机系统的数学模型为G(s)=2/(s+1)和Ks=10/(s+1)，并根据模型参数设计了PID控制器的参数为Kp=1,Ki=0.01和Kd=0.5、实验中设置了不同的目标速度和目标位置，通过对比实际速度和位置与目标值的差异，评价了控制器的性能。
pid控制实验报告[最新版]
![pid控制实验报告[最新版]](https://img.taocdn.com/s3/m/ab6b0408640e52ea551810a6f524ccbff121ca13.png)
pid控制实验报告pid控制实验报告篇一：PID控制实验报告实验二数字PID控制计算机控制是一种采样控制，它只能根据采样时刻的偏差值计算控制量。
因此连续PID控制算法不能直接使用，需要采用离散化方法。
在计算机PID控制中，使用的是数字PID控制器。
一、位置式PID控制算法按模拟PID控制算法，以一系列的采样时刻点kT代表连续时间t，以矩形法数值积分近似代替积分，以一阶后向差分近似代替微分，可得离散PID位置式表达式：Tu T ?kpeu=para; J=0.0067;B=0.1; dy=zeros= y= -+ = k*ts; %time中存放着各采样时刻rineu_1=uerror_1=error;%误差信号更新图2-1 Simulink仿真程序其程序运行结果如表2所示。
Matlab输出结果errori = error_1 = 表2 例4程序运行结果三、离散系统的数字PID控制仿真1．Ex5 设被控对象为G?num 仿真程序：ex5.m%PID Controller clear all; close all;篇二：自动控制实验报告六-数字PID控制实验六数字PID控制一、实验目的1．研究PID控制器的参数对系统稳定性及过渡过程的影响。
2．研究采样周期T对系统特性的影响。
3．研究I型系统及系统的稳定误差。
二、实验仪器1．EL-AT-III型自动控制系统实验箱一台 2．计算机一台三、实验内容1．系统结构图如6-1图。
图6-1 系统结构图图中 Gc（s）=Kp（1+Ki/s+Kds） Gh（s）=（1－e）/s Gp1（s）=5/（（0.5s+1）（0.1s+1）） Gp2（s）=1/（s（0.1s+1））-TS 2．开环系统（被控制对象）的模拟电路图如图6-2和图6-3，其中图6-2对应GP1（s）,图6-3对应Gp2（s）。
图6-2 开环系统结构图1 图6-3开环系统结构图2 3．被控对象GP1（s）为“0型”系统，采用PI控制或PID控制，可使系统变为“I型”系统，被控对象Gp2（s）为“I型”系统，采用PI控制或PID控制可使系统变成“II型”系统。
PID控制实验报告
实验二 数字PID 控制计算机控制是一种采样控制，它只能根据采样时刻的偏差值计算控制量。
因此连续PID 控制算法不能直接使用，需要采用离散化方法。
在计算机PID 控制中，使用的是数字PID 控制器。
一、位置式PID 控制算法按模拟PID 控制算法，以一系列的采样时刻点kT 代表连续时间t ，以矩形法数值积分近似代替积分，以一阶后向差分近似代替微分，可得离散PID 位置式表达式：∑∑==--++=⎪⎪⎭⎫ ⎝⎛--++=k j di p k j D I p T k e k e k T j e k k e k k e k e T T j e T T k e k k u 00)1()()()())1()(()()()( 式中，D p d I pi T k k T k k ==,，e 为误差信号（即PID 控制器的输入），u 为控制信号（即控制器的输出）。
在仿真过程中，可根据实际情况，对控制器的输出进行限幅。
二、连续系统的数字PID 控制仿真连续系统的数字PID 控制可实现D/A 及A/D 的功能，符合数字实时控制的真实情况，计算机及DSP 的实时PID 控制都属于这种情况。
1．Ex3 设被控对象为一个电机模型传递函数BsJs s G +=21)(，式中J=0.0067,B=0.1。
输入信号为)2sin(5.0t π，采用PD 控制，其中5.0,20==d p k k 。
采用ODE45方法求解连续被控对象方程。
因为Bs Js s U s Y s G +==21)()()(，所以u dt dy B dty d J =+22，另y y y y ==2,1，则⎪⎩⎪⎨⎧+-==/J)*u ((B/J)y y y y 12221 ，因此连续对象微分方程函数ex3f.m 如下 function dy = ex3f(t,y,flag,para)u=para;J=0.0067;B=0.1;dy=zeros(2,1);dy(1) = y(2);dy(2) = -(B/J)*y(2) + (1/J)*u;控制主程序ex3.mclear all;close all;ts=0.001; %采样周期xk=zeros(2,1);%被控对象经A/D转换器的输出信号y的初值e_1=0;%误差e(k-1)初值u_1=0;%控制信号u(k-1)初值for k=1:1:2000 %k为采样步数time(k) = k*ts; %time中存放着各采样时刻rin(k)=0.50*sin(1*2*pi*k*ts); %计算输入信号的采样值para=u_1; % D/AtSpan=[0 ts];[tt,xx]=ode45('ex3f',tSpan,xk,[],para); %ode45解系统微分方程%xx有两列，第一列为tt时刻对应的y，第二列为tt时刻对应的y导数xk = xx(end,:); % A/D，提取xx中最后一行的值，即当前y和y导数yout(k)=xk(1); %xk(1)即为当前系统输出采样值y(k)e(k)=rin(k)-yout(k);%计算当前误差de(k)=(e(k)-e_1)/ts; %计算u(k)中微分项输出u(k)=20.0*e(k)+0.50*de(k);%计算当前u(k)的输出%控制信号限幅if u(k)>10.0u(k)=10.0;endif u(k)<-10.0u(k)=-10.0;end%更新u(k-1)和e(k-1)u_1=u(k);e_1=e(k);endfigure(1);plot(time,rin,'r',time,yout,'b');%输入输出信号图xlabel('time(s)'),ylabel('rin,yout');figure(2);plot(time,rin-yout,'r');xlabel('time(s)'),ylabel('error');%误差图程序运行结果显示表1所示。
数字PID控制实验
4.5.1数字PID 控制实验 1 标准PID 控制算法1. 一. 实验要求2. 了解和掌握连续控制系统的PID 控制的原理。
3. 了解和掌握被控对象数学模型的建立。
4. 了解和掌握数字PID 调节器控制参数的工程整定方法。
观察和分析在标准PID 控制系统中, P.I.D 参数对系统性能的影响。
二. 实验内容及步骤 ⑴ 确立模型结构本实验采用二个惯性环节串接组成实验被控对象, T1=0.2S, T2=0.5S Ko=2。
S e T K s G τ-+⨯≈+⨯+=1S 110.2S 21S 5.01)(000⑵ 被控对象参数的确认被控对象参数的确认构成如图4-5-10所示。
本实验将函数发生器（B5）单元作为信号发生器, 矩形波输出（OUT ）施加于被测系统的输入端R, 观察矩形波从0V 阶跃到+2.5V 时被控对象的响应曲线。
图4-5-10 被控对象参数的确认构成实验步骤: 注: 将‘S ST ’用‘短路套’短接!① 在显示与功能选择（D1）单元中, 通过波形选择按键选中‘矩形波’（矩形波指示灯亮）。
② B5的量程选择开关S2置下档, 调节“设定电位器1”, 使之矩形波宽度>2秒（D1单元左显示）。
③ 调节B5单元的“矩形波调幅”电位器使矩形波输出电压= 2.5V 左右（D1单元右显示）。
④ 构造模拟电路: 按图4-5-10安置短路套及测孔联线, 表如下。
（a ）安置短路套 （b ）测孔联线⑤ 运行、观察、记录:A)先运行LABACT 程序, 选择界面的“工具”菜单选中“双迹示波器”（Alt+W ）项, 弹出双迹示波器的界面, 点击开始, 用虚拟示波器观察系统输入信号。
图4-5-11 被控对象响应曲线B) 在图4-5-112被控对象响应曲线上测得t1和t2。
通常取 , 要求从图中测得 ； 通常取 , 要求从图中测得 。
计算 和 : 0.84730.3567t -1.204t )]t (y 1[ln -)]t (y 1[ln )]t (y 1[ln t )]t (y 1[n t 0.8473t t )]t (y 1[ln -)]t (y 1[ln t t T 212010201102122010120==-----=-=---=τC) 求得数字PID 调节器控制参数P K 、I T 、D T （工程整定法）)/0.2(1)/0.37()/0.6(1)/0.5()/2.5(]27.0)/(35.1[10000200000T T T T T T T T T T K K D I P ττττττ+⨯=++⨯=+=据上式计算数字PID 调节器控制参数P K 、I T 、D T⑶ 数字PID 闭环控制系统实验数字PID 闭环控制系统实验构成见图4-5-12, 观察和分析在标准PID 控制系统中, P.I.D 参数对系统性能的影响, 分别改变P.I.D 参数, 观察输出特性, 填入实验报告,模块号 跨接座号 1 A5 S5, S7, S102 A7 S2, S7, S9, P3 B5‘S-ST ’1 输入信号R B5（OUT ）→A5（H1）2 运放级联 A5A （OUTA ）→A7（H1）3 示波器联接 ×1档B5（OUT ）→B3（CH1） 4A7A （OUTA ）→B3（CH2）图4-5-12 数字PID 闭环控制系统实验构成实验步骤: 注: 将‘S ST ’用‘短路套’短接!① 在显示与功能选择（D1）单元中, 通过波形选择按键选中‘矩形波’（矩形波指示灯亮）。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
实验二 数字PID 控制计算机控制是一种采样控制，它只能根据采样时刻的偏差值计算控制量。
因此连续PID 控制算法不能直接使用，需要采用离散化方法。
在计算机PID 控制中，使用的是数字PID 控制器。
一、位置式PID 控制算法按模拟PID 控制算法，以一系列的采样时刻点kT 代表连续时间t ，以矩形法数值积分近似代替积分，以一阶后向差分近似代替微分，可得离散PID 位置式表达式： 式中，D p d I pi T k k T k k ==,，e 为误差信号（即PID 控制器的输入），u 为控制信号（即控制器的输出）。
在仿真过程中，可根据实际情况，对控制器的输出进行限幅。
二、连续系统的数字PID 控制仿真连续系统的数字PID 控制可实现D/A 及A/D 的功能，符合数字实时控制的真实情况，计算机及DSP 的实时PID 控制都属于这种情况。
1．Ex3 设被控对象为一个电机模型传递函数BsJs s G +=21)(，式中J=,B=。
输入信号为)2sin(5.0t π，采用PD 控制，其中5.0,20==d p k k 。
采用ODE45方法求解连续被控对象方程。
因为Bs Js s U s Y s G +==21)()()(，所以u dt dy B dty d J =+22，另y y y y ==2,1，则⎪⎩⎪⎨⎧+-==/J)*u ((B/J)y y y y 12221 ，因此连续对象微分方程函数如下 function dy = ex3f(t,y,flag,para)u=para;J=;B=;dy=zeros(2,1);dy(1) = y(2);dy(2) = -(B/J)*y(2) + (1/J)*u;控制主程序clear all;close all;ts=; %采样周期xk=zeros(2,1);%被控对象经A/D 转换器的输出信号y 的初值e_1=0;%误差e(k-1)初值u_1=0;%控制信号u(k-1)初值for k=1:1:2000 %k 为采样步数time(k) = k*ts; %time中存放着各采样时刻rin(k)=*sin(1*2*pi*k*ts); %计算输入信号的采样值para=u_1; % D/AtSpan=[0 ts];[tt,xx]=ode45('ex3f',tSpan,xk,[],para); %ode45解系统微分方程%xx有两列，第一列为tt时刻对应的y，第二列为tt时刻对应的y导数xk = xx(end,:); % A/D，提取xx中最后一行的值，即当前y和y导数yout(k)=xk(1); %xk(1)即为当前系统输出采样值y(k)e(k)=rin(k)-yout(k);%计算当前误差de(k)=(e(k)-e_1)/ts; %计算u(k)中微分项输出u(k)=*e(k)+*de(k);%计算当前u(k)的输出%控制信号限幅if u(k)>u(k)=;endif u(k)<u(k)=;end%更新u(k-1)和e(k-1)u_1=u(k);e_1=e(k);endfigure(1);plot(time,rin,'r',time,yout,'b');%输入输出信号图xlabel('time(s)'),ylabel('rin,yout');figure(2);plot(time,rin-yout,'r');xlabel('time(s)'),ylabel('error');%误差图程序运行结果显示表1所示。
表1 程序运行结果分析:输出跟随输入，PD 控制中，微分控制可以改善动态特性，调节时间缩短，允许加大比例控制，使稳态误差减小，提高了控制精度.2．Ex4 被控对象是一个三阶传递函数ss s 1047035.8752350023++，采用Simulink 与m 文件相结合的形式，利用ODE45方法求解连续对象方程，主程序由Simulink 模块实现，控制器由m 文件实现。
输入信号为一个采样周期1ms 的正弦信号。
采用PID 方法设计控制器，其中05.0,2,5.1===d i p k k k 。
误差初始化由时钟功能实现，从而在m 文件中实现了误差的积分和微分。
控制主程序：控制子程序：function [u]=ex4f(u1,u2)%u1为Clock ，u2为图2-1中Sum 模块输出的误差信号e 的采样值persistent errori error_1if u1==0 %当Clock=0时，即初始时，e(k)=e(k-1)=0errori=0error_1=0endts=;kp=;ki=;kd=;error=u2;errord=(error-error_1)/ts;%一阶后向差分误差信号表示的误差微分errori=errori+error*ts;%累积矩形求和计算的误差的积分u=kp*error+kd*errord+ki*errori;%由PID算式得出的当前控制信号u(k) error_1=error;%误差信号更新图2-1 Simulink仿真程序其程序运行结果如表2所示。
Matlab输出结果errori =error_1 =表2 例4程序运行结果三、离散系统的数字PID 控制仿真1．Ex5 设被控对象为s s s s G 1047035.87523500)(23++=，采样时间为1ms ，对其进行离散化。
针对离散系统的阶跃信号、正弦信号和方波信号的位置响应，设计离散PID 控制器。
其中S 为信号选择变量，S=1时是阶跃跟踪，S=2时为方波跟踪，S=3时为正弦跟踪。
求出G(s)对应的离散形式)()()(z U z Y z G =，其中Y(z)和U(z)是关于z 的多项式，则可以得到其对应的差分表达式仿真程序：%PID Controllerclear all;close all;ts=;%采样周期sys=tf,[1,,,0]);%被控对象连续传递函数dsys=c2d(sys,ts,'z');%转换成离散z传递函数的形式[num,den]=tfdata(dsys,'v');%提取z传递函数中的分子和分母多项式系数u_1=;u_2=;u_3=;%u(k-1)、u(k-2)、u(k-3)的初值y_1=;y_2=;y_3=; %y(k-1)、y(k-2)、y(k-3)的初值x=[0,0,0]'; %比例、微分、积分项的初值error_1=0;%e(k-1)的初值disp('S=1--step,S=2--sin,S=3--square')% S=1阶跃，S=2方波，S=3正弦S=input('Number of input signal S:')%接收输入信号代号for k=1:1:1500time(k)=k*ts;%各采样时刻if S==1 %阶跃输入时kp=;ki=;kd=; %各项PID系数rin(k)=1; %阶跃信号输入elseif S==2kp=;ki=;kd=; %各项PID系数rin(k)=sign(sin(2*2*pi*k*ts)); %方波信号输入elseif S==3kp=;ki=;kd=; %各项PID系数rin(k)=*sin(2*2*pi*k*ts); %正弦信号输入endu(k)=kp*x(1)+kd*x(2)+ki*x(3); %PID控制信号输出u(k)%控制信号输出限幅if u(k)>=10u(k)=10;endif u(k)<=-10u(k)=-10;end%根据差分方程计算系统当前输出y(k)yout(k)=-den(2)*y_1-den(3)*y_2-den(4)*y_3+num(2)*u_1+num(3)*u_2+num( 4)*u_3;error(k)=rin(k)-yout(k);%当前误差%更新u(k-1)、u(k-2)、u(k-3)、y(k-1)、y(k-2)、y(k-3) u_3=u_2;u_2=u_1;u_1=u(k);y_3=y_2;y_2=y_1;y_1=yout(k);x(1)=error(k); %比例输出x(2)=(error(k)-error_1)/ts; %微分输出x(3)=x(3)+error(k)*ts; %积分输出error_1=error(k); %更新e(k-1)endfigure(1); %作图plot(time,rin,'r',time,yout,'b');xlabel('time(s)'),ylabel('rin,yout');其程序运行结果如表3所示。
2．Ex6针对于Ex5被控对象所对应的离散系统，设计针对三角波、锯齿波和随机信号的位置式响应。
仿真程序：。
程序中当S=1时为三角波，S=2时为锯齿波，S=3时为随机信号。
如果D=1，则通过pause命令实现动态演示仿真。
%PID Controllerclear all;close all;ts=;sys=tf,[1,,,0]);dsys=c2d(sys,ts,'z');[num,den]=tfdata(dsys,'v');u_1=;u_2=;u_3=;r_1=rand;y_1=0;y_2=0;y_3=0;x=[0,0,0]';error_1=0;disp('S=1--Triangle,S=2--Sawtooth,S=3--Random')% S=1三角，S=2锯齿，S=3随机S=input('Number of input signal S:')%接收输入信号代号disp('D=1--Dynamic display,D~=1--Direct display')%D=1动画显示，D~=1直接显示D=input('D=')for k=1:1:3000time(k)=k*ts;kp=;ki=;kd=;if S==1 %Triangle Signalif mod(time(k),2)<1rin(k)=mod(time(k),1);elserin(k)=1-mod(time(k),1);endrin(k)=rin(k);endif S==2 %Sawtooth Signalrin(k)=mod(time(k),;endif S==3 %Random Signalrin(k)=rand;vr(k)=(rin(k)-r_1)/ts; %Max speed iswhile abs(vr(k))>=rin(k)=rand;vr(k)=abs((rin(k)-r_1)/ts);endendu(k)=kp*x(1)+kd*x(2)+ki*x(3); %PID Controller%Restricting the output of controllerif u(k)>=10u(k)=10;endif u(k)<=-10u(k)=-10;end%Linear modelyout(k)=-den(2)*y_1-den(3)*y_2-den(4)*y_3+num(2)*u_1+num(3)*u_2+num( 4)*u_3;error(k)=rin(k)-yout(k);r_1=rin(k);u_3=u_2;u_2=u_1;u_1=u(k);y_3=y_2;y_2=y_1;y_1=yout(k);x(1)=error(k); %Calculating P x(2)=(error(k)-error_1)/ts; %Calculating D x(3)=x(3)+error(k)*ts; %Calculating I xi(k)=x(3);error_1=error(k);if D==1 %Dynamic Simulation Displayplot(time,rin,'b',time,yout,'r');pause;endendplot(time,rin,'r',time,yout,'b');xlabel('time(s)');ylabel('rin,yout');Matlab运行结果为：S=1--Triangle,S=2--Sawtooth,S=3--RandomNumber of input signal S:1（2、3）S =1D=1--Dynamic display,D=0--Direct display %D=1动画显示，D~=1直接显示D=0D =0 %D=0直接显示，如果D=1，则通过pause命令实现动态演示仿真。