„Labor ZH feladatai témakörök szerint csoportosítva” változatai közötti eltérés

Harapeti (vitalap | szerkesztései)
EREDMÉNY kommentben
Kosa333 (vitalap | szerkesztései)
Feladat szöveg fix, kérésre hozzáadott megyarázattal
 
(19 közbenső módosítás, amit 4 másik szerkesztő végzett, nincs mutatva)
18. sor: 18. sor:
    
    
   [Ad,bd,cd,dd]=canon(A,b,c,d)
   [Ad,bd,cd,dd]=canon(A,b,c,d)
  eig(A)


Eredmény:
Eredmény:
41. sor: 42. sor:
==== b./ Irányítható-e, megfigyelhető-e a rendszer? (2 pont) ====
==== b./ Irányítható-e, megfigyelhető-e a rendszer? (2 pont) ====
--> irányítható, megfigyelhető
--> irányítható, megfigyelhető
rank(ctrb(A,b))
--> 2, tehát irányítható (megj.: ha a rang megegyezik A mx dimenziójával akkor irányítható)
rank(obsv(A,c))
--> 2, tehát megfigyelhető (megj.: ha a rang megegyezik A mx dimenziójával akkor megfigyelhető, itt n=2)


==== b./ Ábrázolja az eredeti rendszer állapottrajektóriáját u(t) = 0 és x(0)=[x_1(0);x_2(0)]=[2;6] felételek mellett. (3 pont) ====
==== b./ Ábrázolja az eredeti rendszer állapottrajektóriáját u(t) = 0 és x(0)=[x_1(0);x_2(0)]=[2;6] felételek mellett. (3 pont) ====
66. sor: 75. sor:
   H=ss(A,b,c,d)
   H=ss(A,b,c,d)
   H=zpk(H)
   H=zpk(H)
  eig(A)


Eredmény:
Eredmény:
91. sor: 101. sor:


Rendszer pólusai: 0, -2
Rendszer pólusai: 0, -2
Átviteli fv. pólusok: 0
 
Labilis az integrátor miatt
Az hogy stabil-e az passz, a 0 miatt a stabilitás határán van.
b(1)=0 miatt nem irányítható, de megfigyelhető --> ??????? b(1) nem 2.8284 ???
 
rank(ctrb(A,b))
 
--> 1, tehát nem irányítható (megj.: ha a rang megegyezik A mx dimenziójával akkor irányítható, itt n=2, 1<2)
 
rank(obsv(A,c))
 
--> 2, tehát megfigyelhető (megj.: ha a rang megegyezik A mx dimenziójával akkor megfigyelhető, itt n=2, 2=2 --> IGEN)  




117. sor: 134. sor:


   rank(ctrb(A,b))
   rank(ctrb(A,b))
--> 3, tehát irányítható
--> 3, tehát irányítható (megj.: ha a rang megegyezik A mx dimenziójával akkor irányítható)


   rank(obsv(A,c))
   rank(obsv(A,c))
--> 2, tehát NEM megfigyelhető
--> 2, tehát NEM megfigyelhető (megj.: ha a rang megegyezik A mx dimenziójával akkor megfigyelhető, itt n=3, 2<3 --> NEM)




264. sor: 281. sor:
--> NEM stabil, mivel a 3. pólus pozitív!
--> NEM stabil, mivel a 3. pólus pozitív!


==== b./ Tervezzen állapot-visszacsatolásos szabályozást úgy. hogy a zárt rendszer egy másodrendű lengő tagból és egy egytárolós tagból álljon. A lengő tag csillapítási tényezője 0.6 és időállandója 0.5 legyen. Határozza meg az alapjelkövetéshez a statikus kompenzációs tényező értékét is.  (4 pont) ====
==== b./ Tervezzen állapot-visszacsatolásos szabályozást úgy. hogy a zárt rendszer egy másodrendű lengő tagból és egy egytárolós tagból álljon. A lengő tag csillapítási tényezője 0.6 és időállandója 0.5 legyen. Az egytárolós tag időállandója legyen 2. Határozza meg az alapjelkövetéshez a statikus kompenzációs tényező értékét is.  (4 pont) ====


   T0=0.5
   T0=0.5
270. sor: 287. sor:
   den=[T0*T0, 2*T0*kszi, 1]
   den=[T0*T0, 2*T0*kszi, 1]
   pc=roots(den)
   pc=roots(den)
   pc(3)=-1/2
   pc(3)=-1/2   %T1=2, pc(3)=roots([T1, 1]) <- Az egytárolós tag gyöke [1/(1+T*s)]
   k=acker(A,b,pc)
   k=acker(A,b,pc)
   kr=1/dcgain(A-b*k,b,c,d)
   kr=1/dcgain(A-b*k,b,c,d)
412. sor: 429. sor:
   A=m*Au
   A=m*Au


   m =
   %    m =
      0.3508
  %        0.3508
    
   %   
   f =
   %    f =
  -105.2551
  %    -105.2551
    
   %   
   fi =        % mo!
   %    fi =        % mo!
  -219.8467
  %    -219.8467
    
   %   
   A =        % mo!
   %    A =        % mo!
      0.7016
  %        0.7016


<hr />
<hr />
436. sor: 453. sor:
   A=10*m
   A=10*m


   m =
   %    m =
      0.1085
  %        0.1085
   f =
   %    f =
    -139.3987
  %      -139.3987
    
   %   
   fid =      % mo!!
   %    fid =      % mo!!
    -253.9903
  %      -253.9903
    
   %   
   A =        % mo!!
   %    A =        % mo!!
    1.0847
  %      1.0847




458. sor: 475. sor:
==== a./  Zérusrendű tartószerv esetén adja meg a szakasz G(z) impulzusátviteli függvényét zérus-pólus alakban.  (3 pont) ====
==== a./  Zérusrendű tartószerv esetén adja meg a szakasz G(z) impulzusátviteli függvényét zérus-pólus alakban.  (3 pont) ====


  s=zpk('s');
   P=2/( s*(1+2*s) )
   P=2/( s*(1+2*s) )
   Ts=0.5
   Ts=0.5
466. sor: 484. sor:
   Gz=G1z/(z^d)
   Gz=G1z/(z^d)


   d=2
   %    d=2
    
   %   
   Zero/pole/gain:        %% mo!   
   %    Zero/pole/gain:        %% mo!   
   G(z) =
   %    G(z) =
  0.1152 (z+0.9201)
  %    0.1152 (z+0.9201)
   --------------------
   %    --------------------
   z^2 (z-1) (z-0.7788)
   %    z^2 (z-1) (z-0.7788)


==== b./  A szabályozó impulzusátviteli függvénye C(z) = 0.5*( (z-z_1)/z ). Póluskiejtéses kompenzáció esetén adja meg z_1, értékét. Milyen típusú szabályozót valósítottunk meg? (2 pont) ====
==== b./  A szabályozó impulzusátviteli függvénye C(z) = 0.5*( (z-z_1)/z ). Póluskiejtéses kompenzáció esetén adja meg z_1, értékét. Milyen típusú szabályozót valósítottunk meg? (2 pont) ====
490. sor: 508. sor:
Stabilis: fázistartalék > 0. (Lz amúgy nem stabil (lásd step(Lz), csak így visszacsatolva lesz.)
Stabilis: fázistartalék > 0. (Lz amúgy nem stabil (lásd step(Lz), csak így visszacsatolva lesz.)


   ud =        % mo!
   %    ud =        % mo!
      0.5000
  %        0.5000
      0.1106
  %        0.1106
      0.1106
  %        0.1106
      0.0818
  %        0.0818
      0.0489
  %        0.0489
      0.0367
  %        0.0367


Érdekes, itt a mintamegoldás szerint ennek kell kijönnie:
Érdekes, itt a mintamegoldás szerint ennek kell kijönnie:
   ud[1:5] = 2.0000, 0.4424, 0.4424, -0.0184, -0.5443
   ud[1:5] = 2.0000, 0.4424, 0.4424, -0.0184, -0.5443


--> ???
--> ???
520. sor: 538. sor:
   Gz=G1z/(z^d)
   Gz=G1z/(z^d)


   d=2
   %    d=2
    
   %   
   Zero/pole/gain:        %% mo!   
   %    Zero/pole/gain:        %% mo!   
   G(z) =
   %    G(z) =
    0.13417 (z+0.8008)
  %      0.13417 (z+0.8008)
   -------------------------
   %    -------------------------
   z^2 (z-0.8465) (z-0.6065)
   %    z^2 (z-0.8465) (z-0.6065)


==== b./  A szabályozó impulzusátviteli függvénye C(z) = 0.25*( (z-z_1)/(z-1) ). Póluskiejtéses kompenzáció esetén adja meg z_1, értékét. Milyen típusú szabályozót valósítottunk meg? (2 pont) ====
==== b./  A szabályozó impulzusátviteli függvénye C(z) = 0.25*( (z-z_1)/(z-1) ). Póluskiejtéses kompenzáció esetén adja meg z_1, értékét. Milyen típusú szabályozót valósítottunk meg? (2 pont) ====
539. sor: 557. sor:
   [gm,pm]=margin(Lz)
   [gm,pm]=margin(Lz)
    
    
   gm =      % mo.!!
   %    gm =      % mo.!!
    3.0568
  %      3.0568
  %
  %    pm =      % mo.!!
  %    52.6390


  pm =      % mo.!!
  52.6390
--> stabilis.      % mo.!!
--> stabilis.      % mo.!!


561. sor: 580. sor:
http://i.imgur.com/h3m8ido.png
http://i.imgur.com/h3m8ido.png


u(0) = 0.25
  %  u(0) = 0.25
u(végtelen) = 0.25
  %  u(végtelen) = 0.25




578. sor: 597. sor:
a./ Adja meg a rendszer vágási körfrekvenciáját, fázistartalékát és erősítési tartalékát. Stabilis-e a zárt rendszer?
a./ Adja meg a rendszer vágási körfrekvenciáját, fázistartalékát és erősítési tartalékát. Stabilis-e a zárt rendszer?
Egységugrás zavarójelre és zérus alapjel esetén:
Egységugrás zavarójelre és zérus alapjel esetén:
b./ Ábrázolja minőségileg helyesen az>’ kimenőjel időbeli lefolyását,
b./ Ábrázolja minőségileg helyesen az y kimenőjel időbeli lefolyását,
c./ Adja meg a kimenőjel és a beavatkozójel állandósult értékét.
c./ Adja meg a kimenőjel és a beavatkozójel állandósult értékét.


592. sor: 611. sor:
http://i.imgur.com/k0MFBzL.png
http://i.imgur.com/k0MFBzL.png


   Gm=15.6dB
   Gm=15.6dB


   [gm,pm,wg,wc]=margin(L)
   [gm,pm,wg,wc]=margin(L)


   gm=6, pm=43.2099, wc=0.7793rad/sec
   gm=6, pm=43.2099, wc=0.7793rad/sec


Mivel pm>0, a szabályozás stabilis.
Mivel pm>0, a szabályozás stabilis.
606. sor: 625. sor:
   grid
   grid


   y_vég=0,
   % y_vég=0,
   u_vég=-1
   % u_vég=-1


http://i.imgur.com/ky0WOL8.png
http://i.imgur.com/ky0WOL8.png
614. sor: 633. sor:


=== II. 1. Adott az alábbi szabályozási kör: ===
=== II. 1. Adott az alábbi szabályozási kör: ===
C(s)=(1+10*s)/(10*s)
 
P(s)=1/(1+10*s)(1+s)(1+0.5*s)
  %  C(s)=(1+10*s)/(10*s)
  %  P(s)=1/(1+10*s)(1+s)(1+0.5*s)


http://i.imgur.com/pnitBve.png
http://i.imgur.com/pnitBve.png
637. sor: 657. sor:
http://i.imgur.com/Ml3h14J.png
http://i.imgur.com/Ml3h14J.png


   gm= 30 (29.5dB), pm=81.48, mt=0.89, stabilis (pm>0)
   gm= 30 (29.5dB), pm=81.48, mt=0.89, stabilis (pm>0)
 


   Tz=P/(1+L)
   Tz=P/(1+L)
646. sor: 665. sor:
   grid
   grid


   y_vég=0,
   y_vég=0,
   u_vég=-1
   u_vég=-1


http://i.imgur.com/p6IXH9U.png
http://i.imgur.com/p6IXH9U.png
673. sor: 692. sor:
   L=minreal(L)
   L=minreal(L)


a./ strukturálisan stabilis, kmax=inf
==== a./ strukturálisan stabilis, kmax=inf ====
b./  
==== b./ ====
   [gm,pm]=margin(L)
   [gm,pm]=margin(L)
   m=bode(L+1);
   m=bode(L+1);
   mt=min(m)
   mt=min(m)


   pm=62, mt=0.76, stabilis
   pm=62, mt=0.76, stabilis


c./
==== c./ ====
   H=minreal(1/(1+L))
   H=minreal(1/(1+L))
   step(H)
   step(H)
   grid on
   grid on
d./  
==== d./ ====
   T=minreal(L/(1+L))
   T=minreal(L/(1+L))
   R=1/(s+2)
   R=1/(s+2)
701. sor: 720. sor:
http://i.imgur.com/pnitBve.png
http://i.imgur.com/pnitBve.png


C(s)=(1+20*s)/(20*s)
  % C(s)=(1+20*s)/(20*s)
P(s)=10/( (1+20*s)*(1+2*s)*(1+s) )
  % P(s)=10/( (1+20*s)*(1+2*s)*(1+s) )


==== a./ Adja meg a rendszer erősítési tartalékát, fázistartalékát és modulus tartalékát. Stabilis-e a zárt rendszer? (3 pont) ====  
==== a./ Adja meg a rendszer erősítési tartalékát, fázistartalékát és modulus tartalékát. Stabilis-e a zárt rendszer? (3 pont) ====  
708. sor: 727. sor:
==== c./ r(t) = 0 és 0<=t<=100 (sebességugrás) alapjel és zérus zavarás esetén ábrázolja minőségileg egy koordináta-rendszerben az alapjelet és a kimenőjelet! Mekkora a statikus hiba? (3 pont) ====  
==== c./ r(t) = 0 és 0<=t<=100 (sebességugrás) alapjel és zérus zavarás esetén ábrázolja minőségileg egy koordináta-rendszerben az alapjelet és a kimenőjelet! Mekkora a statikus hiba? (3 pont) ====  


==== a ====  
==== a./ ====  


   s=zpk('s')
   s=zpk('s')
721. sor: 740. sor:
   mt=min(m)
   mt=min(m)
    
    
   gm=3 (9.5dB), pm = 32.6, mt=0.43, stabilis
   % gm=3 (9.5dB), pm = 32.6, mt=0.43, stabilis


==== b./ ====
==== b./ ====
729. sor: 748. sor:
   grid
   grid


   u_kezd = -1
   % u_kezd = -1
   u_vég = -0.1
   % u_vég = -0.1


==== c./ ====
==== c./ ====
749. sor: 768. sor:
mego.:
mego.:


   es=1/K=1/0.5=2
   % es=1/K=1/0.5=2


<hr />
<hr />
1 121. sor: 1 140. sor:


<hr />
<hr />
[[Category:Infoalap]]