clf()//löschen aller vorrigen Graphikfenster
scf(0)//Öffnen eines neuen Graphik Fensters
s=poly(0,'s')//Definition von s als Polynom
G1=syslin('c',7/(2*s))
frq=logspace(-3,2,100000); //frq=frequencies in Hz
[rep] =repfreq(G1,frq);
[db,phi]=dbphi(rep);
// draw frequency characteristics
subplot(3,3,1); 
plot2d(frq*2*%pi,db,style=5,logflag="ln");
xtitle('I-Element','w[rad/sec]','dB')//Achsenbeschriftung
h1=legend("G1=7/(2*s)");//Beschriftung der Legende
subplot(3,3,4);
plot2d(frq*2*%pi,phi,style=5,logflag="ln");
xtitle(' ','w[rad/sec]','°')//Achsenbeschriftung
h1=legend("G1=7/(2*s)");//Beschriftung der Legende


t=[0:0.001:2];
y=csim('step',t,G1);
subplot(3,3,7); 
plot2d(t,y)//Sprungantwort
xtitle(' ','A','E')//Achsenbeschriftung


G2=syslin('c',(0.75*(1+3*s))/(1+s*0.01))
frq=logspace(-3,2,100000); //frq=frequencies in Hz
[rep] =repfreq(G2,frq);
[db,phi]=dbphi(rep);
// draw frequency characteristics
subplot(3,3,2);
plot2d(frq*2*%pi,db,style=5,logflag="ln");
xtitle('PD-Element','w[rad/sec]','dB')//Achsenbeschriftung
h2=legend("G2=0.75*(1+3*s))/(1+s*0.01)");//Beschriftung der Legende
subplot(3,3,5);
plot2d(frq*2*%pi,phi,style=5,logflag="ln");
xtitle(' ','w[rad/sec]','°')//Achsenbeschriftung
h2=legend("G2=0.75*(1+3*s))/(1+s*0.01)");//Beschriftung der Legende


t=[0:0.001:2];
y=csim('step',t,G2);
subplot(3,3,8); 
plot2d(t,y)//Sprungantwort
xtitle(' ','A','E')//Achsenbeschriftung

G3=syslin('c',(5/s)*(0.75*(1+3*s))/(1+s*0.01))
frq=logspace(-3,2,100000); //frq=frequencies in Hz
[rep] =repfreq(G3,frq);
[db,phi]=dbphi(rep);
// draw frequency characteristics
subplot(3,3,3);
plot2d(frq*2*%pi,db,style=5,logflag="ln");
xtitle('I-Element*PD-Element','w[rad/sec]','dB')//Achsenbeschriftung
h3=legend("G2=(5/s)*(0.75*(1+3*s))/(1+s*0.01)");//Beschriftung der Legende
subplot(3,3,6);
plot2d(frq*2*%pi,phi,style=5,logflag="ln");
xtitle(' ','w[rad/sec]','°')//Achsenbeschriftung
h3=legend("G2=(5/s)*(0.75*(1+3*s))/(1+s*0.01)");//Beschriftung der Legende

t=[0:0.001:2];
y=csim('step',t,G3);
subplot(3,3,9);
plot2d(t,y)//Sprungantwort
xtitle(' ','A','E')//Achsenbeschriftung