实验四、最少拍控制算法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,'红线:系统输入;蓝线:系统输出;黑线:控制器输出');  

相关内容

热门资讯

哪个安卓机系统好用,探索安卓系... 你有没有想过,手机里的安卓系统就像是个大厨,不同的系统就像不同的烹饪手法,有的让你吃得津津有味,有的...
安卓如何控制苹果系统,从安卓到... 你知道吗?在这个科技飞速发展的时代,安卓和苹果两大操作系统之间的较量从未停歇。虽然它们各自有着忠实的...
安卓原生系统文件夹,安卓原生系... 你有没有发现,每次打开安卓手机,里面那些文件夹就像是一个个神秘的宝箱,里面藏着各种各样的宝贝?今天,...
基于安卓系统的游戏开发,从入门... 你有没有想过,为什么安卓手机上的游戏总是那么吸引人?是不是因为它们就像是你身边的好朋友,随时随地都能...
安卓系统怎样装驱动精灵,安卓系... 你那安卓设备是不是突然间有点儿不给力了?别急,今天就来手把手教你如何给安卓系统装上驱动精灵,让你的设...
如何本地安装安卓系统包,详细步... 你有没有想过,把安卓系统装在你的电脑上,是不是就像给电脑穿上了时尚的新衣?想象你可以在电脑上直接玩手...
安卓12卡刷系统教程,体验全新... 你有没有发现,你的安卓手机最近有点儿不给力了?运行速度慢得像蜗牛,是不是也想给它来个“换血大法”,让...
安卓系统无法打开swf文件,安... 最近是不是发现你的安卓手机有点儿不给力?打开SWF文件时,是不是总是出现“无法打开”的尴尬局面?别急...
鸿蒙系统依赖于安卓系统吗,独立... 你有没有想过,我们手机里的那个鸿蒙系统,它是不是真的完全独立于安卓系统呢?这个问题,估计不少手机控都...
适合安卓系统的图片软件,精选图... 手机里堆满了各种美美的照片,是不是觉得找起来有点头疼呢?别急,今天就来给你安利几款超级适合安卓系统的...
阴阳师安卓系统典藏,探寻阴阳师... 亲爱的阴阳师们,你是否在安卓系统上玩得如痴如醉,对那些精美的典藏式神们垂涎欲滴?今天,就让我带你深入...
安卓系统有碎片化缺点,系统优化... 你知道吗?在手机江湖里,安卓系统可是个响当当的大侠。它那开放、自由的个性,让无数手机厂商和开发者都为...
安卓4系统手机微信,功能解析与... 你有没有发现,现在市面上还有很多安卓4系统的手机在使用呢?尤其是那些喜欢微信的朋友们,这款手机简直就...
鸿蒙系统是安卓的盗版,从安卓“... 你知道吗?最近在科技圈里,关于鸿蒙系统的讨论可是热闹非凡呢!有人说是安卓的盗版,有人则认为这是华为的...
安卓系统怎么剪辑音乐,轻松打造... 你是不是也和我一样,手机里存了超多好听的歌,但是有时候想给它们来个变身,变成一段专属的旋律呢?别急,...
怎么把安卓手机系统变为pc系统... 你有没有想过,把你的安卓手机变成一台PC呢?听起来是不是有点酷炫?想象你可以在手机上玩电脑游戏,或者...
手机怎么装安卓11系统,手机安... 你有没有想过,让你的手机也来个“青春焕发”,升级一下系统呢?没错,就是安卓11系统!这个新系统不仅带...
安卓系统如何拼网络,构建高效连... 你有没有想过,你的安卓手机是怎么和网络“谈恋爱”的呢?没错,就是拼网络!今天,就让我带你一探究竟,看...
安卓系统怎么看小说,轻松畅享电... 你有没有发现,手机里装了那么多应用,最离不开的竟然是那个小小的小说阅读器?没错,就是安卓系统上的小说...
车载安卓系统12v几安,12V... 你有没有想过,你的车载安卓系统升级到12V几安,竟然能带来如此翻天覆地的变化?没错,今天咱们就来聊聊...