„SzabtechLabZH” változatai közötti eltérés

275. sor:

'''Matlab:'''

~~<pre>~~

s=zpk('s');

z=zpk('z',0.5);

P=2/((1+s)*(1+4*s)*(1+8*s)) % átviteli függvény

ts=0.5; % mintaveteli ido

Pz=c2d(P,ts); % diszkrét átviteli függvény

Pz % zérus-pólus alak

kc=1; %egység körerősítés először

[zd,pd,kd] = zpkdata(Pz,'v');%diszkrét folyamat zérusai pólusai

Cz=kc*(((z-pd(1))*(z-pd(2)))/((z-1)*z)); %diszkrét szab. átv. fv.

Lz = Cz*Pz

Lz=minreal(Lz,0.001);

w=logspace(-1,1,500); %saját logaritmikus skála beállítása

[mag,phase]=bode(Lz,w); %diszkrétből számolt bode diagram

[kc,pm,wg,wc]=margin(mag,phase-60,w); %60 fokos fázistöbblet

Cz=Cz*kc;

Lz=Cz*Pz; %a nyitott kör diszkrét átv. fv.

Lz=minreal(Lz,0.001);

Tz= Lz/(1+Lz); %diszkrét zárt rendszer átv. fv.

Ts=d2c(Tz,'tustin') %folytonos zárt rendszer átv. fv.

figure(2),step(Ts),grid %zart rendszer kimeneti jele

</~~pre>~~

Eredmény:

Zero/pole/gain:

0.0625

------------------------

(s+1) (s+0.25) (s+0.125)

Zero/pole/gain:

0.0011005 (z+3.157) (z+0.2247)

--------------------------------

(z-0.9394) (z-0.8825) (z-0.6065)

Sampling time (seconds): 0.5

Zero/pole/gain:

0.0011005 (z-0.9394) (z-0.8825) (z+3.157) (z+0.2247)

----------------------------------------------------

z (z-1) (z-0.9394) (z-0.8825) (z-0.6065)

Sampling time (seconds): 0.5

Zero/pole/gain:

0.0076717 s (s-7.709) (s-4) (s+4) (s+6.318) (s+0.9797)

------------------------------------------------------

s (s+4) (s+4.126) (s+0.9797) (s^2 + 0.7745s + 0.3623)

http://i.imgur.com/tq7wtWn.png

==== c/ ====

@@ 275. sor: / 275. sor: @@
 '''Matlab:'''
-<pre>
+  s=zpk('s');
-s=zpk('s');
+  z=zpk('z',0.5);
-z=zpk('z',0.5);
+  P=2/((1+s)*(1+4*s)*(1+8*s))  % átviteli függvény
-P=2/((1+s)*(1+4*s)*(1+8*s))  % átviteli függvény
+  ts=0.5;                      % mintaveteli ido
-ts=0.5;                      % mintaveteli ido
+  Pz=c2d(P,ts);                % diszkrét átviteli függvény
-Pz=c2d(P,ts);                % diszkrét átviteli függvény
+  Pz                           % zérus-pólus alak
-Pz                           % zérus-pólus alak
+  kc=1;                        %egység körerősítés először
-kc=1;                        %egység körerősítés először
+  [zd,pd,kd] = zpkdata(Pz,'v');%diszkrét folyamat zérusai pólusai
-[zd,pd,kd] = zpkdata(Pz,'v');%diszkrét folyamat zérusai pólusai
+  Cz=kc*(((z-pd(1))*(z-pd(2)))/((z-1)*z));   %diszkrét szab. átv. fv.
-Cz=kc*(((z-pd(1))*(z-pd(2)))/((z-1)*z));   %diszkrét szab. átv. fv.
+  Lz = Cz*Pz
-Lz = Cz*Pz
+  Lz=minreal(Lz,0.001);
-Lz=minreal(Lz,0.001);
+  w=logspace(-1,1,500);        %saját logaritmikus skála beállítása
-w=logspace(-1,1,500);        %saját logaritmikus skála beállítása
+  [mag,phase]=bode(Lz,w);      %diszkrétből számolt bode diagram
-[mag,phase]=bode(Lz,w);      %diszkrétből számolt bode diagram
+  [kc,pm,wg,wc]=margin(mag,phase-60,w);     %60 fokos fázistöbblet
-[kc,pm,wg,wc]=margin(mag,phase-60,w);     %60 fokos fázistöbblet
+  Cz=Cz*kc;
-Cz=Cz*kc;
+  Lz=Cz*Pz;                    %a nyitott kör diszkrét átv. fv.
-Lz=Cz*Pz;                    %a nyitott kör diszkrét átv. fv.
+  Lz=minreal(Lz,0.001);
-Lz=minreal(Lz,0.001);
+  Tz= Lz/(1+Lz);               %diszkrét zárt rendszer átv. fv.
-Tz= Lz/(1+Lz);               %diszkrét zárt rendszer átv. fv.
+  Ts=d2c(Tz,'tustin')          %folytonos zárt rendszer átv. fv.
-Ts=d2c(Tz,'tustin')          %folytonos zárt rendszer átv. fv.
+  figure(2),step(Ts),grid      %zart rendszer kimeneti jele
-figure(2),step(Ts),grid      %zart rendszer kimeneti jele
-</pre>
+Eredmény:
+  Zero/pole/gain:
+.0625
+  ------------------------
+  (s+1) (s+0.25) (s+0.125)
+  Zero/pole/gain:
+.0011005 (z+3.157) (z+0.2247)
+  --------------------------------
+  (z-0.9394) (z-0.8825) (z-0.6065)
+  Sampling time (seconds): 0.5
+  Zero/pole/gain:
+.0011005 (z-0.9394) (z-0.8825) (z+3.157) (z+0.2247)
+  ----------------------------------------------------
+        z (z-1) (z-0.9394) (z-0.8825) (z-0.6065)
+  Sampling time (seconds): 0.5
+  Zero/pole/gain:
+.0076717 s (s-7.709) (s-4) (s+4) (s+6.318) (s+0.9797)
+  ------------------------------------------------------
+  s (s+4) (s+4.126) (s+0.9797) (s^2 + 0.7745s + 0.3623)
+http://i.imgur.com/tq7wtWn.png
 ==== c/ ====