Advertisement
Guest User

Untitled

a guest
Oct 21st, 2018
67
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 3.26 KB | None | 0 0
  1. 1.
  2. %nr = 6
  3.  
  4. %Olekumudeli algolek:
  5.  
  6. X0=[-0.02;0;0;0]; %olekumudeli algolek
  7.  
  8. X0 =
  9.  
  10.  
  11.  
  12. -0.0200
  13.  
  14. 0
  15.  
  16. 0
  17.  
  18. 0
  19.  
  20. %Olekumudeli soovitud lõppolek
  21.  
  22. Xs=[0;0;nr/10;0]; %olekuvektori seadesuurused
  23.  
  24. Xs =
  25.  
  26.  
  27.  
  28. 0
  29.  
  30. 0
  31.  
  32. 0.6000
  33.  
  34. 0
  35.  
  36.  
  37.  
  38. %Ülesandest lähtuvad nõuded:
  39.  
  40. Täpsuskoridor - Tapsus +-5%
  41.  
  42. Umax =
  43.  
  44.  
  45.  
  46. 4.7000
  47.  
  48. 2.
  49. Diskreetimise sammu valisime siis kui süsteem oli juba häälestatud ja võimalikult kiire pidevajasüsteem.
  50.  
  51. Milline on aeg mille jooksul süsteem läheb täpsuskoridori sisse ja jääb sellesse piiri.
  52.  
  53. t_reg=1.5; - iseloomustab süsteemi kiirust.
  54.  
  55. Diskreetimise samm peab olema väiksem või võrdne kui reguleerimisaeg jagatud 3 süsteemi järguga.
  56.  
  57. Antud juhul siseolekuid on 4. Nurk, nurkkiirus, asend ja liikumisekiirus.
  58.  
  59. td = t_reg/3^4
  60.  
  61.  
  62.  
  63. td=4 %td<=tregf(3-n) [Ad. Bd]=c2d(A.B.td) %diskreetaja mudel saadakse pidevaja süsteemi koventeerimisl. kasutades pidevat A ja B maatriksit ning diskreetimissammu [Ad. Gd]=c2d(A.G.td) Diskreetimistakt määrab ära. kui tihedalt antakse meile siirete muutujate kohta infot. Pigem olgu see väike kui suur. muidu muutub sys ebastabiilseks. Diskreetimissammuks on valitud 6 s. sest selle perioodi jooksul jõuab mudel toimida samamoodi nagu pidevaja mudel ning diskreetaja mudel jääb lubatud vigade ja piirangute piiressse. Diskreetaja ja pidevaja mudelite õigsust saab hinnata simulinki skeemida ja graafikute põhjal. Meie süsteemile arvutatud graafik iseloomustab näitajate muutumist nii nagu oleme ette seadnud. seega sobib see kirjeldama meie süsteemi. Diskreetaja ja pidevaja süsteemi graafikud on küllaltki samased. seega võib eeldada. et toimivad õigesti.
  64.  
  65. Z=exp([P: p3]'td) %omaväärtused teisendatakse diskreetaja tasapinnale Kd=place(Ad. Bd. Z) %disk.aja tagasiside maatriks
  66.  
  67. 3.
  68. ksii=0.9;
  69.  
  70. t_reg=1.5;
  71.  
  72.  
  73.  
  74. wn=5/(ksii*t_reg);
  75.  
  76. kar_pol=[1 2*ksii*wn wn*wn];
  77.  
  78. P=roots(kar_pol);
  79.  
  80. p3=-abs(P(1))*1.5;
  81.  
  82. p4=p3*1.05;
  83.  
  84. K=place(A,B,[P;p3;p4]);
  85.  
  86.  
  87.  
  88. K =
  89.  
  90.  
  91.  
  92. -18.7762 -1.8780 -3.8283 -3.7559
  93.  
  94.  
  95.  
  96. Alustame soovitud reguleerimisajast. See on soovitav järk. Ksii 0-1ni. Kui Ksii on 1 lähedal, siis on süsteem sujuvam. Kui 0le lähemal, siis tekivad võnked.
  97.  
  98. Leiame teist järku karakteristliku polünoomi juured. Need on juhtimissüteemi poolused. Kuna süsteem on neljandat järku, peab sell olema 4 poolust. Peame valima veel kaks.
  99.  
  100. Arvutame tagasiside maatriksi. Sünteesime diskreetaja olekuregulaatori.
  101.  
  102. td=0.1; %diskreetaja periood
  103.  
  104. [Ad, Bd]=c2d(A, B, td); %süsteemi diskreetimine
  105.  
  106. [Ad, Gd]=c2d(A, G, td);
  107.  
  108.  
  109.  
  110. Z=exp([P; p3; p4]*td); % diskreetaja süsteemi omaväärtus
  111.  
  112. Kd=place(Ad, Bd, Z); %diskreetaja olekuregulaator
  113.  
  114.  
  115.  
  116. Kd =
  117.  
  118.  
  119.  
  120. -13.4043 -1.3698 -1.5381 -1.9371
  121.  
  122.  
  123.  
  124. ksii=0.7 %sumbuvus. Mida suurem see on, seda vähem tekib võnkeid. Suurte vöngetega jö'utakse tasakaaluasendisse aeglasemalt. wn=5/(ksirtreg) %omavõnkesagedus nim=[1 2`ksirwn wn*wn] %karakteristlik polünoom %olekumudel: A1-11-1-1 0 0;
  125.  
  126. sigma 0:
  127.  
  128. 0 1 0]
  129.  
  130. B=[ 0: 0] P=roots(nim) %soovitud suletud süsteemi pooluste paiknews p3=-wn %miinusmärgiga kindlasti: esialgne valik: muudab süsteemi kas kiiremaks völ aeglasemaks K=place(ABIP: p3]) %olekuregulaatori süntees
  131.  
  132. K =
  133.  
  134. 0.1774 1 139 0 0666
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement