实验四、最少拍控制算法matlab仿真实验
创始人
2024-06-02 16:57:12
0

实验、最少拍控制算法matlab仿真实验

一、实验目的

1.掌握最少拍有纹波、无纹波系统的设计方法;

2.学会最少拍控制系统的分析方法;

3.了解输入信号对最少拍控制系统的影响及改进措施

二、实验工具

MATLAB 软件( 2022a)  

、实验要求

1、写出广义被控对象的脉冲传递函数G(z)。并求出广义被控对象的差分方程形式。

2、画出未加控制器时系统的单位阶跃响应图,分析系统是否稳定?若稳定,求稳态误差。

3、根据最少拍设计步骤,写出针对单位阶跃输入时的闭环脉冲传递函数和误差脉冲脉冲传递函数1-。给出设计出的最少拍控制器的脉冲传递函数D(z)。并求出最少拍控制器的差分方程形式。

4、画出串上最少拍控制器后的数字控制系统框图。

5、将上述所求的最少拍控制器代入系统,画出系统在单位阶跃输入信号作用下的响应图,以及数字控制器的控制量变化曲线。分析加入控制器后系统能否达到最少拍控制的要求。给出系统输出序列和控制器输出序列,并求出调节时间,讨论输出是否存在纹波。

6、若存在纹波,请设计最少拍无纹波控制器,并给出串上最少拍无纹波控制器后的闭环脉冲传递函数和误差脉冲脉冲传递函数1-。最少拍无纹波控制器的脉冲传递函数D(z)。

7、将得到的最少拍无纹波控制器代入系统,画出系统在单位阶跃输入信号作用下的响应图,以及数字控制器的控制量变化曲线。给出系统输出序列和控制器输出序列,分析加入最少拍无纹波控制器后系统能否达到消除纹波的目的。为消除纹波,调节时间变为多少。

8、写出针对单位斜坡信号输入时的闭环脉冲传递函数和误差脉冲脉冲传递函数1-。给出设计出的最少拍控制器的脉冲传递函数D(z)

、实验内容

设单位反馈线性定常系统的连续部分和零阶保持器的传递函数分别为:

 

,采样周期T=1秒。要求系统在单位阶跃输入时实现最少拍控制,求数字控制器的脉冲传递函数。讨论加上控制器后的系统输出在采样点之间是否存在纹波,若存在纹波,请设计最少拍无纹波控制器。系统结构如图:

 

 

  1. 对未加控制器时的闭环系统分析其是否满足最少拍要求

 

  

由开环传递函数可知,开环脉冲传递函数有z=1的极点个数为1,系统为I型系统。

对于单位速度输入,稳态误差不为零。

即原来的闭环系统达不到最少拍系统中对某个输入稳定无静差的要求。

(3)通过matlab对函数进行绘制并得到具体的输出序列进行对比得到准确的结果。

 

原闭环系统输出序列为:

 

可知,系统输出在第十三拍以后,基本稳定为1,满足稳定性要求。调节时间为13s。动态过程为13拍,达不到最少怕的要求。

原闭环系统偏差序列为:

 

可知,系统偏差输出在第十三拍(系统稳定)以后,基本保持0大小,无静差。

  所以可以从程序输出的内容分析可看出,原系统不满足最少拍要求。需要串入最少拍调节器。

 

 

 

 

 

可知,加数字控制器后,系统输出与系统输入在1s以后相等,说明系统是稳定的且稳态无静差。调节时间为1s,但观察其数字控制器输出曲线及输出序列发现,其控制器输出在28拍(28s)后才变为恒定值,说明输出在第3拍至第28拍间,仅能保证在采样点上与输入相同,在采样点间不相同,即输出有纹波。

 

 

 

 

 

可以看出,此时调节时间变为2s,调节时间比原来延长1拍,但其控制器输出也在2拍后变为恒定值,可以保证2拍后,输出不仅在采样点上与输入相同,还在采样点间保证相同,实现输出无纹波。

实验思考

1. 最少拍受什么限制而使调整节拍增加?

答:调整节拍增加最少拍受输入函数R(z)的阶数限制。R(z)的阶数越高,调整时间越长。

2.分析不同输入信号对最少拍控制系统的影响。

答:与输入信号有关。输入信号的阶数越高,其调整时间越长。

3.无纹波系统对控制器有何要求。

答:有两个要求。

第一,被控对象G(s)要有足够的的积分环节,第二,D(z) 必须包含G(z)中的圆内圆外的全部零点N(z);

六、代码汇总

对未加控制器时的闭环系统分析其是否满足最少拍要求

  1. clear;  
  2. clc;  
  3. G=zpk([],[0,-1],1);           %零极点形式的被控对象  
  4.    
  5. T=1;                         %采样周期  
  6. Gd=c2d(G,T,'zoh')           %广义被控对象脉冲传函  
  7. [num1,den1]=tfdata(Gd,'v');  
  8. GG=feedback(Gd,1)           %原系统闭环传递函数  
  9.    
  10.  yd_1=0;yd_2=0;  
  11.  ed_1=0;ed_2=0;  
  12.    
  13. for k=1:1:35  
  14.     time(k)=k*T;  
  15.     rin(k)=1;      %单位阶跃信号  
  16.    %闭环系统的差分方程  
  17. yd(k)=num1(2)*ed_1+num1(3)*ed_2-den1(2)*yd_1-den1(3)*yd_2;  
  18.  ed(k)=rin(k)-yd(k);      %求偏差  
  19.    
  20. %数据更新  
  21.  ed_2=ed_1;ed_1=ed(k);  
  22. yd_2=yd_1;yd_1=yd(k);  
  23.    
  24. end  
  25. yd   %原闭环系统的输出序列  
  26. ed   %原闭环系统的偏差序列  
  27. figure(1);  
  28. plot(time,rin,'r',time,yd,'g',time,ed,'b');  
  29. xlabel('time');  
  30. ylabel('rin,yd');  
  31. title('闭环系统在单位速度信号下的响应');  
  32. text(3,23,'红线:输入信号;绿线:系统输出响应;蓝线:偏差');  
  33.   
  34. %-------------------------------设计最少拍控制器  -------------------------------
  35. clear all;  
  36. clc;  
  37. G=zpk([],[0,-1],1);           %零极点形式的被控对象  
  38.    
  39. T=1;                         %采样周期  
  40. Gd=c2d(G,T,'zoh');           %广义被控对象脉冲传函  
  41.    
  42. phi=tf([1],[1 0],T);     %闭环脉冲传函  
  43. phie=tf([1 -1],[1 0],T);    %误差脉冲传函  
  44. D=phi/(Gd*phie);               %数字控制器脉冲传函D(z)  
  45. phie,phi,Gd,D  
  46. [num1,den1]=tfdata(Gd,'v');  
  47. [num2,den2]=tfdata(D,'v');  
  48.    
  49.  ud_1=0;ud_2=0;ud_3=0;  
  50.  yd_1=0;yd_2=0;  
  51.  ed_1=0;ed_2=0;ed_3=0;  
  52.    
  53. for k=1:1:35  
  54.     time(k)=k*T;  
  55.     rin(k)=1;      %单位阶跃信号  
  56.    %加数字控制器的差分方程  
  57. yd(k)=num1(2)*ud_1+num1(3)*ud_2-den1(2)*yd_1-den1(3)*yd_2;  
  58.     ed(k)=rin(k)-yd(k);      %求偏差  
  59.    
  60. %数字控制器的差分方程  
  61. ud(k)=(num2(1)*ed(k)+num2(2)*ed_1+num2(3)*ed_2+num2(4)*ed_3-den2(2)*ud_1-den2(3)*ud_2-den2(4)*ud_3)/den2(1);  
  62.    
  63. %数据更新  
  64.     ud_3=ud_2;ud_2=ud_1;ud_1=ud(k);  
  65.     ed_3=ed_2;ed_2=ed_1;ed_1=ed(k);  
  66. yd_2=yd_1;yd_1=yd(k);  
  67.    
  68. end  
  69. yd,ud  
  70. figure(2);  
  71. plot(time,rin,'r',time,yd,'b',time,ud,'k');  
  72. xlabel('time');  
  73. ylabel('rin,yd,ud');  
  74. title('加数字控制器后单位阶跃信号的闭环响应及数字控制器的输出');  
  75. text(3,1.5,'红线:系统输入;蓝线:系统输出;黑线:控制器输出');  
  76. %-------------------------无纹波控制器设计-------------------------------  
  77. clear all;  
  78. close all;  
  79. clc;  
  80. G=zpk([],[0,-1],1);           %%零极点形式的被控对象  
  81. T=1;                         %采样周期  
  82. Gd=c2d(G,T,'zoh');           %广义被控对象脉冲传函  
  83.    
  84. %求误差脉冲传函及闭环脉冲传函  
  85. phi=tf([1/1.7183 0.7183/1.7183],[1 0 0],T);     %闭环脉冲传函  
  86. phie=tf([1 -1/1.7183 -0.7183/1.7183],[1 0 0],T);    %误差脉冲传函  
  87. D=phi/(Gd*phie);               %数字控制器脉冲传函D(z)  
  88. Gd,D,phie,phi  
  89.    
  90. [num1,den1]=tfdata(Gd,'v');%求脉冲传函的分子分母多项式  
  91. [num2,den2]=tfdata(D,'v');  
  92.    
  93.  ud_1=0;ud_2=0;ud_3=0;ud_4=0;ud_5=0;  
  94.  yd_1=0;yd_2=0;  
  95.  ed_1=0;ed_2=0;ed_3=0;ed_4=0;ed_5=0;  
  96.    
  97. for k=1:1:25  
  98.     time(k)=k*T;  
  99.     rin(k)=1;      %单位阶跃信号  
  100.    
  101.      %加数字控制器的差分方程  
  102. yd(k)=num1(2)*ud_1+num1(3)*ud_2-den1(2)*yd_1-den1(3)*yd_2;  
  103.     ed(k)=rin(k)-yd(k);  
  104.    
  105. %数字控制器的差分方程  
  106. ud(k)=(num2(1)*ed(k)+num2(2)*ed_1+num2(3)*ed_2+num2(4)*ed_3+num2(5)*ed_4+num2(6)*ed_5-den2(2)*ud_1-den2(3)*ud_2-den2(4)*ud_3-den2(5)*ud_4-den2(6)*ud_5)/den2(1);  
  107.    
  108.     ud_5=ud_4;ud_4=ud_3;ud_3=ud_2;ud_2=ud_1;ud_1=ud(k);  
  109.     ed_5=ed_4;ed_4=ed_3;ed_3=ed_2;ed_2=ed_1;ed_1=ed(k);  
  110. yd_2=yd_1;yd_1=yd(k);  
  111.    
  112. end  
  113. yd,ud  
  114.    
  115. figure(1);  
  116. plot(time,rin,'r',time,yd,'b',time,ud,'k');  
  117. xlabel('time');  
  118. ylabel('rin,yd,ud');  
  119. title('加数字控制器后单位阶跃信号的闭环响应及数字控制器的输出');  
  120. text(1.5,1.5,'红线:系统输入;蓝线:系统输出;黑线:控制器输出');  

相关内容

热门资讯

安卓系统账号怎么改,轻松实现账... 你是不是也和我一样,在使用安卓手机的时候,突然觉得账号名字太老土了,想要来个焕然一新的改变呢?别急,...
离线语音系统安卓版下载,随时随... 你有没有想过,在手机上也能实现语音助手的功能,而且完全不需要联网?没错,就是那种离线语音系统,听起来...
华为os系统是基于安卓系统吗,... 你有没有想过,华为的手机里那个神秘的OS系统,它是不是就是安卓系统呢?别急,今天就来揭开这个谜底,让...
各品牌安卓系统多大,探索各大品... 你有没有想过,那些我们每天离不开的安卓手机,它们背后的系统到底有多大呢?这可不是一个小问题哦,因为系...
不是安卓系统电视机,探索非安卓... 你有没有想过,家里的电视是不是安卓系统呢?现在市面上,安卓系统电视机可是越来越流行了。但是,你知道吗...
谷歌安卓系统开源免费用,免费体... 你知道吗?在科技的世界里,有时候最让人惊喜的就是那些免费又好用的东西。今天,就让我来给你揭秘一个超级...
安卓电脑版怎么装系统,轻松实现... 你有没有想过,你的安卓电脑版突然间卡得像蜗牛一样,慢得让人抓狂?别急,今天就来教你怎么给它来个焕然一...
安卓系统有几种语音,揭秘多样化... 你知道吗?安卓系统里的语音功能可真是让人爱不释手呢!想象你只需要动动嘴,就能完成各种操作,是不是觉得...
扫码系统推荐安卓游戏,畅享指尖... 你有没有发现,最近手机里的游戏推荐越来越贴心了?没错,就是那个神奇的扫码系统!今天,就让我带你一探究...
手机系统安卓做的最好,探索最佳... 你有没有发现,现在手机市场上,安卓系统简直是要风得风,要雨得雨,简直是手机系统中的佼佼者!不信?那就...
安卓系统越拉越难用,越拉越难用... 你有没有发现,最近安卓系统越来越难用了?是不是每次打开手机,都感觉卡得要命,各种广告弹窗让人头疼不已...
卓安系统创始人,引领科技浪潮的... 你有没有想过,在这个信息爆炸的时代,谁才是那个站在风口浪尖,引领科技潮流的先锋?今天,我要带你认识一...
安卓10系统经常断网,安卓10... 最近是不是你也遇到了安卓10系统经常断网的问题?这可真是让人头疼啊!手机断网,就像失去了方向感的指南...
开发安卓系统用什么,从工具到实... 你有没有想过,为什么你的手机里装了那么多应用,却还是觉得不够用呢?其实,这一切都离不开一个强大的操作...
安卓系统如何双开抖音,安卓系统... 如何在安卓系统上双开抖音:让你的娱乐生活更加丰富多彩在数字化时代,手机已经成为我们生活中不可或缺的一...
安卓8.0系统内测地址,探索系... 你有没有听说最近安卓8.0系统内测的消息?这可是科技圈的大事啊!作为一个紧跟潮流的数码爱好者,我可是...
安卓系统可以账号同步吗,安卓系... 你有没有想过,你的安卓手机里的各种账户信息,比如微信、QQ、支付宝,能不能像变魔术一样,一键同步到另...
安卓系统的平板怎么下载,畅享应... 你有没有想过,拥有一款安卓系统的平板电脑,那感觉简直就像拥有了整个世界?想象随时随地都能畅游网络海洋...
win平板上装安卓系统咋样,开... 亲爱的读者们,你是否曾想过在Windows平板上体验安卓系统的魅力?想象那双熟悉的Windows界面...
给安卓车机写系统,安卓车机系统... 你有没有想过,你的安卓车机系统其实可以焕然一新?想象当你驾驶着爱车,车机界面流畅得像丝滑的巧克力,功...