《计算机控制—史密斯预估器编程》由会员分享，可在线阅读，更多相关《计算机控制—史密斯预估器编程（13页珍藏版）》请在金锄头文库上搜索。
1、1东南大学能源与环境学院实 验 报 告课程名称： 计算机控制概论 实验名称： Smith 预估控制实验 院 （系）： 能源与环境学院 专 业： 热能与动力专业 姓 名： 李栋 学 号： 03012117 实 验 室： 实验组别： 同组人员： 实验时间： 2015 年 04 月 10 日评定成绩： 审阅教师： 2目 录一、实验目的3二、实验内容3三、实验步骤3四、实验分析123Smith 预估控制实验一 实验目的通过实验掌握 Smith 预估控制的方法及程序编制及调试。二 实验内容 Smith 预估控制系统如图所示， 图一对象 G(S)= Ke-s / (1+T1S) ,K = 1, T1 =
2、10 s , = 5 s ,Wc(z)采用数字 PI 控制规律。对象扰动实验画出 U(t) = u01(t)时，y(t)曲线。Smith 预估控制(1) 构造 W (S)，求出 W (Z)。(2) 整定 Wc(s)(按什么整定？)(3) 按图仿真，并打印曲线。（4） 改变 W (S)中 K，（对象不变） ，进行仿真比较，观察它们对调节过程的影响。三 实验步骤1、对象扰动实验源程序如下：#includeiostream.h#includemath.h#includefstream.h4void main()fstream outfile(data1.xls,ios:out);double t;d
3、ouble u0;coutt;coutu0;double ee=pow(2.718,(-t/10.0);int N;int i;double u100,y100;for(i=0;ik;7coutt1;coutt;coutkp;coutki;double ee=pow(2.718,(-t/10.0);int N,N1;int i;double r100,e1100,e2100,cm100,q100,u100,y100;for(i=0;i0&i=N)yi=(1-ee)*ui-N+yi-1*ee;if(i=N1)e1i=ri-yi-1;cmi=ee*cmi-1+k*(1-ee)*ui-1;qi=cm
4、i-cmi-N1;e2i=e1i-qi;ui=ui-1+kp*(e2i-e2i-1)+ki*e2i;if(i=N)yi=(1-ee)*ui-N+yi-1*ee;9for(i=0;i*t100;i+)coutyit;for(i=0;i*t100;i+)outfilei*tt;outfilen;for(i=0;i*t100;i+)outfileyit;outfile.close();（3）输出结果：以下所涉及到的采样周期均为 T=1，PI 控制器的参数均为 Kp=1，Ki=1；当 Smith 预估器中的 K=1，延迟时间 =5 时（即与对象的特性完全符合）：Y（t）输出数据：0 0 0 0 0 0
5、 0.190306 0.421441 0.663641 0.8917551.08676 1.23639 1.37128 1.47104 1.5311 1.549551.52761 1.46956 1.38931 1.29344 1.18983 1.085670.987246 0.89981 0.828799 0.776983 0.745653 0.734524100.741955 0.765251 0.801257 0.846217 0.896223 0.947450.996402 1.04011 1.07631 1.1035 1.1209 1.12848 1.126831.11708 1.1
6、0079 1.07973 1.05581 1.03093 1.00680.984919 0.966463 0.952253 0.942744 0.938032 0.937890.941816 0.949101 0.958895 0.970279 0.982333 0.9941951.00511 1.01448 1.02186 1.02698 1.02978 1.030321.02882 1.02561 1.02108 1.01569 1.00987 1.004060.998627 0.993893 0.990086 0.98735 0.985745 0.9852490.985771 0.987
7、163 0.989238 0.991783 0.994581 0.997421.00011 1.0025 1.00445 1.0059 1.0068 1.00715 1.0071.00641 1.00547 1.00428 1.00293 1.00155 1.000220.999027 0.998028 0.997269 0.996773扰动曲线如下：图三当 Smith 预估器中的 K=1，延迟时间 =2 时（即与对象的特性不完全符合）：Y（t）输出数据如下：0 0 0 0 0 0 0.190306 0.421441 0.663641 0.927971111.21095 1.50619 1.8
8、1053 2.08577 2.31463 2.489892.60123 2.63889 2.59562 2.46564 2.25095 1.958931.59989 1.18774 0.740093 0.277571 -0.176632 -0.598368 -0.963966 -1.25121 -1.44044 -1.51579 -1.4662 -1.28642 -0.977633 -0.547714 -0.01125320.610765 1.29164 1.999962.70093 3.358 3.93455 4.39588 4.71103 4.854644.80862 4.56351 4.
9、11952 3.48712 2.68715 1.750360.716479 -0.367272 -1.44817 -2.47036 -3.37751 -4.11571 -4.63639 -4.89916 -4.87439 -4.54543 -3.91026 -2.98249 -1.79168 -0.38278 1.18524 2.8415 4.5062 6.09408 7.518558.69603 9.55045 10.0176 10.0494 9.61689 8.713477.35632 5.58704 3.47109 1.09587 -1.43244 -3.99312 -6.45626 -
10、8.68888 -10.5616 -11.9554 -12.7687 -12.9234 -12.3704 -11.0941 -9.11507 -6.49149 -3.31832 0.2752394.13026 8.06445 11.88 15.3731 18.3435扰动曲线如下：图四当 Smith 预估器中的 K=2，延迟时间 =2 时（即与对象的特性不完全符合）：12Y（t）输出数据如下：0 0 0 0 0 0 0.190306 0.385225 0.546344 0.7250840.920371 1.11455 1.30834 1.46909 1.59338 1.692661.7608
11、1.79027 1.78227 1.73766 1.66147 1.560211.43778 1.29949 1.15302 1.00558 0.863901 0.7341210.621319 0.529913 0.463425 0.423874 0.411896 0.4269230.467201 0.529943 0.611457 0.707298 0.812552 0.9221031.03084 1.13389 1.22683 1.30585 1.36793 1.410941.4337 1.43598 1.41848 1.38278 1.33121 1.266721.19274 1.112
12、98 1.03127 0.951381 0.876845 0.8108160.75594 0.714253 0.687116 0.675179 0.67838 0.6959770.726605 0.768367 0.818936 0.875681 0.935797 0.9964341.05484 1.10845 1.15505 1.19281 1.22037 1.236891.24206 1.23609 1.21971 1.19405 1.16064 1.12131.07804 1.03296 0.988182 0.945705 0.907359 0.8747110.849012 0.8311
13、46 0.82161 0.820506 0.82755 0.8421020.863208 0.889656 0.920041 0.952835 0.986462 1.01937扰动曲线如下：图五13四 实验分析当系统是特征方程中含有纯迟延项的时候，系统的闭环稳定性事下降的，当迟延时间 比较大的时候，系统就会不稳定。因此采用常规的控制是难以使系统获得满意的控制性能的。理论上，一个被控对象的过程可分为纯迟延环节和 Gp（s） （不含有纯迟延项） ，如果虚拟变量 C 可用某种方法测量，并作为反馈量连接到控制器，就可以把纯迟项移到闭环的外面。因为在反馈信号中没有迟延，系统的响应将大大得到改善，同时在外
14、回路用第二个反馈构成 Smith 预估器控制系统，当中的D（s）控制器采用常规的 PI 或 PID 控制器。采用 Smith 预估器的闭环传递函数为Y（s）/R（s）=D(S)G(S)/(1+D(s)G（s）)* e -s在迟延项从闭环特征方程中去掉后，稳定性将得到改善，控制器可整定的更好。实际中，Smith 预估器主要用作对过程纯迟延的补偿。但是，在实际中，对象的特性是很难精确标示出来的，这样就决定了 Smith 预估器的参数是同对象本身有差别的。若差别很小或没有差别，则对控制效果比较好，没有产生多大的影响如图三所示。但是如果差别很大，预估器反而会使得系统稳定性大大降低，甚至很容易出现不稳定，如图四、五所示，就是这种情况。
