Ich hab mal versucht, den zweiten Beleg komplett von MATLab rechnen zu lassen. Die Ergebnisse von Aufgabe 1 sehen (für mein beschränktes Verständnis) recht vernünftig aus, Aufgabe 2 is aber irgendwie schwer, weil da Polynome dritten Grades zu lösen sind. Hat da jemand ne Idee?
Code: Alles auswählen
%Bitte ändern:
Kw=10^4;
Kb=0.9;
Nr1=18; %1. Buchstabe Familienname
Nr2=17; %2. Buchstabe Vorname
Nr3=0; %3. Buchstabe Vorname
Nr4=4; %4. Buchstabe Familienname
%aus Aufgaben:
Kp=1;
Ki=1;
%Berechnung/Ausgabe Systemparameter:
Kt=0.9^Nr1
Tr=0.95^Nr2 * 100
Tb=1.05^Nr3 * 2
v0=1.1^Nr4 * 2
%1. Bodediagramm proportional:
Z1=[Kb*(Kp*Kw*Kt+v0)];
N1=[Tr*Tb Tr+Tb 1];
g1=tf(Z1,N1);
%Durchtrittsfrequenz
wd1=3^(1/2)/(6*Tr*Tb)*(Tr+Tb+(12*Tr*Tb+(Tr+Tb)^2)^(1/2))
%Proportionalregler (zur Kontrolle mal in dB umrechnen)
Kp1=1/(Kb*Kw*Kt)*(((Tr+Tb)^2*wd1^2+(Tr*Tb*wd1^2-1)^2)^(1/2)-Kb*v0)
%Anstiegszeit
tr1=1.5/wd1
%Regelabweichung
e_infty1=1/(1+Kb*Kp*Kw*Kt+Kb*v0)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5
%%%%%%Aufgabe 2%%%%%%%%%%%%%%%%5
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5
%2. Bodediagramm integral:
Z2=[Kb*v0 Kb*Ki*Kw*Kt];
N2=[Tr*Tb Tr+Tb 1 0];
g2=tf(Z2,N2);
%Durchtrittsfrequenz (wahrscheinlich falsch)
p=[6*Kb*v0*Tr*Tb 2*Kb*v0*3^(1/2)*(Tr+Tb)-3*Kb*v0*pi*(Tr+Tb)-2*Kb*Ki*Kw*Kt*Tr*Tb*3^(1/2)+3*pi*Kb*Ki*Kw*Kt*Tr*Tb 6*Kb*Ki*Kw*Kt-6*Kb*v0 -3*Kb*Ki*Kw*Kt*pi+2*Kb*Ki*Kw*3^(1/2)];
wd2=real(roots(p))
%Proportionalregler (keine Ahnung)
%Anstiegszeit
tr2=1.5/wd2
%Regelabweichung
e_infty2=0
%Ausgabe der Bodediagramme
Bode(g1,g2)