Facebook
From janusze igrazyny, 6 Years ago, written in Plain Text.
Embed
Download Paste or View Raw
Hits: 238
  1. echo on
  2. clc
  3. clear
  4. close all
  5. %dane
  6. n= [2.5 0 1];
  7. m=[ 3 4.2 2];
  8. sys = tf(n,m)
  9. pause
  10.  
  11. %transmitancja G(s)
  12. disp('G(s)')
  13. printsys(n,m), sys
  14. pause
  15.  
  16. %funkcja tf2zp
  17. [z,b,w] =tf2zp(n,m)
  18. [nzp,mzp]=zp2tf(z,b,w)
  19. printsys(nzp,mzp)
  20. pause
  21.  
  22. %dunkcja residue()
  23. [r,p,k]=residue(n,m)
  24. [nres,mres]= residue(r,p,k)
  25. printsys(nres,mres)
  26. pause
  27.  
  28. %zestawienie
  29. [z,b]
  30. [r,p]
  31. [w,k]
  32. pause
  33.  
  34. %pd[pwoedzo ma step i pulse
  35. t = [0:0.01:10];
  36. yi=impulse(n,m,t);
  37. ys=step(n,m,t);
  38. ysres=step(nres,mres,t);
  39. figure(1), plot(t,yi,t,ys),
  40. grid, axis([0 10 -1.5 1.5])
  41. legend('impulse()','step()')
  42. echo off
  43.  
  44. /////////////////////////////////////////////////
  45.  
  46. echo on
  47. clc
  48. clear
  49. close all
  50.  
  51. %--------------------------------------
  52. %dane
  53. n1 = [4 5 6]; m1 = [1 5 4]
  54. printsys(n1,m1)
  55. n2=1;m2=[1 2.5];
  56. printsys(n2,m2);
  57. pause
  58.  
  59. %--------------------------------------
  60. %połączenie szeregowe
  61. [nsz,msz]= series9n1,m1,n2,m2);
  62. disp('gszer(s)')
  63. printsys(nsz,msz)
  64. pause
  65.  
  66. %--------------------------------------
  67. %połaczenie rownoległe
  68. [nr,mr]= parallel(n1,m1,n2,m2);
  69. disp('Grow(s)')
  70. printsys(nr,mr)
  71. pause
  72.  
  73. %--------------------------------------
  74. %uklad z ujemnym sprzezeniem zwrotnym
  75.  
  76. [nusz,musz]=feedback(n1,m1,n2,m2);
  77. disp('Gusz(s)')
  78. printsys(nusz,musz)
  79. pause
  80.  
  81. %--------------------------------------
  82. %residua
  83. [rsz,bsz,wsz] = residue(nsz,msz);
  84. [rr,br,wr] = residue(nr,mr);
  85. [rusz,busz,wusz]= residue(nusz,muszz);
  86. pause
  87.  
  88. %--------------------------------------
  89. %tf2zp
  90. [zsz,psz,ksz] = tf2zp(nsz,msz);
  91. [zr,pr,kr] = tf2zp(nr,mr);
  92. [zusz,pusz,kusz]= tf2zp(nusz,muszz);
  93. pause
  94.  
  95. tf2zp
  96. %odpowiedz na step
  97. t=[0:.01:6];
  98. ysz=step(nsz,msz,t);
  99. yr=step(nr,mr,t);
  100. yusz=step(nusz,musz,t);
  101. figure(1)
  102. plot(t,ysz,t,yr,t,yusz) ,grid, axis ([ 0 6 0 2.5])
  103. legend('ysz','yr','yusz')