% Serieschwingkreis
% Regelungstechnik
% M.Hess 2.11.2000

%Benutzerabfrage
clc;                             %Löschen des command window
clear;                           %Löschen allfälliger Variabeln
%R = input('R [Ohm]= ');         %Variableneingabe
%L = input('L [uH]= ');
%C = input('C [uF]= ');
%A = input('Amplitude [V] = ');
R=50;                            %zum debuggen!
L=2000;
C=0.1;
a=10;
Kp = 1;

%Berechnungen
L = L*10^-6;                    %Multiplizieren mit u-Faktor
C = C*10^-6;                    %dito für C
T = sqrt(L*C);                  %Zeitkonstante
D = (R*C)/(2*T)                 %Dämpfungsgrad
w0 = 1/T                        %ungedämpfte Eigenfrequenz
we = w0*sqrt(1-D*D)             %gedämpfte Eigenfrequenz

%Grafik
figure('Name','Serieschwingkreis','NumberTitle','off');% Festlegen der Grafik
clf;                            %Löschen des Grafikfensters
grid;                           %Gitternetz
hold;                           %Plot einfrieren
t=0:(T/10):20*T;                %t von 0 bis 20*T, 200 Schritte
y=Kp*(1-(exp((-1).*D.*w0.*t)./sqrt(1-D*D)).*sin(w0*sqrt(1-D*D).*t+acos(D)));
plot(t,y);                      %plotten
title('Serieschwingkreis');     %Titel, Achsen inkl Werte, etc.
xlabel('Zeit [s]');
ylabel('Spannung [V]');