% Messpraktikum
% Regelungstechnik
% M.Hess 14.11.2000

% Definition der Konstanten und Variablen
Ue = 1;						    % EingangsSpannung in Volt 
Ua = 0;      				    % AusgangsSpannung in Volt 
R1 = 10e3;						 % EingangsWiderstand vor OP-Amp in Ohm
R2 = 75e3;						 % RückkoppelungsWiderstand in Ohm
C1 = 0.1e-6;					 % EingangsKondensator in Farad

%Messwerte in Sekunden und Volt
TAnstieg = [0 0.0005 0.001 0.0015 0.002 0.0025 0.003 0.0035 0.004 0.005 0.006 0.007 0.008];
Anstieg = [0.33 0.215 0.122 0.037 -0.062 -0.103 -0.15 -0.210 -0.221 -0.274 -0.29 -0.320 -0.34];
TSprung = [0 0.0005 0.001 0.0015 0.002 0.0025 0.003 0.004 0.005 0.006];
Sprung = [-6.6 -5.5 -4.4 -3.5 -2.8 -2.3 -1.8 -1.15 -0.75 -0.47];

% Berechnung Zeitkonstanten
T1 = R1 * C1                %T1, T2 sowie Kp werden ausgegeben
T2 = R2 * C1
t1 = linspace(0,0.01,100);  %bildet Array t1
t2 = linspace(0,0.005,100); %bildet Array t2

% Berechnung VerstärkungFaktor Kp
Kp = - R2/R1

% Festlegungen für das Grafikfenster Sprung- Anstiegsantwort
figure('Name','Testfunktion Messpraktikum','NumberTitle','off'); %Anpassen des Figure-windows
clf;
hold; 

% Zeichnen der Sprungantwort
subplot(2,1,1);
plot(t1, Kp *(Ue*(T2/T1)*exp(-t1./T1)),TSprung, Sprung);
axis([0 0.01 -40 2]);       %Achsen der Übersichtlichkeit halber fix
grid;
title('Sprungantwort','FontWeight','bold');
xlabel('Zeit [s]');
ylabel('Ua [V]');
legend('berechnet','gemessen')

% Zeichnen der Anstiegsantwort
subplot(2,1,2);
plot(t2, Kp*(Ue*T2*(1-exp(-t2./T1))),TAnstieg, Anstieg);
axis([0 0.005 -0.4 0.4]);       %Achsen der Übersichtlichkeit halber fix
grid; 
title('Anstiegsantwort','FontWeight','bold');
xlabel('Zeit [s]');
ylabel('Ua [V]');
legend('berechnet','gemessen');