% Berechnung der Ausgangsspannung des LCR-Spannungsteilers zu Aufg. 7.1

clear;
L  = 0.1;  % Induktivitaet in Henry
R  = 116;  % Innenwiderstand der Kapazitaet in Ohm
C  = 1E-6; % Kapazitaet in Farad
N  = 4*500;% Abtastwerte je Periode (Sample Rate)
NP = 3;    % Anzahl der simulierten Perioden
TP = 2E-2; % Periodenlaenge in s
dt=TP/N;

% Erzeugung einer Periode der Eingangsspannungen
ue1P = 5*[ones(1, N/2) zeros(1, N/2)];
ue2P = 4*[(-N/4:N/4-1) (N/4:-1:-N/4+1)]/N;
ue3P = sin(2*pi*(0:N-1)/N);

% Zusammensetzten von NP Perioden des Eingabesignals
ueP = ue2P;% Signalauswahl
ue = [];
for idx=1:NP	 
  ue = [ue ueP];
end;

% Erzeugung der Zeitfolge
t=dt*(1:length(ue));


iL=0;     % Anfangswerte
uC=0;

for idx=1: length(ue)
    uL      = ue(idx) - uC - R*iL;
    ua(idx) = ue(idx)-uC;
    uC      = uC+dt/C*iL;
    iL      = iL+dt/L*uL;
end;

plot(t, ue, 'LineWidth', 3, t, ua, 'LineWidth', 3);
set(gca,'FontSize',24);
xlabel('t in s'); ylabel('u in V'); 
title('Aufgabe 7.1: LCR-Spannungsteilers');