%Plotten normierter Systemantwort eines Einfreiheitsgradsystem
clc,
color=['y','m','c','r','g','b','y'];
% Normierte ungedaempfte Eigenfrequenz
OMEGA=1.0;
% Daempfungswerte
zeta=[0.1,0.2,0.4,0.6,0.8,1.0];
% Normierte Frequenzen von 0 bis 6
rinit = 0;
rfinal = 6;
r = linspace(rinit,rfinal,512);
for ii=1:6
% Zaehler
Num=r.*r;
omega=OMEGA*sqrt(1.0-zeta(ii)*zeta(ii));
% Nenner
Den=(1-r.*r)+j*(2*zeta(ii).*r);
% Verhaeltnis
TR=Num./Den;
% Aufzeichnen des Amplitudenganges
fig1=figure(1);
if ii==1,clc,end;
semilogy(r,abs(TR),color(ii))
pause
hold on
% Aufzeichnen des Phasenganges
fig2=figure(2);
if ii==1,clc,end;
plot(r,angle(TR),color(ii))
pause
hold on
end
fig1=figure(1);
xlabel('Normierte Frequenz (r)')
ylabel('Amplitudengang')
title('Normierte Eigenschaften des Einfreiheitsgrad-Systems')
grid
fig2=figure(2);
xlabel('Normierte Frequenz (r)')
ylabel('Phasengang')
title('Normierte Eigenschaften des Einfreiheitsgrad-Systems')
grid
